General elastic beam with an elastic foundation
|
|
- Brandon Bailey
- 5 years ago
- Views:
Transcription
1 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 spring has a known displacement, v f, which is almost always zero. The differential equation of equilibrium of an initially straight beam, of flexural stiffness EI, resting on a Winkler foundation, of stiffness k per unit length, with a transverse load per unit length of w, subjected to a tensile axial load N acting along the x-axis is: d 2 [EI d2 v dx 2 dx 2] N d2 v dx2 + kv w(x) =. (1) where v(x) is the transverse displacement of the beam. If the axial force, N, is present, then this is usually called a Beam-Column model. If the axial load is an unknown constant, then it is a buckling load that has to be found by an eigen-solution rather than by solving a linear algebraic system. Beams on elastic foundations are fairly common (like pipelines, rail lines or drill strings). However, most introductions to beams assume no axial load, nor any foundation support. Then, the above general beam equation reduces to that covered in an introduction to solid mechanics: d 2 [EI d2 v dx 2 dx2] = w(x). (2) The homogeneous solution of this equation is simply a cubic polynomial with four constants. Either beam equation is a fourth-order ordinary differential equation. Therefore it will generally need four boundary conditions. Figure 1 A beam-column on an elastic foundation Related physical quantities are the slope, θ(x) = v (x), the bending moment, M(x) = EIv (x), and the transverse shear force, V(x) = EIv (x). The sign conventions are that the position, x, is positive to the right, the deflection, v, point forces, P, and the line load, w, are positive in the y-direction (upward), and the slopes and moments are positive in the counter-clockwise direction. Since the foundation springs are restrained, there can be no rigid body motions of the beam. In other words, the assembled equations should never be singular and it is acceptable to just have non-essential boundary conditions applied to the beam. A simple Winkler foundation model like this one can push or pull on the beam as needed and no gaps can occur. Copyright J.E. Akin. All rights reserved. Page 1 of 11
2 The distributed load per unit length can include point transverse shear loads, V, by using the Dirac Delta distribution. Likewise, employing a doublet distribution in defining w(x) allows for the inclusion of point couples, or moments, M. Engineers designing beams are usually interested in the local moment and transvers shear force since they define the stress levels and the material failure criteria. The exact solution of the homogeneous (w==) general form given above is given in terms of hyperbolic sines and cosines. Based on Tong s Theorem, exact solutions at the nodes are obtained if such functions are used as the interpolation functions in a finite element model. Advanced elements of that type have been applied with excellent results. To apply the Galerkin weak form the governing ODE is multiplied by v(x) and the integral over the length of the beam is set to zero. The highest derivative term is always present, and needs to be integrated twice by parts to reduce the inter-element continuity requirement. That integral becomes L I E = v d2 dx 2 (E(x)I(x) d2 v dx2) dx I E = [v d dx (E(x)I(x) d2 L v dx 2)] [ dv dx (E(x)I(x) d2 L v dx 2)] + d2 v dx 2 Copyright J.E. Akin. All rights reserved. Page 2 of 11 L (E(x)I(x) d2 v dx2) dx which brings the non-essential boundary conditions into the integral form. The first term is the product of the displacement and the point shear force at both ends, while the second term is the product of the slope and the point moment at both ends. Thus, the contributions from the highest derivative term in the ODE are I E = [v V(x)] L [ dv L dx M(x)] L + d2 v dx 2 (E(x)I(x) d2 v dx2) dx Therefore, the integral form contains second derivatives (d 2 v dx 2 ) as its highest derivative term. The presence of second derivatives in the integral form means that calculus requires the elements to have inter-element continuity of the deflection and the slope (v(x) and θ(x) = v (x)). Such elements are said to have an inter-element continuity of C 1. Therefore, each shared beam node must have two degrees of freedom (dof), at least. The presence of the second derivatives in the integral form also means that the essential boundary conditions (EBC) will specify v or θ or both at a point. The non-essential boundary conditions (NBC) are that the transverse shear force, V = EIv, and/or the moment, M = EIv, is given at a point. When the deflection is given, then the shear force is a reaction quantity. When the slope is specified, then the moment is a reaction. Calculus continuity requires the elements share their slope and deflection values. Thus, classic beam elements approximated by at least a cubic polynomial having C 1 continuity can give nodally exact deflections and slopes. One-dimensional polynomial interpolation with C 1 or C 2 continuity usually is
3 based on Hermite interpolation polynomials. The classic beam element is often said to be a Hermite cubic. Some typical Hermite polynomials are given in the Appendix. A line element with only two nodes will provide the four constants needed to define a cubic polynomial. Such elements are referred to here as classic beam elements. But since designers are interested in the second and third derivatives of the solution, classic beam elements give very poor design results unless a large number of them are used. All of the matrices for a cubic beam are well known and are just summarized here for comparison to the three-node quintic beam. For a typical beam (k = ) with a cross-sectional moment of inertia, I, length, L, a depth of h, and material with an elastic modulus of E and a coefficient of thermal expansion of α, the corresponding flexural stiffness transverse load and thermal load matrices for the equilibrium of a single classic beam element, are EI [ L L 12 6L v 1 V 1 6L 4L 2 6L 2L 2 θ ] { 1 M 12 6L 12 6L v } = { 1 } + 2 V 2 6L 2L 2 6L 4L 2 θ 2 M L 3L 6 9 2L 21 2L 3L Copyright J.E. Akin. All rights reserved. Page 3 of 11 ] { w 1 w } + α TEI { 2 h 1 }, (3) where w 1 and w 2 are the load per unit length at the first (left) and second end, respectively, and T is the temperature increase, over h, from the top to the bottom of the beam. Also, the V j and M j terms are point shear forces and point moments applied to the nodes as external loadings and/or unknown reactions. In matrix symbol notation this equation is written as K e v e = F P e + F w e + F α e. (4) The classic beam element has a deflection given by (for r = x L): A slope of v(x) = v 1 (1 3r 2 + 2r 3 ) + θ 1 (r 2r 2 + r 3 )L + v 2 (3r 2 2r 3 ) + θ 2 (r 3 r 2 )L (5) θ(x) = v 1 (6r 2 6r)/L + θ 1 (1 4r + 3r 2 ) + v 2 (6r 6r 2 )/L + θ 2 (3r 2 2r) (6) A linear bending moment of M(x) = [v 1 (12r 6) + θ 1 (6r 4)L + v 2 (6 12r) + θ 2 (6r 2)L]/L 2 (7) and a constant transverse shear force of V(x) = [v 1 (12) + θ 1 (6)L + v 2 ( 12) + θ 2 (6)L]/L 3. (8) If a Winkler foundation is present then the foundation stiffness is For the classic cubic beam this becomes 1 K e f = k e H e dl (9) L e H e T
4 156 22L e K e f = ρe A e L e 22L e L e 2 13L e 13L e 3L e L e 13L e 156 3L e 2 22L e ] (1) 22L e 4L e 2 If an axial load is present, the geometric stiffness matrix after integration by parts is K e = H e T L e N e d2 e H dx 2 dl K e N = Ne 3L e [ 36 3L e 3L e 4L e L e 3L e L e L e 3L e 36 L e 2 3L e ]. (11) 3L e 4L e 2 For dynamic response or natural frequency calculations it is necessary to include the physical mass matrix defined as For the classic cubic beam this becomes m e = H e T L e ρ e A e H e dl L e m e = ρe A e L e 22L e L e 2 13L e 13L e 3L e L e 13L e 156 3L e 2 22L e ]. (12) 22L e 4L e 2 Many exact solutions for beams not on an elastic foundation are tabulated in typical structural design manuals. They show that for the most common load and support conditions the exact solution is a third, fourth or fifth degree polynomial. Thus, a fifth degree polynomial will generally give a very accurate or exact beam solution with only a few elements, perhaps with only one element. However, if a foundation or an axial load is present polynomial approximations will no longer be exact at the nodes. The Hermite family of interpolation polynomials can be increased in order by increasing the number of nodes and/or by increasing the continuity level at each node. Here, a more accurate three-node C 1 beam element will be presented. This will be created by adding a third (middle) node to the element. That adds two more degrees of freedom and raises the polynomial to a fifth degree (a quintic polynomial). For simplicity, all of the coefficients in Eq. (1) are taken as constants evaluated at the center of the element length. That allows the corresponding finite element matrices of the fifth-degree element to be expressed in closed form in terms of the element deflection and slope at each of its three nodes [v 1 θ 1 v 2 θ 2 v 3 θ 3 ]. The element interpolation relations and the corresponding slope, moment and shear distributions are: Copyright J.E. Akin. All rights reserved. Page 4 of 11
5 v(r) = [v 1 (1 23r r 3 68r r 5 ) + θ 1 (r 6r r 3 12r 4 + 4r 5 )L +v 2 (16r 2 32r r 4 ) + θ 2 ( 8r r 3 4r r 5 )L +v 3 (7r 2 34r r 4 24r 5 ) + θ 3 ( r 2 + 5r 3 8r 4 + 4r 5 )L] (13) θ(r) = v (r) = [v 1 ( 46r + 198r 2 272r r 4 )/L + θ 1 (1 12r + 39r 2 48r 3 + 2r 4 )] +v 2 (32r 96r r 3 )/L + θ 2 ( 16r + 96r 2 16r 3 + 8r 4 ) +v 3 (14r 12r r 3 12r 4 )/L + θ 3 ( 2r + 15r 2 32r 3 + 2r 4 )] (14) M(r) = EIv (r) = EI[v 1 ( r 816r r 3 )/L + θ 1 ( r 144r 2 + 8r 3 ) +v 2 (32 192r + 192r 2 )/L + θ 2 ( r 48r r 3 ) +v 3 (14 24r + 624r 2 48r 3 )/L + θ 3 ( 2 + 3r 96r 2 + 8r 3 )]/L (15) V(r) = EIv (r) = EI[v 1 (396 1,632r + 1,44r 2 )/L + θ 1 (78 288r + 24r 2 ) +v 2 ( r)/L + θ 2 (192 96r + 96r 2 ) +v 3 ( ,248r 1,35r 2 )/L + θ 3 (3 192r + 24r 2 )]/L 2 (16) The 6 by 6 element matrix equations of equilibrium of a single element are: (K E e + K N e + K k e ){v e } = {F e } (17) where {v e } T = [v 1, θ 1, v 2, θ 2, v 3, θ 3 ] are the generalized nodal displacements and rotations, and the {F e } are the generalized resultant nodal loading forces and moments. The K e terms are stiffness matrices arising from the first three terms in the differential equation. Usually, only the first is nonzero. This system is singular until enough boundary conditions are applied to prevent rigid body translations and rotations. Since this element has three nodes, the loading per unit length can be input at each of those nodes. The resultant element external loading vectors due to point forces, V k, and couples, M k, and the quadratic distributed load values, w k, at each node (or constant w 1 = w 2 = w 3, or linear w 2 = (w 1 + w 3 ) 2 values) are F P e = V 1 M 1 V 2 M 2 V 3 { M 3 }, F w e = H et H w w e L e dl = L 14,7 [ 1,995 1, L 14L 56 6, L 28L 15 1,54 1,995 14L 15L] w 1 { w 2 } (18) w 3 where the H e are the six Hermite beam interpolations, and the H w are the three quadratic Lagrangian interpolations for the distributed loads. The integral of their product forms a rectangular six by three array to convert distributed loads to concentrated shears and moments at the beam nodes. Likewise, it can be shown that a temperature difference through the depth of a beam causes only end moment loadings given by F e α = α T EI h [ 1 1]. A constant transverse load per unit length reduces the resultant load and moment vector to Copyright J.E. Akin. All rights reserved. Page 5 of 11
6 F w T = wl 14,7 [3,43 245L 7,84 3,43 245L], for w 1 = w 2 = w 3 = w (19) The symmetric flexural stiffness matrix for the three noded quintic element is K E = EI 35L 3 5,92 1,138L 3,584 1,138L 332L 2 896L 3, L 7,168 1,92L 32L 2 1,58 242L 3,584 [ 242L 38L 2 896L 1,92L 1,58 242L 32L 2 242L 38L 2 3, L 1,28L 2 1,92L 32L 2 1,92L 5,92 1,138L 32L 2 1,138L 332L 2 ] (2) The symmetric beam-column (or geometric stiffness) matrix due to any axial load is K N = N 63 L 1, L 1, L L 39 L 28 L 2 48 L 8 L 2 9 L 5 L 2 1, L 3,72 1, L 24 L 8 L L 2 24 L 8 L L 1, L 1, L [ 9 L 5 L 2 48 L 8 L 2 39 L 28 L 2 ] (21) The geometric stiffness matrix always depends on the existing stresses in the element. Here that is just a uniform axial stress, N / A. There is a general equation for the geometric stiffness matrix, involving all the stresses at a point, which comes from a non-linear analysis. It will not be given here. The element foundation stiffness is K k = k L 13,86 2, L L L 114 L 8 L 2 88 L 12 L 2 29 L 3 L L 5, L 16 L 12 L L 2 16 L 12 L L L 2, L [ 29 L 3 L 2 88 L 12 L L 8 L 2 ] (22) If the mass matrix, m e, is included for dynamic response or vibration analysis then the above matrix is the element mass matrix when k is replaced by ρa where ρ is the mass density and A is the crosssectional area of the beam (so ρa is the mass-per unit length). Since it is common for beams to have piecewise constant properties, the classic beam element and the more accurate quintic beam element can be programmed in closed form. Using either the classic or quintic beam element some small beam analysis problems can also be solved in closed analytic form. As an example of the three-node beam, consider a fixed-fixed beam with a constant line load. The essential boundary conditions are v 1 =, θ 1 =, v 3 =, θ 3 =. The external point force and moment at center node 2 are zero (V 2 =, M 2 = ). The middle two rows define the remaining unknown center point deflection and slope: Copyright J.E. Akin. All rights reserved. Page 6 of 11
7 EI [7,168 35L 3 1,28L 2] {v 2 θ } = { 2 } + wl {7,84 14,7 } + { }. (23) Multiplying by the inverse of the square matrix gives the middle node solutions: { v 2 θ 2 } = 35wL4 14,7EI 7,168 [1 1 1,28L 2 ] { 7,84 } = wl4 {1 384EI } (24) which are exact. The slope, θ 2, was expected to be zero due to symmetry. It could have been used as a boundary condition to compute v 2 from one equation. Here, the reactions on the left are found from the first two rows of the equilibrium equations (since all the displacements are now known): EI [ 3,584 35L 3 896L 1,92L 32L 2 ] wl4 {1 384EI } = { V 1 } + wl M {3, ,7 245L } (25) { V 1 } = wl M { L 6 }. Likewise, at the right end, utilizing the last two rows of the equilibrium equations { V 3 } = wl M { 1 }. (26) 3 2 L 6 Note that the total resultant force is wl, which is equal and opposite to the applied transverse load. The net external moment is zero. Usually static equilibrium is taught in undergraduate classes with Newton s Laws. They can used to verify that the computed reactions do indeed satisfy that the sum of the forces is zero, and that the sum of the moments, taken at any reference point, is zero. Since a 5-th degree element was used this problem is exact everywhere. That includes the moment and shear diagrams. The exact shear diagrams of EIv can have a discontinuity (due to point loads P) and the exact moment diagram of EIv can have a discontinuity (due to point moments M). The mesh should be constructed such that any point source occurs at an interface node between elements. Those items can be discontinuous at element interfaces because the calculus requirement for splitting integrals on says that v and v must be continuous there. That explains why a two node C 2 Hermite interpolation was not chosen to create a quintic beam. That choice could be made only if the beam was not allowed to have point sources. Inside any element the solution is continuous and thus is C at any non-interface node. In other words, inside any standard element you can calculate an infinite number of continuous derivatives (but the vast majority is zero). If you were to apply a point source at an interior node then the second and third derivatives would still be continuous and the shear and moment diagrams in that element would be wrong. Next consider the involvement of the foundation in the post-processing needed to recover the reactions on each beam element. Before the deflections of a beam on an elastic foundation (BOEF) Copyright J.E. Akin. All rights reserved. Page 7 of 11
8 are computed, a typical element free body diagram would be similar to Figure 1. As before, any external point load, P, and/or any point couple, M, should be placed at an element interface node and they go directly into the first or second equilibrium row of the node at that point. A distributed line load is integrated to become a consistent load vector. Those two effects are detailed in Eq. 14. Both the beam stiffness and the foundation stiffness resist the external applied loads. The element equilibrium equation (for no axial force) is (K E e + K k e ){v e } = {F P e } + {F w e } (27) The assembled equations have a similar form and are solved for the system displacements, v. Returning to the element to get the beam member reactions, the now known element degrees of freedom, v e, are gathered. The free body diagram of the element now has a significant change because the foundation is now applying a known pressure that opposes the displacement of the beam. That pressure, p(x) = k f v(x), is a polynomial of the same degree as the assumed displacements. That state is shown in Figure 2, which defines the sought element reaction vector, {F R e }. Now, there is a new consistent load vector to evaluate during post processing to find the element reactions. It comes from the foundation pressure and is proportional to the deflection: F e p = H e T p(x)dx = L e H e T k L e f H e dx v e = K e f v e. (28) The foundation stiffness caused the new load, but that load and others must be resisted by only the beam stiffness. Now the beam member equilibrium equations define the reactions on the nodes of the beam, {F R e }: or K E e {v e } = {F R e } + {F p e } + {F P e } + {F w e } (29) K E e {v e } {F p e } {F P e } {F w e } = {F R e }. (3) Now that there is more than one dof per node the scatter/gather operations previously used for axial models have become more generalized. The usual convention for numbering the unknowns is to count all of the unknowns at a node before moving on to another node. Therefore, the displacements are associated with the odd number of rows and columns in a matrix, while the slopes contributions are occurring in the even ones. For previous one-dimensional scalar unknowns scatter simply meant adding a scalar term at any node into the system equilibrium matrix and source vector. Now, that changes to adding a two by two sub-matrix from each element node into the system square matrix and a two by one column matrix into the system source vector. The prior algorithms for doing the scatter and gather still work because they had the number of degrees of freedom per node (n g ) as an input argument. It is just for hand solutions that you need to clearly understand this change. Figure 3 Copyright J.E. Akin. All rights reserved. Page 8 of 11
9 graphically illustrates the sub-matrices scatter operations for a line element with two unknowns per node. To determine the natural frequencies of a beam-column on an elastic foundation the eigen-problem is (K E + K N + K k ) ω 2 M = (31) To determine the natural frequencies of a beam-column the eigen-problem is (K E + K N ) ω 2 M = (32) When the axial load is an unknown constant, say N B, then it factors out of each element matrix and becomes a global constant: where K n is formed with a unit positive load. (K E + N B K n ) ω 2 M = (33) This shows that any axial load, N B, has an influence on the natural frequency, ω. In general, a tensile force (+) increases the natural frequency while a compression force (-) lowers the natural frequency. When the global constant, N B, is unknown it is called the buckling factor of safety (BFS) and the equations lead to a buckling eigen-problem to determine the BFS. The BFS is the factor of safety against buckling, or the ratio of the buckling loads to the applied loads. The following table gives the interpretation of possible values for the buckling factor of safety: Since the geometric stiffness matrix depends on the stress that means that it depends on all of the assembled loads, say F ref. In general, a linear static analysis is completed and the stresses are determined: (K E + K k ){v } = F ref (34) Then the geometric stiffness matrix is calculated from those stresses. The stress, and thus the geometric stiffness matrix, is proportion to the resultant load, so as the loads are scaled by the buckling factor the geometric stiffness increases by the same amount F BFS F ref, K N BFSK ref The scaling value that renders the combined stiffness to have a zero determinant (that is to become unstable) is calculated from a buckling eigen-problem as: K E + K k + BFS K ref = (35) This is solved for the value of the BFS. Then the load that would theoretically cause buckling is BFS F ref. The solution of the eigen-problem also yields the relative buckling mode shape (eigen- Copyright J.E. Akin. All rights reserved. Page 9 of 11
10 vector), v BFS. The magnitude of the mode shape displacements are arbitrary and most commercial software normalizes them to range from to 1. BFS Status Note >1 Buckling not predicted The applied loads are less than the theoretical critical loads. 1 Buckling predicted The applied loads exceed the theoretical critical loads. 1 Buckling predicted The applied loads are exactly equal to the theoretical critical loads. -1 Buckling not predicted The buckling occurs when the directions of the applied loads are all reversed. -1 Buckling not predicted The buckling occurs when the directions of the applied loads are all reversed. < -1 Buckling not predicted Buckling is not predicted even if you reverse all loads. Figure 2 Beam element reaction vector, including foundation pressure Copyright J.E. Akin. All rights reserved. Page 1 of 11
11 Figure 3 Scattering sub-matrices for beam elements with two dof per node Copyright J.E. Akin. All rights reserved. Page 11 of 11
Quintic beam closed form matrices (revised 2/21, 2/23/12) 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 information10. Applications of 1-D Hermite elements
10. Applications of 1-D Hermite elements... 1 10.1 Introduction... 1 10.2 General case fourth-order beam equation... 3 10.3 Integral form... 5 10.4 Element Arrays... 7 10.5 C1 Element models... 8 10.6
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationDiscontinuous Distributions in Mechanics of Materials
Discontinuous Distributions in Mechanics of Materials J.E. Akin, Rice University 1. Introduction The study of the mechanics of materials continues to change slowly. The student needs to learn about software
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationFinite 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 informationBEAM DEFLECTION THE ELASTIC CURVE
BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment
More informationChapter 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 informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationCOPYRIGHTED 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 informationk 21 k 22 k 23 k 24 k 31 k 32 k 33 k 34 k 41 k 42 k 43 k 44
CE 6 ab Beam Analysis by the Direct Stiffness Method Beam Element Stiffness Matrix in ocal Coordinates Consider an inclined bending member of moment of inertia I and modulus of elasticity E subjected shear
More informationChapter 4 Deflection and Stiffness
Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationLecture 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 informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationInstitute 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 informationCHAPTER -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 informationTransactions of the VŠB Technical University of Ostrava, Mechanical Series. article No Roland JANČO *
Transactions of the VŠB Technical University of Ostrava, Mechanical Series No. 1, 013, vol. LIX article No. 1930 Roland JANČO * NUMERICAL AND EXACT SOLUTION OF BUCKLING LOAD FOR BEAM ON ELASTIC FOUNDATION
More informationFINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES
th International Conference on Earthquake Geotechnical Engineering June 5-8, 7 Paper No. 11 FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES A.Halim KARAŞĐN 1, Polat GÜLKAN ABSTRACT Plates on elastic foundations
More informationMECE 3321: Mechanics of Solids Chapter 6
MECE 3321: Mechanics of Solids Chapter 6 Samantha Ramirez Beams Beams are long straight members that carry loads perpendicular to their longitudinal axis Beams are classified by the way they are supported
More information[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 informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationStructural 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 informationUnit - 7 Vibration of Continuous System
Unit - 7 Vibration of Continuous System Dr. T. Jagadish. Professor for Post Graduation, Department of Mechanical Engineering, Bangalore Institute of Technology, Bangalore Continuous systems are tore which
More informationIntroduction 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 informationChapter 12 Elastic Stability of Columns
Chapter 12 Elastic Stability of Columns Axial compressive loads can cause a sudden lateral deflection (Buckling) For columns made of elastic-perfectly plastic materials, P cr Depends primarily on E and
More informationChapter Objectives. Copyright 2011 Pearson Education South Asia Pte Ltd
Chapter Objectives To generalize the procedure by formulating equations that can be plotted so that they describe the internal shear and moment throughout a member. To use the relations between distributed
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationMechanics 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 information9.1 Introduction to bifurcation of equilibrium and structural
Module 9 Stability and Buckling Readings: BC Ch 14 earning Objectives Understand the basic concept of structural instability and bifurcation of equilibrium. Derive the basic buckling load of beams subject
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
More information6. 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 information3. 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 informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationMechanics 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 informationLecture 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 informationMaterials: engineering, science, processing and design, 2nd edition Copyright (c)2010 Michael Ashby, Hugh Shercliff, David Cebon.
Modes of Loading (1) tension (a) (2) compression (b) (3) bending (c) (4) torsion (d) and combinations of them (e) Figure 4.2 1 Standard Solution to Elastic Problems Three common modes of loading: (a) tie
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationApplication of Finite Element Method to Create Animated Simulation of Beam Analysis for the Course of Mechanics of Materials
International Conference on Engineering Education and Research "Progress Through Partnership" 4 VSB-TUO, Ostrava, ISSN 156-35 Application of Finite Element Method to Create Animated Simulation of Beam
More informationCHAPTER 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 informationLecture 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 informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationConsider 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 informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationLaboratory 4 Topic: Buckling
Laboratory 4 Topic: Buckling Objectives: To record the load-deflection response of a clamped-clamped column. To identify, from the recorded response, the collapse load of the column. Introduction: Buckling
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationGeometric 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 informationFree vibration analysis of elastically connected multiple-beams with general boundary conditions using improved Fourier series method
Free vibration analysis of elastically connected multiple-beams with general boundary conditions using improved Fourier series method Jingtao DU*; Deshui XU; Yufei ZHANG; Tiejun YANG; Zhigang LIU College
More informationConverting Plane Stress Statics to 2D Natural Frequencies, changes in red
The following illustrates typical modifications for converting any plane stress static formulation into a plane stress natural frequency and mode shape calculation. The changes and additions to a prior
More informationModule 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method
Module 2 Analysis of Statically Indeterminate Structures by the Matrix Force Method Lesson 8 The Force Method of Analysis: Beams Instructional Objectives After reading this chapter the student will be
More informationstructural analysis Excessive beam deflection can be seen as a mode of failure.
Structure Analysis I Chapter 8 Deflections Introduction Calculation of deflections is an important part of structural analysis Excessive beam deflection can be seen as a mode of failure. Extensive glass
More informationPart 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 informationChapter 7: Internal Forces
Chapter 7: Internal Forces Chapter Objectives To show how to use the method of sections for determining the internal loadings in a member. To generalize this procedure by formulating equations that can
More informationLATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS
LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS By John J. Zahn, 1 M. ASCE ABSTRACT: In the analysis of the lateral buckling of simply supported beams, the ends are assumed to be rigidly restrained
More informationChapter 2: Deflections of Structures
Chapter 2: Deflections of Structures Fig. 4.1. (Fig. 2.1.) ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 1 (2.1) (4.1) (2.2) Fig.4.2 Fig.2.2 ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 2
More informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 008/9 Dr. Colin Caprani 1 Contents 1. Introduction... 3 1.1 General... 3 1. Background... 4 1.3 Discontinuity Functions...
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 009/10 Dr. Colin Caprani 1 Contents 1. Introduction... 4 1.1 General... 4 1. Background... 5 1.3 Discontinuity Functions...
More informationChapter 3 Higher Order Linear ODEs
Chapter 3 Higher Order Linear ODEs Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 3.1 Homogeneous Linear ODEs 3 Homogeneous Linear ODEs An ODE is of
More informationChapter 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 information1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load
1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical
More informationIntroduction 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 informationSection 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 informationStrength of Materials Prof. S.K.Bhattacharya Dept. of Civil Engineering, I.I.T., Kharagpur Lecture No.26 Stresses in Beams-I
Strength of Materials Prof. S.K.Bhattacharya Dept. of Civil Engineering, I.I.T., Kharagpur Lecture No.26 Stresses in Beams-I Welcome to the first lesson of the 6th module which is on Stresses in Beams
More informationModule 3. Analysis of Statically Indeterminate Structures by the Displacement Method
odule 3 Analysis of Statically Indeterminate Structures by the Displacement ethod Lesson 21 The oment- Distribution ethod: rames with Sidesway Instructional Objectives After reading this chapter the student
More informationCHAPTER 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 informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationMODULE C: COMPRESSION MEMBERS
MODULE C: COMPRESSION MEMBERS This module of CIE 428 covers the following subjects Column theory Column design per AISC Effective length Torsional and flexural-torsional buckling Built-up members READING:
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More information1 Static Plastic Behaviour of Beams
1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial
More information2 marks Questions and Answers
1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and
More informationChapter 8 Supplement: Deflection in Beams Double Integration Method
Chapter 8 Supplement: Deflection in Beams Double Integration Method 8.5 Beam Deflection Double Integration Method In this supplement, we describe the methods for determining the equation of the deflection
More informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More informationNON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil
Abstract: NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil The geometric non-linear total potential energy equation is developed and extended to study the behavior of buckling
More informationEML4507 Finite Element Analysis and Design EXAM 1
2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 17
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras Module - 01 Lecture - 17 In the last class, we were looking at this general one-dimensional
More informationPrepared 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 informationCH. 4 BEAMS & COLUMNS
CH. 4 BEAMS & COLUMNS BEAMS Beams Basic theory of bending: internal resisting moment at any point in a beam must equal the bending moments produced by the external loads on the beam Rx = Cc + Tt - If the
More informationME 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 informationLecture 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 informationCHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions
1. Deflections of Beams and Shafts CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 1. Integration method. Discontinuity functions 3. Method
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationJEPPIAAR 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 informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationAdvanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras
Advanced Structural Analysis Prof. Devdas Menon Department of Civil Engineering Indian Institute of Technology, Madras Module No. # 5.1 Lecture No. # 27 Matrix Analysis of Beams and Grids Good morning,
More informationDynamic 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 informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationFree vibration and stability of axially functionally graded tapered Euler-Bernoulli beams
Shock and Vibration 18 211 683 696 683 DOI 1.3233/SAV-21-589 IOS Press Free vibration and stability of axially functionally graded tapered Euler-Bernoulli beams Ahmad Shahba a,b, Reza Attarnejad a,b, and
More informationVariational Formulation of Plane Beam Element
13 Variational Formulation of Plane Beam Element IFEM Ch 13 Slide 1 Beams Resist Primarily Transverse Loads IFEM Ch 13 Slide 2 Transverse Loads are Transported to Supports by Flexural Action Neutral surface
More information7.4 The Elementary Beam Theory
7.4 The Elementary Beam Theory In this section, problems involving long and slender beams are addressed. s with pressure vessels, the geometry of the beam, and the specific type of loading which will be
More informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationMECHANICS OF STRUCTURES SCI 1105 COURSE MATERIAL UNIT - I
MECHANICS OF STRUCTURES SCI 1105 COURSE MATERIAL UNIT - I Engineering Mechanics Branch of science which deals with the behavior of a body with the state of rest or motion, subjected to the action of forces.
More informationAnalysis III for D-BAUG, Fall 2017 Lecture 11
Analysis III for D-BAUG, Fall 2017 Lecture 11 Lecturer: Alex Sisto (sisto@math.ethz.ch) 1 Beams (Balken) Beams are basic mechanical systems that appear over and over in civil engineering. We consider a
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
More informationSTRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS
1 UNIT I STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define: Stress When an external force acts on a body, it undergoes deformation. At the same time the body resists deformation. The
More information