CONTROL BY INTERCONNECTION OF THE TIMOSHENKO BEAM. DEIS, University of Bologna, viale Risorgimento 2, Bologna (Italy)

Size: px
Start display at page:

Download "CONTROL BY INTERCONNECTION OF THE TIMOSHENKO BEAM. DEIS, University of Bologna, viale Risorgimento 2, Bologna (Italy)"

Transcription

1 CONTROL BY INTERCONNECTION OF THE TIMOSHENKO BEAM Alessandro Macchelli Claudio Melchiorri EIS, Universiy of Bologna, viale Risorgimeno 2, 4036 Bologna Ialy {amacchelli, Asrac: In his paper, he dynamical conrol of a mixed finie and infinie dimensional mechanical sysem is approached wihin he framework of por Hamilonian sysems. As an applicaive example of he presened mehodology, a flexile eam, modeled according o he Timoshenko heory, wih a mass under graviy field conneced o a free end, is considered. Afer he disriued por Hamilonian dph model of he eam is inroduced, he conrol prolem is discussed. The concep of srucural invarian Casimir funcion is generalized o he infinie dimensional case and he so-called conrol y inerconnecion conrol echnique is exended o he infinie dimensional case. In his way, finie dimensional passive conrollers can sailize disriued parameer sysems y shaping heir oal energy, i.e. y assigning a new minimum in he desired equilirium configuraion ha can e reached if a dissipaion effec is inroduced. Keywords: disriued por Hamilonian sysems, energy shaping, Timoshenko eam, conrol. INTROUCTION Flexile eams are generally modeled according o he classical Euler-Bernoulli heory: his formulaion provides a good descripion of he dynamical ehavior of he sysem if he eam s cross secional dimension are small in comparison of is lengh. In his case, he effecs of he roary ineria of he eam are no considered. A more accurae eam model is provided y he Timoshenko heory, according o which he roary ineria and also he deformaion due o shear are considered. The resuling Timoshenko model of he eam is generally more accurae in predicing he eam s response han he Euler-Bernoulli one, u, on he oher hand, i is more difficul o uilize for conrol purposes ecause of is complexiy. The dph van der Schaf and Maschke, 2002 formulaion of he Timoshenko eam Macchelli and Melchiorri, 2003 provides a eer undersanding of he model from a physical poin of view, oviously wihou reducing he complexiy of he model iself. I makes clear how he kineic and he poenial elasic energy domains inerac and how he sysem can exchange power wih he environmen hrough is order and/or a disriued por. Moreover, from he conrol poin of view, he dph represenaion of he sysem allows o easily exend well-esalished passive conrol sraegies, ha were originally developed for finie dimensional por Hamilonian sysems, o deal wih disriued sysems. This paper is organized as follows: in Sec. 2, he dph model of he Timoshenko eam is presened y inroducing he underlying irac srucure. Then, he mixed finie and infinie dimensional por Hamilonian model of he plan is presened in Sec. 3 and, in Sec. 4, some consideraions on he srucural invarians Casimir funcions of he closed-loop sysem are developed. This is he

2 saring poin in he definiion of he conroller, ha is widely discussed in Sec. 5. Finally, some concluding remarks and some suggesions for fuure work are illusraed in Sec PH MOEL OF THE TIMOSHENKO BEAM 2. Background. The classical formulaion According o he Timoshenko heory, he moion of a eam can e descried y he following sysem of PE: ρ 2 w 2 I ρ 2 φ 2 EI 2 φ x 2 K K 2 w x 2 K φ x = 0 φ w = 0 x where is he ime and x 0, L] is he spaial coordinae along he eam in is equilirium posiion. Then, wx, is he deflecion of he eam from he equilirium configuraion, denoed y w 0, and φx, is he roaion of he eam s cross secion due o ending. Assume ha he moion akes place in he wx-plane. The coefficiens ρ, I ρ, E and I are he mass per uni lengh, he mass momen of ineria of he cross secion, Young s modulus and he momen of ineria of he cross secion, respecively. The parameer K is equal o kga, where G is he modulus of elasiciy in shear, A is he cross secional area and k is a consan depending on he shape of he cross secion. The oal energy is given y he following relaion kineic energy {}}{ L 2 ] 2 w φ H := ρ I ρ dx 2 0 L K φ w 2 ] 2 φ EI dx 2 0 x x }{{} poenial elasic energy 2 ha poins ou ha, as every infinie dimensional sysem, he eam is characerized y a spaial domain := 0, L], wih order = {0, L}, and y he presence of wo ineracive energy domains, he kineic and he poenial elasic. 2.2 Timoshenko eam in dph form The saring poin in he definiion of a por Hamilonian sysem oh finie and infinie dimensional is he idenificaion of a suiale space of power or energy variales, sricly relaed o he geomery of he sysem, and he definiion of a irac srucure on his space of power variales, in order o descrie he inernal and exernal inerconnecion of he sysem. The poenial elasic energy in 2 is a funcion of he shear and of he ending, given y he following -forms: ɛ, x = ] w, x φ, x dx x ɛ r, x = φ, xdx x 3 The associaed co-energy variales are he shear force and he ending momenum, given y he 0-forms funcions σ, x = K ɛ, x and σ r, x = EI ɛ r, x, where is he Hodge sar operaor, Araham and Marsden, 978. Besides, he kineic energy is funcion of he ranslaional and roaional momena, i.e. of he following -form: p, x = ρ w, xdx φ p r, x = I ρ, xdx 4 and he associaed co-energy variales are he ranslaional and roaional momena, given y he 0-forms funcions v, x = ρ p, x and v r, x = I ρ p r, x. If N is an n-dimensional Riemannian manifold, he space of k-forms on N, i.e. he space of k- linear alernaing funcions, is given y Ω k N. So, we have ha p, p r, ɛ, ɛ r Ω and ha w, φ Ω 0. The oal energy 2 ecomes he following quadraic funcional: Hp, p r, ɛ, ɛ r = Hp, p r, ɛ, ɛ r 5 wih H : Ω Ω Ω, Hp, p r, ɛ, ɛ r = 2 ρ p p I ρ p r p r K ɛ ɛ EI ɛ r ɛ r he energy densiy. Consider a ime funcion p, p r, ɛ, ɛ r Ω Ω wih R, and evaluae he energy H along his rajecory. A any ime, he variaion of inernal energy, ha is he power exchanged wih he environmen, is given y dh d = δ p H p δ p r H p r δ ɛ H ɛ δ ɛ r H ɛ r 6 The differenial forms p, pr, ɛ and ɛr are he ime derivaives of he energy variales p, p r, ɛ and ɛ r, and represen he generalized velociies flows, while δ p H, δ pr H, δ ɛ H, δ ɛr H are he variaional derivaive of he oal energy H : Ω Ω R. They are relaed o he

3 rae of change of he sored energy and represen he generalized forces effors. The dph formulaion of he Timoshenko eam can e oained or y expressing in erms of p, p r, ɛ and ɛ r 3 and 4, or, in a more rigorous way, y revealing he underlying irac srucure of he model. Firs of all, i is necessary o define he space of power variales. The space of flows is given y F := Ω Ω Ω Ω Ω 0 Ω 0, and he space of effor E, he dual of F, is given y E := Ω 0 Ω 0 Ω 0 Ω 0 Ω 0 Ω 0, since he concep of dualiy over he space of forms can e given y he following: Proposiion. Assume ha N is an n-dimensional manifold. Then, he dual space Ω k N of Ω k N can e idenified wih Ω n k N and he dualiy produc eween Ω k N and Ω k N y β, α := β α 7 wih α Ω k N and β Ω n k N. The same resul holds for Ω k N. Proof. See van der Schaf and Maschke, Given he dualiy produc defined in 7, i is possile o inroduce he pairing operaor on he space of power variales y means of he following definiion, ha can e easily specialized for he Timoshenko eam case. efiniion 2. Suppose ha α i, β i Ω k N Ω n k N, i =, 2. Then α, β, α 2, β 2 := β 2, α β, α 2 = β 2 α β α 2 N Since he space of power variales and he pairing operaor are specified, i is possile o give he following: efiniion 3. Suppose ha F and E are linear spaces wih a dual pairing,. A irac srucure is a linear suspace F E such ha =, wih denoing he orhogonal complemen wih respec o he pairing,. The previous definiion is quie general and can e reformulaed as follows. Given f i, e i F E, i =, 2, hen f, e, f 2, e 2 = 0. I is immediae ha, if f, e, hen N 0 = f, e, f, e = 2 e, f 8 In oher words, if f, e F E is a couple of power conjugaed variales, he fac ha hey elong o he irac srucure implies power conservaion, i.e. he dual produc is equal o 0. The irac srucure is he geomerical ool y means of which i is possile o deal wih power conserving inerconnecion in physical sysems. Once a proper inerconnecion srucure is defined, he por Hamilonian model of a physical sysem, oh finie and infinie dimensional, follows auomaically. Wih he following proposiion, he main resul of his secion is presened. Proposiion 4. Consider he space of power variales F E and he ilinear form pairing operaor, inroduced y ef. 2; hen, define he following linear suspace of F E: = { f p,..., f r, e p,..., e r F E f p 0 0 d 0 e p f pr f ɛ = 0 0 d e pr d 0 0 e ɛ, f ɛr 0 d 0 0 e ɛr f f r e e r = e p e pr e ɛ e ɛr } 9 where denoes he resricion on he order of he spaial domain. Then =, or, equivalenly, is a irac srucure. Proof. As in van der Schaf and Maschke, 2002, he proof can e divided in wo seps. In he firs one, i is verified ha, while, in. Suppose ha ω i, ei p,..., e r,i F E, wih i =, 2; hen, if ω, ω 2, i happens ha ω, ω 2 = 0. From he definiion 9 of he irac srucure, we have he second one, ha = fp i,..., f r,i ω, ω 2 = de 2 p e ɛ e 2 p de ɛ de p e 2 ɛ e p de 2 ɛ de 2 p r e ɛ r e 2 p r de ɛ r de p r e 2 ɛ r e p r de 2 ɛ r f, e,2 f,2 e, f r, e r,2 f r,2 e r, Since de e 2 = de e 2 e de 2, and, for he Sokes heorem, if α Ω 0, hen dα = α, we deduce ha ω, ω 2 = 0 and. In order o prove ha, consider ω 2. From he definiion of irac srucure we have ha ω, ω, ω 2 = 0. Since ω, from 9 we have:

4 0 = ω, ω 2 = de ɛ e 2 p fp 2 e p e ɛ de ɛ r e 2 p r fp 2 r e ] p r de p e p r e 2 ɛ fɛ 2 e ] ɛ de p r e 2 ɛ r fɛ 2 r e ɛ r e p e,2 f,2 e ɛ e p r e r,2 f r,2 e ɛ r From he Sokes Theorem and he properies of he exerior derivaive, is is possile o oain ha 0 = e ɛ de 2 p e 2 p r fɛ 2 ] e p fp 2 de 2 ɛ e ɛ r de 2 p r fɛ 2 r ] e pr fp 2 r e 2 ɛ de 2 ɛ r ] ] e ɛ e 2 p f,2 ] e ɛ r e 2 p r f r,2 ] e p e 2 ɛ e,2 ] e p r e 2 ɛ r e r,2 for every ω. Since he exerior produc of forms is non degeneraive, we deduce ha he previous relaion holds if and only if ω 2. So, and his complees he proof. Consider he oal energy 5 as he Hamilonian of he sysem, i.e. a quadraic funcional of he energy variales p, p r, ɛ and ɛ r ounded from elow. The rae of change of hese energy variale generalized velociies can e conneced o he irac srucure 9 y seing f p = p, f ɛ = ɛ, f p r = pr and f ɛr = ɛr, where he minus sign is necessary in order o have a consisen energy flow descripion. Moreover, he rae of change of he Hamilonian wih respec o he energy variales, ha is is variaional derivaives, can e relaed o he irac srucure y seing e p = δ p H, e ɛ = δ ɛ H, e pr = δ pr H and e ɛr = δ ɛr H. Then, he disriued Hamilonian formulaion wih oundary energy flow of he Timoshenko eam can e given y he following: efiniion 5. The dph model of he Timoshenko eam wih irac srucure 9 and Hamilonian H 5 is given y p 0 0 d 0 δ p H p r ɛ = 0 0 d δ pr H d 0 0 δ ɛ H ɛ r 0 d 0 0 δ ɛr H f δ p H f r e = δ pr H δ ɛ H e r δ ɛr H 0 Noe ha 0 coincides, afer some manipulaions, wih. The power conserving propery 8 implies ha for every f p,..., e r e p f p e pr f pr e ɛ f ɛ e ɛr f ɛr e f e r f r = 0 Then, from 6, he following proposiion can e easily proved Proposiion 6. Consider he dph model of he Timoshenko eam 0. Then dh d = e f e r f r = e, xf, x e r, xf r, x ] x=l x=0 or, in oher words, he increase of energy kineic/poenial energy of he eam is equal o he power supplied hrough he order. I is imporan o noe ha he Timoshenko eam has een already reaed wihin he framework of dph sysems also in Golo e al., This approach differs from he one presened in his paper in he mahemaical descripion of he space of energy/power variales L 2 funcions and differenial forms and in he definiion of dualiy over he space of power variales, Golo, Moreover, in he proposed approach, he Sokes heorem is he key poin in undersanding he relaion eween power variales on he order and evoluion of he energy variales inside he domain. 3. MOEL OF THE PLANT Consider he mechanical sysem of Fig., in which a flexile eam is conneced o a rigid ody wih mass m and ineria momenum J in x = L and o a conroller in x = 0. The conroller acs on he sysem wih a force f c and a orque τ c. Since he Timoshenko model of he eam is valid only for small deformaions, i is possile o assume ha he moion of he mass is he cominaion of a roaional and of a ranslaional moion along x = L. The por Hamilonian model of he mass

5 PSfrag replacemens { f L f r f L ] T = e f c τ c 0 L q 2 m, J Fig.. Flexile link wih mass in x = L. Plan overview. is given y: q = p H ṗ = q H p H f e = p H q 2 where q = q, q 2 ] T Q are he generalized coordinaes, wih q he disance from he equilirium configuraion and q 2 he roaion angle, p T Q are he generalized momena and Hp, q := 2 pt M p mgq = p 2 2 m p2 2 J mgq is he oal energy Hamilonian, wih g he graviy acceleraion. Moreover, he por variales are given y f, e R 2. As regard he conroller, we assume ha i can e modeled y means of he following finie dimensional por Hamilonian sysems q c = pc H c ṗ c = qc H c c pc H c G c f c 3 e c = G T c pc H c where q c Q c are he generalized coordinaes, wih dimq c = 2, p c T Q c are he generalized momena and f c, e c R 2 are he power conjugaed por variales. Moreover, H c q c, p c is he Hamilonian and i will e specified in he remaining par of his secion in order o drive he whole sysem in a desired equilirium configuraion The por causaliy of oh he mass and he conroller is assumed o e wih flows as inpu and effors as oupus. As poined ou in Sramigioli, 200, i is possile o inerconnec wo por Hamilonian sysems only if a por dualizaion is applied on one of he sysem. In his way, a sysem can have an effor as inpu and a flow as oupu. The ond graph represenaion of he closed loop sysem made of he Timoshenko eam, he mass in x = L and he finie dimensional rag replacemens por Hamilonian conroller acing in x = 0 is given in Fig. 2, where SGY is he symplecic gyraor ha implemens he por dualizaion. Then, he inerconnecions consrains eween H c SGY H SGY Fig. 2. Bond graph represenaion of he closedloop sysem. he por variales of he susysems are given y he following power-preserving relaions: H e L e r L ] T = f { f 0 ff r 0 ] 4 T = ec e 0 e r 0 ] T = fc From 0, 2, 3 and 4 i is possile o oain he mixed finie and infinie dimensional por Hamilonian represenaion of he closedloop sysem. The oal energy H cl is defined in he exended space χ := q, p, q c, p c, p, p r, ɛ r, ɛ ] T and i is given y he sum of he energy funcions of he susysems, ha is H cl := H H c H 5 Moreover, i is easy o verify ha he energy rae is equal o dh cl T H cl = d p H cl p T H cl H cl c p c p c where c and H c have o e designed in order o drive he sysem o he desired equilirium posiion, in general differen from he rivial one, for which only some damping injecion, i.e. c > 0, is enough. The asic idea is o shape he oal energy H cl y properly choosing he conroller Hamilonian H c in order o have a new minimum of energy in he desired configuraion ha can e reached if a dissipaive effec is inroduced. 4. STRUCTURAL INVARIANTS IN THE CLOSE-LOOP SYSTEM In general, he energy variales of he conroller are no relaed o he energy variales of he plan and, as a consequence, i is no clear how he oal energy can e shaped and he desired minimum assigned. The firs sep is o relae he energy variales of he conroller o he energy variales of he plan independenly from he Hamilonian of he conroller and, if i is possile, from he Hamilonian of he plan. These invarians are called, in he finie dimensional framework, Casimir funcions, van der Schaf, 999; Orega e al., 999. The Casimir funcions are srucural invarians in he sense ha hey are invarian independenly from he Hamilonian of he sysem. In he case of mixed finie and infinie dimensional sysem, i is possile o give he following: efiniion 7. Consider a mixed finie and infinie dimensional por Hamilonian sysem wih sae space given y X X, where X is he sae space of he finie dimensional par and X he sae space of he infinie dimensional one, and Hamilonian H : X X R. A funcional C : X X R is a Casimir funcional for he sysem if and only if dc d = 0, H Hamilonian of he sysem

6 ha is C is consan independenly of he Hamilonian of he sysem. For he sysem of Fig., we have ha X = Q T Q Q c T Q c and X = Ω Ω Ω Ω, wih Hamilonian given y 5. Given χ := q, p, q c, p c, p, p r, ɛ r, ɛ ] T X X and a funcional C : X X R, we have ha dc C d = T q q T C p ṗ T C q c T C ṗ c q c p c δ p C p δ p r C p r δ ɛ C ɛ δ ɛ r C ɛ r ha has o e equal o zero for every Hamilonian H, H c and H. I is possile o prove, ha his is rue if and only if dδ ɛ C = 0 dδ p C δ pr C = 0 dδ ɛr C δ ɛ C = 0 dδ pr C = 0 C p = C ] ] δp C = L δp C = 0 = 0 p c δ pr C L δ pr C 0 ] ] C q = δɛ C L C δɛ C = G 0 δ ɛr C c L q c δ ɛr C 0 In oher words, he following proposiion holds. 6 Proposiion 8. Consider he mixed finie and infinie dimensional por Hamilonian sysem of Fig. 2, ha is he resul of he power conserving inerconnecion 4 of he susysems 0, 2 and 3. If X X is he exended sae space of he sysem, hen a funcional C : X X R is a Casimir for he closed-loop sysem if and only if relaions 6 hold. In order o conrol he flexile eam wih he finie dimensional conroller 3, he firs sep is o find Casimir funcionals for he closed-loop sysem ha can relae he sae variales of he conroller q c o he sae variales ha descrie he configuraion of he flexile eam and of he mass conneced o is exremiy. In paricular, we are looking for some funcionals C i, i =, 2, such ha C i q, p, q c, p c, p, p r, ɛ, ɛ r := q c,i C i q, p, p c, p, p r, ɛ, ɛ r, wih i =, 2, are Casimir funcionals for he closed loop sysem, ha is hey saisfies he condiions of Prop. 8. Firs of all, from 6, i is immediae o noe ha every Casimir funcional canno depend on p and p c. Moreover, since i is necessary ha dδ ɛ C i = 0 and dδ pr C i = 0, we deduce ha δ ɛ C i and δ pr C i have o e consan as funcion on x on and heir value will e deermined y he oundary condiions on C i. Since, from 6, δ pr C i = 0, we deduce ha δ pr C i = 0 on. Bu dδ p C i = δ pr C i = 0, hen, from he oundary condiions, we deduce ha also δ p C i = 0 on. As a consequence, all he admissile Casimir funcionals are also independen from p and p r, ha is C i q, q c, ɛ, ɛ r := q c,i C i q, ɛ, ɛ r, wih i =, 2. Assuming G c = I, we have ha ] ] C δɛ C = = x=0 q c 0 δ ɛr C x=0 7 and, consequenly, δ ɛ C = on. From 6, we have ha dδ ɛr C = δ ɛ C = = dx; hen, δ ɛr C = x c, where c is deermined y he oundary condiions. Since, from 7, δ ɛr C x=0 = 0, hen c = 0; moreover, we deduce ha δ ɛr C x=l = L, relaion ha inroduces a new oundary condiion in x = L. A consequence is ha ] ] C q = δɛ C x=l = δ ɛr C x=l L The firs conclusion is ha C q,..., ɛ r = q c, Lq 2 q xɛ r ɛ 8 is a Casimir for he closed loop sysem. Following he same procedure, i is possile o calculae C 2. From 6, we have ha C 2 q c = 0 ] ] δɛ C = 2 x=0 δ ɛr C 2 x=0 9 and hen δ ɛ C 2 = 0 on ; moreover, dδ ɛr C 2 = 0 and, consequenly, δ ɛr C 2 = on since 9 holds. Again from 6, we deduce ha ] ] δɛ C 2 x=l C 2 q = δ ɛr C 2 x=l So we can sae ha C 2 q, q c, ɛ, ɛ r = q c,2 q 2 = 0 ɛ r 20 is anoher Casimir funcional for he closed loop sysem. In conclusion, he following proposiion has een proved. Proposiion 9. Consider he mixed finie and infinie dimensional por Hamilonian sysem of Fig. 2, ha is he resul of he power conserving inerconnecion 4 of he susysems 0, 2 and 3. Then 8 and 20 are Casimir funcionals for his sysem. Since C i, i =, 2, are Casimir funcionals, hey are invarian for he sysem of Fig. 2. Then, for every energy funcion H c of he conroller, we have ha q c, = Lq 2 q xɛ r ɛ C 2 q c,2 = q 2 ɛ r C 2 where C and C 2 depend on he iniial condiions. If he iniial configuraion of he sysem is known,

7 hen i is possile o assume hese consans equal o zero. Since H c is an arirary funcion of q c, i is possile o shape he oal energy funcion of he closed-loop sysem in order o have a minimum of energy in a desired configuraion: if some dissipaion effec is presen, he new equilirium configuraion will e reached. 5. CONTROL BY ENERGY SHAPING In his secion, he sailizaion of he mechanical sysem of Fig. y energy shaping is presened. The desired equilirium configuraion expressed in erms of energy variales can e calculaed as he soluion of 0 wih f,r p = p r = ɛ = ɛ r = 0 0,L = 0, e r 0,L = e 0 = 0, e L = mg ha leads o he following equilirium soluion: p, x = 0 p r, x = 0 ɛ, x = mg ɛ r, x = mg K EI x L 22 If α and q are he desired roaional angle and posiion of he mass, ha is φ, L = α and w, L = q, he desired configuraion 22 of he eam can e easily expressed in erms of is verical displacemen and cross secion roaion: w, x = mg x L3 6EI α mg x L q 23 K φ, x = mg 2EI x L2 α The energy funcion H c of he conroller 3 will e developed in order o regulae he closed-loop sysem in he configuraion χ = q, p, p c, p, p r, ɛ, ɛ r In he remaining par of his secion i will proved ha, y choosing he conroller energy as H c q c,p c := 2 pt c Mc p c 2 K c,q c, q c, 2 2 K c,2q c,2 q c,2 2 Ψ q c, q c, Ψ 2 q c,2 q c,2 24 wih M c = Mc T > 0, K c,, K c,2 > 0, qc, and qc,2 given y 2 evaluaed in he equilirium configuraion, Ψ and Ψ 2 o e properly chosen, he configuraion χ previously inroduced is gloally asympoically sale. As in he case of finie dimensional Hamilonian sysem, he sailiy of a mixed finie and infinie dimensional Hamilonian sysem can e proved if i can e shown ha he equilirium is a sric exremum of he oal energy of he closed-loop sysem. The only difference is ha, in order o prove he sailiy for he infinie dimensional par, i is necessary o fix a norm: i is imporan o noe ha he sailiy wih respec o his given norm will no assure he sailiy wih respec o a differen one. The sailiy definiion in he sense of Lyapunov for mixed finie and infinie dimensional sysem can e given y he following: efiniion 0. The equilirium configuraion χ for a mixed finie and infinie dimensional sysem is said o e sale in he sense of Lyapunov wih respec o he norm if for every ɛ > 0 here exiss δ ɛ > 0 such ha χ0 χ < δ ɛ χ χ < ɛ for all > 0, where χ0 is he iniial configuraion of he sysem, Swaers, As proposed in Rodriguez e al., 200; Swaers, 2000, in order o verify he sailiy of χ, i is necessary o show ha χ is an exremum of he closed-loop energy funcion H cl inroduced in 22, wih H c given y 24, ha is he condiion H cl χ = 0 has o hold. Moreover, if χ is he displacemen from he equilirium configuraion χ, inroduce he non linear funcional N χ := H cl χ χ H cl χ 25 ha is proporional o he second variaion of H cl. Then, he configuraion χ is sale if i is possile o find γ, Γ > 0 such ha Swaers, 2000: 0 N χ Γ χ γ 26 In his case, he norm can e convenienly chosen as he Euclidean norm on he exended space of he closed-loop sysem. I easy o prove ha H cl χ = 0 if and only if Ψ = mg and Ψ 2 = mgl in 24. Moreover, since H cl χ = 0, he linear funcional inroduced in 25 is given y N χ = 2 pt M p 2 pt c Mc p c H p, p r, ɛ, ɛ r 2 K c, L q 2 q 2 2 K c,2 q 2 ɛ r 2 x ɛ r ɛ Since N χ > 0 if χ 0, he sailiy proof is compleed if proper consans γ, Γ saisfying 26 are found. If Γ = 2 max{ M, Mc, K c,, K c,2, 2 ρ, 2 I ρ, 2K, 2EI}, hen

8 N χ Γ p 2 p c 2 χ 2] Γ L q 2 q x ɛ r ɛ 2 Γ q 2 ɛ r 2 If γ = 2 and Γ = Γ max{2, 2L 2 2L }, hen condiion 26 is saisfied and he configuraion χ is sale. Summarizing, he following proposiion has een proved. Proposiion. Consider he mixed finie and infinie dimensional por Hamilonian sysem of Fig. 2, ha is he resul of he power conserving inerconnecion 4 of he susysems 0, 2 and 3. If in 3 i is assumed ha G c = I and H c is chosen according o 24, hen he configuraion descried in 22, 23 is sale in he sense of Lyapunov. 6. CONCLUSIONS AN FUTURE WORK Saring from he dph formulaion of he Timoshenko eam, i has een shown ha i is possile o exend he conrol y inerconnecion and energy shaping o rea mixed finie and infinie dimensional sysems. In paricular, he conrol of a mechanical sysem made of a flexile eam wih a mass conneced a one of is exremiy under graviy field has een presened. The finie dimensional conroller, acing on he sysem hrough he oher exremiy, is developed y properly exending he concep of Casimir funcions o he infinie dimensional case. Fuure work will deal wih he exension of hese conceps o he modeling and conrol of simple kinemaic chains wih flexile links. Acknowledgmen. This work has een suppored y he European projec GeoPlex, IST Furher informaions a hp:// Orega, R., A.J. van der Schaf, B. Maschke and G. Escoar 999. Energy-shaping of porconrolled hamilonian sysems y inerconnecion. In: Proc. of he 38h CC Conference. Rodriguez, H., A.J. van der Schaf and R. Orega 200. On sailizaion of non-linear disriued parameer por-conrolled hamilonian sysems via energy shaping. In: Proc. of he 40h CC Conference. Sramigioli, S Modeling and IPC conrol of ineracive mechanical sysems: a coordinae free approach. Springer-Verlag. London, UK. Swaers, G.E Inroducion o Hamilonian fluid dynamics and sailiy heory. Chapman & Hall / CRC. van der Schaf, A.J L 2 -gain and Passiviy Techniques in Nonlinear Conrol. Springer- Verlag. London, UK. van der Schaf, A.J. and B.M. Maschke Hamilonian formulaion of disriued parameer sysems wih oundary energy flow. J. of Geom. and Ph. 7. REFERENCES Araham, R. and J.E. Marsden 978. Foundaion of Mechanics. Benjamin/Cummings Pulishing Company. Golo, G Inerconnecion srucures in por-ased modeling: ools for analysis and simulaion. Ph hesis. Universiy of Twene. Golo, G., V. Talasila and A.J. van der Schaf A hamilonian formulaion of he imoshenko eam model. In: Proc. of Mecharonics Macchelli, A. and C. Melchiorri isriued por hamilonian formulaion of he imoshenko eam: Modeling and conrol. In: Proc. of 4h MATHMO, Vienna.

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/

More information

Flow-Induced Vibration Analysis of Supported Pipes with a Crack

Flow-Induced Vibration Analysis of Supported Pipes with a Crack Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Phys 331: Ch 4. Conservative Forces & 1-D Systems 1. 1 W b

Phys 331: Ch 4. Conservative Forces & 1-D Systems 1. 1 W b Phys 33: Ch 4. Conservaive orces & -D Sysems Thurs. 9/3 ri., / Mon. /4 Wed. /6 Thurs. /7 ri. /8 4.4-.6 Curl of Conservaive orce, Varying Poenial, -D sysems 4.7-.8 Curvilinear -D, Cenral orce 4.9 Energy

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Finite Element Analysis of Structures

Finite Element Analysis of Structures KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Basilio Bona ROBOTICA 03CFIOR 1

Basilio Bona ROBOTICA 03CFIOR 1 Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Finite element method for structural dynamic and stability analyses

Finite element method for structural dynamic and stability analyses Finie elemen mehod for srucural dynamic and sabiliy analyses Module- Nonlinear FE Models Lecure-39 Toal and updaed Lagrangian formulaions Prof C Manohar Deparmen of Civil Engineering IIc, Bangalore 56

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Q.1 Define work and its unit?

Q.1 Define work and its unit? CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar

More information

Generalized Chebyshev polynomials

Generalized Chebyshev polynomials Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT

More information

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

RC, RL and RLC circuits

RC, RL and RLC circuits Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie e Quanum eory of Aoms and Molecules: e Scrodinger equaion Hilary erm 008 Dr Gran Ricie An equaion for maer waves? De Broglie posulaed a every paricles as an associaed wave of waveleng: / p Wave naure of

More information

Single and Double Pendulum Models

Single and Double Pendulum Models Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double

More information

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product 11.1 APPCATON OF AMPEE S AW N SYMMETC MAGNETC FEDS - f one knows ha a magneic field has a symmery, one may calculae he magniude of by use of Ampere s law: The inegral of scalar produc Closed _ pah * d

More information

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

4. Advanced Stability Theory

4. Advanced Stability Theory Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium

More information

Chapter 10 INDUCTANCE Recommended Problems:

Chapter 10 INDUCTANCE Recommended Problems: Chaper 0 NDUCTANCE Recommended Problems: 3,5,7,9,5,6,7,8,9,,,3,6,7,9,3,35,47,48,5,5,69, 7,7. Self nducance Consider he circui shown in he Figure. When he swich is closed, he curren, and so he magneic field,

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Kinematics and kinematic functions

Kinematics and kinematic functions Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Lab 10: RC, RL, and RLC Circuits

Lab 10: RC, RL, and RLC Circuits Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7]. Chaper Numerical Tess We here collec a few numerical ess, in order o pu ino evidence he poenialiies of HBVMs [4, 6, 7]. Tes problem Le us consider he problem characerized by he polynomial Hamilonian (4.)

More information

Acceleration. Part I. Uniformly Accelerated Motion: Kinematics & Geometry

Acceleration. Part I. Uniformly Accelerated Motion: Kinematics & Geometry Acceleraion Team: Par I. Uniformly Acceleraed Moion: Kinemaics & Geomery Acceleraion is he rae of change of velociy wih respec o ime: a dv/d. In his experimen, you will sudy a very imporan class of moion

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables

More information

Oscillations. Periodic Motion. Sinusoidal Motion. PHY oscillations - J. Hedberg

Oscillations. Periodic Motion. Sinusoidal Motion. PHY oscillations - J. Hedberg Oscillaions PHY 207 - oscillaions - J. Hedberg - 2017 1. Periodic Moion 2. Sinusoidal Moion 3. How do we ge his kind of moion? 4. Posiion - Velociy - cceleraion 5. spring wih vecors 6. he reference circle

More information

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

AP Calculus BC Chapter 10 Part 1 AP Exam Problems AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

Unsteady Mass- Transfer Models

Unsteady Mass- Transfer Models See T&K Chaper 9 Unseady Mass- Transfer Models ChEn 6603 Wednesday, April 4, Ouline Conex for he discussion Soluion for ransien binary diffusion wih consan c, N. Soluion for mulicomponen diffusion wih

More information

Rhythm of the Heart. Megan Levine. Abstract

Rhythm of the Heart. Megan Levine. Abstract Rhyhm of he Hear Megan Levine Asrac An overlook of he hear and how i works as a pump wih an eplanaion of how o uild a mahemaical model using differenial equaions. The Hear. The Hear Bea The hear ea is

More information

Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced

Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced Lecure 1 Inroducion Why is i imporan o sudy waves in plasma? Waves are naurally found in plasmas and have o be deal wih. This includes insabiliies, flucuaions, waveinduced ranspor... Waves can deliver

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Fall 200 Dep. of Chemical and Biological Engineering Korea Universiy CHE302 Process Dynamics and Conrol Korea Universiy

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

IMPACT OF AN OBLIQUE BREAKING WAVE ON A WALL

IMPACT OF AN OBLIQUE BREAKING WAVE ON A WALL Source: Physics of Fluids Vol 6 No pp 6-64 4 DOI: 6/64445 IMPACT OF AN OIQUE REAKING WAVE ON A WA Jian-Jun SHU School of Mechanical & Aerospace Engineering Nanyang Technological Universiy 5 Nanyang Avenue

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN

More information

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3 The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped

More information

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Failure of the work-hamiltonian connection for free energy calculations. Abstract

Failure of the work-hamiltonian connection for free energy calculations. Abstract Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,

More information