Towards effective shell modelling with the FEniCS project
|
|
- Lee Grant
- 5 years ago
- Views:
Transcription
1 Towards effective shell modelling with the FEniCS project *, P. M. Baiz Department of Aeronautics 19th March 2013 Shells and FEniCS - FEniCS Workshop 13 1
2 Outline Introduction Shells: chart shear-membrane-bending and membrane-bending models example forms Two proposals for discussion: geometry: chart object for describing shell geometry discretisation: projection/reduction operators for implementation of generalised displacement methods Summary Shells and FEniCS - FEniCS Workshop 13 2
3 So far... dolfin manifold support already underway, merged into trunk [Marie Rognes, David Ham, Colin Cotter] I have already implemented locking-free (uncurved) beams and plate structures using dolfin manifold Next step: curved surfaces, generalised displacement methods (?) Aim of my talk is to start discussion on the best path Shells and FEniCS - FEniCS Workshop 13 3
4 Why shells? The mathematics: Shells are three-dimensional elastic bodies which occupy a thin region around a two-dimensional manifold situated in three-dimensional space The practical advantages: Shell structures can hold huge applied loads over large areas using a relatively small amount of material. Therefore they are used abundantly in almost all areas of mechanical, civil and aeronautical engineering. The computational advantages: A three-dimensional problem is reduced to a two-dimensional problem. Quantities of engineering relevance are computed directly. Shells and FEniCS - FEniCS Workshop 13 4
5 Figure : British Museum Great Court. Source: Wikimedia Commons. Shells and FEniCS - FEniCS Workshop 13 5
6 Figure : Specialized OSBB bottom bracket. Source: bikeradar.com Shells and FEniCS - FEniCS Workshop 13 6
7 A huge field There are many different ways of: obtaining shell models representing the geometry of surfaces on computers discretising shell models successfully And therefore we need suitable abstractions to ensure generality and extensibility of any shell modelling capabilities in FEniCS. Shells and FEniCS - FEniCS Workshop 13 7
8 Two methodologies Mathematical Model approach: 1. Derive a mathematical shell model. 2. Discretise that model using appropriate numerical method for description of geometry and fields. Degenerated Solid approach: 1. Begin with a general 3D variational formulation for the shell body. 2. Degenerate a solid 3D element by inferring appropriate FE interpolation at a number of discrete points. 3. No explicit mathematical shell model, one may be implied. Shells and FEniCS - FEniCS Workshop 13 8
9 Mathematical model Shells and FEniCS - FEniCS Workshop 13 9
10 Mathematical model shear-membrane-bending (smb) model Find U V smb : h 3 A b (U, V) + ha s (U, V) + ha m (U, V) = F(V) V V smb (1) membrane-bending (mb) model Find U V mb : h 3 A b (U, V) + ha m (U, V) = F(V) V V mb (2) Shells and FEniCS - FEniCS Workshop 13 10
11 Discretisation smb models vs mb models smb model takes into account the effects of shear; closer to the 3D solution for thick shells, matches the mb model for thin shells. Boundary conditions are better represented in smb model; hard and soft supports, boundary layers. smb U H 1 (Ω) vs mb U H 2 (Ω) Shells and FEniCS - FEniCS Workshop 13 11
12 Mathematical model Let s just take a look at the bending bilinear form A b for the mb model: A b (U, V) = ρ αβ H αβγδ b ρ γδ da (3a) Ω ρ αβ := ϕ,αβ t 1 j (u,1 (ϕ,2 t) u,2 (ϕ,1 t)) + 1 j (u,1 (ϕ,αβ ϕ,2 u,2 (ϕ,αβ ϕ,1 )) u,αβ t H αβγδ b := Eh 3 ( ) 12(1 ν 2 ν(ϕ,α ϕ,β ) +... ) (3b) (3c) Shells and FEniCS - FEniCS Workshop 13 12
13 Geometry Continuous model Terms describing the differential geometry of the shell mid-surface. The mid-surface is defined by the chart function. Discrete model We do not (usually) have an explicit representation of the chart. It must be constructed implicitly from the mesh and/or data from a CAD model. There are many different ways of doing this. Shells and FEniCS - FEniCS Workshop 13 13
14 Geometry Proposal 1 A base class Chart object which exposes various new symbols describing the geometry of the shell surface. Specific subclasses of Chart will implement a particular computational geometry procedure. The user can then express their mathematical shell model independently from the underlying geometrical procedure using the provided high-level symbols. Shells and FEniCS - FEniCS Workshop 13 14
15 shell_mesh = mesh (" shell. xml ") normals = MeshFunction (...) C = FunctionSpace ( shell_mesh, " CG", 2) chart = Chart ( shell_mesh, C, method =" patch_averaged ") chart = Chart ( shell_mesh, C, method =" CAD_normals ", normals = normals )... b_cnt = chart. contravariant_ basis () b_cov = chart. covariant_ basis () da = chart. measure () a = chart. first_ fundamental_ form ()... A_b =... Shells and FEniCS - FEniCS Workshop 13 15
16 Current discretisation options mb model: H 2 (Ω) conforming finite elements DG methods smb model: straight H 1 (Ω) conforming finite elements mixed finite elements (CG, DG) generalised displacement methods Shells and FEniCS - FEniCS Workshop 13 16
17 A quick note For simplicity I will just talk about the smb model reduced to plates, the chart function is the identity matrix; considerably simpler asymptotic behaviour but concepts apply to shells also. Shells and FEniCS - FEniCS Workshop 13 17
18 Locking 10 0 e L t = 0.1 t = 0.01 t = dofs Locking Inability of the basis functions to represent the limiting Kirchhoff mode. Shells and FEniCS - FEniCS Workshop 13 18
19 Move to a mixed formulation Treat the shear stress as an independent variational quantity: γ h = λ t 2 ( z 3h θ h ) S h Discrete Mixed Weak Form Find (z 3h, θ h, γ h ) (V 3h, R h, S h ) such that for all (y 3h, η, ψ) (V 3h, R h, S h ): a b (θ h ; η) + (γ h ; y 3 η) L 2 = f (y 3 ) t ( z 3h θ h ; ψ) L 2 2 λ (γ h; ψ) L 2 = 0 Shells and FEniCS - FEniCS Workshop 13 19
20 Move back to a displacement formulation Linear algebra level: Eliminate the shear stress unknowns a priori to solution [ ] { } { } A B u f B T = (5) C γ 0 To do this we can rearrange the second equation and then if and only if C is diagonal/block-diagonal we can invert cheaply giving a problem in original displacement unknowns: (A + BC 1 B T )u = f (6) Shells and FEniCS - FEniCS Workshop 13 20
21 Currently, this can be done with CBC.Block TRIA0220, TRIA1B20 [Arnold and Brezzi][Boffi and Lovadina] David Ham, Kent Andre-Mardal, Anders Logg, Joachim Haga and myself A, B, BT, C = [ assemble (a), assemble (b), assemble (bt), assemble (c)] K = collapse ( A - B * LumpedInvDiag ( C) * BT) Shells and FEniCS - FEniCS Workshop 13 21
22 Move back to a displacement formulation Variational form level: γ h = λ t 2 Π h ( z 3h θ h ) (7) Π h : V 3h R h S h (8) giving: a b (θ h ; η) + (Π h ( z 3 θ); y 3 η) L 2 = f (y 3 ) (9) Shells and FEniCS - FEniCS Workshop 13 22
23 Proposal 2 A new class Projection in UFL that signals to FFC that a projection between FunctionSpace objects is required. This requires additions in DOLFIN, UFL, FFC and FIAT. Shells and FEniCS - FEniCS Workshop 13 23
24 MITC7 Shells and FEniCS - FEniCS Workshop 13 24
25 ... V_3 = FunctionSpace (mesh, "CG", 2) R = VectorFunctionSpace ( mesh, " CG", 2, dim =2) + VectorFunctionSpace (mesh,"b", 3, dim =2) S = FunctionSpace ( mesh, " N1curl ", order =2) Pi_h = Projection ( from =R, to =S)... U = MixedFunctionSpace ([R, V_3 ]) theta, z_3 = TrialFunctions ( U) eta, y_3 = TestFunctions ( U) a_s = inner ( grad ( z_3 ) - Pi_h ( theta ), grad ( y_3 ) - Pi_h ( eta ))*dx Shells and FEniCS - FEniCS Workshop 13 25
26 Summary A big field with lots of approaches; need appropriate abstractions (inc. other PDEs on surfaces) Proposal 1: Expression of geometric terms in shell models using a natural form language which reflects the underlying mathematics Proposal 2: Effective discretisation options for the implementation of generalised displacement methods Shells and FEniCS - FEniCS Workshop 13 25
27 Thanks for listening. Shells and FEniCS - FEniCS Workshop 13 26
FEniCS Course. Lecture 8: A posteriori error estimates and adaptivity. Contributors André Massing Marie Rognes
FEniCS Course Lecture 8: A posteriori error estimates and adaptivity Contributors André Massing Marie Rognes 1 / 24 A priori estimates If u H k+1 (Ω) and V h = P k (T h ) then u u h Ch k u Ω,k+1 u u h
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationThe FEniCS Project. Anders Logg. Simula Research Laboratory / University of Oslo. Det Norske Veritas
The FEniCS Project Anders Logg Simula Research Laboratory / University of Oslo Det Norske Veritas 2011 05 11 The FEniCS Project Free Software for Automated Scientific Computing C++/Python library Initiated
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 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 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 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 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 informationFEniCS Course. Lecture 1: Introduction to FEniCS. Contributors Anders Logg André Massing
FEniCS Course Lecture 1: Introduction to FEniCS Contributors Anders Logg André Massing 1 / 15 What is FEniCS? 2 / 15 FEniCS is an automated programming environment for differential equations C++/Python
More informationFEniCS Course Lecture 13: Introduction to dolfin-adjoint. Contributors Simon Funke Patrick Farrell Marie E. Rognes
FEniCS Course Lecture 13: Introduction to dolfin-adjoint Contributors Simon Funke Patrick Farrell Marie E. Rognes 1 / 25 dolfin adjoint automatic adjoint models for FEniCS P. E. Farrell, S. W. Funke, D.
More informationFEniCS Course. Lecture 1: Introduction to FEniCS. Contributors Anders Logg André Massing
FEniCS Course Lecture 1: Introduction to FEniCS Contributors Anders Logg André Massing 1 / 30 What is FEniCS? 2 / 30 FEniCS is an automated programming environment for differential equations C++/Python
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 informationContinuum mechanics of beam-like structures using one-dimensional finite element based on Serendipity Lagrange cross-sectional discretisation, Mayank Patni, Prof. Paul Weaver, Dr Alberto Pirrera Bristol
More informationGETTING STARTED WITH FENICS
GETTING STARTED WITH FENICS DOUGLAS N. ARNOLD 1. A first program in FEniCS 1.1. The boundary value problem. As a first problem we consider the Neumann boundary value problem: (1) u + u = f in, u = 0 on,
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationSolving an elasto-plastic model using DOLFIN
Solving an elasto-plastic model using DOLFIN TTI 2005 Johan Jansson johanjan@math.chalmers.se Chalmers University of Technology Solving an elasto-plastic model using DOLFIN p. Overview Motivation Previous
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 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 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 Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
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 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 informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
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 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 informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
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 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 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 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 informationPlates and Shells: Theory and Computation. Dr. Mostafa Ranjbar
Plates and Shells: Theory and Computation Dr. Mostafa Ranjbar Outline -1-! This part of the module consists of seven lectures and will focus on finite elements for beams, plates and shells. More specifically,
More informationINTRODUCTION TO THE FINITE ELEMENT SOFTWARE FEniCS ROLAND HERZOG
INTRODUCTION TO THE FINITE ELEMENT SOFTWARE FEniCS ROLAND HERZOG 1 Introduction FEniCS is free software for automated solution of differential equations. As a particular feature, it allows the formulation
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 informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
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 informationEfficient multigrid solvers for mixed finite element discretisations in NWP models
1/20 Efficient multigrid solvers for mixed finite element discretisations in NWP models Colin Cotter, David Ham, Lawrence Mitchell, Eike Hermann Müller *, Robert Scheichl * * University of Bath, Imperial
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 informationSolving Non-linear Elasticity BVPs in FEniCS
Solving Non-linear Elasticity BVPs in FEniCS Bifurcations and Buckling of Beams Christiana Mavroyiakoumou University of Oxford A special topic report submitted for the degree of M.Sc. in Mathematical Modelling
More informationMulti-scale approach of the mechanical behavior of RC structures Application to nuclear plant containment buildings
Multi-scale approach of the mechanical behavior of RC structures Application to nuclear plant containment buildings Martin DAVID June 19, 2012 1 Industrial context EDF is responsible for numerous reinforced
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 informationTopology Optimization Mathematics for Design
Topology Optimization Mathematics for Design Homogenization Design Method (HMD) Why topology? Change in shape & size may not lead our design criterion for reduction of structural weight. 1 Structural Design
More informationModule-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination
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 informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationExpressive Environments and Code Generation for High Performance Computing
Expressive Environments and Code Generation for High Performance Computing Garth N. Wells University of Cambridge SIAM Parallel Processing, 20 February 2014 Collaborators Martin Alnæs, Johan Hake, Richard
More informationFEA CODE WITH MATLAB. Finite Element Analysis of an Arch ME 5657 FINITE ELEMENT METHOD. Submitted by: ALPAY BURAK DEMIRYUREK
FEA CODE WITH MATAB Finite Element Analysis of an Arch ME 5657 FINITE EEMENT METHOD Submitted by: APAY BURAK DEMIRYUREK This report summarizes the finite element analysis of an arch-beam with using matlab.
More 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 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 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 informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
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 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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMAE 323: Chapter 6. Structural Models
Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need
More informationAalto University School of Engineering
Aalto University School of Engineering Kul-4.4 Ship Structural Design (P) ecture 6 - Response of Web-frames, Girders and Grillages Kul-4.4 Ship Structures Response ecture 5: Tertiary Response: Bending
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 informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
More informationA NEW CLASS OF MIXED FINITE ELEMENT METHODS FOR REISSNER MINDLIN PLATES
A NEW CLASS OF IXED FINITE ELEENT ETHODS FOR REISSNER INDLIN PLATES C. LOVADINA Abstract. A new class of finite elements for Reissner indlin plate problem is presented. The family is based on a modified
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 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 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 informationThe FEniCS Project. Anders Logg. Simula Research Laboratory University of Oslo NOTUR / 34
The FEniCS Project Anders Logg Simula Research Laboratory University of Oslo NOTUR 2011 2011 05 23 1 / 34 What is FEniCS? 2 / 34 FEniCS is an automated programming environment for differential equations
More informationMITOCW MITRES2_002S10linear_lec07_300k-mp4
MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More informationACCELERATION WAVES IN THE VON KARMAN PLATE THEORY 1
ACCELERATION WAVES IN THE VON KARMAN PLATE THEORY V. VASSILEV and P. DJONDJOROV Institute of Mechanics Bulgarian Academy of Sciences Acad. G. Bontchev St., Block 4, 3 Sofia, BULGARIA arxiv:math-ph/000606v2
More informationFinite element modelling of structural mechanics problems
1 Finite element modelling of structural mechanics problems Kjell Magne Mathisen Department of Structural Engineering Norwegian University of Science and Technology Lecture 10: Geilo Winter School - January,
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationA free-vibration thermo-elastic analysis of laminated structures by variable ESL/LW plate finite element
A free-vibration thermo-elastic analysis of laminated structures by variable ESL/LW plate finite element Authors: Prof. Erasmo Carrera Dr. Stefano Valvano Bologna, 4-7 July 2017 Research group at Politecnico
More informationParametric study on the transverse and longitudinal moments of trough type folded plate roofs using ANSYS
American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-4 pp-22-28 www.ajer.org Research Paper Open Access Parametric study on the transverse and longitudinal moments
More 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 informationNew Model Stability Criteria for Mixed Finite Elements
New Model Stability Criteria for Mixed Finite Elements Andrew Gillette Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin http://www.math.utexas.edu/users/agillette
More informationMechanics of Irregular Honeycomb Structures
Mechanics of Irregular Honeycomb Structures S. Adhikari 1, T. Mukhopadhyay 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Singleton Park, Swansea SA2 8PP, UK Sixth International
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 informationComputational Reality XV. Rotary viscometer. Determining material coefficients I. B. Emek Abali, Wilhelm LKM - TU Berlin
Computational Reality XV Rotary viscometer Determining material coefficients I B. Emek Abali, Wilhelm Hübner @ LKM - TU Berlin Abstract A non-newtonian fluid can be modeled with three-parameter-constitutive
More informationThe FEniCS Project Philosophy, current status and future plans
The FEniCS Project Philosophy, current status and future plans Anders Logg logg@simula.no Simula Research Laboratory FEniCS 06 in Delft, November 8-9 2006 Outline Philosophy Recent updates Future plans
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 informationFinite Elements for Large Strains - A double mixed (M 2 ) Formulation
Finite Elements for Large Strains - A double mixed (M 2 ) Formulation Motivation Development of user friendly elements robustness simple treatment of incompressible materials complex geometries geometrical
More informationFinite Element Code Generation: Simplicity, Generality, Efficiency
: Simplicity, Generality, Efficiency Simula Research Laboratory SCSE 07 Uppsala, August 15 2007 Acknowledgments: Martin Sandve Alnæs, Johan Hoffman, Johan Jansson, Claes Johnson, Robert C. Kirby, Matthew
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 informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More 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 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 informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 008/9 Dr. Colin Caprani 1 Contents 1. Introduction... 3 1.1 General... 3 1. Background... 4 1.3 Discontinuity Functions...
More informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 009/10 Dr. Colin Caprani 1 Contents 1. Introduction... 4 1.1 General... 4 1. Background... 5 1.3 Discontinuity Functions...
More informationBECAS - an Open-Source Cross Section Analysis Tool
BECAS - an Open-Source Cross Section Analysis Tool José P. Blasques and Robert D. Bitsche Presented at DTU Wind Energy stand at the EWEA 2012 conference, Copenhagen, 16.4.2012 BECAS-DTUWind@dtu.dk Motivation
More informationANSYS Explicit Dynamics Update. Mai Doan
ANSYS Explicit Dynamics Update Mai Doan Mai.Doan@ansys.com +1 512 687 9523 1/32 ANSYS Explicit Dynamics Update Outline Introduction Solve Problems that were Difficult or Impossible in the Past Structural
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 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 informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationTHICK SHELL FINITE ELEMENTS IN THE ANALYSIS OF HELICOIDAL STAIR SLABS
00328 ivil and Environmental Engineering onference New Frontiers and hallenges 8-12 November 1999 Bangkok Thailand THK SHELL FNTE ELEMENTS N THE ANALYSS OF HELODAL STAR SLABS A. F. M. S. Amin & S. Ahmad
More informationSpinor Formulation of Relativistic Quantum Mechanics
Chapter Spinor Formulation of Relativistic Quantum Mechanics. The Lorentz Transformation of the Dirac Bispinor We will provide in the following a new formulation of the Dirac equation in the chiral representation
More informationFinite Elements with Symbolic Computations and Code Generation
Finite Elements with Symbolic Computations and Code Generation August 15, 2007 Outline We will focus on creating/assembling matrices based on FEM Manual implementation of the integrand/variational form
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More informationFEM Validation. 12th January David Schmid Teamleader Structural Analysis
FEM Validation 12th January 2012 David Schmid Teamleader Structural Analysis FEM Validation and Verification Each FE model which is used to substantiate flight material must be verified Depending on the
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationDESIGN OF LAMINATES FOR IN-PLANE LOADING
DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
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 information