In-Plane Responses of Axially Moving Plates Subjected to Arbitrary Edge Excitations

Size: px
Start display at page:

Download "In-Plane Responses of Axially Moving Plates Subjected to Arbitrary Edge Excitations"

Transcription

1 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 In-Plane Responses of Axially oing Plates Sjected to Aritrary Edge Excitations T. H. Yong and Y. S. Cio International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 Astract The free and forced in-plane irations of axially moing plates are inestigated in this work. The plate possesses an internal damping of which the constittie relation oeys the Kelin-Voigt model, and the excitations are aritrarily distrited on two opposite edges. First, the eqations of motion and the ondary conditions of the axially moing plate are deried. Then, the extended Ritz method is sed to otain discretized system eqations. Finally, nmerical reslts for the natral freqencies and the mode shapes of the in-plane iration and the in-plane response of the moing plate sjected to aritrary edge excitations are presented. It is osered that the symmetry class of the mode shapes of the in-plane iration disperses gradally as the moing speed gets higher, and the - and -components of the mode shapes elong to different symmetry class. In addition, large response amplitdes haing shapes similar to the mode shapes of the plate can e excited y the edge excitations at the resonant freqencies and with the same symmetry class of distrition as the -components. Keywords Aritrary edge excitations, axially moing plates, in-plane iration, extended Ritz method. A I. ITRODUCTIO XIALLY moing wes are freqently encontered applications in indstries, sch as magnetic tapes, driing elts, paper sheets, and plastic wes, etc. In order to raise prodction efficiency, the moing speed of the we gets faster and faster, which worsens the iration prolem of the we. Therefore, a ast nmer of researches concerning the dynamic ehaiors of moing wes can e fond in the literatre [1]-[6]. Becase the lowest natral freqency and critical speed of the ot-of-plane iration of the moing we are lower than those of the in-plane iration, most of the references in literatre deal with the ot-of-plane iration of the moing wes. In literatre, the references concerning the in-plane iration of rectanglar plates are mch lesser than those pertaining to the ot-of-plane iration of rectanglar plates. Hyde et al. [7] stdied the free in-plane iration of rectanglar plates ndergoing plane stress deformation throgh the Ritz discretization. The natral freqencies and mode shapes of rectanglar plates with three different kinds of ondary conditions were presented in this work. The free in-plane irations of a rectanglar plate nder comined static pressre/tension edge loads were inestigated y Waer [8]. In this paper, the in-plane deformation de to static edge loads is T. H. Yong and Y. S. Cio are with ational Taiwan Uniersity of Science and Technology, Taipei 16, Taiwan (phone: ; fax: ; ( thyong@mail.ntst.ed.tw, 131@ mail.ntst.ed. tw). determined first y soling the plane stress prolem, and the free in-pane iration is analyzed y the Rayleigh-Ritz method. A mathematical model was deeloped y Farag and Pan [9] for the prediction of the forced response of rectanglar plates with all edges clamped to in-plane point force excitations. All the aoe three references dealt with the in-plane iration of a stationary rectanglar plate, and the references dealing with the in-plane iration of axially moing rectanglar plates are een less. Shin et al. [1] stdied the free in-plane iration of an axially moing plate y the Galerkin method. The elastic plate is sjected to a constant tension along two mass-transport ondaries, and two sets of ondary conditions, which are free and fixed constraints in the lateral direction at two mass-transport ondaries are discssed in this stdy. To athors knowledge, no research work on the in-plane response of axially moing plates sjected to aritrary edge excitations can e fond in the literatre. Therefore, this work inestigates the free and forced in-plane irations of axially moing plates. The plate possesses an internal damping of which the constittie relation oeys the Kelin-Voigt model, and the excitations are aritrarily distrited on two opposite edges. First, the eqations of motion and the ondary conditions of the axially moing plate are deried. Then, the extended Ritz method is sed to otain discretized system eqations. Finally, nmerical reslts for the natral freqencies and the mode shapes of the in-plane iration and the in-plane response of the moing plate sjected to aritrary edge excitations are presented. II. EQUATIOS OF OTIO Consider a rectanglar plate of dimension a moing in the x-direction with a constant speed V. The plate is simply spported on two opposite edges x = and a, and is free on the other two edges y = and. The plate is sjected to an excitation f ( y, t) aritrarily distrited on two simply spported edges, as shown in Fig. 1. The plate possesses an internal damping of which the constittie relation oeys the Kelin-Voigt model, i.e., ε σ = E( ε + c ) (1) where σ and ε are the stress and strain of the plate, respectiely, and t is a temporal ariale. E and c are Yong s modls and internal damping coefficient of the plate, International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

2 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 respectiely. Therefore, the stress-strain relations for the plane stress state are gien y E ε ε x y σ x = [( ε ) c( )] x + νε y + + ν 1 ν (a) E ε y ε x σ y = [( ε ) c( )] y + νε x + + ν 1 ν () E ε xy σ xy = [ ε xy + c ] (c) 1+ ν where ν is Poisson s ratio of the plate. In terms of the in-plane displacement components, the strain energy of the plate can e expressed as Eh U = (1 ν ) a [( ) + ( ) 1 ν + ν( )( ) + ( + ) ] dxdy where h is the thickness of the plate, and and are the in-plane displacement components of the plate in the x- and y-directions, respectiely. The kinetic energy of the axially moing plate is gien y ρh T = a [( V + + V ) + ( + V ) ] dxdy where ρ is the mass density of the plate. The irtal work done y the edge excitation is written as δw K = x= a δ x= ] (3) (4) f ( y, t)[ δ dy (5) The irtal work done y the internal damping force is gien y δw d Ehc = 1 ν a [( + ν ) δ( ) + ( + 1 ν ν ) δ( ) + ( + ) δ( + )] dxdy The irtal work done y the momentm of the plate moing across the ondaries is fond to e δ = (6) a a ρhv[( V + + V ) δ + ( + V ) δ ] dy (7) The extended Hamilton principle for non-conseratie systems reqires t t1 ( δt δu + δwk + δwd δ ) dt = (8) where t 1 and t are two time instants at which the irtal displacements δ and δ anish. Sstitting (3) (7) into (8) and going throgh ariational calcls maniplations yield the eqations of motion and ondary conditions for the axially moing plate, ρh( + V + V Eh ) 1 ν Eh Ehc ( + ) (1 + ν) x y 1 ν ( + ν ) ( + ν ) Ehc ( + ) = f ( y, t)[ δ( x a) δ( x)] (1 + ν) (9a) ρh( + V + V ) Eh ( + ν ) 1 ν where Eh Ehc ( + ) (1 + ν ) x y 1 ν Ehc ( + ) = (1 + ν ) x, ( + ν ) (9) at x = and a, = and = (1a) x at y = and, = and = (1) y, and y xy are the in-plane stress resltants of the plate. Eqations (1a) and (1) can e redced to those otained y Shin et al. [1] if the internal damping force terms are deleted. Eqations (1a) and (1) constitte a system of the second order partial differential eqations in terms of the displacement components and. Since the exact soltion of them is not feasile, to otain approximate soltions in the finite dimensional space, discretization of the system eqations is carried ot. De to the ondary conditions of the prolem, it is ery difficlt to find the comparison fnctions, the ones that satisfy oth the geometric and natral ondary conditions, as the trial fnctions. Therefore, the extended Ritz method, which admits the admissile fnctions satisfying the geometric ondary conditions only as the trial fnctions, is adopted in conjnction with the extended Hamilton principle. Assme the displacement components and to e of the forms, x m y n ( x, y, t) = A m n mn ( t)( ) ( ), = = a (11a) ( m + 1) πx y n ( x, y, t) = B m n mn ( t) sin ( ) = = a (11) where A mn (t) and B mn (t) are ndetermined fnctions. Sstitting (11a) and (11) into (8) and non- dimensionalizing all parameters and ariales yield the following eqation, xy International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

3 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 [ ] U + ( V [ G] + c[ K ]) U + ([ K ] + V [ K ]) U F( τ ) (1) e e g = (1 ν ρa where the non-dimensional parameters V = V, E E E c = c, τ = t, and a prime denotes (1 ν ) ρa ( 1 ν ) ρa a differentiation with respect to τ. The mass matrix [ ], the moing-speed dependent matrix [G], the elastic stiffness matrix [ K e ], and the geometric stiffness matrix [ K g ] hae the following forms, [ ] [] [] [ ] =, [ ] = [ [ ] = ] [ G ] [] G, [] [ G ] [ K ] [ K ] [ K g ] [] K e, [ K g ] =, [ K ] [ K ] [] [ K g ] in which the lock matrices [ ], [ ], [ K ], and [ K ] T are symmetric, and [ K ] = [ K ]. Therefore, the mass matrix [ ] and the elastic stiffness matrix K ] are ) [ e symmetric, the moing-speed dependent matrix [G] and the geometric stiffness matrix [ K g ] are neither symmetric nor skew-symmetric. The response ector U and the forcing ector F (τ) in (1) hae the following forms, A f ( τ) U =, F ( τ) = B where A and B are colmn matrices formed y all and B mn (τ) /, and f (τ) is a colmn matrix. A mn (τ) / Eqation (1) goerns the in-plane forced response of an axially moing plate. To find the natral freqencies and the mode shapes of the free in-plane iration, neglect the nonhomogeneos forcing term on the right-hand side of the eqation, and rewritten the set of the second order differential eqations into a set of first order differential eqations, which yield the following eigenale prolem, s [ ] [ ] [ ] [ I ] ( V Z + [ G] + c [ K ]) ([ K ] + V [ K ]) e e g Z = (13) [ I ] [ ] where [I ] is an identity matrix. The eigenales of the aoe eqation are in pairs of complex conjgates, s j 1, j = λ j ± iω j, j=1,,, L, where L = ( + 1)( + 1) is the total degrees of freedom of the discretized system. The eigenectors of the aoe eqation are also in pairs of complex conjgates, Z j 1, j = X j ± iy j, j=1,,, L. The imaginary part of the eigenale ω j is the non- dimensionalized natral freqency of the jth mode of the axially moing plate, while ς = λ / ω is the damping ratio of the jth mode of the axially j j j moing plate. If the real parts of all eigenales are smaller than or eqal to, the axially moing plate is stale; otherwise, the axially moing plate is nstale. Therefore, the critical speeds for the in-plane iration of the plate can e otained when the real part of one eigenale ecomes positie. To calclate the in-plane forced response of the axially moing plate, assme the edge excitation is simply harmonic in time, i.e., f ( y, t) = q( y) sin Ωt, where Ω is the exciting freqency. The forcing ector F (τ) in (1) will hae the form F ( τ) = Q sin Ωτ, where Q is a constant ector, and Ω = Ω (1 ν ) ρa / E. To match the form of the forcing ector, assme the response ector to e of the form, U τ) = U s sin Ωτ + U cosωτ (14) ( c where U s and U c are the sine and cosine components of the response ector, respectiely. Sstitting (14) into (1) yields the following eqation, ( Q = [ K ] + [ ] Ω [ ] Ω [ ] + e V K g ) ( V G c[ Ke ]) Ω( [ ] + V G c[ K ]) ([ K ] + V [ K ] Ω [ ] e e g U ) U s c (15) Once U s and U s are soled, the fnctions A ) and B mn (τ) are determined, from which the in-plane responses and are calclated. Since the response ector contains the sine and cosine components, the fnctions A ) and B ) mn (τ consist of the sine and cosine components also, i.e., Amn ( τ ) s c = A sin Ωτ + cosωτ mn Amn Bmn ( τ ) s c = B sin Ωτ + cosωτ mn Bmn mn (τ mn (τ (16a) (16) s c s c where A mn, A mn, B mn, and B mn are constants. The responses and can e written in the amplitde-phase angle form, where the amplitdes ( ξ, η, τ) = ( ξ, η)sin( Ωτ φ ) (17a) C ( ξ, η, τ) = ( ξ, η)sin( Ωτ φ ) (17) C International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

4 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 C ( ξ, η) = C ( ξ, η) = [ [ and the phase angles φ φ = tan = tan 1 A. Free Viration [ [ m= n= m= n= m= n= m= n= 1 B B s mn c mn A A m= n= m= n= m= n= m= n= s mn c mn x ( ) a x ( ) a y n ( ) ] y ( ) ] m m n + ( m + 1) πx y n sin ( ) ] a ( m + 1) πx y n sin ( ) ] a s x m y n Amn ( ) ( ) a c x m y n Amn ( ) ( ) a s ( m + 1) πx y Bmn sin ( ) a c ( m + 1) πx y Bmn sin ( ) a III. UERICAL RESULTS In order to erify the correctness of the reslts otained in this work, comparison is carried ot with the reslts y Shin et al. [1]. Tale I presents the lowest six natral freqencies of an axially moing rectanglar plate withot internal damping. It is osered that with ==8, the nmerical reslts otained in this work agree excellently with those reported in the reference [1]. The maximm difference etween these two sets of reslts sing the same degrees of freedom is aot.5%. Fig. illstrates the lowest six mode shapes of an axially moing rectanglar plate with internal damping. The -component and the -component of each mode shape are drawn in different plots, respectiely. Since the moing speed of the plate is ery low, the mode shapes almost presere the symmetry class that is owned only y stationary plates. The figre shows that the -components of the first and the forth mode shapes are anti-symmetric in oth the x- and y-directions. The -components of the second and the fifth mode shapes are symmetric in the x-direction and are anti- symmetric in the y-direction. The -component of the third mode shape is anti-symmetric in the x-direction and is symmetric in the y-direction, while that of the sixth mode shape is symmetric in oth the x- and y-directions. The figre shows the opposite trend for the -components. oreoer, the fifth mode shape of the axially moing rectanglar plate is a prely -component mode. The symmetry class of the mode shapes disperses gradally as the moing speed of the plate ecomes higher ecase the moing speed dependent matrix and the geometric stiffness matrix are neither symmetric nor skew-symmetric. Tale II shows the lowest six natral freqencies and modal damping ratios of an axially moing rectanglar plate with internal damping. It is fond that the natral freqencies of the n n + plate with internal damping are slightly lower than those of the plate withot internal damping, and the modal damping ratio is monotonically increasing for higher modes. Fig. 3 depicts the real and the imaginary parts of the eigenales of the lowest seen modes of an axially moing rectanglar plate with internal damping. In Fig. 3 (a), all the imaginary parts of the eigenales, which are the natral freqencies of the axially moing plate, drop as the moing speeds of the plate increase. At V =.486, the imaginary part of the lowest mode redces to zero. oreoer, eering phenomena occr when the imaginary part of the seenth mode meets that of the sixth mode at V =.1587; when the imaginary part of the sixth mode meets that of the fifth mode at V =.1637, and when the imaginary part of the sixth mode meets that of the seenth mode at V = In Fig. 3 (), all the real parts of the eigenales are negatie and hold nearly constant except the followings. The real part of the lowest mode rises and ecomes positie at the speed which the imaginary part of the lowest mode redces to zero. That is the moing plate ecomes nstale at this speed, and this speed is called the first critical speed of the moing plate. The real parts of the fifth, sixth and seenth modes jmp p and down respectiely, and cross each other at the speeds where eering phenomena occr. Becase the internal damping of the plate in this figre is small, Fig. 3 (a) looks ery similar to that otained y Shin et al. [1] in which no internal damping is considered. B. Forced Viration Fig. 4 illstrates the forced response of an axially moing plate sjected to excitations niformly distrited on two simply spported edges at zero and the lowest fie resonant freqencies. All -components of the response plots in the figre are nearly anti-symmetric in the axial direction and are nearly symmetric in the lateral direction, whereas all -components of the response plots show the opposite symmetry class, that is, nearly symmetric in the axial direction and nearly anti-symmetric in the lateral direction. This trend loses gradally for higher resonant freqencies. Becase the response plots possess the same symmetry class of the third mode, the amplitdes of the response plot at the third resonant freqency are one-order larger than those at the other resonant freqencies. Except for that, the response amplitdes generally redce gradally as the resonant freqencies get higher, and the static deformation of the plate sjected to static edge loads are close to the forced response amplitdes at the first resonant freqency. Fig. 5 shows the forced response of an axially moing plate sjected to excitations linearly distrited on two simply spported edges sch that the distrition of the excitations are anti-symmetric in the lateral direction. All -components of the response plots in the figre are almost anti-symmetric in oth axial and lateral directions, whereas all -components of the response plots are almost symmetric in oth directions. Again, this trend loses gradally for higher resonant freqencies. Becase the response plots possess the same symmetry class of International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

5 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 the first and the forth modes, the amplitdes of the response plot at the first resonant freqency are at least two-order larger than those at the other resonant freqencies, and the amplitdes of the response plot at the forth resonant freqency are larger than those at the second and the third resonant freqencies. The response plot at the first resonant freqency is similar to the static deformation plot sjected to static edge loads, and the response amplitdes at the first resonant freqency are aot one hndred times of the static deformation sjected to static edge loads, which is the inerse of the internal damping factor c. ote that the amplitdes of the response plot at the first resonant freqency in this figre are een larger than those of the plate sjected to niformly distrited edge excitations at the third resonant freqency, as shown in Fig. 3 (c). The forced response of an axially moing plate sjected to excitations linearly, anti-symmetrically distrited on the middle hales of two simply spported edges is depicted in Fig. 6. The distrition of the excitations is again anti- symmetric in the lateral direction, and the intensity of the excitation is twice as that in Fig. 5. All corresponding plots in oth figres look similar, t the amplitdes in this figre are a little it smaller than those in the preios figre. The effect of the internal damping on the forced response of an axially moing plate is shown in Fig. 7. In this figre, the forced response of the plate considered in Fig. 5 t with a larger internal damping factor is illstrated. Comparing etween oth figres, it is osered that the response plots in oth figres look similar, and the response plots at the first, the second, and the forth resonant freqencies in this figre deiate more from the mode shapes of the plate than the corresponding plots in Fig. 5. Becase the response plots possess the same symmetry class of the first and the forth modes, the amplitde ratios etween Figs. 5 and 7 at the first and the forth resonant freqencies are aot 1:1, which is approximately inersely proportional to the internal damping ratio etween these two figres. Bt the amplitdes etween these two figres at the other resonant freqencies are of the same orders. The effect of the aspect ratio on the forced response of an axially moing plate is presented in Fig. 8. In this figre, the forced response of a sqare plate is considered. Becase the aspect ratio of the plate is smaller, the natral freqencies of the plate rise and the mode shapes of the plate change also. The -component of the first mode shape is anti-symmetric in oth the x- and y-directions. The -components of the second and the forth mode shapes are symmetric in the x-direction and are anti-symmetric in the y-direction. The -components of the third and the fifth mode shapes are anti-symmetric in the x-direction and are symmetric in the y-direction. The -components of the mode shapes hae the opposite class of symmetry. Becase the distrition of the edge excitation is linear, anti-symmetric in the lateral direction, the amplitdes of the response at the first resonant freqency is the highest among those at the lowest fie resonant freqencies. In addition, the amplitdes of the response at the first resonant freqency in this figre are mch smaller than those at the first resonant freqency in Fig. 5 (a) ecase the aspect ratio of the plate in this figre is smaller than that of the plate considered in Fig. 5, and hence the plate considered in this figre is more stiffer. IV. COCLUSIOS The free and forced in-plane irations of axially moing plates sjected to aritrarily distrited edge excitations on two opposite, simply-spported edges are stdied y the extended Ritz method in this work. merical reslts for the natral freqencies, mode shapes and the critical speed of the free in-plane iration and the response amplitdes of the forced in-plane iration are presented for different system parameters and the distrition of the edge excitations. The nmerical reslts for the natral freqencies, mode shapes and the critical speed of the free in-plane iration are compared with [1], and excellent agreement of the nmerical reslts etween these two works is osered. From these nmerical reslts, the following conclsions can e drawn: 1) The mode shapes of the corresponding stationary plate presere symmetry class, and the -components and -components of the mode shapes elong to the opposite symmetry class. The symmetry class of the mode shapes of the free in-plane iration disperses gradally as the moing speed of the plate ecomes higher. ) The natral freqencies of the plate with internal damping are slightly lower than those of the plate withot internal damping, and the modal damping ratio is monotonically increasing for higher modes. 3) At the first critical speed, the imaginary part of the lowest mode redces to zero, and the real part of the lowest mode ecomes positie. Veering phenomena occr among the fifth, sixth and seenth modes when two imaginary parts of these modes meet each other, and the real parts of these modes jmp p and down respectiely and cross each other at the speeds where eering phenomena occr. 4) For an axially moing plate sjected to excitations symmetrically distrited in the lateral direction on two simply spported edges, all -components of the response plots are nearly anti-symmetric in the axial direction and are nearly symmetric in the lateral direction. For an axially moing plate sjected to excitations anti-symmetrically distrited in the lateral direction on two simply spported edges, all -components of the response plots are nearly anti-symmetric in oth the axial and the lateral directions. The -components of the response plots show the opposite symmetry class. This trend loses gradally for higher resonant freqencies and faster moing speeds. 5) If the response plot possesses the same symmetry class as the mode shape of a mode, the amplitdes of the response plot at that resonant freqency will e mch larger than those at the other resonant freqencies. Except for that, the response amplitdes generally redce gradally as the resonant freqencies get higher. 6) The effect of the internal damping is ery effectie only when the symmetry class of the edge excitation meets the International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

6 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 same symmetry class of the mode shape at that resonant freqency. Otherwise, its effect is insignificant. TABLE I THE ATURAL FREQUECIES OF A UDAPED AXIALLY OVIG PLATE ode no. Present work a / =3, V =.1, c = [1] Error %.131 %.34 %.54 % % % International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 TABLE II THE ATURAL FREQUECIES AD ODAL DAPIG RATIOS OF A DAPED AXIALLY OVIG PLATE a / =3, V =.1, c =.1 ode no ω j ς j Fig. 1 An axially moing plate sjected to aritrarily distrited excitations () nd mode (c) 3 rd mode (d) 4 th mode (a) 1 st mode International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

7 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (e) 5 th mode (f) 6 th mode Fig. The mode shapes of a damped axially moing plate a / =3,V =.1, c =.1 (a) Imaginary part () Real part Fig. 3 The eigenales cres of the lowest seen modes of an axially moing plate a / =3, c =.1 (a) Ω / ω1 = () Ω / ω1 =.9999 International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

8 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (c) Ω / ω1 =3.749 (d) Ω / ω1 =3.837 (f) Ω / ω1 = Fig. 4 The forced response plots of an axially moing plate sjected to niformly distrited excitations (a) Ω / ω1 = a / =3, V =.1, c =.1 (e) Ω / ω1 =5.418 () Ω / ω1 =.9999 International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

9 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (c) Ω / ω1 =3.749 (d) Ω / ω1 =3.837 (e) Ω / ω1 =5.418 (f) Ω / ω1 = Fig. 5 The forced response plots of an axially moing plate sjected to linearly, anti-symmetrically distrited excitations a / =3,V =.1, c =.1 (a) Ω / ω1 = () Ω / ω1 =.9999 International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

10 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (c) Ω / ω1 =3.749 (d) Ω / ω1 =3.837 (e) Ω / ω1 =5.418 (f) Ω / ω1 = Fig. 6 The forced response plots of an axially moing plate sjected to excitations linearly, anti-symmetrically distrited on the middle hales of the edges a / =3,V =.1, c =.1 (a) Ω / ω1 = () Ω / ω1 =.999 International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

11 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (c) Ω / ω1 =3.59 (d) Ω / ω1 = (f) Ω / ω1 =6.865 Fig 7 The forced response plots of an axially moing plate sjected to linearly, anti-symmetrically distrited excitations a / =3,V =.1, c =.1 (a) Ω / ω1 = () Ω / ω1 =1. (e) Ω / ω1 =5.899 (c) Ω / ω1 =1.31 International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

12 International Jornal of echanical and echatronics Engineering Vol:7, o:7, 13 [9]. H. Farag, and J. Pan, Free and forced in-plane iration of rectanglar plates, Jornal of Acostical Society of America, ol. 13, pp , [1] C. Shin, W. Kim, and J. Chng, Free in-plane iration of an axially moing memrane, Jornal of Sond and Viration, ol. 7, pp , 4. (d) Ω / ω1 = International Science Index, echanical and echatronics Engineering Vol:7, o:7, 13 waset.org/plication/16414 (e) Ω / ω1 =.379 (f) Ω / ω1 =.391 Fig. 8 The forced response plots of an axially moing plate sjected to linearly, anti-symmetrically distrited excitations a / =1,V =.1, c =.1 REFERECES [1] F. Fng, J. S. Hang, and Y. C. Chen, The transient amplitde of the iscoelastic traelling string: An integral constittie law, Jornal of Sond and Viration, ol. 1, pp , [] L. Q. Chen, and X. D. Yang, Steady-state response of axially moing iscoelastic eams with plsating speed: comparison of two nonlinear models, International Jornal of Solids and Strctres, ol. 4, pp. 35-5, 5. [3] C. Shin, J. Chng, and W. Kim, Dynamic characteristics of the ot-of-plane irations for an axially moing memrane, Jornal of Sond and Viration, ol. 97, pp , 6. [4] K. arynowski, Two-dimensional rheological element in modeling of axially moing iscoelastic we, Eropean Jornal of echanics A/Solids, ol. 5, pp , 6. [5] C. C. Lin, Staility and iration characteristics of axially moing plates, International Jornal of Solids and Strctres, ol. 34, pp , [6] S. Hatami, H. R. Ronagh, and. Azhari, Exact free iration analysis of axially moing iscoelastic plates, Compters & Strctres, ol. 86, pp , 8. [7] K. Hyde, J. Y. Chang, C. Bacca, and J. A. Wickert, Parameter stdies for plane stress in-plane iration of rectanglar plates, Jornal of Sond and Viration, ol. 47, pp , 1. [8] J. Waer, In-plane irations of rectanglar plates nder inhomogeneos static edge load, Eropean Jornal of echanics A/Solids, ol. 11, pp , 199. International Scholarly and Scientific Research & Innoation 7(7) scholar.waset.org/ /16414

CHAPTER 8 ROTORS MOUNTED ON FLEXIBLE BEARINGS

CHAPTER 8 ROTORS MOUNTED ON FLEXIBLE BEARINGS CHAPTER 8 ROTORS MOUNTED ON FLEXIBLE BEARINGS Bearings commonly sed in heavy rotating machine play a significant role in the dynamic ehavior of rotors. Of particlar interest are the hydrodynamic earings,

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b)

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b) . Sole Problem.-(c). Sole Problem.-(b). A two dimensional trss shown in the figre is made of alminm with Yong s modls E = 8 GPa and failre stress Y = 5 MPa. Determine the minimm cross-sectional area of

More information

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL BBAA VI International Colloqim on: Blff Bodies Aerodynamics & Applications Milano, Italy, Jly, 2-24 28 FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

ON THE PERFORMANCE OF LOW

ON THE PERFORMANCE OF LOW Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

NON-LINEAR DYNAMICS OF CABLE-STAYED MASTS

NON-LINEAR DYNAMICS OF CABLE-STAYED MASTS Copyright 29 by ABCM Janary 4-8, 21, Foz do Igaç, R, Brazil NON-LINEAR DYNAMICS OF CABLE-STAYED MASTS Zenon J. G. N. del rado, zenon@eec.fg.br School of Ciil Engineering, Federal Uniersity of Goiás, UFG.

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

Xihe Li, Ligong Wang and Shangyuan Zhang

Xihe Li, Ligong Wang and Shangyuan Zhang Indian J. Pre Appl. Math., 49(1): 113-127, March 2018 c Indian National Science Academy DOI: 10.1007/s13226-018-0257-8 THE SIGNLESS LAPLACIAN SPECTRAL RADIUS OF SOME STRONGLY CONNECTED DIGRAPHS 1 Xihe

More information

Evolution of N-gram Frequencies under Duplication and Substitution Mutations

Evolution of N-gram Frequencies under Duplication and Substitution Mutations 208 IEEE International Symposim on Information Theory (ISIT) Evoltion of N-gram Freqencies nder Dplication and Sstittion Mtations Hao Lo Electrical and Compter Engineering University of Virginia Charlottesville,

More information

A Model-Free Adaptive Control of Pulsed GTAW

A Model-Free Adaptive Control of Pulsed GTAW A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,

More information

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qalitatie featres of a discrete dynamical system of homogeneos difference

More information

The Open Civil Engineering Journal

The Open Civil Engineering Journal Send Orders for Reprints to reprints@benthamscience.ae 564 The Open Ciil Engineering Jornal, 16, 1, 564-57 The Open Ciil Engineering Jornal Content list aailable at: www.benthamopen.com/tociej/ DOI: 1.174/187414951611564

More information

We automate the bivariate change-of-variables technique for bivariate continuous random variables with

We automate the bivariate change-of-variables technique for bivariate continuous random variables with INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,

More information

The spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam

The spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College

More information

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently Relatiity II I. Henri Poincare's Relatiity Principle In the late 1800's, Henri Poincare proposed that the principle of Galilean relatiity be expanded to inclde all physical phenomena and not jst mechanics.

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

Experiment and mathematical model for the heat transfer in water around 4 C

Experiment and mathematical model for the heat transfer in water around 4 C Eropean Jornal of Physics PAPER Experiment and mathematical model for the heat transfer in water arond 4 C To cite this article: Naohisa Ogawa and Fmitoshi Kaneko 2017 Er. J. Phys. 38 025102 View the article

More information

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION APPENDIX B MAIX NOAION he Deinition o Matrix Notation is the Deinition o Matrix Mltiplication B. INODUCION { XE "Matrix Mltiplication" }{ XE "Matrix Notation" }he se o matrix notations is not necessary

More information

MEG 741 Energy and Variational Methods in Mechanics I

MEG 741 Energy and Variational Methods in Mechanics I MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS 7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,

More information

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element

More information

Flexure of Thick Simply Supported Beam Using Trigonometric Shear Deformation Theory

Flexure of Thick Simply Supported Beam Using Trigonometric Shear Deformation Theory International Jornal of Scientific and Research Pblications, Volme, Isse 11, November 1 1 ISSN 5-15 Flere of Thick Simply Spported Beam Using Trigonometric Shear Deformation Theory Ajay G. Dahake *, Dr.

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

Development of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model

Development of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti

More information

Bounded UDE-based Controller for Systems with Input Constraints

Bounded UDE-based Controller for Systems with Input Constraints Bonded UDE-ased Controller for Systems with Inpt Constraints Yeqin Wang, Beiei Ren, and Qing-Chang Zhong Astract The ncertainty and distrance estimator (UDE- ased control is a rost control method, which

More information

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007 Graphs and Networks Lectre 5 PageRank Lectrer: Daniel A. Spielman September 20, 2007 5.1 Intro to PageRank PageRank, the algorithm reportedly sed by Google, assigns a nmerical rank to eery web page. More

More information

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model . ozdog and W. W. Olson Complex Tire-Grond Interaction Simlation: ecent eelopments Of n danced Shell Theory ased Tire odel EFEECE: ozdog. and Olson W. W. Complex Tire-Grond Interaction Simlation: ecent

More information

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary

More information

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS COLLEEN ACKERMANN AND WILL HARDESTY Abstract. The Cahn-Morral System has often been sed to model the dynamics of phase separation in mlti-component

More information

Numerical Model for Studying Cloud Formation Processes in the Tropics

Numerical Model for Studying Cloud Formation Processes in the Tropics Astralian Jornal of Basic and Applied Sciences, 5(2): 189-193, 211 ISSN 1991-8178 Nmerical Model for Stdying Clod Formation Processes in the Tropics Chantawan Noisri, Dsadee Skawat Department of Mathematics

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

Study on the Mathematic Model of Product Modular System Orienting the Modular Design

Study on the Mathematic Model of Product Modular System Orienting the Modular Design Natre and Science, 2(, 2004, Zhong, et al, Stdy on the Mathematic Model Stdy on the Mathematic Model of Prodct Modlar Orienting the Modlar Design Shisheng Zhong 1, Jiang Li 1, Jin Li 2, Lin Lin 1 (1. College

More information

Study on Combustion Characteristics of the Blast Furnace Gas in the Constant Volume Combustion Bomb

Study on Combustion Characteristics of the Blast Furnace Gas in the Constant Volume Combustion Bomb Proceedings of the 7th WEA International Conference on Power ystems, Beijing, China, eptemer 15-17, 2007 115 tdy on Comstion Characteristics of the Blast Frnace Gas in the Constant Volme Comstion Bom IU

More information

System identification of buildings equipped with closed-loop control devices

System identification of buildings equipped with closed-loop control devices System identification of bildings eqipped with closed-loop control devices Akira Mita a, Masako Kamibayashi b a Keio University, 3-14-1 Hiyoshi, Kohok-k, Yokohama 223-8522, Japan b East Japan Railway Company

More information

Exact Free Vibration of Webs Moving Axially at High Speed

Exact Free Vibration of Webs Moving Axially at High Speed Eact Free Viration of Wes Moving Aially at High Speed S. HATAMI *, M. AZHARI, MM. SAADATPOUR, P. MEMARZADEH *Department of Engineering, Yasouj University, Yasouj Department of Civil Engineering, Isfahan

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

Chaotic and Hyperchaotic Complex Jerk Equations

Chaotic and Hyperchaotic Complex Jerk Equations International Jornal of Modern Nonlinear Theory and Application 0 6-3 http://dxdoiorg/0436/ijmnta000 Pblished Online March 0 (http://wwwscirporg/jornal/ijmnta) Chaotic and Hyperchaotic Complex Jerk Eqations

More information

Professor Terje Haukaas University of British Columbia, Vancouver The M4 Element. Figure 1: Bilinear Mindlin element.

Professor Terje Haukaas University of British Columbia, Vancouver   The M4 Element. Figure 1: Bilinear Mindlin element. Professor Terje Hakaas University of British Colmbia, ancover www.inrisk.bc.ca The M Element variety of plate elements exist, some being characterized as Kirchhoff elements, i.e., for thin plates, and

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 5 (22) 265 292 NONLINEARITY PII: S95-775(2)349-4 A scalar nonlocal bifrcation of solitary waes for copled nonlinear Schrödinger systems Alan R Champneys and

More information

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L. .4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

6.4 VECTORS AND DOT PRODUCTS

6.4 VECTORS AND DOT PRODUCTS 458 Chapter 6 Additional Topics in Trigonometry 6.4 VECTORS AND DOT PRODUCTS What yo shold learn ind the dot prodct of two ectors and se the properties of the dot prodct. ind the angle between two ectors

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) MAE 5 - inite Element Analysis Several slides from this set are adapted from B.S. Altan, Michigan Technological University EA Procedre for

More information

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks Nmerical verification of the existence of localization of the elastic energy for closely spaced rigid disks S. I. Rakin Siberian State University of transport Rssia, 6349, Novosibirsk, Dsy Kovalchk street,

More information

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE 13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth

More information

Differential Geometry of Surfaces

Differential Geometry of Surfaces Differential Geometry of urfaces Jordan mith and Carlo équin C Diision, UC Berkeley Introduction These are notes on differential geometry of surfaces ased on reading Greiner et al. n. d.. Differential

More information

Similarity Solution for MHD Flow of Non-Newtonian Fluids

Similarity Solution for MHD Flow of Non-Newtonian Fluids P P P P IJISET - International Jornal of Innovative Science, Engineering & Technology, Vol. Isse 6, Jne 06 ISSN (Online) 48 7968 Impact Factor (05) - 4. Similarity Soltion for MHD Flow of Non-Newtonian

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Interrogative Simulation and Uncertainty Quantification of Multi-Disciplinary Systems

Interrogative Simulation and Uncertainty Quantification of Multi-Disciplinary Systems Interrogative Simlation and Uncertainty Qantification of Mlti-Disciplinary Systems Ali H. Nayfeh and Mhammad R. Hajj Department of Engineering Science and Mechanics Virginia Polytechnic Institte and State

More information

08.06 Shooting Method for Ordinary Differential Equations

08.06 Shooting Method for Ordinary Differential Equations 8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method

More information

Magnitude-Recurrence Relationship and its Effect on Uniform Hazard Spectra: A Current Assessment

Magnitude-Recurrence Relationship and its Effect on Uniform Hazard Spectra: A Current Assessment Magnitde-Recrrence Relationship and its Effect on Uniform Hazard Spectra: A Crrent Assessment T. K. Sen PhD (London), CEng., FICE Granherne, UK, a ssidiary of Kellogg Brown & Root Ltd SUMMARY: Gtenerg

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion rocedre (demonstrated with a -D bar element problem) MAE - inite Element Analysis Many slides from this set are originally from B.S. Altan, Michigan Technological U. EA rocedre for Static Analysis.

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information

Mixed Type Second-order Duality for a Nondifferentiable Continuous Programming Problem

Mixed Type Second-order Duality for a Nondifferentiable Continuous Programming Problem heoretical Mathematics & Applications, vol.3, no.1, 13, 13-144 SSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Mied ype Second-order Dality for a Nondifferentiable Continos Programming Problem.

More information

MEC-E8001 Finite Element Analysis, Exam (example) 2017

MEC-E8001 Finite Element Analysis, Exam (example) 2017 MEC-E800 Finite Element Analysis Eam (eample) 07. Find the transverse displacement w() of the strctre consisting of one beam element and po forces and. he rotations of the endpos are assmed to be eqal

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

Modelling, Simulation and Control of Quadruple Tank Process

Modelling, Simulation and Control of Quadruple Tank Process Modelling, Simlation and Control of Qadrple Tan Process Seran Özan, Tolgay Kara and Mehmet rıcı,, Electrical and electronics Engineering Department, Gaziantep Uniersity, Gaziantep, Trey bstract Simple

More information

NEW PERIODIC SOLITARY-WAVE SOLUTIONS TO THE (3+1)- DIMENSIONAL KADOMTSEV-PETVIASHVILI EQUATION

NEW PERIODIC SOLITARY-WAVE SOLUTIONS TO THE (3+1)- DIMENSIONAL KADOMTSEV-PETVIASHVILI EQUATION Mathematical and Comptational Applications Vol 5 No 5 pp 877-88 00 Association for Scientific Research NEW ERIODIC SOLITARY-WAVE SOLUTIONS TO THE (+- DIMENSIONAL KADOMTSEV-ETVIASHVILI EQUATION Zitian Li

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the

More information

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas

More information

Gravitational Instability of a Nonrotating Galaxy *

Gravitational Instability of a Nonrotating Galaxy * SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy

More information

Stability And Unbalance Response Of Rotor Bearing System

Stability And Unbalance Response Of Rotor Bearing System Stability And Unbalance Response Of Rotor Bearing Syste T.V.V..N. Rao and. Athre Mechanical Engineering Grop Departent of Mechanical Engineering Birla Institte of Technology & Science Indian Institte of

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

Lewis number and curvature effects on sound generation by premixed flame annihilation

Lewis number and curvature effects on sound generation by premixed flame annihilation Center for Trblence Research Proceedings of the Smmer Program 2 28 Lewis nmber and crvatre effects on sond generation by premixed flame annihilation By M. Talei, M. J. Brear AND E. R. Hawkes A nmerical

More information

Low-emittance tuning of storage rings using normal mode beam position monitor calibration

Low-emittance tuning of storage rings using normal mode beam position monitor calibration PHYSIAL REVIEW SPEIAL TOPIS - AELERATORS AND BEAMS 4, 784 () Low-emittance tning of storage rings sing normal mode beam position monitor calibration A. Wolski* Uniersity of Lierpool, Lierpool, United Kingdom

More information

SOIL NON-LINEAR BEHAVIOR AND HYSTERETIC DAMPING IN THE SPRING-DASHPOT ANALOG

SOIL NON-LINEAR BEHAVIOR AND HYSTERETIC DAMPING IN THE SPRING-DASHPOT ANALOG SOIL NON-LINEAR BEHAVIOR AND HYSTERETIC DAMPING IN THE SPRING-DASHPOT ANALOG Nikolaos OROLOGOPOULOS 1 and Dimitrios LOUKIDIS 2 ABSTRACT This paper presents reslts from nmerical simlations of footing vibration

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

JEPPIAAR ENGINEERING COLLEGE

JEPPIAAR ENGINEERING COLLEGE JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III

More information

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory Entrop 3, 5, -3; doi:.339/e54 Article OPEN ACCESS entrop ISSN 99-43 www.mdpi.com/jornal/entrop Deriation of D Power-Law Velocit Distribtion Using Entrop Theor Vija P. Singh,, *, stao Marini 3 and Nicola

More information

ON THE LIMITS OF THE BETZ S EFFICIENCY COEFFICIENT OF WIND TURBINES WITH THE HORIZONTALLY SHAFT

ON THE LIMITS OF THE BETZ S EFFICIENCY COEFFICIENT OF WIND TURBINES WITH THE HORIZONTALLY SHAFT ON THE LIMITS OF THE BETZ S EFFICIENCY COEFFICIENT OF WIND TURBINES WITH THE HORIZONTALLY SHAFT Prof.dr. ing. Petre Terzi SOCIETATEA PENTRU PROMOVAREA ENERGIILOR REGENERABILE, INEPUIZABILE SI NOI SPERIN

More information

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL 8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a

More information

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS

AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS th Eropean Conference on Comptational Mechanics (ECCM ) 7th Eropean Conference on Comptational Flid Dynamics (ECFD 7) 5 Jne 28, Glasgow, UK AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD

More information

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY APPLIED PHYSICS MEDICAL WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY L. CÃPITANU, A. IAROVICI, J. ONIªORU Institte of Solid Mechanics, Romanian Academy, Constantin Mille 5, Bcharest Received

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Lecture 7 Waveguides. TC 412 Microwave Communications

Lecture 7 Waveguides. TC 412 Microwave Communications Lectre 7 Wavegides TC 41 Microwave Commnications RS 1 Review Impedance matching to minimie power relection rom load Lmped-element tners Single-stb tners Microstrip lines The most poplar transmission line

More information

Nonlinear Dynamics of Thick Composite Laminated Plates Including the Effect of Transverse Shear and Rotary Inertia

Nonlinear Dynamics of Thick Composite Laminated Plates Including the Effect of Transverse Shear and Rotary Inertia The st Regional Conference of Eng. Sci. NUCEJ Spatial ISSUE ol.no. 8 pp 5- Nonlinear Dnamics of Thick Composite Laminated Plates Inclding the Effect of Transerse Shear and Rotar Inertia Mhsin J. Jeeg Sameera

More information

Diffraction of light due to ultrasonic wave propagation in liquids

Diffraction of light due to ultrasonic wave propagation in liquids Diffraction of light de to ltrasonic wave propagation in liqids Introdction: Acostic waves in liqids case density changes with spacing determined by the freqency and the speed of the sond wave. For ltrasonic

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

Motion in an Undulator

Motion in an Undulator WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Sven Reiche :: SwissFEL Beam Dnamics Grop :: Pal Scherrer Institte Motion in an Undlator CERN Accelerator School FELs and ERLs On-ais Field of Planar Undlator For planar

More information

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 2015 The athors and IOS Press. This article is pblished online with Open Access by IOS Press and distribted nder the terms of the Creative

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

1 Differential Equations for Solid Mechanics

1 Differential Equations for Solid Mechanics 1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to

More information

LPV control of an active vibration isolation system

LPV control of an active vibration isolation system 29 American Control Conference Hyatt Regency Rierfront, St. Lois, MO, USA Jne 1-12, 29 ThC15.1 LPV control of an actie ibration isolation system W.H.T.M. Aangenent, C.H.A. Criens, M.J.G. an de Molengraft,

More information

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica Visalisations of Gssian and Mean Cratres by Using Mathematica and webmathematica Vladimir Benić, B. sc., (benic@grad.hr), Sonja Gorjanc, Ph. D., (sgorjanc@grad.hr) Faclty of Ciil Engineering, Kačićea 6,

More information

University of Bristol

University of Bristol Uniersity of Bristol AEROGUST Workshop 27 th - 28 th April 2017, Uniersity of Lierpool Presented by Robbie Cook and Chris Wales Oeriew Theory Nonlinear strctral soler copled with nsteady aerodynamics Gst

More information