Discontinuous Galerkin methods for nonlinear elasticity
|
|
- Bernice Lora Green
- 5 years ago
- Views:
Transcription
1 Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008
2 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity and demonstrate some of their advantages over more traditional methods. Although the use of DG for elliptic problems may not readily appear advantageous since we typically expect solutions to be smooth, we shall see that when the constraint of inter-element continuity is relaxed not only are some problems easier to solve, but, in some cases, we may be able to solve them more accurately. For instance, consider a nearly incompressible neo-hookean material. It is well known that low order conforming elements will lock and one must use either a mixed method or high order elements. As we will see, when continuity is relaxed, a DG displacement-only formulation with low order elements does not lock. To see an example more clearly, consider a very poorly chosen mesh such as in Fig. 1. Due to the structure, a conforming method will produce asymmetric deformations even when symmetric loads are applied. However, when we relax the continuity constraint, we see that DG does not suffer this pathology. In addition to reporting some of the advantages, we also hope to refute some of the common thought disadvantages of DG for elliptic problems. There is often a predominant concern over the excessive costs stemming from extra degrees of freedom (i.e. the discontinuities). We will see that in some cases, the DG approximation to the derivative is slightly more accurate than a highly-refined conforming approximation with the same number of degrees of freedom. This means that a given level of strain accuracy may be just as expensive for a conforming method as it is for DG. Another concern involves the loops over interior element faces which adds a significant cost in three-dimensions. As we will see, DG can be formulated such that no loops over interior faces are needed to assemble the linearized elasticity problem. In this case, the assembly routine is virtually the same as that of a conforming method with a slight modification to the element matrices. Another common concern is the resulting illconditioned linearized problem. This could make it impossible to solve large three-dimensional problems in parallel. We found that the traditional Additive Schwarz preconditioner seemed (a) DG approximation (b) Conforming approximation Fig. 1. The reference configuration is a compressible neo-hookean square with unit length. The upper surface is translated down via Dirichlet boundary conditions with 100 increments. The left figure shows an intermediate configuration, the right figure shows the final configuration. For the DG method we set β = 0.01k where k is the increment number. to work reasonably well for a worst case scenario in three-dimensions where DG is used to model artery micromechanics. The linear system of equations has 4.4(10 6 ) unknowns and is extremely ill-conditioned since the material coefficients rapidly oscillate by a factor of Despite this rather arduous challenge, we are able to solve the linearized system in parallel using 343 processors with a total solve time of 30 minutes. For many supercomputing clusters, this is not an unreasonable request. Another concern is that most DG methods are sensitive to the method of stabilization which almost always is a penalty on discontinuities. Penalizing discontinuities can cause problems, since these are the very degrees of freedom in which many of DG s benefits arise. Several authors have shown that it can be easy to over penalize and substantially degrade the performance of the DG method. We will see that DG can be stabilized in such a way that the performance is virtually unaffected by the user s choice of penalty. Based on our research thus far, we see no reason to discount DG as a viable tool to solve real world nonlinear elasticity problems. 2
3 Let us begin by defining the nonlinear elasticity problem for a body that occupies a connected open domain B 0 R d, with d = 2, 3. The body deforms due to the presence of applied forces or displacements on the surface, B 0. The deformation of the reference configuration is described by the deformation mapping, ϕ : B 0 R d. The mechanical behavior of the material is assumed to be hyperelastic and defined by a strain energy density W : R d d R, W (F), which depends on the deformation gradient, F(X) = ϕ(x). Mechanical loads in the form of surface tractions T are applied on τ B 0 B 0. The potential energy of the body when deformed with a deformation mapping ϕ is computed as I[ϕ] = W (F) dv T ϕ ds. (1) B 0 τb 0 The nonlinear elasticity problem consists in finding the deformations in which the body is in mechanical equilibrium under the imposed boundary conditions. We seek these configurations within a space of admissible deformations V. The set of all admissible variations of this space consists of all functions of the form δϕ = ϕ ɛ / ɛ ɛ=0, where ϕ ɛ is a smooth one-parameter family of functions in V. An equilibrium configuration of the body is a deformation mapping ϕ V that renders the potential energy stationary with respect to all admissible variations. For the discrete problem, we begin by defining Vh, a discrete vector space of smooth scalar functions on an element. The DG finite element space is then V h = Π Vh, and hence the elements of V h are allowed to have discontinuities across element boundaries. For the discrete nonlinear elasticity problem we shall use two types of such spaces, V h and W h. The DG space Vh d will be used to approximate the deformation mapping ϕ, and the space Wh d d will be used to approximate its gradient (rather than Vh d ). In order to define the discrete approximation to the gradient we use the lifting operator R : L 2 (Γ) d Wh d d which requires R(ϕ) z dv = B 0 Γ ϕ N {z } ds, (2) Traction (Pa) e-04 1e-05 1e-06 xact-ν=0.5 Conforming-ν= DG-ν= Conforming-ν=0.499 DG-ν= r 0 /R 0 Fig. 2. Tractions measured on the inner wall of an expanding cylinder as a function of the ratio of the deformed to reference radius using 10,000 increments (1.1 represents a 10% increase in radius). Stabilization parameter β = 10 4 k where k is the increment number. As measured, the traction represents the internal pressure needed to expand the cylinder to a given radius. Notice that for the conforming elements the tractions are several orders of magnitude larger than necessary. However, for DG we see a good approximation to the exact solution. for all z Wh d d, where Γ is the set of all element boundaries and {ϕ } = 1 2 (ϕ+ + ϕ ). The calculation of the lifting operator is the only operation that involves a loop over element faces. However, this is a preprocessing step that must be done once for a given mesh. Since its calculation involves a local element mass matrix inversion (and not a global one) the cost is trivial compared to solving the linearized elasticity problem. Once the lifting operator is computed the DG derivative D DG : Vh d W d d, is simply h D DG ϕ = ϕ + R([ϕ]) (3) where we have used the jump operator, [ϕ] = ϕ ϕ +. For this particular DG method (often called the Bassi-Rebay method) we clearly see a distinction, in the discrete derivative, between DG and conforming. The DG derivative uses the additional degrees of freedom (i.e. the discontinuities) via the lifting operator which takes functions in Γ and returns an enhancement in Wh d d. This DG method yields a particularly simple formulation for the discrete linear elasticity problem. In 3
4 the discrete problem we again look for stationary points in the total potential energy, however now in the discontinuous space, Vh d. Here, we use the same total potential energy as shown in q. 1 now approximating ϕ with ϕ h, and F with D DG ϕ h. The resulting problem is nonlinear, with the nonlinearity embedded in the constitutive relation, W (D DG ϕ h ). It can be solved by taking a second variation, δ 2, and using Newton type iteration. In the absence of external tractions, the linearized problem becomes W F (D DG ϕ h) : δd DG ϕ h dv = δd DG ϕ h : 2 W F 2 (D DG ϕ h): δ 2 D DG ϕ h dv. (4) Most DG methods will display numerical instabilities without proper stabilization. In nonlinear elasticity, numerical instabilities can be an especially challenging task due to the loss of positive definiteness of the elasticity tensor, A = 2 W/ F 2, in selected regions of B 0. Typically, a stable method is achieved through a penalty on the discontinuities. Using the traditional approach one would add to the total potential energy a term such as β [ϕ h h ] 2 ds, (5) Γ where β > 0 and h is a measure of the mesh size. In Fig. 2 we see an example where traditional stabilization is used to approximate the deformations of a nearly incompressible neo-hookean material. Despite having to satisfy the nonlinear constraint, det(f) = 1, we see that DG does not lock even though a purely displacement formulation is used with linear triangles. 1 Deformation Gradient rror Uniformly stabilized Adaptive ε=0.0 Adaptive ε= Fig. 3. (Left) Highly refined DG approximation using adaptive stabilization on 125 processors with trilinear hexahedra. The material is a compressible neo-hookean. The DG approximation has h = 1/36 resulting in 1.1(10 6 ) dof. Colored regions show where A is indefinite and thus where discontinuities are penalized. (Right) Comparison of the deformation gradient errors for adaptive stabilization and traditional stabilization with uniform penalty. The deformation gradient error is measured as F h F exact 0,B0 where F exact is the exact deformation gradient, and F h is its numerical approximation. For DG, the deformation gradient approximation is F h = D DG ϕ h. The mesh refinement is fixed with h = 1/12, and the errors are scaled by the maximum error observed. For each method, a decrease in β beyond what is shown resulted in a loss of stability. As β increases the uniform method quickly becomes sensitive whereas the adaptive method with ɛ = 0 is virtually insensitive and with ɛ = 0.01 a moderate reduction in error is seen for small values of β Although, traditional stabilization may perform well, several authors have shown that many DG methods are rather sensitive to the penalization parameter, β which is user defined. The methods displayed large variations in the accuracy of the discrete approximation to both ϕ and F when β is changed by very little. In an attempt to automatically determine β and optimize the accuracy of the approximation for a given h, we propose a method of adaptive stabilization. In this method, we take advantage of the indefiniteness β 4
5 of A. For a point X B 0 we set λ X = λ(a(f(x))), with λ(a) = max { } g : A : g ɛ, min 0 g R d d g : g (6) where ɛ R + and ɛ << 1 so that if A(F(X)) is positive definite then λ X is nearly zero, and if A(F(X)) has one or more negative eigenvalues then λ X is nearly equal to the absolute value of the smallest eigenvalue of A(F(X)). The newly stabilized method is formulated in terms of an incremental variational principle. At an increment k + 1 we add a stabilization term of the form β 2 T h λ k X R([ϕ h ]) 2 dv (7) where β > 0 is the stabilization parameter and λ k X is computed from the deformation obtained at increment k. Several authors have proved the equivalence of this term to the one shown in q. 5. In this case, the method is guaranteed to be stable provided that β is large enough. However, we have found for this method a suitable value of β is generally O(1). 0.1 DG Linear conforming 1 DG Linear conforming Displacement rror Deformation Gradient rror 0.1 1e Total DOF Total DOF Fig. 4. Convergence plot for the displacements (left) and deformation gradient (right) as a function of the number of degrees of freedom. The reference configuration is a neo-hookean square with unit length. Here we use ɛ = 0 and β = 1 with linear triangles. The displacement error is ϕ h ϕ exact 0,B0, the deformation gradient error is F h F exact 0,B0 where ϕ exact and F exact are calculated from a highly-refined conforming solution. In this case, the linear conforming approximation is more efficient than the proposed method for the computation of displacements, while the situation is reversed in the case of deformation gradient. In Fig. 3 the accuracy in the approximation to the deformation gradient is compared for the adaptive DG method and one using a traditional stabilization mechanism of the form β B 0 R([ϕ h ]) 2 dv. The standard stabilization uniformly penalizes discontinuities on each of the faces in the mesh, while the adaptive stabilization adds additional penalization to those areas in B 0 where A is indefinite. In this study we evaluate the effect of β on the accuracy of the DG approximation. Here reference configuration is a neo-hookean cube with unit length. A highly refined DG approximation to the deformation is shown in Fig. 3. For this particular deformation, simulations were done in parallel. To solve the linear system of equations we used Additive Schwarz preconditioning. For the preconditioner, an overlap of exactly one element (h) was implemented, and a direct LU factorization was used to solve the local problem. The linear solver and preconditioner were implemented within the PTSc framework. In the figure we clearly see the robustness and improved accuracy of the adaptive method for both ɛ = 0 and ɛ = In Fig. 4 we compare the accuracy as a function of the number of degrees of freedom in the mesh for the deformation gradient and the displacements. A DG and conforming approximation to the final deformation is shown in Fig.5. We measure the computational cost as the total number of degrees of freedom in the discrete problem. In the figure we see that the DG approximation to the deformation gradient is competitive with the conforming one. To more clearly see the advantage of using Wh d d to approximate derivatives, we 5
6 have plotted contours of D DG ϕ h for one DG approximation and ϕ h for a conforming approximation with the same number of degrees of freedom. In Fig. 5 we see that Wh d d yields piecewise linear interpolation while Vh d yields piecewise constant interpolation. Fig. 5. A coarse DG approximation (h = 1/9, left) and a refined conforming approximation (h = 1/24, right). The contours show the F 11 component of the approximation to the deformation gradient. Notice that the DG deformation gradient are piecewise linear inside each element, in contrast to the piecewise constant approximations in the conforming case. Both solutions have nearly the same error in deformation gradient computed with nearly the same number of degrees of freedom, DG has 972 and conforming In Fig. 6 we use DG with adaptive stabilization to model the elastin structure within a small cubic sample of arterial tissue with length 25(10 6 )m. The elastin structure was constructed from images obtained by Mary K. O Connell and Charles A. Taylor in the Cardiovascular Biomechanics Research Laboratory at Stanford. Due to the accuracy of DG and its non-locking behavior it is an ideal candidate for investigating arterial micromechanics. A highly refined DG mesh with 4.4(10) 6 dof on 343 processors is shown in Fig. 6. Fig. 6. (Left) lastin (red) suspended in soft extracellular matrix (blue). Both the elastin and extracellular matrix are modeled as a nearly incompressible linear elastic solid. (Right) Contours of the shear stress are shown when the sample is subjected to a transverse displacement of the upper surface. We see the DG approximation is able to capture some of the localizations that occur due to the rapidly oscillating material coefficients. The ability to capture stress localizations will be key when investigating phenomena such as aneurysm. 6
Some applications of discontinuous Galerkin methods in solid mechanics
Some applications of discontinuous Galerkin methods in solid mechanics Adrian Lew, Alex Ten Eyck, and Ramsharan Rangarajan Mechanical Engineering, Stanford University, USA. Email: {lewa,alabute,rram}@stanford.edu
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 informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More 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 informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
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 informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
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 informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (E-post:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
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 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 informationConcepts. 3.1 Numerical Analysis. Chapter Numerical Analysis Scheme
Chapter 3 Concepts The objective of this work is to create a framework to implement multi-disciplinary finite element applications. Before starting, it is necessary to explain some basic concepts of the
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 informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
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 informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
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 informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationNumerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions
TECHNISCHE MECHANIK, 33, 2, (2013), 97 103 submitted: December 11, 2012 Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions R. Glüge, M. Weber
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
More informationSimulation of elasticity, biomechanics and virtual surgery. Joseph M. Teran
Contents Simulation of elasticity, biomechanics and virtual surgery Joseph M. Teran 1 Simulation of elasticity, biomechanics and virtual surgery 3 Introduction 3 Real-time computing 4 Continuum mechanics
More informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
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 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 methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More information2.1 Strain energy functions for incompressible materials
Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to
More informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationA Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements
A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationPerformance evaluation of different model mixers by numerical simulation
Journal of Food Engineering 71 (2005) 295 303 www.elsevier.com/locate/jfoodeng Performance evaluation of different model mixers by numerical simulation Chenxu Yu, Sundaram Gunasekaran * Food and Bioprocess
More informationPreprint Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction
Fakultät für Mathematik und Informatik Preprint 2015-15 Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction ISSN 1433-9307
More informationA note on finite elastic deformations of fibre-reinforced non-linearly elastic tubes
Arch. Mech., 67, 1, pp. 95 109, Warszawa 2015 Brief Note A note on finite elastic deformations of fibre-reinforced non-linearly elastic tubes M. EL HAMDAOUI, J. MERODIO Department of Continuum Mechanics
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationUsing Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses
Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses Brian Baillargeon, Ramesh Marrey, Randy Grishaber 1, and David B. Woyak 2 1 Cordis Corporation,
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 informationIntroduction to Finite Element computations
Non Linear Computational Mechanics Athens MP06/2012 Introduction to Finite Element computations Vincent Chiaruttini, Georges Cailletaud vincent.chiaruttini@onera.fr Outline Continuous to discrete problem
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 informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
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 information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More 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 note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationLecture 1 NONLINEAR ELASTICITY
Lecture 1 NONLINEAR ELASTICITY Soft-Matter Engineering: Mechanics, Design, & Modeling Mechanics Rigid/Hard Systems bending or torsion (but not stretching) small strains (~0.1% metals; ~1% plastics) linearized
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
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 informationFor an imposed stress history consisting of a rapidly applied step-function jump in
Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationin this web service Cambridge University Press
CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationCreasing Critical Strain Dependence on Surface Defect Geometry. EN234 Final Project
Creasing Critical Strain Dependence on Surface Defect Geometry EN234 Final Project A Landauer Dec 16, 2015 1 Problem Description In elastic soft homogeneous materials that admit large compressive deformations
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
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 informationNumerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35
Numerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35 Mirjam Yvonne van Leeuwen Supervisor: Dr. Ir. M.C.M. Rutten Ir. N.J.B. Driessen TUE Eindhoven, The Netherlands
More informationCoupled Thermomechanical Contact Problems
Coupled Thermomechanical Contact Problems Computational Modeling of Solidification Processes C. Agelet de Saracibar, M. Chiumenti, M. Cervera ETS Ingenieros de Caminos, Canales y Puertos, Barcelona, UPC
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationNonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec
Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling
More informationModel-independent approaches for the XFEM in fracture mechanics
Model-independent approaches for the XFEM in fracture mechanics Safdar Abbas 1 Alaskar Alizada 2 and Thomas-Peter Fries 2 1 Aachen Institute for Computational Engineering Science (AICES), RWTH Aachen University,
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 informationFinite Element Modeling and Analysis. CE 595: Course Part 2 Amit H. Varma
Finite Element Modeling and Analysis CE 595: Course Part 2 Amit H. Varma Discussion of planar elements Constant Strain Triangle (CST) - easiest and simplest finite element Displacement field in terms of
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationNonlinear analysis in ADINA Structures
Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent
More informationA FACE-BASED SMOOTHED FINITE ELEMENT METHOD FOR HYPERELASTIC MODELS AND TISSUE GROWTH
th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
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 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 informationMITOCW MITRES2_002S10nonlinear_lec05_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationComputing Periodic Orbits and their Bifurcations with Automatic Differentiation
Computing Periodic Orbits and their Bifurcations with Automatic Differentiation John Guckenheimer and Brian Meloon Mathematics Department, Ithaca, NY 14853 September 29, 1999 1 Introduction This paper
More informationDiscontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, 22 24 March 2016 Structure
More information