A Tool for Converting FEM Models into Representations Suitable for Control Synthesis

Size: px
Start display at page:

Download "A Tool for Converting FEM Models into Representations Suitable for Control Synthesis"

Transcription

1 Proceedings o the 17th Word Congress The Internationa Federation o Automatic Contro Seou, Korea, Juy 6-11, 8 A Too or Converting FEM Modes into Representations Suitabe or Contro Synthesis Syed Fawad Raza Ai Bokhari, Sauat S. Chughtai, Herbert Werner Institute o Contro Systems, Hamburg University o Technoogy, Germany Abstract: Finite Eement Methods (FEM) pay an important roe in modeing o compex systems. Modes generated rom FEM are o very high dimension and are diicut to hande or contro system design. Inthis paper, an agorithm is presentedthat converts asystem generated rom FEM into the state space mode o an interconnected system, thus vasty reducing the compexity o synthesizing a distributed contro strategy. A simpe cantiever beam probem is used as an exampe to iustrate the proposed method. A state space representation o this system is obtained through FEM anaysis. This mode is converted rom umped state space orm to interconnected orm. Adecentraized controer is then designed usinga homotopy approach. Keywords: Sotware too, Interconnected Systems, Finite Eement Methods (FEM). 1. INTRODUCTION Engineering systems in which subsystems interact with each other are o considerabe practica interest. With advances in design techniques and computationa power compex system modes o very high dimensions can be deveoped, incuding umped approximations o partia dierentia equations (PDEs) exampes incudethedeection obeams, pates andmembranes, andthetemperature distribution o thermay conductive materias. Soutions to these types o probems are possibe or very simpe systems; or more reaistic systems the Finite Eement Method (FEM)(Becker, 4) is a very eective too or soving such probems numericay (Huebner et a., 1). Currenty FEM is used in structura anaysis, heat transer, uid mechanics etc. In FEM a compete system is divided into sma eements which are spatiay distributed and interact with each other, based on their spatia ocation. In structura anaysis behavior a system is described by the dispacement o eements satisying materia aws (constitutive equation) (Becker, 4). A eements are assembed together and the requirements o continuity and equiibrium are satisied between neighboring eements. A unique soution can be ound or a given boundary conditions. A state space representation obtained rom FEM generay has a arge number o states aong with arge number o inputs andoutputs. Forsuchsystems moderncontro design techniques may ai due to the size o the probem. One way o deaing with such systems is to reduce the size o system using standard mode order reduction techniques. However, this reduction may eiminate some aspects o the dynamic behavior which may get exited in the cosed oop, thus degrading the achieved perormance. In this paper we propose a dierent approach, which is to decompose the mode used or contro synthesis into subsystems. Each subsystem has ony a sma number o states and its own sensor and actuator. This decomposition can aciitate the use o modern controer synthesis techniques to arge systems, using the ramework o interconnected systems. A genera interconnected systems is composed o N subsystems interconnected with each other. The state space reaization o i th subsystem is given by ẋ i (t) = A i x i (t) + N j=1,j i B vij x j (t) + B i w i (t) z i (t) = C i x i (t) + D i w i (t) (1) where x i R ni, w i R ki, z i R oi are the state vector, disturbance inputs and outputs, respectivey. The matrices A i, B vij, B i, C i and D i are constant matrices o appropriate dimensions. Here, the matrix B vij represents the interconnection o the i th subsystem with the j th subsystem. Modes in the orm o (1) have been extensivey studied, and eicient synthesis toos have been proposed to design controers or such systems. Some earier work on such systems is presented in (Sijak, 1978), whie (Sijak and Zecevic, 5) presents some new approaches based on inear matrix inequaities. In (Chen et a., 5) the authors have presented a homotopy based approach to design robust controers or (1). In (D Andrea and Duerud, 3) the authors have shown that the Kaman- Yaqubovich-Popov emma can be generaized or acass o interconnected systems with identica subsystems having interconnections ony with their near neighbors, which can be considered as a specia case o system (1). In this paper we present a sotware too that can decompose an FEM generated mode into the orm o (1), /8/$. 8 IFAC /876-5-KR

2 17th IFAC Word Congress (IFAC'8) Seou, Korea, Juy 6-11, 8 thus making it suitabe or contro system design. The paper is organized as oows. A genera matrix decomposition unction is presented in Section. This unction is appied in Section 3 to a FEM generated system and a compete agorithm to decompose the given system is discussed. Insection4auniorm cantiever exampe is used to demonstrate the appication o the too. In section 5 the homotopy approach o (Chen et a., 5) is used to designdecentraized controer orthedecomposedsystem. Concuding remarks and uture directions are given in Section 5.. THE BASIC FUNCTION In this section the basic unction decomp() is introduced. The unction decomposes any matrix G R p q into the bock structure G 11 G 1 G 1n G 1 G G n G n1 G n G nn. () In the oowing G i wi represent the i th diagona bock G(i,i), where G i R ( hi ji) such that, n p = hi (3) q = i=1 n ji (4) i=1 The unction requires two urther inputs rom the user. These are T h = [ h1 h hn ] T (5) T j = [ j1 j jn ] T (6) These vectors speciythedimension oeach bock, andor the genera case ( hi ji ). Let us deine G vi as G vi = [ G (i,1) G (i,i 1), G (i,i+1) G (i,n) ] In order to extract G i and G vi rom G we need upper and ower row and coumn indices in terms o the eements o vectors h and j. Next we wi ind these indices or G i and G vir, where r = 1,...,i 1,i + 1,...,n..1 Indices o G i Let e h be the ower row index and e j be the ower coumn index o G i respectivey, then these can be ound rom vectors h and j as Then we have i 1 e h = hk + 1 (7) k=1 i 1 e j = jk + 1 (8) k=1 G i = G(e h+ i,e j+ g i ), where. Indices o G vir i =,1,...,( hi 1) g i =,1,...,( ji 1) (9) Let G vir = G( vir,g vir ), then vir = e h + i. For the case when r < i, i 1 g vir = e j jk + {,1,..., jr 1} k=r and or the case when r > i, r 1 g vir = e j + jk + {,1,..., jr 1} k=i The above approach is coded into a MATLAB unction decomp(g, h, j,i,m) where the binary input m can be used to get G i i m = or G vi i m. 3. APPLICATION TO FEM GENERATED SYSTEMS A system mode generated via FEM is given as M q + C q + Kq = F d + F (1) where q represents the generaized co-ordinates and M = Mass matrix K = Stiness matrix C = Damping matrix F d = Appied noda orces (externa orces) F = Boundary conditions These matrices are obtained by combining the equations o motion o each eement o the FEM anaysis. Foowing this the continuity constraints between neighboring eements are appied. The above equation can be written in state space orm as [ ] [ ] I x = M 1 K M 1 x + C M 1 F(t) }{{}}{{} Ã y = Cx + Du (11) where x R n, u R k, y R o. The state vector in (11) is arranged as x T = [ x T 1 x T x T n ẋ T 1 ẋ T ẋ T n ]T This can be transormed into x T = [ x T 1 ẋ T 1 x T ẋ T x T n ẋ T n ]T (1) Deining a unitary transormation matrix T r such that x = T r x we obtain A = T r ÃTr 1, B = T r B, C = 1 CT r and D = D. Next we show how to decompose the transormed system into the rom (1) with N subsystems. Let or the i th subsystem ẋ i,x i, R xi, z i R zi and w i R wi. Deine the vectors 667

3 17th IFAC Word Congress (IFAC'8) Seou, Korea, Juy 6-11, 8 F x = ( x1,..., xn ) (13) F z = ( z1,..., zn ) (14) F w = ( w1,..., wn ) (15) The unction decomp() canthen be appied toa, B, C and D to decompose them to get the state space representation (1). This is expressed in the oowing agorithm. 3.1 Agorithm For an FEM generated state space in the rom (11) Step 1 Find the transormation matrix T r and appy this simiarity transormation to get A, B, C, D. Step 3 Deine the vectors F x, F z and F w. Step Obtain the other matrices by A i = decomp(a,f x,f x,i,); B vi = decomp(a,f x,f x,i,1); B i = decomp(b,f x,f w,i,); C i = decomp(c,f z,f x,i,); D i = decomp(d,f z,f w,i,); 4. EXAMPLE: CANTILEVER BEAM This section briey presents the FEM treatment o a uniorm cantiever (Juang and Phan, 4). Fig. 1. Beam cantiever FEM o singe cantiever beam eement Firstconsider abeam with camped boundary condition at one end and ree boundary condition at the other as shown in Fig 1. In order to iustrate the generation o an FEM mode, a singe eement with two nodes o a camped-ree uniorm beam is taken. Let an impuse act downwards at the node as shown in the Fig 1. This orce resuts in a dispacement d(x,t) and an anguar rotation, θ = d(x,t) x. The equation o motion can be written as ρ d(x,t) t + EM I 4 d(x,t) x 4 = (16) Here ρ, E, M I are density, Young moduus o easticity and second moment o inertia, respectivey. The behavior o dispacement and rotation can be approximated interms o interpoation unctions as d 1 (t) θ d(x,t) = [ N 1 (x) N (x) N 3 (x) N 4 (x) ] 1 (t) d (t) θ (t) where the subscripts 1 and denote the irst and second node o the eement. The N i (x) are the interpoation unctions o the corresponding variabe, which can be approximated as N k (x) = c + c 1 x + c x + c 3 x 3 k = 1,,3,4 Vaues o these coeicients can be determined by appying geometric boundary conditions on d and θ at node 1 and. Let us deine N T (x) = [ N 1 (x) N (x) N 3 (x) N 4 (x) ] Substituting in (16) we obtain ρn T (x) q(t) 4 N T (x)q(t) t + EM I x 4 = (x,t) where q T = [d 1,θ 1,d,θ ] T. Soving this equation resuts in which is equivaent to ρn(x)n T (x)dx }{{} M q(t) t + N(x) N T (x) EM I x x q(t) }{{} = K N(x)(x,t)dx + F (17) } {{ } F d M q + Kq = F d + F F in terms o shear orce and moment m is given as 1 m F = 1 = m EM 3 d(x,t) x= I x 3 EM d(x,t) x= I x EM 3 d(x,t) x= I x 3 EM d(x,t) = I x = For the above interpoation unction the oowing oca mass and stiness matrices can be cacuated: 668

4 17th IFAC Word Congress (IFAC'8) Seou, Korea, Juy 6-11, 8 M = ρ K = EM I Generaization to cantiever with a arge number o eements The above expressions or K and M represent oca eements in FEM which are assembed in a goba mass and stiness matrix to generate umped matrices. The goba equation or the whoe assemby can be obtained by combining the matrix contribution o a individua eements, suchthatcoeicients beongingtocommon nodes are added together. Let us deine k (i) = [ k (i) 11 k(i) 1 k (i) 1 k(i) ] which is the stiness matrix o i th eement, then k (1) 11 k (1) 1 k (1) 1 k(1)... K = + k() 11 k (i 1) + k (i) 11 k (i) 1 k (i) 1 k (i) + k(i+1) 11 Simiary the goba mass matrix M can aso be ound. These goba mass and stiness matrices are used in (11) to generate a state space representation. In this paper a cantiever beam consisting o 1 nodes, as shown in Fig, is considered. Each node has our states (d, θ, d, θ), and the overa system has ive inputs and outputs at every second node aso shown in Fig. The inputs are the orces, and the outputs are the dispacements d o these eements. It is desired to convert the mode into an interconnected system with ive subsystems, each consisting o two nodes (shown by dashed boxes in Fig ). Thus we have xi = 8 and zi, wi = 1 i = 1,...,5. The numerica vaues o the physica parameters used or the simuation are isted in Tabe 1. By appying the too presented here to this umped system, an interconnected system is created. Both umped and interconnected systems are simuated. The resuting dispacement and ange o the ast node are shown in Fig 3 and 4. It is observed that the error between umped and interconnected system remain beow the numerica precision, which shows that the decomposition retains a... w 1 w z 1 z Fig.. Beam cantiever Tabe 1. Physica parameters o the cantiever exampe Fied w 3 z 3 w 4 z 4 Vaues Area o crossection m Length 1 m Length o eement.1 m Moduus o Easticity N m the dynamics o the system. Note that each subsystem has ony 8 states, as compared to a tota o 4 states. Dispacement [m] Times [sec] Fig. 3. Dispacement (d) o the 1 th node Anguar dispacement [rad] Times [sec] Fig. 4. Anguar dispacement (θ) o the 1 th node 5. CONTROLLER SYNTHESIS The active vibration contro o acantiever can be considered as an input disturbance rejection probem as shown w 5 z 5 669

5 17th IFAC Word Congress (IFAC'8) Seou, Korea, Juy 6-11, 8 in Figure 5, where, W 1 and W are the weighting iters. Let the state space o the weighted generaized pant be given as, ẋ = Ax + Bw + Bu z = C 1 x + D 1 u (18) y = C x + D 1 w where, d,z 1,z are the disturbance,inputs and outputs acting on every node. In the cantiever exampe we have considered state eedback case and d incude both orce and torque disturbances acting on each node. Using the d K(s) u G(s) Fig. 5. Generaized pant used or oopshaping. approach presented in section the umped system (18) is decomposed, where dynamics o the i th subsystem can be written as ẋ i (t) = A ii x i (t) + B 1i w i (t) + B i u i (t) + z i (t) = C 1i x i (t) + D 1 u i (t) y i (t) = C i x i (t) + D 1i w i (t) y W1 W N j=1,j i z1 z A ij x j (t) i = 1,,...N (19) x i R ni, w i R ri, u i R mi, z i R (i+mi) and y i R i are the states, disturbance input, controed input, controed output and measured output o the i th subsystem. The matrices A ii, B 1i, B i, C 1i, C i, D 1i and D 1i are o appropriate dimensions ound by decomp(). A stricty proper i th output eedback controer is deined by ẋ ci = A ci x ci + B ci y i u i = C ci x ci () i = 1,,...,N From (19) and () the cosed oop system rom w to z is given by where ẋ c = A c x c + B c w [ ] A B C A c = c B c C A c z = C c x c (1) B c = [ B1 B c D 1 ] C c = [C 1 D 1 C c ] () Decentraized contro o the system (18) requires that A c = diag{a c1,...a cn }, B c = diag{b c1,...b cn }, C c = diag{c c1,...,c cn }. In order to obtain a decentraized contro aw one needs to impose structura constraints on the bounded rea emma. In (Chen et a., 5) it has been proposed that or a given constant γ >, the system (18) is stabiizabe with the disturbance attenuation eve γ by a decentraized controer (), i there exist positive deinite bockdiagonamatrices X,Y andbock-diagona matrices F,L,Q such that T <, where T = and J 11 J1 T B 1 XC 1 J 1 J Y B 1 + LD 1 C1 T B1 T B1 T Y + D1L T T γi C 1 X + D 1 F C 1 γi [ ] X I > I Y (3) Positive deinite bock diagona matrices X and Y and diagona matrices F, L, Q are obtained by soving the above biinear matrix inequaity (BMI). Bock diagona coeicient matrices are then given by A c = V 1 QU T, B c = V 1 L, C c = FU T (4) where J 11 = AX + XA T + B F + F T B T, J 1 = A T + Y AX + LC X + Y B F + Q, J = Y A + A T Y + LC + C T L T. The size o sub-matrices in the bock diagona are compatibe with the dimensions o the state, input, and output vectors o the subsystems. The approach presented in (Chen et a., 5) invoves ixing a group o variabes to convert it into an LMI. First the variabes L, and Y and then the variabes X, F are ixed. Soving these two LMIs aternatey is then a way o inding a soution. In the proposed approach the path rom centraized controer to decentraized controer is divided into M steps and the structura constraints are graduay appied, by deining a rea number λ such that λ = k M, where k is graduay increased rom to M. This deines a matrix unction H as H(X,Y,F,L,Q,λ) = T(X,Y,F,L,(1 λ)q F + λq) < This is same as (3) except that J 1 is repaced by A T + Y AX + LC X + Y B F + (1 λ)q F + λq. Then { T(X,Y,F,L,QF ), λ = H(X,Y,F,L,Q,λ) = T(X,Y,F,L,Q), λ = 1 By using the approach o (Chen et a., 5) a state eedback decentraized controer or the active vibration contro o exibe beam is designed or W 1 = 1 and W =.1. The requency domain cosed-oop response between theorce input, ineardispacement andtheorce disturbance or the ast node are as shown in Figure 6 and 7. These igures show that the open oop system wi exhibit an osciatory response whereas in the cosed oop system the disturbance is rejected, whie the contro gain remains in reasonabe imits. One point which must be kept in mind is that or simpicity these requency domain resuts are generated rom the state space representation without imposing the appropriate boundary conditions. 67

6 17th IFAC Word Congress (IFAC'8) Seou, Korea, Juy 6-11, 8 Singuar Vaues (db) Singuar Vaues 1 Frequency 1 1 (rad/sec) Fig. 6. Frequencyresponse o the open-oop system (upper curve) and the cosed-oop system (ower curve) Singuar Vaues (db) N. Chen, M. Ikeda, and W. Gui. Design o robust H contro or interconnected systems: A homotopy method. Int. J. o Contro, Automation and Systems, 3(), 5. R. D Andrea and G. E. Duerud. Distributed contro design or spatiay interconnected systems. IEEE Transactions on Automatic Contro, 48(9): , 3. K. Huebner, D. Dewhirst, D. Smith, and T. Byrom. The Finite Eement Method or Engineers. A Wiey- Interscience Pubication, 1. ISBN J. Juang and M. Phan. Identiication and Contro o Mechanica Systems. Cambridge Universitypress, Cambridge, 4. ISBN D. Sijak. Large Scae Dynamic Systems - Stabiity and Structure. The Netherands: North Hoand, Amsterdam, D. Sijak and A. Zecevic. Contro o arge-scae systems: Beyond decnetraized eedback. Annua Reviews in Contro, 9(): , 5. G. Zhai, M. Ikeda, and Y. Fujisaki. Decentraized H controer design: A matrix inequaity approach using a homotopy method. Automatica, 37(4), Frequency 1 4 (rad/sec) Fig. 7. Contro sensitivity 6. CONCLUSION The main idea presented in this paper is that instead o reducing the mode order o a high order arge system one can decompose the overa system into sma subsystems, where each o these subsystems interacts with its neighboring eements. This approach is suitabe or FEM generated modes since these modes are oten o very high order and are generated rom a combination o sma eements. A simpe yet eicient too is presented here to carry out such a decomposition. This aows the appication o anaysis and synthesis toos or arge interconnected systems, and enabes the designer to synthesize controers or such systems in a systematic way. As an iustrative exampe, the too presented in this paper is appied to a cantiever mode. Resuts shows that the time responses o the origina system and the one decomposed into interconnected orm give identica resuts, up to numerica precision. Furthermore a previousy presented approach o designingdecentraized controer is appied to the decomposed system to demonstrate that by using the decomposition one can systematicay design a distributed contro strategy or FEM-generated systems without reducing the mode order. REFERENCES A. Becker. An introductory Guide to inite eement anaysis. Proessiona Engineering Pubishing, 4. ISBN

PERFORMANCE COMPARISON OF OPEN AND CLOSED LOOP OPERATION OF UPFC

PERFORMANCE COMPARISON OF OPEN AND CLOSED LOOP OPERATION OF UPFC OL. 3, NO. 5, OCTOE 008 ISSN 89-6608 AN Journa o Engineering and Appied Sciences 006-008 Asian esearch ubishing Networ (AN. A rights reserved. www.arpnournas.com EFOMANCE COMAISON OF OEN AND CLOSED LOO

More information

A Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem

A Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

Interactive Fuzzy Programming for Two-level Nonlinear Integer Programming Problems through Genetic Algorithms

Interactive Fuzzy Programming for Two-level Nonlinear Integer Programming Problems through Genetic Algorithms Md. Abu Kaam Azad et a./asia Paciic Management Review (5) (), 7-77 Interactive Fuzzy Programming or Two-eve Noninear Integer Programming Probems through Genetic Agorithms Abstract Md. Abu Kaam Azad a,*,

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Modal analysis of a multi-blade system undergoing rotational motion

Modal analysis of a multi-blade system undergoing rotational motion Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

COMPENSATION FOR DISTORTION OF DYNAMIC SENSOR DATA USING RECEPTANCE COUPLING

COMPENSATION FOR DISTORTION OF DYNAMIC SENSOR DATA USING RECEPTANCE COUPLING COMPENSATION FO DISTOTION OF DYNAMIC SENSO DATA SIN ECEPTANCE COPIN Tony. Schmitz niversity o Forida, Department o Mechanica and Aerospace Engineering, ainesvie, F 36. INTODCTION Sensors can be used to

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Mode in Output Participation Factors for Linear Systems

Mode in Output Participation Factors for Linear Systems 2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong

More information

Available online at ScienceDirect. IFAC PapersOnLine 50-1 (2017)

Available online at   ScienceDirect. IFAC PapersOnLine 50-1 (2017) Avaiabe onine at www.sciencedirect.com ScienceDirect IFAC PapersOnLine 50-1 (2017 3412 3417 Stabiization of discrete-time switched inear systems: Lyapunov-Metzer inequaities versus S-procedure characterizations

More information

SECTION A. Question 1

SECTION A. Question 1 SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second

More information

Recursive Asymptotic Hybrid Matrix Method for Acoustic Waves in Multilayered Piezoelectric Media

Recursive Asymptotic Hybrid Matrix Method for Acoustic Waves in Multilayered Piezoelectric Media Open Journa o Acoustics, 2,, 27-33 doi:.4236/oja.2.4 Pubished Onine September 2 (http://.scirp.org/journa/oja) Recursive Asymptotic Hybrid Matrix Method or Acoustic Waves in Mutiayered Piezoeectric Media

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

Fitting affine and orthogonal transformations between two sets of points

Fitting affine and orthogonal transformations between two sets of points Mathematica Communications 9(2004), 27-34 27 Fitting affine and orthogona transformations between two sets of points Hemuth Späth Abstract. Let two point sets P and Q be given in R n. We determine a transation

More information

Categories and Subject Descriptors B.7.2 [Integrated Circuits]: Design Aids Verification. General Terms Algorithms

Categories and Subject Descriptors B.7.2 [Integrated Circuits]: Design Aids Verification. General Terms Algorithms 5. oward Eicient arge-scae Perormance odeing o Integrated Circuits via uti-ode/uti-corner Sparse Regression Wangyang Zhang entor Graphics Corporation Ridder Park Drive San Jose, CA 953 wangyan@ece.cmu.edu

More information

A New Backpressure Algorithm for Joint Rate Control and Routing with Vanishing Utility Optimality Gaps and Finite Queue Lengths

A New Backpressure Algorithm for Joint Rate Control and Routing with Vanishing Utility Optimality Gaps and Finite Queue Lengths A New Backpressure Agorithm or Joint Rate Contro and Routing with Vanishing Utiity Optimaity Gaps and Finite Queue Lengths Hao Yu and Michae J. Neey Abstract The backpressure agorithm has been widey used

More information

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Simulation of single bubble rising in liquid using front tracking method

Simulation of single bubble rising in liquid using front tracking method Advances in Fuid Mechanics VI 79 Simuation o singe bubbe rising in iquid using ront tracking method J. Hua & J. Lou Institute o High Perormance Computing, #01-01 The Capricorn, Singapore Abstract Front

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Finite element method for structural dynamic and stability analyses

Finite element method for structural dynamic and stability analyses Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

Strain Energy in Linear Elastic Solids

Strain Energy in Linear Elastic Solids Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua

More information

SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT

SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT 8 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT P. Wang, N. Hamia *, P. Boisse Universite de Lyon, INSA-Lyon,

More information

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS

More information

CE601-Structura Anaysis I UNIT-IV SOPE-DEFECTION METHOD 1. What are the assumptions made in sope-defection method? (i) Between each pair of the supports the beam section is constant. (ii) The joint in

More information

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

ISSN: [Alkhatib* et al., 5(10): October, 2016] Impact Factor: 4.116

ISSN: [Alkhatib* et al., 5(10): October, 2016] Impact Factor: 4.116 [Akhatib* et a., 5(10): October, 2016] Impact Factor: 4.116 IC Vaue: 3.00 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY VIBRATION ANALYSIS, SHAPE AND DESIGN OF MOTORCYCLE MOUNTING

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

Algorithms to solve massively under-defined systems of multivariate quadratic equations

Algorithms to solve massively under-defined systems of multivariate quadratic equations Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

Leader-Follower Consensus Modeling Representative Democracy

Leader-Follower Consensus Modeling Representative Democracy Proceedings o the Internationa Conerence o Contro, Dynamic Systems, and Robotics Ottawa, Ontario, Canada, May 7-8, 2015 Paper o. 156 Leader-Foower Consensus Modeing Representative Democracy Subhradeep

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

More information

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2 MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis

More information

A nodal collocation approximation for the multidimensional P L equations. 3D applications.

A nodal collocation approximation for the multidimensional P L equations. 3D applications. XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

arxiv: v1 [cs.sy] 31 May 2014

arxiv: v1 [cs.sy] 31 May 2014 Cooperative Contro of Linear Muti-Agent Systems via Distributed Output Reguation and Transient Synchronization Georg S. Seyboth, Wei Ren, Fran Agöwer arxiv:46.85v cs.sy 3 May 24 Abstract A wide range of

More information

UNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS

UNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS 11th Word Congress on Computationa Mechanics WCCM XI 5th European Conference on Computationa Mechanics ECCM V 6th European Conference on Computationa Fuid Dynamics ECFD VI E. Oñate J. Oiver and. Huerta

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence

More information

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe

More information

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR 1 Maximizing Sum Rate and Minimizing MSE on Mutiuser Downink: Optimaity, Fast Agorithms and Equivaence via Max-min SIR Chee Wei Tan 1,2, Mung Chiang 2 and R. Srikant 3 1 Caifornia Institute of Technoogy,

More information

VTU-NPTEL-NMEICT Project

VTU-NPTEL-NMEICT Project MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid

More information

Input-to-state stability for a class of Lurie systems

Input-to-state stability for a class of Lurie systems Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems

More information

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water

More information

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com XI PHYSICS M. Affan Khan LECTURER PHYSICS, AKHSS, K affan_414@ive.com https://promotephysics.wordpress.com [TORQUE, ANGULAR MOMENTUM & EQUILIBRIUM] CHAPTER NO. 5 Okay here we are going to discuss Rotationa

More information

Substructuring Preconditioners for the Bidomain Extracellular Potential Problem

Substructuring Preconditioners for the Bidomain Extracellular Potential Problem Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di

More information

Optimum Design Method of Viscous Dampers in Building Frames Using Calibration Model

Optimum Design Method of Viscous Dampers in Building Frames Using Calibration Model The 4 th Word Conference on Earthquake Engineering October -7, 8, Beijing, China Optimum Design Method of iscous Dampers in Buiding Frames sing Caibration Mode M. Yamakawa, Y. Nagano, Y. ee 3, K. etani

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Instructional Objectives:

Instructional Objectives: Instructiona Objectives: At te end of tis esson, te students soud be abe to understand: Ways in wic eccentric oads appear in a weded joint. Genera procedure of designing a weded joint for eccentric oading.

More information

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

3.10 Implications of Redundancy

3.10 Implications of Redundancy 118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

Distributed average consensus: Beyond the realm of linearity

Distributed average consensus: Beyond the realm of linearity Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,

More information

HILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT

HILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

Distributed Cross-Layer Optimization of Wireless Sensor Networks: A Game Theoretic Approach

Distributed Cross-Layer Optimization of Wireless Sensor Networks: A Game Theoretic Approach Distributed Cross-Layer Optimization o Wireess Sensor Networks: A Game Theoretic Approach Jun Yuan and Wei Yu Eectrica and Computer Engineering Department, University o Toronto {steveyuan, weiyu}@commutorontoca

More information

Adjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models

Adjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models IO Conference Series: Earth and Environmenta Science AER OEN ACCESS Adjustment of automatic contro systems of production faciities at coa processing pants using mutivariant physico- mathematica modes To

More information

Path planning with PH G2 splines in R2

Path planning with PH G2 splines in R2 Path panning with PH G2 spines in R2 Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru To cite this version: Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru. Path panning with PH G2 spines

More information

SPRINGS - Functions. Spring Classifications. Spring Performance Characteristics DEFINITION FUNCTIONS C R I T E R I A. Spring rate (spring constant) k

SPRINGS - Functions. Spring Classifications. Spring Performance Characteristics DEFINITION FUNCTIONS C R I T E R I A. Spring rate (spring constant) k SPRINGS - unctions DEINITION machine parts that are designed and constructed to give a reativey arge eastic deection when oaded UNCTIONS To appy oad and to contro motion (brakes & cutches, cam oowers,

More information

VI.G Exact free energy of the Square Lattice Ising model

VI.G Exact free energy of the Square Lattice Ising model VI.G Exact free energy of the Square Lattice Ising mode As indicated in eq.(vi.35), the Ising partition function is reated to a sum S, over coections of paths on the attice. The aowed graphs for a square

More information

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

Technical Data for Profiles. Groove position, external dimensions and modular dimensions Technica Data for Profies Extruded Profie Symbo A Mg Si 0.5 F 25 Materia number.206.72 Status: artificiay aged Mechanica vaues (appy ony in pressing direction) Tensie strength Rm min. 245 N/mm 2 Yied point

More information

Proceedings of Meetings on Acoustics

Proceedings of Meetings on Acoustics Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.

More information

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY

More information

Volume 13, MAIN ARTICLES

Volume 13, MAIN ARTICLES Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi

More information

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

The influence of temperature of photovoltaic modules on performance of solar power plant

The influence of temperature of photovoltaic modules on performance of solar power plant IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance

More information

CALIBRATION OF RIVER BED ROUGHNESS

CALIBRATION OF RIVER BED ROUGHNESS CALIBRATION OF RIVER BED ROUGHNESS A. M. Wasantha La 1, M. ASCE Singuar vaue decomposition (SVD) is used to caibrate the Manning s roughness coefficients in a 1-D unsteady fow mode of the Upper Niagara

More information

Identification of macro and micro parameters in solidification model

Identification of macro and micro parameters in solidification model BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7 6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the

More information