General elastic beam with an elastic foundation

Size: px
Start display at page:

Download "General elastic beam with an elastic foundation"

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

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 information

10. Applications of 1-D Hermite elements

10. 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 information

Contents as of 12/8/2017. Preface. 1. Overview...1

Contents 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 information

Discontinuous Distributions in Mechanics of Materials

Discontinuous 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 information

3. Numerical integration

3. 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 information

Finite 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 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 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

BEAM DEFLECTION THE ELASTIC CURVE

BEAM 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 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

CIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives

CIVL 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 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

k 21 k 22 k 23 k 24 k 31 k 32 k 33 k 34 k 41 k 42 k 43 k 44

k 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 information

Chapter 4 Deflection and Stiffness

Chapter 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 information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY 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 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

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen

CIV-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 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

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

Transactions 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. 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 information

FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES

FINITE 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 information

MECE 3321: Mechanics of Solids Chapter 6

MECE 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 [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

Eigenvalues of Trusses and Beams Using the Accurate Element Method

Eigenvalues 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 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

Unit - 7 Vibration of Continuous System

Unit - 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 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

Chapter 12 Elastic Stability of Columns

Chapter 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 information

Chapter Objectives. Copyright 2011 Pearson Education South Asia Pte Ltd

Chapter 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 information

M5 Simple Beam Theory (continued)

M5 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 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

9.1 Introduction to bifurcation of equilibrium and structural

9.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 information

2. Polynomial interpolation

2. 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 information

Computational Stiffness Method

Computational 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 information

Methods of Analysis. Force or Flexibility Method

Methods 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 information

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14

Chapter 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 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

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

Mechanics of Inflatable Fabric Beams

Mechanics 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 information

Structural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.

Structural 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 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

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

Materials: engineering, science, processing and design, 2nd edition Copyright (c)2010 Michael Ashby, Hugh Shercliff, David Cebon.

Materials: 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 information

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS

Bilinear 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 information

Application of Finite Element Method to Create Animated Simulation of Beam Analysis for the Course of Mechanics of Materials

Application 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 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

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

Interpolation Functions for General Element Formulation

Interpolation 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 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

Stress analysis of a stepped bar

Stress 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 information

Laboratory 4 Topic: Buckling

Laboratory 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 information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL 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 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

Free 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 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 information

Converting Plane Stress Statics to 2D Natural Frequencies, changes in red

Converting 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 information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method

Module 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 information

structural analysis Excessive beam deflection can be seen as a mode of failure.

structural 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 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

Chapter 7: Internal Forces

Chapter 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 information

LATERAL STABILITY OF BEAMS WITH ELASTIC END RESTRAINTS

LATERAL 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 information

Chapter 2: Deflections of Structures

Chapter 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 information

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering

Deflection 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 information

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.

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. 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 information

Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering

Deflection 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 information

Chapter 3 Higher Order Linear ODEs

Chapter 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 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

1859. 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 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 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

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

Strength 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 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 information

Module 3. Analysis of Statically Indeterminate Structures by the Displacement Method

Module 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 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

Mechanical Design in Optical Engineering

Mechanical 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 information

MODULE C: COMPRESSION MEMBERS

MODULE 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 information

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations

ENGI 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 information

1 Static Plastic Behaviour of Beams

1 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 information

2 marks Questions and Answers

2 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 information

Chapter 8 Supplement: Deflection in Beams Double Integration Method

Chapter 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 information

Mechanics of Materials Primer

Mechanics 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 information

NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil

NON 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 information

EML4507 Finite Element Analysis and Design EXAM 1

EML4507 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 information

Finite 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 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 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

CH. 4 BEAMS & COLUMNS

CH. 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 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

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

CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions

CHAPTER 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 information

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).

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). 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 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

Vibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1

Vibration 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 information

Advanced 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 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 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

Chapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship

Chapter 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 information

Free vibration and stability of axially functionally graded tapered Euler-Bernoulli beams

Free 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 information

Variational Formulation of Plane Beam Element

Variational 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 information

7.4 The Elementary Beam Theory

7.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 information

Chapter 12. Static Equilibrium and Elasticity

Chapter 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 information

MECHANICS OF STRUCTURES SCI 1105 COURSE MATERIAL UNIT - I

MECHANICS 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 information

Analysis III for D-BAUG, Fall 2017 Lecture 11

Analysis 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 information

FEA 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 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 information

STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS

STRESS 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