Structural Element Stiffness, Mass, and Damping Matrices

Size: px
Start display at page:

Download "Structural Element Stiffness, Mass, and Damping Matrices"

Transcription

1 Structural Element Stiffness, Mass, and Damping Matrices CEE 541. Structural Dynamics Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall 18 1 Preliminaries This document describes the formulation of stiffness and mass matrices for structural elements such as truss bars, beams, plates, and cables?). The formulation of each element involves the determination of gradients of potential and kinetic energy functions with respect to a set of coordinates defining the displacements at the ends, or nodes, of the elements. The potential and kinetic energy of the functions are therefore written in terms of these nodal displacements i.e., generalized coordinates). To do so, the distribution of strains and velocities within the element must be written in terms of nodal coordinates as well. Both of these distributions may be derived from the distribution of internal displacements within the solid element. 1.1 Displacements Figure 1. Displacements within a solid continuum.

2 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin A component of a time-dependent displacement u i x, t), i = 1,, 3) in a solid continuum can be expressed in terms of the displacements of a set of nodal displacements, ū n t) n = 1,, N) and a corresponding set of shape functions ψ in. u i x, t) = N ψ in x 1, x, x 3 ) ū n t) 1) = Ψ i x) ūt) ) ux, t) = [Ψx)] 3 N ūt) 3) Engineering strain, axial strain ɛ ii, shear strain γ ij. Displacement gradient ɛ ii x, t) = u ix, t) x i u i,i x, t) 4) γ ij x, t) = u ix, t) x j + u jx, t) x i u j,i x, t) + u i,j x, t) 5) 6) u i x) x j = u i,j x) = N N x j ψ in x 1, x, x 3 ) ū n t) 7) ψ in,j x) ū n t) 8) Strain-displacement relations ɛ ii x, t) = γ ij x, t) = N ψ in,i x) ū n t) 9) N ψ in,j x) + ψ jn,i x)) ū n t) 1) Strain vector ɛ T x, t) = { ɛ 11 ɛ ɛ 33 γ 1 γ 3 γ 13 } 11) ɛx, t) = [ Bx) ] 6 N ūt) 1)

3 Structural Element Stiffness, Mass, and Damping Matrices 3 1. Geometric Deformation In frame elements, nonuniform transverse displacements u yx) ) induce longitudinal deformation, u xx). These longitudinal deformation are called geometric deformations. dx u yx) + u x) y dx u x) y u x)dx x Figure. Transverse deformation u yx) and geometric deformation u xx). dx Geometric deformation can be derived from the Pythagorean theorem and is quadratic in u xx) and u yx). dx u xx)dx) + u yx)dx) = dx) dx) u xx)dx) + u xx)) dx) = dx) u yx)dx) u xx)dx) u xx)) dx) = u yx)) dx) u xx) = 1 u xx)) + 1 u yx)) Since slender structural elements are much more stiff in axial deformation than in transverse deformation, and since elastic tensile strains are limited to within. for most structural materials, u x) <.) u y), the geometric portion of the deformation for structural frame elements is approximated as u xx) 1 u yx)). 13) The accuracy of the approximation 13) can be quantified in terms of the strains in a displaced bar. Figure 3 shows a bar with four displacement coordinates and Figure 4 shows the errors associated with the finite strain approximation. In this figure, 1 is the true stretched length of the truss bar, 1 = o + u 3 u 1 ) + u 4 u ), uses the finite-strain approximation shown above, = o + u 3 u 1 ) + 1/ o ))u 4 u ), and the logarithmic strain is ln 1 / o ). For a bar stretched to a level of linear strain u 3 u 1 )/ o that is significant for most metallic materials, the approximation 13) is accurate to within.1% for rotations u y <.1, at which point the total strain would be approximately.5, which is about twice the yield strain for most metals. For analyses of structures incorporating geometric nonlinearity but within the range of elastic behavior of metals, the approximation 13) for the geometric portion of the deformation is accurate to within one part per million.

4 4 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin u 4 u u 1 o u 3 Figure 3. A bar with four displacement coordinates in its original and displaced configurations. strain u 4 -u )/ o =.5 1 Pythagorean Thm.)/ o -1 approx)/ o -1 log 1 / o ) relative error, / u 3 -u 1 )/ o Figure 4. Finite-strain approximations and their associated errors. 1.3 Green-agrange Strain Following the development from the previous section, we can consider the planar displacement of a short line segment of length dl o to a line segment of length dl. In the original configuration, the line segment has a squared length of dl o = dx + dy and in the displaced configuration it has a squared length of dl = dx + u x,x dx + u x,y dy) + dy + u y,x dx + u y,y dy) = dx) + u x,x ) dx) + u x,y ) dy) + u x,x )dx) + u x,y )dx)dy) + u x,x u x,y )dx)dy) + dy) + u y,x ) dx) + u y,y ) dy) + u y,x )dx)dy) + u y,y )dy) + u y,x u y,y )dx)dy)

5 Structural Element Stiffness, Mass, and Damping Matrices 5 dl u y x,y) + u y,xx,y) dx + u x,y) dy y,y u x,y) y dl o x+dx, y+dy) u xx,y) + u x,xx,y) dx + u x,yx,y) dy x, y) u x,y) x Figure 5. Displacement of line segment dl o to line segment dl through planar displacements u x x, y) and u y x, y). The difference between the squared lengths, dl dl o, is quadratic in dx and dy and can be expressed in a matrix form. l lo = [ ] T [ dx u x,x + 1 u x,x + 1 u 1 y,x x,y + 1u y,x + 1u x,yu y,x 1 dy x,y + 1u y,x + 1u x,yu y,x u y,y + 1 u x,y + 1 u y,y = [ ] T [ dx 1 ɛxx ] [ ] xy dx 1 dy yx ɛ yy dy ] [ dx dy This matrix is called the Green-agrange strain tensor, or the agrangian strain tensor, or the Green-Saint Venant strain tensor for plane strain. This relationship may be confirmed in two simple cases: dx = dl o, dy =, u x,x = /l o, u x,y =, u y,y =, dl dl o = l o +. dx = dl o, dy =, u y,x = /l o, u x,x =, u y,y =, dl dl o =. ] 14) Generalizing this result to three dimensions, the strain-displacement equations for linear strain and geometric deformation are: ɛ ii = u i,i + 1 u i,i + 1 u j,i + 1 u k,i = u i x i + 1 ui x i ) + 1 uj x i ) + 1 ) uk 15) γ ij = u i,j + u j,i + u i,j u j,j + u j,i u i,i = u i x j + u j x i + u i x j u j x j + u j x i u i x i 16) In the following models, finite deformation effects are limited to ɛ xx u x,x + 1/)u y,x for the exclusive purpose of deriving geometric stiffness matrices. For metalic materials in which strains are typically around.1, this approximation does not compromise accuracy. x i

6 6 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 1.4 Stress-strain relationship isotropic elastic solid) σ 11 σ σ 33 τ 1 τ 3 τ 13 = Stress vector E 1 + ν)1 ν) 1.5 Potential Energy and Stiffness 1 ν ν ν ν 1 ν ν ν ν 1 ν 1 ν 1 ν 1 ν ɛ 11 ɛ ɛ 33 γ 1 γ 3 γ 13 17) σ T x, t) = { σ 11 σ σ 33 τ 1 τ 3 τ 13 } 18) σ = [ S e E, ν) ] 6 6 ɛ 19) Consider a system comprising an assemblage of linear springs, with stiffness k i, each with an individual stretch, d i. The total potential energy in the assemblage is U = 1 k i d i i If displacements of the assemblage of springs is denoted by a vector u, not necessarily equal to the stretches in each spring, then the elastic potential energy may also be written Uu) = 1 ut Ku = 1 n u i f i i=1 = 1 n n u i K ij u j i=1 j=1 where K is the stiffness matrix with respect to the coordinates u. The stiffness matrix K relates the elastic forces f i to the collocated displacements, u i. f 1 = K 11 u K 1j u j + + K 1N u N f i = K i1 u K ij u j + + K in u N f N = K N1 u K Nj u j + + K NN u N A point force f i acting on an elastic body is the gradient of the elastic potential energy U with respect to the collocated displacement u i f i = u i U The i, j term of the stiffness matrix may therefore be found from the potential energy function Uu), K ij = Uu) ) u i u j

7 Structural Element Stiffness, Mass, and Damping Matrices Strain Energy and Stiffness in inear Elastic Continua Elastic element stiffness matrix 1.7 Kinetic Energy and Mass Uū) = 1 σx, t) T ɛx, t) dω 1) Ω = 1 ɛx, t) T S e E, ν) ɛx, t) dω Ω = 1 ūt) T Bx) T S e E, ν) Bx)ūt) dω Ω = 1 [ ūt)t Bx) T S e E, ν) Bx) ] dω ūt) ) N N Ω fe = U ū = K e ū K e = The impulse-momentum relationship states that f dt = δm u) Ω [ Bx) T S e E, ν) Bx) ] dω 3) N N f = d m u) dt f = d ) 1 dt u m u f = d ) dt u T, where T is the kinetic energy, and is assumed to be independent of ut). Consider a system comprising an assemblage of point masses, m i, each with an individual velocity, v i. The total kinetic energy in the assemblage is T = 1 m i vi i If displacements of the assemblage of masses are defined by a generalized coordinate vector u, not necessarily equal to the velocity coordinates, above, then the kinetic energy may also be written T u) = 1 ut M u = 1 n i=1 u i n j=1 M ij u j

8 8 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin where M is the constant mass matrix with respect to the generalized coordinates u. The mass matrix M relates the inertial forces f i to the collocated accelerations, ü i. f 1 = M 11 ü M 1j ü j + + M 1N ü N f i = M i1 ü M ij ü j + + M in ü N f N = M N1 ü M Nj ü j + + M NN ü N The i, j term of the constant mass matrix may therefore be found from the kinetic energy function T, M ij = T u) = T u) 4) ü i t u j u i u j 1.8 Inertial Energy and Mass in Deforming Continua Consistent mass matrix T ū) = 1 ρ ux, t) dω 5) Ω = 1 ρ ux, t) T ux, t) dω Ω = 1 ρ ūt) T Ψx) T Ψx) ūt) dω Ω = 1 ūt) T ρ [ Ψx) T Ψx) ] dω ūt) 6) N N Ω T ū Ω = ρ [ Ψx) T Ψx) ] dω ūt) 7) N N fi = d ) T = ρ [ Ψx) dt ū T Ψx) ] dω ūt) 8) Ω N N M = ρ [ Ψx) T Ψx) ] dω 9) N N Ω

9 Structural Element Stiffness, Mass, and Damping Matrices 9 Bar Element Matrices D prismatic homogeneous isotropic truss bar. Uniform uni-axial stress σ T = {σ xx,,,,, } T Corresponding strain ɛ T = σ xx /E) {1, ν, ν,,, } T. Incremental strain energy du = 1 σt ɛ dω = 1σ xxɛ xx dω = 1 Eɛ xx dω.1 Bar Displacements Figure 6. Truss bar element coordinates and displacements. u x x, t) = 1 x ) x ū 1 t) + ū 3 t) 3) ) = ψ x1 x) ū 1 t) + ψ x3 x) ū 3 t) 31) u y x, t) = 1 x ) x ū t) + ū 4 t) 3) ) = ψ y x) ū t) + ψ y4 x) ū 4 t) 33) Ψx) = [ 1 x x 1 x x ] 34) [ ux x, t) u y x, t) ] = Ψx) ūt) 35)

10 1 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin. Bar Strain Energy and Elastic Stiffness Matrix Strain-displacement relation ɛ xx = u x x + 1 ) uy 36) x = ψ x1,x ū 1 + ψ x3,x ū {ψ y,xū + ψ y4,x ū 4 } 37) = = 1 ) 1 ū 1 + ū 3 + ) 1 { 1 [ 1 1 ] ū + 1 { 1 = B ū + 1 ) } 1 ū + ū 4 38) ) ) } 1 ū + ū 4 39) ) { 1 ) } 1 ū + ū 4 4) ) B = [ 1 1 ]. 41) Strain energy and elastic stiffness U = 1 K e = Ω x= = EA = EA ɛ xx E ɛ xx dω 4) [ B T E B ] A dx 43) x= 1/ 1/ 1/ 1/ dx 44) 45)

11 Structural Element Stiffness, Mass, and Damping Matrices 11.3 Bar Kinetic Energy and Mass Matrix T = 1 ū T Ω ρ [ Ψx) T Ψx) ] N N dω ūt) 46) M = = ρa x= x= = 1 6 ρa ρ [ Ψx) T Ψx) ] A dx 47) 1 x ) 1 x ) x ) 1 x ) 1 x ) x ) x )1 x ) x ) x )1 x ) x ) Bar Stiffness Matrix with Geometric Strain Effects dx 48) 49) Substitute to obtain So, U ū = EA U = 1 = EA = EA ɛ xx E ɛ xx A dx 5) U = EA u x x + 1 ux x u x x u y x uy x ) + u x x ) dx 51) ) uy + 1 x 4 ) 4 uy dx 5) x = 1 ū1 + 1 ū3 53) = 1 ū + 1 ū4 54) ū 3 ū 1 ) + 1 ū 3 ū 1 )ū 4 ū ) ) ū 1 ū ū 3 ū 4 + EAū 3 ū 1 ) ū 1 ū ū 3 ū 4 55) 56) = K e ū + N K g ū 57)

12 1 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 3 Bernoulli-Euler Beam Element Matrices D prismatic homogeneous isotropic beam element, neglecting shear deformation and rotatory inertia. 3.1 Bernoulli-Euler Beam Coordinates and Internal Displacements Consider the geometry of a deformed beam. The functions u x x) and u y x) describe the translation of points along the neutral axis of the beam as a function of the location along the un-stretched neutral axis. Figure 7. Beam element coordinates and displacements. We will describe the deformation of the beam as a function of the end displacements ū 1, ū, ū 4, ū 5 ) and the end rotations ū 3, ū 6 ). In a dynamic context, these end displacements will change with time. u x x, t) = u y x, t) = ψ xn x) ū n t) ψ yn x) ū n t) The functions ψ xn x) and ψ yn x) satisfy the boundary conditions at the end of the beam and the differential equation describing bending of a Bernoulli-Euler beam loaded statically at the nodal coordinates. In such beams the effects of shear deformation and rotatory inertia are neglected. For extension of the neutral axis, ψ x1 x) = 1 x ψ x4 x) = x

13 Structural Element Stiffness, Mass, and Damping Matrices 13 and ψ x = ψ x3 = ψ x5 = ψ x6 = along the neutral axis. For bending of the neutral axis, ) x x 3 ψ y x) = ) ) ) x 3 ψ y3 x) = x ψ y5 x) = 3 ψ y6 x) = x x + ) ) x 3 ) ) x ) x 3 + ) and ψ y1 = ψ y4 =. Ψx) = [ 1 x x 1 3 x ) + x ) 3 x x ) + x ) ) 3 3 x ) x ) 3 x ) + x ) 3 ) ] 58) [ ux x, t) u y x, t) ] = Ψx) ūt) 59) These expressions are analytical solutions for the displacements of Bernoulli-Euler beams loaded only with concentrated point loads and concentrated point moments at their ends. Internal bending moments are linear within beams loaded only at their ends, and the beam displacements may be expressed with cubic polynomials. 3. Bernoulli-Euler Beam Strain Energy and Elastic Stiffness Matrix In extension, the elastic potential energy in a beam is the strain energy related to the uniform extensional strain, ɛ xx. If the strain is small, then the extensional strain within the cross section is equal to an extension of the neutral axis, u x / x), plus the bending strain, u y / x )y. ɛ xx = u x x u y x y = x ψ xnx) ū n = = ψ xnx) ū n B n x, y) ū n x ψ ynx) y ū n 6) ψ ynx) y ū n = Bx, y) ū 61) where Bx, y) = [ 1, 6y 1xy, 4y 3 6xy, 1 3, 6y + 1xy, y 3 6xy ]. 6)

14 14 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin The elastic stiffness matrix can be found directly from the strain energy of axial strains ɛ xx. U = 1 K e = Ω x= ɛ xx E ɛ xx dω 63) A [ Bx, y) T E Bx, y) ] da dx. 64) Note that this integral involves terms such as A y da and A yda in which the origin of the coordinate axis is placed at the centroid of the section. The integral A y da is the bending moment of inertia for the cross section, I, and the integral A yda is zero. It is also important to recognize that the elastic strain energy may be evaluated separately for extension effects and bending effects. For extension, the elastic strain energy is U = 1 EA ɛ xx ) dx x= = 1 ) EA ψ xnx) ū n x= dx and the ij stiffness coefficient for indices 1 and 4) is K ij = 1 ) EA ψ ū i ū j xnx) ū n x= = x= EA ψ xix) ψ xjx) dx. 65) In bending, the elastic potential energy in a Bernoulli-Euler beam is the strain energy related to the curvature, κ z. κ z = u y x = x ψ ynx) ū n = The elastic strain energy for pure bending is dx ψ ynx) ū n U = 1 EI κ z ) dx x= = 1 EI ψ ynx) ū n) dx x= and the ij stiffness coefficient for indices,3,5 and 6) is K ij = 1 EI ψ ū i ū j ynx) ū n) dx x= = x= EI ψ yix) ψ yjx) dx. 66)

15 Structural Element Stiffness, Mass, and Damping Matrices Bernoulli-Euler Beam Kinetic Energy and Mass Matrix The kinetic energy of a particle within a beam is half the mass of the particle, ρadx, times its velocity, u, squared. For velocities along the direction of the neutral axis, u x x) = ψ xn x) ū n, The kinetic energy function and the mass matrix may be by substituting equation 58) into equations 6) and 9). T = 1 ū T Ω ρ [ Ψx) T Ψx) ] N N dω ūt) 67) M = x= ρ [ Ψx) T Ψx) ] A dx 68) It is important to recognize that kinetic energy and mass associated with extensional velocities may be determined separately from those associated with transverse velocities. The kinetic energy for extension of the neutral axis is T = 1 ρa u x ) dx x= = 1 ) ρa ψ xn x) ū n dx x= and the ij mass coefficient for indices 1 and 4) is M ij = 1 ) ρa ψ xn x) ū n dx ū i ū j x= = x= For velocities transverse to the neutral axis, ρa ψ xi x) ψ xj x) dx. 69) u y x) = ψ yn x) ū n, the kinetic energy for velocity across the neutral axis is T = 1 ρa u y ) dx x= = 1 ) ρa ψ yn x) ū n dx x= and the ij mass coefficient for indices,3,5 and 6) is M ij = 1 ) ρa ψ yn x) ū n dx ū i ū j x= = x= ρa ψ yi x) ψ yj x) dx. 7)

16 16 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 3.4 Bernoulli-Euler Stiffness Matrix with Geometric Strain Effects The axial strain in a Bernoulli-Euler beam including the geometric strain is ɛ xx = u x x u y x y + 1 The potential energy with geometric strain effects is U = 1 = 1 = 1 x= A E E A A ) uy 71) x ɛ xx E ɛ xx da dx 7) ux x u y x y + 1 ) ) uy dx 73) x u x,x u x,x u y,xx y + u x,x u y,x + u y,xxy u y,xx u y,xy u4 y,x Note that A yda = and A y da = I and neglect u 4 y,x so that ) dadx 74) U = 1 EA u x,x ) 1 dx + EI ) u y,xx dx + EA u x,x u y,x) dx. 75) Substitute u y,x = u y,xx = u x,x = ψ ynx) ū n 76) ψ ynx) ū n 77) ψ xn ū n = N 78) EA and differentiate with respect to ū i and ū j to obtain, K ij = EA ψ xiψ xj dx + EI ψ yix)ψ yjx) dx + N ψ yix)ψ yjx) dx 79) so that, K = K e + N K g 8)

17 Structural Element Stiffness, Mass, and Damping Matrices Bernoulli-Euler Beam Element Stiffness and Mass Matrices For prismatic homogeneous isotropic beams, substituting the expressions for the functions ψ xn and ψ yn into equations 65) - 7), or substituting equation 6) into equation 64) and 58) to equation 68) results in element stiffness matrices K e, M, and K g. K e = M K g = ρa 4 = N EA EA 1EI 3 sym 6EI 1EI 3 4EI 6EI 6EI EI EA 1EI 3 6EI EI sym sym ) 8) 83)

18 18 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 4 Timoshenko Beam Element Matrices D prismatic homogeneous isotropic beam element, including shear deformation and rotatory inertia Consider again the geometry of a deformed beam. When shear deformations are included sections that are originally perpendicular to the neutral axis may not be perpendicular to the neutral axis after deformation. The functions u x x) and u y x) describe the translation of Figure 8. Deformation of beam element including shear deformation. points along the neutral axis of the beam as a function of the location along the un-stretched neutral axis. If the beam is not slender length/depth < 5), then shear strains will contribute significantly to the strain energy within the beam. The deformed shape of slender beams is different from the deformed shape of stocky beams. The beam carries a bending moment Mx) related to axial strain ɛ xx and a shear force, S related to shear strain γ xy. The potential energy has a bending strain component and a shear strain component. U = 1 σ T ɛ dω Ω = 1 σ xx ɛ xx dω + 1 τ xy γ xy dω Ω Ω = 1 Mx)y Mx)y da dx + 1 SQy) SQy) da dx A I EI A Iby) GIby) = 1 = 1 Mx) y da dx + 1 EI A Mx) dx + 1 EI S GA/α) S GI A Qy) by) da dx dx 84) where the shear area coefficient α reduces the cross section area to account for the non-uniform distribution of shear stresses in the cross section, α = A Qy) I A by) da. For solid rectangular sections α = 6/5 and for solid circular sections α = 1/9 [, 3, 4, 5, 8].

19 Structural Element Stiffness, Mass, and Damping Matrices Timoshenko Beam Coordinates and Internal Displacements including shear deformation effects) The transverse deformation of a beam with shear and bending strains may be separated into a portion related to shear deformation and a portion related to bending deformation, u y x, t) = u b)y x) + u s)y x) 85) where EIu b)yx) = Mx) = M 1 1 x ) ) x + M 86) GA/α)u s)yx) = Sx) = 1 M 1 + M ) 87) It can be shown that the following shape functions satisfy the Timoshenko beam equations equations 85), 86) and 87)) for transverse displacements. ψ y x) = ψ y3 x) = ψ y5 x) = ψ y6 x) = [ ) 1 x x Φ ) x ) ] Φ [ ) x x x 3 ) ) ] 1 + Φ 1 x x + + ) Φ [ ) 1 x ) ] x 3 x Φ Φ [ ) x ) x 3 ) ) ] 1 x x Φ Φ The term Φ gives the relative importance of the shear deformations to the bending deformations, Φ = 1EI ) r GA/α) = 4α1 + ν), 88) where r is the radius of gyration of the cross section, r = I/A, ν is Poisson s ratio. Shear deformation effects are significant for beams which have a length-to-depth ratio less than 5. To neglect shear deformation, set Φ =. These displacement functions are exact for frame elements with constant shear forces S and linearly varying bending moment distributions, Mx), in which the strain energy has both a shear stress component and a normal stress component, U = 1 ) EI ψ b)ynx)ū n dx + 1 ) GA/α) ψ s)ynx)ū n dx 89) where the bending and shear components of the shape functions, ψ b)yn x) and ψ s)yn x) are:

20 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin [ ) 1 x ] x 3 ψ b)y x) = Φ ) [ Φ ψ s)y x) = 1 x ] 1 + Φ [ ) x x x 3 ψ b)y3 x) = 1 + Φ ) x ) ) ] x Φ ψ s)y3 x) = Φ [ ] 1 x 1 + Φ [ ) 1 x ] x 3 ψ b)y5 x) = Φ ) [ Φ x ψ s)y5 x) = ] ψ b)y6 x) = 1 + Φ [ 1 + Φ ψ s)y6 x) = 1 + Φ ) x x + ] x Φ 4. Timoshenko Beam Element Stiffness Matrices [ 1 ) ) ) ] x Φ The geometric stiffness matrix for a Timoshenko beam element may be derived as was done with the Bernoulli-Euler beam element from the potential energy of linear and geometric strain, K ij = EA + EI + GA/α) + N ψ xix)ψ xjx) dx ψ b)yix)ψ b)yjx) dx ψ s)yix)ψ s)yjx) dx ψ yix)ψ yjx) dx 9) where the displacement shape functions ψx) are provided in section 4.1.

21 Structural Element Stiffness, Mass, and Damping Matrices Timoshenko Beam Element Stiffness and Mass Matrices, including shear deformation effects but not rotatory inertia) For prismatic homogeneous isotropic beams, substituting the previous expressions for the functions ψ xn x) and ψ b)yn x), and ψ s)yn x) into equation 9) and 68), results in the Timoshenko element elastic stiffness matrices K e, mass matrix M, and geometric stiffness matrix K g EA EA M = ρa 84 K e = 1 EI 6 EI 1 EI 6 EI 1+Φ 3 1+Φ 1+Φ 3 1+Φ sym 4+Φ EI 6 EI 1+Φ 1+Φ Φ EI 1+Φ EA 1 EI 6 EI 1+Φ 3 1+Φ 4+Φ EI 1+Φ Φ + 8Φ Φ + 35Φ ) Φ + 175Φ Φ + 35Φ ) Φ + 7Φ ) Φ + 35Φ ) Φ + 7Φ ) 9) 8 sym Φ + 8Φ Φ + 35Φ ) K g = N /15+ Φ/6+ Φ /1 1+Φ) sym 6/5+Φ+Φ 1+Φ) 6/5+Φ+Φ 1+Φ) /1 1+Φ) 6/5 Φ Φ 1+Φ) /1 1+Φ) Φ + 7Φ ) /1 /3 Φ/6 Φ /1 1+Φ) 1+Φ) /1 1+Φ) /15+ Φ/6+ Φ /1 1+Φ) 91) 93)

22 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 4.4 Timoshenko Beam Element Mass Matrix including rotatory inertia but not shear deformation effects) Consider again the geometry of a deformed beam with linearly-varying axial beam displacements outside of the neutral axis. The functions u x x, y) and u y x, y) now describe the Figure 9. Deformation of beam element showing axial-direction displacements u x x, y, t) outside the neutral axis. translation of points anywhere within the beam, as a function of the location within the beam. We will again describe these displacements in terms of a set of shape functions, ψ xn x, y) and ψ yn x), and the end displacements ū 1,, ū 6. u x x, y, t) = u y x, t) = ψ xn x, y) ū n t) ψ yn x) ū n t) The shape functions for transverse displacements ψ yn x) are the same as the shape functions ψ yn x) used previously. The shape functions for axial displacements along the neutral axis, ψ x1 x, y) and ψ x4 x, y) are also the same as the shape functions ψ x1 x) and ψ x4 x) used previously. To account for axial displacements outside of the neutral axis, four new shape functions are derived from the assumption that plane sections remain plane, u x x, y) = u b)y x)y. ψ x x, y) = ψ y y = x x y 6 )) ψ x3 x, y) = ψ y3 y = x ) ) x 3 y ψ x5 x, y) = ψ y5 y = 6 x ) ) x y + ψ x6 x, y) = ψ y6 y = x ) ) x 3 y

23 Structural Element Stiffness, Mass, and Damping Matrices 3 Because ψ yn, ψ x1 and ψ x4 are unchanged, the stiffness matrix is also unchanged. The kinetic energy of the beam, including axial and transverse effects is now, T = 1 h/ ρby) ψ xn x, y) ū n) dy dx + 1 ) ρa ψ yn x) ū n x= y= h/ x= dx and the mass matrix coefficients are found from M ij = ū i ū j T ū) Evaluating equation 4) using the new shape functions ψ x, ψ x3, ψ x5, and ψ x6, results in a mass matrix incorporating rotatory inertia. M = ρa r r 9 6 r r r 4 1 sym r r 1 r r 94) Beam element mass matrices including the effects of shear deformation on rotatory inertia are more complicated. Refer to p 95 of Theory of Matrix Structural Analysis, by J.S. Przemieniecki Dover Pub., 1985).

24 4 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 5 Coordinate Transformations for Bars and Beams 5.1 Beam Element Stiffness Matrix in ocal Coordinates, K EI θ 3EI 3 3EI θ 6EI θ 1EI 3 3EI 3EI θ 4EI θ 6EI EI θ EI θ 3EI 3 3EI θ 6EI θ 1EI 3 6EI N 1 V 1 M 1 N V M = EA EA 1EI 3 sym 6EI 1EI 3 4EI 6EI 6EI EI EA 1EI 3 6EI 4EI ū 1 ū ū 3 ū 4 ū 5 ū 6 f = K ū

25 Structural Element Stiffness, Mass, and Damping Matrices 5 5. Beam Element Stiffness Matrix in Global Coordinates, K Geometric relationship between ū and u: ū = T u ū 1 = u 1 cos θ + u sin θ ū = u 1 sin θ + u cos θ ū 3 = u 3 where T = c s s c 1 c s s c 1 c = cos θ = x x 1 s = sin θ = y y 1 The stiffness matrix in global coordinates is K = T T K T K = EA c + 1EI EA cs 3 s 1EI 3 EA s EA c cs 6EI s 1EI + 1EI 3 c 6EI sym c 4EI EA cs 3 s + 1EI cs 6EI 3 EA s cs EA s + 1EI cs 1EI 3 c 6EI 3 6EI s EA c 6EI c EA c EI cs + 1EI s 1EI 6EI 3 cs 3 s EA s + 1EI 3 c 6EI c 4EI f = K u

26 6 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 5.3 Beam Element Consistent Mass Matrix in ocal Coordinates, M N 1 V 1 M 1 N V M = ρa sym ū 1 ū ū 3 ū 4 ū 5 ū 6 f = M ū

27 Structural Element Stiffness, Mass, and Damping Matrices Beam Element Consistent Mass Matrix in Global Coordinates, M Geometric relationship between ū and u: ū = T u ū 1 = u 1 cos θ + u sin θ ū = u 1 sin θ + u cos θ ū 3 = u 3 where T = c s s c 1 c s s c 1 c = cos θ = x x 1 s = sin θ = y y 1 The consistent mass matrix in global coordinates is M = T T M T M = ρa 4 14c 16cs s 7c 16cs 13s +15s +54s 14s c 16cs 7s 13c +156c +54c 4 13s 13c 3 14c 16cs s sym +156s 14s c +156c 4 f = M u

28 8 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 6 D Plane-Stress and Plane-Strain Rectangular Element Matrices D, isotropic, homogeneous element, with uniform thickness h. Approximate element stiffness and mass matrices based on assumed distribution of internal displacements. 6.1 D Rectangular Element Coordinates and Internal Displacements Consider the geometry of a rectangle with edges aligned with a Cartesian coordinate system. x a, y b) The functions u x x, y, t) and u y x, y, t) describe the in-plane displacements as a function of the location within the element. Figure 1. D rectangular element coordinates and displacements. Internal displacements are assumed to vary linearly within the element. x u x x, y, t) = c 1 a + c x y a b + c y 3 b + c 4 x u y x, y, t) = c 5 a + c x y 6 a b + c y 7 b + c 8

29 Structural Element Stiffness, Mass, and Damping Matrices 9 The eight coefficients c 1,, c 8 coordinates at the corners. may be found uniquely from matching the displacement resulting in internal displacements u x a, b) = ū 1, u y a, b) = ū u x, b) = ū 3, u y, b) = ū 4 u x, ) = ū 5, u y, ) = ū 6 u x a, ) = ū 7, u y a, ) = ū 8 u x x, y, t) = ˆxŷ ū 1 t) + 1 ˆx)ŷ ū 3 t) + 1 ˆx)1 ŷ) ū 5 t) + ˆx1 ŷ) ū 7 t) 95) u y x, y, t) = ˆxŷ ū t) + 1 ˆx)ŷ ū 4 t) + 1 ˆx)1 ŷ) ū 6 t) + ˆx1 ŷ) ū 8 t) 96) where ˆx = x/a ˆx 1) and ŷ = y/b ŷ 1) so that [ ˆxŷ 1 ˆx)ŷ 1 ˆx)1 ŷ) ˆx1 ŷ) Ψˆx, ŷ) = ˆxŷ 1 ˆx)ŷ 1 ˆx)1 ŷ) ˆx1 ŷ) 97) and [ ] ux x, y, t) = Ψˆx, ŷ) ūt) 98) u y x, y, t) ] Strain-displacement relations so that or ɛ xx ɛ yy γ xy = ɛ xx = u x x = 1 u x a ˆx ɛ yy = u y y = 1 u y b ŷ γ xy = u x y + u y x = 1 b u x ŷ + 1 u y a ˆx ŷ/a ŷ/a 1 ŷ)/a 1 ŷ)/a ˆx/b 1 ˆx)/b 1 ˆx)/b ˆx/b ˆx/b ŷ/a 1 ˆx)/b ŷ/a 1 ˆx)/b 1 ŷ)/a ˆx/b 1 ŷ)/a ɛx, y, t) = Bx, y) ūt) ū 1 ū ū 3 ū 4 ū 5 ū 6 ū 7 ū 8

30 3 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 6. Stress-Strain relationships 6..1 Plane-Stress In-plane behavior of thin plates, σ zz = τ xz = τ yz = For plane-stress elasticity, the stress-strain relationship simplifies to σ xx σ yy = E 1 ν ɛ xx 1 ν τ ɛ yy xy γ xy or ν ν) 99) σ = S pσ ɛ 1) 6.. Plane-Strain In-plane behavior of continua, ɛ zz = γ xz = γ yz = For plane-strain elasticity, the stress-strain relationship simplifies to σ xx 1 ν ν σ yy = E 1 + ν)1 ν) ν 1 ν 1 τ xy ν or ɛ xx ɛ yy γ xy 11) σ = S pɛ ɛ 1) 6.3 D Rectangular Element Strain Energy and Elastic Stiffness Matrix V = 1 σx, y, t) T ɛx, y, t) h dx dy 13) A = 1 [ ūt)t Bx, y) T S e E, ν) Bx, y) ] h dx dy ūt) 14) 8 8 A Elastic element stiffness matrix [ K e = Bx, y) T S e E, ν) Bx, y) ] h dx dy 15) 8 8 A 6.4 D Rectangular Element Kinetic Energy and Mass Matrix T ū) = 1 ρ ux, y, t) h dx dy 16) A = 1 ūt) T ρ [ Ψx, y) T Ψx, y) ] h dx dy ūt) 17) 8 8 A Consistent mass matrix M = A ρ [ Ψx, y) T Ψx, y) ] 8 8 h dx dy 18)

31 Structural Element Stiffness, Mass, and Damping Matrices D Rectangular Plane-Stress and Plane-Strain Element Stiffness and Mass Matrices Plane-Stress stiffness matrix K e = Eh 11 ν ) 4c + k A k B 4c + k A / k C c k A / k B c k A k C k B 4/c + k D k C /c k D k B /c k D / k C 4/c + k D / 4c + k A / k C 4c + k A k B c k A k C c k A / k B k C /c k D k B 4/c + k D k C 4/c + k D / k B /c k D / c k A / k B c k A k C 4c + k A k B 4c + k A / k C k B /c k D / k C 4/c + k D / k B 4/c + k D k C /c k D c k A k C c k A / k B 4c + k A / k C 4c + k A k B k C 4/c + k D / k B /c k D / k C /c k D k B 4/c + k D where c = b/a and k A = /c)1 ν) k B = 3/)1 + ν) k C = 3/)1 3ν) k D = c)1 ν) 6.5. Plane-Strain stiffness matrix K e = Eh 11+ν)1 ν) k A + k B 3/ k A + k B / 6ν 3/ k A / k B / 3/ k A / k B 3/ 6ν 3/ k C + k D 3/ 6ν k C / k D 3/ k C / k D / 6ν 3/ k C + k D / k A + k B / 3/ 6ν k A + k B 3/ k A / k B 6ν 3/ k A / k B / 3/ 6ν 3/ k C / k D 3/ k C + k D 3/ 6ν k C + k D / 3/ k C / k D / k A / k B / 3/ k A / k B 3/ 6ν k A + k B 3/ k A + k B / 6ν 3/ 3/ k C / k D / 6ν 3/ k C + k D / 3/ k C + k D 3/ 6ν k C / k D k A / k B 6ν 3/ k A / k B / 3/ k A + k B / 3/ 6ν k A + k B 3/ 3/ 6ν k C + k D / 3/ k C / k D / 6ν 3/ k C / k D 3/ k C + k D where c = b/a and k A = 4c)1 ν) k B = /c)1 ν) k C = 4/c)1 ν) k D = c)1 ν)

32 3 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin Mass matrix The element mass matrix for the plane-stress and plane-strain elements is the same M = ρabh ) Note, again, that these element stiffness matrices are approximations based on an assumed distribution of internal displacements. 7 Element damping matrices Damping in vibrating structures can arise from diverse linear and nonlinear phenomena. If the structure is in a fluid liquid or gas), the motion of the structure is resisted by the fluid viscosity. At low speeds low Reynolds numbers), this damping effect can be taken to be linear in the velocity, and the damping forces are proportional to the total rate of displacement not the rate of deformation). If the fluid is flowing past the structure at high flow rates high Reynolds numbers), the motion of the structure can interact with the flowing medium. This interaction affects the dynamics natural frequencies and damping ratios) of the coupled structure-fluid system. Potentially, at certain flow speeds, the motion of the structure can increase the transfer of energy from the flow into the structure, giving rise to an aero-elastic instabililty. Damping can also arise within structural systems from friction forces internal to the structure the micro-slip within joints and connections) inherent material viscoelasticity, and inelastic material behavior. In many structural systems, a type of damping in which damping stresses are proportional to strain and in-phase with strain-rate are assumed. Such so-called complex-stiffness damping or structural damping is commonly used to model the damping in soils. Fundamentally, this kind of damping is neither elastic nor viscous. The forcedisplacement behavior does not follow the same path in loading and unloading, behavior but instead follows a butterfly shaped path. Nevertheless, this type of damping is commonly linearized as linear viscous damping, in which forces are proportional to the rate of deformation. In materials in which stress depends on strain and strain rate, a Voigt viscoelasticity model may be assumed, in which stress is proportional to both strain ɛ and strain-rate ɛ, σ = [ S e E, ν) ] ɛ + [ S v η) ] ɛ

33 Structural Element Stiffness, Mass, and Damping Matrices 33 The internal virtual work of real viscous stresses S v ɛ moving through virtual strains δɛ is δw ū) = σx, t) T δɛx, t) dω 11) Ω = ɛx, t) T S v η) δɛx, t) dω Ω = ūt) T Bx) T S v η) Bx)δūt) dω Ω [ = ūt) T Bx) T S v η) Bx) ] dω δūt) 111) N N Ω Given a material viscous damping matrix, S v, a structural element damping matrix can be determined for any type of structural element, through the integral in equation 111), as has been done for stiffness and mass element matrices earlier in this document. In doing so, it may be assumed that the internal element displacements u i x, t) and the matrices [Ψ] and [B]) are unaffected by the presence of damping, though this is not strictly true. Further, the parameters in S v η) are often dependent of the frequency of the strain and the strain amplitude. Damping behavior that is amplitude-dependent is outside the domain of linear analysis. 7.1 Rayleigh damping matrices for structural systems In an assembled model for a structural system, a damping matrix that is proportional to system s mass and stiffness matrices is called a Rayleigh damping matrix. R T C s R = ζ 1 ω n1... ζ N ω nn C s = αm s + βk s 1 = α β ω n 1... ω n N 11) where ω nj is an eigen-value squared natural frequency) and the columns of R are massnormalized eigen-vectors modal vectors) of the generalized eigen-problem [K s ω n j M s ] r j =. 113) From equations 11) it can be seen that the damping ratios satisfy ζ j = α 1 ω nj + β ω nj and the Rayleigh damping coefficients α and β) can be determined so that the damping ratios ζ j have desired values at two frequencies. The damping ratios modeled by Rayleigh damping can get very large for low and high frequencies. Rayleigh damping grows to as ω and increases linearly with ω for large values of ω. Note that the Rayleigh damping matrix has the same banded form as the mass and stiffness matrices. In other words, with Rayleigh damping, internal damping forces are applied only between coordinates that are connected by structural elements.

34 34 CEE 541. Structural Dynamics Duke University Fall 18 H.P. Gavin 7. Caughey damping matrices for structural systems The Caughey damping matrix is a generalization of the Rayleigh damping matrix. Caughey damping matrices can involve more than two parameters and can therefore be used to provide a desired amount of damping over a range of frequencies. The Caughey damping matrix for an assembled model for a structural system is j=n C s = M s α j M 1 s K s ) j j=n 1 where the index range limits n 1 and n can be positive or negative, as long as n 1 < n. As with the Rayleigh damping matrix, the Caughey damping matrix may also be diagonalized by the real eigen-vector matrix R. The coefficients α j are related to the damping ratios, ζ k, by ζ k = 1 1 ω k j=n α j ω j k j=n 1 The coefficients α j may be selected so that a set of specified damping ratios ζ k are obtained at a corresponding set of frequencies ω k. If n 1 = and n = 1, then the Caughey damping matrix is the same as the Rayleigh damping matrix. For other values of n 1 and n the Caughey damping matrix loses the banded structure of the Rayleigh damping matrix, implying the presence of damping forces between coordinates that are not connected by structural elements. Structural systems with classical damping have real-valued modes r j that depend only on the system s mass and stiffness matrices equation 113)), and can be analyzed as a system of uncoupled second-order ordinary differential equations. The responses of the system coordinates can be approximated via a modal expansion of a select subset of modes. The convenience of the application of modal-superpostion to the transient response analysis of structures is the primary motivation 7.3 Rayleigh damping matrices for structural elements An element Rayleigh damping matrix may be easily computed from the element s mass and stiffness matrix C = αm+βk and assembled into a damping matrix for the structural system C s. The element damping is presumed to increases linearly with the mass and the stiffness of the element; larger elements will have greater mass, stiffness, and damping. System damping matrices assembled from such element damping matrices will have the same banding as the mass and stiffness matrices; internal damping forces will occur only between coordinates connected by a structural element. However, such an assembled damping matrix will not be diagonizeable by the real eigenvectors of the structural system mass matrix M s and stiffness matrix K s.

35 Structural Element Stiffness, Mass, and Damping Matrices inear viscous Damping elements Some structures incorporate components designed to provide supplemental damping. These supplemental damping components can dissipate energy through viscosity, friction, or inelastic deformation. In a linear viscous damping element a dash-pot), damping forces are linear in the velocity across the nodes of the element and the forces act along a line between the two nodes of the element. The element node damping forces f d are related to the element node velocities v d through the damping coefficient c d [ fd1 f d ] = [ c d c d c d c d ] [ vd1 The damping matrix for a linear viscous damper connecting a node at x 1, y 1 ) to a node at x, y ) is found from the element coordinate transformation, v d ] C 6 6 = [ c s c s ] T [ c d c d c d c d ] [ c s c s ] where c = x x 1 )/ and s = y y 1 )/. Structural systems with supplemental damping components generally have non-classical system damping matrices. References [1] Clough, Ray W., and Penzien, Joseph, Dynamics of Structures, nd ed. revised), Computers and Structures, 3. [] Cowper, G.R., Shear Coefficient in Timoshenko Beam Theory, J. Appl. Mech., 33)1966): [3] Dong, S.B, Alpdogan, C., and Taciroglu, E. Much ado about shear correction factors in Timoshenko beam theory, Int. J. Solids & Structures, 471): [4] Gruttman, F., and Wagner, W., Shear correction factors in Timoshenko s beam theory for arbitrary shaped cross-sections, Comp. Mech. 71): [5] Kaneko, T., An experimental study of the Timoshenko s shear coefficient for flexurally vibrating beams, J. Phys. D: Appl. Phys ): ; [6] Paz, Mario, Structural Dynamics Theory and Computation, Chapman & Hall,. [7] Przemieniecki, J.S., Theory of Matrix Structural Analysis, Dover, [8] Rosinger, H.E., and Ritchie, I.G., On Timoshenko s correction for shear in vibrating isotropic beams, J. Phys. D: Appl. Phys., ):

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations

More information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly .3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original

More information

Geometric Stiffness Effects in 2D and 3D Frames

Geometric Stiffness Effects in 2D and 3D Frames Geometric Stiffness Effects in D and 3D Frames CEE 41. Matrix Structural Analsis Department of Civil and Environmental Engineering Duke Universit Henri Gavin Fall, 1 In situations in which deformations

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

Lecture 7: The Beam Element Equations.

Lecture 7: The Beam Element Equations. 4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite

More information

Aircraft Structures Kirchhoff-Love Plates

Aircraft Structures Kirchhoff-Love Plates University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

Comb resonator design (2)

Comb resonator design (2) Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory

More information

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:

More information

VIBRATION PROBLEMS IN ENGINEERING

VIBRATION PROBLEMS IN ENGINEERING VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor

More information

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness

More information

MEC-E8001 FINITE ELEMENT ANALYSIS

MEC-E8001 FINITE ELEMENT ANALYSIS MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,

More information

Lecture 8: Tissue Mechanics

Lecture 8: Tissue Mechanics Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials

More information

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,

More information

FINAL EXAMINATION. (CE130-2 Mechanics of Materials)

FINAL EXAMINATION. (CE130-2 Mechanics of Materials) UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,

More information

CO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS

CO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL

More information

[8] Bending and Shear Loading of Beams

[8] Bending and Shear Loading of Beams [8] Bending and Shear Loading of Beams Page 1 of 28 [8] Bending and Shear Loading of Beams [8.1] Bending of Beams (will not be covered in class) [8.2] Bending Strain and Stress [8.3] Shear in Straight

More information

CHAPTER 5. Beam Theory

CHAPTER 5. Beam Theory CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e

More information

NONLINEAR CONTINUUM FORMULATIONS CONTENTS

NONLINEAR CONTINUUM FORMULATIONS CONTENTS NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell

More information

Table of Contents. Preface... 13

Table of Contents. Preface... 13 Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...

More information

Lecture 27: Structural Dynamics - Beams.

Lecture 27: Structural Dynamics - Beams. 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

More information

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may

More information

ME FINITE ELEMENT ANALYSIS FORMULAS

ME FINITE ELEMENT ANALYSIS FORMULAS ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness

More information

NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS

NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.

More information

ENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1

ENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1 ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at

More information

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION

More information

Part D: Frames and Plates

Part D: Frames and Plates Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with

More information

Finite Element Method in Geotechnical Engineering

Finite Element Method in Geotechnical Engineering Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps

More information

Symmetric Bending of Beams

Symmetric Bending of Beams Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications

More information

Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum

Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure

More information

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition

Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition Fluid Structure Interaction and Moving Boundary Problems IV 63 Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition K.-H. Jeong, G.-M. Lee, T.-W. Kim & J.-I.

More information

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Mathematical Properties of Stiffness Matrices

Mathematical Properties of Stiffness Matrices Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the

More information

Comb Resonator Design (2)

Comb Resonator Design (2) Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory

More information

Toward a novel approach for damage identification and health monitoring of bridge structures

Toward a novel approach for damage identification and health monitoring of bridge structures Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:

More information

JEPPIAAR ENGINEERING COLLEGE

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

More information

The CR Formulation: BE Plane Beam

The CR Formulation: BE Plane Beam 6 The CR Formulation: BE Plane Beam 6 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction..................... 6 4 6.2 CR Beam Kinematics................. 6 4 6.2. Coordinate

More information

Finite element modelling of structural mechanics problems

Finite element modelling of structural mechanics problems 1 Finite element modelling of structural mechanics problems Kjell Magne Mathisen Department of Structural Engineering Norwegian University of Science and Technology Lecture 10: Geilo Winter School - January,

More information

3D Elasticity Theory

3D Elasticity Theory 3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.

More information

ME 475 Modal Analysis of a Tapered Beam

ME 475 Modal Analysis of a Tapered Beam ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory

More information

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.

More information

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.

More information

CHAPTER -6- BENDING Part -1-

CHAPTER -6- BENDING Part -1- Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and

More information

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics

More information

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix

More information

202 Index. failure, 26 field equation, 122 force, 1

202 Index. failure, 26 field equation, 122 force, 1 Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic

More information

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric

More information

Gyroscopic matrixes of the straight beams and the discs

Gyroscopic matrixes of the straight beams and the discs Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation

More information

CIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass

CIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional

More information

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60. 162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides

More information

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING CHETTINAD COLLEGE OF ENGINEERING AND TECHNOLOGY-KARUR FINITE ELEMENT ANALYSIS 2 MARKS QUESTIONS WITH ANSWER Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING FINITE ELEMENT ANALYSIS

More information

Basic Energy Principles in Stiffness Analysis

Basic Energy Principles in Stiffness Analysis Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting

More information

Bending of Simply Supported Isotropic and Composite Laminate Plates

Bending of Simply Supported Isotropic and Composite Laminate Plates Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,

More information

Section 6: PRISMATIC BEAMS. Beam Theory

Section 6: PRISMATIC BEAMS. Beam Theory Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam

More information

3. The linear 3-D elasticity mathematical model

3. The linear 3-D elasticity mathematical model 3. The linear 3-D elasticity mathematical model In Chapter we examined some fundamental conditions that should be satisfied in the modeling of all deformable solids and structures. The study of truss structures

More information

Dynamic Model of a Badminton Stroke

Dynamic Model of a Badminton Stroke ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815

More information

Aircraft Stress Analysis and Structural Design Summary

Aircraft Stress Analysis and Structural Design Summary Aircraft Stress Analysis and Structural Design Summary 1. Trusses 1.1 Determinacy in Truss Structures 1.1.1 Introduction to determinacy A truss structure is a structure consisting of members, connected

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how

More information

NONLINEAR WAVE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES

NONLINEAR WAVE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES NONLINEAR WAE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES DONGMING WEI Department of Mathematics, University of New Orleans, 2 Lakeshore Dr., New Orleans, LA 7148,USA E-mail: dwei@uno.edu

More information

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts

More information

Introduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil

Introduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Introduction to Waves in Structures Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Waves in Structures Characteristics of wave motion Structural waves String Rod Beam Phase speed, group velocity Low

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,

More information

The Kinematic Equations

The Kinematic Equations The Kinematic Equations David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 0139 September 19, 000 Introduction The kinematic or strain-displacement

More information

Chapter 6: Cross-Sectional Properties of Structural Members

Chapter 6: Cross-Sectional Properties of Structural Members Chapter 6: Cross-Sectional Properties of Structural Members Introduction Beam design requires the knowledge of the following. Material strengths (allowable stresses) Critical shear and moment values Cross

More information

University of Groningen

University of Groningen University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published

More information

MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM

MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,

More information

3. BEAMS: STRAIN, STRESS, DEFLECTIONS

3. BEAMS: STRAIN, STRESS, DEFLECTIONS 3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More information

Geometry-dependent MITC method for a 2-node iso-beam element

Geometry-dependent MITC method for a 2-node iso-beam element Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh

More information

Free Vibration Analysis of Kirchoff Plates with Damaged Boundaries by the Chebyshev Collocation Method. Eric A. Butcher and Ma en Sari

Free Vibration Analysis of Kirchoff Plates with Damaged Boundaries by the Chebyshev Collocation Method. Eric A. Butcher and Ma en Sari Free Vibration Analysis of Kirchoff Plates with Damaged Boundaries by the Chebyshev Collocation Method Eric A. Butcher and Ma en Sari Department of Mechanical and Aerospace Engineering, New Mexico State

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

General elastic beam with an elastic foundation

General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Dynamic analysis of railway bridges by means of the spectral method

Dynamic analysis of railway bridges by means of the spectral method Dynamic analysis of railway bridges by means of the spectral method Giuseppe Catania, Silvio Sorrentino DIEM, Department of Mechanical Engineering, University of Bologna, Viale del Risorgimento, 436 Bologna,

More information

Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section

Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,

More information

A Boundary Integral Formulation for the Dynamic Behavior of a Timoshenko Beam

A Boundary Integral Formulation for the Dynamic Behavior of a Timoshenko Beam A Boundary Integral Formulation for the Dynamic Behavior of a Timoshenko Beam M. Schanz and H. Antes Institute of Applied Mechanics Technical University Braunschweig, Germany e-mail: m.schanz@tu-bs.de,

More information

Outline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices

Outline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices Outline in MDOF Systems Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano May 8, 014 Additional Today we will study the properties of structural matrices, that is the operators that

More information

Structural Matrices in MDOF Systems

Structural Matrices in MDOF Systems in MDOF Systems http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 9, 2016 Outline Additional Static Condensation

More information

A Timoshenko beam theory with pressure corrections for plane stress problems

A Timoshenko beam theory with pressure corrections for plane stress problems A Timoshenko beam theory with pressure corrections for plane stress problems Graeme J. Kennedy a,,, Jorn S. Hansen a,2, Joaquim R.R.A. Martins b,3 a University of Toronto Institute for Aerospace Studies,

More information

7. Hierarchical modeling examples

7. Hierarchical modeling examples 7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in

More information

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017 CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and

More information

Continuum mechanism: Stress and strain

Continuum mechanism: Stress and strain Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the

More information

FREE VIBRATIONS OF FRAMED STRUCTURES WITH INCLINED MEMBERS

FREE VIBRATIONS OF FRAMED STRUCTURES WITH INCLINED MEMBERS FREE VIBRATIONS OF FRAMED STRUCTURES WITH INCLINED MEMBERS A Thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Technology in Civil Engineering By JYOTI PRAKASH SAMAL

More information

Mechanics in Energy Resources Engineering - Chapter 5 Stresses in Beams (Basic topics)

Mechanics in Energy Resources Engineering - Chapter 5 Stresses in Beams (Basic topics) Week 7, 14 March Mechanics in Energy Resources Engineering - Chapter 5 Stresses in Beams (Basic topics) Ki-Bok Min, PhD Assistant Professor Energy Resources Engineering i Seoul National University Shear

More information

Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces

Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces Mechanics and Mechanical Engineering Vol. 18, No. 1 2014) 37 51 c Lodz University of Technology Free Vibrations of Timoshenko Beam with End Mass in the Field of Centrifugal Forces A. I. Manevich V. Yu.

More information

Nonconservative Loading: Overview

Nonconservative Loading: Overview 35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three

More information

Lecture 8: Assembly of beam elements.

Lecture 8: Assembly of beam elements. ecture 8: Assembly of beam elements. 4. Example of Assemblage of Beam Stiffness Matrices. Place nodes at the load application points. Assembling the two sets of element equations (note the common elemental

More information

Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams

Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body

More information

6. Bending CHAPTER OBJECTIVES

6. Bending CHAPTER OBJECTIVES CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where

More information

Steps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen

Steps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,

More information

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ. Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,

More information

Analytical Mechanics: Elastic Deformation

Analytical Mechanics: Elastic Deformation Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda

More information