Finite Element Simulation of Smart Lightweight Structures for Active Vibration and Interior Acoustic Control

Size: px
Start display at page:

Download "Finite Element Simulation of Smart Lightweight Structures for Active Vibration and Interior Acoustic Control"

Transcription

1 ECHNISCHE MECHANIK, Bnd 3, Heft 1, (3), Mnuskripteingng: 3. Februr 3 Finite Element Simultion of Smrt Lightweight Strutures for Ative Vibrtion nd Interior Aousti Control J. Lefèvre, U. Gbbert Dedited to Dr.-Ing. Friedrih Whl on the osion of his 65 th birthdy. he pper presents numeril pproh to tive vibrtion nd noise ontrol of smrt lightweight strutures. he struture is provided with thin piezoeletri wfers s tutors nd sensors to ontrol vibrtions of the struture. Fully oupled eletromehnil field equtions were tken into ount where model bsed ontrollers re pplied for design purposes. he objetive of vibrtion ontrol of elsti strutures is to redue interior noise levels. Hene, the mehnil field is lso oupled with the ousti field, nd onsequently fully oupled eletro-mehnil-oustil problem needs to be solved. he numeril solution is bsed on the finite element method, introduing veloity potentil for the ousti fluid to reeive overll symmetri system mtries of the semi-disrete form of the eqution of motion. It is shown tht the vibro-ousti oupling n be negleted for ontroller design purposes, nd onsequently the modl truntion tehnique onsidering only the unoupled struturl modes, n be dpted to vibro-ousti systems. he behviour of smrt plte struture oupled with n ousti vity is studied s referene exmple. 1 Introdution Over the pst few yers smrt struturl onepts hve met with growing interest in mny engineering brnhes nd severl new pprohes nd tehnologies hve been developed (zou&gurn, 1998; Gbbert, ). Smrt strutures, or strutroni systems, re hrterized by synergisti integrtion of tive mterils into pssive struture onneted by ontrol system to filitte utomti djustment to hnging environmentl onditions. Piezoeletri mterils re widely used s distributed sensors nd tutors in smrt strutures. Commerilly vilble piezoeletri wfers re ommonly used s tive mterils for ontrolling struture vibrtions. Suh thin wfers my be glued on the surfe of the bse struture, or embedded in omposite mteril during the mnufturing proess. Intensified tivities in developing nd pplying piezoeletri smrt strutures require effetive nd relible design tools. he finite element method (FEM) is n exellent bsis to develop powerful softwre tools. Being widely spred, it hs beome theoretilly nd prtilly estblished method for solving oupled piezoeletri field problems. Bsed on the generl purpose finite element softwre COSAR (COSAR, 199), tool hs been designed nd grdully implemented over the pst few yers (Berger et l., ; Gbbert et l., ). he softwre pkge ontins n extended librry of multi-field finite elements for 1D, D nd 3D ontinu s well s for shell-type thin-wlled strutures (Seeger et l. ; Gbbert et l., b). he COSAR FEM softwre is linked with ontroller design tools, suh s MtLb/Simulink, through generl dt interfe, designing the ontroller on the bsis of finite element models (Gbbert et l. 1; Gbbert et l., ). Sound rdition is nother mjor problem in designing engineering strutures. Whenever n elsti struture gets into ontt with surrounding fluid, struturl vibrtions nd the ousti pressure field will influene eh other. his vibro-ousti oupling results in dditionl pressure lods on the struture-fluid interfe used by fluid pressure, s well s unwnted noise rdition used by struturl vibrtions. Strong vibro-ousti oupling effets our in thin-wlled lightweight strutures of lrge surfe. Ative noise ontrol by pplying piezoeletri pth tutors to the struture is n lterntive wy to redue noise rdition of elsti strutures (see Blhndrn et l., 1996; Ro&Bz, 1999; Kim et l., 1999; Gopinthn et l.,, nd others). Only reently, we extended our softwre tool COSAR by inluding new brik-type finite elements for disretizing ousti volumes (Lefèvre, ). he veloity potentil of the fluid ws onsidered s n dditionl nodl degree of freedom in order to reeive symmetri system mtries. Also, the vibro-ousti oupling effet ws tken into ount. Now, it is possible to solve the fully oupled eletro-mehnil-ousti field problem on the bsis of finite element disretiztion of the eletromehnil struture nd the fluid. hermopiezoeletri effets s disussed by Görnndt et l.,, re negleted. Further below, the theoretil bsis of this finite 59

2 element pproh is presented, followed by disussion of the ontroller design for vibro-ousti problems. As lrge-sle finite element models nnot be used for ontroller design purposes, n pproprite model redution is required to redue the number of the degrees of freedom. Here, modl truntion seems to be the best suited method for estblishing the ontroller design of strutures bsed on the finite element disretiztion s flexible strutures exhibit low-pss hrteristis, whih llow to neglet high-frequeny dynmis (Strßberger, 1997). Bsed on seleted eigenmodes the finite element model is redued nd trnsformed into the stte spe form. All required dt of the modl stte spe model n be trnsferred to the ontroller design tool (MtLb/Simulink) vi generl dt exhnge interfe to design suitble model-bsed ontroller. o study the ontrolled behviour of the struture, the ontroller mtries n be retrnsferred to the originl finite element model. As shown in the pper, dequte results n be reeived for ontroller design purposes by tking into ount only the eigenmodes of the unoupled struture. his proedure hs been pplied to interior ousti problems with strong fluid-struture intertion. he simultion of smrt retngulr elsti plte struture oupled with n ousti vity is investigted s referene exmple. Governing Equtions nd Finite Element Anlysis his setion presents the theoretil bkground of our finite element softwre tool for simulting smrt strutures oupled with n ousti fluid. First, the finite element nlysis of oupled eletromehnil strutures nd ousti fluids re presented seprtely, followed by demonstrtion how both pprohes re oupled to reeive n overll vibro-ousti finite element model. he following equtions re bsed on Crtesin (x 1,x,x 3 )-oordinte system..1 Finite Element Model of Piezoeletri Smrt Strutures he oupled eletromehnil behviour of polrizble piezoeletri mteril in low voltge pplitions n be modelled with suffiient ury by mens of linerized onstitutive equtions. Furthermore, smll displements re onsidered. he omplete derivtion of the finite element model is disussed in detil in Gbbert et l.(). Hene, only brief summry of the equtions is presented here. In three-dimensionl ontinuum the finite element equtions re bsed on the mehnil equilibrium D p = ρ u, (1) u the eletri equilibrium (4 th Mxwell eqution) D ϕ D = () nd the liner oupled eletromehnil onstitutive equtions = C ee, (3) D = e (, (4) with the stress vetor = [ σ11 σ σ 33 σ1 σ 3 σ 31], the body fore vetor p = [ p1 p p3] displement vetor u = [ u1 u u3] = [ D1 D D3] [3 3], the strin vetor = [ ε11 ε ε33 ε1 ε 3 ε 31] E = [ E E ]. D u nd D ϕ desribe the following differentil mtrix opertors, the, the mss density ρ, the vetor of eletril displements D nd the elstiity mtrix C[6 6], the piezoeletri mtrix e [6 3], the dieletri mtrix 1 E3 nd the eletri field vetor nd x1 x x3 D u = (5) x x1 x3 x3 x x1 6

3 D ϕ =. (6) x1 x x3 ogether with stress boundry onditions nd the hrge boundry onditions, the mehnil nd the eletri blne equtions n be written in wek form s δu ( D ) dv ( u p ρu δϕ Dϕ )dv χ = D V δ ( q q) δϕ( Q Q) do = V O q u do (7) where δu is virtul displement, δϕ is virtul eletri potentil, q is given trtion vetor on the surfe O q nd Q is given hrge on the surfe O Q. In eqution (7) the strin-displement reltion = D u u (8) nd the reltion between the eletri field nd the eletri potentil E = D ϕ ϕ, (9) re introdued. he displements u 1, u, u 3 nd the eletri potentil ϕ re pproximted elementwise by shpe funtions, ontining the element nodl prmeters s unknown prmeters. Inluding the shpe funtions of finite element in the mtries G u nd G ϕ nd inorporting the unknown nodl prmeters in the vetors w nd, the pproximtion of the displement vetor u nd the eletri potentil ϕ n be written s ( x x, x ) G ( x, x x ) w u 1 3 u 1, 3, =, (1) ( x, x ) ( x, x x ) ϕ x1, 3 = Gϕ 1, 3. (11) Following the stndrd finite element proedure, the semi-disrete form of the equtions of motion of generl element (e) n be derived from eqution (7) s follows M uu w Kuu Kuϕ w fu () () () () =. (1) e e e e Kuϕ Kϕϕ fϕ he vetor w ontins the nodl displements nd the vetor ontins the nodl eletri potentils; M uu is the element mss mtrix, K uu is the element stiffness mtrix, K ϕϕ is the element eletri mtrix, K uϕ is the piezoeletri oupling mtrix, f u is the lod vetor resulting from mehnil lods nd f ϕ is the lod vetor resulting from eletri hrges.. Finite Element Model of the Aousti Fluid Aousti responses in fluid re usully regrded s smll perturbtions relted to n mbient referene stte. Hene, the finite element model of the fluid is derived from the liner ousti wve eqution (Fhy, 1994) 1 p = p, (13) onsidering the dibti wve propgtion in homogeneous, invisid fluid. In eqution (13) is the Lplin opertor, p is the ousti pressure nd is the sound speed. Although we re minly interested in sound pressure distribution, the use of pressure p s nodl vrible is rther disdvntgeous s in this se the finite element formultion of the oupled vibro-ousti problem results in un-symmetri system mtries (Desmet&Vndepitte, 1999). Everstine (1997) reommended the introdution of the veloity potentil Φ of the fluid s new degree of freedom in order to obtin symmetri mtries. he veloity potentil is slr field relted to the veloity v of the fluid prtiles by O Q 61

4 v = D Φ (14) nd to the sound pressure by p = ρ Φ, (15) with the differentil mtrix opertor D, hving the sme form s the opertor in eqution (6), nd the fluid density ρ (Kollmnn, ). Introduing equtions (14) nd (15) into eqution (13) results in the new form of the ousti wve eqution 1 Φ = Φ. (16) he next step onsiders n interior ousti problem, i.e. the whole fluid volume is enlosed by boundries in ll diretions with finite distne from ommon referene point. o obtin the wek form of the ousti wve eqution, eqution (16) is multiplied by ny test funtion τ nd integrted over the fluid volume V i. Using the Gussin integrl theorem the wve eqution is derived in wek form s 1 τ D D Φ dv τ Φ dv = ( τ D Φ) n do Vi Vi Oi, (17) where n is the outwrd direted norml vetor of the surfe O i of the volume V i. he right hnd side of eqution (17) orresponds to the boundry onditions. O i n be divided into the surfes O v with given norml veloity v n nd the surfes O z with n imposed impedne funtion Z n. he impedne funtion desribes the reltionship between the ousti pressure nd the norml omponent of the veloity p Z n =. (18) v n Following this proedure, it is possible to model boundry dmping effets (Kim et l., 1999) nd with the equtions (14) nd (15) eqution (17) n be written s 1 ρ D τ Φ do. (19) Z τ DΦ dv τ Φ dv = τ vn do Vi Vi O n v Oz In finite element the veloity potentil Φ nd the test funtion τ re pproximted by using the vetors of the nodl degrees of freedom, nd the mtrix G ontining the pproximte funtions s nd ( ) ( ), x, x = x, x x Φ x G 3 () 1 3 1, (, x, x ) = ( x, x x ) τ x G 3. (1) 1 3 1, With the bbrevition B = D G, () nd the equtions () nd (1) the eqution (19) n be written s 1 G G dv ρ Z () () () () G e e e e G G do B B dv = n Vi Oz Vi Ov v do. (3) n As this eqution must to be fulfilled for ny prmeter of the test funtion from eqution (3) we reeive M C K = f, (4) 6

5 with the ousti element mss mtrix 1 M = V i the ousti element dmping mtrix ρ = Z C G n O z the ousti element stiffness mtrix G G dv, (5) G do, (6) K = B B dv (7) V i nd the ousti element lod vetor f = G v do. (8) O v.3 Vibro-Aousti Coupling n First, it is ssumed tht the ousti pressure represents n dditionl lod with respet to the struture. Developing our finite element model, we must lso onsider n dditionl lod vetor f () u e resulting from the ousti pressure. Consequently, eqution (1) n be ugmented to M uu w K uu Kuϕ w fu fu () () () () =. (9) e e e e Kuϕ Kϕϕ fϕ With eqution (15) the vetor f n be lulted s u fu = Gu npdo = Gu ng do = C ρ u, (3) O O S S where O s is the fluid struture interfe nd C u is the oupling mtrix with regrd to the ousti pressure. Furthermore the veloity of the vibrting struture ts on O s s new ousti lod of the fluid M C K = f f. (31) Considering tht the norml veloity of the struture n be expressed by u = n G w, (3) n u we reeive n expression similr to eqution (8) f = G n G do w = C w, (33) O S u with C s the oupling mtrix regrding struturl vibrtions. Compring eqution (3) nd eqution (33), symmetri system mtries re reeived by multiplying ll lines relted to fluid degrees of freedom by (-ρ ). herefore nd when the element mtries re inluded in globl mtrix nd the oupling mtries re inserted 63

6 in the left hnd side, the semi-disrete system of equtions of the eletro-mehnil-ousti field problem n be written s M uu w Cuu Cu w M ρ C C u ρ Kuu Kuϕ w fu = Kuϕ Kϕϕ fϕ ρ K ρ f. (34) In eqution (34) the Ryleigh dmping hs been introdued s the struturl dmping C uu. Mention should be mde of the ft tht the vetor w lso represents the nodl displements of the pssive prt of the struture. Bsed on eqution (34) two ousti brik-type elements with 8 nd nodes of the Serendipity element fmily hve been developed nd implemented in our COSAR softwre tool. he isoprmetri element onept ws used to properly pproximte the element geometry. Furthermore, lso the vibro-ousti oupling proedures s presented in equtions (9)-(3) were inorported in the COSAR softwre tool. ogether with initil onditions eqution (34) n be integrted numerilly, e.g. by using the Newmrk formuls. 3 Control of Smrt Strutures Numeril simultions of smrt strutures under the finite element onept require n overll finite element model omprising the pssive struture, the tive sensor nd the tutor elements s well s suitble model of the ontroller. ody, omprehensive design tools suh s MtLb/Simulink re vilble to support the designing proess. A generl dt exhnge interfe is required to exhnge dt nd informtion between the finite element model nd the ontroller design tool. Only reently, suh dt interfe ws developed to ouple our finite element COSAR softwre with MtLb/Simulink (see Gbbert et l., 1, ). 3.1 Controller Design For ontroller design purposes the vibro-ousti oupling effet is negleted in the tive struture. Combining the displement nd the eletri potentil degrees of freedom in one vetor w = [ w ] eqution (1) n be written s, (35) () t Bu () t M w Cw Kw = Ef, (36) r where C is n dditionl dmping mtrix nd u r (t) is the vetor of the ontroller influene of the struture. he mtries E nd B desribe the positions of the externl fores nd the ontroller prmeters in the finite element model, respetively. o use the modl truntion tehnique the liner eigenvlue problem ( K M) wˆ = ω, (37) k k =, where m is the number of eigen- needs to be solved. he result is the modl mtrix Q [ wˆ wˆ... wˆ 1 m ] modes onsidered. Ortho-normlizing Q with Q MQ = I = dig() 1 nd Q KQ = = dig( ω ) Q CQ = = dig ( δ ω ) k k, onsidering proportionl dmping. Introduing modl oordintes qˆ s k, we lso obtin w = Qqˆ (38) into eqution (36), we obtin the modl trunted system of (m m) differentil equtions () t Q Bu () t q ˆ qˆ qˆ = Q Ef r. (39) 64

7 Eqution (39) is trnsformed into the stte spe form, whih is more onvenient for the ontrol theory. Introduing the stte spe vetor qˆ z = q ˆ (4) in eqution (39) results in I z = z ur() t f() t = Aqz() t Bqur() t Eqf() t. (41) Q B Q E In smrt strutures vriety of signls n be mesured, first of ll eletri potentils by using piezoeletri sensors; but lso other signls, suh s exiting fores or displements, n be mesured in order to hrterize the stte spe of the struture. Hene, the mesurement mtrix is used in generl form s () t C z D u () t F f() t y =, (4) q q r q where C q, D q nd F q re mpping mtries desribing the reltions between the mesured quntities of the vetors z, u r nd f with respet to the mesuring vetor y. he mtries A q, B q, E q, C q, D q, F q re trnsferred to MtLb/Simulink where the ontroller is designed. Designing time independent LQ ontroller, the tutor signl is generted by the ontroller mtrix R s u r () t Rz() t 3. Solution Conept =. (43) he semi-disrete form of the oupled vibro-ousti eqution of motion inluding ontrol re s follows M w C Cu M ρ Cu ρc K w w Ef = ρ K ρ () t Bur () t f () t. (44) C u is the oupling mtrix with respet to the vetor w. o solve eqution (44), the tutor signl (43) must be expressed in terms of w nd w, respetively. If only few seleted eigenmodes re onsidered, Q is not squre mtrix. o get the modl oordintes from eqution (38), the pseudo-inverse mtrix of Q, 1 ( Q Q) Q Q =, (45) is required. With eqution (45) the stte spe vetor hs the following form Q z = Q w. (46) w Mention should be mde tht the use of Q in eqution (46) does not provide n ext solution with respet to the modl oordintes. However, if the motion of the struture is dominted by the seleted eigenmodes, there is very good greement between the pproh in eqution (46) nd the ext solution s we ould estblish. Following the integrtion proedure of the initil finite element model, the tutor signl u r n be inluded s n dditionl fore vetor in the right hnd side of eqution (44). 4 Exmples his setion ontins numeril demonstrtion of the ft tht ontroller designed with modl redued system s bove n be pplied with suffiient ury to fully oupled vibro-ousti system. he finite element 65

8 model of lmped retngulr plte struture oupled with n ousti vity is used s test exmple (Fig. 1). Four piezoeletri pthes tthed to the plte s two olloted tutor/sensor pirs (tutors on the top nd sensors on the bottom of the plte) re employed. he struture is exited by hrmoni fore t given point. All dimensions nd mteril properties of the plte, the ousti vity nd the tutors/sensors re given in bles 1-4. he plte is pproximted by 96 pssive hexhedron elements ontining only mehnil degrees of freedom; the vity is disretized with 384 ousti hexhedron elements. Eh tutor or sensor is pproximted by oupled piezoeletri hexhedron element ontining mehnil s well s eletril degrees of freedom. In ll these finite elements qudrti shpe funtions re used. he Newmrk formuls re employed for numeril time integrtion of the equtions of motion. sensor/tutor pirs F(t) x h l 3 x 1 x 3 l 1 l Figure 1. Smrt Plte Struture Coupled with n Aousti Cvity. Plte/Cvity System Atutors/Sensors l 1 6mm Length in x 1 -diretion 1mm l 4mm Length in x -diretion 5mm l 3 4mm Pth thikness.mm h (Plte thikness) mm ble 1. Geometry of the Vibro-ousti System Inluding Atutors nd Sensors. Elsti Constnts Piezoeletri Constnts Dieletri Constnts N/mm² 33 14N/mm² e N/(mV)mm κ N/(mV)² 1 631N/mm² N/mm² e N/(mV)mm κ N/(mV)² N/mm² 66 N/mm² e N/(mV)mm Density ρ Ns²/mm 4 ble. Mteril Properties of the Piezoeletri Atutors nd Sensors. Young s modulus E 1N/mm² Poisson s rtio ν.3 Density ρ P Ns²/mm 4 ble 3. Mteril Properties of the Elsti Plte. Speed of sound 34mm/s Fluid density ρ Ns²/mm 4 ble 4. Mteril Properties of the Aousti Fluid. Number Frequeny Hz 9.9Hz Hz Hz Hz ble 5. Eigenfrequenies of the Unoupled Elsti Plte. 66

9 In ll simultions fore exittion with n mplitude of 1N nd frequeny ontining the first eigenfrequeny (see b. 5) ws pplied, whih is very importnt exittion regrding the influene to the ousti pressure in the enlosure. In the simultions the ontroller ws swithed on fter period of 1s. A LQ ontroller for the elsti plte ws designed on the bsis of redued model ontining the first five struturl modes (see b. 5). o verify the reeived ontroller mtrix, the response of the plte ws simulted without ny vibro-ousti oupling. Fig. displys the resulting signl of the sensor. 8.E3 6.E3 4.E3.E3.E M-.E3-4.E3-6.E3-8.E sensor signl t [s] Figure. Sensor Signl of the Unoupled Plte Resulting from Hrmoni Exittion. [mv] As n be seen in Fig., the designed ontroller redues struturl vibrtions to pproximtely qurter of its mximum mplitude. he seond simultion dditionlly onsidered the vibro-ousti oupling. For this simultion the pressure relese ondition (p=) ws pplied to ll boundries of the fluid in ddition to oupling with n elsti plte, s inluded in the oupling equtions (3) nd (33). Fig. 3 displys the sensor signl. [mv] 8.E3 6.E3 4.E3.E3.E M-.E3-4.E3-6.E3-8.E sensor signl t [s] Figure 3. Resulting Sensor Signl with Respet to the Coupled Vibro-ousti System under the Assumption of Boundry Pressure Relese Conditions. A omprison between Fig. nd Fig. 3 revels tht sound pressure does not exert strong influene on the elsti plte when pressure relese onditions re ssumed t the boundries. For the lst simultion we ssumed tht the pressure relese ondition only ts on the boundry x 3 =l 3. Aprt from the plte ll other four boundries of the vity were onsidered s oustilly hrd surfes (v n =). Fig. 4 displys the signl reeived by the sensors. Fig. 5 shows the orresponding sound pressure in the entre of the vity. Fig. 4 nd Fig. 5 disply strong oupling effet. he oupled fluid serves s n dditionl dmping with respet to the vibrting plte. Furthermore, the motion of the plte is no longer only dominted by the first unoupled eigenmode. he reson is the ousti fluid ltering the eigenfrequenies nd eigenmodes of the whole system. [mv] 8.E3 6.E3 4.E3.E3.E -.E3-4.E3-6.E3-8.E sensor signl t [s] M Figure 4. Resulting Sensor Signl of the Coupled Vibro-ousti System Considering Aoustilly Hrd Boundries. 67

10 p [N/mm²] 5.E-6 4.E-6 3.E-6.E-6 1.E-6.E -1.E-6 -.E-6-3.E-6-4.E-6-5.E t [s] sound pressure Figure 5. Resulting Sound Pressure in the Cvity Centre when Considering Aoustilly Hrd Boundries. As n be seen in ll simultions, the ontroller redues the sensor signls to lmost the sme level. In our opinion the fluid-struture intertion n be negleted for ontroller design purposes if there is no strong oupling between the struture nd the fluid. However, we often enounter problems where strong vibro-ousti oupling ours. Bsilly, the intertion of fluid nd struture needs to be onsidered for ontroller design whenever lrge surfe re of thin-wlled lightweight struture is overed by fluid. his is lso of gret importne when suh models re used for lulting optiml tutor nd sensor positions t the struture. In suh ses the modl truntion should be bsed on the fully oupled eletro-mehnil-oustil eigenvlue problem. 5 Conlusions he pper presents the theoretil bkground of new finite element-bsed softwre tool for solving threedimensionl eletro-mehnil-oustil field problems, inluding ontrol lgorithms. his tool ws inorported in our finite element nlysis softwre COSAR, filitting simultions of the ontrolled behviour of oupled vibro-ousti systems with distributed piezoeletri tutors nd sensors with regrd to the interior ousti rdition problem. A generl dt exhnge interfe is used to design ontrollers by MtLb/Simulink. Furthermore, numeril investigtions were performed whether the modl truntion tehnique using the unoupled struturl modes for ontroller design n be dpted to vibro-ousti systems. If there is no strong oupling between the fluid nd the struture this method yields suffiient results. Otherwise, the eigenmodes of the ompletely oupled system need to be onsidered. Aknowledgement his work hs been supported by the postgrdute progrm of the federl stte of Shsen-Anhlt. his support is grtefully knowledged. Referenes Blhndrn, B., Smprth, A., Prk, J.: Ative ontrol of interior noise in three-dimensionl enlosure, J. of Smrt Mterils nd Strutures, Vol. 5, (1996), pp Berger, H., Gbbert, U., Köppe, H., Seeger, F.: Finite Element Anlysis nd Design of Piezoeletri Controlled Smrt Strutures, J. of heoretil nd Applied Mehnis, 3, 38, (), pp COSAR Generl Purpose Finite Element Pkge: Mnul, FEMCOS mbh Mgdeburg (see lso (199). Desmet, W., Vndepitte, D.: Finite Element Method in Aoustis, Pro. of the 1 th Interntionl Seminr on Advned ehniques in Applied nd Numeril Aoustis, Leuven, Belgium, (1999). Everstine, G.C.: Finite Element Formultions of Struturl Aousti Problems, J. Computers & Strutures, Vol. 65, No. 3, (1997), pp

11 Fhy, F.: Sound nd Struturl Vibrtion. Rdition, rnsmission nd Response, Ademi Press, London, (1994). Gbbert, U., Berger, H., Köppe, H., Co, X.: On Modelling nd Anlysis of Piezoeletri Smrt Strutures by the Finite Element Method, J. of Applied Mehnis nd Engineering, Vol. 5, No. 1, (), pp Gbbert, U., Köppe, H., Seeger, F.: Overll design of tively ontrolled smrt strutures by the finite element method, SPIE Proeedings Series, Vol. 436, (1), pp *DEEHUW8 SSH1HVWRURYLü7UDMNRY7(QWZXUILQWHOOLJHQWHU6WUukturen unter Einbeziehung der Regelung. t Automtisierungstehnik 5, 9, (), pp Gbbert, U., Köppe, H., Seeger, F., Berger, H.: Modeling of smrt omposite shell strutures, J. of heoretil nd Applied Mehnis, 3, 4, (b), pp Gbbert, U.: Reserh tivities in smrt mterils nd strutures nd expettions to future developments, J. of heoretil nd Applied Mehnis, 3, 4, (), pp Görnndt, A., Gbbert, U. (): Finite Element Anlysis of hermopiezoeletri Smrt Strutures. At Mehni, 154, (), pp Gopinthn, S. V., Vrdn, V. V., Vrdn, V. K.: Finite Element/Boundry Element Simultion of Interior Noise Control Using Ative-Pssive Control ehnique, SPIE Proeedings Series, Vol. 3984, (), pp. -3. Kim, J., Ko, B., Lee, J.-K., Cheong, C.-C.: Finite element modeling of piezoeletri smrt struture for the bin noise problem, J. of Smrt Mterils nd Strutures, Vol. 8, (1999), pp Kim, J., Lee, J.-K., Im, B.-S., Chung, C.-C.: Modeling of Piezoeletri Smrt Strutures Inluding Absorbing Mterils for Cbin Noise Problems, SPIE Proeedings Series, Vol. 3667, (1999), pp Kollmnn, F. G.: Mshinenkustik. Grundlgen, Meßtehnik, Berehnung, Beeinflussung, Springer-Verlg, Berlin-Heidelberg, (). Lefèvre, J.: Modellierung und Simultion geregelter Leihtbustrukturen zur ktiven Reduktion der Shllbstrhlung, Diplomrbeit, Otto-von-Guerike-Universität Mgdeburg, Institut für Mehnik, (). Ro, J., Bz, A.: Control of sound rdition from plte into n ousti vity using tive onstrined lyer dmping, J. of Smrt Mterils nd Strutures, Vol. 8, (1999), pp Seeger, F., Gbbert, U., Köppe, H., Fuhs, K.: Anlysis nd Design of hin-wlled Smrt Strutures in Industril Applitions, SPIE Proeedings Series, Vol. 4698, (), pp Strßberger, M.: Aktive Shllreduktion durh digitle Zustndsregelung der Strukturshwingungen mit Hilfe piezo-kermisher Aktoren, Mitteilungen us dem Institut für Mehnik, Ruhr-Universität Bohum, (1997). zou, H.-S., Gurn, A. (Eds.): Strutroni Systems: Smrt Strutures, Devies nd Systems, Prt 1: Mterils nd Strutures, Prt : Systems nd Control, World Sientifi, (1998). Address: Dipl.-Ing. Jen Lefèvre, Prof. Dr.-Ing. hbil. Dr. h.. Ulrih Gbbert, Otto-von-Guerike-Universität Mgdeburg, Universitätspltz, D-3916 Mgdeburg, e-mil: jen.lefevre@mb.uni-mgdeburg.de, ulrih.gbbert@mb.uni-mgdeburg.de 69

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction ECHNISCHE MECHANIK, Band 25, Heft 3-4, (25), 241-247 Manuskripteingang: 11. Dezember 25 Finite Element Modelling of Vibro-Aousti Systems for Ative Noise Redution J. Lefèvre, U. Gabbert Over the past years

More information

Electromagnetism Notes, NYU Spring 2018

Electromagnetism Notes, NYU Spring 2018 Eletromgnetism Notes, NYU Spring 208 April 2, 208 Ation formultion of EM. Free field desription Let us first onsider the free EM field, i.e. in the bsene of ny hrges or urrents. To tret this s mehnil system

More information

Lecture Notes No. 10

Lecture Notes No. 10 2.6 System Identifition, Estimtion, nd Lerning Leture otes o. Mrh 3, 26 6 Model Struture of Liner ime Invrint Systems 6. Model Struture In representing dynmil system, the first step is to find n pproprite

More information

Lecture Summaries for Multivariable Integral Calculus M52B

Lecture Summaries for Multivariable Integral Calculus M52B These leture summries my lso be viewed online by liking the L ion t the top right of ny leture sreen. Leture Summries for Multivrible Integrl Clulus M52B Chpter nd setion numbers refer to the 6th edition.

More information

Dorf, R.C., Wan, Z. T- Equivalent Networks The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

Dorf, R.C., Wan, Z. T- Equivalent Networks The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000 orf, R.C., Wn,. T- Equivlent Networks The Eletril Engineering Hndook Ed. Rihrd C. orf Bo Rton: CRC Press LLC, 000 9 T P Equivlent Networks hen Wn University of Cliforni, vis Rihrd C. orf University of

More information

ANALYSIS AND MODELLING OF RAINFALL EVENTS

ANALYSIS AND MODELLING OF RAINFALL EVENTS Proeedings of the 14 th Interntionl Conferene on Environmentl Siene nd Tehnology Athens, Greee, 3-5 Septemer 215 ANALYSIS AND MODELLING OF RAINFALL EVENTS IOANNIDIS K., KARAGRIGORIOU A. nd LEKKAS D.F.

More information

Some Aspects of Non-Orthogonal Stagnation-Point Flow towards a Stretching Surface

Some Aspects of Non-Orthogonal Stagnation-Point Flow towards a Stretching Surface Engineering, 00,, 705-709 doi:0.436/eng.00.909 Published Online September 00 (http://www.sirp.org/journl/eng) Some Aspets of Non-Orthogonl Stgntion-Point Flow towrds Strething Surfe Abstrt Mothr Rez, Andi

More information

Electromagnetic-Power-based Modal Classification, Modal Expansion, and Modal Decomposition for Perfect Electric Conductors

Electromagnetic-Power-based Modal Classification, Modal Expansion, and Modal Decomposition for Perfect Electric Conductors LIAN: EM-BASED MODAL CLASSIFICATION EXANSION AND DECOMOSITION FOR EC 1 Eletromgneti-ower-bsed Modl Clssifition Modl Expnsion nd Modl Deomposition for erfet Eletri Condutors Renzun Lin Abstrt Trditionlly

More information

University of Sioux Falls. MAT204/205 Calculus I/II

University of Sioux Falls. MAT204/205 Calculus I/II University of Sioux Flls MAT204/205 Clulus I/II Conepts ddressed: Clulus Textook: Thoms Clulus, 11 th ed., Weir, Hss, Giordno 1. Use stndrd differentition nd integrtion tehniques. Differentition tehniques

More information

First compression (0-6.3 GPa) First decompression ( GPa) Second compression ( GPa) Second decompression (35.

First compression (0-6.3 GPa) First decompression ( GPa) Second compression ( GPa) Second decompression (35. 0.9 First ompression (0-6.3 GP) First deompression (6.3-2.7 GP) Seond ompression (2.7-35.5 GP) Seond deompression (35.5-0 GP) V/V 0 0.7 0.5 0 5 10 15 20 25 30 35 P (GP) Supplementry Figure 1 Compression

More information

Lecture 1 - Introduction and Basic Facts about PDEs

Lecture 1 - Introduction and Basic Facts about PDEs * 18.15 - Introdution to PDEs, Fll 004 Prof. Gigliol Stffilni Leture 1 - Introdution nd Bsi Fts bout PDEs The Content of the Course Definition of Prtil Differentil Eqution (PDE) Liner PDEs VVVVVVVVVVVVVVVVVVVV

More information

Solutions to Assignment 1

Solutions to Assignment 1 MTHE 237 Fll 2015 Solutions to Assignment 1 Problem 1 Find the order of the differentil eqution: t d3 y dt 3 +t2 y = os(t. Is the differentil eqution liner? Is the eqution homogeneous? b Repet the bove

More information

Bridging Methods for Atomistic-to-Continuum Coupling and Their Implementation

Bridging Methods for Atomistic-to-Continuum Coupling and Their Implementation Commun. Comput. Phys. doi:.428/ip.29.9.53 Vol. 7, No. 4, pp. 83-876 April 2 Bridging Methods for Atomisti-to-Continuum Coupling nd Their Implementtion Pblo Seleson nd Mx Gunzburger Deprtment of Sientifi

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

A Matlab/Simulink Model of a Langevin s Ultrasonic Power Transducers

A Matlab/Simulink Model of a Langevin s Ultrasonic Power Transducers Mtl/Simulink Model of Lngevin s Ultrsoni Power Trnsduers Igor Jovnović, Ugleš Jovnović nd Drgn Mnčić strt Ultrsoni sndwih trnsduer, lso known s Lngevin s trnsduer, is hlf-wve resonnt struture tht osilltes

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx,

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx, MATH3403: Green s Funtions, Integrl Equtions nd the Clulus of Vritions 1 Exmples 5 Qu.1 Show tht the extreml funtion of the funtionl I[y] = 1 0 [(y ) + yy + y ] dx, where y(0) = 0 nd y(1) = 1, is y(x)

More information

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4.

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4. Mth 5 Tutoril Week 1 - Jnury 1 1 Nme Setion Tutoril Worksheet 1. Find ll solutions to the liner system by following the given steps x + y + z = x + y + z = 4. y + z = Step 1. Write down the rgumented mtrix

More information

NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE

NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE V.S. Gordeev, G.A. Myskov Russin Federl Nuler Center All-Russi Sientifi Reserh Institute of Experimentl Physis (RFNC-VNIIEF)

More information

The study of dual integral equations with generalized Legendre functions

The study of dual integral equations with generalized Legendre functions J. Mth. Anl. Appl. 34 (5) 75 733 www.elsevier.om/lote/jm The study of dul integrl equtions with generlized Legendre funtions B.M. Singh, J. Rokne,R.S.Dhliwl Deprtment of Mthemtis, The University of Clgry,

More information

Exercise 3 Logic Control

Exercise 3 Logic Control Exerise 3 Logi Control OBJECTIVE The ojetive of this exerise is giving n introdution to pplition of Logi Control System (LCS). Tody, LCS is implemented through Progrmmle Logi Controller (PLC) whih is lled

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

A Mathematical Model for Unemployment-Taking an Action without Delay

A Mathematical Model for Unemployment-Taking an Action without Delay Advnes in Dynmil Systems nd Applitions. ISSN 973-53 Volume Number (7) pp. -8 Reserh Indi Publitions http://www.ripublition.om A Mthemtil Model for Unemployment-Tking n Ation without Dely Gulbnu Pthn Diretorte

More information

Hyers-Ulam stability of Pielou logistic difference equation

Hyers-Ulam stability of Pielou logistic difference equation vilble online t wwwisr-publitionsom/jns J Nonliner Si ppl, 0 (207, 35 322 Reserh rtile Journl Homepge: wwwtjnsom - wwwisr-publitionsom/jns Hyers-Ulm stbility of Pielou logisti differene eqution Soon-Mo

More information

THE DESIGN OF ADAPTIVE CONTROLLER AND SYNCHRONIZER FOR QI-CHEN SYSTEM WITH UNKNOWN PARAMETERS

THE DESIGN OF ADAPTIVE CONTROLLER AND SYNCHRONIZER FOR QI-CHEN SYSTEM WITH UNKNOWN PARAMETERS THE DESIGN OF ADAPTIVE CONTROLLER AND SYNCHRONIZER FOR QI-CHEN SYSTEM WITH UNKNOWN PARAMETERS Sundrpndin Vidynthn 1 1 Reserh nd Development Centre, Vel Teh Dr. RR & Dr. SR Tehnil University Avdi, Chenni-600

More information

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1) Green s Theorem Mth 3B isussion Session Week 8 Notes Februry 8 nd Mrh, 7 Very shortly fter you lerned how to integrte single-vrible funtions, you lerned the Fundmentl Theorem of lulus the wy most integrtion

More information

Table of Content. c 1 / 5

Table of Content. c 1 / 5 Tehnil Informtion - t nd t Temperture for Controlger 03-2018 en Tble of Content Introdution....................................................................... 2 Definitions for t nd t..............................................................

More information

Line Integrals and Entire Functions

Line Integrals and Entire Functions Line Integrls nd Entire Funtions Defining n Integrl for omplex Vlued Funtions In the following setions, our min gol is to show tht every entire funtion n be represented s n everywhere onvergent power series

More information

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix tries Definition of tri mtri is regulr rry of numers enlosed inside rkets SCHOOL OF ENGINEERING & UIL ENVIRONEN Emple he following re ll mtries: ), ) 9, themtis ), d) tries Definition of tri Size of tri

More information

1.3 SCALARS AND VECTORS

1.3 SCALARS AND VECTORS Bridge Course Phy I PUC 24 1.3 SCLRS ND VECTORS Introdution: Physis is the study of nturl phenomen. The study of ny nturl phenomenon involves mesurements. For exmple, the distne etween the plnet erth nd

More information

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

SECTION A STUDENT MATERIAL. Part 1. What and Why.? SECTION A STUDENT MATERIAL Prt Wht nd Wh.? Student Mteril Prt Prolem n > 0 n > 0 Is the onverse true? Prolem If n is even then n is even. If n is even then n is even. Wht nd Wh? Eploring Pure Mths Are

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah 1. Born-Oppenheimer pprox.- energy surfces 2. Men-field (Hrtree-Fock) theory- orbitls 3. Pros nd cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usully does HF-how? 6. Bsis sets nd nottions 7. MPn,

More information

Learning Objectives of Module 2 (Algebra and Calculus) Notes:

Learning Objectives of Module 2 (Algebra and Calculus) Notes: 67 Lerning Ojetives of Module (Alger nd Clulus) Notes:. Lerning units re grouped under three res ( Foundtion Knowledge, Alger nd Clulus ) nd Further Lerning Unit.. Relted lerning ojetives re grouped under

More information

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know

More information

Development of Failure Probability Analysis Method for. Concrete Piers of Multi-span Continuous Bridges using

Development of Failure Probability Analysis Method for. Concrete Piers of Multi-span Continuous Bridges using Development o Filure Probbility Anlysis Method or Conrete Piers o Multi-spn Continuous Bridges using the Probbilisti Cpity Spetrum Method Je Shin CHOI, Je Kwn KIM ABSTRACT When erthqukes our, strutures

More information

, g. Exercise 1. Generator polynomials of a convolutional code, given in binary form, are g. Solution 1.

, g. Exercise 1. Generator polynomials of a convolutional code, given in binary form, are g. Solution 1. Exerise Genertor polynomils of onvolutionl ode, given in binry form, re g, g j g. ) Sketh the enoding iruit. b) Sketh the stte digrm. ) Find the trnsfer funtion T. d) Wht is the minimum free distne of

More information

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 7. Waveguides Part 4: Rectangular and Circular Waveguide

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 7. Waveguides Part 4: Rectangular and Circular Waveguide ECE 5317-6351 Mirowve Engineering Fll 01 Prof. Dvid R. Jkson Dept. of ECE Notes 7 Wveguides Prt 4: Retngulr nd Cirulr Wveguide 1 Retngulr Wveguide One of the erliest wveguides. Still ommon for high power

More information

Plates on elastic foundation

Plates on elastic foundation Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016) Introduction Winkler

More information

Explain shortly the meaning of the following eight words in relation to shells structures.

Explain shortly the meaning of the following eight words in relation to shells structures. Delft University of Technology Fculty of Civil Engineering nd Geosciences Structurl Mechnics Section Write your nme nd study number t the top right-hnd of your work. Exm CIE4143 Shell Anlysis Tuesdy 15

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P.

(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P. Chpter 7: The Riemnn Integrl When the derivtive is introdued, it is not hrd to see tht the it of the differene quotient should be equl to the slope of the tngent line, or when the horizontl xis is time

More information

2. Topic: Summation of Series (Mathematical Induction) When n = 1, L.H.S. = S 1 = u 1 = 3 R.H.S. = 1 (1)(1+1)(4+5) = 3

2. Topic: Summation of Series (Mathematical Induction) When n = 1, L.H.S. = S 1 = u 1 = 3 R.H.S. = 1 (1)(1+1)(4+5) = 3 GCE A Level Otober/November 008 Suggested Solutions Mthemtis H (970/0) version. MATHEMATICS (H) Pper Suggested Solutions. Topi: Definite Integrls From the digrm: Are A = y dx = x Are B = x dy = y dy dx

More information

Can one hear the shape of a drum?

Can one hear the shape of a drum? Cn one her the shpe of drum? After M. K, C. Gordon, D. We, nd S. Wolpert Corentin Lén Università Degli Studi di Torino Diprtimento di Mtemti Giuseppe Peno UNITO Mthemtis Ph.D Seminrs Mondy 23 My 2016 Motivtion:

More information

MATH Final Review

MATH Final Review MATH 1591 - Finl Review November 20, 2005 1 Evlution of Limits 1. the ε δ definition of limit. 2. properties of limits. 3. how to use the diret substitution to find limit. 4. how to use the dividing out

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

Physics 505 Homework No. 11 Solutions S11-1

Physics 505 Homework No. 11 Solutions S11-1 Physis 55 Homework No 11 s S11-1 1 This problem is from the My, 24 Prelims Hydrogen moleule Consider the neutrl hydrogen moleule, H 2 Write down the Hmiltonin keeping only the kineti energy terms nd the

More information

8 THREE PHASE A.C. CIRCUITS

8 THREE PHASE A.C. CIRCUITS 8 THREE PHSE.. IRUITS The signls in hpter 7 were sinusoidl lternting voltges nd urrents of the so-lled single se type. n emf of suh type n e esily generted y rotting single loop of ondutor (or single winding),

More information

VIBRATION ANALYSIS OF AN ISOLATED MASS WITH SIX DEGREES OF FREEDOM Revision G

VIBRATION ANALYSIS OF AN ISOLATED MASS WITH SIX DEGREES OF FREEDOM Revision G B Tom Irvine Emil: tom@virtiondt.om Jnur 8, 3 VIBRATION ANALYSIS OF AN ISOLATED MASS WITH SIX DEGREES OF FREEDOM Revision G Introdution An vionis omponent m e mounted with isoltor grommets, whih t s soft

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

Math 124A October 04, 2011

Math 124A October 04, 2011 Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model

More information

Restraint of purlins for various roof systems

Restraint of purlins for various roof systems NS009 Restrint of purlins for vrious roof systems T. Vrny, M. Brhm & A. Beli ulty of ivil Engineering, zeh Tehnil University, Prh, zehi Astron Buildings S.A., Diekirh, Luxemourg, A memer of the Lind Group

More information

Chapter Gauss Quadrature Rule of Integration

Chapter Gauss Quadrature Rule of Integration Chpter 7. Guss Qudrture Rule o Integrtion Ater reding this hpter, you should e le to:. derive the Guss qudrture method or integrtion nd e le to use it to solve prolems, nd. use Guss qudrture method to

More information

Factors affecting the phonation threshold pressure and frequency

Factors affecting the phonation threshold pressure and frequency 3SC Fctors ffecting the phontion threshold pressure nd frequency Zhoyn Zhng School of Medicine, University of Cliforni Los Angeles, CA, USA My, 9 57 th ASA Meeting, Portlnd, Oregon Acknowledgment: Reserch

More information

Finite Element Simulation on Frictional and Brittle Preseismic fault slip

Finite Element Simulation on Frictional and Brittle Preseismic fault slip Finite Element Simultion on Fritionl nd Brittle Preseismi fult slip Zhishen Wu (1) Yun Go (1) Yutk Murkmi (2) (1) Deprtment of Urn & Civil Engineering. Irki University, Jpn (e-mil: zswu@ip.irki..jp; goyun@hs.irki..jp,

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

Introduction to Olympiad Inequalities

Introduction to Olympiad Inequalities Introdution to Olympid Inequlities Edutionl Studies Progrm HSSP Msshusetts Institute of Tehnology Snj Simonovikj Spring 207 Contents Wrm up nd Am-Gm inequlity 2. Elementry inequlities......................

More information

A Non-parametric Approach in Testing Higher Order Interactions

A Non-parametric Approach in Testing Higher Order Interactions A Non-prmetri Approh in Testing igher Order Intertions G. Bkeerthn Deprtment of Mthemtis, Fulty of Siene Estern University, Chenkldy, Sri Lnk nd S. Smit Deprtment of Crop Siene, Fulty of Agriulture University

More information

Problems set # 3 Physics 169 February 24, 2015

Problems set # 3 Physics 169 February 24, 2015 Prof. Anhordoqui Problems set # 3 Physis 169 Februry 4, 015 1. A point hrge q is loted t the enter of uniform ring hving liner hrge density λ nd rdius, s shown in Fig. 1. Determine the totl eletri flux

More information

Section 3.6. Definite Integrals

Section 3.6. Definite Integrals The Clulus of Funtions of Severl Vribles Setion.6 efinite Integrls We will first define the definite integrl for funtion f : R R nd lter indite how the definition my be extended to funtions of three or

More information

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets.

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets. I MATRIX ALGEBRA INTRODUCTION TO MATRICES Referene : Croft & Dvison, Chpter, Blos, A mtri ti is retngulr rr or lo of numers usull enlosed in rets. A m n mtri hs m rows nd n olumns. Mtri Alger Pge If the

More information

The Emission-Absorption of Energy analyzed by Quantum-Relativity. Abstract

The Emission-Absorption of Energy analyzed by Quantum-Relativity. Abstract The mission-absorption of nergy nlyzed by Quntum-Reltivity Alfred Bennun* & Néstor Ledesm** Abstrt The uslity horizon llows progressive quntifition, from n initil nk prtile, whih yields its energy s blk

More information

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS Throughout, we let [, b] be bounded intervl in R. C 2 ([, b]) denotes the spce of functions with derivtives of second order continuous up to the endpoints. Cc 2

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Math 270A: Numerical Linear Algebra

Math 270A: Numerical Linear Algebra Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner

More information

TOPIC: LINEAR ALGEBRA MATRICES

TOPIC: LINEAR ALGEBRA MATRICES Interntionl Blurete LECTUE NOTES for FUTHE MATHEMATICS Dr TOPIC: LINEA ALGEBA MATICES. DEFINITION OF A MATIX MATIX OPEATIONS.. THE DETEMINANT deta THE INVESE A -... SYSTEMS OF LINEA EQUATIONS. 8. THE AUGMENTED

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

Co-ordinated s-convex Function in the First Sense with Some Hadamard-Type Inequalities

Co-ordinated s-convex Function in the First Sense with Some Hadamard-Type Inequalities Int. J. Contemp. Mth. Sienes, Vol. 3, 008, no. 3, 557-567 Co-ordinted s-convex Funtion in the First Sense with Some Hdmrd-Type Inequlities Mohmmd Alomri nd Mslin Drus Shool o Mthemtil Sienes Fulty o Siene

More information

International Industrial Informatics and Computer Engineering Conference (IIICEC 2015)

International Industrial Informatics and Computer Engineering Conference (IIICEC 2015) Interntionl Industril Informtis nd Computer Engineering Conferene (IIICEC 5) Simultion Implementtion nd Applition Reserh of High Speed Milling Regenertion ype Chtter Bsed on MALAB Yun Lingling, Mei Wentob,

More information

Section 4.4. Green s Theorem

Section 4.4. Green s Theorem The Clulus of Funtions of Severl Vriles Setion 4.4 Green s Theorem Green s theorem is n exmple from fmily of theorems whih onnet line integrls (nd their higher-dimensionl nlogues) with the definite integrls

More information

On the Scale factor of the Universe and Redshift.

On the Scale factor of the Universe and Redshift. On the Sle ftor of the Universe nd Redshift. J. M. unter. john@grvity.uk.om ABSTRACT It is proposed tht there hs been longstnding misunderstnding of the reltionship between sle ftor of the universe nd

More information

THE PYTHAGOREAN THEOREM

THE PYTHAGOREAN THEOREM THE PYTHAGOREAN THEOREM The Pythgoren Theorem is one of the most well-known nd widely used theorems in mthemtis. We will first look t n informl investigtion of the Pythgoren Theorem, nd then pply this

More information

A Study on the Properties of Rational Triangles

A Study on the Properties of Rational Triangles Interntionl Journl of Mthemtis Reserh. ISSN 0976-5840 Volume 6, Numer (04), pp. 8-9 Interntionl Reserh Pulition House http://www.irphouse.om Study on the Properties of Rtionl Tringles M. Q. lm, M.R. Hssn

More information

THE ANALYSIS AND CALCULATION OF ELECTROMAGNETIC FIELD AROUND OVERHEAD POWER LINE HongWang Yang

THE ANALYSIS AND CALCULATION OF ELECTROMAGNETIC FIELD AROUND OVERHEAD POWER LINE HongWang Yang 5th Interntionl Conferene on Advned Mterils nd Computer Siene (ICAMCS 6) THE ANALYSIS AN CALCULATION OF ELECTROMAGNETIC FIEL AROUN OVERHEA POWER LINE HongWng Yng eprtment of eletril engineering, North

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

Hybrid Systems Modeling, Analysis and Control

Hybrid Systems Modeling, Analysis and Control Hyrid Systems Modeling, Anlysis nd Control Rdu Grosu Vienn University of Tehnology Leture 5 Finite Automt s Liner Systems Oservility, Rehility nd More Miniml DFA re Not Miniml NFA (Arnold, Diky nd Nivt

More information

Part 4. Integration (with Proofs)

Part 4. Integration (with Proofs) Prt 4. Integrtion (with Proofs) 4.1 Definition Definition A prtition P of [, b] is finite set of points {x 0, x 1,..., x n } with = x 0 < x 1

More information

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides. 1 PYTHAGORAS THEOREM 1 1 Pythgors Theorem In this setion we will present geometri proof of the fmous theorem of Pythgors. Given right ngled tringle, the squre of the hypotenuse is equl to the sum of the

More information

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS IN-SITU PROBING OF DOMAIN POLING IN Bi 4 Ti 3 O 12 THIN FILMS BY OPTICAL SECOND HARMONIC GENERATION YANIV BARAD, VENKATRAMAN GOPALAN Mterils Reserh Lortory

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

Matrices, Moments and Quadrature, cont d

Matrices, Moments and Quadrature, cont d Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl

More information

Parabola and Catenary Equations for Conductor Height Calculation

Parabola and Catenary Equations for Conductor Height Calculation ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ, 6 (), nr. 3 9 Prbol nd Ctenr Equtions for Condutor Height Clultion Alen HATIBOVIC Abstrt This pper presents new equtions for ondutor height lultion bsed on the

More information

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0. STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

More information

Lecture 6. CMOS Static & Dynamic Logic Gates. Static CMOS Circuit. PMOS Transistors in Series/Parallel Connection

Lecture 6. CMOS Static & Dynamic Logic Gates. Static CMOS Circuit. PMOS Transistors in Series/Parallel Connection NMOS Trnsistors in Series/Prllel onnetion Leture 6 MOS Stti & ynmi Logi Gtes Trnsistors n e thought s swith ontrolled y its gte signl NMOS swith loses when swith ontrol input is high Peter heung eprtment

More information

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums Green s Theorem If f is funtion of one vrible x with derivtive f x) or df dx to the Fundmentl Theorem of lulus, nd [, b] is given intervl then, ording This is not trivil result, onsidering tht b b f x)dx

More information

Generalization of 2-Corner Frequency Source Models Used in SMSIM

Generalization of 2-Corner Frequency Source Models Used in SMSIM Generliztion o 2-Corner Frequeny Soure Models Used in SMSIM Dvid M. Boore 26 Mrh 213, orreted Figure 1 nd 2 legends on 5 April 213, dditionl smll orretions on 29 My 213 Mny o the soure spetr models ville

More information

Engr354: Digital Logic Circuits

Engr354: Digital Logic Circuits Engr354: Digitl Logi Ciruits Chpter 4: Logi Optimiztion Curtis Nelson Logi Optimiztion In hpter 4 you will lern out: Synthesis of logi funtions; Anlysis of logi iruits; Tehniques for deriving minimum-ost

More information

Note 16. Stokes theorem Differential Geometry, 2005

Note 16. Stokes theorem Differential Geometry, 2005 Note 16. Stokes theorem ifferentil Geometry, 2005 Stokes theorem is the centrl result in the theory of integrtion on mnifolds. It gives the reltion between exterior differentition (see Note 14) nd integrtion

More information

Ordinary differential equations

Ordinary differential equations Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties

More information

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL P3.1 Kot Iwmur*, Hiroto Kitgw Jpn Meteorologil Ageny 1. INTRODUCTION Jpn Meteorologil Ageny

More information

Structural Systems. Structural Engineering

Structural Systems. Structural Engineering 101-305 Theor of Strutures 1-1 Instrutor: ST, P, J nd TPS Struturl Engineering 101-305 Theor of Strutures 1 - Instrutor: ST, P, J nd TPS Struturl Sstems Struture the tion of building: ONSTUTION Struturl

More information

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable INTEGRATION NOTE: These notes re supposed to supplement Chpter 4 of the online textbook. 1 Integrls of Complex Vlued funtions of REAL vrible If I is n intervl in R (for exmple I = [, b] or I = (, b)) nd

More information

The Role of Transmitting Boundaries in Modeling Dynamic Soil-Structure Interaction Problems

The Role of Transmitting Boundaries in Modeling Dynamic Soil-Structure Interaction Problems Interntionl Journl of Engineering nd Tehnology Volume No., Ferury, 1 The Role of Trnsmitting Boundries in Modeling Dynmi Soil-Struture Intertion Prolems Mohmmed Y. Ftth *, T. Shnz **, Shth H. Dwood ***

More information

Probability Distributions for Gradient Directions in Uncertain 3D Scalar Fields

Probability Distributions for Gradient Directions in Uncertain 3D Scalar Fields Technicl Report 7.8. Technische Universität München Probbility Distributions for Grdient Directions in Uncertin 3D Sclr Fields Tobis Pfffelmoser, Mihel Mihi, nd Rüdiger Westermnn Computer Grphics nd Visuliztion

More information

Distributed Generation Placement in Unbalanced Distribution System with Seasonal Load Variation

Distributed Generation Placement in Unbalanced Distribution System with Seasonal Load Variation Distriuted Genertion Plement in Unlned Distriution System with Sesonl Lod Vrition Rvi Tej Bhimrsetti Dept. of Eletril Engg., NT Kurukshetr Kurukshetr, ndi svrtej@gmil.om Ashwni Kumr, Memer, EEE Dept. of

More information

COMPARISON BETWEEN TWO FRICTION MODEL PARAMETER ESTIMATION METHODS APPLIED TO CONTROL VALVES. Rodrigo Alvite Romano and Claudio Garcia

COMPARISON BETWEEN TWO FRICTION MODEL PARAMETER ESTIMATION METHODS APPLIED TO CONTROL VALVES. Rodrigo Alvite Romano and Claudio Garcia 8th Interntionl IFAC Symposium on Dynmis nd Control of Proess Systems Preprints Vol., June 6-8, 007, Cnún, Mexio COMPARISON BETWEEN TWO FRICTION MODEL PARAMETER ESTIMATION METHODS APPLIED TO CONTROL VALVES

More information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

NANO-SCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE

NANO-SCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE Nno-Sle Effets in Clindril Contts Sri et l. NANO-SCALE EFFECTS IN THE ADHERENCE, SLIDING AND ROLLING OF A CYLINDER ON A SUBSTRATE Ö. T. Sri, G. G. Adms, S. Müftü Mehnil Engineering Deprtment Northestern

More information

1 Bending of a beam with a rectangular section

1 Bending of a beam with a rectangular section 1 Bending of bem with rectngulr section x3 Episseur b M x 2 x x 1 2h M Figure 1 : Geometry of the bem nd pplied lod The bem in figure 1 hs rectngur section (thickness 2h, width b. The pplied lod is pure

More information

T b a(f) [f ] +. P b a(f) = Conclude that if f is in AC then it is the difference of two monotone absolutely continuous functions.

T b a(f) [f ] +. P b a(f) = Conclude that if f is in AC then it is the difference of two monotone absolutely continuous functions. Rel Vribles, Fll 2014 Problem set 5 Solution suggestions Exerise 1. Let f be bsolutely ontinuous on [, b] Show tht nd T b (f) P b (f) f (x) dx [f ] +. Conlude tht if f is in AC then it is the differene

More information

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 H 3 CH3 C. NMR spectroscopy. Different types of NMR

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 H 3 CH3 C. NMR spectroscopy. Different types of NMR 6.. Spetrosopy NMR spetrosopy Different types of NMR NMR spetrosopy involves intertion of mterils with the lowenergy rdiowve region of the eletromgneti spetrum NMR spetrosopy is the sme tehnology s tht

More information

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials: Summry of equtions chpters 7. To mke current flow you hve to push on the chrges. For most mterils: J E E [] The resistivity is prmeter tht vries more thn 4 orders of mgnitude between silver (.6E-8 Ohm.m)

More information