A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
|
|
- Damon Holt
- 5 years ago
- Views:
Transcription
1 European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES Gordan Jelenić 1, Maja Gaćeša 1, and Miran Saje 2 1 University of Rijeka, Faculty of Civil Engineering Radmile Matejčić 3, 51 Rijeka, Croatia {gordan.jelenic,maja.gacesa@gradri.hr 2 University of Ljubljana, Faculty of Civil and Geodetic Engineering Jamova cesta 2, 11 Ljubljana, Slovenia msaje@fgg.uni-lj.si Keywords: 3D beams, non-linear analysis, configuration-dependent interpolation, fixed-pole approach Abstract. Fixed-pole approach in the non-linear analysis of beam structures offers a different and promising framework for development of new element types with particular potential in the area of conservative time integrators. The finite elements formulated using this approach, however, utilise non-standard degrees of freedom which makes them difficult to combine with standard finite elements in existing codes. A modified approach which combines the benefits of the fixed-pole approach with the use of standard kinematic unknowns from the existing (movingpole) approaches is presented in this work along with a family of solution procedures for static analysis of geometrically non-linear spatial beam structures. Such a hybrid formulation opens a number of questions which are specifically addressed.
2 1 INTRODUCTION 1.1 Notation L, x initial length of the beam, co-ordinate of a point along the centroid axis of the beam E i orthonormal base vectors of the material co-ordinate system (i 1, 2, 3) I, 3 3 unit and null matrix, respectively r, Λ position vector of a point along the beam, rotation matrix of a triad rigidly attached to the cross-section u displacement vector N, M vectors of stress and stress-couple resultants in material co-ordinate system n, m vectors of stress and stress-couple resultants in spatial co-ordinate system n, m fixed-pole vectors of stress and stress-couple resultants in spatial co-ordinate system Γ, K translational and rotational strain measures in material co-ordinate system γ, κ translational and rotational strain measures in spatial co-ordinate system γ, κ fixed-pole translational and rotational strain measures in spatial co-ordinate system n e, m e vectors of applied distributed forces and torques n e, m e fixed-pole vectors of applied distributed forces and torques F, F L,M, M L applied concentrated forces at ends of the beam F, F L, M, M L fixed pole applied concentrated forces at ends of the beam C N, C M translational and rotational constitutive matrices E, G Young s (elastic) modulus, shear modulus A 1, A 2, A 3 cross-sectional and shear areas J, I 2, I 3 torsional and cross-sectional moments of inertia A diag [A 1, A 2, A 3 ] J ρ diag [J, I 2, I 3 ] ρ material density k, π vectors of specific translational momentum and angular momentum δr, δϑ variation of the position vector, spin variable N number of interpolation points (nodes) of a finite element I i (x) i-th Lagrangian interpolation polynomial of order N 1 G (if not in C N ) weak form g i, q i i, qi m, q i e residual and nodal vectors of internal, inertial and external forces at node i K ij stiffness matrix relating the trial functions at node j to the residual at node i ( ), ( ) ( ), ( ) t ( ) i, ( ) j δ( ), ( ) derivatives of ( ) with respect to x, and time t cross-product operator ( ), transposed quantity ( ) ( ) at nodes i and j variation of ( ), incremental change of ( ) The geometrically exact beam theory provided by Simo in 1985 [1] makes one of the milestones in the development of the non-linear finite-element method in introducing a non-linear manifold as the solution space and the related differential geometry tools in order to implement a robust solution procedure [2]. This approach, however, turns out to be incompatible with algorithmic preservation of strain invariance with respect to a rigid-body motion, or simultaneous conservation of energy and the momentum vectors during a free motion of an unloaded beam. Borri and Bottasso [3, 4] propose the fixed-pole approach in which they introduce new stresscouple resultant and angular momentum - defined with respect to a unique point for the whole structure - the fixed pole. Different implementations of this general concept result in the algo- 2
3 rithms which naturally inherit the strain-invariance of the underlying formulation with respect to a rigid body motion [3] or are capable of simultaneous conservation of both energy and momentum vectors [4]. The main problem with the fixed-pole finite element implementation is that it uses nonstandard problem unknowns, and therefore is not easily combinable with existing finite element codes. For this reason, we derive and present a formulation which uses standard kinematic unknowns but is inherently based on the fixed-pole approach, along with with the algorithmic benefits of this approach. 2 KINEMATIC AND CONSTITUTIVE EQUATIONS Kinematic equations follow from the Reissner-Simo theory [1]. The material translational and rotational strain measures are given as follows: Γ Λ t γ Λ t r E 1 ˆK Λ t κλ Λ t Λ, (1) while the spatial translational and angular velocity vectors, v and w, are defined as v ṙ and ŵ Λ t Λ. Varying (1) leads to δγ δr + r δϑ δκ δϑ (2). Let us consider a linear elastic material with the material cross-sectional stress and stresscouple resultants defined as [ ] N CN Γ, (3) M C M K with EA GI t C N GA 1, C M EI 1. (4) GA 2 EI 2 3 THE STANDARD APPROACH The standard approach [2] given by Simo and Vu-Quoc follows naturally from the principle of virtual work: V i + V m V e, (5) where V i, V m and V e are the virtual works of internal, inertial and external forces, respectively, which are defined as follows: n L [ ] I V i δγ t δκ t dx δr m t δϑ t d dx n r I d dx, (6) m dx k V m δr t δϑ t dx, (7) π V e δr t δϑ t ne F dx + δr m t δϑ t F + δr e M t L L δϑt L. (8) M L 3
4 A certain irregularity in the integrals should be noted while the virtual quantities simply dot-multiply the specific vectors of inertial and external forces, an additional operator occurs between the virtual quantities and the stress resultants: [ ] I d dx r I d dx. (9) After interpolating virtual quantities δr, δϑ using Lagrangian polynomials I i (x): δr I i (x)δr i, δϑ i1 I i (x)δϑ i, (1) i1 we have the approximate weak form: G h t δr i δϑ t i g i, (11) i1 from which, because of the arbitrariness of δr i and δϑ i, follows the vector of nodal dynamic residuals: g i q i i + q i m q i e, (12) with q i i, q i m and q i e as the nodal vectors of internal, inertial and external forces, respectively, given as follows: [ I q i i ] I n i I i r I i dx, (13) I m k q i m I i dx, (14) π q i e I i ne F dx + δ1 i F + δn i L. (15) m e 4 THE FIXED-POLE APPROACH In [3, 4] Borri and Bottasso thoroughly investigated the idea of replacing the stress-couple resultant m and the specific angular momentum π, which are defined with respect to the beam reference axis at a cross-section, with another stress-couple resultant m and specific angular momentum π, which are to be defined with respect to a unique point for the whole structure the fixed pole naturally taken to be the origin of the spatial frame, ie: M M L m r n + m, (16) π r k + π. (17) We also note that n n and k k. In order to obtain a relationship between standard and fixed-pole strain measures, we request the strain energy density to remain invariant to the change of virtual quantities: 1 2 ( γ n + κ m) 1 (γ n + κ m) (18) 2 4
5 By substituting the relationship (16) between standard and fixed pole stress resultants, we obtain the relationship between standard and fixed-pole spatial strain measures: γ γ + r κ, (19) κ κ. (2) It is also useful to note the relationship between the time derivatives of standard and fixedpole specific momenta: k k, π ṙ {{ k + r k + π, (21) ṙ Aρṙ and considering virtual work of the inertial forces (7) to be an objective quantity gives us the fixed-pole kinematics which is virtual-work conjugate to the time-derivatives of the specific fixed-pole momenta as δ r δr + r δϑ, δ ϑ (22) δϑ, from where their respective derivatives follow as δ r δr + r δϑ + r δϑ, δ ϑ δϑ. (23) Analogously to (16), the fixed-pole distributed moment is introduced: m e m e + r n e, (24) noting that n e n e. We also introduce the relationship between standard and fixed-pole concentrated moments at ends of the beam: M M + r F, M L M L + r L F L, (25) noting that F F and F L F L. 4.1 Original formulation We rewrite the expression for the virtual work of internal forces (6) and substitute (23) and (16) into it: n L V i (δr + r δϑ) t δϑ t dx ((δr + r δϑ) n + δϑ m) dx m ((δ r r δϑ r δϑ + r δϑ) n + δϑ ( m r n)) dx ( δ r n + δ ϑ ) m dx, (26) 5
6 and obtain the virtual work of internal forces using the fixed-pole virtual quantities and stress and stress-couple resultants: ( ) { d n m V i δ r dx t δ ϑ t dx. (27) Note that there is no irregularity (9) in this term! In a like manner, we rewrite the expression for virtual work of the inertial forces (7) and substitute (22) and (21) into it: V m ( (δ r r δϑ) k + δϑ π r k ) (δ r r δϑ) t δϑ t to get the fixed-pole virtual work of inertial forces as V m δ r t δ ϑ t dx { k π r k dx ( δ r k ) + δ ϑ π dx, (28) k dx. (29) π Finally, the virtual work of external forces (8) may be treated in the same way, so that after substituting (25) and (22) into it: { V e (δ r r δϑ) t δϑ t n e dx+ m e r n e { + (δ r r δϑ ) t δϑ t F + M r F { + (δ r L r L δϑ L ) t δϑ t L F L, (3) M L r L F L we arrive at the following fixed-pole virtual work of external forces: V e δ r t δ ϑ t ne F dx + δ r m t δ ϑ t F + δ r e M t L δ ϑ t L L. M L Deciding to interpolate the fixed-pole virtual quantities δ r, δ ϑ using Lagrangian polynomials I i (x), δ r I i (x)δ r i, δ ϑ i1 we obtain the approximate weak form: I i (x)δ ϑ i, (31) i1 G h δ r t i δ ϑ t i ḡi, (32) i1 from which follows, because of the arbitrariness of δ r i and δ ϑ i, ḡ i q i i + q i m q i e, (33) 6
7 with ḡ i as the fixed-pole vector of nodal dynamic residuals and q i i, q i m and q i e as the fixed-pole nodal vectors of internal, inertial and external forces, respectively, given as q i i I i dx, (34) n m k q i m I i dx, (35) π q i e 4.2 Modified formulation I i { ne m e dx + δ i 1 F F + δn M i L. (36) M L As shown in (32), the fixed-pole nodal dynamic residual ḡ i is virtual-work conjugate to the fixed-pole nodal virtual displacements and rotations which are non-standard and prevent the elements from being combined with standard elements in a finite-element mesh. This problem, however, is easily overcome by noting that from (22) the nodal fixed-pole virtual quantities are related to the standard nodal virtual quantities via: [ ] δ ri δri δ ϑ + r i δϑ i I ri δri, (37) i δϑ i I δϑ i so that the fixed-pole weak form (32) becomes with G h [ ] I δr t i δϑ t i ḡ r i i I i1 δr t i δϑ t i g i, (38) i1 g i q i i + q i m q i e, (39) and q i i q i m q i e [ ] I n I i r dx, (4) r i I m [ ] I k I i r dx, (41) r i I π [ ] I I i ne F r dx + δ i F 1 + δ i L N. (42) r i I m e M M L After substituting (4), (41) and (42) in g i and then back to (32), it can be shown that, using a non-linear interpolation of the type: δr δϑ I i [δr i + (r i r) δϑ i ] i1 (43) I i δϑ i i1 is equivalent to the standard (Lagrangian) interpolation of the fixed-pole virtual quantities (31). 7
8 5 SOLUTION PROCEDURE The non-linear equation (39) may be solved using the Newton-Raphson solution procedure. In order to concentrate on the issues that arise when attempting to utilise a fixed-pole approach in the hybrid manner presented here, in which the standard displacement and rotation quantities have been kept as the problem unknowns, we proceed by considering only the static case. When linearising (39), account has to be taken of the fact that Λ ϑλ and Γ Λ t ( r + r ϑ) as well as [1] ( ( ) ) t ( ( ) ) K new Λ t newλ new exp ϑ Λ old exp ϑ Λ old ( ( )) t ( ( )) ( ( )) t ( ) Λ t old exp ϑ exp ϑ Λold + Λ t old exp ϑ exp ϑ Λ old {{ I K + K old, (44) where and ( ) exp ϑ I + sin ϑ ϑ ϑ 1 cos ϑ + ϑ 2, (45) ϑ 2 with H ϑ sin ϑ ϑ ϑ t + ϑ 3 K Λ t oldh ϑ, (46) sin ϑ ϑ The strain measure updates for K and Γ then follow as 1 cos ϑ I ϑ. (47) ϑ 2 K new K old + K, (48) Γ new : Λ t newr new E 1. (49) After expanding equation (39) into a Taylor series around a known configuration and omitting higher-order terms we have where g i + g i, (5) g i q i i q i e. (51) After the linearisation of the modified nodal internal force vector (4) we have q i i [ ] I i ri r L [ ] [ ] I n r dx + I n i dx+ ϑ i ϑ r r i I m ϑ [ ] [ ] [ ] [ ] I Λ + I i CN Λ t r + r ϑ r r i I Λ C M Λ t ϑ dx. (52) 8
9 Linearizing the modified nodal external force vector (42) we have: [ ] q i e I i ri r dx. (53) n e ϑ i ϑ An interesting phenomenon occurs when linearising the external force vector, which is a direct consequence of the use of the fixed-pole virtual quantities and their conversion back to standard virtual quantities: even though the external loading is configuration-independent, the modified nodal external loading vector (42) is not. As a result of this, the stiffness matrix will have an additional part as a result of the linearisation of q i e. In order to complete the solution procedure, the unknown functions and/or their iterative changes must be interpolated in some way. Within this paper we consider three different interpolation options, which arise as a consequence of the fact that the virtual position vector has been interpolated in a non-linear manner, dependent on the actual position vector not only at the nodal points, but also at x (i.e. the integration point). 5.1 Interpolation option 1 Within this option r and ϑ are interpolated in the same way as δr and δϑ (43): r ϑ I j [ r j + (r j r) ϑ j ] j1 I j ϑ j, while the unknown displacement functions are interpolated as: j1 r (54) I k r k. (55) k1 It should be noted that this interpolation is in contradiction with the interpolation for r in (54), which obviously does not follow as a linearisation of (55). The above contradiction necessarily results in loss of quadratic convergence of the Newton-Raphson solution process, and has been introduced as a simplest means of providing r(x) in (54). After introducing (54) into (52) and (53) we have: q i i K ij rj int, (56) ϑ j q i e K ij rj ext, (57) ϑ j where [ ] int I i L [ ] n (δ ij I j ) n I j n( r j dx + I i I j r) ( r dx r i ) n m [ ] [ ] [ ] [ ] [ ] I Λ + I i I j CN Λ t I rj r r r i I Λ C M Λ t dx, I (58) K ij 9
10 and and finally K ij ext [ ] I i (I j δ ij ) n e I j n e ( r j dx, (59) r) K ij K ij int Kij ext. (6) 5.2 Interpolation option 2 Within this option, r is interpolated using (55). Interpolation of the displacement increments follows consistently by linearising r and ϑ: r I j r j, ϑ j1 I j ϑ j. (61) After introducing (61) into (52) and (53) we have q i i K ij rj int, (62) ϑ j q i e K ij rj ext, (63) ϑ j where and K ij int + [ ] I i (δ n ij I j ) dx + ] [ Λ [ I I i r r i I K ij ext j1 [ ] [ ] I n I i I j dx+ r r i I m ] [ ] [ ] [ CN Λ t I j ] I I j r Λ C M Λ t I j dx, I I i (I j δ ij ) [ ] dx, (64) n e and finally K ij K ij int Kij ext. (65) In this interpolation option, different interpolations have been used for the test functions and the trial functions, which is bound to make the tangent stiffness matrix strongly non-symmetric. 5.3 Interpolation option 3 Within this option r is not interpolated, but only updated as follows: r new (x) r old (x) + r(x), (66) where r old (x) is the last known value for r(x), not necessarily associated with an equilibrium state and r(x) following interpolation of the Newton-Raphson changes given in (54). This is a consistent choice which yields a tangent stiffness matrix. In this case the values of r at both integration and nodal points must be saved. The earlier results (56) (6) can still be used, but now noting that: r r g r(x g ). (67) 1
11 6 CONCLUDING REMARKS The fixed-pole approach, first introduced by Borri and Bottasso [3] has been analysed and discussed. The approach uses a new set of stress-couple resultants and specific angular momentum, defined with respect to the fixed-pole rather than the reference line at the cross-section, which in turn implies that r(x), the standard position vector, ceases to be the problem unknown. The presented modified fixed-pole approach makes it easy to combine the standard choice of the problem unknowns with the benefits of the original fixed-pole approach. The solution procedure for a spatial beam under static loading has been presented and the complexities of the modified approach discussed. In particular, the approach is based on a non-linear interpolation of the virtual displacements, which is dependent on the position vector field itself, which in turn opens a number of combinations in interpolating the position vector field and its iterative change. The solution procedure for a spatial beam under dynamic loading will be derived and the numerical examples shown in future works. ACKNOWLEDGEMENT Research resulting with this paper was made within the scientific project No : Improved accuracy in non-linear beam elements with finite 3D rotations financially supported by the Ministry of Science, Education and Sports of the Republic of Croatia. The second author additionally acknowledges the support of the Croatian Science Foundation which funded the bilateral project No 3.1/129: Preservation of mechanical constants in numerical time integration of non-linear beam equations of motion. REFERENCES [1] J.C. Simo. A finite strain beam formulation. The three-dimensional dynamic problem. Part I. Computer Methods in Applied Mechanics and Engineering, 49:55 7, [2] J.C. Simo and L. Vu-Quoc. On the dynamics in space of rods undergoing large motions A geometrically exact approach. Computer Methods in Applied Mechanics and Engineering, 66: , [3] M. Borri and C. Bottasso. An intrinsic beam model based on a helicoidal approximation Part I: Formulation. International Journal for Numerical Methods in Engineering, 37(13): , [4] C. Bottasso and M. Borri. Integrating finite rotations. Computer Methods in Applied Mechanics and Engineering, 164(3-4):37 331,
University of Rijeka, Faculty of Civil Engineering, R. Matejčić 3, Rijeka,
Proceedings of the XII Finnish Mechanics Days R. Kouhia, J. Mäkinen, S. Pajunen and T. Saksala (Eds.) c The Authors, 215. Open access under CC BY-SA 4. license. Finite elements on non-linear manifolds
More informationFinite Elements in Analysis and Design
Finite Elements in Analysis and Design 99 (15) 39 48 Contents lists available at ScienceDirect Finite Elements in Analysis and Design journal homepage: www.elsevier.com/locate/finel Modified fixed-pole
More informationModeling of 3D Beams
Modeling of 3D Beams Friday, July 11 th, (1 Hour) Goran Golo Goran Golo Preliminaries Goran Golo Preliminaries Limitations of liner beam theory Large deviations Equilibrium is not a straight line 3 Preliminaries
More informationMODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS
PROCEEDINGS Proceedings of the 25 th UKACM Conference on Computational Mechanics 12-13 April 217, University of Birmingham Birmingham, United Kingdom MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN
More informationGeometrically exact beam dynamics, with and without rotational degree of freedom
ICCM2014 28-30 th July, Cambridge, England Geometrically exact beam dynamics, with and without rotational degree of freedom *Tien Long Nguyen¹, Carlo Sansour 2, and Mohammed Hjiaj 1 1 Department of Civil
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 informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationNUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS
NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering
More informationJamova cesta 2 SI 1000 Ljubljana, Slovenia Jamova cesta Ljubljana, Slovenija
Univerza v Ljubljani Fakulteta za gradbeništvo in geodezijo University of Ljubljana Faculty of Civil and Geodetic Engineering Jamova cesta 2 1 Ljubljana, Slovenija http://www3.fgg.uni-lj.si/ Jamova cesta
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 informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
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 informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
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 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 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 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 informationINTERPOLATION OF ROTATIONAL VARIABLES IN NONLINEAR DYNAMICS OF 3D BEAMS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 43, 1193 1222 (1998) INTERPOLATION OF ROTATIONAL VARIABLES IN NONLINEAR DYNAMICS OF 3D BEAMS G. JELENI C AND M. A.
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 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 informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More 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 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 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 informationLarge deflection analysis of planar solids based on the Finite Particle Method
yuying@uiuc.edu 10 th US National Congress on Computational Mechanics Large deflection analysis of planar solids based on the Finite Particle Method 1, 2 Presenter: Ying Yu Advisors: Prof. Glaucio H. Paulino
More informationUNIT IV FLEXIBILTY AND STIFFNESS METHOD
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : SA-II (13A01505) Year & Sem: III-B.Tech & I-Sem Course & Branch: B.Tech
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part 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 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 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 informationChapter 4 Statics and dynamics of rigid bodies
Chapter 4 Statics and dynamics of rigid bodies Bachelor Program in AUTOMATION ENGINEERING Prof. Rong-yong Zhao (zhaorongyong@tongji.edu.cn) First Semester,2014-2015 Content of chapter 4 4.1 Static equilibrium
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
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 informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
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 informationNonlinear analysis in ADINA Structures
Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent
More informationResearch Article A Hybrid Interpolation Method for Geometric Nonlinear Spatial Beam Elements with Explicit Nodal Force
Mathematical Problems in Engineering Volume 216, Article ID 898676, 16 pages http://dx.doi.org/1.1155/216/898676 Research Article A Hybrid Interpolation Method for Geometric Nonlinear Spatial Beam Elements
More informationStructural Dynamics Prof. P. Banerji Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 1 Introduction
Structural Dynamics Prof. P. Banerji Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 1 Introduction Hello, I am Pradipta Banerji from the department of civil engineering,
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationarxiv: v1 [cs.ce] 1 Sep 2016
Geometrically Exact Finite Element Formulations for Curved Slender Beams: Kirchhoff-Love Theory vs. Simo-Reissner Theory Christoph Meier a,, Wolfgang A. Wall b, Alexander Popp b a Mechanosynthesis Group,
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 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 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 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 informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationCOROTATIONAL METHOD AND DIFFERENTIAL GEOMETRY
www.arpapress.com/volumes/vol15issue1/ijrras_15_1_07.pdf COROAIONAL MEHOD AND DIFFERENIAL GEOMERY Mohammed Khalil Ferradi SEEC-PI, 42/52 Quai de la Rapée 75012, Paris, France E-mail: mohammed-khalil.ferradi@tpi.setec.fr
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More 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 informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
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 informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
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 informationGeneral method for simulating laboratory tests with constitutive models for geomechanics
General method for simulating laboratory tests with constitutive models for geomechanics Tomáš Janda 1 and David Mašín 2 1 Czech Technical University in Prague, Faculty of Civil Engineering, Czech Republic
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 informationOn characterising fracture resistance in mode-i delamination
9 th International Congress of Croatian Society of Mechanics 18-22 September 2018 Split, Croatia On characterising fracture resistance in mode-i delamination Leo ŠKEC *, Giulio ALFANO +, Gordan JELENIĆ
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 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 informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
More informationUniversity of Groningen
University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published
More informationNUMERICAL SIMULATION OF STEEL CATENARY RISER
SIMMEC/EMMCOMP 214 XI Simpósio de Mecânica Computacional II Encontro Mineiro de Modelagem Computacional Juiz De Fora, MG, 28-3 de Maio De 214 NUMERICAL SIMULATION OF STEEL CATENARY RISER Marcus V. S. Casagrande
More informationMITOCW MITRES2_002S10nonlinear_lec05_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec05_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 informationMODAL DERIVATIVES BASED REDUCTION METHOD FOR FINITE DEFLECTIONS IN FLOATING FRAME
Modal derivatives based reduction method for finite deflections in floating frame 11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th
More informationPhysics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics
Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES
ANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES Jerome Merel formerly at Hutchinson Corporate Research Center Israël Wander Apex Technologies Pierangelo Masarati, Marco Morandini Dipartimento
More informationThe CR Formulation: BE Plane Beam
6 The CR Formulation: BE Plane Beam 6 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction..................... 6 4 6.2 CR Beam Kinematics................. 6 4 6.2. Coordinate
More informationDynamic Model of a Badminton Stroke
ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815
More informationGEO E1050 Finite Element Method Autumn Lecture. 9. Nonlinear Finite Element Method & Summary
GEO E1050 Finite Element Method Autumn 2016 Lecture. 9. Nonlinear Finite Element Method & Summary To learn today The lecture should give you overview of how non-linear problems in Finite Element Method
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationNeuro -Finite Element Static Analysis of Structures by Assembling Elemental Neuro -Modelers
Neuro -Finite Element Static Analysis of Structures by Assembling Elemental Neuro -Modelers Abdolreza Joghataie Associate Prof., Civil Engineering Department, Sharif University of Technology, Tehran, Iran.
More informationLecture 4: PRELIMINARY CONCEPTS OF STRUCTURAL ANALYSIS. Introduction
Introduction In this class we will focus on the structural analysis of framed structures. We will learn about the flexibility method first, and then learn how to use the primary analytical tools associated
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 informationDynamic analysis of 3D beams with joints in presence of large rotations
Comput. Methods Appl. Mech. Engrg. 19 1) 4195±43 www.elsevier.com/locate/cma Dynamic analysis of 3D beams with joints in presence of large rotations G. Jelenic, M.A. Cris eld * Department of Aeronautics,
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
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 informationSensitivity and Reliability Analysis of Nonlinear Frame Structures
Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,
More informationVibration analysis of circular arch element using curvature
Shock and Vibration 15 (28) 481 492 481 IOS Press Vibration analysis of circular arch element using curvature H. Saffari a,. Tabatabaei a, and S.H. Mansouri b a Civil Engineering Department, University
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 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 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 informationMATHEMATICAL ANALYSIS OF FLY FISHING ROD STATIC AND DYNAMIC RESPONSE
MATHEMATICAL ANALYSIS OF FLY FISHING ROD STATIC AND DYNAMIC RESPONSE DER-CHEN CHANG, GANG WANG, AND NORMAN M. WERELEY We develop two mathematical models to study the fly fishing rod static and dynamic
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_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 informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
More informationTheory & Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory & Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 5 Torsional Vibrations Lecture - 4 Transfer Matrix Approach
More informationBiaxial Analysis of General Shaped Base Plates
Biaxial Analysis of General Shaped Base Plates R. GONZALO ORELLANA 1 Summary: A linear model is used for the contact stresses calculation between a steel base plate and a concrete foundation. It is also
More informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationIntroduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0
Introduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0 Randy H. Shih Oregon Institute of Technology SDC PUBLICATIONS Schroff Development Corporation www.schroff.com Better Textbooks.
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 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 informationInvestigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
Print ISSN: 2322-2093; Online ISSN: 2423-6691 DOI: 10.7508/ceij.2016.02.011 Technical Note Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More information4 NON-LINEAR ANALYSIS
4 NON-INEAR ANAYSIS arge displacement elasticity theory, principle of virtual work arge displacement FEA with solid, thin slab, and bar models Virtual work density of internal forces revisited 4-1 SOURCES
More information