Chapter #16: Structural Dynamics and Time Dependent Heat Transfer. Lectures #1-6 have discussed only steady systems. There has been no time dependence in any problems. We will investigate beam dynamics and show the additional steps in the numerical solution of a time dependent problem. There are many time dependent problems discussed in Chapter #16. 16.1 1 DOF spring mass system. 16., 16.4 & 16.5 bar element dynamic response. 16.7 truss and plane frame analysis. 16.8 time dependent (unsteady) heat transfer.
16.6 Beam Element Mass Matrices and Natural Frequencies. We will revisit the formation of the governing DE for the equilibrium of the beam element. Equilibrium now involves an inertial force acting on the differential beam element. The approximate displacement field v(x) will be revised to include the time dependence. We will see the generation of the consistent mass matrix through the evaluation of the weighted residuals. We will look at the assembly process for a beam dynamics problem and discuss the computational advantages of a lumped mass approximation. The solution process will be expanded to include the integration of the state vector ( 16.3). Define the phrase initial conditions.
Beam Dynamics: y ( A dx) ρ v x The differential element is now accelerating in the positive y direction in reaction to the external forces.
Consider what happens as a beam element moves (vibrates or translates in space). The profile of our element is defined by node coordinates and node rotations. The nodal values (the state vector d) is blended by the shape function matrix. For the moving beam the profile is fluctuating. In the case shown the slopes at both nodes 1 and and the vertical displacement at node 1 have changed over an interval Δt.
The approximation must be time dependent. vx ( ) = N d = [ N ]{ } 1 N N3 N 4 d1y φ1 dy φ The shape functions are dependent on space. They are set to satisfy element boundary conditions. 1 3 3 N1 = ( x 3x L+ L 3 ) L 1 N = x 3 L x L + xl L 1 3 N3 = 3 ( x + 3x L) L 1 3 N4 = 3 ( x L x L ) L 3 ( 3) T So, the node state vector becomes the time dependent quantity. vx (, t) = Nx ( ) dt ( ) vx ( ) = Nx ( ) dt ( ) vx ( ) = N( x) d () t
A force balance on the differential element gives. continuum mechanics V = EI ( ρ ) dv w( x) dx = Adx vx ( ) ( ) 3 dx 3 dvx We can evaluate our approximate inertial term. 4 vxt (, ) vxt (, ) EI + w + 4 ( ρ A) = x t A new term in our residual equations for the beam element. Note that the shape functions have not changed.
Since the shape functions have not changed, the contributions of the first two terms to the element equations will not change. We need only evaluate the residual term produced by the acceleration term in the motion equation. We do this in a consistent manner: we try to drive four weighted residuals to zero using Galerkin s choice of weighting functions. This leads to the phrase consistent mass matrix. + + N A dx = x= L 4 x= L x = L T vxt (, ) T T vxt (, ) N EI dx N 4 { w} dx ( ρ ) x x x x t = = =
Following the MWR procedure: x = L T vxt (, ) I ( ) a = N ρ A dx x t = x = L ( ) { } T ρ I = A N Nd dx a x = x = L T I ( ) a = ρ A N N dx d x = { } 1y φ1 y φ d = d d T The consistent mass matrix.
When we append this inertial term onto the results of our previous weighted residual evaluations: d m m L L L L L L m w f 1 1 1 1 1 6 1 6 156 54 13 y y f y L L d y L L w 1z 1z EI 4 6 φ1 ρ AL 4 13 3 φ1 + = + 3 w f L 1 6L 4 156 L y f y d y w d y 4L 4L z mz i φ i φ F k m Our approximate displacement field keeps the element equations linear (we do not see products of state variables). In MECH 33 you learned how to acquire the natural frequencies and mode shapes for a multi-dof undamped system.
An alternative to the consistent mass matrix is the lumped mass matrix. The lumped mass representation is a heuristic discretization of mass. The consistent mass matrix is formed by using the actual mass distribution within the residual evaluations. The lumped approximation assumes that the beam is very slender and it lumps the mass at node points. There is no inertia associated with changing node slopes. The lumped approximatio is used because it diagonalizes the mass matrix. Lumping is not a mathematical process it is a idealization that is executed on the element equations.
The fluctuations of the curvilinear profile do not accelerate any mass. Half of the beam s total mass is at node 1. Half of the beam s total mass is at node.
When we apply the lumped mass approximation: w 1 1 1 1 1 6 1 6 1 d f y y f y L L d y w m 1z m 1z 4L 6L L EI φ 1 AL φ ρ 1 + = 3 + w f L 1 6L 1 y f y d y d y w m 4L z m z i φ i φ F k m We will see that the lumped mass approximation eases the computational burden of modal analysis and time domain simulation (Lecture #8).
We know that in the undamped free vibration case, the transverse vibration at any point along the element will be harmonic. Consider free (natural) vibration. = k d + m d ( k D ω n ) ( iωnt k ω n m) De = m D e = iω t Must be a singular matrix. n d d 1y () t D1 y φ () t Φ = = e 1 1 d y () t D y φ () t Φ iω t n An assumed harmonic variation of the state variables. Magnitudes of the 4 vibration signals. The solution to the undamped free vibration problem is the same as the eigen problem for the system dynamic matrix.
( ω ) n = k m De ( ω ) n = 1 m k m D ( 1 ω ) n iω t = m k I D 1 ω n D= m k D n The ω values are the eigenvalues of the dynamic matrix. Dynamic Matrix The D vectors are the eigenvectors of the dynamic matrix.
Example problem 16.5:
Forming the consistent element equations. (1) d f 1y 1y 1 6L -1 6L d1y 156 L 54 13L (1) m 1z EI 4L 6L L φ1 ρ AL 4L 13L 3L φ 1 (1) = 3 + f y L 1 6L d y 4 156 L d y (1) m 4L φ 4L z i i φ Element #1 () d f y y 1 6L -1 6L d y 156 L 54 13L () m z EI 4L 6L L φ ρ AL 4L 13L 3L φ () = 3 + f3 y L 1 6L d3y 4 156 L d 3 y () m 4L φ 4L 3z i 3 i φ 3 Element #
Assembling and applying the homogeneous boundary conditions. (1) f1y F1 y 1 6L -1 6L d1y (1) m m 1z 1z 6L 4L 6L L φ1 (1) () fy + f y F y EI 1 6L 1 + 1 6L+ 6L 1 6L d y (1) () = = 3 m m 6L L 6L 6L 4L 4L 6L L z + mz z L + + φ () f F 1 6L 1 6L 3 y 3 y d3 y () m m 6 3z 3z L L 6L 4L φ3 d 156 L 54 13L 1y L 4L 13L 3L φ1 ρ AL 54 13L 156 + 156 L+ L 54 13L d y + 4 13L 3L L+ L 4L + 4L 13L 3L φ 54 13L 156 L d 3 y 13L 3L L 4L φ3
The application of the BC s is what specifies the system nature. The natural frequencies we calculate are now specific to the fixedfixed supports. EI 4 d y ρ AL 31 d y = 3 L 8L + φ 8 4 L φ The associated eigen problem is: EI 3 L 4 ρ AL 31 ( ω ) 8L 4 = 8L The natural frequencies (x) are given by: 4EI 5.68 EI 4EI.49 EI ω1 = ω 4 1 = ; ω = ω 4 = 13ρ AL L ρa ρal L ρa