The CR Formulation: BE Plane Beam
|
|
- Corey Caldwell
- 6 years ago
- Views:
Transcription
1 6 The CR Formulation: BE Plane Beam 6
2 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction CR Beam Kinematics Coordinate Systems Degrees of Freedom Partial Derivatives Arbitrary Base Configuration Stress Resultants The Deformational Strain Energy Internal Force Vector and Tangent Stiffness Matrix Internal Force Vector Material Stiffness Matrix Geometric Stiffness Matrix Notes and Bibliography Exercises
3 TABE OF CONTENTS Page 6. Introduction CR Beam Kinematics Coordinate Systems Degrees of Freedom Partial Derivatives Arbitrary Base Configuration Stress Resultants The Deformational Strain Energy Internal Force Vector and Tangent Stiffness Matrix Internal Force Vector Material Stiffness Matrix Geometric Stiffness Matrix Notes and Bibliography Exercises
4 Chapter 6: THE CR FORMUATION: BE PANE BEAM 6.. Introduction This Chapter uses the CR kinematic description to construct a geometrically nonlinear, 2-node, Bernoulli-Euler plane beam. Unlike the T beam element presented in Chapters 0-, we can do a C (Hermitian) beam from scratch, without need of MacNeal s RBF trick to finish the job. The reason is that with the CR approach it is as easy to do C (Bernoulli-Euler) or C 0 (Timoshenko) beams, and the former has a better performing geometric stiffness matrix CR Beam Kinematics The CR formulation of the beam motion is quite similar to that of the bar element in many respects, and much of the development can be reused. Only the major differences will be noted here Coordinate Systems As in Chapter 0, we consider a plane, straight, prismatic beam element with two end nodes. The element is initially aligned with the global X axis in the initial configuration C 0, with the origin O 0 located at the element midpoint. This configuration is assumed to be straight and undeformed although it may be under initial uniform axial stress with resultant N 0. The bar properties include the elastic modules, E, the cross section area A 0 and the moment of inertia I 0 about the neutral axis. The reference length is 0. The motion on the {X, Y } plane carries it to the current configuration C. The corotated configuration C R is selected as depicted in Figure 6... The longitudinal axis passes through the current position of the end nodes. This defines the local axis x e. The origin of {x e, y e } is placed halfway between the nodes. This forms an angle ψ with X. 2. The C R nodes are placed at an equal distance from the C nodes. Hence the corotated axes {xr e, ye R }, including origin, coincide with {x e, y e }. The new ingredient is the rotation angle θ about Z or z. With C R chosen as indicated, the deformation part of these rotations is easily extracted: θ = θ ψ. Other choices for C R are possible. The foregoing one has the advantages of simplicity, and of being compatible with that of the bar element discussed in the previous Chapter Degrees of Freedom The CR plane beam element has six degrees of freedom, which are placed in the vectors u X ū e X u Y ū e 2 d Y 0 θ u =, ū e θ e θ = u X2 ū e = X2 u Y 2 ū e 2 d. (6.) Y 2 0 θ 2 θ 2 e θ 2 See Figure 6.2 for a picture of the global displacements and Figure 6.3 for the deformational displacements. 6 4
5 Current configuration // Y _ C // Y _ θ CR BEAM KINEMATICS φ = ψ+ϕ ψ // X _ ϕ 2(x 2,y 2) // X u Y2 Y, y Base or initial configuration C 0 X, x Y _ γ (X,Y ) θ (x,y ) u X u Y Corotated configuration X _ 2(X 2,Y 2) ϕ u X2 C R Reference length 0, cross section area A 0, second moment of inertia I 0, and elastic modulus E Figure 6.. Kinematics of corotational C plane beam element. Deformations in current configuration grossly exaggerated for visibility. 2 θ 2 //X, x C u X θ //X, x u Y2 u Y Y, y X, x C u X2 Figure 6.2. Global displacements of CR plane beam element initially aligned with X axis. y e d/2 θ _ C R d/2 2 _ θ 2 x e 0 C Figure 6.3. Deformational displacements of CR plane beam element in its local (CR) system. Deformations grossly exaggerated for visibility. 6 5
6 Chapter 6: THE CR FORMUATION: BE PANE BEAM Proceeding as in the general formulation specialized to the 2D case, we can obtain the following relation: ū e x c ψ s ψ u X u X0 2 0( c ψ ) ū e y s ψ c ψ u Y u Y 0 2 ū e θ e θ 0s ψ ψ = ū e = + x c ψ s ψ 0 u X2 u X0 2 ū e s y2 ψ c ψ 0 u Y 2 u Y 0 0(c ψ ) (6.2) 2 0s ψ θ 2 e θ 2 ψ Here c ψ and s ψ and the angle ψ are implicitly defined by the displacements through the trigonometric relations s ψ = sin ψ = y, c ψ = cos ψ = x, ψ = arctan y (6.3) x where x = 0 + u X2 u X, y = u Y 2 u Y, and = 2 x + 2 y (6.4) is the bar length in the current configuration, ignoring the bending deformation. We note the following relations = = c ψ, u X2 u X c ψ u X2 = c ψ u X = s2 ψ, = = s ψ, u Y 2 u Y c ψ u Y 2 = c ψ u Y = s ψc ψ, c ψ θ s ψ = s ψ = s ψc ψ u X2 u X, s ψ = s ψ = c2 ψ u Y 2 u Y, ψ = ψ = s ψ u X2 u X, ψ = ψ = c ψ u Y 2 u Y, which are useful in the calculations that follow Partial Derivatives = = 0, θ θ 2 = c ψ θ 2 = 0, s ψ = s ψ = 0, θ θ 2 ψ = ψ = 0. θ θ 2 (6.5) The first and second partial derivatives of the deformations d, θ and θ 2 with respect to the node displacements are necessary for the computations of internal forces and stiffness matrices. Using (6.5) and Mathematica, one obtains for the first derivatives: δū e x δū e 2 c ψ 2 s ψ 0 2 c ψ 2 s ψ 0 y s ψ c ψ 0 / cψ 2 δ θ e 0/ 0 s ψ c ψ 0 / cψ 2 0/ 0 s ψ / c ψ / s ψ / c ψ / 0 δū e = x2 δū e 2 c ψ 2 s ψ 0 2 c ψ 2 s ψ 0 y2 s ψ c ψ 0 / cψ 2 0/ 0 s ψ c ψ 0 / cψ 2 0/ 0 s ψ / c ψ / 0 s ψ / c ψ / δ θ e 2 δu X δu Y δθ δu X2 δu Y 2 δθ 2 (6.6) 6 6
7 6.2 CR BEAM KINEMATICS θ 2 // Y _ φ = ψ+ϕ ψ // X _ ϕ 2(x 2,y 2) // X // Y _ C u Y 2 θ (x,y ) _ X Y, y X, x _ Y C 0 u Y 2(X,Y ) 2 2 u X2 ϕ u X (X,Y ) Figure 6.4. Beam element with arbitrarily oriented base configuration C 0, forming an angle ϕ with X. Corotated configuration not shown to reduce clutter. Deformations in current configuration grossly exaggerated for visibility. Since u e x = ue x2 = 2 d, θ e = θ, and θ e 2 = θ 2 we get δd c ψ s ψ 0 c ψ s ψ 0 δ θ = s ψ / c ψ / s ψ / c ψ / 0 δ θ 2 s ψ / c ψ / 0 s ψ / c ψ / δu X δu Y δθ δu X2 δu Y 2 δθ 2 The second derivatives of deformation variables are sψ 2 s ψ c ψ 0 s 2 ψ s ψ c ψ 0 2 d u u = s ψ c ψ cψ 2 0 s ψ c ψ cψ 2 0 sψ 2 s ψ c ψ 0 sψ 2 s ψ c ψ 0 s ψ c ψ cψ 2 0 s ψ c ψ cψ 2 0 2s ψ c ψ cψ 2 s2 ψ 0 2s ψ c ψ sψ 2 c2 ψ 0 2 θ u u = cψ 2 s2 ψ 2s ψ c ψ 0 sψ 2 c2 ψ 2s ψ c ψ 0 2 2s ψ c ψ sψ 2 c2 ψ 0 2s ψ c ψ cψ 2 s2 ψ 0 sψ 2 c2 ψ 2s ψ c ψ 0 cψ 2 s2 ψ 2s ψ c ψ 0 (6.7) (6.8) (6.9) 6 7
8 Chapter 6: THE CR FORMUATION: BE PANE BEAM 2 θ 2 u u = Arbitrary Base Configuration 2s ψ c ψ cψ 2 s2 ψ 0 2s ψ c ψ sψ 2 c2 ψ 0 cψ 2 s2 ψ 2s ψ c ψ 0 sψ 2 c2 ψ 2s ψ c ψ 0 2s ψ c ψ sψ 2 c2 ψ 0 2s ψ c ψ cψ 2 s2 ψ 0 sψ 2 c2 ψ 2s ψ c ψ 0 cψ 2 s2 ψ 2s ψ c ψ 0 (6.0) The foregoing relations can be generalized to the case of a base (initial) configuration C 0 not aligned with the X axis as shown in Figure 6.4. Given the node coordinates and displacements shown in that figure, it is easily shown (cf. Chapter 0) that cos ϕ = X 2 / 0, sin ϕ = Y 2 / 0, cos φ = cos(ψ + ϕ) = x 2 /, sin φ = sin(ψ + ϕ) = y 2 /, cos ψ = (X 2 x 2 + Y 2 y 2 )/( 0 ) and sin ψ = (X 2 y 2 Y 2 x 2 )/( 0 ). The preceding transformation rules remain correct if ψ is replaced by φ = ϕ + ψ, except for the deformation angle computation, which remain θ = θ ψ and θ 2 = θ 2 ψ because the θs are measured from X. The relation between deformational and global displacements become d = 0 = u X2 c φ + u Y 2 s φ + 0 ( c φ ) θ = θ ψ (6.) θ 2 = θ 2 ψ The first derivatives of deformation variables, expressed as variations to get a neast matrix form, are δu X δd c φ s φ 0 c φ s φ 0 δu Y δ θ = s φ / c φ / s φ / c φ / 0 δθ (6.2) δu δ θ 2 s φ / c φ / 0 s φ / c φ / X2 δu Y 2 δθ 2 The second derivatives of deformation variables are sφ 2 s φ c φ 0 s 2 φ s φ c φ 0 2 d u u = s φ c φ cφ 2 0 s φ c φ cφ 2 0 sφ 2 s φ c φ 0 sφ 2 s φ c φ 0 (6.3) s φ c φ cφ 2 0 s φ c φ cφ 2 0 2s φ c φ cφ 2 s2 φ 0 2s φ c φ sφ 2 c2 φ 0 2 θ u u = cφ 2 s2 φ 2s φ c φ 0 sφ 2 c2 φ 2s φ c φ 0 2 2s φ c φ sφ 2 c2 φ 0 2s φ c φ cφ 2 s2 φ 0 (6.4) sφ 2 c2 φ 2s φ c φ 0 cφ 2 s2 φ 2s φ c φ 0 6 8
9 6.4 INTERNA FORCE VECTOR ANDTANGENT STIFFNESS MATRIX 2 θ 2 u u = Stress Resultants 2s φ c φ cφ 2 s2 φ 0 2s φ c φ sφ 2 c2 φ 0 cφ 2 s2 φ 2s φ c φ 0 sφ 2 c2 φ 2s φ c φ 0 2s φ c φ sφ 2 c2 φ 0 2s φ c φ cφ 2 s2 φ 0 sφ 2 c2 φ 2s φ c φ 0 cφ 2 s2 φ 2s φ c φ 0 (6.5) The stress resultants in the reference configuration (either C 0 or C R ) are N 0, M 0 and M0 2. The initial shear force is V 0 = (M 0 M0 2 )/ 0. See Figure 6.5 for sign conventions. Denote by N, V and M the stress resultants in the current configuration. Whereas N and V are constant along the element, M = M(x e ) varies linearly along the length because this is a Hermitian or model, which relies on cubic transverse displacements. Consequently we will define its variation by the two node values M and M 2. The shear V is recovered from equilibrium as V = (M M 2 )/, which is also constant. The stress resultants can be obtained from the deformations as N = N 0 + EA 0 0 d, M = M 0 2EI 0 0 (2 θ + θ 2 ), M 2 = M EI 0 ( θ + 2 θ 2 ), V = M M 2 0 = V EI 0 ( θ θ 2 ). (6.6) 6.3. The Deformational Strain Energy The next step in the CR formulation is to work out the deformational strain energy of the beam. The basic choices are:. A linear beam 2. A nonlinear T beam In the interest of resuing the linear element, we select the first choice. The strain energy of the beam for small strains can be written U = U a + U b + U g (6.7) where U a, U b and U g are the energy taken by axial (bar) deformation, bending deformation, and initial-stress geometric effects, respectively. We adopt the following energy expressions: U a = N 0 d + 2 (N N 0 )d 2 = N 0 0 e + 2 EA 0 0 e 2, [ ] T [ ][ ] θ U b = M2 0 θ 2 M 0 θ + EI θ, 2 θ θ 2 U g = 2 [ θ θ 2 ] T N [ ][ ] 4 θ. 4 θ 2 (6.8) The 2 2 matrices appearing in U b and U g may be derived from those given in Chapters 5 and 5, respectively, of Przemieniecki s book. This book, howevr, omits the initial stress terms. J. S. Przemieniecki, Theory of Matrix Structural Analysis, Dover, New York,
10 Chapter 6: THE CR FORMUATION: BE PANE BEAM M 2 V N C M V N C 0 V 0 M 0 2 N 0 N 0 V 0 M 0 Figure 6.5. Plane beam stress resultants displaying positive sign conventions. Axial forces N and transverse shear forces V are constant along the length, but bending moments M vary linearly. Hence two nodal values of M are required Internal Force Vector and Tangent Stiffness Matrix The internal force vector and tangent stiffness matrix of the corrotational element are then obtained by the usual formulas: p = U u, K = p u = K M + K G (6.9) To develop these quantities it is necessary to find the first and second partial derivatives of d, θ and θ 2 in terms of the node displacements Internal Force Vector Using the partial derivatives compiled above and Mathematica, one obtains the following expression for the internal forces. p = p a + p b + p g (6.20) where p a = U a u = N [ c φ s φ 0 c φ s φ 0 ] T p b = U b u = [ Vs φ Vc φ M Vs φ Vc φ M 2 ] T p g = U g u = N [ 3s φ( θ + θ 2 )/ 3c φ ( θ + θ 2 )/ 4 θ θ 2 3s φ ( θ + θ 2 )/ 3c φ ( θ + θ 2 )/ 4 θ 2 θ ] T (6.2) 6 0
11 6. Notes and Bibliography Material Stiffness Matrix Carrying out the computations one obtains the following compact expression for the material stiffness: K M = T T K M0 T (6.22) where EA 0 0 EA EI 6EI 2 0 2EI 6EI 2 K M0 = 0 6EI 4EI 0 6EI 2EI, (6.23) EA 0 0 EA EI 2 6EI 0 2EI 2 6EI 0 6EI 2EI 0 6EI 4EI is the stiffness matrix of the linear beam-column element, and T is the transformation matrix c φ s φ s φ c φ T =, (6.24) c φ s φ s φ c φ which introduces the effect of finite rigid body motions Geometric Stiffness Matrix The expression for the geometric stiffness is a bit more complicated. It can be presented in a compact form as follows: K G = T T KG N T + KV G (6.25) where T is the transformation matrix (6.24), KG N is the well known geometric stiffness for a Hermitian beam element under axial force: KG N = N (6.26) and the remaining term introduces the effect of varying moments through the transverse shear force in C: sin 2φ cos 2φ 0 sin 2φ cos 2φ 0 cos 2φ sin 2φ 0 cos 2φ sin 2φ 0 KG V = V (6.27) sin 2φ cos 2φ 0 sin 2φ cos 2φ 0 cos 2φ sin 2φ 0 cos 2φ sin 2φ 0 in which sin 2φ = 2s φ c φ and cos 2φ = cφ 2 s2 φ. Note that K G is symmetric. 6
12 Chapter 6: THE CR FORMUATION: BE PANE BEAM Notes and Bibliography Most of the foregoing expressions were obtained in the writer s 966 Berkeley thesis [208], but not using the CR approach (which was not fully developed until the mid 980s). One thesis result omitted here is the more detailed expression of the contribution of the varying bending moment to the geometric stiffness. For most practical problems these refinements are not important. Similar results can be found in the 99 book by Crisfield [72, Ch. 7]. 6 2
13 Exercises Homework Exercises for Chapter 6 The Corotational Description: 2D C Beam EXERCISE 6. Complete the derivation of p for the 2-node C beam element and implement in Mathematica, using the same inputs as in Chapter 9 Addendum. (Implemented and posted on Web) EXERCISE 6.2 Complete the derivation of K for the 2-node C beam element and implement in Mathematica, using the same inputs as in Chapter 9 Addendum. (Implemented and posted on Web) EXERCISE 6.3 A plane 2-node C beam element has properties 0 = 6, E = 3000, A 0 = 2, I 0 = 2, N 0 = 5 in the initial state C 0 along X, with node at (0,0) and node 2 at ( 0, 0). The beam rotates by 45 about the origin so that at the current configuration C node stays at {0, 0} while node 2 moves to {( 0 + d)/ 2,( 0 + d)/ 2}, where d = 0 /000. The rotational freedoms at C are θ = θ 2 = 45 = π/4 radians. Compute p, K M and K G at the current configuration, and compare those quantities with those of the C 0 beam element presented in Chapter 9, using RBF for the latter. Note: AMathematica implementation of this C element has been posted on the Web as a Mathematica 4. Notebook PlaneBeamC.nb. The element checks out when moving about the reference configuration C 0. It gives excellent buckling values for the problem of Exercise 9.3. More tests are needed, however, for an arbitrary configuration to make sure the internal force vector and the tangent stiffness are consistent. EXERCISE 6.4 Confirm the previous statement by repeating the buckling calculations of Exercise 9.3 using the CR beam element provided in the Mathematica Notebook mentioned above (extract the material and stiffness matrices, ignore the rest). Compare the speed of convergence of the CR and T element for the cantilever buckling problem. 6 3
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 informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationThe TL Timoshenko Plane Beam Element
. 9 The TL Timoshenko Plane Beam Element 9 1 Chapter 9: THE TL TIMOSHENKO PLANE BEAM ELEMENT 9 2 TABLE OF CONTENTS Page 9.1. Introduction 9 3 9.2. Beam Models 9 3 9.2.1. Basic Concepts and Terminology...........
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 informationLaboratory 4 Topic: Buckling
Laboratory 4 Topic: Buckling Objectives: To record the load-deflection response of a clamped-clamped column. To identify, from the recorded response, the collapse load of the column. Introduction: Buckling
More informationLecture notes Models of Mechanics
Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /
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 informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationVariational Formulation of Plane Beam Element
. Variational Formulation of Plane Beam Element Chapter : VARIATIONAL FORMULATION OF PLANE BEAM ELEMENT TABLE OF CONTENTS Page.. Introduction 3.. What is a Beam? 3... Terminology................. 3...
More information9.1 Introduction to bifurcation of equilibrium and structural
Module 9 Stability and Buckling Readings: BC Ch 14 earning Objectives Understand the basic concept of structural instability and bifurcation of equilibrium. Derive the basic buckling load of beams subject
More informationENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1
ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at
More 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 informationCO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS
COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL
More informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More informationGyroscopic matrixes of the straight beams and the discs
Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation
More 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 informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More 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 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 informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
More informationCHAPTER -6- BENDING Part -1-
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 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 informationChapter 2: Deflections of Structures
Chapter 2: Deflections of Structures Fig. 4.1. (Fig. 2.1.) ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 1 (2.1) (4.1) (2.2) Fig.4.2 Fig.2.2 ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 2
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationGeneral 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 informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn
ENGN234 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the Euler-Bernuolli
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More 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 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 informationComb resonator design (2)
Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
More informationVariational Formulation of Plane Beam Element
13 Variational Formulation of Plane Beam Element IFEM Ch 13 Slide 1 Beams Resist Primarily Transverse Loads IFEM Ch 13 Slide 2 Transverse Loads are Transported to Supports by Flexural Action Neutral surface
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More information4.5 The framework element stiffness matrix
45 The framework element stiffness matri Consider a 1 degree-of-freedom element that is straight prismatic and symmetric about both principal cross-sectional aes For such a section the shear center coincides
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationA study of the critical condition of a battened column and a frame by classical methods
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School 003 A study of the critical condition of a battened column and a frame by classical methods Jamal A.H Bekdache
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
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 beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationMITOCW MITRES2_002S10linear_lec07_300k-mp4
MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More information7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment
7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment à It is more difficult to obtain an exact solution to this problem since the presence of the shear force means that
More informationDue Monday, September 14 th, 12:00 midnight
Due Monday, September 14 th, 1: midnight This homework is considering the analysis of plane and space (3D) trusses as discussed in class. A list of MatLab programs that were discussed in class is provided
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationWorkshop 8. Lateral Buckling
Workshop 8 Lateral Buckling cross section A transversely loaded member that is bent about its major axis may buckle sideways if its compression flange is not laterally supported. The reason buckling occurs
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationComb Resonator Design (2)
Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationChapter 3. Load and Stress Analysis. Lecture Slides
Lecture Slides Chapter 3 Load and Stress Analysis 2015 by McGraw Hill Education. This is proprietary material solely for authorized instructor use. Not authorized for sale or distribution in any manner.
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationDue Tuesday, September 21 st, 12:00 midnight
Due Tuesday, September 21 st, 12:00 midnight The first problem discusses a plane truss with inclined supports. You will need to modify the MatLab software from homework 1. The next 4 problems consider
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 informationParametric study on the transverse and longitudinal moments of trough type folded plate roofs using ANSYS
American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-4 pp-22-28 www.ajer.org Research Paper Open Access Parametric study on the transverse and longitudinal moments
More informationNonconservative Loading: Overview
35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed 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 informationTheories of Straight Beams
EVPM3ed02 2016/6/10 7:20 page 71 #25 This is a part of the revised chapter in the new edition of the tetbook Energy Principles and Variational Methods in pplied Mechanics, which will appear in 2017. These
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 informationDecember 10, PROBLEM NO points max.
PROBLEM NO. 1 25 points max. PROBLEM NO. 2 25 points max. B 3A A C D A H k P L 2L Given: Consider the structure above that is made up of rod segments BC and DH, a spring of stiffness k and rigid connectors
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 informationAE3160 Experimental Fluid and Solid Mechanics
AE3160 Experimental Fluid and Solid Mechanics Cantilever Beam Bending Claudio Di Leo 1 Learning Objectives 1. On Structural Mechanics: a) Mechanics of Slender Beams b) Strain Transformation Theory c) Principal
More informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationThermomechanical Effects
3 hermomechanical Effects 3 Chapter 3: HERMOMECHANICAL EFFECS ABLE OF CONENS Page 3. Introduction..................... 3 3 3.2 hermomechanical Behavior............... 3 3 3.2. hermomechanical Stiffness
More informationNumerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies
Numerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies F. D. Sorokin 1, Zhou Su 2 Bauman Moscow State Technical University, Moscow,
More information14. *14.8 CASTIGLIANO S THEOREM
*14.8 CASTIGLIANO S THEOREM Consider a body of arbitrary shape subjected to a series of n forces P 1, P 2, P n. Since external work done by forces is equal to internal strain energy stored in body, by
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 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 informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More 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 arch-beam with using matlab.
More informationVariational Formulation of Plane Beam Element
Variational Formulation of Plane Beam Element Chapter : VARIATIONAL FORMULATION OF PLANE BEAM ELEMENT TABLE OF CONTENTS Page. Introduction..................... 3. What is a Beam?................... 3..
More informationLecture Slides. Chapter 4. Deflection and Stiffness. The McGraw-Hill Companies 2012
Lecture Slides Chapter 4 Deflection and Stiffness The McGraw-Hill Companies 2012 Chapter Outline Force vs Deflection Elasticity property of a material that enables it to regain its original configuration
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 informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationMoment Area Method. 1) Read
Moment Area Method Lesson Objectives: 1) Identify the formulation and sign conventions associated with the Moment Area method. 2) Derive the Moment Area method theorems using mechanics and mathematics.
More informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationOn Nonlinear Buckling and Collapse Analysis using Riks Method
Visit the SIMULIA Resource Center for more customer examples. On Nonlinear Buckling and Collapse Analysis using Riks Method Mingxin Zhao, Ph.D. UOP, A Honeywell Company, 50 East Algonquin Road, Des Plaines,
More informationModule 3. Analysis of Statically Indeterminate Structures by the Displacement Method
odule 3 Analysis of Statically Indeterminate Structures by the Displacement ethod Lesson 14 The Slope-Deflection ethod: An Introduction Introduction As pointed out earlier, there are two distinct methods
More informationε t increases from the compressioncontrolled Figure 9.15: Adjusted interaction diagram
CHAPTER NINE COLUMNS 4 b. The modified axial strength in compression is reduced to account for accidental eccentricity. The magnitude of axial force evaluated in step (a) is multiplied by 0.80 in case
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More informationME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft.
ME 323 - Final Exam Name December 15, 2015 Instructor (circle) PROEM NO. 4 Part A (2 points max.) Krousgrill 11:30AM-12:20PM Ghosh 2:30-3:20PM Gonzalez 12:30-1:20PM Zhao 4:30-5:20PM M (x) y 20 kip ft 0.2
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 10 Columns
EMA 370 Mechanics & Materials Science (Mechanics of Materials) Chapter 10 Columns Columns Introduction Columns are vertical prismatic members subjected to compressive forces Goals: 1. Study the stability
More informationPURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.
BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally
More informationExercise: concepts from chapter 5
Reading: Fundamentals of Structural Geology, Ch 5 1) Study the oöids depicted in Figure 1a and 1b. Figure 1a Figure 1b Figure 1. Nearly undeformed (1a) and significantly deformed (1b) oöids with spherulitic
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More informationFree vibration and stability of axially functionally graded tapered Euler-Bernoulli beams
Shock and Vibration 18 211 683 696 683 DOI 1.3233/SAV-21-589 IOS Press Free vibration and stability of axially functionally graded tapered Euler-Bernoulli beams Ahmad Shahba a,b, Reza Attarnejad a,b, and
More information7.4 The Elementary Beam Theory
7.4 The Elementary Beam Theory In this section, problems involving long and slender beams are addressed. s with pressure vessels, the geometry of the beam, and the specific type of loading which will be
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More 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 informationMulti Linear Elastic and Plastic Link in SAP2000
26/01/2016 Marco Donà Multi Linear Elastic and Plastic Link in SAP2000 1 General principles Link object connects two joints, i and j, separated by length L, such that specialized structural behaviour may
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 informationIndeterminate Analysis Force Method 1
Indeterminate Analysis Force Method 1 The force (flexibility) method expresses the relationships between displacements and forces that exist in a structure. Primary objective of the force method is to
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
More informationModule 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy
The Lecture Contains: Torsion of Rods Torsional Energy This lecture is adopted from the following book 1. Theory of Elasticity, 3 rd edition by Landau and Lifshitz. Course of Theoretical Physics, vol-7
More information