# Ch3 - CHAPTER 3 LONGITUDINAL WAVES IN RODS 3.1 EQUATION OF...

This preview shows pages 1–13. Sign up to view the full content.

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: CHAPTER 3 LONGITUDINAL WAVES IN RODS 3.1 EQUATION OF MOTION We consider the longitudinaﬂ wave motion in an elastic body whose E length is much greater than its cross—sectional dimensions, i.e., a rod. We neglect any {arenaﬂ coninaciion 0d expanéion of the rod, and hence this theory is an approximate one based on a Athengih 06 matehiaﬁA theony. The approximation is very good if the wavelengths of the motion are long compared with the rod's cross—sectional dimensions. In addition, it is also assumed that plane cross—sections of the rod remain plane and parallel during the deformation associated with the wave motion. Consider a section of an elastic rod, A+dA Let 0 be the normal stress in the X-direction and u be the longitudinal displacement of any cross—section in the X~direction. Both 0 and u depend I only on the variables x and t. Summing forces in the x—direction gives i 2 (O+d0)(A+dA) — (G)A = OK dx §—% , (3.1) at where p is the mass density of the material, and K'is the average cross— sectional area. Equation (3.1) reduces to 2 odA+dOA= pA dXB—l‘zl , at or odA + doA = 6* azu dX 8t2 Taking the limit as dx + 0 yields 2 3 _ 3 u SE-(Ao) — pA 3;? . (3-2) For the case A = constant, we would have swat. , (3.3) 2 u 31:2 Now from strength of materials theory with small deformations O = E8 (3.4) 8:22 3X ’ where 8 is the normal strain in the x—direction and E is modulus of elasticity. Combining (3.3) and (3.4) gives 32u azu 9 3X St or (3.5) =__l_u uxx 2 tt CL Where CL = XE , is the speed of propagation of the longitudinal disturbance. * We then see that the dtépﬁaccmcnt u(x,t) is governed by the one—dimensional wave equation. * It can also be shown that the Aineéé o(x,t), must also satisfy the one— dimensional wave equation (see homework exercise). ’1 1 l i l < E l l 3—3 Common boundary conditions which are used in conjunction with (3.5) are at a ﬁﬂee and, G = O e> EB = 0, while at a ﬁixed and, u = 0. 3x 3-2 TRAVELING WAVE SOLUTION RESULTS Recall the traveling wave solution to (3.5) was u(X,t) = f(x—ct) + g(x+ct) , where c = c . L Consider the case of an impact of the rod with velocity v to produce a wave traveling in the +x—direction. We therefore would have i u(X,t) = f(x~ct) . i l I Note that g E Bu , E =—-—-'=— i V at cf i and 39-9-. 8X _ E _ f ' Consequently, ,9 = _ X E c or E O = —pvc . (3.6) I Note for a wave traveling in the —x—direction we would find I i O = pvc . Equation (3.6) then gives the magnitude of the impact stress in terms of i the impact velocity. One can show the various reflection phenomena using the traveling wave solution. The results are that i g i ‘ 3—4 a u \$+F€\$S i.) at a ﬁnea end, 0 = 0, -—-= O; a reflected wave has the same 3x A magnitude as the incident wave except opposite in sign, and at the end the displacement is doubled, sirgse ‘ ii.) at a ﬁéxed and, u = O; a reflected wave has the same magnitude /\ and Sign as the incident wave, and at the end the stress is doubled. 3.3 §§3ARAT10N OF VARIABLES SOLUTION Consider the free vibrations of an unrestrained rod of length R with a prescribed initial displacement rm“ ‘; —-~ m % The problem is formulated as O<X<5Z, u =~l——u XX 2 tt ’ —- —- c O(O,t) = o(£,t) = O ut(x,0) = 0 I i i l E E (3.7) 2 u(x,0) = f(X) To separate the variables use the solution form l u(x,t) = X(X)T(t) g i ! This produces the forms X = A sinkx + B coskx l H H C sinkct + D coslct . Now from conditions (3.7)2 3u(0,t) = sumo = 3x 3x 0 ’ O(0,t) = o(5L,t) = 0 => 3—5 so A = 0 sinAR = O é> A = ——'; n = l, 2, Note that the n = 0 case corresponds to the A=O case, which gives XT = constant. Relation (3.7)3 implies that Using superposition our solution may then be written as 8 a o nﬁx nﬁct u(x,t) = Er‘+ 1 an cos ~Ef-cos g . (3.8) “M 11 Finally applying the last condition (3.7)4 yields that a 00 n'ITX f(x) = 352 + i an cos -7Z- , (3.9) n=l which is the Fourier cosine series for f(x); hence 2 an = %-f0f(£) cos 3%; dE . (3.10) Equations (3.8) and (3.10) then solve the problem for an arbitrary function f. 3'4 MW Consider now the longitudinal wave motion along a composite rod composed of two semi—infinite rods of different material and cross—section. Assume a perfectly rigid bonding of the two rods. 3—6 Working in terms of the stress wave motion we write oi = incident stress pulse or = reflected stress pulse 0t = transmitted stress pulse Now at the discontinuity of cross-section we must have the following matching conditions holding, (oi + or)Al = o A t 2 (3.11) vi + vr = vt where vi r are the incident, reflected or transmitted velocities, respectively. ! ’ Using (3.6) in (3.11)2 yields _ 0i + gr _ Ot _' " a plcl plcl pzcz or O =.l (g _ 0 ) (3.12) t u i r ’ p C . - . Where u = l l . prodigch P4; 13 m/fccl wafer/cl Impedance pzcz Solving (3.11)l and (3.12) for or and o in terms of oi produces t A 1—14-—; A2 0 = U. r Al 1 l+u K— 2 (3.13) A1 2 _ A2 0t Al Oi . l+Ll K“ 2 Note it follows from (3.13) that at a fixed end 3—7 while at a free end A + o :> 0r = —o. 3 o§;:c3 2 1 g 3-5 MW We now want to consider the case of variable cross—section where A = A(X). The rod or wave guide will now affect the wave amplitude because of the geometrical variation. x i l 5 Recall the governing equation (3.2) read 2 gym) =pA8—% 8t Using (3.4) in the above yields 2 3 Bu 3 u —~ [EA ——] = pA ——— 3X 3X atZ Assuming E = constant, we get 11mg =_1_82u 5 A 8X 3x C2 at2 or (3.14) Bzu + g;_gy£ = g;_ 82u 3x2 A 3X c2 Btz Consider next the case of a harmonic time dependent solution form, i.e., let u(x,t) = U(x)ei‘*’t . (3.15) Using (3.15) in (3.14) gives 2 A' w n _ v ____ = U + A U + 2 U 0 . (3.16) c For the case where A = constant (i.e., A' = 0), we get (1) cos —-x c U: . w . (3.17) Sln'— x c Now assume A' is small, so that we are investigating slightly varying rods. In addition try the solution form U(X) = f(X) sin g(X) . (3.18) Substituting form (3.18) in (3.16) produces 2 A' w2 [f" — f(g') + Kf~f' +-*§ f] sin g(X) C (3.19) I + [Zf'g' + fg” + %— fg'] cos g(X) = O For (3.19) to be zero for all x we must have 2 2 A' w n_ l _____ v __ = f f(g) +A f +c2f o A! 2f'g' + fg" + K—~fg' = O Rewriting equation (3.20)2 gives 1 u I 2.1+g_f.+é._=0 f g A Integrating the above yields 2 1n f + 1n g' + 1n A = Cl , where C1 = arbitrary constant. or 1n[f2Ag'] = Cl C L szg' = e l = C2 , say . (3.21) I Now in (3.20)l, assuming A' is small the term %— f' is negligible in comparison to the other terms. In addition, we take f" << 1, so the term f" is also negligible. Therefore (3.20)l becomes approximately 2 —f(g')2+%f = 0 . C Using (3.21) in the above yields Solving for f gives cC f=V ——2=KA_l/2 , (3.22) (DA /CzcI where K = —~—-. U) Hence, if we have a gradual variation in the cross—sectional area, then we will have a small variation in the amplitude f. Using (3.21) and (3.22) gives §=C2=9 sz C Integrating the above g=gx+c3, 02% where C “ constant. 3 l Combining results (3.22) and (3.23) in (3.18) gives U(x) ='E— sin(%'x + C3) . (3.24) JR Note the comparison of (3.24) with the A' = 0 case (3.17). 3—10 3-6 VIBRATION OF RODS WITH RIGID END MASSES Consider the case of the free vibration of a finite uniform rod with a rigid end mass. The problem is formulated as _ L 1 uxx — 2 utt , O f X_: £ C 1 u(0,t) = 0 (3.25) : AB R — W i uX( \$t) “ _ g utt( at) Separating the variables yields u(x,t) = X(X)T(t) X(x) = Cl cosXx + C2 sinAx (3.26) T(t) II C cosACt + C4 Sinkct 3 Condition (3.25)2 gives while (3.25)3 implies that ' -- W.— 2 2 ‘ o AE(C3 cosAct + C4 Slnkct)CZX coslﬂ — g A c C2 Sinl£ ' (C3 coskct + C4 sinAct) The previous relation may be written as AEg = tanAI . (3.27) Wk c2 3—ll Letting CL = AEgSL Q m mrmm‘s of rod Wcz MM ) M: was: of rigid wag/174 B = Xi , relation (3.27) then becomes a = B tanB . (3.28) Relation (3.28) is a thanccndentaﬂ cquaiion for the eigenvalues X. It is sometimes called the ﬁhcqucncy equaiion of the system. The solution to (3.28) can be graphically illustrated by the following figure. g i 3 E i K . I E Hence from the figure for a given a, there will be k—roots (k = l, 2, ...) to equation (3.28). Call the roots Bk. Hence we write ka Bkct Bkct uk(x,t) = sin —E—~(C3k cos % + 04k sin z ) , (3.29) which represents the motion for the k—th phincépaﬁ mode, The complete solution is the superposition of all modes, i.e., 3—12 ka Bkct . Bkct 3k cos + C Sln u(X,t) = f sin ——— (C k=l £ z 4k 2 ) ’ (3'30) where 03k and C4k would be determined from the initial conditions. 3-7 L_A_TERAL INERTIA EFFECTS (Lame Thea/3]) We mentioned in the derivation of the governing equation (3.5), that the effects of lateral contractionior expansion would be neglected. We now re—examine this point and reformulate the analysis to include this effect. This modified analysis still follows a strength of materials approach II and initially assumes that 0y = O 0. So from Hooke's law, the transverse Z strains are given by _ .2 _ _ _ _ .Ee €y—€Z—-EOX- vex—vax, 7 (3.31) where v is POM/son’xs We. The lateral displacements in the y and 2 directions are taken to be distributed linearly and so we write v=€y=-\)y-82 y 3X (3.32) w=€z=—\)ziq-. 2 8X With all the variables now expressed solely in terms of u(X,t), the governing equation of motion can be derived by using the eneagy method. 3—13 This method incorporates Hamiﬂton’é pntncépﬁc, and states that a system moves in such a way as to make the total kinetic energy and system work an * exiaemum. The details of this process are involved , so we give here only the end result 82u vZEZ 34u _ l 32u ——+ ——-—————-————- , (3.33) 3x2 c 2 3x28t2 2 3t2 L CL where k is the poﬁan aadtué 06 gynaiion of the cross—section, i.e., Ip = sz with IP being the poﬁad moment 06 anatia of the cross—sectional area. This particular longitudinal wave theory is sometimes called Love's theory or Rayleigh's theory. Equation (3.33) is not the one—dimensional wave equation, and so we can expect significantly different propagational characteristics for this case. Attempting a harmonic traveling wave solution form produces the fact that, when lateral inertia is included, we have a ﬁaequencg dependent phaée ucﬂocity, i.e. dibpehAion (see homework exercise). 7': For details see Graff or Love. E g E 3 § s ...
View Full Document

## This note was uploaded on 10/03/2011 for the course MCE 565 taught by Professor Staff during the Spring '11 term at Rhode Island.

### Page1 / 13

Ch3 - CHAPTER 3 LONGITUDINAL WAVES IN RODS 3.1 EQUATION OF...

This preview shows document pages 1 - 13. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online