Static and dynamic rectangular finite elements for plate vibration analysis
|
|
- Christiana Erin Berry
- 5 years ago
- Views:
Transcription
1 UDC 34.: ]:9.6 Originl scientific pper Received: Sttic nd dynmic rectngulr finite elements for plte virtion nlysis Ivo Senjnović, Mrko Tomić, Neven Hdžić, Nikol Vldimir University of Zgre, Fculty of Mechnicl Engineering nd Nvl Architecture, Zgre, CROATIA e-mil: SUMMARY Nowdys dynmic finite elements with frequency dependent unified stiffness mtrix re considered in order to reduce finite element mesh density, nd, consequently, system of lgeric equtions in eigenvlue formultion. In this pper, n ttempt to develop 4-noded rectngulr finite element for thin plte ending is presented. Since there is no simple generl nlyticl solution of plte virtions for ritrry oundry conditions, dynmic shpe functions re ssumed s products of em shpe functions in longitudinl nd trnsverse direction. Dynmic finite element with frequency dependent stiffness nd mss mtrix is derived y energy pproch. Convergence of results is checked nd it is found tht such finite element formultion is not suitle for wider use. Therefore, strip finite element for virtion nlysis of rectngulr plte with simply supported two opposite edges is developed, which gives exct results. Appliction of oth finite elements is illustrted y severl numericl exmples. KEY WORDS: thin plte, dynmic finite elements, rectngulr element, strip element, dynmic stiffness nd mss mtrix.. INTRODUCTION The finite element method (FEM) is numericl method in which physicl properties from the element re re condensed into its nodes. In tht wy, the governing differentil equtions descriing the ehviour of n element re trnsformed into the system of lgeric equtions sed on energy equivlence []. FEM is widely used for structurl nlysis of engineering structures due to the possiility of detiled modelling of complex structure topology nd reltively esy computing. Depending on the type of structure, different finite elements hve een developed for strength nlysis, virtion nlysis, stility nlysis, etc. []. At the very eginning of computer ppliction, mtrix methods hve een developed s sis for lter development of nowdys generlly ccepted FEM [, 3]. Two FEM formultions hve een developed, i.e. the displcement finite element method nd force finite element method. In the former nd ltter cse, nodl ENGINEERING MODELLING 9 (06) -4, -
2 displcements nd nodl forces re unknown, respectively. Nowdys, the displcement FEM is mostly used s more prcticl solution []. This pper dels with thin plte virtion nlysis [4]. For this purpose, two types of 4-noded rectngulr finite elements re used; one with deflection comptiility (non-conforming) nd the other with deflection nd slope comptiility (conforming) []. In oth cses, shpe functions re presented in polynomil form. Sttic stiffness nd mss mtrix re developed. Numericl solution converges to the exct one y incresing finite element mesh density. Nowdys, there is tendency to develop dynmic finite elements with dynmic stiffness. For this purpose, nlyticl solution of differentil eqution of motion is used. However, tht is reltively esy tsk if D prolems s rs nd ems re considered, or if D prolem cn e reduced to D prolem, s in the cse of plte with simply supported two opposite edges []. An dvntge of dynmic finite elements is the pplicility of corser finite element mesh nd high ccurcy. In this pper, n ttempt to develop dynmic 4-noded rectngulr finite element with frequency dependent stiffness nd mss mtrix is descried. Shpe functions re ssumed in the form of products of the nlyticl em shpe functions. Since such plte shpe functions re n pproximtion of the rel ones, the usul energy pproch is pplied. Also, dynmic strip finite element for virtion nlysis of pltes with simply supported two opposite edges is presented. Its ppliction is illustrted in severl numericl exmples. Depending on oundry conditions, oth displcement nd force finite element methods re used in order to void the ill-conditioned formultion of eigenvlue prolem nd rther complicted procedure for its solving [6].. STATE-OF-THE ART. GENERAL There re two sic 4-noded rectngulr plte ending finite elements. One stisfies only deflection comptiility conditions long the edges, nd the other oth deflection nd slope comptiility conditions. The former is clled non-conforming element nd my yield over or under-stiff results. The ltter is clled conforming element nd demonstrtes monotonic convergence from the over-stiff side to the exct solution [, ]. Shpe functions of the non-conforming element re expressed in polynomil terms sed on sttic plte deformtion. Shpe functions of conforming element re defined s product of em shpe functions in longitudinl nd trnsversl directions. The element properties, i.e. stiffness mtrix, mss mtrix, geometric stiffness mtrix nd lod vector for sttic, dynmic nd stility nlyses, re sed on the vrition principle or energy pproch, nd re minly presented in symolic form. ENGINEERING MODELLING 9 (06) -4, -
3 . FINITE ELEMENT WITH DEFLECTION COMPATIBILITY The 4-noded rectngulr finite element is shown in Figure in Crtesin (0, x, y) nd nturl (0, ξ, η) coordinte systems. Fig. Rectngulr plte ending finite element with origin in node Nodl displcements, i.e. nodl deflection nd rottion ngles re indicted. Deflection function is ssumed in the polynomil form involving terms equl to the numer of degrees of freedom (DOF): w = + ξ + η + ξ + ξη + η + ξ + ξ η + ξη + η + ξ η + ξη, () where ξ = x /, η = y / nd nd re finite element length nd redth, respectively. In Eq. () ll polynomil terms of the third order nd lower re included. Among terms of the 3 3 fourth order only two mixed terms i.e., ξ η nd ξη re tken into ccount s the est choice from the ccurcy point of view. Rottion ngles re specified s: ϕ w w w w = =, ψ = =. () y η x ξ Shpe functions, i.e. deflection functions due to unit vlue of one nodl displcement (deflection or rottion) nd zero vlue of the others, re otined in the following form: ENGINEERING MODELLING 9 (06) -4, - 3
4 { Φ} ξη (3 ξ )ξ ( η) ( ξ )(3 η)η ( ξ )η( η) ξ( ξ ) ( η) (3 ξ )ξ ( η) + ξη( η)( η) ξη( η) ( ξ )ξ ( η) = (3 ξ )ξ η ξη( η)( η) ξ( η)η ( ξ )ξ η ( ξ )(3 η)η + ξ( ξ )( ξ )η ( ξ )( η)η ξ( ξ ) η Shpe functions due to unit nodl displcements w i re of order for unit nodl rottion ngle ϕ nd ψ, i.e. 3 ξ η nd (3) 3 ξη, s well s those 3 ξη ndξ 3 η, respectively. The ove shpe functions nd corresponding stiffness mtrix in symolic form re shown in Refs. [, 7]. Moreover, stiffness nd mss mtrix in symolic form re presented in Refs. [4, 8, 9]..3 FINITE ELEMENT WITH DEFLECTION AND SLOPE COMPATIBILITY Finite element shown in Figure is considered. Shpe functions Φ j(ξ,η), j =,..., re specified s products of em shpe functions in ξ nd η directions, X i(ξ ) nd Y i(η), i =,, 3, 4 []: Φ = X Y = (+ ξ )( ξ ) (+ η)( η) Φ = X Y = (+ ξ )( ξ ) η( η) Φ = X Y = ξ( ξ ) (+ η)( η) Φ = X Y = (3 ξ )ξ (+ η)( η) Φ = X Y = (3 ξ )ξ η( η) Φ = X Y = ( ξ )ξ (+ η)( η) Φ = X Y = (3 ξ )ξ (3 η)η Φ = X Y In this cse, ll shpe functions re of = (3 ξ )ξ ( η)η Φ = X Y = ( ξ )ξ (3 η)η Φ = X Y = (+ ξ )( ξ )(3 η)η Φ = X Y = (+ ξ )( ξ ) ( η)η Φ = X Y = ξ( ξ ) (3 η)η. 3 3 ξ η order. Stiffness mtrix is shown in the symolic form in Ref. [], while oth stiffness nd mss mtrices re specified in Ref. [9]. (4) 4 ENGINEERING MODELLING 9 (06) -4, -
5 3. FINITE ELEMENT BASED ON FREE PLATE NATURAL MODES Finite element of virting plte cn e considered s free structurl element in spce. Therefore, it is resonle to ssume its deflection function s set of free plte nturl modes. Only nturl virtions of pltes with simply supported two opposite edges nd ritrry oundry conditions on the remining two edges cn e determined nlyticlly in reltively simple wy. Strictly speking, some other oundry vlue prolems re lso solved nlyticlly ut in very complicted mnner [0, ]. Therefore, for ritrry oundry conditions, numericl methods like Ryleigh-Ritz method, finite difference method nd finite element method (FEM) re preferred. However, nturl virtions of free plte cn lso e nlysed pproximtely in simple nlyticl wy y the Ryleigh quotient, if the mode shpes re ssumed in the form of products of free em nturl modes, W(ξ,η) = X(ξ )Y(η) []. Bem nturl modes re divided into two groups, i.e. symmetric nd ntisymmetric. Therefore, the origin of the coordinte system is locted in the middle of the em length, l =. Free em nturl modes re expressed y the nturl coordinte ξ = x /, s follows []: symmetric: ntisymmetric: chγiξ cosγiξ X0 =, X i = +, i =,4... γ =.360, γ4 =.4978,... chγi cos γi shγiξ sinγiξ X = ξ, X i = +, i = 3,... γ3 = 3.966, γ = ,... shγi sinγi The sme expressions re vlid for em nturl modes in the cross-direction η = y /, denoted s Y i(η). The first three plte nturl modes re due to rigid ody motion with zero nturl frequencies since there is no restoring stiffness. Finite element virtion nlysis of free squre plte shows tht elstic nturl modes re of ordinry form []: W (ξ,η) X (ξ )Y (η) ij = i j, (7) () (6) or extrordinry form: W ± (ξ,η) = X (ξ )Y (η) ± X (ξ )Y (η). (8) mn kl m n k l The first elstic nturl modes re shown in Figure nd the corresponding formule re specified in Tle. Shpes of the nd nd 3 rd modes re hyperolic nd n ellipticl proloid, respectively []. Dimensionless frequency prmeters determined y oth FEM nlysis nd the Ryleigh quotient (RQ) re listed in Tle. Their ccurcy is estimted with respect to the relile results determined y the Ryleigh-Ritz method (RRM) in Ref. [3]. ENGINEERING MODELLING 9 (06) -4, -
6 Tle Frequency prmeter ( ) Ω = ω π ρh D of free squre plte elstic modes, Figure 3 Mode no. Mode type FEM Discrepncy Discrepncy RRM, RQ % % Liew [3] XY XY0 X0Y XY0 + X0Y XY XY X3Y0 X0Y X3Y0 + X0Y XY X3Y XY X3Y + XY XY (.0)* X3Y (.0)* * With respect to FEM Fig. Free squre plte nturl modes, FEM 6 ENGINEERING MODELLING 9 (06) -4, -
7 Fig. 3 Rectngulr plte ending finite element in Crtesin nd nturl coordinte systems It is necessry to point out tht extrordinry modes for rectngulr plte move to higher frequencies y incresing vlue of the plte spect rtio /. Therefore, deflection of rectngulr finite element shown in Figure 3 is ssumed s set of three rigid ody modes nd 9 squre plte elstic modes. Ordinry modes nd constitutive prts of extrordinry modes re tken into ccount. Symmetric mode X Y is excluded. Hence, one cn write: where { } is vector of unknown coefficients nd: { } W (ξ,η) = P(ξ,η), (9) P(ξ,η) = ξ η ξη X Y X3 Xη ξy Y3 X3η ξy3. (0) Structure of P(ξ,η) is similr to the polynomil terms in the cse of finite element with deflection comptiility conditions, Eq. (). For rottion ngles, the following is otined: w P ϕ = = η η { } w P ψ = = { }, ξ ξ () ENGINEERING MODELLING 9 (06) -4, - 7
8 P η P ξ = = ' ' ' ' ξ 0 Y 0 X ξy Y X ξy ' ' ' ' ' η X 0 X X η Y 0 X η Y. () Displcement field is specified s: w ψ { f } = ϕ = [ Z(ξ,η)]{ } [ Z(ξ,η)] Displcement vector of the i-th node ccording to Eq. (3) reds: { } ϕ [ ] { } i nd the finite element displcement vector:, (3) P P =. (4) η P ξ wi δ = i i = Z, i =,,3,4 i ψ { } { } { } [ ]{ } i () δ = δ = Z, (6) [ Z] [ Z] [ Z] [ Z] [ Z] =. (7) 3 4 By sustituting {} from Eq. (6) into Eq. (9) one cn express deflection function y the nodl displcement vector: is vector of shpe functions. { δ} w(ξ,η) = Φ, (8) Φ [ ] Deformtion (curvture) vector is otined in the form: = P(ξ,η) Z (9) 8 ENGINEERING MODELLING 9 (06) -4, -
9 w x w y w x y { κ} = = [ B]{ δ} [ B] [ A][ Z], (0) = () [ A] P ξ p P = q = η P r ξ η " " " " 3 3 p = X 0 X X η 0 0 X η 0 " " " " 3 3 q = Y 0 0 ξy Y 0 ξy () (3) ' ' ' ' 3 3 r = X Y 0 X Y. According to the definition of stiffness mtrix []: T [ ] [ ] [ ][ ] K = B D B d A, (4) A 3 ν 0 Eh D = D d, D =, d = ν 0 ( ν ) ν 0 0 [ ] [ ] [ ] () is mtrix of flexurl rigidity, the following is otined y employing Eq. (): [ ] [ ] T ( ) [ ][ ] =, (6) K D Z k Z T [ ] [ ] [ ][ ] k A d A dξ dη =. (7) Tking into ccount Eq. (), the element of mtrix [k] cn e written in index nottion: ENGINEERING MODELLING 9 (06) -4, - 9
10 i j = i j + j + i j + j + i j, (8) k p p νq q q νp ( ν)r r where p i, q i nd r i re defined y Eq. (3). Mss mtrix is expressed with shpe functions: [ ] { Φ} M = m Φ d A, (9) A where m = ρh is mss per unit plte re. By sustituting Eq. (9) into Eq. (9), the following is otined: The finite element eqution reds: [ ] [ ] T ( ) [ ][ ] =, (30) M m Z J Z [ J] { } = P P dξ dη. (3) { F} ([ K ] ω [ M] ){ δ} =, (3) where { F } nd { δ } re vectors of nodl forces nd nodl displcements, respectively. The first two free em nturl modes, Eqs. () nd (6), re shown in Figure 4. Virtion nlysis of free squre plte modelled y the presented finite element shows tht ll nturl frequencies determined y only one finite element gree very well with the exct solution, wheres for dense finite element mesh results converge to slightly higher vlues of the exct ones, Tle. Therefore, em nturl modes, Eqs. () nd (6), re pproximted y polynomils: = ( ), X3 ξ ( ξ 3) X 3ξ = (33) nd compred in Figure 4. Discrepncies of nturl frequencies of free squre plte determined y one finite element re very lrge for higher modes. However, the results rpidly converge to the exct solution y incresing the finite element mesh density, Tle. 0 ENGINEERING MODELLING 9 (06) -4, -
11 Tle Frequency prmeter Ω = ω( / π ) ρh/ D of free squre plte, Figure 3 Free plte Polynomil, mesh Mode No. Liew [3] modes, mesh x x x 4x4 6x6 8x8 0x Fig. 4 The first two free em elstic modes nd their polynomil pproximtion Bsed on the foresid numericl test, the polynomilly pproximted em nturl modes, Eq. (33), re used for the derivtion of finite element properties. Tht lso enles their symolicl presenttion. Mtrix of the polynomil coefficients nd its inversion red: ENGINEERING MODELLING 9 (06) -4, -
12 [ Z] [ Z ] = = Shpe functions, y using Eq. (9), cn e presented in index nottion for ech node i =,, 3, 4: Φi = + ξiξ(3 ξ ) + ηiη(3 η ) + ξiηiξη(4 ξ η ) 8 Φ i = (+ ξiξ + ηiη + ξiηiξη)(η ) 8 Φ i3 = (+ ξiξ + ηiη + ξiηiξη)( ξ ). 8 Shpe functions of the finite element re numered s follows: where k =,.... k i j (34) (3) (36) Φ = Φ, i =,,3,4, j =,,3, (37) Stiffness mtrix [K] is determined ccording to Eqs. (6), (7) nd (8), nd mss mtrix [M] ccording to Eqs. (30) nd (3). Mtrix [J] in [M] is digonl since polynomil terms in Eq. (0), where X (ξ ) nd X 3(ξ ) re specified y Eq. (33), nd Y (η) nd Y 3(η) in n nlogous wy, re orthogonl, i.e. ENGINEERING MODELLING 9 (06) -4, -
13 digonl [ ] J = (38) Stiffness nd mss mtrix cn e derived in symolic form. They re the sme s those presented in Refs. [4, 9], since polynomil Eq. (0) with functions Eq. (33) cn e rerrnged into the form Eq. (). The dvntge of the present pproch is the physicl explntion of the 3 polynomil terms s pproximte free plte nturl modes. As result, polynomil terms ξ η nd 3 ξη re chosen insted of 4 ξ nd 4 η due to lower nturl frequencies, i.e. strin energy, Tle. Another dvntge of the present formultion is the specifiction of the shpe functions in index nottion, which fcilittes the progrmming of the finite element. 4. FINITE ELEMENT BASED ON BEAM NATURAL MODES Finite element with the origin of the coordinte system plced in node is considered, Figure. Its shpe functions cn e specified s the product of em shpe functions X i(ξ ) nd Y i(η), i =,, 3, 4, s follows: Plte deflection is presented in the form: Node Node Node 3 Node4 Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y Φ = X Y { }. (39) w(ξ,η) = Φ δ, (40) where Φ is the vector of shpe functions nd { δ } is the vector of nodl displcements. Curvture vector is defined s: w x w y w x y { κ} = = [ B]{ δ} [ B] Φ ξ p Φ = q = η Φ r ξ η, (4) (4) ENGINEERING MODELLING 9 (06) -4, - 3
14 " " " " " " " " " " " " p = X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y " " " " " " " " " " " " q = X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' r = X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y. (43) Stiffness mtrix is defined ccording to Eq. (4): where referring to Eq. (8):, (44) 0 0 T [ K ] = [ B] [ D][ B] dξ dη = D[ k] i j = i j + i j + ( i j + i j ) + i j. (4) 0 0 k p p q q ν p q q p ( ν )r r dξdη Prmeters p i, q i nd r i re specified y Eq. (43). For mss mtrix Eq. (9), one cn write: [ ] { Φ} M m Φ dξ dη =, (46) 0 0 In order to determine em shpe functions, the deflection of nturlly virting em with length l nd nturl coordinte ξ = x / l, 0 ξ, is expressed y Krylov functions U i(ξ ) due to simplicity: X = AU + A U + A3 U3 + A4 U4, (47) U = ( chγξ + cosγξ ) U = ( shγξ + sinγξ ) U3 = ( chγξ cosγξ ) U4 = ( shγξ sinγξ ). An dvntge of Krylov functions is their chrcteristic vlues t ξ = 0 : U (0) =, 3 4 nd reltion etween derivtives nd functions: (48) U (0) = U (0) = U (0) = 0 (49) ' ' ' ' = 4 = 3 = 4 = 3 " " " " = 3 = 4 3 = 4 = U γu, U γu, U γu, U γu U γ U, U γ U, U γ U, U γ U,etc. Hence, derivtives of Krylov functions cn e expressed y the sme functions with the chnge of indexes in the nticlockwise direction. Furthermore, one otins for the ngle of rottion: (0) ϕ = dx γ ( A U A U A U A U ). dx = () 4 ENGINEERING MODELLING 9 (06) -4, -
15 The integrtion constnts A i cn e expressed in terms of nodl displcements. For the st node ξ = 0 nd one finds from Eqs. (47) nd (0) A = X(0) = W, nd A = ϕ(0 ) = φ γ γ. Constnts A 3 nd A 4 re expressed y A nd A from the sme equtions setting for the nd node X() = W nd ϕ () = φ. Hence, deflection function Eq. (47) cn e written in the form: X(ξ ) = X (ξ )W + X (ξ )φ + X 3(ξ )W + X 4(ξ )φ, () where X i(ξ ) refers to shpe functions with unit vlue of deflection nd rottion ngle t ξ = 0 nd ξ =, respectively: with the following coefficients: X (ξ ) = U (ξ ) + U (ξ ) + U (ξ ) 3 4 X (ξ ) = [ U (ξ ) + U 3(ξ ) + U 4(ξ )] γ X (ξ ) = U (ξ ) + U (ξ ) X 4(ξ ) = [ 4U 3(ξ ) + 4U 4(ξ )] γ shγ sinγ = U ()U 3() U 4 () =, C chγcosγ chγ sinγ shγcosγ = [ U ()U 4() U ()U 3() ] =, C chγcosγ chγ cos γ shγ sinγ 3 = U 3() =, 4 = U 4() =, C chγcosγ C chγcosγ shγcosγ + chγ sinγ = [ U 3()U 4() U ()U () ] =, C chγcosγ shγ sinγ = U ()U 3() U () =, C chγcosγ shγ + sinγ chγ cos γ 3 = U () =, 4 = U 3() =, C chγcosγ C chγcosγ C = U ()U 4() U 3 () = ( chγcosγ ). The chrcteristic feture of the first shpe function is tht sher force: tkes zero vlue t ξ 0 3 d X D 3 Q( ξ ) = D = γ U ( ξ ) + U ( ξ ) + U ( ξ ), 3 3 dx =, i.e. ( ) 3 3 Q 0 Dγ 0 (3) (4) 4 () = =. This condition is relted to ccording to Eq. (4), one otins eigenvlue eqution thγ + tgγ = 0, with roots Eq. (). In similr wy, ending moment of the second shpe function: tkes zero vlue t ξ 0 d X D M ( ξ ) = D = γ U ( ξ ) + U ( ξ ) + U ( ξ ), dx =, i.e. ( ) M 0 Dγ 0 = 0 nd 4 (6) = =. Tking into ccount tht eigenvlue function thγ tgγ = 0 is otined with eigenvlues Eq. (6). = 0, Eq. (4), ENGINEERING MODELLING 9 (06) -4, -
16 Finite element is constructed with em shpe functions X ( ξ ) nd X3 ( ) =, nd X ( ξ ) nd ( ) γ γ.360 ξ with rgument X4 ξ with rgument 3 =. Similr expressions to Eq. = with corresponding (3) cn e written for em functions in η direction Yi ( η ), i,,3,4 Krylov function Vj ( η ), j,,3,4 =. For derivtives of Krylov functions, reltions in Eq. (0) re ville. Stiffness nd mss mtrix re specified y Eqs. (44) nd (46), respectively. The first two em shpe functions expressed y Krylov functions re shown in Figures nd 6 nd compred with ordinry em shpe functions known s Hermite polynomils []. The " difference etween X i(ξ ) vlues is smll, wheres the one etween X i (ξ ) vlues is lrger. Polynomil " i X (ξ ) function is stright line, nd tht expressed y Krylov functions is of " prolic shpe. Curvture functions X i (ξ ) influence the stiffness mtrix nd deflection functions X i(ξ ) influence mss mtrix. If shpe function X (ξ ) is defined with γ, it is very close to the polynomil shpe function, Figure 6. For γ 0., shpe functions in terms of Krylov functions give the sme numericl vlues s Hermite polynomils, Figures nd 6. If trigonometric nd hyperolic functions in Krylov functions Eq. (48) nd coefficients Eq. (4) re expnded into power series nd setting γ 0, shpe functions Eq. (3) re reduced to polynomil em shpe functions in Eq. (4). Fig. Comprison of the first em shpe function determined y Krylov functions nd polynomil 6 ENGINEERING MODELLING 9 (06) -4, -
17 Fig. 6 Comprison of the second em shpe function determined y Krylov functions nd polynomil In order to investigte the ccurcy of different finite element formultions, virtion nlysis of clmped squre plte is performed. Plte is modelled y different mesh density, nd corresponding frctions of prmeters γ nd γ 3. The otined results for the first 8 nturl modes re shown in Tle 3, nd compred with the Liew vlues [3]. Mesh x with γ nd γ 3 gives only the first three nturl frequencies due to only one centrl free node with 3 d.o.f., which re highly ccurte. For mesh 4x4 discrepncies re incresed, ut vlues of frequency prmeter converge to the ccurte vlues, Tle 3. Tle 3 Frequency prmeter Ω = ω( / π ) functions, γ =.360, γ3 = ρh/ D of clmped squre plte, dynmic shpe Mesh Mode Mode Mode 3 Mode 4 Mode Mode 6 Mode 7 Mode 8 x γ, γγ (0.34%) (0.47%) (0.47%) 4x4 γ/, γγ3/ x6 γ/3, γγ3/ x8 γ/4, γγ3/ x0 γ/, γγ3/ x0 γ/0, γγ3/ x30 γ/, γγ3/ x40 γ/0, γγ3/ Liew [3] ENGINEERING MODELLING 9 (06) -4, - 7
18 If γ = 0. is tken into ccount, vlues of frequency prmeter re equl to those determined y the sttic finite element with polynomil shpe functions, Eqs. (4), Tle 4. Vlues of frequency prmeter converge to the exct solution y incresing finite element mesh density. Polynomil shpe functions cn e presented in simpler form, if the coordinte system of the finite element is plced in the middle of the element, Figure 3. In tht cse one otins: where ξ nd η. Xi = + ξiξ(3 ξ ), i =,3 4 Xi = ξ i(4 + 3ξiξ )(ξ ), i =,4 6 Yi = + ηiη(3 η ), i =,3 4 Yi = η i(4 + 3ηiη)(η ), i =,4, 6 Stiffness nd mss mtrix, Eqs. (44) nd (46), cn e generted in symolic form. They re the sme s those in Ref. [], since polynomils Eq. (7) cn e trnsformed into form Eq. (4) y chnging the coordinte system. Tle 4 Frequency prmeter Ω = ω( / π ) ρh/ D of clmped squre plte, shpe functions in terms of Krylov functions ( γ = 0.) nd Hermite polynomils (7) Mesh Mode Mode Mode 3 Mode 4 Mode Mode 6 Mode 7 Mode 8 x x x x x x x x Liew [3] STRIP FINITE ELEMENT FOR SIMPLY SUPPORTED PLATE AT TWO OPPOSITE EDGES. THEORETICAL CONSIDERATION Differentil eqution of plte virtion reds [4]: W W W m + + ω W = x x y y D In cse of simply supported longitudinl edges one cn ssume deflection in the form: (8) W ( ξ,η) = X ( ξ ) sinγη, γ = nπ, n =,, K (9) 8 ENGINEERING MODELLING 9 (06) -4, -
19 By sustituting Eq. (9) into Eq. (8), ordinry differentil eqution of Lévy form is otined [4]: m X γ X + γ ω X = 0. D Solution of Eq. (60) is ssumed in the exponentil form X = e nd four roots re otined, i.e. λ, = ± α, λ3,4 = ± iβ, λξ (60) m m α = ω + γ, β ω γ. D = D (6) Hence, deflection function Eq. (9) nd rottion ngle Eq. () red: W W ( ξ,η) = X ( ξ ) sinγη, ψ( ξ,η) = = Ψ ( ξ ) sinγη, (6) x ( ) X ξ = A shγη + A chγη + A sinβξ + A cosβξ, 3 4 (63) Ψ ( ξ ) = ( A αchαξ + A αshαξ + A3 βcosβξ A4 β sinβξ ) sinγη. Bending moment nd totl shering force, ccording to Ref. [4], red: W W 0 Mx ( ξ,η) = D + ν = Mx ( ξ ) sinγη, x y (64) 3 3 W W 0 Qx ( ξ,η) = D + ( ν) Qx ( ξ ) sinγη, 3 = x x y nd: ( ) = ( + ) ( ) = ( + + ) 0 x 3 4 M ξ D A c shαξ A c chαξ A c sinβξ A c cosβξ, 0 x 3 4 Q ξ D A d chαξ A d shαξ A d cosβξ A d sinβξ, α γ β γ c = ν, c ν, = + α γ α β γ β d = ( ν ), d ( ν ). = + Nodl displcements X ( 0) = X, Ψ ( 0) = Ψ, X ( ) = X nd Ψ ( ) Ψ (6) (66) = cn e written in the mtrix nottion: 0 0 X α β A 0 0 Ψ A =. (67) X shα chα sinβ cos β A 3 Ψ α α β β A 4 chα shα cosβ sinβ ENGINEERING MODELLING 9 (06) -4, - 9
20 On the other side nodl forces Q ( 0) = Q, M ( 0) = M, Q ( ) = Q nd ( ) lso presented in the mtrix nottion: x x x Q d 0 d 0 A M 0 c 0 c A = D. Q dchα dshα d cosβ d sinβ A3 M cshα cchα csinβ ccosβ A 4 M = M re x Equtions (67) nd (68) cn e written in condensed mtrix form { δ} [ B]{ A} { F} = [ C]{ A}, nd y eliminting vector { A }, the strip finite element eqution is otined: { } ( ) { } (68) = nd F = K ω δ, (69) ( ) [ ][ ] K ω = C B, (70) is 4x4 symmetric dynmic stiffness mtrix with 6 independent terms shown in the Appendix. Eqution (69) is derived y the displcement finite element method. In some cses, it is more convenient to use the finite element eqution otined y the force finite element pproch, i.e. { } ( ) { } δ = L ω F, (7) ( ) = ( ) = [ ][ ] L ω K ω B C, (7) is dynmic flexiility mtrix []. Its terms re lso given in the Appendix. A plte with uniform stiffness nd mss distriution is modelled y only one or two strip elements, depending on oundry conditions. For ech vlue of γ = nπ, n =,, K infinite numer of exct nturl frequencies nd modes cn e theoreticlly determined. If nturl virtions of plte with vrile properties in longitudinl direction is nlysed, plte is modelled y set of strip finite elements. Nturl frequencies re determined from condition Det Ks ( ω) = 0.. ILLUSTRATIVE EXAMPLES Simply supported squre plte, SSSS Only one strip finite element is sufficient to solve this tsk with comined geometric nd physicl oundry conditions. By tking into ccount oundry conditions X = X = 0, the finite element eqution is reduced to two d.o.f. system: M k k = M k k Ψ. Ψ Since k44 = k nd k4 = k4 one cn write for the determinnt of Eq. (73): 4 4 (73) Det = k + k k k = 0. (74) Ech of the ove expressions in the rckets cn e equl to zero. By pplying formule (A) from the Appendix, the following is otined: 0 ENGINEERING MODELLING 9 (06) -4, -
21 k ± k = ( c + c ) α sinβ( chαm ) βshα ( m cos β) = 0. (7) 4 Dk Deflection of simply supported plte is descried only y trigonometric functions, nd excluding hyperolic functions from Eq. (7), one finds sinβ = 0, i.e. β = mπ. Finlly, y employing the second formul of Eq. (6) for β, the well-known formul for nturl frequencies of simply supported plte is otined: Dimensionless frequency prmeter: tkes the simple form: mπ nπ D ω = +. (76) m Ω = ω π Ω = m + n m, (77) D. (78) Vlues of Ω for squre plte re compred with those determined numericlly y Ryleigh- Ritz method in Ref. [3], Tle. Tle Frequency prmeter ( ) Ω = ω π ρh D of squre plte simply supported t longitudinl edges Mode no. SSSS CSCS FSFS Liew, m x n Exct [3] n Liew, Exct [3] n Liew, Exct [3] x x x x x x x x x x Clmped squre plte t trnsverse edges, CSCS This prolem with geometric oundry conditions t trnsverse edges cn e solved with t lest two strip finite elements. By tking into ccount oundry conditions, the reduced ssemly eqution reds: ( ) ( ) ( ) ( ) Q k33 + k k34 + k X =, M ( ) ( ) ( ) ( ) k 43 k k44 k Ψ + + (79) ENGINEERING MODELLING 9 (06) -4, -
22 where ( p) k, p =, re relted to strip elements, nd Q, M, X nd Ψ re relted to the ij middle node of the ssemly. Since two finite elements re equl, one cn write ccording to Eqs. (A): Determinnt of Eq. (79) reds: ( ) ( ) ( ) ( ) ( ) ( ) k = k = k, k = k = k, k = k. (80) Det = 4k k = d α + d β c + c DK αshα cos β + βchα sinβ αchα sinβ βshα cos β = 0. Two trnscendentl frequency equtions re otined from Eq. (8), i.e. (8) αthα + βtgβ=0, βthα αtgβ=0. (8) The former nd the ltter re relted to symmetric nd ntisymmetric modes, respectively. Tking into ccount Eqs. (6), vlues of frequency prmeter Ω, Eq. (77), re determined for squre plte nd listed in Tle. Tht exct solutions re compred with numericl vlues clculted y the Ryleigh-Ritz method in Ref. [3]. Free squre plte t trnsverse edges, FSFS Due to physicl oundry conditions t the trnsverse edges, it is more convenient to solve this tsk y using the forced finite element formultion nd t lest two strip elements. Tking oundry conditions into ccount, the reduced ssemly eqution tkes the form: ( ) ( ) ( ) ( ) X l33 + l l34 + l Q =, Ψ ( ) ( ) ( ) ( ) l M 43 l l44 l + + ( p) where l, p =, re elements of the dynmic flexiility mtrix Eq. (79) for two elements. In ij similrity to the previous exmple, ccording to Eqs. (A): nd the determinnt of Eq. (84) reds: (83) ( ) ( ) ( ) ( ) ( ) ( ) l = l = l, l = l = l, l = l, (84) Det = 4l l = d α + d β c + c DL c d chα sinβ c d shαcosβ c d chα sinβ + c d shαcosβ = 0. Two trnscendentl frequency equtions re otined from Eq. (8), i.e. (8) cdthα cdtgβ=0, cdthα + cdtgβ=0. (86) The former nd ltter re relted to ntisymmetric nd symmetric modes, respectively. By employing reltions (6), vlues of the frequency prmeter Ω, Eq. (77), re determined for the squre plte nd included in Tle. They re compred with numericlly determined vlues from Ref. [3]. In the literture, the ove tsk is solved y employing ordinry used displcement finite element method nd only one finite element. However, in tht cse one is fced with the prolem of nturl frequencies determintion since Det K ( ω) 0. This prolem is overcome ENGINEERING MODELLING 9 (06) -4, -
23 y specil nd rther complex Wittrick-Willims lgorithm [, 6]. However, y using the force finite element formultion for similr prolems, the solution is otined stright wy nd unnecessry complictions re voided. 6. CONCLUSION Nowdys, dynmic finite elements with frequency dependent properties re more nd more used, due to course finite element mesh nd high ccurcy. If nlyticl solution of the considered issue is pplied, unique dynmic stiffness mtrix is constructed nd exct results re otined. If free em nturl modes re used for descriing plte deflection, the energy pproch is employed nd sttic stiffness nd mss mtrix re derived. The ltter pproch is used for the construction of non-conforming 4-noded rectngulr plte finite element (Section 3). Course mesh gives pproximte results for the spectrum of nturl frequencies quite close to the exct ones. In cse of conforming finite element with shpe functions s product of dynmic em shpe functions, only the first few nturl frequencies re close to the exct vlues (Section 4). By incresing finite element mesh density nd decresing vlues of rgument of trigonometric nd hyperolic functions ccordingly, the otined results converge to the exct ones. If very smll vlue of tht rgument is used, dynmic shpe functions ecome equl to sttic shpe functions, nd sttic finite element is otined with ordinry convergence of numericl results. In cse of plte simply supported t two opposite edges, D prolem is reduced to D, i.e. strip element, with closed-form nlyticl solution. Dynmic stiffness mtrix sed on tht solution gives exct vlues of nturl frequencies for ny comintion of oundry condition t the remining two plte edges. The developed 4-noded dynmic finite elements re not relile for wider prcticl use. On the contrry, ppliction of the -noded strip dynmic finite element is very effective. In future reserch, the sme pproch cn e used for the development of -noded finite element for thick pltes, y employing the modified Mindlin theory [4]. 7. ACKNOWLOGEMENT This work ws supported y the Ntionl Reserch Foundtion of Kore (NRF) grnt funded y the Koren Government (MSIP) through GCRC-SOP (Grnt No ). ENGINEERING MODELLING 9 (06) -4, - 3
24 8. APPENDIX: DYNAMIC STIFFNESS AND FLEXIBILITY MATRIX OF THE PLATE SIMPLY SUPPORTED AT TWO OPPOSITE EDGES Elements of dynmic stiffness mtrix ccording to Eqs. (67), (68) nd (70): k = k = d α + d β αshαcosβ + βchα sinβ, 33 DK k = k = ( d α d β)( chαcosβ ) + ( d α + d β) shα sinβ, 34 DK 3 DK 4 3 DK 44 DK K k = d α + d β αshα + βsinβ, k = k = d α + d β chα cosβ, k = k = c + c αchα sinβ βshα cos β, D k = = c + c βshα α sinβ, 4 DK ( ) ( ) = αβ chα cos β + α β shα sinβ. Elements of dynmic flexiility mtrix ccording to Eqs. (67), (68) nd (7): l = l = c + c c d chα sinβ c d shαcosβ, 33 DL l = l = d d c c chαcosβ + c d + c d shα sinβ, ( ) 34 D L l = c + c c d shα c d sinβ, 3 DL l = l = d d c + c chα cosβ, 4 3 DL l = l = d α + d β c d chα sinβ + c d shαcosβ, 44 DL l = = d α + d β c d shα + c d sinβ, 4 DL ( ) ( ) D = c c d d chαcosβ + c d c d shα sinβ. L (A) (A) 9. REFERENCES [] O.C. Zienkiewicz, The Finite Element Method in Engineering Science, McGrw-Hill, New York, 97. [] J.S. Przemieniecki, Theory of Mtrix Structurl Anlysis, McGrw-Hill, New York, 968. [3] D.J. Dwe, Mtrix nd Finite Elements Displcement Anlysis of Structures, Clrendon Press, Oxford, 984. [4] R. Szilrd, Theory nd Applictions of Plte Anlysis, John Wiley & Sons, Hooken, ENGINEERING MODELLING 9 (06) -4, -
25 [] M. Boscolo, J.R. Bnerjee, Dynmic stiffness formultion for pltes using first order sher deformtion theory, st AIAA/ASME/ASCE/AHS/ASC Structures, Structurl Dynmics, nd Mterils Conference <BR> 8 th - April 00, Orlndo, Florid. [6] F.W. Willims, S. Yun, K. Ye, D. Kennedy, M.S. Djoudi, Towrds deep nd simple understnding of trnscendentl eigenprolem of structurl virtions, Journl of Sound nd Virtion, Vol. 6, No. 4, pp , 00. [7] J. Sorić, The Finite Element Method, Golden mrketing Tehničk knjig, Zgre, 003. (in Crotin). [8] I. Holnd, K. Bell, Finite Element Method in Stress Anlysis, Tpir-Trykk, 970. [9] I. Senjnović, The Finite Element Method in Ship Structure Anlysis, University of Zgre, Zgre, 998. (In Crotin). [0] Y. Xing, B. Liu, Closed form solutions for free virtions of rectngulr Mindlin pltes, Act Mechnic Sinic, Vol., No., pp , 009. [] I. Senjnović, M. Tomić, N. Vldimir, D.S. Cho, Anlyticl solution for free virtion of modertely thick rectngulr plte, Hindwi Pulishing Corportion, Mthemticl Prolems in Engineering, Volume 03, Article ID [] I. Senjnović, M. Tomić, N. Vldimir, An pproximte nlyticl procedure for nturl virtion nlysis of free rectngulr pltes, Thin-Wlled Structures, Vol. 9, pp. 0-4, 0. [3] K.M. Liew, Y. Xing, S. Kitipornchi, Trnsverse Virtion of Thick Rectngulr Pltes I. Comprehensive Set of Boundry Conditions, Computers & Structures, Vol. 49, No., pp. - 9, 993. [4] I. Senjnović, N. Vldimir, M. Tomić, An dvnced theory of modertely thick plte virtions, Journl of Sound nd Virtions, Vol. 33, pp , 03. ENGINEERING MODELLING 9 (06) -4, -
Review of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationPhysics 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 informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationNew 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 informationPlates 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 informationPlate Bending Analysis by using a Modified Plate Theory
Copyright c 2006 Tech Science Press CMES, vol.11, no.3, pp.103-110, 2006 Plte Bending Anlysis y using Modified Plte Theory Y. Suetke 1 Astrct: Since Reissner nd Mindlin proposed their clssicl thick plte
More information2. VECTORS AND MATRICES IN 3 DIMENSIONS
2 VECTORS AND MATRICES IN 3 DIMENSIONS 21 Extending the Theory of 2-dimensionl Vectors x A point in 3-dimensionl spce cn e represented y column vector of the form y z z-xis y-xis z x y x-xis Most of the
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationSection 4: Integration ECO4112F 2011
Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic
More informationA Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation
Americn Journl of Engineering Reserch (AJER) 13 Americn Journl of Engineering Reserch (AJER) e-issn : 3-847 p-issn : 3-936 Volume-, Issue-1, pp-388-393 www.jer.org Reserch Pper Open Access A Brief Note
More informationMAC-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 informationIntroduction to Finite Element Method
Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............
More informationKirchhoff and Mindlin Plates
Kirchhoff nd Mindlin Pltes A plte significntly longer in two directions compred with the third, nd it crries lod perpendiculr to tht plne. The theory for pltes cn be regrded s n extension of bem theory,
More informationPolynomials and Division Theory
Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the
More informationAn Alternative Approach to Estimating the Bounds of the Denominators of Egyptian Fractions
Leonrdo Journl of Sciences ISSN 58-0 Issue, Jnury-June 0 p. -0 An Alterntive Approch to Estimting the Bounds of the Denomintors of Egyptin Frctions School of Humn Life Sciences, University of Tsmni, Locked
More informationPatch Antennas. Chapter Resonant Cavity Analysis
Chpter 4 Ptch Antenns A ptch ntenn is low-profile ntenn consisting of metl lyer over dielectric sustrte nd ground plne. Typiclly, ptch ntenn is fed y microstrip trnsmission line, ut other feed lines such
More informationAdomian Decomposition Method with Green s. Function for Solving Twelfth-Order Boundary. Value Problems
Applied Mthemticl Sciences, Vol. 9, 25, no. 8, 353-368 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/.2988/ms.25.486 Adomin Decomposition Method with Green s Function for Solving Twelfth-Order Boundry
More informationThin plates with large deflection - Buckling instability and response
Chpter 2 Thin pltes with lrge deflection - Buckling instility nd response In the present chpter some issue on the theory of thin pltes with lrge deflections will e ddressed. This llows to study oth structurl
More informationExplain 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 informationSCHOOL OF ENGINEERING & BUILT ENVIRONMENT
SCHOOL OF ENGINEERING & BUIL ENVIRONMEN MARICES FOR ENGINEERING Dr Clum Mcdonld Contents Introduction Definitions Wht is mtri? Rows nd columns of mtri Order of mtri Element of mtri Equlity of mtrices Opertions
More informationApplication Chebyshev Polynomials for Determining the Eigenvalues of Sturm-Liouville Problem
Applied nd Computtionl Mthemtics 5; 4(5): 369-373 Pulished online Septemer, 5 (http://www.sciencepulishinggroup.com//cm) doi:.648/.cm.545.6 ISSN: 38-565 (Print); ISSN: 38-563 (Online) Appliction Cheyshev
More informationA finite thin circular beam element for out-of-plane vibration analysis of curved beams
Journl of Mechnicl Science nd echnology (009) 196~1405 Journl of Mechnicl Science nd echnology www.springerlink.com/content/178-494x DOI 10.1007/s106-008-11- A finite thin circulr bem element for out-of-plne
More informationLecture Solution of a System of Linear Equation
ChE Lecture Notes, Dept. of Chemicl Engineering, Univ. of TN, Knoville - D. Keffer, 5/9/98 (updted /) Lecture 8- - Solution of System of Liner Eqution 8. Why is it importnt to e le to solve system of liner
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationModule 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur
Module Anlysis of Stticlly Indeterminte Structures by the Mtrix Force Method Version CE IIT, Khrgpur esson 8 The Force Method of Anlysis: Bems Version CE IIT, Khrgpur Instructionl Objectives After reding
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
More informationMethods of Analysis and Selected Topics (dc)
Methods of Anlysis nd Selected Topics (dc) 8 Ojectives Become fmilir with the terminl chrcteristics of current source nd how to solve for the voltges nd currents of network using current sources nd/or
More informationCOSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature
COSC 336 Numericl Anlysis I Numericl Integrtion nd Dierentition III - Guss Qudrture nd Adptive Qudrture Edgr Griel Fll 5 COSC 336 Numericl Anlysis I Edgr Griel Summry o the lst lecture I For pproximting
More information10 Vector Integral Calculus
Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve
More informationLecture 4 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell
Lecture 4 Notes, Electromgnetic Theory I Dr. Christopher S. Bird University of Msschusetts Lowell 1. Orthogonl Functions nd Expnsions - In the intervl (, ) of the vrile x, set of rel or complex functions
More informationMethod of Localisation and Controlled Ejection of Swarms of Likely Charged Particles
Method of Loclistion nd Controlled Ejection of Swrms of Likely Chrged Prticles I. N. Tukev July 3, 17 Astrct This work considers Coulom forces cting on chrged point prticle locted etween the two coxil,
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationLINEAR ALGEBRA APPLIED
5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nth-order
More informationI1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3
2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is
More informationMATH34032: 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 informationNumerical Linear Algebra Assignment 008
Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com
More informationM344 - ADVANCED ENGINEERING MATHEMATICS
M3 - ADVANCED ENGINEERING MATHEMATICS Lecture 18: Lplce s Eqution, Anltic nd Numericl Solution Our emple of n elliptic prtil differentil eqution is Lplce s eqution, lso clled the Diffusion Eqution. If
More informationConducting Ellipsoid and Circular Disk
1 Problem Conducting Ellipsoid nd Circulr Disk Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 (September 1, 00) Show tht the surfce chrge density σ on conducting ellipsoid,
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More information13: 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 informationEdexcel GCE Core Mathematics (C2) Required Knowledge Information Sheet. Daniel Hammocks
Edexcel GCE Core Mthemtics (C) Required Knowledge Informtion Sheet C Formule Given in Mthemticl Formule nd Sttisticl Tles Booklet Cosine Rule o = + c c cosine (A) Binomil Series o ( + ) n = n + n 1 n 1
More information1 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 informationOrdinary 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 informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationSeptember 13 Homework Solutions
College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationCHAPTER 1 PROGRAM OF MATRICES
CHPTER PROGRM OF MTRICES -- INTRODUCTION definition of engineering is the science y which the properties of mtter nd sources of energy in nture re mde useful to mn. Thus n engineer will hve to study the
More informationThe practical version
Roerto s Notes on Integrl Clculus Chpter 4: Definite integrls nd the FTC Section 7 The Fundmentl Theorem of Clculus: The prcticl version Wht you need to know lredy: The theoreticl version of the FTC. Wht
More informationu( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.
nlyzing Dmped Oscilltions Prolem (Medor, exmple 2-18, pp 44-48): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $
More informationCHAPTER 4a. ROOTS OF EQUATIONS
CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment
More informationDETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE
Determintion RevAdvMterSci of mechnicl 0(009) -7 properties of nnostructures with complex crystl lttice using DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING
More informationMatrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24
Mtrix lger Mtrix ddition, Sclr Multipliction nd rnsposition Mtrix lger Section.. Mtrix ddition, Sclr Multipliction nd rnsposition rectngulr rry of numers is clled mtrix ( the plurl is mtrices ) nd the
More informationThings to Memorize: A Partial List. January 27, 2017
Things to Memorize: A Prtil List Jnury 27, 2017 Chpter 2 Vectors - Bsic Fcts A vector hs mgnitude (lso clled size/length/norm) nd direction. It does not hve fixed position, so the sme vector cn e moved
More information12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION
12. DYNAMIC ANALYSIS Force Equilibrium is Fundmentl in the Dynmic Anlysis of Structures 12.1 INTRODUCTION { XE "Newton's Second Lw" }All rel physicl structures behve dynmiclly when subjected to lods or
More informationMath 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 informationBend Forms of Circular Saws and Evaluation of their Mechanical Properties
ISSN 139 13 MATERIALS SCIENCE (MEDŽIAGOTYRA). Vol. 11, No. 1. 5 Bend Forms of Circulr s nd Evlution of their Mechnicl Properties Kristin UKVALBERGIENĖ, Jons VOBOLIS Deprtment of Mechnicl Wood Technology,
More informationSUPPLEMENTARY INFORMATION
DOI:.38/NMAT343 Hybrid Elstic olids Yun Li, Ying Wu, Ping heng, Zho-Qing Zhng* Deprtment of Physics, Hong Kong University of cience nd Technology Cler Wter By, Kowloon, Hong Kong, Chin E-mil: phzzhng@ust.hk
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationLinear Systems with Constant Coefficients
Liner Systems with Constnt Coefficients 4-3-05 Here is system of n differentil equtions in n unknowns: x x + + n x n, x x + + n x n, x n n x + + nn x n This is constnt coefficient liner homogeneous system
More informationSome basic concepts of fluid dynamics derived from ECE theory
Some sic concepts of fluid dynmics 363 Journl of Foundtions of Physics nd Chemistry, 2, vol. (4) 363 374 Some sic concepts of fluid dynmics derived from ECE theory M.W. Evns Alph Institute for Advnced
More informationChapter 6 Techniques of Integration
MA Techniques of Integrtion Asst.Prof.Dr.Suprnee Liswdi Chpter 6 Techniques of Integrtion Recll: Some importnt integrls tht we hve lernt so fr. Tle of Integrls n+ n d = + C n + e d = e + C ( n ) d = ln
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationModule 1. Energy Methods in Structural Analysis
Module 1 Energy Methods in Structurl Anlysis Lesson 4 Theorem of Lest Work Instructionl Objectives After reding this lesson, the reder will be ble to: 1. Stte nd prove theorem of Lest Work.. Anlyse stticlly
More informationdx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.
Chpter 8 Stility theory We discuss properties of solutions of first order two dimensionl system, nd stility theory for specil clss of liner systems. We denote the independent vrile y t in plce of x, nd
More information1 ELEMENTARY ALGEBRA and GEOMETRY READINESS DIAGNOSTIC TEST PRACTICE
ELEMENTARY ALGEBRA nd GEOMETRY READINESS DIAGNOSTIC TEST PRACTICE Directions: Study the exmples, work the prolems, then check your nswers t the end of ech topic. If you don t get the nswer given, check
More informationThe Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5
The Dirichlet Prolem in Two Dimensionl Rectngle Section 13.5 1 Dirichlet Prolem in Rectngle In these notes we will pply the method of seprtion of vriles to otin solutions to elliptic prolems in rectngle
More informationA - INTRODUCTION AND OVERVIEW
MMJ5 COMPUTATIONAL METHOD IN SOLID MECHANICS A - INTRODUCTION AND OVERVIEW INTRODUCTION AND OVERVIEW M.N. Tmin, CSMLb, UTM MMJ5 COMPUTATIONAL METHOD IN SOLID MECHANICS Course Content: A INTRODUCTION AND
More informationThe 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 informationJack 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 informationCreating A New Planck s Formula of Spectral Density of Black-body Radiation by Means of AF(A) Diagram
nd Jogj Interntionl Physics Conference Enhncing Network nd Collortion Developing Reserch nd Eduction in Physics nd Nucler Energy Septemer 6-9, 007, Yogykrt-Indonesi Creting A New Plnck s Formul of Spectrl
More informationPlate Theory. Section 11: PLATE BENDING ELEMENTS
Section : PLATE BENDING ELEMENTS Plte Theor A plte is structurl element hose mid surfce lies in flt plne. The dimension in the direction norml to the plne is referred to s the thickness of the plte. A
More informationFarey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University
U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions
More informationStatically indeterminate examples - axial loaded members, rod in torsion, members in bending
Elsticity nd Plsticity Stticlly indeterminte exmples - xil loded memers, rod in torsion, memers in ending Deprtment of Structurl Mechnics Fculty of Civil Engineering, VSB - Technicl University Ostrv 1
More informationINEQUALITIES FOR TWO SPECIFIC CLASSES OF FUNCTIONS USING CHEBYSHEV FUNCTIONAL. Mohammad Masjed-Jamei
Fculty of Sciences nd Mthemtics University of Niš Seri Aville t: http://www.pmf.ni.c.rs/filomt Filomt 25:4 20) 53 63 DOI: 0.2298/FIL0453M INEQUALITIES FOR TWO SPECIFIC CLASSES OF FUNCTIONS USING CHEBYSHEV
More informationDerivations for maximum likelihood estimation of particle size distribution using in situ video imaging
2 TWMCC Texs-Wisconsin Modeling nd Control Consortium 1 Technicl report numer 27-1 Derivtions for mximum likelihood estimtion of prticle size distriution using in situ video imging Pul A. Lrsen nd Jmes
More informationLinear Inequalities. Work Sheet 1
Work Sheet 1 Liner Inequlities Rent--Hep, cr rentl compny,chrges $ 15 per week plus $ 0.0 per mile to rent one of their crs. Suppose you re limited y how much money you cn spend for the week : You cn spend
More informationShear and torsion interaction of hollow core slabs
Competitive nd Sustinble Growth Contrct Nº G6RD-CT--6 Sher nd torsion interction of hollow core slbs HOLCOTORS Technicl Report, Rev. Anlyses of hollow core floors December The content of the present publiction
More informationGeneralizations of the Basic Functional
3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............
More information5: The Definite Integral
5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description
More informationThe Algebra (al-jabr) of Matrices
Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense
More informationOPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS
IJRET: Interntionl Journl of Reserch in Engineering nd Technology eissn: 9-6 pissn: -78 OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS T. Rmchndrn R.Priml
More informationCBE 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 informationThomas Whitham Sixth Form
Thoms Whithm Sith Form Pure Mthemtics Unit C Alger Trigonometry Geometry Clculus Vectors Trigonometry Compound ngle formule sin sin cos cos Pge A B sin Acos B cos Asin B A B sin Acos B cos Asin B A B cos
More informationChapter 5 Bending Moments and Shear Force Diagrams for Beams
Chpter 5 ending Moments nd Sher Force Digrms for ems n ddition to illy loded brs/rods (e.g. truss) nd torsionl shfts, the structurl members my eperience some lods perpendiculr to the is of the bem nd will
More informationTHEORY OF VIBRATIONS OF TETRA-ATOMIC SYMMETRIC BENT MOLECULES
terils Physics nd echnics (0 9- Received: rch 0 THEORY OF VIRTIONS OF TETR-TOIC SYETRIC ENT OLECLES lexnder I eler * ri Krupin Vitly Kotov Deprtment of Physics of Strength nd Plsticity of terils Deprtment
More informationVibration Analysis for Rectangular Plate Having a Circular Central Hole with Point Support by Rayleigh-Ritz Method
Journl of Solid Mechnics Vol. 6, No. 1 (014) pp. 8-4 Vibrtion Anlysis for ectngulr Plte Hving Circulr Centrl Hole with Point Support by yleigh-itz Method K. Torbi *, A.. Azdi Deprtment of Mechnicl Engineering,
More informationPLEASE SCROLL DOWN FOR ARTICLE
This rticle ws downloded y:[knt, Trun] On: 8 Jnury 2008 Access Detils: [suscription numer 789308335] Pulisher: Tylor & Frncis Inform Ltd Registered in Englnd nd Wles Registered Numer: 1072954 Registered
More informationDEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.
398 CHAPTER 11 ORTHOGONAL FUNCTIONS AND FOURIER SERIES 11.1 ORTHOGONAL FUNCTIONS REVIEW MATERIAL The notions of generlized vectors nd vector spces cn e found in ny liner lger text. INTRODUCTION The concepts
More informationPlate Theory. Section 13: PLATE BENDING ELEMENTS
Section : PLATE BENDING ELEENTS Wshkeic College of Engineering Plte Theor A plte is structurl element hose mid surfce lies in flt plne. The dimension in the direction norml to the plne is referred to s
More informationGreen function and Eigenfunctions
Green function nd Eigenfunctions Let L e regulr Sturm-Liouville opertor on n intervl (, ) together with regulr oundry conditions. We denote y, φ ( n, x ) the eigenvlues nd corresponding normlized eigenfunctions
More informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.
Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil
More informationLecture 2: January 27
CS 684: Algorithmic Gme Theory Spring 217 Lecturer: Év Trdos Lecture 2: Jnury 27 Scrie: Alert Julius Liu 2.1 Logistics Scrie notes must e sumitted within 24 hours of the corresponding lecture for full
More informationFinite Element Determination of Critical Zones in Composite Structures
Finite Element Determintion of Criticl Zones in Composite Structures Alexey I. Borovkov Dmitriy V. Klimshin Denis V. Shevchenko Computtionl Mechnics Lb., St. Petersburg Stte Polytechnicl University, Russi
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More information