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


 Stuart Fletcher
 1 years ago
 Views:
Transcription
1 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the workequivalence method for replacing distributed loading by a set of discrete loads To introduce the general formulation for solving beam problems with distributed loading acting on them To analyze beams with distributed loading acting on them To compare the finite element solution to an exact solution for a beam To derive the stiffness matrix for the beam element with nodal hinge To show how the potential energy method can be used to derive the beam element equations To apply Galerkin s residual method for deriving the beam element equations General Formulation We can account for the distributed loads or concentrated loads acting on beam elements by considering the following formulation for a general structure: F=KdF where F are the equivalent nodal forces, expressed in terms of the globalcoordinate components. These forces would yield the same displacements as the original distributed load. If we assume that the global nodal forces are not present (F = ) then: F =Kd
2 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 General Formulation We now solve for the displacements, d, given the nodal forces F. Next, substitute the displacements and the equivalent nodal forces F back into the original expression and solve for the global nodal forces. F=KdF This concept can be applied on a local basis to obtain the local nodal forces in individual elements of structures as: f=kdf Example 5  oad Replacement Consider the beam shown below; determine the equivalent nodal forces for the given distributed load. The work equivalent nodal forces are shown above. Using the beam stiffness equations: w f y 6 6 v w m fy 6 6v w m w
3 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Example 5  oad Replacement Apply the boundary conditions: v 6 6 v v We can solve for the displacements w 6v w w v 8 w 6 Example 5  oad Replacement In this case, the method of equivalent nodal forces gives the exact solution for the displacements and rotations. To obtain the global nodal forces, we will first define the product of Kd to be F e, where F e is called the effective global nodal forces. Therefore: w e F 6 6 y e 5w M e 4 F 6 6 w y 8 w e M w 6 w
4 CIV 7/87 Chapter 4  Development of Beam Equations  Part 4/4 Example 5  oad Replacement Using the above expression and the fixend moments in: w w w F y 5w w w M F=KdF F y w w M w w Fy M Fy w w w Example 6  Cantilever Beam Consider the beam, shown below, determine the vertical displacement and rotation at the freeend and the nodal forces, including reactions. Assume is constant throughout the beam. We will use one element and replace the concentrated load with the appropriate nodal forces.
5 CIV 7/87 Chapter 4  Development of Beam Equations  Part 5/4 Example 6  Cantilever Beam The beam stiffness equations are: P fy 6 6v P m fy 6 6v P m P 8 Apply the boundary conditions: 6 6 v v v Example 6  Cantilever Beam The beam stiffness equations become: P 6v P P v 48 P 8 To obtain the global nodal forces, we begin by evaluating the effective nodal forces. e F 6 6 y e M e F P y 48 e M P 8 P P 8 P P 8
6 Example 6  Cantilever Beam Using the above expression in the following equation, gives: F=KdF 8 8 P P P P e Kd = F F F P P 8 8 P P P P y y F M F M F Example 6  Cantilever Beam In general, for any structure in which an equivalent nodal force replacement is made, the actual nodal forces acting on the structure are determined by first evaluating the effective nodal forces F e for the structure and then subtracting off the equivalent nodal forces F for the structure. Similarly, for any element of a structure in which equivalent nodal force replacement is made, the actual local nodal forces acting on the element are determined by first evaluating the effective local nodal forces for the element and then subtracting off the equivalent local nodal forces associated only with the element. f=kdf CIV 7/87 Chapter 4  Development of Beam Equations  Part 6/4
7 CIV 7/87 Chapter 4  Development of Beam Equations  Part 7/4 Comparison of FE Solution to Exact Solution We will now compare the finite element solution to the exact classical beam theory solution for the cantilever beam shown below. Both one and twoelement finite element solutions will be presented and compared to the exact solution obtained by the direct doubleintegration method. et E = x 6 psi, I = in 4, = in, and uniform load w = Ib/in. Comparison of FE Solution to Exact Solution To obtain the solution from classical beam theory, we use the doubleintegration method: Mx ( ) y where the double prime superscript indicates differentiation twice with respect to x and M is expressed as a function of x by using a section of the beam as shown: F y V( x) w wx M wx w Mx ( ) wx
8 CIV 7/87 Chapter 4  Development of Beam Equations  Part 8/4 Comparison of FE Solution to Exact Solution To obtain the solution from classical beam theory, we use the doubleintegration method: Mx ( ) y w x y x dxdx w x x x y C dx Boundary Conditions 6 y() y() 4 w x x x y Cx C w x x x y Comparison of FE Solution to Exact Solution Recall the oneelement solution to the cantilever beam is: 4 w v 8 w 6 Using the numerical values for this problem we get: 4 lb in in 6 4 v 8 psi in.8in lb in in. rad psi in
9 CIV 7/87 Chapter 4  Development of Beam Equations  Part 9/4 Comparison of FE Solution to Exact Solution The slope and displacement from the oneelement FE solution identically match the beam theory values evaluated at x =. The reason why these nodal values from the FE solution are correct is that the element nodal forces were calculated on the basis of being energy or work equivalent to the distributed load based on the assumed cubic displacement field within each beam element. Comparison of FE Solution to Exact Solution Values of displacement and slope at other locations along the beam for the FE are obtained by using the assumed cubic displacement function. vx ( ) x xv xx The value of the displacement at the midlength v(x = 5 in) is: v( x 5 in).78in Using beam theory, the displacement at v(x = 5 in) is: v( x 5 in).95in
10 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Comparison of FE Solution to Exact Solution In general, the displacements evaluated by the FE method using the cubic function for v are lower than by those of beam theory except at the nodes. This is always true for beams subjected to some form of distributed load that are modeled using the cubic displacement function. The exception to this result is at the nodes, where the beam theory and FE results are identical because of the workequivalence concept used to replace the distributed load by workequivalent discrete loads at the nodes. Comparison of FE Solution to Exact Solution The beam theory solution predicts a quartic (fourthorder) polynomial expression for a beam subjected to uniformly distributed loading, while the FE solution v(x) assumes a cubic (thirdorder) displacement behavior in each beam all load conditions. The FE solution predicts a stiffer structure than the actual one. This is expected, as the FE model forces the beam into specific modes of displacement and effectively yields a stiffer model than the actual structure. However, as more elements are used in the model, the FE solution converges to the beam theory solution.
11 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Comparison of FE Solution to Exact Solution For the special case of a beam subjected to only nodal concentrated loads, the beam theory predicts a cubic displacement behavior. The FE solution for displacement matches the beam theory solution for all locations along the beam length, as both v(x) and y(x) are cubic functions. Comparison of FE Solution to Exact Solution Under uniformly distributed loading, the beam theory solution predicts a quadratic moment and a linear shear force in the beam. However, the FE solution using the cubic displacement function predicts a linear bending moment and a constant shear force within each beam element used in the model.
12 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Comparison of FE Solution to Exact Solution We will now determine the bending moment and shear force in the present problem based on the FE method. Mx ( ) Mx ( ) d y dx B d Nd d dx Nd 6 x 4 6x 6 x 6x B 6 x 4 6x M v 6 x 6x v Comparison of FE Solution to Exact Solution We will now determine the bending moment and shear force in the present problem based on the FE method. Position M FE M exact X = 8, lbin , lbin X = 5 in , lbin 5, lbin X = in 6,667 lbin
13 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Comparison of FE Solution to Exact Solution The plots below show the displacement, bending moment, and shear force over the beam using beam theory and the oneelement FE solutions. Comparison of FE Solution to Exact Solution The FE solution for displacement matches the beam theory solution at the nodes but predicts smaller displacements (less deflection) at other locations along the beam length.
14 CIV 7/87 Chapter 4  Development of Beam Equations  Part 4/4 Comparison of FE Solution to Exact Solution The bending moment is derived by taking two derivatives on the displacement function. It then takes more elements to model the second derivative of the displacement function. Comparison of FE Solution to Exact Solution The shear force is derived by taking three derivatives on the displacement function. For the uniformly loaded beam, the shear force is a constant throughout the singieelement model.
15 CIV 7/87 Chapter 4  Development of Beam Equations  Part 5/4 Comparison of FE Solution to Exact Solution To improve the FE solution we need to use more elements in the model (refine the mesh) or use a higherorder element, such as a fifthorder approximation for the displacement function. Beam Element with Nodal Hinge Consider the beam, shown below, with an internal hinge. An internal hinge causes a discontinuity in the slope of the deflection curve at the hinge and the bending moment is zero at the hinge.
16 CIV 7/87 Chapter 4  Development of Beam Equations  Part 6/4 Beam Element with Nodal Hinge For a beam with a hinge on the right end: [ K] [ K] fy 6 6v m fy 6 6v m [ ] [ K ] K The moment m is zero and we can partition the matrix to eliminate the degree of freedom associated with. f K Kd f KdKd f KK d f K d K d Beam Element with Nodal Hinge For a beam with a hinge on the right end: fy 6 6v m fy 6 6v m d v v d f f y m f y f m d K f K d K d K K f K d f Kd Kd f c kcd fc fkk f m c k K K K K
17 CIV 7/87 Chapter 4  Development of Beam Equations  Part 7/4 Beam Element with Nodal Hinge We can condense out the degree of freedom by using the partitioning method discussed earlier. k K K K K c 6 6 kc Therefore, the condensed stiffness matrix is: kc Beam Element with Nodal Hinge The element forcedisplacement equations are: f y v m f y v Expanding the element forcedisplacement equations and maintaining m = gives: fy v m fy v m
18 CIV 7/87 Chapter 4  Development of Beam Equations  Part 8/4 Beam Element with Nodal Hinge Once the displacements d are found, d may be computed: d K f K d m v K m K v In this case, d is (the rotation at right side of the element at the hinge. Beam Element with Nodal Hinge For a beam with a hinge on the left end: fy 6 6v m fy 6 6v m Rewriting the equations to move and m to the first row: m f y 6 6 v fy 6 6 v m 6 6 4
19 CIV 7/87 Chapter 4  Development of Beam Equations  Part 9/4 Beam Element with Nodal Hinge For a beam with a hinge on the left end: [ K] [ K] m f y 6 6 v fy 6 6 v m K [ ] [ K ] The moment m is zero and we can partition the matrix to eliminate the degree of freedom associated with. f K Kd f KdKd f KK d f K d K d Beam Element with Nodal Hinge For a beam with a hinge on the right end: m f 6 6 v f 6 6 v m d K f K d d f m d v v f fy f m y f K d K d K d K K f K d m c kcd fc f KK f kc K KK K f
20 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Beam Element with Nodal Hinge We can condense out the degree of freedom by using the partitioning method discussed earlier. k K K K K c 6 6 kc Therefore, the condensed stiffness matrix is: k c Beam Element with Nodal Hinge The element forcedisplacement equations are: fy v f y v m Expanding the element forcedisplacement equations and maintaining m = gives: fy v m fy v m
21 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Beam Element with Nodal Hinge Once the displacements d are found, d may be computed: d K f K d m v K m K v In this case, d is (the rotation at right side of the element at the hinge. Beam Element with Nodal Hinge For a beam element with a hinge at its left end, the element forcedisplacement equations are: fy v f y v m Expanding the element forcedisplacement equations and maintaining m = gives: fy v m fy v m
22 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Example 7  Beam With Hinge In the following beam, shown below, determine the vertical displacement and rotation at node and the element forces for the uniform beam with an internal hinge at node. Assume is constant throughout the beam. Example 7  Beam With Hinge The stiffness matrix for element (with hinge on right) is: fy v m fy v m v v fy a v m a a a fy a a v m
23 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Example 7  Beam With Hinge The stiffness matrix for element (with no hinge) is: f 6 6v m f 6 6v m v v f 6b 6bv m 6b 4b 6b b f b 6b 6b v m 6b b 6b 4b Example 7  Beam With Hinge The assembled equations are: Element f a a a y v m a a a 6 6 f y a a b a b b b v m b b b b f 6 6 y v b b b b m b b b b Element
24 CIV 7/87 Chapter 4  Development of Beam Equations  Part 4/4 Example 7  Beam With Hinge The boundary conditions are: v v f a a a y v m a a a 6 6 f y a a a b b b b v m b b b b f 6 6 y v b b b b m b b b b Example 7  Beam With Hinge After applying the boundary conditions the global beam equations reduce to: 6 a b b v P 6 4 b b abp v b a abp b a
25 CIV 7/87 Chapter 4  Development of Beam Equations  Part 5/4 Example 7  Beam With Hinge The element forcedisplacement equations for element are: f y a m a a a a f y a abp b a bp b a f y ab P m b a f y bp b a Example 7  Beam With Hinge The slope on element may be found using the condensed matrix: m v d K f Kd K m K element v a 4 a abp b a 6a a 6a abp b a
26 CIV 7/87 Chapter 4  Development of Beam Equations  Part 6/4 Example 7  Beam With Hinge The element forcedisplacement equations for element are: abp f y 6b 6b b a m 6b 4b 6b b abp f y b 6b 6b b a m 6b b 6b 4b ap b a fy m f y ap m b a ba P b a Example 7  Beam With Hinge Displacements and rotations on each element. v abp v v b a abp b a abp b a
27 CIV 7/87 Chapter 4  Development of Beam Equations  Part 7/4 Potential Energy Approach to Derive Beam Element Equations et s derive the equations for a beam element using the principle of minimum potential energy. The procedure for applying the principle of minimum potential energy is similar to that used for the bar element. The total potential energy p is defined as the sum of the internal strain energy U and the potential energy of the external forces : U p Potential Energy Approach to Derive Beam Element Equations The differential internal work (strain energy) du in a onedimensional beam element is: U x x dv For a single beam element, shown below, subjected to both distributed and concentrated nodal forces, the potential energy due to forces (or the work done these forces) is: S V TvdS Pv m y iy i i i i i
28 CIV 7/87 Chapter 4  Development of Beam Equations  Part 8/4 Potential Energy Approach to Derive Beam Element Equations If the beam element has a constant crosssectional area A, then the differential volume of the beam is given as: dv da dx The differential element where the surface loading acts is given as: ds = b dx (where b is the width of the beam element). Therefore the total potential energy is: da dx T vb dx P v m p x x y iy i i i xa x i Potential Energy Approach to Derive Beam Element Equations dv The straindisplacement relationship is: x y dx We can express the strain in terms of nodal displacements and rotations as: x 6 6x4 x 6 6x x y d [ ] x yb d x 6 6x4 x 6 6x [ B]
29 CIV 7/87 Chapter 4  Development of Beam Equations  Part 9/4 Potential Energy Approach to Derive Beam Element Equations The stressstrain relationship in onedimension is: [ E] x x where E is the modulus of elasticity. Therefore: [ ][ ] x ye B d [ ] The total potential energy can be written in matrix form as: T T T da dx bt v dx d P p x x y xa x x yb d Potential Energy Approach to Derive Beam Element Equations If we define, w bt y as a line load (load per unit length) in the y direction and the substitute the definitions of x and x the total potential energy can be written in matrix form as: E p y d B B d dadx w d N dx d P T T T T T [ ] [ ] [ ] A Use the following definition for moment of inertia: I Then the total potential energy expression becomes: p T T T T T d [ B] [ B] ddx wd [ N] dx d P A y da
30 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Potential Energy Approach to Derive Beam Element Equations Differentiating the total potential energy with respect to the displacement and rotations (v, v, and ) and equating each term to zero gives: The nodal forces vector is: T T [ B] [ B] dx d w[ N] dx P T f w[ N] dx P The elemental stiffness matrix is: T k [ B] [ B] dx Potential Energy Approach to Derive Beam Element Equations Integrating the previous matrix expression gives: T k [ B] [ B] dx k
31 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Galerkin s Method to Derive Beam Element Equations The governing differential equation for a onedimensional beam is: 4 dv w 4 dx We can define the residual R as: 4 dv w N,,,and4 4 i dx i dx If we apply integration by parts twice to the first term we get: xxxx i xx i, xx i xxx i, x xx v N dx v N dx N v N v where the subscript x indicates a derivative with respect to x Galerkin s Method to Derive Beam Element Equations Since v N d to x is: or v xx xx [ ], then the second derivative of v with respect x 6 6x4 x 6 6x v B d [ ] Therefore the integration by parts becomes: Nixx, B dxd Nw i dx NV i Nix, m d d [ ] i,,, and 4
32 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Galerkin s Method to Derive Beam Element Equations The above expression is really four equations (one for each N i ) and can be written in matrix form as: x T T T T [ B] [ B] dx d [ N] w dx [ N] m [ N] V The interpolation function in the last two terms can be evaluated: [ N] x [ ] [ N] x [ ] x [ N] x [ ] [ N] x [ ] x Galerkin s Method to Derive Beam Element Equations Therefore, the last two terms of the matrix form of the Galerkin formulation become (see the figure below): i V() i m() i V( ) i 4 m( )
33 CIV 7/87 Chapter 4  Development of Beam Equations  Part /4 Galerkin s Method to Derive Beam Element Equations When beam elements are assembled, as shown below: Two shear forces and two moments form adjoining elements contribute to the concentrated force and the concentrated moment at the node common to both elements. Problems: 7. Verify the four beam element equations are contained in the following matrix expression. T T [ B] [ B] dx d w[ N] dx P 8. Do problems 4., 4., 4.8, 4., 4.4 and 4.47 in your textbook.
34 CIV 7/87 Chapter 4  Development of Beam Equations  Part 4/4 Problems: 9. Work problem 4.6 in your using the SAP. Attempt to select the lightest standard W section to support the loads for the beam. The bending stress must not exceed 6 MPa and the allowable deflection must not exceed ( = 6 m)/6 End of Chapter 4b
CIVL 8/7117 Chapter 12  Structural Dynamics 1/75. To discuss the dynamics of a singledegreeof freedom springmass
CIV 8/77 Chapter  /75 Introduction To discuss the dynamics of a singledegreeof freedom springmass system. To derive the finite element equations for the timedependent stress analysis of the onedimensional
More informationCIVE1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIVE16 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 informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beamcolumn on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
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 informationQuintic 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 beamcolumn on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
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 crosssectional area of a beam. Supports that apply a moment
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 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 informationEML4507 Finite Element Analysis and Design EXAM 1
21715 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 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 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 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 informationFinite Element MethodPart II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element MethodPart II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationDevelopment of Beam Equations
CHAPTER 4 Development of Beam Equations Introduction We begin this chapter by developing the stiffness matrix for the bending of a beam element, the most common of all structural elements as evidenced
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 informationENGN2340 Final Project: Implementation of a EulerBernuolli Beam Element Michael Monn
ENGN234 Final Project: Implementation of a EulerBernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the EulerBernuolli
More informationFINAL EXAMINATION. (CE1302 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130 Mechanics of Materials) Problem 1: (15 points) pinned bar structure is shown in Figure 1. There is an external force, W = 5000N,
More 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 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 informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture  06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture  06 In the last lecture, we have seen a boundary value problem, using the formal
More informationIntroduction to Structural Member Properties
Introduction to Structural Member Properties Structural Member Properties Moment of Inertia (I): a mathematical property of a crosssection (measured in inches 4 or in 4 ) that gives important information
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 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 informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More information10. Applications of 1D Hermite elements
10. Applications of 1D Hermite elements... 1 10.1 Introduction... 1 10.2 General case fourthorder beam equation... 3 10.3 Integral form... 5 10.4 Element Arrays... 7 10.5 C1 Element models... 8 10.6
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationChapter 11. Displacement Method of Analysis Slope Deflection Method
Chapter 11 Displacement ethod of Analysis Slope Deflection ethod Displacement ethod of Analysis Two main methods of analyzing indeterminate structure Force method The method of consistent deformations
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 crosssectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
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 informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3  Truss Equations  Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and ComputerAided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More information= 50 ksi. The maximum beam deflection Δ max is not = R B. = 30 kips. Notes for Strength of Materials, ET 200
Notes for Strength of Materials, ET 00 Steel Six Easy Steps Steel beam design is about selecting the lightest steel beam that will support the load without exceeding the bending strength or shear strength
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 informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline EulerBernoulli Beams Assumptions
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 informationSoftware Verification
EXAMPLE 1026 FRAME MOMENT AND SHEAR HINGES EXAMPLE DESCRIPTION This example uses a horizontal cantilever beam to test the moment and shear hinges in a static nonlinear analysis. The cantilever beam has
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 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 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 informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationFinite element modelling of structural mechanics problems
1 Finite element modelling of structural mechanics problems Kjell Magne Mathisen Department of Structural Engineering Norwegian University of Science and Technology Lecture 10: Geilo Winter School  January,
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More 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 informationGeometrydependent MITC method for a 2node isobeam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 33 Geometrydependent MITC method for a node isobeam element PhillSeung Lee Samsung Heavy Industries, Seocho, Seoul 37857, Korea HyuChun Noh
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 58, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationHIGHERORDER THEORIES
HIGHERORDER THEORIES Thirdorder Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
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 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 informationHIGHERORDER THEORIES
HIGHERORDER THEORIES THIRDORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 ThirdOrder Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
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 informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
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 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 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 informationTruss Structures: The Direct Stiffness Method
. Truss Structures: The Companies, CHAPTER Truss Structures: The Direct Stiffness Method. INTRODUCTION The simple line elements discussed in Chapter introduced the concepts of nodes, nodal displacements,
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 informationChapter 3 Higher Order Linear ODEs
Chapter 3 Higher Order Linear ODEs Advanced Engineering Mathematics WeiTa 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 informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November  2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN
More informationTwo Tier projects for students in ME 160 class
ME 160 Introduction to Finite Element Method Spring 2016 Topics for Term Projects by Teams of 2 Students Instructor: Tai Ran Hsu, Professor, Dept. of Mechanical engineering, San Jose State University,
More informationVibration of a Euler Bernoulli uniform beam carrying a rigid body at each end
Vibration of a Euler Bernoulli uniform beam carrying a rigid body at each end S. Naguleswaran Department of Mechanical Engineering, University of Canterbury, Christchurch, New Zealand EMail: sivapatham.naguleswaran@canterbury.ac.nz
More informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic StressStrain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic StressStrain Relationship A stress in the xdirection causes a strain in the xdirection by σ x also causes a strain in the ydirection & zdirection
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 informationBending Load & Calibration Module
Bending Load & Calibration Module Objectives After completing this module, students shall be able to: 1) Conduct laboratory work to validate beam bending stress equations. 2) Develop an understanding of
More informationPhysics 8 Wednesday, November 29, 2017
Physics 8 Wednesday, November 29, 2017 HW11 due this Friday, Dec 1. After another day or two on beams, our last topic of the semester will be oscillations (a.k.a. vibration, periodic motion). Toward that
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pinjointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationCHAPTER 6 BENDING Part 1
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER 6 BENDING Part 11 CHAPTER 6 Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationPhysics 8 Monday, November 23, 2015
Physics 8 Monday, November 23, 2015 Handing out HW11, due Friday, December 4. One or two more beamrelated examples, then we ll move on to oscillations ( periodic motion ). This week s reading is Mazur
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T  I INTRODUCTION
More informationMECE8001 FINITE ELEMENT ANALYSIS
MECE800 FINIE EEMEN ANAYSIS 07  WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationStrength of Materials Prof: S.K.Bhattacharya Dept of Civil Engineering, IIT, Kharagpur Lecture no 28 Stresses in Beams III
Strength of Materials Prof: S.K.Bhattacharya Dept of Civil Engineering, IIT, Kharagpur Lecture no 28 Stresses in Beams III Welcome to the 3 rd lesson of the 6 th module which is on Stresses in Beams part
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
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 informationTypes of Structures & Loads
Structure Analysis I Chapter 4 1 Types of Structures & Loads 1Chapter Chapter 4 Internal lloading Developed in Structural Members Internal loading at a specified Point In General The loading for coplanar
More information5. What is the moment of inertia about the x  x axis of the rectangular beam shown?
1 of 5 Continuing Education Course #274 What Every Engineer Should Know About Structures Part D  Bending Strength Of Materials NOTE: The following question was revised on 15 August 2018 1. The moment
More informationUniversity of Illinois at UrbanaChampaign College of Engineering
University of Illinois at UrbanaChampaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationName :. Roll No. :... Invigilator s Signature :.. CS/B.TECH (CENEW)/SEM3/CE301/ SOLID MECHANICS
Name :. Roll No. :..... Invigilator s Signature :.. 2011 SOLID MECHANICS Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks. Candidates are required to give their answers
More informationUNSYMMETRICAL BENDING
UNSYMMETRICAL BENDING The general bending stress equation for elastic, homogeneous beams is given as (II.1) where Mx and My are the bending moments about the x and y centroidal axes, respectively. Ix and
More informationSolution: The moment of inertia for the crosssection is: ANS: ANS: Problem 15.6 The material of the beam in Problem
Problem 15.4 The beam consists of material with modulus of elasticity E 14x10 6 psi and is subjected to couples M 150, 000 in lb at its ends. (a) What is the resulting radius of curvature of the neutral
More information3.5 Reinforced Concrete Section Properties
CHAPER 3: Reinforced Concrete Slabs and Beams 3.5 Reinforced Concrete Section Properties Description his application calculates gross section moment of inertia neglecting reinforcement, moment of inertia
More informationCOLUMNS: BUCKLING (DIFFERENT ENDS)
COLUMNS: BUCKLING (DIFFERENT ENDS) Buckling of Long Straight Columns Example 4 Slide No. 1 A simple pinconnected truss is loaded and supported as shown in Fig. 1. All members of the truss are WT10 43
More informationErrata Sheet for S. D. Rajan, Introduction to Structural Analysis & Design (1 st Edition) John Wiley & Sons Publication
S D Rajan, Introduction to Structural Analsis & Design ( st Edition) Errata Sheet for S D Rajan, Introduction to Structural Analsis & Design ( st Edition) John Wile & Sons Publication Chapter Page Correction
More informationFEA CODE WITH MATLAB. Finite Element Analysis of an Arch ME 5657 FINITE ELEMENT METHOD. Submitted by: ALPAY BURAK DEMIRYUREK
FEA CODE WITH MATAB Finite Element Analysis of an Arch ME 5657 FINITE EEMENT METHOD Submitted by: APAY BURAK DEMIRYUREK This report summarizes the finite element analysis of an archbeam with using matlab.
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 informationPart 1 is to be completed without notes, beam tables or a calculator. DO NOT turn Part 2 over until you have completed and turned in Part 1.
NAME CM 3505 Fall 06 Test 2 Part 1 is to be completed without notes, beam tables or a calculator. Part 2 is to be completed after turning in Part 1. DO NOT turn Part 2 over until you have completed and
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More information3 2 6 Solve the initial value problem u ( t) 3. a If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b True or false and why 1. if A is
More informationPhysics 8 Monday, November 20, 2017
Physics 8 Monday, November 20, 2017 Pick up HW11 handout, due Dec 1 (Friday next week). This week, you re skimming/reading O/K ch8, which goes into more detail on beams. Since many people will be traveling
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 informationNumerical and Analytical Approach of ThermoMechanical Stresses in FGM Beams
Numerical and Analytical Approach of ThermoMechanical Stresses in FGM Beams Fatemeh Farhatnia, GholamAli Sharifi and Saeid Rasouli Abstract In this paper, thermomechanical stress distribution has been
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 informationCIVL4332 L1 Introduction to Finite Element Method
CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such
More informationTRIANGULAR, NINEDEGREESOFFREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY*
QUARTERLY OF APPLIED MATHEMATICS 303 OCTOBER, 1973 TRIANGULAR, NINEDEGREESOFFREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY* By ISAAC FRIED and SHOK KENG YANG (Boston University) Introduction.
More information(2) ANALYSIS OF INDETERMINATE STRUCTRES
Chapter (2) ANALYSIS OF INDETERMINATE STRUCTRES 1.1 Statically Indeterminate Structures A structure of any type is classified as statically indeterminate when the number of unknown reaction or internal
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 informationApplications in Solid Mechanics
Companies, 4 CHAPTER 9 Applications in Solid 9. INTRODUCTION The bar and beam elements discussed in Chapters 4 are line elements, as only a single coordinate axis is required to define the element reference
More information