Node-to-segment and node-to-surface interface finite elements for fracture mechanics

The topologies of existing interface elements used to discretize cohesive cracks are such that they can be used to compute the relative displacements (displacement discontinuities) of two opposing segments (in 2D) or of two opposing facets (in 3D) belonging to the opposite crack faces and enforce the cohesive traction-separation relation. In the present work we propose a novel type of interface element for fracture mechanics sharing some analogies with the node-to-segment (in 2D) and with the node-to-surface (in 3D) contact elements. The displacement gap of a node belonging to the finite element discretization of one crack face with respect to its projected point on the opposite face is used to determine the cohesive tractions, the residual vector and its consistent linearization for an implicit solution scheme. The following advantages with respect to classical interface finite elements are demonstrated: (i) non-matching finite element discretizations of the opposite crack faces is possible; (ii) easy modelling of cohesive cracks with non-propagating crack tips; (iii) the internal rotational equilibrium of the interface element is assured. Detailed examples are provided to show the usefulness of the proposed approach in nonlinear fracture mechanics problems.Comment: 37 pages, 17 figure

Hybrid coupling of CG and HDG discretizations based on Nitsche’s method

This is a post-peer-review, pre-copyedit version of an article published in Computational mechanics. The final authenticated version is available online at: http://dx.doi.org/10.1007/s00466-019-01770-8A strategy to couple continuous Galerkin (CG) and hybridizable discontinuous Galerkin (HDG) discretizations based only on the HDG hybrid variable is presented for linear thermal and elastic problems. The hybrid CG-HDG coupling exploits the definition of the numerical flux and the trace of the solution on the mesh faces to impose the transmission conditions between the CG and HDG subdomains. The con- tinuity of the solution is imposed in the CG problem via Nitsche’s method, whereas the equilibrium of the flux at the interface is naturally enforced as a Neumann con- dition in the HDG global problem. The proposed strategy does not affect the core structure of CG and HDG discretizations. In fact, the resulting formulation leads to a minimally-intrusive coupling, suitable to be integrated in existing CG and HDG libraries. Numerical experiments in two and three dimensions show optimal global convergence of the stress and superconvergence of the displacement field, locking-free approximation, as well as the potential to treat structural problems of engineering interest featuring multiple materials with compressible and nearly incompressible behaviors.Peer ReviewedPostprint (author's final draft

Stabilized MorteX method for mesh tying along embedded interfaces

We present a unified framework to tie overlapping meshes in solid mechanics applications. This framework is a combination of the X-FEM method and the mortar method, which uses Lagrange multipliers to fulfill the tying constraints. As known, mixed formulations are prone to mesh locking which manifests itself by the emergence of spurious oscillations in the vicinity of the tying interface. To overcome this inherent difficulty, we suggest a new coarse-grained interpolation of Lagrange multipliers. This technique consists in selective assignment of Lagrange multipliers on nodes of the mortar side and in non-local interpolation of the associated traction field. The optimal choice of the coarse-graining spacing is guided solely by the mesh-density contrast between the mesh of the mortar side and the number of blending elements of the host mesh. The method is tested on two patch tests (compression and bending) for different interpolations and element types as well as for different material and mesh contrasts. The optimal mesh convergence and removal of spurious oscillations is also demonstrated on the Eshelby inclusion problem for high contrasts of inclusion/matrix materials. Few additional examples confirm the performance of the elaborated framework

Singular function mortar finite element methods

This is the published version, also available here: http://dx.doi.org/10.2478/cmam-2003-0014.We consider the Poisson equation with Dirichlet boundary conditions on a polygonal domain with one reentrant corner. We introduce new nonconforming finite element discretizations based on mortar techniques and singular functions. The main idea introduced in this paper is the replacement of cut-off functions by mortar element techniques on the boundary of the domain. As advantages, the new discretizations do not require costly numerical integrations and have smaller a priori error estimates and condition numbers. Based on such an approach, we prove optimal accuracy error bounds for the discrete solution. Based on such techniques, we also derive new extraction formulas for the stress intensive factor. We establish optimal accuracy for the computed stress intensive factor. Numerical examples are presented to support our theory

Isogeometric dual mortar methods for computational contact mechanics

International audienceIn recent years, isogeometric analysis (IGA) has received great attention in many fields of computational mechanics research. Especially for computational contact mechanics, an exact and smooth surface representation is highly desirable. As a consequence, many well-known finite e lement m ethods a nd a lgorithms f or c ontact m echanics h ave b een t ransferred t o I GA. I n t he present contribution, the so-called dual mortar method is investigated for both contact mechanics and classical domain decomposition using NURBS basis functions. In contrast to standard mortar methods, the use of dual basis functions for the Lagrange multiplier based on the mathematical concept of biorthogonality enables an easy elimination of the additional Lagrange multiplier degrees of freedom from the global system. This condensed system is smaller in size, and no longer of saddle point type but positive definite. A very simple and commonly used element-wise construction of the dual basis functions is directly transferred to the IGA case. The resulting Lagrange multiplier interpolation satisfies discrete inf–sup stability and biorthogonality, however, the reproduction order is limited to one. In the domain decomposition case, this results in a limitation of the spatial convergence order to O(h 3 /2) in the energy norm, whereas for unilateral contact, due to the lower regularity of the solution, optimal convergence rates are still met. Numerical examples are presented that illustrate these theoretical considerations on convergence rates and compare the newly developed isogeometric dual mortar contact formulation with its standard mortar counterpart as well as classical finite elements based on first and second order Lagrange polynomials