The Plane Stress Problem
|
|
- Marjorie Grant
- 5 years ago
- Views:
Transcription
1 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, / 24
2 Outline Plane stress Discrete interpolation Triangular elements Rectangular elements The algebraic problem Solution, postprocessing Accuracy of FEM Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
3 Plane stress A 2D problem: plane stress Assumptions: thin specimen free to move normal to the plane (x 3 direction, i.e. direction of unit vector e 3 ) only loaded in the plane (x 1, x 2 direction) Image from Thus no stress normal to the plane, but specimen typically bulges ([sv. bukta]) in the normal direction so displacements is typically not zero normal to the plane (u 3 0). Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
4 Plane stress Plane stress modeling No stress (force) normal to the plane σ e 3 = 0 third column of σ is zero, and, because of σ = σ T, also the third row σ 11 σ 12 0 ( ) σ = σ 12 σ 22 0 ˆσ 0 = T. 0 Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
5 Plane stress Plane stress modeling No stress (force) normal to the plane σ e 3 = 0 third column of σ is zero, and, because of σ = σ T, also the third row σ 11 σ 12 0 ( ) σ = σ 12 σ 22 0 ˆσ 0 = T. 0 Constitutive law: σ = λi tr ɛ + 2µɛ = E [ 1 + ν ( ) ˆɛ 0 ɛ = 0 T ɛ 33, ν 1 2ν I tr ɛ ɛ, where ɛ 33 = ν tr ˆɛ. 1 ν ] Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
6 Plane stress Plane stress modeling Express constitutive law in terms of upper 2 2 block, ˆσ = E [ ] ν I tr ˆɛ + ˆɛ or 1 + ν 1 ν ˆσ = ˆλI tr ˆɛ + 2µˆɛ, ˆλ Eν = (1 + ν)(1 ν) Conclusion: for plane stress, compute in 2D with modified constitutive law Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
7 Plane stress Plane strain Adjoint/dual problem to plane stress (observe: not the mathematical adjointness ) Assumptions: no displacements orthogonal to the plane (a thick specimen or a restricted movement) loads constant in direction orthogonal to the plane Assumptions ɛ 11 ɛ 12 0 ɛ = ɛ 12 ɛ 12 0 = ( ) ˆɛ 0 0 T. 0 Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
8 Finite Element Discretization for Plane Stress Problems Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
9 Finite Element (FE) discretiztation for 2D linear elasticity Consider the equilibrium equation ) ) ] [ˆλ tr (ˆɛ(v) tr (ˆɛ(u) + 2µˆɛ(v) : ˆɛ(u) dω = Ω or, in abbreviated form, Ω v f dω + v g ds, Γ 1 a(v, u) = l(v) for all admissible v Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
10 Finite Element (FE) discretiztation for 2D linear elasticity Consider the equilibrium equation ) ) ] [ˆλ tr (ˆɛ(v) tr (ˆɛ(u) + 2µˆɛ(v) : ˆɛ(u) dω = Ω or, in abbreviated form, Ω v f dω + v g ds, Γ 1 a(v, u) = l(v) for all admissible v FEM aims at satisfying this equation for a subset of admissible displacements work on a discretization associated with a triangulation the domain, consisting of nonoverlapping triangles (quadrilaterals also common) An example triangulation Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
11 Discrete interpolation Discrete interpolation Admissible displacements u h are continuous on Ω and polynomials on each triangle (or quadrilateral). An example of one component of an admissible displacement when u h is a linear function (first-order polynomial) at each triangle. Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
12 Discrete interpolation Discrete interpolation Admissible displacements u h are continuous on Ω and polynomials on each triangle (or quadrilateral). An example of one component of an admissible displacement when u h is a linear function (first-order polynomial) at each triangle. For piecewise-linear functions, the nodes are the vertices of the triangulations lowest-order linear Lagrange element or the three-node triangle Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
13 Discrete interpolation Each admissible displacement u h is expressed as a weighted sum of basis functions φ j (x): N u h = φ j (x)u j. j=1 Basis function for continuous, piecewise-linear functions The sum interpolates discrete displacement vectors u j located at nodes (interpolation points) A nodal basis function is a continuous function in the finite element space that is one at one node, and zero at all other nodes Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
14 Discrete interpolation Lagrange elements P p Interpolation properties of nodal basis functions ( ): Polynomials of degree p on each triangle (element) Continuous across each edge Linear: P 1 The three-node triangle u(x, y) = a 1 + a 2x + a 3y Quadratic: P 2 The six-node triangle u(x, y) = a 1 + a 2x + a 3y + a 4x 2 + a 5xy + a 6y 2 Cubic; P 3 The ten-node triangle Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
15 Discrete interpolation Example basis functions Functions within a single element: For quadratic Lagrangian elements, there are two distinct types of nodal basis functions: Associated with mesh vertex Associated with edge midpoints Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
16 Discrete interpolation Example basis functions II Basis functions on their support: The basis function φ i(x) when i corresponds to a corner node The basis function φ i(x) when i corresponds to an edge-midpoint node Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
17 Discrete interpolation Rectangular elements A common element type for 2D elasticity. The displacements at the nodes ( ) are interpolated by a product of one-dimensional polynomials in each coordinate direction The four-node rectangular element Q 1. Bilinear, 4 coefficients u(x, y) = a 1 + a 2 x + a 3 y + a 4 xy (obtained as product of φ(x) = s 1 +s 2 x and ψ(y) = t 1 + t 2 y) Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
18 Discrete interpolation Rectangular elements, II The nine-node rectangular element Q 2. Biquadratic, 9 coefficients u(x, y) = a 1 + a 2 x + a 3 y + a 4 xy + a 5 x 2 + a 6 y 2 + a 7 x 2 y + a 8 xy 2 + a 9 x 2 y 2 Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
19 Discrete interpolation The four-node Q 1 on general quadrilaterals: Linear on each edge Not in general of form q(x, y) = a 0 + a 1 x + a 2 y + a 3 xy! A quadrilateral is not the image of an affine map of a rectangle! Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
20 The algebraic problem The algebraic problem Problem (FE equation for linear elasticity and plane stress) Find displacement field u h, such that a(v h, u h ) = l(v h ) for all v h V h. (2) Form of an admissible displacement in the finite element subspace: N N u h (x) = φ j (x)u j = φ j (x)(u 1,j e 1 + u 2,j e 2 ) = 2N j=1 j=1 j=1 where, for j = 1,..., N, φ j (x) u j (3) u 2j 1 = u 1,j u 2j = u 2,j φ 2j 1 (x) = e 1 φ j (x) φ 2j (x) = e 2 φ j (x). Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
21 The algebraic problem Substitute (3) into FE equation and test with v h = φ k, k = 1,..., 2N. Hence, (2) is equivalent to requiring 2N j=1 a(φ k, φ j ) u j = l(φ k ) for k = 1,..., 2N, which is a linear system Au = b, where [ ( A kj = a(φ k, φ j ) = λ tr ɛ(φk ) ) tr ( ɛ(φ j ) ) + 2µɛ(φ k ) ɛ(φ j ) ] dω, Ω b k = l(φ k ) = φ k f dω + φ k g ds, Ω Γ 1 u = (u 1, u 2,..., u 2N ) T = (u 1,1, u 2,1, u 1,2, u 2,2,..., u 1,N, u 2,N ) T Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
22 The algebraic problem Practical implementation of elastic equations Assembly of matrix A by a loop over all cells; deal.ii open source finite element library based on C++, authors W. Bangerth, R. Hartmann, G. Kanschat, tutorial step-8 for (CellIterator cell = cell_begin; cell!= cell_end; ++cell) { phi.reinit (cell); cell_matrix = 0; for (int q_point = 0; q_point < n_q_points; ++q_point) for (int i = 0; i < dofs_per_cell; ++i) for (int j = 0; j < dofs_per_cell; ++j) cell_matrix(i,j) += (lambda * trace(phi[u].gradient(i,q_point)) * trace(phi[u].gradient(j,q_point)) + 2 * mu * phi[u].symmetric_gradient(i,q_point) * phi[u].symmetric_gradient(j,q_point) ) * phi.jxw(q_point); cell->get_dof_indices(dof_indices); global_matrix.add (dof_indices, cell_matrix); } Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
23 Solution, postprocessing Solution and postprocessing Solve linear system Au = b Direct: Gauss/Cholesky decomposition, COMSOL: UMFPACK (sparse direct solver) Iteratively: for larger size (> degrees of freedom in 2D, > in 3D) CG with preconditioning, multigrid Postprocessing: Translate components in u to displacements at certain points Get stresses and/or strains from derivatives of (3), using kinematic and constitutive relation Computer lab Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
24 Accuracy of FEM Accuracy of FEM On the accuracy of the discrete FEM solutions Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
25 Accuracy of FEM Error estimates Discretization errors measured in integral norms. If everything is nice, we have for triangular or rectangular elements ( u u h L = 2 (Ω) Ω ) u u h 2 1/2 dω Ch p+1. h: largest edge in mesh; p: element polynomial order; C depends on the (p + 1)-th derivatives of u. Estimate requires some niceness conditions, typically Nondegenerate mesh refinements: r outer /r inner bounded as h 0 (or that the minimum (or maximum) angle is bounded away from 0 (or 180 ) Smooth solutions. Accuracy typically reduced in vicinity of reentrant corners on the boundary (ω > 180 ) Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24 r inner r outer!
26 Accuracy of FEM Remarks Stresses are functions of derivatives of the displacement (constitutive & kinematic relation) Therefore: accuracy of stress approximations typically one order lower than the displacement approximations In particular: stresses are piecewise constant for the P 1 element (3-node triangle) poor stress approximation Rule of thumb: P 1 too inaccurate (P 2 much better), see computer lab, task 1.1. Adaptive grid refinement (e.g. around reentrant corners) and error estimates central aspects in CSM Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
27 Accuracy of FEM P 1 versus Q 1 Advantage Q 1 Q 1 (4-node rectangle) performs usually much better than P 1 (both have error term h 2!) Reason: the presence of the quadratic xy term for Q 1 (stresses not constant in element) Rectangular elements can give higher accuracy for the same number of degrees of freedom, particularly in 3D Disadvantage quadrilaterals Mesh generation more difficult and less automatic for rectangular meshes compared to triangles Martin Kronbichler (TDB) The Plane Stress Problem February 2, / 24
The Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
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 informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
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 informationShape Function Generation and Requirements
Shape Function Generation and Requirements Requirements Requirements (A) Interpolation condition. Takes a unit value at node i, and is zero at all other nodes. Requirements (B) Local support condition.
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
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 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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
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 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 informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
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 informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationPreconditioned space-time boundary element methods for the heat equation
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 Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
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 informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in 2-D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf 2 Matthias K. Gobbert 2 1 Department of Mechanical Engineering 2 Department of Mathematics and Statistics
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 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 informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
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 informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationParallel Discontinuous Galerkin Method
Parallel Discontinuous Galerkin Method Yin Ki, NG The Chinese University of Hong Kong Aug 5, 2015 Mentors: Dr. Ohannes Karakashian, Dr. Kwai Wong Overview Project Goal Implement parallelization on Discontinuous
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
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 informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
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 informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationBuilding a finite element program in MATLAB Linear elements in 1d and 2d
Building a finite element program in MATLAB Linear elements in 1d and 2d D. Peschka TU Berlin Supplemental material for the course Numerische Mathematik 2 für Ingenieure at the Technical University Berlin,
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 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 informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
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 informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
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 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 informationFormulation of the displacement-based Finite Element Method and General Convergence Results
Formulation of the displacement-based Finite Element Method and General Convergence Results z Basics of Elasticity Theory strain e: measure of relative distortions u r r' y for small displacements : x
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationOverlapping Domain Decomposition and Multigrid Methods for Inverse Problems
Contemporary Mathematics Volume 8, 998 B -88-988--35- Overlapping Domain Decomposition and Multigrid Methods for Inverse Problems Xue-Cheng Tai, Johnny Frøyen, Magne S. Espedal, and Tony F. Chan. Introduction
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 informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
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 informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationIV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80
www..com www..com Code No: M0322/R07 Set No. 1 IV B.Tech. I Semester Supplementary Examinations, February/March - 2011 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
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 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 informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
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 information