SOLID M! Al DYNAMICS TUTORIAL -NATURAL VIBRATIONS — ONE DEGREE OF FREEDOM This work covers elements of the syllabus for the Engineering Council Exam D225 - Dynamics of Mechanical Systems, C105 Mechanical and Structural Engineering and the Edexcel HNC/D module Mechanical Science. Oscillations with two degrees of freedom are covered in other tutorials and these are vital for those sitting exam D225. On completion of this tutorial you should be able to do the following. Explain the meaning of “degrees of freedom”. Define and explain Simple Harmonic Motion Explain the meaning of free vibrations. Show how various systems vibrate with simple harmonic motion. Determine the natural frequency and periodic time for simple systems. Determine the displacement, velocity and acceleration of bodies vibrating with simple harmonic motion. Explain the use of energy to analyse vibrations (Rayleigh’s Method) Explain how to combine several natural frequencies (Dunkerley’s Method) This tutorial covers the theory of natural vibrations and brings together other areas of studies in the process. To do the tutorial fully you must be familiar with the following concepts, Spring Force. Inertia. Moment of inertia. Moment of force. Law’s of motion. Basic calculus. The deflection of beams. The twisting of shafts. Strain energy. Buoyancy. All these topics may be found in other tutorials. © D.LDUNNContents 1. Free Vibrations 1.1 Degrees of Freedom 1.2 Free Vibrations - Examples 1.3 Simple Harmonic Motion. 1.4 Angular Frequency, Frequency and Periodic Time 1.5 Equations of Simple harmonic Motion. . Analysis of Natural Vibrations. Simple Pendulum. Linear Elastic Vibrations. 4.1 Mass-Spring System 4.2 Transverse Vibrations (of beams) 4.3 Energy Methods (Rayleigh) 4.4 Transverse Vibrations duc to the distributed mass. 4.5 Combination of Distributed and Point Loads (Dunkerley) Torsional Oscillations. Oscillation of a Floating Body © D.LDUNN 2REE VIBRATION: 1.1 DEGREES OF FREEDOM The numbers of degrees of freedom that ° ‘a body possesses are those necessary to 7 completely define its position and = orientation in space. This is useful in several fields of study such as robotics One degree and vibrations. Two degrees Figure | Consider a spherical object that can only be positioned somewhere on the x axis. This needs only one dimension, ‘x’ to define the position to the centre of gravity so it has one degree of freedom. If the object was a cylinder, we also need an angle ‘@” to define the orientation so it has two degrees of freedom. Now consider a sphere that can be positioned in Cartesian coordinates anywhere on the z plane. This needs two coordinates ‘x’ and “y’ to define the position of the centre of gravity so it has two degrees of freedom. A cylinder, however, needs the angle “0” also to define its orientation in that plane so it has three degrees of freedom, 4 Two degrees Three degrees Figure 2 “ In order to completely specify the position and orientation of a cylinder in Cartesian space, we would need three coordinates x, y and z and three angles relative to cach angle. This makes six degrees of freedom. In the study of free vibrations, we will be constrained to one degree of freedom, Six degrees Figure 3 © D.LDUNN 31.2 FREE VIBRATIONS- EXAMPL! A free vibration is one that occurs naturally with no energy being added to the vibrating system. The vibration is started by some input of energy but the vibrations die away with time as the energy is dissipated. In each case, when the body is moved away from the rest position, there is a natural force that tries to return it to its rest position. Here are some examples of vibrations with one degree of freedom. A pendulum swings Dak: aud forth Amass goes up down on a spuing Figure 4 Note that the mass on the spring could be made to swing like a pendulum as well as bouncing up and down and this would be a vibration with two degrees of freedom. The motion that all these examples perform is called SIMPLE HARMONIC MOTION (S.H.M.). This motion is characterised by the fact that when the displacement is plotted against time, the resulting graph is basically sinusoidal. Displacement can be linear (e.g. the distance moved by the mass on the spring) or angular (e.g. the angle moved by the simple pendulum). Although we are studying natural vibrations, it will help us understand S.H.M. if we study a forced vibration produced by a mechanism such as the Scotch Yoke. 1.3 SIMPLE HARMONIC MOTION = eR ina) 7 sin = Rsiniot) . ig eee . acceleration a Figure 5 Figure 6 The wheel revolves at @ radians/sec and the pin forces the yoke to move up and down. The pin slides in the slot and Point P on the yoke oscillates up and down as it is constrained to move only in the vertical direction by the hole through which it slides. The motion of point P is simple harmonic motion. Point P moves up and down so at any moment it has a displacement x, velocity v and an acceleration a © D.LDUNN 4The pin is located at radius R from the centre of the wheel. The vertical displacement of the pin from the horizontal centre line at any time is x. This is also the displacement of point P. The yoke reaches a maximum displacement equal to R when the pin is at the top and-R when the pin is at the bottom. This is the amplitude of the oscillation, If the wheel rotates at « radian/sec then after time t seconds the angle rotated is @ = ct radians. From the right angle triangle we find x = R sin(ot) and the graph of x - 6 is shown on figure 7a. Velocity is the rate of change of distance with time and in calculus form v = dx/dt, If we differentiate x we get v = dx/dt = OR cos(at). The plot is also shown on figure 7a, The maximum velocity or amplitude is @R and this occurs as the pin passes through the horizontal position and is plus on the way up and minus on the way down. This makes sense since the tangential velocity of a point moving in a circle is v= @R and at the horizontal point they are the same thing, Acceleration is the rate of change of velocity with time and in calculus form a = dvi Differentiating v we get a = dv/dt = ~@?R sin(ot). The plot is also shown on figure 7a. The amplitude is R and this is positive at the bottom and minus at the top (when the yoke is about to change direction) ox Since R sin(«t) =x then substituting x we find This is the usual definition of S.H.M. The equation tells us that any body that performs sinusoidal motion must have an acceleration that is directly proportional to the displacement and is always directed to the point of zero displacement. The constant of proportionality is @”, Any vibrating body that has a motion that can be described in this way must vibrate with S.H.M. and have the same equations for displacement, velocity and acceleration. 1.4 ANGULAR FREQUENCY, FREQUENCY AND PERIODIC TIME is the angular velocity of the wheel but in any vibration such as the mass on the spring, it is called the angular frequency as no physical wheel exists The frequency of the wheel in revolutions/second is equivalent to the frequency of the vibration. If the wheel rotates at 2 rev/s the time of one revolution is 1/2 seconds. If the wheel rotates at 5 rev/s the time of one revolution is 1/5 second, If it rotates at f rev/s the time of one revolution is I/F. This formula is important and gives the periodic time. Periodic Time T = time needed to perform one cycle. fis the frequency or number of cycles per second, It follows that Wf and f= 1/T Each cycle of an oscillation is equivalent to one rotation of the wheel and 1 revolution is an angle of 2n radians. When 0=2n, t=T. Tt follows that since @ = wt then 2x = oT Rearrange and @ = 2n/T. Substitute T= 1/f and @ =2nf © D.LDUNN 51.5 EQUATIONS OF S.H.M, Consider the three equations derived earlier. Displacement x=Rsin(wt). Velocity v= dx/dt = oR cos(at) Acceleration a= dv/dt =-o2R sin(ot) The plots of x, v and a against angle 0 are shown on figure 7a, In the analysis so far made, we measured angle 0 from the horizontal position and arbitrarily decided that the time was zero at this, point. Suppose we start the timing after the angle has reached a value of $ from this point. In these cases, $ is called the phase angle, The resulting equations for displacement, velocity and acceleration are then as follows. Displacement xR sin(ot + 9). Velocity v= dwidt = oR cos(at + 6) Acceleration a=dv/dt=-o2R sin(ot + $) The plots of x, v and a are the same but the vertical axis is displaced by ¢ as shown on figure 7b. A point to note on figure 7a and 7b is that the velocity graph is shifted ‘4 cycle (90°) to the left and the acceleration graph is shifted a further “ cycle making it ‘4 cycle out of phase with x. YH. ’ velocity a) \ ee ’ ( Figure 7a Figure 7b © D.LDUNN 6WORKED EXAMPLE No.1 The displacement of a body performing simple harmonic motion is described by the following equation x=A sin (ot + 4) where A is the amplitude, @ is the natural frequency and @ is the phase angle. Given A = 20 mm, @ = 50 rad/s and $ = 7/8 radian, calculate the following. i. The frequency. ii, The periodic time. iii. The displacement, velocity and acceleration when t = 1/4, Sketch the raphs of x, v and a and confirm your answers. SOLUTION First deduce the frequency. f = @/2m = S0/2n = 7.96 Hz. Next deduce the periodic time. T = 1/f = 0.126 s Next deduce the time t. t= T/4 = 0.0314 s Next write out the equation for displacement and solve x at t= 0.0314 s x=20 sin{ 60 x 0.314) +2) \ 8) rg) x= 20sin| 1.57+ ry = 20sin1.963 = 18.48 mm \ Next write down the equations for v and a x= 20sin(ot+ 9) Vv =20e cos(t +) =20x 50x cos(1.963) 382.2 mm/s a =-20’sin(ot + @) =-20 x 50* sin(1.963) =-46203 mm/s* The plots of x, v and a confirm these answers. we 2 = > =m ay |\ / « i \ E fim § Hime 3 ip P wT e wah wh Foon / i i ac] ™ —_ sooo Figure & © D.LDUNN 72 NATURAL VIBRATION: EXPLANATION In the following work we will show how some simple cases of natural vibrations are examples of simple harmonic motion, Remember that one important point common to all of them is that there must be a natural force that makes the body move to the rest position. Another point common to all the following examples is that the body has mass (inertia) and that in order to accelerate, there must be an inertia force or torque present. A free vibration has no external energy added after it starts moving so it follows that all the forces and all the moment of force acting on the body must add up to zero. This is the basis of the analysis. 3 SIMPLE PENDULUM The restoring force in this case is gravity. When the pendulum is displaced through an angle 0, the weight tries to restore it to the rest position, This analysis is based on moment of force (torque). RESTORING TORQUE Weight = mg To find any moment we must always use the distance to the centre of rotation measured at 90° (normal) to the direction of the force. In this case the distance is L sind Denote the torque as T, .This tries to return the mass to the rest position. T,= weight x mg (Lind) ‘Weight=mg Figure 9 INERTIA TORQU! Since the pendulum has angular acceleration as it slows down and speeds up, it requires an inertia torque to produce it. Denote this torque as T;. From Newton’s second law for angular motion T, = Ia «a. is the angular acceleration and I is the moment of inertia. The mass is assumed to be concentrated at radius L. (If it was not, the problem would be more complicated), The moment of inertia is then simply given as I= mL? BALANCE OF MOMENTS If there is no applied torque from any external source, then the total torque on the body must be zero, Te tT, mg LsinO + mL2a=0 mg LsinO =- mL2a. g Lsin0 =-L2a, The sin of small angles is very similar to the angle itself in radians (try it on your calculator). The smaller the angle, the truer this becomes. In such cases sin(0) = 0 radians and so we may simplify the equation to 20=-La a=-(g/L)0 This meets the requirements for SHM since the acceleration « is directly proportional to the displacement 0 and the minus sign indicates that it is always accelerating towards the rest point, It follows that the constant of proportionality is so (g/L). = g/L = (g/L)" © D.LDUNN 8If the displacement Q was plotted against time the same graphs as shown in figure 7 would result. Note that the displacement in this example is angle and should not be confused with the angle on the Scotch Yoke. The frequeney of oscillation is obtained from f= e/2n. 1 jg 2nVL Note that the mass makes no difference to the frequency. On earth we can only change the frequency by altering the length L. If we took the pendulum to the moon, it would oscillate more slowly since gravity is smaller. In outer space where g is very close to zero, the pendulum would have no weight and would not swing at all if moved to the side. Hene Remember also that the above equation is only true if the pendulum swings through a small angle If the angle is large, the motion is not perfect S.H.M. WORKED EXAMPLE No.2 ‘A mass is suspended from a string 60 mm long. It is nudged so that it makes a small swinging oscillation, Determine the frequency and periodic time. SOLUTION = (12m) (g/L) #5 (1/2n) (9.81/0.06) 4 = 2.035 Hz, = 1/f=0.491 s © D.LDUNN 94 LINEAR ELASTIC OSCILLATIONS Most natural oscillations occur because the restoring force is due to a spring. A spring is any elastic body which when stretched, compressed, bent or twisted, will produce a force or torque directly proportional to displacement, Examples range from the oscillation of a mass on the end of a spring to the motion of a tree swaying in the wind. Let's start with a simple mass suspended on a spring. 4.1 MASS ON SPRING Consider the mass is pulled down with a force F as shown. The spring is stretched a distance xq. This is called the initial displacement When the mass is the released, it oscillates up and down with simple harmonic motion. Let's analyse the forces involved, < Rest positon xe F is the applied force in Newtons. Mass Me ~ x is the displacement from the rest position at any time and k is the spring stiffness. Figure 10 F, = Spring force that tries to return the mass to the rest position. From spring theory we know that F, =k x Since the motion of the mass clearly has acceleration then there is an inertia force F. From Newton's second law of motion we know that Fi = mass x acceleration = M a Balancing forces gives F=F+F,=Matkx If the mass is disturbed and released so that it is oscillating, the applied force must be zero and this, is the requirement for it to be a free natural oscillation. 0-Ma+kx Rearrange to make a the subject a=-(h/m) x This is the equation for SHM and tells us the acceleration is directly proportional to displacement and it is directed towards the rest position. If x was plotted against time, a sinusoidal graph would k result, The constant of proportionality is 4 and this must be the angular frequency squared so k ® 1 1 [k The frequency of oscillation is f = 2. = — J ~ Qn 2nVM Because this is a natural oscillation the frequencies are often denoted as «, and f,. This equation is true for all elastic oscillations. © D.LDUNN 10WORKED EXAMPLE No.3 ‘A spring of stiffness 20 kN/m supports a mass of 4 kg. The mass is pulled down 8 mm and] released to produce linear oscillations. Calculate the frequency and periodic time. Sketch the} graphs of displacement, velocity and acceleration. Calculate the displacement, velocity and| acceleration 0.05 s after being released. The oscillation starts at the bottom of the cycle so xo =-8 mm. The resulting graph of x against time will be a negative cosine curve with an amplitude of 8 mm, The equations describing the motion are as follows. X=Xgcoset When t = 0.05 seconds x = -8 cos(70.71 x 0.05) x= 7.387 mm, (Note angles are in radian) This is confirmed by the graph. If we differentiate once we get the equation for velocity. V=-oxosin ot v= -oxgsin ct = -70.71 (-8)sin(70.71 x 0.05) v=-217 mm/s This is confirmed by the graph. Differentiate again to get the acceleration. a=-o'xcose! and since x = xgcos@t a= -0'x a=-70.71" x 7.387 = -36 934 mmis* This is confirmed by the graph, u i sooo vat 8 ‘ \ a {\ \ F 7 we fam § \ / : “ win al oom Figure 11 © D.LDUNN unSELF ASSESSMENT EXERCISE No.1 1, Calculate the frequency and periodic time for the oscillation produced by a mass — spring system given that the mass is 0.5 kg and the spring stiffiness is 3 N/mm. (12.3 Hz, 0.081 s). 2. A mass of 4 kg is suspended from a spring and oscillates up and down at 2 Hz. Determine the stiffness of the spring. (631.6 Nim), The amplitude of the oscillation is 5 mm. Determine the displacement, velocity and acceleration 0,02 s after the mass passes through the mean or rest position in an upwards direction, (1.243 mm, 60.86 mm/s and -196.4 mm/s’) 3. From recordings made of a simple harmonic motion, it is found that at a certain point in the motion the velocity is 0.3 mis and the displacement is 20 mm, both being positive downwards in direction. Determine the amplitude of the motion and the maximum velocity and acceleration. Write down the equations of motion. Note that the data given is at time t= 0, You will have to assume that X—x9cos(at + §) at time (-0 Ans. x= 0.0311 cos(wt - 50°) 3914 sin(at - 50°) a=-157.9x © D.LDUNN 1242 TRANSVERSE VIBRATIONS A transverse oscillation is one in which the motion is sideways to the length. This occurs in shafts and beams, They may be solved as either angular or linear motion and linear is chosen here. CANTILEVER \ Fore Force ‘Mass a SS AN / us Figure 12 Consider a mass on the end of a cantilever beam. If the mass is moved sideways to the beam, the beam bends and it will always be found that the force is directly proportional to displacement. In other words, it is simple transverse spring. The stiffness of the beam may be found from beam theory. For example in the tutorial on beams it was shown that the deflection of a cantilever due to the point load on the end is given by FL* 3EI y is the symbol for the deflection of a beam but since we have denoted the displacement as x in the studies so far so we will use x 3EI The stiffness is hence E_ 31 x LU ‘The theory is the same as a mass on the spring sof = |X -1 3B 2nYM 2 \ ML’ SIMPLY SUPPORTED BEAM 2 For a simply supported beam with a point load at the middle, the deflection is x= 7-7 and it follows that k= 48"! an Figure 13 © D.LDUNN 13sTA DEFLECTION ‘The stiffness may also be found by measuring the static deflection. Suppose the mass deflects a distance xg under its own weight. The force is the weight so F = Mg. It follows that This formula works for both a cantilever and a simply supported beam. Note that none of the forgoing work takes into account the mass of the beam or shaft and unless this is small in ‘comparison to the concentrated mass, the formulae does not give accurate answers. This is covered later WORKED EXAMPLE No.4 ‘A rod 20 mm diameter and 1.2 m long is rigidly fixed at one end and has a mass of 2 kg concentrated at the other end. Ignoring the weight of the rod, calculate the frequency of| transverse vibrations. Take E = 200 GPa. SOLUTION First calculate the second moment of area. For a circular section : aD! _nx0.02" sey 610 m SEL _ 1 [3x200x10°x7.854x10" _ 5 og MLE 2eY xt? WORKED EXAMPLE No.5 A horizontal shaft may be treated as a simply supported beam and has a mass of 40 kg placed at the middle and it deflects | mm under the weight. Ignore the mass of the shaft and calculate the frequency of transverse oscillations. LUTION rt fe 1 [st 2n\x, 2x V0.001 © D.LDUNN 14SELF ASSESSMENT EXERCISE 1. A shaft is mounted horizontally between two bearings. A mass is hung from the centre point and the shaft deflects 2 mm, Calculate the natural frequency of the shaft due to this weight only. (11.145 Hz) 2, A beam is 4 m long and rests on simple supports at each end. The beam has a second moment of area I of 800 x 10° m‘ and a modulus of elasticity of 180 GPa. It has a mass of 30 kg placed at the middle. Calculate the natural frequency of the shaft due to this weight only. (9.549 Hz) 43. ENERGY METHODS (RAYLEIGH) Rayleigh came up with a method of solving complex oscillations for mass — spring systems based on the fact during an oscillation the maximum kinetic energy of the oscillating mass is equal to the maximum strain (spring) energy. This may be illustrated at a basic level by considering a simple mass on a spring, Let the maximum deflection of the mass be x. The spring, force is F =k xo. The work done is F x= kx” Let the displacement at any time be > oping sttesk Nm X= X,8iN Of The velocity is V= 0X, cos ot The maximum value is @ x. The maximum kinetic energy = Mv” mar K.E.(max) = % Ma x Equate energy and Ys kx,?= Y% Moo? xo? = k/M and this is the same result found earlier. Figure 14 An important point to note is that the result would be the same whatever value is chosen for displacement. The deflection due to the static mass could be used and this would lead to the same result so it is arguable that the forgoing including the static deflection method, is a result of applying Rayleigh’s method, The method is particularly useful for transverse oscillations. © D.LDUNN 1sWORKED EXAMPLE No.6 Calculate the natural frequency of the system shown ignoring the mass of the shaft. The flexural stiffness is EI = 20 300 N m’. 580 mm 420 mam, 500 mm Figure 15 SOLUTION ‘We can use two different but related strategies. We could determine the strain energy or the static deflection. Let's use the strain energy method. First determine the reactions at A and C. 1m 05m f 20k Ra Re Figure 16 Moments about the left end yields 20 x 1.5 = Rex 1 Re=30 kg and hence Ra = 10 kg down. Ra=-98.1N Ry = 294.3 N Up and the load on the end is 196.2 N down, STRAIN ENERGY FOR SECTION Ato wax = Sone | 98.1x)'dx 0.079 Joules 2EI 2x 20500 STRAIN ENERGY FOR SE! TION CoB This is simplified if we measure x from the free end and then M = -196,2 x 1 1 os 1962? _[x3]"* U =< J Mdx = ————. [(- 196.2x)Pdx =§ <= 2ET 2x 20300 5, 2x 20300 | 3 ‘The total strain energy is U = 0.0395 + 0.079 = 0.1185 J 0.0395 Joules The deflection produced at the free end is Yn ‘The strain energy =" F ym nd F = weight = 20g = 196.2 N Hence yq,=2 x 0.1185/196.2 = 0.0012 m or 1.2 mm. Ifa mass is oscillating up and down sinusoidally y ~ ym sin wt The velocity is © ym 0s (ot) © D.LDUNN 16The kinetic energy is ¥ Mv"/2 = %& Mf @ yu cos (aot)}? Maximum value of ¥% M/2){@ Ym} Ym is the deflection at the end = 0.0021 m M = 20 kg Equating energies we have Mio Ym}? = 0.1185 10 w* (0.00217) = 0.1185 = 0.1185/{10 (0.00217)} = 8229 @=90.7 rad/s f= 90.7/(2n) = 14.43 Hz Note we could have found the static deflection by beam theory and then used £, = (/2n)V(glym) = (1/2m)V(g/0.0012) = 14.4 He, Also note that if the disc on the shaft is large, the disc will rotate about its diameter as the shat deflects and so angular kinetic energy will also be involved MA: Beams and shafis of many possible configurations (cantilever, simply supported, encastré and so on) will vibrate in much the same manner as the string on a musical instrument, Each configuration will produce a natural frequency depending on the mass, how it is distributed and the transverse spring stiffness, These oscillations can be approached by consideration of the energy changes involved. When a beam is bent, longitudinal layers become stretched and compressed so that strain energy is stored in it. When the beam oscillates it has a velocity at any moment in time and hence has kinetic energy. The strain energy and kinetic energy are converted back and forth from one to the other. The energy trapped in the system cannot get out and cannot stay static so this interchange of energy occurs with the total energy being conserved (assuming no damping). STRAIN ENERGY Consider a simply supported beam with negligible mass. Suppose a weight hanger is suspended from it at distance x metres from the end. Further suppose that weights are added to the hanger in small increments causing the beam to deflect a distance y. x mean or t rest position ———— weights Figure 17 © D.LDUNN 7wergnt If we conducted this experiment and plotted the weights W * Newtons against deflection y we would get a straight line graphs as shown. The maximum weight is W and the maximum deflection is ym The work done is the area under the graph (a triangle). Work Done=%W Yu Yon Figure 18 w Since an equal amount of energy has been used to deflect the beam, this energy must be stored in the beam as mechanical energy due to the stretching and compression of the material along the length. This is strain energy and we normally denote it as U. The strain energy stored in the beam is hence U=% W Ym ‘Now consider a beam that does not have a negligible mass and is uniform along its length with a distributed weight of w N/m. Consider a small length of the beam 8x at distance x metres from the left end, rest position = Z Figure 19 The weight of the element is w 5x Newtons. This is equivalent to the weight W in figure 17 so the strain energy is a small part of the total 6U = 4 w 8x yy The strain energy is due to the weight of the small length only so if we want the strain energy for the whole length of the beam we must integrate along the length. Reverting to calculus 8x—rdx and SUsdU then w FJ vadx KINETIC ENERGY If at any moment in time the velocity of the element is v then the kinetic energy of the element is “4 m v2 where Sm is the small part of the total mass. The mass of the element is In order to understand the velocity of the clement we need to consider the S.H.M., which it performs. ym is the amplitude of the motion, The deflection is zero at horizontal point so we can describe the displacement (y) with the following equation. The deflection at any moment is y= ymsinot he velocity is v=dy/dt= © ym cost The maximum velocity is hence @ ym. © D.LDUNN 18The time plots of y and v are as shown. 7 + The kinetic energy of the small f i bs an element at any time is then y y > 5 1 wé \ / 8(K.E.) = ~~ (ay,,cos ot)? -0y,g— 28 Yin Figure 20 ‘The maximum value is when the cosine has a value of 1. The maximum K.E. of the element is 8(K.E)max = 1" oy? x 2g ‘The maximum K.E. for the whole beam is found by integrating with respect to x. Again revert to caleulus 8x—odx and 8(K.E.)>d(K.F.) then wor, K.E(max) = ee 2 . setae wo? 2 Equating U and K.E. we have FJ ynte ay dont 2 _ 8) ¥ndx vind In order to evaluate this expression, the function relating the deflection yq, and position x must be known, In the case of a simply supported beam we know this is as follows. andf = 2 From this the angular frequency of the vibration is «@ =F za SIMPLY SUPPORTED BEAM [atx —x‘—L'x] and the solution (which is hard to do) is Ke , a a aaa 3 Figure 21 ELXED ENDS (Encastre Beam) The lowest critical speed is f = 3.562, et and the higher critical speeds are given by 1) [eer ney) Epa where m= 2,3, « Note that the lowest speed almost corresponds to n= 1 ) Vw ] | st moae second mode Encastre Figure 22 © D.LDUNN 19CANTILEVER, The lowest critical speed is f =0.565,, wi 1y = {n ——| 22) where n=2,3 ...... Note that the lowest speed almost corresponds to n= 1 20" 2) Vw Figure 23 45 COMBINATION OF DISTRIBUTED AND POINT LOADS. When a beam or shaft has one or more concentrated masses as well as a uniformly distributed mass, the frequency of oscillation may be found using Dunkerley’s method. This is a method, which enables the frequency of an oscillation to be deduced when it is due to two or more masses. For example a cantilever oscillates because of its distributed mass and because of any point loads on it. We know how to work these out separately. Suppose a beam oscillates a frequency f] due to one load on its own, £2 due to another load on its own and f3 due to a third load on its own. When all three loads vibrate together, the resulting frequency fis found using the reciprocal rule below, 1 1 1 1 This formula is based on observation and works very well. WORKED EXAMPLE No.7 A cantilever beam is 2 m long and has a cross section 60 mm wide and 40 mm deep. It is made from a material with a density of 2800 kg/m3 and a modulus of elasticity of 78 GPa. It has a mass of 5 kg attached on the end. Calculate the following. i, The natural frequency due to its own weight only. ii, The natural frequency due to the point load only. iii, The combined frequency. LUTION ‘The natural frequency due to its own weight is f= 0.56, wit First calculate the weight of | metre length (w). Volume = sectional area x length = 0.06 x 0.04 x 1 = 0.0024 m* Next change this to mass. m = density x volume = 2800 x 0.0024 = 6.72 kg/m. ext change this to weight. w = mg = 6.72 x 9.81 = 65.92 Nim Next calculate the second moment of area about the horizontal centre line. BD’ _ 0.06 x0.04* 2 2 320x10°m* © D.LDUNN 20Now calculate f gET os 28s ate X320X10" 9 aoe £=0.56,/2— wh? 65.92 x2" ‘The natural frequency due to the mass at the end only is f = Evaluate this. [as x 78x10" 2n\ 5x2" Now find the combined frequency. 1 1 8532? 27.54 1° =66.42 £=8.15Hz WORKE! EXAMPLE No.8 A light shaft carries a pulley at the centre with its centre of gravity on the centre line. The shaft is supported in self aligning bearing at the ends. The shaft deflects 0.5 mm under the static weight of the shaft, Determine the lowest natural frequency. UTION ‘The frequency corresponds to that of a simply supported beam so fe at Je at | 8! Loo si 2m Vx, 2m V 0.0005 WORKED EXAMPLE No.9 The shaft in the previous example has a distributed weight of 40 N/m and the bearings are 2 m apart. The flexural stiffness is 4500 N m’. What is the natural frequency? SOLUTION The frequency due to the distributed weight only is * [SEL _) 579 Bx 4500 463 40x12 2 VwL* The critical speed is N is found from 1 1 1 1 223° 36.3 Hence f= 19 Hz © D.LDUNN 2WORKED EXAMPLE No.10 A steel wire 2 mm diameter is held between chucks Im apart. The wire weighs 0.241 N/m. The flexural stiffness is 0.157 N m’. Calculate the first and second mode frequencies. SOLUTION The lowest frequency due to the distributed weight only is BET _ 5 56p [S1XO.IST _ 9, 0.241x1 wL ray ray (04) [8 #1) PER oan 2) Vwit 20 2) V 0241x1' £=3,562 The second mode frequency is f LF. ENT EXERCISE Ja. A shaft is 20 mm diameter and 3 m long and simply supported at each end. The density is 7 830 kgim3, E = 205 GPa, Calculate the natural frequency. (Ans. 4.47 Tz) The same shaft has a flywheel mounted at the middle of mass 10 kg 1b. Calculate the new natural frequency due to the point load only. (Ans. 2.69 Hz) 1c. Calculate the combined natural frequency for the uniform and point load. (Ans. 2.3 Hz) 2. A shaft 1.8 m long and 30 mm diameter has an moment of inertia of 2.5 kg m? at one end and and 4 kg m? at the other end, The shaft rests in bearings. The modulus of rigidity is 80 GN/m’. Calculate the natural frequency of torsional vibrations and the position of the node. (Answers 7.628 Hz and 1.108 m from the smaller inertia) 3. A shaft 30 mm diameter is made from steel with density 7830 kg/m*. Calculate the weight per metre (54.3 N/m) The shaft runs between self aligning bearings 2.2 m apart. The modulus of elasticity E is 200GN/m*, Calculate the lowest natural frequeney. (12.3 Tz) A pulley is mounted at the middle and this makes the shaft deflect 1.5 mm. Calculate the lowest natural frequency. (12.7 Hz) 5. A thin rod 5 mm diameter is held between two chucks 0.8 m apart. The wire weighs 1.508 N/m and the flexural stiffiness is 6.136 N m”, Calculate the first and second mode frequency. (35.2 Hz and 96.9 Hz) 6. A shaft is 50 mm diameter and 8 m long and may be regarded as simply supported. The density is 7 830 kg/m3, E = 205 GPa. Calculate the first three critical frequencies. (1.571 , 6.279 and 14.13 Hz) 7. An aluminium rod is held in a chuck with the other end unsupported. It is 12 mm diameter and 400 mm long. The density of aluminium is 2710 kg/m3 and the modulus of elasticity E is 71 GPa, Calculate the first two critical frequencies. (54.2 Hz and 340 Hz) © D.LDUNN 25. TORSIONAL OSCILLATIONS, A torsional oscillation occurs when a restoring torque acts on a body which is displaced by turning it about its axis, The restoring torque may be caused by twisting a shaft or by some other form of spring. A torsional spring produces a torque Fixed end (1) directly proportional to the angle of twist (Q) sueh that T = k; 0. , is the torsional spring constant and for a simple shaft this may be derived in terms of the material dimensions and properties. From the tutorial on torsion we found that: T _ © tom which the torsional spring constant is k, T_GS JL 6 ob Remember that J is the polar second moment of area, L is the length and G is the modulus of rigidity for the material. For a solid round shaft J = xD¥/32. Consider a flywheel of moment of inertia Ton the end of a vertical shaft. Figure 24 Ifa torque T is applied to the flywheel it must 1) Overcome the inertia (inertia torque) by the law T= Ia 2) Twist the shaft (spring torque) by the law T= ky @ Hence T=la+ Ko ‘When the flywheel is released, the torque is zero so =a. + Ky @ k +9 Rearrange to make the acceleration (a) the subject @ = — This shows that the angular acceleration is directly proportional to displacement angle and always towards the rest position so the motion is S.H.M. The constant of proportionality o” hence: o 1 an In (ORKED EXAMPLE No.11 ‘A vertical shaft 20 mm diameter and 0.8 m long is made from metal with a modulus of rigidity G of 80 GPa. The shaft has a flywheel of moment of inertia 30kg m2 mounted on the end. Determine the frequency of the torsional oscillations. SOLUTION nD‘ _2x0.024 J 5.7.x 10° m* 32 32 k, = = G4 S00 ISTO" _ 1579.8 Nmiradian oo 08 .152Hz ¢— 1 [ke fis70.8 2mVi IV 30 © D.LDUNN 23WORKE! AMPLE No.12 A large drum of radius 0.5 m is mounted on a horiz with a mass of 400 kg on one end. The other end is ‘Nimm. The drum has a moment of inertia of 60 kg m system. The belt does not slip on the drum, Also dete in the belt if the maximum displacement of the mass i Figure 25 T. SOLUTION a The solution must be done analytically since it is a combination of torsional and linear oscillations. Suppose the mass is pulled downwards with a force F. This must overcome the inertia of the mass, the inertia of the drum and stretch the spring. Inertia force to accelerate the drum The Torque required to overcome the inertia of the drum is T = Ta. Torque = Force x radius or T= F R and the force is F =T/R Substitute T = lo. Fu = [a/R where a: is the angular acceleration of the drum. Inertia force to accelerate the mass Fa = ma where ais the linear acceleration, Force to stretch the spring F, =k x where k is the spring stiffness. Force balance F=Fu+Fatk F=Ia/R + ma + kx The angular acceleration is linked to the linear acceleration by q = a/R where R is the drum radius, F = Ia/R2 + ma + kx = a(m + Ir?) + kx Fora free oscillation F = 0 hence ( ) la (1...) getmte Makea thesubject a=-| Ix R AR I ' sim KR?) Chis shows that the acceleration is directly proportional to displacement so the motion must be simple harmonic. The constant of proportionality is the angular frequency squared so: © D.LDUNN 24Evaluate f using k = 300 000 N/m, m =400 kg ,1= 6 0 kg m? and R = 0.5 m 1 |} 300000 | 60 3.44 Hz +400 Consider the force in the belt on the side with the spring, ‘The static force when at rest = weight = 400 x 9.81 = 3 924N When the oscillation changes the length of the spring by 5 mm this kx = 300 x 5=1500N is reduced or increased by ‘The minimum force is hence 3 924 - 1500 = 2.424 N ‘The maximum is 3 924 + 1500 = 5 424.N ‘Now consider the side with the mass. Static force when at rest = 3 924 N The maximum acceleration occurs when the mass is about to reverse direction and this is when x=Smm, From the studies of s.h.m, we know a= 2x = 21.652 x 0.005 = 2.34 mis? Force in belt to accelerate the mass = ma = 400 x 2.34 = 937 N Force in belt when maximum acceleration is upwards is 3 924 + 937 =4 861 N Force in belt when acceleration is downwards is 3 924 - 937= 2987 N ‘The maximum and minimum forces in the belt are 5 424 N and 2 424 N respectively. SELF ASSESSMENT EXERCISE No.4 1. A large drum of radius 0.12 m is mounted on a horizontal shafi. A belt runs over it as shown in fig.19 with a mass of 10 kg on one end. The other end is restrained with a spring of stiffness 788 Nim. ‘The drum has a moment of inertia of 0.04 kg m2. Determine the following, i. The natural frequency of the system assuming the belt does not slip on the drum. ii, The maximum and minimum forces in the belt on each side given the maximum displacement is 0.1 m. Ans 1.25 Hz, 160 N and 36.4 N, 176.9 N and 19.3 N. 2. A wheel 0.6 m diameter has a mass of 8 kg and is mounted on a horizontal axle. The wheel is perfectly balanced about the axle, A small mass of 200 g is fixed to the outer edge of the wheel. The wheel is the made to oscillate with a small amplitude and the periodic time is found to be 5.9 s. Determine the radius of gyration of the wheel. (Ans. 0.255 m) ‘The maximum amplitude of oscillation is 10°, Find the maximum angular velocity of the wheel as it oscillates. (Ans. 0.186 rad/s) © D.LDUNN 253. A pendulum is shown in fig.20. It consists of a mass of 5 kg on a weightless rod. The centre of gravity G is 400 mm below the pivot. The radius of gyration about G is 0.05 m. Two light springs are attached to the rod at point A as shown each with stiffness 750 N/m. Point A is 120 mm below the pivot. springs Mass Figure 26 The rod is at rest in the vertical position. Determine the frequency of oscillation if it is disturbed slightly. Note that you need the parallel axis theorem to do this. (Ans 1.13 Hz). The angular displacement of the rod is 2° when rel acceleration. (Ans, .77 rad/s?) sed, Determine the maximum angular 4, A mass of 4 kg is suspended from a spring and oscillates up and down at 2 Hz, Determine the stiffness of the spring, (Ans, 631.6 N/m) From recordings made of the oscillation it is found that at a certain point in the motion the velocity is 0.3 m/s and the displacement is 20 mm, both being positive downwards in direction. Determine the amplitude of the motion and the maximum velocity and acceleration, Write down the equations of motion. Note that the data given is at time t = 0. You will have to assume that X= xgcos(at + 6) at time 0 Ans. x= 0.0311 cos(wt + 50°) v=0.3914 sin(wt + 50°) a=-157.9x © D.LDUNN 266 OSCILLATION OF A FLOA’ ING BODY You may have observed that some bodies floating in water bob up and down. This is another example of simple harmonic motion and the Force restoring force in this case is buoyaney. { Consider a floating body of mass M kg. Initially it is at rest and all the forces acting on it add up to zero. Suppose a force F is applied to the top to push it down a distance x. The applied force F must overcome this buoyancy force and also overcome the inertia of the body. Area position Figure 27 BUOYANCY FORCE The pressure on the bottom increases by Ap = p gx. The buoyancy force pushing it up on the bottom is F, and this increases by Ap A. Substitute for Ap and F,=p gx A INERTIA FORCE The inertia force acting on the body is Fi = Ma BALANCE OF FORCES The applied force must be F =F; + Fy and this must be zero if the vibration is free. O=F, +R, =Ma+pgxA Make the acceleration the subject. pAg This shows that the acceleration is directly proportional to displacement and is always directed towards the rest position so the motion must be simple harmonic and the constant of proportionality must be the angular frequency squared. > _ PAY ot - PAS SELF A! MENT EXERCISE. 1. A cylindrical rod is 80 mm diameter and has a mass of 5 kg. It floats vertically in water of density 1036 kg/m’, Calculate the frequency at which it bobs up and down. (Ans, 0.508 Hz) © D.LDUNN 21