VTU-NPTEL-NMEICT Project
|
|
- Theodore Potter
- 5 years ago
- Views:
Transcription
1 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid NMEICT, MHRD, New Dehi SME Name : Course Name: Type of the Course Modue Subject Matter Expert Detais Dr.MOHAMED HANEEF PRINCIPA, VTU SENATE MEMBER Vibration engineering web X VTU-NPTE-NMEICT Project DEPARTMENT OF MECHANICA ENGINEERING, GHOUSIA COEGE OF ENGINEERING, RAMANARAM Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 1 of 4
2 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 CONTENTS S. No. DISCRETION 1. ecture Notes (Continuous System: Approximate method).. Quadrant a. Animations. b. Videos. c. Iustrations. a. Wikis. b. Open Contents a. Probems. Quadrant -3 VTU-NPTE-NMEICT Project Quadrant -4 b. Sef Assigned Q & A.. Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page of 4
3 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Modue-X Continuous system : Approximate Soution ECTURE NOTES 1. Rayeigh-Ritz method Rayeigh-Ritz method is an extension of Rayeigh s method. It not ony provides a means of obtaining a more accurate vaue for the fundamenta frequency, but it aso gives approximations to the higher frequencies and mode shapes. In Rayeigh s method a singe shape function is used to approximate the natura mode of vibration but in Rayeigh-Ritz s method a series of shape functions mutipied by constant coefficients are used. The coefficients are adjusted by minimizing the frequency with respect to each of the coefficients, which resuts in n agebraic equations in ω. The soution of these equations then gives the natura frequencies and mode shapes of the system. The number frequencies cacuated is equa to the number arbitrary functions used. If n functions are seected for approximating the defection y(x) for the transverse vibration of a uniform beam, then. y(x) = c 1 y 1 (x) + c y (x) + c 3 y 3 (x) c n y n (x) (1) where y 1 (x), y (x), y 3 (x), y n (x) are the constant coefficients to be determined. We have, ω n = Maximum potentia energy (U max) Maximum kinetic energy (T max ) = EI d i y (x) dx dx = d EI i y dx dx ρa{y(x)} dx my dx () Substituting equation (1) in equation () and equating the first partia derivative of the resuting expression with respect to each of the constant coefficients C i to zero wi make the natura frequency stationary. i.e., ω n c 1 = ; ω n c VTU-NPTE-NMEICT Project = ; ω n c 3 = In genera ω n c i = where I = 1,, n (3) Equation 3 represents a set of n agebraic equation in the n coefficients c 1, c, c 3, ---- c n and ω n to be determined. Since the right hand side of equation is a function of c 1, c, c 3, ---- c i, its partia derivatives with respect to c 1, c, c 3, ---- c i wi a be zero for minimum vaue of ω n. Hence we have Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 3 of 4
4 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 my dx EI d y c i dx dx EI d y dx dx my c i dx = where i = 1,, 3, n (4) From equation () Therefore equation (4) becomes where A = EI d y EI d y dx EI d y c i dx dx dx dx dx = ω n my dx ω n my dx = (5) c i B i.e., A ω c n = (6) i c i and B = my dx Equation 6 are set of I homogeneous inear equations in unknowns c 1, c, c 3, ---- c n. It wi have a non-trivia soution if the determinant of its coefficients is equa to zero. It gives an i th degree equation in ω n, the owest root of which is equa to ω n1 and the next higher root ω n and so on.. Rayeigh s energy method: Using the energy method to find the natura frequencies of singe degree of freedom system, the principe of conservation of energy, in the context of an undamped vibrating system, can be restated as T 1 + U 1 = T + U (1) Where the subscripts 1 and denote two different instants of time; Specificay, we use the subscript 1 to denote the time when the mass is passing through its static equiibrium position and choose u 1 = as reference for the potentia energy. If we et the subscript indicate the time corresponding to the maximum dispacement of the mass, we have T =. Thus Eq. [1] Becomes. VTU-NPTE-NMEICT Project T 1 +O = O + U --- () If the system is undergoing harmonic motion, then T1 and U denote the maximum vaues of T and U, respectivey, and Eq. [] becomes. T max = U max --- (3) Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 4 of 4
5 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING Ritz-Gaerkin method In the Ritz-Gaerkin method, an approximate soution of the probem is found by satisfying the governing noninear equation in the average. To see how the method works, et the noninear differentia equations be represented as E[x] = (1) An approximate soution of Eq. (1)is assumed as x (t) = a 1 1 (t) + a (t) + + a n n (t) () Where 1 (t), (t). n (t) are prescribed functions of time and a 1, a,. a n are weighting factors to be determined. If Eq. (), we get a function E x (t). since x (t) is not, in genera, the exact soution of Eq. (1), E (t) = E x (t)wi not be zero. However, the vaue of E (t)wi serve as a measure of the accuracy of the approximation; in fact, E (t) as x x. The weighting factor a i are determined by minimizing the integra τ E [t]dt (3) Where τ denoted as the period of the motion. The minimization of the function of Eq.(3) requires τ E τ E [t] [t]dt = E [t] dt = (4) a i a i Eq. (4) represents a system of n agebraic equations can be soved simutaneousy to find vaues of a 1, a. a n. VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 5 of 4
6 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Animations QUADRANT- (Animation inks reated, Continuous System: Approximate method) gation.htm Videos (Animation inks reated, Continuous System: Approximate method) rayeigh-s-energy-method gaerkin-s-method-1d-finite-eement-method VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 6 of 4
7 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 IUSTRATIONS 1. Expain Rayeigh Energy Method. So: Rayeigh-Ritz method is an extension of Rayeigh s method. It not ony provides a means of obtaining a more accurate vaue for the fundamenta frequency, but it aso gives approximations to the higher frequencies and mode shapes. In Rayeigh s method a singe shape function is used to approximate the natura mode of vibration but in Rayeigh-Ritz s method a series of shape functions mutipied by constant coefficients are used. The coefficients are adjusted by minimizing the frequency with respect to each of the coefficients, which resuts in n agebraic equations in ω. The soution of these equations then gives the natura frequencies and mode shapes of the system. The number frequencies cacuated is equa to the number arbitrary functions used. If n functions are seected for approximating the defection y(x) for the transverse vibration of a uniform beam, then. y(x) = c 1 y 1 (x) + c y (x) + c 3 y 3 (x) c n y n (x) (1) where y 1 (x), y (x), y 3 (x), y n (x) are the constant coefficients to be determined. We have, ω n = Maximum potentia energy (U max) Maximum kinetic energy (T max ) = EI d i y (x) dx dx = d EI i y dx dx ρa{y(x)} dx my dx () Substituting equation (1) in equation () and equating the first partia derivative of the resuting expression with respect to each of the constant coefficients C i to zero wi make the natura frequency stationary. i.e., ω n c 1 = ; ω n c = ; ω n = c 3 In genera ω n c i = where I = 1,, n (3) Equation 3 represents a set of n agebraic equation in the n coefficients c 1, c, c 3, ---- c n and ω n to be determined. Since the right hand side of equation is a function of c 1, c, c 3, ---- c i, its partia derivatives with respect to c 1, c, c 3, ---- c i wi a be zero for minimum vaue of ω n. Hence we have my dx EI d y c i dx dx EI d y dx dx my c i dx = where i = 1,, 3, n (4) From equation () Therefore equation (4) becomes VTU-NPTE-NMEICT Project EI d y dx dx = ω n my dx Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 7 of 4
8 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 where A = EI d y EI d y c i dx dx dx dx ω n my dx = (5) c i i.e., A ω B c n = (6) i c i and B = my dx Equation 6 are set of I homogeneous inear equations in unknowns c 1, c, c 3, ---- c n. It wi have a non-trivia soution if the determinant of its coefficients is equa to zero. It gives an i th degree equation in ω n, the owest root of which is equa to ω n1 and the next higher root ω n and so on.. Expian the Concept of Rayeigh Energy Method. Soution) Using the energy method to find the natura frequencies of singe degree of freedom system, the principe of conservation of energy, in the context of an undamped vibrating system, can be restated as T 1 + U 1 = T + U (1) Where the subscripts 1 and denote two different instants of time; Specificay, we use the subscript 1 to denote the time when the mass is passing through its static equiibrium position and choose u 1 = as reference for the potentia energy. If we et the subscript indicate the time corresponding to the maximum dispacement of the mass, we have T =. Thus Eq. [1] Becomes. T 1 +O = O + U --- () If the system is undergoing harmonic motion, then T1 and U denote the maximum vaues of T and U, respectivey, and Eq. [] becomes. T max = U max --- (3) VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 8 of 4
9 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING Expain the concept of Ritz Gaerkins Method. Soution) In the Ritz-Gaerkin method, an approximate soution of the probem is found by satisfying the governing noninear equation in the average. To see how the method works, et the noninear differentia equations be represented as E[x] = (1) An approximate soution of Eq. (1)is assumed as x (t) = a 1 1 (t) + a (t) + + a n n (t) () Where 1 (t), (t). n (t) are prescribed functions of time and a 1, a,. a n are weighting factors to be determined. If Eq. (), we get a function E x (t). since x (t) is not, in genera, the exact soution of Eq. (1), E (t) = E x (t)wi not be zero. However, the vaue of E (t)wi serve as a measure of the accuracy of the approximation; in fact, E (t) as x x. The weighting factor a i are determined by minimizing the integra τ E [t]dt (3) Where τ denoted as the period of the motion. The minimization of the function of Eq.(3) requires τ E τ E [t] [t]dt = E [t] dt = (4) a i a i VTU-NPTE-NMEICT Project Eq. (4) represents a system of n agebraic equations can be soved simutaneousy to find vaues of a 1, a. a n. 4. Determine the first two natura frequencies of fixed-fixed uniform string of mass density P per unit ength stretched between x = and x = with an initia tension T. The defection functions are given by y1(x) = x( x) and y(x) = x ( x) Soution : We have, y (x) = c 1 y 1 (x) + c y (x) = c 1 x( x) + c x ( x) = c 1 x c 1 x +c x + c x 4 c x 3 dy = c 1 c 1 x + c x + 4c x 3 6c x dx Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 9 of 4
10 c c MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 For a string, Max. PE = U max = 1 Max. KE = T max = ω n T dy et A = T c 1 3 dx Tρy 3 B = ρ c 1 5 Equating the equations (1) and () we get Considering, ω n c 1 = ω n A = Bω n ω n = A B c = A A B = A c 1 B c 1 A A B = A c B c Now, A dx = c dx = c c 7 + c 1c c 9 + c 1c c 7 + c 1c 5 T c 9 + c 1c 7 ρω n (1) () (3) (4) ω B c n = (5) 1 c 1 ω B c n = (6) c = T c 1 c c B = ρ c c c A = T c 4c 5 + c B = ρ c c 9 + c Therefore A ω B c n = gives 1 c 1 T c c 5 3 i.e., c 1 T3 ρω n i.e., c 1 T5 VTU-NPTE-NMEICT Project 15 ω n ρ c R 15 + c 7 7 = T5 ρω n 7 15 Simiary, A c ω n B c = gives T 4c 5 + c 1 5 ω n ρ c 9 + c 1 7 = ρω n R T7 ρω n = (7) 63 = (8) Equations (7) and (8) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 1 of 4
11 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 1 μ 1 μ μ μ = (9) where µ = ρ ω n T Expanding (9) we get, µ - 11µ + 18 = ω n1 = 3.14 T rad ρ sec From equation (1), µ = +11± µ 1 = 9.87 and µ = 1.13 and ω n = 1.11 T rad ρ sec 5. Find the fundamenta natura frequency of transverse vibration of a tapered bar fixed at its base as shown in figure beow. Consider the width of this bar as unity. Soution : Thickness of bar at distance x from base t = T 1 x VTU-NPTE-NMEICT Project Moment of inertia at distance x from base I(x) = 1 bt3 Mass per unit ength at distance x from base. = 1 1 T3 1 x 3 (since b = 1) m(x) = ρt = ρt1 x we have y(x) = c 1 y 1 (x) + c y (x) c n y n (x) x Assume y = c 1 + c x 3 + c x in order that the foowing 4 boundary conditions are satisfied. Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 11 of 4
12 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Tria : 1 We know, y(o) = ; dy () dx Consider y = c 1 x Soving the equation (1) we get, Tria : We know, dy dx = c 1x and d y dx = c 1 ω n = EI(x)dy ω n1 = 1.58 T rad = ; EI () d y () dx = ; EI () d3 y () dx 3 = dx dx = 1 E 1 T3 1 x 3 c 1 m(x)y dx ρt1 x x c 1 dx sec To get more accurate resuts consider first two terms. i.e., y = c 1 x + c dy = c 1x + 3c x dx 3 d y = c 1 + 6c x dx 3 EI(x) d y dx c i dx i.e., E 1 c i 1 T 1 x 3 c 1 + 6c x 3 x 3 3 ω n m(x)y dx = c i dx ω n ρt 1 x c i x x c c Integrating first and then differentiating with the respect to c 1 and c separatey, we get, (.8µ)c 1 + (1..57µ)c = (1..57µ) + (1..43µ)c = VTU-NPTE-NMEICT Project where + µ = ρ4 ωn T E (1) 3 dx = Equations () and (3) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence.8µ 1..57µ 1..57µ 1..43µ = Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 1 of 4
13 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 i.e., (.8µ) (1..43µ) - (1..57µ) (1..57µ) = i.e.,.4.86µ -.96µ +.344µ µ µ = i.e.,.191µ.45µ +.96 = i.e., µ 3.665µ = For the fundamenta natura frequency, i.e., µ = 3.665± µ 1 =.36 and µ = 1.31 µ 1 =.36 = ρ4 ωn T E ω n1 = T rad/sec It is more accurate than the tria (1) vaue. 6. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the torsiona system of figure beow. Use cubic poynomias as tria functions. J = m 4 G = N m The Rayeigh-Ritz method can be appied using tria functions satisfying ony the geometric boundary conditions (i.e., boundary conditions deveoped soey from geometric considerations). Poynomias satisfying the boundary conditions for a fixed-free shaft are φ i (x) = x 3 3 x, φ (x) = x x The coefficients used in Equation (7.) are α ij = JG d i d j dx + k dx dx t i () j () VTU-NPTE-NMEICT Project β ij = ρj i (x) j (x)dx Kt = N m r Evauation of these coefficients eads to α 11 = JG(3x 3 ) dx + k t ( ) = α 1 = JG(3x 3 )(x ) dx + k t ( )( ) = α 13 = JG(x ) dx + t ( ) = β 11 = ρj (x ) dx =.589 β 1 = ρj (x 3 3 x)(x x)dx =.93 Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 13 of 4
14 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 β = ρj (x x) dx =.3159 Equation (7.) become ( ω )C 1 + ( ω )C = ( ω )C 1 + ( ω )C = A nontrivia soution is obtained if and ony if the determinant of the coefficient matrix is set to zero, eading to.1369ω = whose soutions are ω 1 = 886 rad/s, ω = 6643 rad/s, 7. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of a uniform fixed-fixed beam. Use the foowing tria functions, which satisfy a boundary conditions: = x 4 x 3 + x, φ (x) = x 5 3 x x The Rayeigh-Ritz approximation to the mode shape is w(x) = C 1 φ 1 (x) + C φ (x) The appropriate form of the coefficients for equation (7.) is α ij = EI d i d j dx dx dx β ij = ρa i (x) j (x)dx Using the suggested tria functions, the coefficients are cacuated as α 11 = EI(1x 1x + ) dx =.8EI 5 α 1 = EI(1x 1x + )(x 3 18 x ) dx = EI 6 α = EI(x 3 18 x ) dx = 5.148EI 7 VTU-NPTE-NMEICT Project β 11 = ρa (x 4 x 3 + x ) dx =.1587ρA 9 β 1 = ρa (x 4 x 3 + x )(x 5 3 x + 3 x )dx =.3968ρA 1 β = ρa (x 5 3 x x ) dx =.99567ρA 11 Substitution into equation (7.) and rearrangement yieds ( φ)C 1 + (.3968φ)C = (.3968φ)C 1 + ( φ)C = where φ = ω ρa4 EI Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 14 of 4
15 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 A nontrivia soution of the above system exists if and ony if the determinant of the coefficient matrix is set to zero. To this end ( φ)( φ) (.3968φ) = φ φ = The soutions of the above equations are 54.9 and 48.3 eading to ω 1 =.4 EI ρa 4, ω = EI ρa 4 8. The exhaust from a singe-cyinder four-stroke diese engine is to be connected to a siencer, and the pressure therein is to be measured with a simpe U-tube manometer [see fig. beow]. Cacuate the minimum ength of the manometer tube so that the natura frequency of osciation of the mercury coumn. Wi be 3.5 times sower than the frequency of the pressure fuctuation in the siencer at an engine speed of 6 rpm. The frequency of pressure fuctuation in the siencer is equa to. Number of cyinders x speed of the engine x Figure: U-tube manometer VTU-NPTE-NMEICT Project x SOUTION 1. Natura frequency of osciation of the iquid coumn: et the datum in fug..18 be taken as the equiibrium position of the iquid. If the dispacement of the iquid coumn from the equiibrium position is denoted by x, the change in potentia energy is given by potentia energy of raised potentia energy of depressed U = + iquid coumn iquid coumn Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 15 of 4
16 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 U = weight mercury raised x dispacement of the C. G of the segment weight of mercury depressed + x dispacement of the C. G. of the segment = [Axγ] x + [Axγ] = [Ax γ] (1) Where A is the cross-sectiona area of the mercury coumn and y is the specific weight of mercury. The change in kinetic energy is given by T = 1 [mass of mercury] [veocity] = 1 Aγ g x () Where is the ength of the mercury coumn. By assuming harmonic motion, we can write where x is the maximum dispacement and ω n is the natura frequency. x(t) = Xcosω n t (3) By substituting Eq. [3] into Eqs. [1] and [], we obtain U = U max cos ω n t T = T max sin ω n t U max = AγX (4) (5) (6) T max = 1 Aγω n X (7) g By equating U max to T max we obtain the norma frequency: VTU-NPTE-NMEICT Project ω n = g (8). ength of mercury coumn: the frequency of pressure fuctuation in the siencer = 1 6 = 3 rpm 3 π = = 1π 6 rad/sec Thus the frequency of osciation of iquid coumn in the manometer is 1 π /3.5 = 9. rad/sec. by using Eq. [8], we obtain Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 16 of 4
17 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 ω n = g = 9. = =. 43 m VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 17 of 4
18 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 QUADRANT-3 Wikis: (This incudes wikis reated to Continuous System: Approximate Method) Open Contents: (This incudes wikis reated to Continuous System: Approximate Method) Mechanica Vibrations, S. S. Rao, Pearson Education Inc, 4 th edition, 3. Mechanica Vibrations, V. P. Singh, Dhanpat Rai & Company, 3 rd edition, 6. Mechanica Vibrations, G. K.Grover, Nem Chand and Bros, 6 th edition, 1996 Theory of vibration with appications,w.t.thomson,m.d.daheh and C Padmanabhan,Pearson Education inc,5 th Edition,8 Theory and practice of Mechanica Vibration : J.S.Rao&K,Gupta,New Age Internationa Pubications,New Dehi,1 VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 18 of 4
19 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 QUADRANT-4 Probems 1. Determine the first two natura frequencies of fixed-fixed uniform string of mass density P per unit ength stretched between x = and x = with an initia tension T. The defection functions are given by y1(x) = x( x) and y(x) = x ( x) Soution : We have, y (x) = c 1 y 1 (x) + c y (x) = c 1 x( x) + c x ( x) For a string, Max. PE = U max = 1 Max. KE = T max = ω n = c 1 x c 1 x +c x + c x 4 c x 3 dy = c 1 c 1 x + c x + 4c x 3 6c x dx T dy dx dx = c 1 3 Tρy dx = c 1 5 et A = T c B = ρ c 1 5 Equating the equations (1) and () we get Considering, ω n c 1 = ω n A = Bω n ω n = A B c = A A B = A c 1 B c 1 A A B = A c B c + c 7 + c 1c c 9 + c 1c c 7 + c 1c 5 T c 9 + c 1c 7 ρω n VTU-NPTE-NMEICT Project Now, A (1) () (3) (4) ω B c n = (5) 1 c 1 ω B c n = (6) c = T c 1 c c 5 15 B = ρ c c c A = T c 4c 5 + c B = ρ c c 9 + c Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 19 of 4
20 c c MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 T c i.e., c 1 T3 i.e., c 1 T5 Therefore A c 1 ω n B c 1 = gives + c 5 15 ω n ρ c c 7 ρω n R 7 = T5 ρω n 7 15 Simiary, A c ω n B c = gives T 4c 5 + c 1 5 ω n ρ c 9 + c 1 7 = ρω n R T7 ρω n = (7) 63 = (8) Equations (7) and (8) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence 1 μ 1 μ μ = (9) μ where µ = ρ ω n Expanding (9) we get, µ - 11µ + 18 = ω n1 = 3.14 T rad ρ sec T From equation (1), µ = +11± and ω n = 1.11 T 1 µ 1 = 9.87 and µ = 1.13 rad ρ sec. Find the fundamenta natura frequency of transverse vibration of a tapered bar fixed at its base as shown in figure beow. Consider the width of this bar as unity. VTU-NPTE-NMEICT Project Soution : Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page of 4
21 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Thickness of bar at distance x from base t = T 1 x Moment of inertia at distance x from base I(x) = 1 bt3 Mass per unit ength at distance x from base. = 1 1 T3 1 x 3 (since b = 1) m(x) = ρt = ρt1 x we have y(x) = c 1 y 1 (x) + c y (x) c n y n (x) x Assume y = c 1 + c x 3 + c x in order that the foowing boundary conditions are satisfied. Tria : 1 We know, y(o) = ; dy () Consider y = c 1 x Soving the equation (1) we get, Tria : We know, dy = c 1x dx dx and d y ω n = EI(x)dy ω n1 = 1.58 T rad = ; EI () d y () dx = ; EI () d3 y () dx 3 = = c 1 dx dx dx = 1 E 1 T3 1 x 3 c 1 m(x)y dx ρt1 x x c 1 dx sec To get more accurate resuts consider first two terms. VTU-NPTE-NMEICT Project x i.e., y = c 1 + c dy dx = c 1x + 3c x 3 d y = c 1 + 6c x dx 3 EI(x) d y dx c i dx x 3 3 ω n m(x)y dx = c i (1) Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 1 of 4
22 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 i.e., E 1 c i 1 T 1 x 3 c 1 + 6c x 3 dx ω n ρt 1 x c i x x c c Integrating first and then differentiating with the respect to c 1 and c separatey, we get, (.8µ)c 1 + (1..57µ)c = (1..57µ) + (1..43µ)c = where + µ = ρ4 ωn T E 3 dx = Equations () and (3) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence.8µ 1..57µ 1..57µ 1..43µ = i.e., (.8µ) (1..43µ) - (1..57µ) (1..57µ) = i.e.,.4.86µ -.96µ +.344µ µ µ = i.e.,.191µ.45µ +.96 = i.e., µ 3.665µ = For the fundamenta natura frequency, i.e., µ = 3.665± µ 1 =.36 and µ = 1.31 µ 1 =.36 = ρ4 ωn T E ω n1 = T rad/sec It is more accurate than the tria (1) vaue. VTU-NPTE-NMEICT Project 3. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the torsiona system of figure beow. Use cubic poynomias as tria functions. J = m 4 G = N m Kt = N m r The Rayeigh-Ritz method can be appied using tria functions satisfying ony the geometric boundary conditions (i.e., boundary conditions deveoped soey from geometric considerations). Poynomias satisfying the boundary conditions for a fixed-free shaft are φ i (x) = x 3 3 x, φ (x) = x x The coefficients used in Equation (7.) are Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page of 4
23 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 α ij = JG d i dx β ij = ρj i (x) j (x)dx d j dx dx + k t i () j () Evauation of these coefficients eads to α 11 = JG(3x 3 ) dx + k t ( ) = α 1 = JG(3x 3 )(x ) dx + k t ( )( ) = α 13 = JG(x ) dx + t ( ) = β 11 = ρj (x ) dx =.589 β 1 = ρj (x 3 3 x)(x x)dx =.93 β = ρj (x x) dx =.3159 Equation (7.) become ( ω )C 1 + ( ω )C = ( ω )C 1 + ( ω )C = A nontrivia soution is obtained if and ony if the determinant of the coefficient matrix is set to zero, eading to.1369ω = whose soutions are ω 1 = 886 rad/s, ω = 6643 rad/s, 4. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of a uniform fixed-fixed beam. Use the foowing tria functions, which satisfy a boundary conditions: = x 4 x 3 + x, φ (x) = x 5 3 x x The Rayeigh-Ritz approximation to the mode shape is w(x) = C 1 φ 1 (x) + C φ (x) The appropriate form of the coefficients for equation (7.) is VTU-NPTE-NMEICT Project α ij = EI d i dx d j dx dx β ij = ρa i (x) j (x)dx Using the suggested tria functions, the coefficients are cacuated as α 11 = EI(1x 1x + ) dx =.8EI 5 α 1 = EI(1x 1x + )(x 3 18 x ) dx = EI 6 Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 3 of 4
24 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 α = EI(x 3 18 x ) dx = 5.148EI 7 β 11 = ρa (x 4 x 3 + x ) dx =.1587ρA 9 β 1 = ρa (x 4 x 3 + x )(x 5 3 x + 3 x )dx =.3968ρA 1 β = ρa (x 5 3 x x ) dx =.99567ρA 11 Substitution into equation (7.) and rearrangement yieds where ( φ)C 1 + (.3968φ)C = (.3968φ)C 1 + ( φ)C = φ = ω ρa4 EI A nontrivia soution of the above system exists if and ony if the determinant of the coefficient matrix is set to zero. To this end ( φ)( φ) (.3968φ) = φ φ = The soutions of the above equations are 54.9 and 48.3 eading to ω 1 =.4 EI ρa 4, ω = EI ρa 4 5. The exhaust from a singe-cyinder four-stroke diese engine is to be connected to a siencer, and the pressure therein is to be measured with a simpe U-tube manometer [see fig. beow]. Cacuate the minimum ength of the manometer tube so that the natura frequency of osciation of the mercury coumn. Wi be 3.5 times sower than the frequency of the pressure fuctuation in the siencer at an engine speed of 6 rpm. The frequency of pressure fuctuation in the siencer is equa to. Number of cyinders x speed of the engine VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 4 of 4
25 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 x x SOUTION Natura frequency of osciation of the iquid coumn: et the datum in fug..18 be taken as the equiibrium position of the iquid. If the dispacement of the iquid coumn from the equiibrium position is denoted by x, the change in potentia energy is given by potentia energy of raised potentia energy of depressed U = + iquid coumn iquid coumn weight mercury raised U = x dispacement of the C. G of the segment Figure: U-tube manometer weight of mercury depressed + x dispacement of the C. G. of the segment = [Axγ] x + [Axγ] = [Ax γ] (1) Where A is the cross-sectiona area of the mercury coumn and y is the specific weight of mercury. The change in kinetic energy is given by VTU-NPTE-NMEICT Project T = 1 [mass of mercury] [veocity] = 1 Aγ g x () Where is the ength of the mercury coumn. By assuming harmonic motion, we can write where x is the maximum dispacement and ω n is the natura frequency. x(t) = Xcosω n t (3) By substituting Eq. [3] into Eqs. [1] and [], we obtain Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 5 of 4
26 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 U = U max cos ω n t T = T max sin ω n t U max = AγX (4) (5) (6) T max = 1 Aγω n X (7) g By equating U max to T max we obtain the norma frequency: ω n = g (8) ength of mercury coumn: the frequency of pressure fuctuation in the siencer = 1 6 = 3 rpm 3 π = = 1π 6 rad/sec Thus the frequency of osciation of iquid coumn in the manometer is 1 π /3.5 = 9. rad/sec. by using Eq. [8], we obtain ω n = g = = 9. =. 43 m VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 6 of 4
27 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Frequenty asked Questions. 1. Determine the first two natura frequencies of fixed-fixed uniform string of mass density P per unit ength stretched between x = and x = with an initia tension T. The defection functions are given by y 1 (x) = x( x) and y (x) = x ( x). Find the fundamenta natura frequency of transverse vibration of a tapered bar fixed at its base as shown in figure beow. Consider the width of this bar as unity. 3. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the torsiona system of figure beow. Use cubic poynomias as tria functions. J = m 4 G = N m VTU-NPTE-NMEICT Project Kt = N m r 4. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of a uniform fixed-fixed beam. Use the foowing tria functions, which satisfy a boundary conditions: = x 4 x 3 + x, φ (x) = x 5 3 x x Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 7 of 4
28 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the system of figure beow. Use the two owest mode shapes for a uniform fixed-free bar as tria functions. 6. Write a short Notes on i) Rayeigh Ritz method ii) iii) Rayeigh Energy Method Ritz Gaerkin Method Sef Question & Answer 1. Determine the first two natura frequencies of fixed-fixed uniform string of mass density P per unit ength stretched between x = and x = with an initia tension T. The defection functions are given by y1(x) = x( x) and y(x) = x ( x) Soution : We have, y (x) = c 1 y 1 (x) + c y (x) = c 1 x( x) + c x ( x) = c 1 x c 1 x +c x + c x 4 c x 3 dy = c 1 c 1 x + c x + 4c x 3 6c x dx For a string, Max. PE = U max = 1 Max. KE = T max = ω n T dy dx dx = c 1 3 Tρy dx = c 1 5 et A = T c B = ρ c 1 5 Equating the equations (1) and () we get A = Bω n + c 7 + c 1c c 9 + c 1c c 7 + c 1c 5 T c 9 + c 1c 7 ρω n VTU-NPTE-NMEICT Project (1) () (3) (4) Considering, ω n c 1 ω n = A B = ω n c = A A B = A c 1 B c 1 A A B = A c B c ω B c n = (5) 1 c 1 ω B c n = (6) c Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 8 of 4
29 c c MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 T c i.e., c 1 T3 i.e., c 1 T5 Now, A = T c 1 c c 5 15 B = ρ c c c A = T c 4c 5 + c B = ρ c c 9 + c Therefore A c 1 ω n B c 1 = gives + c 5 15 ω n ρ c c 7 ρω n R 7 = T5 ρω n 7 15 Simiary, A ω B c n = gives c T 4c 5 + c ρω n ω n ρ c R 14 = (7) c 1 7 T7 ρω n = 63 = (8) Equations (7) and (8) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence 1 μ 1 μ μ = (9) μ where µ = ρ ω n Expanding (9) we get, µ - 11µ + 18 = T From equation (1), µ = +11± VTU-NPTE-NMEICT Project 1 µ 1 = 9.87 and µ = 1.13 ω n1 = 3.14 T rad ρ sec and ω n = 1.11 T rad ρ sec. Find the fundamenta natura frequency of transverse vibration of a tapered bar fixed at its base as shown in figure beow. Consider the width of this bar as unity. Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 9 of 4
30 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Soution : Thickness of bar at distance x from base t = T 1 x Moment of inertia at distance x from base I(x) = 1 bt3 Mass per unit ength at distance x from base. = 1 1 T3 1 x 3 (since b = 1) m(x) = ρt = ρt1 x we have y(x) = c 1 y 1 (x) + c y (x) c n y n (x) x Assume y = c 1 + c x 3 + c x in order that the foowing boundary conditions are satisfied. Tria : 1 We know, y(o) = ; dy () Consider y = c 1 x dx dy dx = c 1x and d y dx = c 1 = ; EI () d y () dx = ; EI () d3 y () dx 3 = VTU-NPTE-NMEICT Project Soving the equation (1) we get, dx dx = 1 E 1 T3 1 x 3 c 1 m(x)y dx ρt1 x x c 1 dx ω n = EI(x)dy (1) ω n1 = 1.58 T rad sec Tria : Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 3 of 4
31 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 We know, To get more accurate resuts consider first two terms. i.e., y = c 1 x + c dy dx = c 1x + 3c x 3 d y = c 1 + 6c x dx 3 EI(x) d y dx c i dx i.e., E 1 c i 1 T 1 x 3 c 1 + 6c x 3 x 3 3 ω n m(x)y dx = c i dx ω n ρt 1 x c i x x c c Integrating first and then differentiating with the respect to c 1 and c separatey, we get, (.8µ)c 1 + (1..57µ)c = (1..57µ) + (1..43µ)c = where + µ = ρ4 ωn T E 3 dx = Equations () and (3) are homogeneous inear equations in unknowns c 1 and c. For a non trivia soution of c 1 and c, the determinant of the coefficients of c 1 and c must be zero. Hence.8µ 1..57µ 1..57µ 1..43µ = i.e., (.8µ) (1..43µ) - (1..57µ) (1..57µ) = i.e.,.4.86µ -.96µ +.344µ µ µ = i.e.,.191µ.45µ +.96 = i.e., µ 3.665µ = VTU-NPTE-NMEICT Project For the fundamenta natura frequency, i.e., µ = 3.665± µ 1 =.36 and µ = 1.31 µ 1 =.36 = ρ4 ωn T E ω n1 = T rad/sec It is more accurate than the tria (1) vaue. Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 31 of 4
32 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the torsiona system of figure beow. Use cubic poynomias as tria functions. J = m 4 G = N m The Rayeigh-Ritz method can be appied using tria functions satisfying ony the geometric boundary conditions (i.e., boundary conditions deveoped soey from geometric considerations). Poynomias satisfying the boundary conditions for a fixed-free shaft are φ i (x) = x 3 3 x, φ (x) = x x The coefficients used in Equation (7.) are α ij = JG d i d j dx + k dx dx t i () j () β ij = ρj i (x) j (x)dx Kt = N m r Evauation of these coefficients eads to α 11 = JG(3x 3 ) dx + k t ( ) = α 1 = JG(3x 3 )(x ) dx + k t ( )( ) = α 13 = JG(x ) dx + t ( ) = β 11 = ρj (x ) dx =.589 β 1 = ρj (x 3 3 x)(x x)dx =.93 β = ρj (x x) dx =.3159 Equation (7.) become ( ω )C 1 + ( ω )C = ( ω )C 1 + ( ω )C = A nontrivia soution is obtained if and ony if the determinant of the coefficient matrix is set to zero, eading to.1369ω = whose soutions are ω 1 = 886 rad/s, ω = 6643 rad/s, VTU-NPTE-NMEICT Project 4. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of a uniform fixed-fixed beam. Use the foowing tria functions, which satisfy a boundary conditions: = x 4 x 3 + x, φ (x) = x 5 3 x x Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 3 of 4
33 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 The Rayeigh-Ritz approximation to the mode shape is w(x) = C 1 φ 1 (x) + C φ (x) The appropriate form of the coefficients for equation (7.) is α ij = EI d i dx d j dx dx β ij = ρa i (x) j (x)dx Using the suggested tria functions, the coefficients are cacuated as α 11 = EI(1x 1x + ) dx =.8EI 5 α 1 = EI(1x 1x + )(x 3 18 x ) dx = EI 6 α = EI(x 3 18 x ) dx = 5.148EI 7 β 11 = ρa (x 4 x 3 + x ) dx =.1587ρA 9 β 1 = ρa (x 4 x 3 + x )(x 5 3 x + 3 x )dx =.3968ρA 1 β = ρa (x 5 3 x x ) dx =.99567ρA 11 Substitution into equation (7.) and rearrangement yieds where ( φ)C 1 + (.3968φ)C = (.3968φ)C 1 + ( φ)C = φ = ω ρa4 EI A nontrivia soution of the above system exists if and ony if the determinant of the coefficient matrix is set to zero. To this end ( φ)( φ) (.3968φ) = φ φ = The soutions of the above equations are 54.9 and 48.3 eading to ω 1 =.4 EI VTU-NPTE-NMEICT Project ρa 4, ω = EI ρa 4 5. The exhaust from a singe-cyinder four-stroke diese engine is to be connected to a siencer, and the pressure therein is to be measured with a simpe U-tube manometer [see fig. beow]. Cacuate the minimum ength of the manometer tube so that the natura frequency of osciation of the mercury coumn. Wi be 3.5 times sower than the frequency of the pressure fuctuation in the siencer at an engine speed of 6 rpm. The frequency of pressure fuctuation in the siencer is equa to. Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 33 of 4
34 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Number of cyinders x speed of the engine x x SOUTION Natura frequency of osciation of the iquid coumn: et the datum in fug..18 be taken as the equiibrium position of the iquid. If the dispacement of the iquid coumn from the equiibrium position is denoted by x, the change in potentia energy is given by potentia energy of raised potentia energy of depressed U = + iquid coumn iquid coumn weight mercury raised U = x dispacement of the C. G of the segment Figure: U-tube manometer weight of mercury depressed + x dispacement of the C. G. of the segment = [Axγ] x + [Axγ] = [Ax γ] (1) VTU-NPTE-NMEICT Project Where A is the cross-sectiona area of the mercury coumn and y is the specific weight of mercury. The change in kinetic energy is given by T = 1 [mass of mercury] [veocity] = 1 Aγ g x () Where is the ength of the mercury coumn. By assuming harmonic motion, we can write where x is the maximum dispacement and ω n is the natura frequency. x(t) = Xcosω n t (3) Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 34 of 4
35 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 By substituting Eq. [3] into Eqs. [1] and [], we obtain U = U max cos ω n t T = T max sin ω n t U max = AγX (4) (5) (6) T max = 1 Aγω n X (7) g By equating U max to T max we obtain the norma frequency: ω n = g (8) ength of mercury coumn: the frequency of pressure fuctuation in the siencer = 1 6 = 3 rpm 3 π = = 1π 6 rad/sec Thus the frequency of osciation of iquid coumn in the manometer is 1 π /3.5 = 9. rad/sec. by using Eq. [8], we obtain = ω n = g = 9. VTU-NPTE-NMEICT Project =. 43 m Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 35 of 4
36 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 Probems 1. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of the torsiona system of figure beow. Use cubic poynomias as tria functions. J = m 4 G = N m The Rayeigh-Ritz method can be appied using tria functions satisfying ony the geometric boundary conditions (i.e., boundary conditions deveoped soey from geometric considerations). Poynomias satisfying the boundary conditions for a fixed-free shaft are φ i (x) = x 3 3 x, φ (x) = x x The coefficients used in Equation (7.) are α ij = JG d i d j dx + k dx dx t i () j () β ij = ρj i (x) j (x)dx Kt = N m r Evauation of these coefficients eads to α 11 = JG(3x 3 ) dx + k t ( ) = α 1 = JG(3x 3 )(x ) dx + k t ( )( ) = α 13 = JG(x ) dx + t ( ) = β 11 = ρj (x ) dx =.589 β 1 = ρj (x 3 3 x)(x x)dx =.93 β = ρj (x x) dx =.3159 Equation (7.) become ( ω )C 1 + ( ω )C = ( ω )C 1 + ( ω )C = A nontrivia soution is obtained if and ony if the determinant of the coefficient matrix is set to zero, eading to.1369ω = whose soutions are ω 1 = 886 rad/s, ω = 6643 rad/s, VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 36 of 4
37 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14. Use the Rayeigh-Ritz method to approximate the two owest natura frequencies of a uniform fixed-fixed beam. Use the foowing tria functions, which satisfy a boundary conditions: = x 4 x 3 + x, φ (x) = x 5 3 x x The Rayeigh-Ritz approximation to the mode shape is w(x) = C 1 φ 1 (x) + C φ (x) The appropriate form of the coefficients for equation (7.) is α ij = EI d i dx d j dx dx β ij = ρa i (x) j (x)dx Using the suggested tria functions, the coefficients are cacuated as α 11 = EI(1x 1x + ) dx =.8EI 5 α 1 = EI(1x 1x + )(x 3 18 x ) dx = EI 6 α = EI(x 3 18 x ) dx = 5.148EI 7 β 11 = ρa (x 4 x 3 + x ) dx =.1587ρA 9 β 1 = ρa (x 4 x 3 + x )(x 5 3 x + 3 x )dx =.3968ρA 1 β = ρa (x 5 3 x x ) dx =.99567ρA 11 Substitution into equation (7.) and rearrangement yieds where ( φ)C 1 + (.3968φ)C = (.3968φ)C 1 + ( φ)C = φ = ω ρa4 EI A nontrivia soution of the above system exists if and ony if the determinant of the coefficient matrix is set to zero. To this end ( φ)( φ) (.3968φ) = φ φ = The soutions of the above equations are 54.9 and 48.3 eading to ω 1 =.4 EI ρa 4, VTU-NPTE-NMEICT Project ω = EI ρa 4 Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 37 of 4
38 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING The exhaust from a singe-cyinder four-stroke diese engine is to be connected to a siencer, and the pressure therein is to be measured with a simpe U-tube manometer [see fig. beow]. Cacuate the minimum ength of the manometer tube so that the natura frequency of osciation of the mercury coumn. Wi be 3.5 times sower than the frequency of the pressure fuctuation in the siencer at an engine speed of 6 rpm. The frequency of pressure fuctuation in the siencer is equa to. Number of cyinders x speed of the engine SOUTION Natura frequency of osciation of the iquid coumn: et the datum in fug..18 be taken as the equiibrium position of the iquid. If the dispacement of the iquid coumn from the equiibrium position is denoted by x, the change in potentia energy is given by potentia energy of raised U = + iquid coumn x Figure: U-tube manometer VTU-NPTE-NMEICT Project x potentia energy of depressed iquid coumn U = weight mercury raised x dispacement of the C. G of the segment weight of mercury depressed + x dispacement of the C. G. of the segment = [Axγ] x + [Axγ] = [Ax γ] (1) Where A is the cross-sectiona area of the mercury coumn and y is the specific weight of mercury. The change in kinetic energy is given by Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 38 of 4
39 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 T = 1 [mass of mercury] [veocity] = 1 Aγ g x () Where is the ength of the mercury coumn. By assuming harmonic motion, we can write where x is the maximum dispacement and ω n is the natura frequency. x(t) = Xcosω n t (3) By substituting Eq. [3] into Eqs. [1] and [], we obtain U = U max cos ω n t T = T max sin ω n t U max = AγX (4) (5) (6) T max = 1 Aγω n X (7) g By equating U max to T max we obtain the norma frequency: ω n = g (8) ength of mercury coumn: the frequency of pressure fuctuation in the siencer = 1 6 = 3 rpm 3 π = = 1π 6 rad/sec VTU-NPTE-NMEICT Project Thus the frequency of osciation of iquid coumn in the manometer is 1 π /3.5 = 9. rad/sec. by using Eq. [8], we obtain ω n = g = 9. = =. 43 m Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 39 of 4
40 MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Dr.MOHAMED HANEEF,PRINCIPA,GHOUSIA COEGE OF ENGINEERING,RAMANAGARAM Page 4 of 4
SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationVibrations of Structures
Vibrations of Structures Modue I: Vibrations of Strings and Bars Lesson : The Initia Vaue Probem Contents:. Introduction. Moda Expansion Theorem 3. Initia Vaue Probem: Exampes 4. Lapace Transform Method
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More informationFinite element method for structural dynamic and stability analyses
Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationVTU-NPTEL-NMEICT Project
VTU-NPTEL-NMEICT Project Progress Report The Project on Development of Remaining Three Quadrants to NPTEL Phase-I under grant in aid NMEICT, MHRD, New Delhi SME Name : Course Name: Type of the Course Module
More informationSECTION A. Question 1
SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationVTU-NPTEL-NMEICT Project
MODULE-II --- SINGLE DOF FREE S VTU-NPTEL-NMEICT Project Progress Report The Project on Development of Remaining Three Quadrants to NPTEL Phase-I under grant in aid NMEICT, MHRD, New Delhi SME Name : Course
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationRadiation Fields. Lecture 12
Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationWork and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame
Work and energy method EI EI T x-axis Exercise 1 : Beam with a coupe Determine the rotation at the right support of the construction dispayed on the right, caused by the coupe T using Castigiano s nd theorem.
More informationStrain Energy in Linear Elastic Solids
Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua
More informationLobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z
Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant
More information1 Equations of Motion 3: Equivalent System Method
8 Mechanica Vibrations Equations of Motion : Equivaent System Method In systems in which masses are joined by rigid ins, evers, or gears and in some distributed systems, various springs, dampers, and masses
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationCE601-Structura Anaysis I UNIT-IV SOPE-DEFECTION METHOD 1. What are the assumptions made in sope-defection method? (i) Between each pair of the supports the beam section is constant. (ii) The joint in
More information2.1. Cantilever The Hooke's law
.1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationCABLE SUPPORTED STRUCTURES
CABLE SUPPORTED STRUCTURES STATIC AND DYNAMIC ANALYSIS OF CABLES 3/22/2005 Prof. dr Stanko Brcic 1 Cabe Supported Structures Suspension bridges Cabe-Stayed Bridges Masts Roof structures etc 3/22/2005 Prof.
More informationLecture 8 February 18, 2010
Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some
More informationModule 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur
odue 2 naysis of Staticay ndeterminate Structures by the atri Force ethod Version 2 E T, Kharagpur esson 12 The Three-oment Equations- Version 2 E T, Kharagpur nstructiona Objectives fter reading this
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationCandidate Number. General Certificate of Education Advanced Level Examination January 2012
entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More information1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2
MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationUI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE
UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationC1B Stress Analysis Lecture 1: Minimum Energy Principles in Mechanics Prof Alexander M. Korsunsky Hilary Term (January 08)
CB Stress Anaysis ecture : Minimum Energy Principes in Mechanics Prof Aexander M. Korsunsky Hiary Term (January 8) http://users.ox.ac.uk/~engs6/4me6.htm This course introduces seected chapters in Stress
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More informationMA 201: Partial Differential Equations Lecture - 11
MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f(q ) = exp ( βv( q )) 1, which is obtained by summing over
More informationProblem Set 6: Solutions
University of Aabama Department of Physics and Astronomy PH 102 / LeCair Summer II 2010 Probem Set 6: Soutions 1. A conducting rectanguar oop of mass M, resistance R, and dimensions w by fas from rest
More informationVibrations of beams with a variable cross-section fixed on rotational rigid disks
1(13) 39 57 Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks Abstract The work is focused on the probem of vibrating beams with a variabe cross-section fixed on a rotationa
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationChapter 5. Wave equation. 5.1 Physical derivation
Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing
More information3.10 Implications of Redundancy
118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationChapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory
Chapter 4 Moving Observer Method 4.1 Overview For a compete description of traffic stream modeing, one woud reuire fow, speed, and density. Obtaining these parameters simutaneousy is a difficut task if
More informationLECTURE 10. The world of pendula
LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive
More informationElements of Kinetic Theory
Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationAPPENDIX C FLEXING OF LENGTH BARS
Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationMONTE CARLO SIMULATIONS
MONTE CARLO SIMULATIONS Current physics research 1) Theoretica 2) Experimenta 3) Computationa Monte Caro (MC) Method (1953) used to study 1) Discrete spin systems 2) Fuids 3) Poymers, membranes, soft matter
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 7
Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationProcess Capability Proposal. with Polynomial Profile
Contemporary Engineering Sciences, Vo. 11, 2018, no. 85, 4227-4236 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.88467 Process Capabiity Proposa with Poynomia Profie Roberto José Herrera
More informationPhysicsAndMathsTutor.com
. Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached
More informationSupporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers
Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,
More informationInterim Exam 1 5AIB0 Sensing, Computing, Actuating , Location AUD 11
Interim Exam 1 5AIB0 Sensing, Computing, Actuating 3-5-2015, 14.00-15.00 Location AUD 11 Name: ID: This interim exam consists of 1 question for which you can score at most 30 points. The fina grade for
More information1. Measurements and error calculus
EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More information2M2. Fourier Series Prof Bill Lionheart
M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationLecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series
Lecture Notes for Math 251: ODE and PDE. Lecture 32: 1.2 Fourier Series Shawn D. Ryan Spring 212 Last Time: We studied the heat equation and the method of Separation of Variabes. We then used Separation
More informationChapter 32 Inductance
Chapter 3 nductance 3. Sef-nduction and nductance Sef-nductance Φ BA na --> Φ The unit of the inductance is henry (H). Wb T H A A When the current in the circuit is changing, the agnetic fux is aso changing.
More informationModal analysis of a multi-blade system undergoing rotational motion
Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa
More informationThe Group Structure on a Smooth Tropical Cubic
The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,
More informationApplication of the Finite Fourier Sine Transform Method for the Flexural-Torsional Buckling Analysis of Thin-Walled Columns
IOSR Journa of Mechanica and Civi Engineering (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 3-334X, Voume 14, Issue Ver. I (Mar. - Apr. 17), PP 51-6 www.iosrjournas.org Appication of the Finite Fourier Sine Transform
More informationAbout the Torsional Constant for thin-walled rod with open. cross-section. Duan Jin1,a, Li Yun-gui1
Internationa Forum on Energy, Environment Science and Materias (IFEESM 17) bout the Torsiona Constant for thin-waed rod with open cross-section Duan Jin1,a, Li Yun-gui1 1 China State Construction Technica
More informationTHE NUMERICAL EVALUATION OF THE LEVITATION FORCE IN A HYDROSTATIC BEARING WITH ALTERNATING POLES
THE NUMERICAL EVALUATION OF THE LEVITATION FORCE IN A HYDROSTATIC BEARING WITH ALTERNATING POLES MARIAN GRECONICI Key words: Magnetic iquid, Magnetic fied, 3D-FEM, Levitation, Force, Bearing. The magnetic
More information14-6 The Equation of Continuity
14-6 The Equation of Continuity 14-6 The Equation of Continuity Motion of rea fuids is compicated and poory understood (e.g., turbuence) We discuss motion of an idea fuid 1. Steady fow: Laminar fow, the
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationNumerical simulation of javelin best throwing angle based on biomechanical model
ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha
More informationMalaysian Journal of Civil Engineering 30(2): (2018)
Maaysian Journa of Ci Engineering 3():331-346 (18) BUBNOV-GALERKIN METHOD FOR THE ELASTIC BUCKLING OF EULER COLUMNS Ofondu I.O. 1, Ikwueze E. U. & Ike C. C. * 1 Dept. of Mechanica and Production Engineering,
More informationModel Solutions (week 4)
CIV-E16 (17) Engineering Computation and Simuation 1 Home Exercise 6.3 Mode Soutions (week 4) Construct the inear Lagrange basis functions (noda vaues as degrees of freedom) of the ine segment reference
More informationImproving the Reliability of a Series-Parallel System Using Modified Weibull Distribution
Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution
More informationOSCILLATIONS. dt x = (1) Where = k m
OSCILLATIONS Periodic Motion. Any otion, which repeats itsef at reguar interva of tie, is caed a periodic otion. Eg: 1) Rotation of earth around sun. 2) Vibrations of a sipe penduu. 3) Rotation of eectron
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationMeasurement of acceleration due to gravity (g) by a compound pendulum
Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationSolution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...
: Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................
More informationUniversity of California, Berkeley Physics 7A Spring 2009 (Yury Kolomensky) SOLUTIONS TO PRACTICE PROBLEMS FOR THE FINAL EXAM
1 University of Caifornia, Bereey Physics 7A Spring 009 (Yury Koomensy) SOLUIONS O PRACICE PROBLEMS FOR HE FINAL EXAM Maximum score: 00 points 1. (5 points) Ice in a Gass You are riding in an eevator hoding
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More information