Some Remarks on the Reissner Mindlin Plate Model
|
|
- Hester Beasley
- 5 years ago
- Views:
Transcription
1 Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10,
2 Outline The 3D Problem and its Modeling Full Discretization Playing Around with the Timoshenko Beam Conclusions 2
3 The 3D Problem and its Modeling 3D linearly elastic plates A Hellinger Reissner Principle Reissner Mindlin Model Derivation Convergence Estimate Full Discretization Playing Around with the Timoshenko Beam Conclusions 3
4 The Domain Ω P ε x = (x, x 3 ) 2ε Middle surface: Ω R 2, 3D Domain: P ε = Ω ( ε, ε), Lateral Boundary: P ε L = Ω ( ε, ε), Top and bottom boundaries: P ε ± = Ω { ε, ε}. 4
5 Bending of Linearly Elastic Plate Find u ε : P ε R 3 and σ ε : P ε R 3 3 sym such that div σ ε = 0, σ ε = Ce(u ε ), in P ε, σ ε n = g ε on P ε ±, u ε = 0 on P ε L, where g ε = (0, 0, g 3 ) : P ε ± R 3 is the traction load. Also, e(u ε ) = 1 2 ( uε + T u ε ) and Cτ = 2ντ + λ tr(τ)δ, where µ and λ are the Lamé coeff., and δ is the 3 3 ident. matrix. Assume g 3 (x, ε) = g 3 (x, ε). 5
6 Dimension Reduction The goal of dimension reduction is to pose a system of equations in the middle surface Ω that is easy to solve and which solution approximates the original 3D solution. The most popular models are variants of Reissner Mindlin Kirchhoff Love (also known as biharmonic model) formal higher order models 6
7 Notation: u = u u 3 R2 R, σ = σ σ σ t σ 33 R2 2 R 2 R 2t R, 7
8 Hellinger Reissner Principle Let S g (P ε ) = { τ H(div, P ε ) : τn = g on P± ε }. Hellinger-Reissner Principle: (σ ε, u ε ) is critical point of L(τ, v) = 1 C τ : τ dx + div τ v dx 2 P ε P ε on S g (P ε ) L 2 (P ε ). So σ ε S g (P ε ) and u ε L 2 (P ε ) satisfy C σ ε : τ dx + P ε u ε div τ dx = 0 P ε for all τ S 0 (P ε ), div σ ε v dx = 0 P ε for all v L 2 (P ε ). 8
9 We now derive a model of Reissner Mindlin type by looking for critical points of L within subspaces of S g (P ε ) L 2 (P ε ). Consider subspaces composed of functions which are polynomial in the transverse direction, with the following degrees: 1, Then our solution will be of the form θ u R = (x )x 3, σ R = ω(x ) + ω 2 (x )r(x 3 ) σ(x)x 3 σ(x )q(x 3 ), g 3 x 3 /ε + σ 33 s(x 3 ) where r(z) = 3(z 2 /ε 2 1/5)/2, q(z) = 3(1 z 2 /ε 2 )/2, and s(z) = 5(z/ε z 3 /ε 3 )/4. 9
10 Consider the eqtn P ε div σ v dx = 0, where v is a test function. Choose v = (ψ (x )x 3, 0), with ψ L 2 (Ω) arbitrary, to obtain ε2 3 Ω div σ ψ dx + 2ε Ω σ ψ dx = 0 for all ψ L 2 (Ω). Choosing v = (0, 0, w(x )), with w L 2 (Ω) arbitrary, it follows that ε Ω div σ w dx = Ω g 3 w dx, for all w L 2 (Ω) Finally, with v = (0, 0, w 2 (x )r(x 3 )), and w 2 L 2 (Ω) arbitrary, Ω σ 33 w 2 dx = 2 5 Ω g 3 w 2 dx, for all w 2 L 2 (Ω) 10
11 Hence ε2 3 div σ ψ dx + 2ε σ ψ dx = 0 for all ψ L2 (Ω) Ω Ω ε div σ w dx = g 3 w dx, for all w L 2 (Ω) Ω Ω σ 33 w 2 dx = 2 g 3 w 2 dx Ω 5, for all w 2 L 2 (Ω) Ω 11
12 Similarly for the stresses, choosing arbitrary test functions with polynomial profile, we gather that A σ : τ 2εθ div τ dx = ν (σ g 3 ) tr(τ E ) dx Ω Ω for all τ H1 (Ω), ( ) 1 6 Ω 2µ 5 σ + θ τ + ω div τ dx = 0 for all τ H 1 (Ω), ω 2 τ 33 dx = ε ( 1 Ω E 3 g σ 33 ν ) 6 tr(σ) τ 33 dx Ω for all τ 33 H 1 (Ω). Integrating by parts, we can write the above equations in terms of the displacement unknowns only. 12
13 Reissner Mindlin Model: Find θ and ω such that 2ε3 3α div e (θ ) + 4εµ5 6 (θ ω) = ε εµ 5 6 div(θ ω) = 2g 3 in Ω, θ = 0, ω = 0 on Ω. β α g 3 in Ω, After that we can perform the following computations. σ = 2εA e (θ) λ 2µ + λ g 3δ, σ 33 = 2 5 g 3, σ = 2µ5 6 ( θ + ω), ω 2 = ε ( 1 E 3 g σ 33 ν ) 6 tr(σ). 13
14 Convergence of Reissner Mindlin Theorem 1 Assume g 3 = ε α ĝ 3, where ĝ 3 is ε-independent and smooth, and α is a nonnegative number. Then u ε u R E(P ε ) u ε E(P ε ) + σε σ R L 2 (P ε ) σ ε L 2 (P ε ) Cε 1/2. The constant C depends on Ω, and ĝ 3 only. 14
15 The proof of this theorem relies on the two-energies principle. This principle says that if, given a stress that is statically admissible, and a displacement that is kinematically admissible, then the sum of the squared energy and complementary energy norms is equal to the complementary energy norm of the constitutive residual. The complementary energy norm is simply [ σ C ε = (C σ) : σ dx P ε ] 1/2 15
16 Theorem 2 (The two energies principle.) Suppose that σ H(div, P ε ), is statically admissible, i.e. divσ = 0 in P ε, σn = g ε on P ε ±, and suppose u H 1 (P ε ) is kinematically admissible, i.e. u = 0 on P ε L. Then u ε u 2 E ε + σε σ 2 C ε = σ Ce(u) 2 C ε. 16
17 Our derivation of the Reissner Mindlin system yields an statically admissible stress field. Also, we add a boundary layer to make the displacement field kinematically admissible. 17
18 The constitutive residual ρ = C σ R e(u R ) is given by ρ = 0, ρ = 5 8µ (x2 3 ε )[µ(θ w)], ρ 33 = x 3 2µ + λ (λ div θ + g 3). We shall need the following a priori estimates: Lemma 3 Let γ = µ(θ w). Then there exists a constant C only dependent on Ω such that θ H 1 + w H 1 + ε γ L 2 C g 3 L 2. 18
19 The 3D Problem and its Modeling Full Discretization Playing Around with the Timoshenko Beam Equations Full Discretization Results for Constant Traction Conclusions 19
20 Consider the beam (2D elasticity) problem of finding such that u ε = u 1 and σ ε = σ 1,1 σ 1,2 u 2 σ 1,2 σ 2,2 div σ ε = 0, σ ε = C e(u ε ), in (, 1) ( ε, ε), σ ε ε n = g on (, 1) { ε, ε}, u ε = 0 on {, 1} ( ε, ε), where g ε = (0, g 2 ) is the traction, and e (u ε ) = 1 2 ( u ε + T u ε ). 20
21 Hellinger-Reissner Principle: (σ ε, u ε ) is critical point of L(τ, v ) = 1 2 on S g (P ε ) L 2 (P ε ). P ε C τ : τ dx + div P ε τ v dx Full discretization by choosing σ = σ 11 σ 12 P 0 2 (, 1) x 2 P0 1 (, 1) q(x 2 ) σ 12 σ 22 P(, 0 1) s(x 2 ) u = u 1 P (, 1 1) x 2 P(, 0 1) {1, r(x 2 )} u 2 P k 0 (, 1): continuous piecewise polynomials of degree k in (, 1) P k (, 1): discontinuous polynomials of degree k in (, 1) 21
22 Our (Timoshenko Beam Model) solution will be of the form σ T (x 1, x 2 ) = σ 11(x 1 )x 2 σ 12 (x 1 )q(x 2 ), g 3 x 3 /ε + σ 22 s(x 2 ) θ(x u T (x 1, x 2 ) = 1 )x 2, ω(x 1 ) + ω 2 (x 1 )r(x 2 ) where σ 11 P0 2 (, 1), σ 12 P0 1 (, 1), σ 22 P(, 0 1), θ P(, 1 1), ω P(, 0 1) 22
23 ε ε 1 dσ 11 dx 1 ψ dx 1 + 2ε dσ 12 dx 1 w dx 1 = σ 22 w 2 dx 1 = σ 12 ψ dx 1 = 0 ψ P 1 (, 1), g 2 w dx 1 w P 0 (, 1), g 2 w 2 dx 1 w 2 P 0 (, 1), ασ 11 τ 11 βσ 22 τ 11 2εθ dτ 11 dx 1 = 2β g 2 τ 11 dx 1, dx 1 1 ( ) 1 6 2µ 5 σ 12 + θ τ 12 + ω dτ 12 dx 1 = 0, dx 1 1 ( )[ 1 1 ω 2 τ 22 dx 1 = ε 2µ β 3 g ] 21 σ 22 τ 22 β ε 6 σ 11τ 22 dx 1, for all τ 11 P 2 0 (, 1), τ 12 P 1 0 (, 1), τ 22 P 0 (, 1). 23
24 Assume g 2 constant by parts. It is possible to solve for the stresses using the equations above: σ h 11 = 6 ε 2 P ε c 0x 1 + c 1, σ h 12 = P ε + c 0, σ h 22 = 2 5 g 2, where P (x) = x g 2(ξ) dξ, P 2 (x) = x P (ξ) dξ, and we compute c 0 and c 1 from (1+ 5µα 3ε 2 )c 0 = 1 2ε 1 P dx 1 5µα 2ε 3 1 P 2 x 1 dx 1 βµ 1 g 2 x 1 dx 1. 2αc 1 = β g 2 dx 1 6α ε 2 1 P 2 dx 1. 24
25 Hence, we have to solve for θ P 1 (, 1) and ω P 0 (, 1), where 1 1 2εθ dτ 11 dx 1 = dx 1 1 ω dτ 12 dx 1 dx 1 = 1 (βσ 22 ασ βg 2 ) τ 11 dx 1, ( 1 2µ for all τ 11 P 2 0 (, 1), τ 12 P 1 0 (, 1). ) 6 5 σ 12 + θ τ 12 dx 1, 25
26 For the plate problem, things are not that easy. Indeed, the planar stress is no longer a scalar, but a 2 2 symmetric matrix with rows in H(div, P ε ). It then becomes nontrivial to find stable elements for stresses and displacements. An option is to solve for the planar stress in terms of the displacement and shear stress. This is the same as start with Reissner Mindlin, and include shear stress, and it s traditionally considered to develop stable elements. 26
27 The 3D Problem and its Modeling Full Discretization Playing Around with the Timoshenko Beam Conclusions 27
28 Conclusions The derivation of the Reissner Mindlin Model using Hellinger Reissner Principle is mathematically consistent. The model obtained is actually the first in a hierarchical sequence. It also makes the derivation of the modeling error through the two energies principle easier, since it yields statically admissible stress field. Reissner Mindlin is better than the biharmonic: if there is shear, the biharmonic fails to deliver a meaningful result. 28
29 Example: General load Domain: (0, 1) ( ε, ε), with ε = 1/40 Traction: g ε = (1, 10 3 ) on the top and g ε = (, 10 3 ) on the bottom Clamped lateral boundary conditions 29
30 Planar elasticity, Reissner Mindlin and Kirchhoff Love solutions:
31 Reissner Mindlin is the same as the biharmonic: for purely transverse load, both models converge at the same rate. I now of no proof indicating that one is better than the other in this case. Both suffer from a quite low convergence rate: O(ε 1/2 ) in relative energy norm. Is there a way to improve this rate? 31
32 Thank You! 32
Questions (And some Answers) (And lots of Opinions) On Shell Theory
1 Questions (And some Answers) (And lots of Opinions) On Shell Theory 2 Shells Q: What is a shell? A: A three-dimensional elastic body occupying a thin neighborhood of a two-dimensional submanifold of
More informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationON THE ACCURACY OF REISSNER MINDLIN PLATE MODEL FOR STRESS BOUNDARY CONDITIONS. Sheng Zhang 1
ESAIM: M2AN Vol. 40, N o 2, 2006, pp. 269 294 DOI: 10.1051/m2an:2006014 ESAIM: Mathematical Modelling and Numerical Analysis www.edpsciences.org/m2an ON THE ACCURACY OF REISSNER MINDLIN PLATE MODEL FOR
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 informationMathematical modelling and numerical analysis of thin elastic shells
Mathematical modelling and numerical analysis of thin elastic shells Zhang Sheng Wayne State University, Michigan Zhang Sheng (szhang@wayne.edu) Thin elastic shell May, 202 / 9 Models for mechanics of
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationOn the Justification of Plate Models
On the Justification of Plate Models Dietrich Braess Stefan Sauter Christoph Schwab October 26, 2009 Abstract In this paper, we will consider the modelling of problems in linear elasticity on thin plates
More informationLocking phenomena in Computational Mechanics: nearly incompressible materials and plate problems
Locking phenomena in Computational Mechanics: nearly incompressible materials and plate problems C. Lovadina Dipartimento di Matematica Univ. di Pavia IMATI-CNR, Pavia Bologna, September, the 18th 2006
More informationA LINEAR SHELL THEORY BASED ON VARIATIONAL PRINCIPLES
The Pennsylvania State University The Graduate School Department of Mathematics A LINEAR SHELL THEORY BASED ON VARIATIONAL PRINCIPLES AThesisin Mathematics by Sheng Zhang c 2001 Sheng Zhang Submitted in
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationHierarchical Modeling of Heterogeneous Plates 1
TEMA Tend. Mat. Apl. Comput., 8, No. 2 (27), 219-228. c Uma Publicação da Sociedade Brasileira de Matemática Aplicada e Computacional. Hierarchical Modeling of Heterogeneous Plates 1 A.C. CARIUS 2, A.L.
More informationThéories des Coques Monique DAUGE Rennes, France http://www.maths.univ-rennes1.fr/~dauge Shell theory attempts the impossible: to provide a two-dimensional representation of an intrinsically three-dimensional
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 informationSEMIGROUP WELL-POSEDNESS OF MULTILAYER MEAD-MARKUS PLATE WITH SHEAR DAMPING
SEMIGROUP WELL-POSEDNESS OF MULTILAYER MEAD-MARKUS PLATE WITH SHEAR DAMPING Scott W. Hansen Department of Mathematics Iowa State University Ames, IA 50011 shansen@iastate.edu Abstract A multilayer plate
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
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 information4. Mathematical models used in engineering structural analysis
4. Mathematical models used in engineering structural analysis In this chapter we pursue a formidable task to present the most important mathematical models in structural mechanics. In order to best situate
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
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 informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationASYMPTOTICS AND HIERARCHICAL MODELING OF THIN DOMAINS
The Pennsylvania State University The Graduate School Department of Mathematics ASYMPTOTICS AND HIERARCHICAL MODELING OF THIN DOMAINS AThesisin Mathematics by Alexandre Loureiro Madureira c 1999 Alexandre
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 informationLinear Cosserat elasticity, conformal curvature and bounded stiffness
1 Linear Cosserat elasticity, conformal curvature and bounded stiffness Patrizio Neff, Jena Jeong Chair of Nonlinear Analysis & Modelling, Uni Dui.-Essen Ecole Speciale des Travaux Publics, Cachan, Paris
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationTangential-Displacement and Normal-Normal-Stress Continuous Mixed Finite Elements for Elasticity
www.oeaw.ac.at angential-displacement and Normal-Normal-Stress Continuous Mixed inite Elements for Elasticity J. Schöberl, A. Sinwel RICAM-Report 2007-10 www.ricam.oeaw.ac.at ANGENIAL-DISPLACEMEN AND NORMAL-NORMAL-SRESS
More informationMATHEMATICAL AND NUMERICAL ANALYSIS OF SOME PARAMETER-DEPENDENT MODELS
MATHEMATICAL AND NUMERICAL ANALYSIS OF SOME PARAMETER-DEPENDENT MODELS A SUMMARY OF MY RECENT WORK SHENG ZHANG 1. A general formulation of the problems Much of my recent work is on mathematical models
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 information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationHybrid equilibrium plate elements of high degree
Hybrid equilibrium plate elements of high degree Edward A.W. Maunder Department of Engineering, School of Engineering and Computer Science University of Exeter, Harrison Building, North Park Road, Exeter,
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationFurther Linear Elasticity
Torsion of cylindrical bodies Further Linear Elasticity Problem Sheet # 1. Consider a cylindrical body of length L, the ends of which are subjected to distributions of tractions that are statically equivalent
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 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 informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT METHOD FOR REISSNER-MINDLIN PLATES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 1, Pages 48 62 c 2017 Institute for Scientific Computing and Information ANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationHIERARCHICAL MODELING OF THE HEAT EQUATION IN A HETEROGENEOUS PLATE
HIERARCHICAL MODELING OF THE HEAT EQUATION IN A HETEROGENEOUS PLATE ANA CAROLINA CARIUS ALEXANDRE L. MADUREIRA Abstract. In this work we investigate the modeling of heterogeneous plates, where the length
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationLinear elastic analysis of thin laminated beams with uniform and symmetric cross-section
Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,
More informationANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES
ANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES F. Auricchio Dipartimento di Meccanica Strutturale Università di Pavia, Italy C. Lovadina Dipartimento di Ingegneria Meccanica e Strutturale
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 informationA-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements Joint work with Lourenco Beirão da Veiga (Milano), Claudia Chinosi (Alessandria) and Carlo Lovadina (Pavia) Previous
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 informationCorrection of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams
Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,
More informationFree vibration analysis of beams by using a third-order shear deformation theory
Sādhanā Vol. 32, Part 3, June 2007, pp. 167 179. Printed in India Free vibration analysis of beams by using a third-order shear deformation theory MESUT ŞİMŞEK and TURGUT KOCTÜRK Department of Civil Engineering,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationTwo Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate
Two Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate Pingbing Ming and Zhong-ci Shi Institute of Computational Mathematics & Scientific/Engineering Computing, AMSS, Chinese Academy of
More informationAsymptotics of Arbitrary Order for a Thin Elastic Clamped Plate, I. Optimal Error Estimates
Asymptotics of Arbitrary Order for a Thin Elastic Clamped Plate, I. Optimal Error Estimates Monique Dauge & Isabelle Gruais Abstract. This paper is the first of a series of two, where we study the asymptotics
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationOn the characterization of drilling rotation in the 6 parameter resultant shell theory
On the characterization of drilling rotation in the 6 parameter resultant shell theory Mircea Birsan and Patrizio Neff Chair for Nonlinear Analysis and Modelling Faculty of Mathematics, University Duisburg-Essen,
More informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationMODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM
MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,
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 informationEML4507 Finite Element Analysis and Design EXAM 1
2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationFinite Element Exterior Calculus and Applications Part V
Finite Element Exterior Calculus and Applications Part V Douglas N. Arnold, University of Minnesota Peking University/BICMR August 15 18, 2015 De Rham complex 0 H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div
More informationMATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command.
MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command Synopsis The MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command is used to specify the parameters associated with an isotropic, linear elastic material idealization
More informationarxiv: v1 [math.na] 20 May 2018
arxiv:1805.07835v1 [math.na] 20 May 2018 An ultraweak formulation of the Kirchhoff Love plate bending model and DPG approximation homas Führer Norbert Heuer Antti H. Niemi Abstract We develop and analyze
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationA GENERAL NONLINEAR THIRD-ORDER THEORY OF FUNCTIONALLY GRADED PLATES
International Journal of Aerospace and Lightweight Structures Vol. 1, No. 1 011 1 1 c Research Publishing Services DOI: 10.3850/S010486110000X A GENERAL NONLINEAR THIRD-ORDER THEORY OF FUNCTIONALLY GRADED
More informationAnalysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load
Analysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load Dr. Mama, Benjamin Okwudili Department of Civil Engineering, University of Nigeria, Nsukka. benjamin.mama@unn.edu.ng
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
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 informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More informationKirchhoff Plates: Field Equations
20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the
More information2008 by authors and 2008 Springer Science+Business Media
Antti H. Niemi, Harri Hakula, and Juhani Pitkäranta. 28. Point load on a shell. In: Karl Kunisch, Günther Of, and Olaf Steinbach (editors). Numerical Mathematics and Advanced Applications. Proceedings
More informationSEVERAL RELATED MODELS FOR MULTILAYER SANDWICH PLATES. Scott W. Hansen Iowa State University, Ames, Iowa, USA
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company SEVERAL RELATED MODELS FOR MULTILAYER SANDWICH PLATES Scott W. Hansen Iowa State University, Ames, Iowa, 500 USA
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 informationTRANSIENT RESPONSE OF SANDWICH AND LAMINATED COMPOSITES WITH DAMPING UNDER IMPULSE LOADING
TRANSIENT RESPONSE OF SANDWICH AND LAMINATED COMPOSITES WITH DAMPING UNDER IMPULSE LOADING Evgeny Barkanov, Andris Chate and Rolands Rikards Institute of Computer Analysis of Structures, Riga Technical
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 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 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 study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation
Multibody Syst Dyn (2014) 31:309 338 DOI 10.1007/s11044-013-9383-6 A study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation Marko K. Matikainen Antti
More informationBEAM DEFLECTION THE ELASTIC CURVE
BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment
More informationA UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE*
SIAM J. NUMER. ANAL. c 989 Society for Industrial and Applied Mathematics Vol. 26, No. 6, pp. 5, December 989 000 A UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE* DOUGLAS N. ARNOLD
More informationHierarchic Isogeometric Large Rotation Shell Elements Including Linearized Transverse Shear Parametrization
Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry October 11-13, 2017 in Stuttgart, Germany Hierarchic Isogeometric Large Rotation Shell
More information1 Bending of beams Mindlin theory
1 BENDNG OF BEAMS MNDLN THEORY 1 1 Bending of beams Mindlin theory Cross-section kinematics assumptions Distributed load acts in the xz plane, which is also a plane of symmetry of a body Ω v(x = 0 m Vertical
More informationDynamics of beams. Modes and waves. D. Clouteau. September 16, Department of Mechanical and Civil Engineering Ecole Centrale Paris, France
Dynamics of and waves Department of Mechanical and Civil Engineering Ecole Centrale Paris, France September 16, 2008 Outline 1 2 Pre-stressed Bernoulli-Euler in moving frames Kinematical assumption Resultant
More informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
More informationMathematical model of static deformation of micropolar elastic circular thin bar
Mathematical model of static deformation of micropolar elastic circular thin bar Mathematical model of static deformation of micropolar elastic circular thin bar Samvel H. Sargsyan, Meline V. Khachatryan
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 informationFinite Element Approximation of Piezoelectric Plates
Finite Element Approximation of Piezoelectric Plates F. Auricchio Department of Structural Mechanics, University of Pavia, 27100 Pavia, Italy P. Bisegna Department of Civil Engineering, University of Roma
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 informationAvailable online at ScienceDirect. Procedia IUTAM 13 (2015 ) 82 89
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 13 (215 ) 82 89 IUTAM Symposium on Dynamical Analysis of Multibody Systems with Design Uncertainties The importance of imperfections
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 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 informationProfessor Terje Haukaas University of British Columbia, Vancouver The M4 Element. Figure 1: Bilinear Mindlin element.
Professor Terje Hakaas University of British Colmbia, ancover www.inrisk.bc.ca The M Element variety of plate elements exist, some being characterized as Kirchhoff elements, i.e., for thin plates, and
More informationAN ITERATION METHOD OF SOLUTION OF A NONLINEAR EQUATION FOR A STATIC BEAM. 2 University Str., Tbilisi 0186, Georgia 2 Georgian Technical University
AN ITERATION METHOD OF SOLUTION OF A NONLINEAR EQUATION FOR A STATIC BEAM T. Davitashvili 1, J. Peradze 1,, Z. Tsiklauri 1 Ivane Javakhishvili Tbilisi State University University Str., Tbilisi 186, Georgia
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More information