Kaczmarczyk, L.

Pearce, C. J.

Ullah, Zahur

English

The mortar contact formulation is a well-established technique to tie non-conforming finite element meshes in domain decomposition and is also the basis of many well-known contact algorithms. Mortar contact formula- tion allows for a variationally consistent treatment of contact conditions including mesh tying, non-penetration, frictionless and frictional sliding leading to satisfaction of contact patch test. Efficient, accurate and robust nu- merical implementation of the interface coupling terms associated with the mortar contact formulation remains challenging, especially in three-dimensional case. The computational contact algorithm presented in this paper is carefully designed for accuracy, efficiency and robustness and making use of the cutting-edge third-party computa- tional tools including Mesh-Oriented datABase (MOAB), Portable, Extensible Toolkit for Scientific Computation (PETSc), Boost and clipper libraries. The computational framework is designed to take advantage of distributed memory high-performance computing and hierarchic basis functions. The numerical implementation is validated with two non-conforming mesh tying examples, which, on the one hand, remove some of the complexities as- sociated with actual unilateral contact formulation but, on the other hand, clarify many of the conceptual and implementational aspects of the contact mechanics

Ulster University's Research Portal

Proceedings of the 25th UKACM Conference on Computational Mechanics12 - 13 April 2017, University of BirminghamBirmingham, United KingdomTHREE-DIMENSIONAL MORTAR CONTACT FORMULATION:AN EFFICIENT AND ACCURATE NUMERICAL IMPLEMENTATION*Zahur Ullah1, Łukasz Kaczmarczyk1 and Chris J. Pearce11School of Engineering, Rankine Building, The University of Glasgow, Glasgow, UK, G12 8LT*Zahur.Ullah@glasgow.ac.ukABSTRACTThe mortar contact formulation is a well-established technique to tie non-conforming finite element meshes indomain decomposition and is also the basis of many well-known contact algorithms. Mortar contact formula-tion allows for a variationally consistent treatment of contact conditions including mesh tying, non-penetration,frictionless and frictional sliding leading to satisfaction of contact patch test. Efficient, accurate and robust nu-merical implementation of the interface coupling terms associated with the mortar contact formulation remainschallenging, especially in three-dimensional case. The computational contact algorithm presented in this paper iscarefully designed for accuracy, efficiency and robustness and making use of the cutting-edge third-party computa-tional tools including Mesh-Oriented datABase (MOAB), Portable, Extensible Toolkit for Scientific Computation(PETSc), Boost and clipper libraries. The computational framework is designed to take advantage of distributedmemory high-performance computing and hierarchic basis functions. The numerical implementation is validatedwith two non-conforming mesh tying examples, which, on the one hand, remove some of the complexities as-sociated with actual unilateral contact formulation but, on the other hand, clarify many of the conceptual andimplementational aspects of the contact mechanics.Key Words: finite element analysis; mortar contact formulation; mesh tying; numerical integration; hierarchicalbasis functions1. IntroductionMechanical interaction between different bodies often referred as "contact or impact", is of great impor-tance in many engineering applications including prosthetics in biomedical engineering, pneumatic tiresin automotive engineering and adhesion or slip between concrete and reinforcing steel in civil engineer-ing. In mechanical engineering applications, contact can be found in gears, bearings, metal forming andcar crash test [1, 2]. At the minimum, contact mechanics involves searching for the contact area betweeninteracting bodies and subsequent prevention of inter-penetration. Relative movement or slip includ-ing both frictionless and frictional, intermittent interaction and wear also comes under the umbrella ofcontact mechanics. Due to the associated nonlinearities and complicated nature of problems involvingcontact, special assumptions were used in the past for their solution. Advances in Computational mod-elling allows solving these problems to be solved numerically with sufficient accuracy for engineeringanalysis/design. On the other hand, as compared to the current state of the art nonlinear finite elementtechnology including finite deformation kinematics, inelastic material behavior and linear and nonlinearequations solvers, contact mechanics is relatively immature. Therefore, the design of efficient, accurateand robust contact algorithms is still a challenge for practical engineering problems.In this paper, one of the most commonly used contact discretisation technique for the solution of contactproblem, i.e. the mortar contact formulation is adopted, which is also referred to as a special type ofsegment-to-segment approach in the literature [3]. Alternative methods, including node-to-segment andenforcing contact condition at specific finite element node, cannot guarantee the satisfaction of contactpatch test and suffers from non-physical oscillation in the contact forces. As compared to the strongor point-wise satisfaction of the interface continuity condition in the aforementioned contact mechan-ics approaches, a weak or integral form is used in the mortar contact formulation. The computationalcontact algorithm presented in this paper is carefully designed for accuracy, efficiency and robustnessand makes use of the cutting-edge third party computational tools including Mesh-Oriented datABase(MOAB), Portable, Extensible Toolkit for Scientific Computation (PETSc), Boost and clipper [4] li-braries. The developed algorithm is implemented within our group’s finite element code, MOFEM [5].Relatively simple, linear-elastic two body problems with non-conforming meshes are used to validatethe numerical implementation. Tetrahedral elements are used to discretised the two contacting bodiesleading to triangular elements on the common interfaces. The clipper library is used to create an inter-secting polygon between every pair of triangles belongs to the two contacting surfaces. These polygonsare then triangulated, with their integration points projected back on to the parent triangles, which areused subsequently for the numerical integration of the mortar contact integral. For the convenient andefficient handling of the mesh data associated with the contacting surfaces, prisms are inserted betweenevery pair of contacting triangles and are stored in the multi-index containers. Furthermore, the compu-tational framework is designed to take advantage of distributed memory high-performance computingand hierarchic basis functions [6].2. Problem formulationIn this paper, we limit ourselves to the application of mortar contact formulation for tying non-conforming meshes in linear-elasticity by removing some of the complexities associated with unilateralcontact mechanics. Nonetheless, this will clarify many of the conceptual and implementational aspectsof the computational contact mechanics [7]. Consider two subdomains Ω(i) ⊂ R3, i = 1,2 bounded byboundaries ∂Ω(i). The subdomain boundary ∂Ω(i) is divided into disjoint sets of Dirichlet boundary Γ(i)u ,Neumann boundary Γ(i)σ and mesh tying interface Γc . The body force acting over the individual subdo-main is b(i). The strong form of governing equations in the case of linear-elasticity for each subdomainis written as:Divσ (i) + b(i) = 0 in Ω(i), (1)where σ (i) is the Cauchy stress tensor. The associated boundary conditions are the Dirichlet, Neumannand mesh tying constraint and are written as:u(i) = u(i) on Γ(i)u , (2a)σ (i) · n(i) = t(i) on Γ(i)σ , (2b)u(1) = u(2) on Γc . (2c)The weak formulation is subsequently derived together with Lagrange multipliers λ to enforce the meshtying constraint. In Figure 1(a), Ω(2) is considered as the slave side and Ω(1) is considered as the masterside. The final discretised set of equations is written as:KNN KNM KNS 0KMN KMM 0 −MTKSN 0 KSS DT0 −M D 0dNdMdSλ=FNFMFS0, (3)where dS and dM are the degrees of freedom associated with slave and master contact surfaces and dNconsists of all the remaining degrees of freedom. K and F are the standard stiffness matrices and forcevectors respectively. Both slave and master surfaces are not connected directly and are tied using themortar contact formulation leading to blocks of zeros for both KMS and KSM . Constraint matrices Dand M are written as:D =∫ΓcNλNSdΓc =npr∑pr=1*.,ngp∑g=1wSg Nλ(ξS (ξg ), ηS (ηg ))NS(ξS (ξg ), ηS (ηg ))JS+/- , (4a)M =∫ΓcNλNMdΓc =npr∑pr=1*.,ngp∑g=1wMg Nλ(ξS (ξg ), ηS (ηg ))NM(ξM (ξg ), ηM (ηg ))JM+/- , (4b)(a) Two body mesh tying problem in 3Dclip polygon(b) Triangles in 3D space (c) Projected triangles and      clip polygon(d) Triangulation and      integration points(e) Projection back on     both original trianglesFigure 1: Mesh tying problem and steps involved in 3D mortar contact formulation for one pair of master andslave triangular elementswhere NS and NM are matrices of shape functions for slave and master sides respectively. Nλ is amatrix of shape functions used for the discretisation of Lagrange multipliers, only exists on the slaveside and is assumed the same as NS . w and J are the weight and Jacobian associated with Gauss points.(ξS (ξg ), ηS (ηg ))and(ξM (ξg ), ηM (ηg ))represents the projection of Gauss points on the original slaveand master triangles in local coordinates, explained below in detail.A step-by-step procedure used for the 3D mortar contact formulation for one pair of slave and mastertriangles is shown in Figures 1(b-e). After projection of both triangles in xy-plane Figure (1(b)), theclipper library is used to determine the clip polygon (Figure 1(c)), which is subsequently triangulated(Figure 1(d)) with their integration points projected back on both the original master and slave triangles(Figure 1(e)). For the easy and efficient data handling, prisms are inserted between each intersectingslave and master triangles. In Equation (4), npr and ngp are the number of prisms and Gauss pointsrespectively.3. Numerical ExampleThe two pairs of both non-conforming and conforming meshes, as shown in Figure 2, are used to validatethe numerical implementation of the mortar contact formulation. The dimensions, loading, boundaryconditions and coordinate system are also shown in Figure 2. Young’s modulus and Poisson’s ratio usedin this case are E = 10 and ν = 0.5 respectively. For both pairs, the conforming meshes, i.e. casesshown in Figures 2(b, d), are used for validation and is solved without mortar contact formulation. InFigure 2(a) and Figure 2(c) non-conforming meshes exists at their curved interfaces. The contours ofz-component of displacement, i.e. uz is also shown for all of the four cases. It can be seen that foreach pair, the contours plots of uz is exactly the same, demonstrating the correct implementation of thecomputational framework.4. ConclusionsIn this study, an efficient, accurate and robust numerical framework is presented for three-dimensionalmortar contact formulation, which allows consistent treatment of the contact conditions. The compu-tational framework is implemented in our group’s finite element code, MOFEM, which uses state of(a) Non-conforming meshes(b) Conforming meshes(c) Non-conforming meshes(d) Conforming meshes102060301010Figure 2: Two pairs of non-conforming and conforming meshes and corresponding contours of uzthe art MOAB, PETSc and Boost libraries. The implementation is validated with two non-conformingmesh tying examples. Although, not demonstrated in this paper, the developed computational frame-work is designed to take advantages of the hierarchic basis functions and high-performance computing.This paper is restricted to linear-elastic problems but will be extended subsequently to include bothmaterial and geometric nonlinearities. The actual contact formulation including efficient contact search,non-penetration, frictionless and frictional sliding will be implemented next.AcknowledgementsThis work was supported by EDF Energy Nuclear Generation Ltd and the Royal Academy of Engineer-ing (RAEng). The views expressed in this paper are those of the authors and not necessarily those ofEDF Energy Nuclear Generation Ltd or RAEng.References[1] P. Wriggers. Computational Contact Mechanics. Springer Berlin Heidelberg, 2006.[2] T. A. Laursen. Computational Contact and Impact Mechanics: fundamentals of modeling interfacialphenomena in nonlinear finite element analysis. Springer Science & Business Media, 2013.[3] A. Popp, M. W. Gee, and W. A. Wall. A finite deformation mortar contact formulation usinga primal–dual active set strategy. International Journal for Numerical Methods in Engineering,79(11):1354–1391, 2009.[4] A. Johnson. Clipper - an open source freeware library for clipping and offsetting lines and polygons.http://www.angusj.com/delphi/clipper.php, 2014.[5] Mesh Oriented Finite Element Method (MOFEM) Version 0.5.38, University of Glasgow, Glasgow,UK. http://mofem.eng.gla.ac.uk/mofem/html/, 2017.[6] M. Ainsworth and J. Coyle. Hierarchic finite element bases on unstructured tetrahedral meshes.International Journal for Numerical Methods in Engineering, 58(14):2103–2130, 2003.[7] A. Popp. Mortar Methods for Computational Contact Mechanics and General Interface Problems.PhD thesis, Technical University of Munich, 2012.

Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation

The mortar contact formulation is a well-established technique to tie non-conforming finite element meshes in domain decomposition and is also the basis of many well-known contact algorithms. Mortar contact formulation allows for a variationally consistent treatment of contact conditions including mesh tying, non-penetration, frictionless and frictional sliding leading to satisfaction of contact patch test. Efficient, accurate and robust numerical implementation of the interface coupling terms associated with the mortar contact formulation remains challenging, especially in three-dimensional case. The computational contact algorithm presented in this paper is carefully designed for accuracy, efficiency and robustness and making use of the cutting-edge third-party computational tools including Mesh-Oriented datABase (MOAB), Portable, Extensible Toolkit for Scientific Computation (PETSc), Boost and clipper libraries. The computational framework is designed to take advantage of distributed memory high-performance computing and hierarchic basis functions. The numerical implementation is validated with two non-conforming mesh tying examples, which, on the one hand, remove some of the complexities associated with actual unilateral contact formulation but, on the other hand, clarify many of the conceptual and implementational aspects of the contact mechanics

Kaczmarczyk, Lukasz

Pearce, Chris J

Queen's University Belfast Research Portal

https://pure.ulster.ac.uk/ws/files/11630584/UKACM2017-Zahur.pdf

Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation

Abstract

Similar works

Full text

Available Versions

Ulster University's Research Portal

Queen's University Belfast Research Portal