A new approach for Kirchhoff-Love plates and shells
|
|
- Maximillian Freeman
- 5 years ago
- Views:
Transcription
1 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 Linz) Kirchhoff-Love plates and shells AANMPDE 10 1 / 30
2 Acknowledgment Joint work with Wolfgang Krendl, Katharina Rafetseder The research was funded by the Austrian Science Fund (FWF): S11702-N23 Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 2 / 30
3 Outline 1 The Kirchhoff plate bending problem 2 Extension to Mindlin-Reissner plates? 3 Extension to Kirchhoff-Love shells 4 Numerical results 5 Conclusion and outlook Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 3 / 30
4 Kirchhoff plate bending problem Figure: Geometry of a shell Unknown: covariant components of the displacement u = (u i ) Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 4 / 30
5 Kirchhoff plate bending problem plate: S is flat, for simplicity: θ = Id pure bending: u 3 = w. dimensional reduction of 3D elasticity for Ω = ( t/2, t/2) find w W H 2 () which minimizes J(w) = 1 C 2 w : 2 w dx 2 with Cε = 2µ ( ε + ν ) (tr ε)i 1 ν g w dx Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 5 / 30
6 Kirchhoff plate bending problem (kinematic) boundary conditions: = γ c γ s γ f if the plate is clamped on γ c simply supported on γ s then free on γ f W = {v H 2 () : v = 0, n v = 0 on γ c, v = 0 on γ s }. variational formulation: find w W such that C 2 w : 2 v dy = g v dx for all v W, Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 6 / 30
7 Kirchhoff plate bending problem quantities of interest: w vertical deflection M = C 2 w bending moments Mixed variational formulation: Find M V and w Q such that C 1 M : L dy + L : 2 w dy = 0 for all L V, M : 2 v dy = g v dx for all v Q, e.g., with the function spaces V = L 2 (Ω) sym, Q = W. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 7 / 30
8 Kirchhoff plate bending problem An different choice for Q with less smoothness: if W = {v H 2 () : v = 0, n v = 0 on γ c, v = 0 on γ s }, then Q = {v H 1 () : v = 0, on γ c, v = 0 on γ s } = H 1 0,γ c γ s (). The choice for V : transposition method b(v, L) = L : 2 v dy Bv, L B L, v, where B L, v is well-defined for v Q. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 8 / 30
9 Kirchhoff plate bending problem The linear operator B: B : D(B) Q Y, here with domain D(B) = W and Y = L 2 (Ω) sym, is a densely defined linear operator and unbounded operator in Q. Its adjoint operator is given by the identity B : D(B ) Y Q and all L from the domain B L, v = Bv, L for all v W D(B ) = {L Y : Bv, L c v Q for all v D(B)} Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE 10 9 / 30
10 Kirchhoff plate bending problem standard integration by parts: L : 2 v dy = Div L v dy + = Div L v dy + (Ln) v ds L nt t v ds + L nn n v ds?? L nn = 0 on γ s γ f (L nn = Ln n, L nt = Ln t) = γ s : D(B ) = {L L 2 (Ω) sym : div Div L H 1 ()} Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
11 Kirchhoff plate bending problem Lemma Let L L 2 (Ω) sym C 1 (). Then L D(B ) if and only if L nn = 0 on γ s γ f and L nt x = 0 on interior corner points x of γ f. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
12 Kirchhoff plate bending problem Notations for B = 2 = Grad grad: with B = div Div, D(B ) = H(div Div, ; Q ) sym. H(div Div, ; Q ) sym = {L L 2 (Ω) sym : L : 2 v dy c v H 1 () for all v W } Observe H 1 0 () sym H(div Div, ; Q ) sym L 2 (Ω) sym Bernardi/Girault/Maday (1992), Z. (2015), Pechstein/Schöberl (2011), Rafetseder/Z. (2017) Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
13 Kirchhoff plate bending problem Mixed variational formulation of Find M V and w Q such that C 1 M : L dy + div Div L, w = 0 for all L V, div Div M, v = g v dx for all v Q, with the function spaces V = H(div Div, ; Q ) sym, Q = H 1 0,γ c γ s (), equipped with norms v Q = v 1 and L V = L div Div;Q, given by L 2 div Div;Q = L 2 L 2 (Ω) + div Div L 2 Q. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
14 Kirchhoff plate bending problem Theorem The mixed problem is well-posed in V Q = H(div Div, ; Q ) sym H 1 0,γ c γ s (). Corollary The primal variational problem and the mixed problem are equivalent: If w W solves the primal problem, then M = C 2 w V and (M, w) solves the mixed problem. If (M, w) V Q solves the mixed problem, then w W and w solves the primal problem. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
15 Kirchhoff plate bending problem Theorem Let be simply connected. M H(div Div, ; Q ) sym if and only if M = p I + symcurl φ with with Ψ p = p Q = H 1 0,γ c γ s () and φ Ψ p { ψ (H 1 ()) 2 : t ψ, v Γ = p n v ds Γ } for all v W ( ) ( 1 φ1 φ =, sym Curl φ = 2 φ 1 2 ( ) 2φ 2 1 φ 1 ) φ ( 2φ 2 1 φ 1 ) 1 φ 2 Krendl/Rafetseder/Z. (2014,2016), Rafetseder/Z. (2017) Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
16 Kirchhoff plate bending problem Equivalent variational formulation in ( 2 p Q = H0,γ 1 c γ s (), φ Ψ p H ()) 1, w Q = H 1 (). 0,γc γs 1 Find p Q such that p v dy = g v dx for all v Q. 2 For given p Q, find φ Ψ p such that C 1 symcurl φ : symcurl ψ dy = G[p], ψ for all ψ Ψ 0. 3 For given M = pi + symcurl φ, find w Q such that w q dy = H[M], q for all q Q. Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
17 Decomposition of the problem 2 for φ = ( φ 2, φ 1 ) T we obtain Ĉ 1 ε(φ ) : ε(ψ) dy = Ĝ[p], ψ for all ψ Ψ 0 with an appropriately rotated material tensor Ĉ 1 and ε αβ (ψ) = 1 2 ( βψ α + α ψ β ) decomposition of the Kirchhoff model into 1 Poisson problem 2 plane linear elasticity problem 3 Poisson problem Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
18 Mindlin-Reissner plates Mindlin-Reissner model J = 1 { Cε(θ) : ε(θ) + µ t 2 w θ 2} dx 2 g w dx min Kirchhoff model: θ = w. J = 1 { } C 2 w : 2 w 2 dx g w dx min decomposition of the Mindlin-Reisser model into 1 Poisson problem 2 Stokes-like problem 3 Poisson problem Brezzi/Fortin 1986 Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
19 Kirchhoff-Love shells variational formulation of the linearized model: find u W such that ] [t Cε(u) : ε(v) + t3 a Cκ(u) : κ(v) dy = F, v for all v W 12 with the membrane strain the bending strain W H 1 () H 1 () H 2 (), ε αβ (u) = 1 2 (u α β + u β α ) b αβ u 3, κ αβ (u) = u 3 αβ b σ αb σβ u 3 + b σ αu σ β + b τ β u τ α + bτ β α u τ = αβ u Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
20 Kirchhoff-Love shells quantities of interest: u displacement M = t3 12 Cκ(u) a bending moments M = Ĉκ(u) Mixed variational formulation: Find M V and u Q such that Ĉ 1 M : L dy div Div L, u 3 L : κ 1 (u) dy = 0 div Div M, v 3 M : κ 1 (v) dy t Cɛ(u) : ɛ(v) a dy = F, v for all L V and v Q, with the function spaces where Q i H 1 (). V = H(div Div, ; Q 3 ) sym, Q = Q 1 Q 2 Q 3, Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
21 Numerical results benchmark examples of the shell obstacle course: Scordelis-Lo roof pinched cylinder pinched hemisphere midsurface represented as B-spline surfaces approximation spaces: equal-order approximation for u i, p, φ single patch: B-splines of maximum smoothness multipatch: C 0 implementation in G+Smo (C++ library for IgA) Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
22 Numerical results Example: Scordelis-Lo roof Figure: Geometry and load Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
23 u Numerical results Example: Scordelis-Lo roof 0.35 Scordelis-Lo roof p = 2 p = 3 p = 4 p = D.O.F. Figure: Results, 1 patch Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
24 u Numerical results Example: Scordelis-Lo roof 0.35 Scordelis-Lo roof (4 patches) p = 2 p = 3 p = 4 p = D.O.F. Figure: Results, 4 patches Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
25 Numerical results Example: pinched Cylinder Figure: Geometry and load Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
26 Numerical results Example: pinched Cylinder 2 Pinched cylinder (4 patches) # u p=2 p=3 p=4 p= D.O.F. Figure: Results, 4 patches Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
27 Numerical results Example: pinched Hemisphere Figure: Geometry and load Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
28 u Numerical results Example: pinched Hemisphere 0.1 Pinched hemisphere (4 patches) p = 2 p = 3 p = 4 p = D.O.F. Figure: Results, 4 patches Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
29 Conclusion and outlook Kirchhoff plate bending problems and similar 4-th order problems can be decomposed in three (consecutively to solve) second-order problems. extension to Kirchhoff-Love shell: formulation in H 1 spaces. work in progress: extend the numerical analysis from Kirchhoff plates to Kirchhoff-Love shells. behavior of equal-order discretization method w.r.t. membrane locking extension to Mindlin-Reissner plates and shells Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
30 References W. Krendl, K. Rafetseder, and W. Z. A decomposition result for biharmonic problems and the Hellan-Herrmann-Johnson method. ETNA, Electron. Trans. Numer. Anal., 45: , K. Rafetseder and W. Z. A decomposition result for Kirchhoff plate bending problems and a new discretization approach. arxiv: ArXiv e-prints, March Walter Zulehner (JKU Linz) Kirchhoff-Love plates and shells AANMPDE / 30
A new mixed isogeometric approach to Kirchhoff-Love shells. Katharina Rafetseder and Walter Zulehner. G+S Report No. 70
A new mixed isogeometric approach to Kirchhoff-Love shells Katharina Rafetseder and Walter Zulehner G+S Report No. 70 March 2018 A new mixed isogeometric approach to Kirchhoff-Love shells Katharina Rafetseder
More informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Decomposition Result for Kirchhoff Plate Bending Problems and a New Discretization Approach Katharina Rafetseder Institute of Computational
More informationSome Remarks on the Reissner Mindlin Plate Model
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, 2003 1 Outline The 3D Problem and its Modeling Full
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
More informationarxiv: v2 [math.na] 20 Dec 2017
arxiv:1703.07962v2 [math.na] 20 Dec 2017 A decomposition result for Kirchhoff plate bending problems and a new discretization approach Katharina Rafetseder and Walter Zulehner October 10, 2018 Abstract
More informationGustafsson, Tom; Stenberg, Rolf; Videman, Juha A posteriori analysis of classical plate elements
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Gustafsson, Tom; Stenberg, Rolf;
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 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 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 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 informationWellposedness for A Thermo-Piezo-Electric Coupling Model.
Wellposedness for A Thermo-Piezo-Electric Coupling Model. AANMPDE 216, Strobl, Austria Rainer Picard (joint work with T. Mulholland, S. Trostor & M. Waurick) Technische Universität Dresden, Germany 1/3
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 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 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 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 informationDevelopment DKMQ Shell Element with Five Degrees of Freedom per Nodal
International Journal of Mechanical Engineering and Robotics Research Vol. 6, No., May 7 Development DKMQ Shell Element with Five Degrees of Freedom per Nodal Herry Irpanni, Irwan Katili, and Imam J. Maknun
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 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 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 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 informationQuestions (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 informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
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 informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
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 informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems Ulrich Langer Institute of Computational Mathematics,
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 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 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 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 informationKirchhoff-Love shell theory based on tangential differential calculus
Kirchhoff-Love shell theory based on tangential differential calculus D. Schöllhammer, T.P. Fries arxiv:1805.11978v2 [cs.ce] 10 Oct 2018 October 11, 2018 Institute of Structural Analysis Graz University
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 informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
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 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 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 informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
More informationA spacetime DPG method for acoustic waves
A spacetime DPG method for acoustic waves Rainer Schneckenleitner JKU Linz January 30, 2018 ( JKU Linz ) Spacetime DPG method January 30, 2018 1 / 41 Overview 1 The transient wave problem 2 The broken
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 informationUNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING
UNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ACKNOWLEDGMENTS This report consists of the dissertation by Ms. Yan Jane Liu, submitted in partial fulfillment
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 informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationMultiscale Analysis of Vibrations of Streamers
Multiscale Analysis of Vibrations of Streamers Leszek Demkowicz Joint work with S. Prudhomme, W. Rachowicz, W. Qiu and L. Chamoin Institute for Computational Engineering and Sciences (ICES) The University
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for Kirchhoff plate elements Jarkko Niiranen Laboratory of Structural
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationFinite Element Exterior Calculus. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016
Finite Element Exterior Calculus Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016 The Hilbert complex framework Discretization of Hilbert complexes Finite element
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Computational Methods in Elasticity - Michel SALAUN
COMPUAIONAL MEHODS IN ELASICIY Michel SALAUN Départment de Mathématiques et Informatique, ENSICA, oulouse, France Keywords Elastodynamics, Energy minimization, Finite element method, Kirchhoff- Love model,
More informationRigorous derivation of a sixth order thin film equation
Rigorous derivation of a sixth order thin film equation Boris Muha Department of Mathematics, Faculty of Science, University of Zagreb Waves in Flows, Prague, August 2018 Joint work with M. Bukal, University
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 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 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 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 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 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 elements for plates and shells. Advanced Design for Mechanical System LEC 2008/11/04
Finite elements for plates and shells Sstem LEC 008/11/04 1 Plates and shells A Midplane The mid-surface of a plate is plane, a curved geometr make it a shell The thickness is constant. Ma be isotropic,
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationNonlinear Thermo- Mechanics of Plates and Shallow Shells
Nonlinear Thermo- Mechanics of Plates and Shallow Shells Payam Khazaeinejad 1, Asif S. Usmani 1, Omar Laghrouche 1 IIE, School of Engineering, The University of Edinburgh IIE, School of the Built Environment,
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationPhase-field description of brittle fracture in plates and shells
Phase-field description of brittle fracture in plates and shells Josef Kiendl a,, Marreddy Ambati a, Laura De Lorenzis a, Hector Gomez b, Alessandro Reali c,d a Institute of Applied Mechanics TU Braunschweig,
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 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 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 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 information3D-SHELL ELEMENTS AND THEIR UNDERLYING MATHEMATICAL MODEL
Mathematical Models and Methods in Applied Sciences Vol. 4, No. (2004) 05 42 c World Scientific Publishing Company 3D-SHELL ELEMENTS AND THEIR UNDERLYING MATHEMATICAL MODEL D. CHAPELLE and A. FERENT INRIA-Rocquencourt,
More informationTowards effective shell modelling with the FEniCS project
Towards effective shell modelling with the FEniCS project *, P. M. Baiz Department of Aeronautics 19th March 2013 Shells and FEniCS - FEniCS Workshop 13 1 Outline Introduction Shells: chart shear-membrane-bending
More informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
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 informationLecture Notes Introduction to Vector Analysis MATH 332
Lecture Notes Introduction to Vector Analysis MATH 332 Instructor: Ivan Avramidi Textbook: H. F. Davis and A. D. Snider, (WCB Publishers, 1995) New Mexico Institute of Mining and Technology Socorro, NM
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
An alternative boundary element formulation for plate bending analysis J.B. Paiva, L.O. Neto Structures Department, Sao Carlos Civil Engineering School, f, BrazzY Abstract This work presents an alternative
More informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
More informationANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company ANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL DOUGLAS N. ARNOLD Department of
More informationExternal Work. When a force F undergoes a displacement dx in the same direction i as the force, the work done is
Structure Analysis I Chapter 9 Deflection Energy Method External Work Energy Method When a force F undergoes a displacement dx in the same direction i as the force, the work done is du e = F dx If the
More informationLecture 3. Vector fields with given vorticity, divergence and the normal trace, and the divergence and curl estimates
Lecture 3. Vector fields with given vorticity, divergence and the normal trace, and the divergence and curl estimates Ching-hsiao (Arthur) Cheng Department of Mathematics National Central University Taiwan,
More informationarxiv: v2 [cs.ce] 22 Feb 2018
Isogeometric analysis of thin Reissner-Mindlin plates and shells: locking phenomena and B-bar method Qingyuan Hu a,b, Yang Xia c, Sundararajan Natarajan d, Andreas Zilian b, Ping Hu c, Stéphane P.A. Bordas
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 informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
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 informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Schur Complement Preconditioners for Multiple Saddle Point Problems of Block Tridiagonal Form with Application to Optimization Problems
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 informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
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 informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationTHE INFLUENCE OF LATERAL BOUNDARY CONDITIONS ON THE ASYMPTOTICS IN THIN ELASTIC PLATES
THE INFLUENCE OF LATERAL BOUNDARY CONDITIONS ON THE ASYMPTOTICS IN THIN ELASTIC PLATES MONIQUE DAUGE, ISABELLE GRUAIS, AND ANDREAS RÖSSLE Abstract. Here we investigate the limits and the boundary layers
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 informationROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS
ROBUST BDDC PRECONDITIONERS FOR REISSNER-MINDLIN PLATE BENDING PROBLEMS AND MITC ELEMENTS L. BEIRÃO DA VEIGA, C. CHINOSI, C. LOVADINA, AND L. F. PAVARINO Abstract. A Balancing Domain Decomposition Method
More informationA Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear Elasticity Equations
www.oeaw.ac.at A Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear lasticity quations B. Ayuso de Dios, I. Georgiev, J. Kraus, L. Zikatanov RICAM-Report 2011-29 www.ricam.oeaw.ac.at
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationNumerical results: Isotropic plates and shells
Chapter 5 Numerical results: Isotropic plates and shells In this chapter numerical results are presented for the isotropic plate and shell elements presented in Chapter 4. In the following, SA denotes
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 informationIsogeometric Kirchhoff Love shell formulations for general hyperelastic materials
Isogeometric Kirchhoff Love shell formulations for general hyperelastic materials Josef Kiendl a,, Ming-Chen Hsu b, Michael C. H. Wu b, Alessandro Reali a,c,d a Dipartimento di Ingegneria Civile ed Architettura,
More informationStress Analysis of Elastic Roof with Boundary Element Formulations
Copyright 211 Tech Science Press CMC, vol.26, no.3, pp.187-2, 211 Stress Analysis of Elastic Roof with Boundary Element Formulations Dan Ma 1,2 and Xianbiao Mao 1 Abstract: Roof is one of the most important
More informationEDGE EFFECTS IN THE REISSNER-MINDLIN PLATE THEORY*
EDGE EFFECTS IN THE REISSNER-MINDLIN PLATE THEORY* DOUGLAS N ARNOLD and RICHARD S FALK Abstract We study the behavior of solutions of five different boundary value problems for the Reissner- Mindlin plate
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More information