An X‐FEM approach for large sliding contact along discontinuities

The extended finite element method (X‐FEM) has been developed to minimize requirements on the mesh design in a problem with a displacement discontinuity. This advantage, however, still remains limited to the small deformation hypothesis when considering sliding discontinuities. The approach presented in this paper proposes to couple X‐FEM with a Lagrangian large sliding frictionless contact algorithm. A new hybrid X‐FEM contact element was developed with a contact search algorithm allowing for an update of contacting surfaces pairing. The stability of the contact formulation is ensured by an algorithm for fulfilling Ladyzhenskaya‐Babuska‐Brezzi (LBB) condition. Several 2D simple examples are presented in this paper in order to prove its efficiency and stability. Copyright © 2008 John Wiley & Sons, Ltd.     

Vertex‐to‐face contact searching algorithm for three‐dimensional frictionless contact problems

This article presents a new vertex‐to‐face contact searching algorithm for the three‐dimensional (3‐D) discontinuous deformation analysis (DDA). In this algorithm, topology is applied to the contact rule when any two polyhedrons are close to each other. Attempt is made to expand the original contact searching algorithm from two‐dimensional (2‐D) to 3‐D DDA. Examples are provided to demonstrate the new contact rule for vertex‐to‐face contacts between two polyhedrons with planar boundaries. Copyright © 2005 John Wiley & Sons, Ltd.     

Extension of the node‐to‐segment contact element for surface‐expansion‐dependent contact laws

A class of friction laws depending on the measure of contact surface expansion is defined in the paper within the continuum contact mechanics framework. The nominal and spatial forms of constitutive relations are discussed, including incremental penalty relations. Further, an extended node‐to‐segment element is derived which is capable of treating surface‐expansion‐dependent contact laws in a consistent way. The approach is suitable for any kind of node‐to‐segment contact elements. Finally, the computational efficiency of the extended element as well as other possible approaches are illustrated by numerical examples relevant to metal forming applications. Copyright © 2001 John Wiley & Sons, Ltd.     

Modelling of reinforced‐concrete structures providing crack‐spacing based on X‐FEM,ED‐FEM and novel operator split solution procedure

In this work, we present a novel approach to the finite element modelling of reinforced‐concrete (RC) structures that provides the details of the constitutive behavior of each constituent (concrete, steel and bond‐slip), while keeping formally the same appearance as the classical finite element model. Each component constitutive behavior can be brought to fully non‐linear range, where we can consider cracking (or localized failure) of concrete, the plastic yielding and failure of steel bars and bond‐slip at concrete steel interface accounting for confining pressure effects. The standard finite element code architecture is preserved by using embedded discontinuity (ED‐FEM) and extended (X‐FEM) finite element strain representation for concrete and slip, respectively, along with the operator split solution method that separates the problem into computing the deformations of RC (with frozen slip) and the current value of the bond‐slip. Several numerical examples are presented in order to illustrate very satisfying performance of the proposed methodology. Copyright © 2010 John Wiley & Sons, Ltd.     

A local multigrid X‐FEM strategy for 3‐D crack propagation

A multiscale method for 3‐D crack propagation simulation in large structures is proposed. The method is based on the extended finite element method (X‐FEM). The asymptotic behavior of the crack front is accurately modeled using enriched elements and no remeshing is required during crack propagation. However, the different scales involved in fracture mechanics problems can differ by several orders of magnitude and industrial meshes are usually not designed to account for small cracks. Enrichments are therefore useless if the crack is too small compared with the element size. To overcome this drawback, a project combining different numerical techniques was started. The first step was the implementation of a global multigrid algorithm within the X‐FEM framework and was presented in a previous paper (Eur. J. Comput. Mech. 2007; 16 :161–182). This work emphasized the high efficiency in cpu time but highlighted that mesh refinement is required on localized areas only (cracks, inclusions, steep gradient zones). This paper aims at linking the different scales by using a local multigrid approach. The coupling of this technique with the X‐FEM is described and computational aspects dealing with intergrid operators, optimal multiscale enrichment strategy and level sets are pointed out. Examples illustrating the accuracy and efficiency of the method are given. Copyright © 2008 John Wiley & Sons, Ltd.     

Stabilized global–local X‐FEM for 3D non‐planar frictional crack using relevant meshes

A stabilized global–local quasi‐static contact algorithm for 3D non‐planar frictional crack is presented in the X‐FEM/level set framework. A three‐field weak formulation is considered and allows an independent discretization of the bulk and the crack interface. Then, a fine discretization of the interface can be defined according to the possible complex contact state along the crack faces independently from the mesh in the bulk. Furthermore, an efficient stabilized non‐linear LATIN solver dedicated to contact and friction is proposed. It allows solving in a unified framework frictionless and frictional contact at the crack interface with a symmetric formulation, no iterations on the local stage (unilateral contact law with/without friction), no calculation of any global tangent operator, and improved convergence rate. 2D and 3D patch tests are presented to illustrate the relevance of the proposed model and an actual 3D frictional crack problem under cyclic fretting loading is modeled. Copyright © 2011 John Wiley & Sons, Ltd.     

Direct estimation of generalized stress intensity factors using a three‐scale concurrent multigrid X‐FEM

A concurrent multigrid method is devised for the direct estimation of stress intensity factors and higher‐order coefficients of the elastic crack tip asymptotic field. The proposed method bridges three characteristic length scales that can be present in fracture mechanics: the structure, the crack and the singularity at the crack tip. For each of them, a relevant model is proposed. First, a truncated analytical reduced‐order model based on Williams' expansion is used to describe the singularity at the tip. Then, it is coupled with a standard extended finite element (FE) method model which is known to be suitable for the scale of the crack. A multigrid solver finally bridges the scale of the crack to that of the structure for which a standard FE model is often accurate enough. Dedicated coupling algorithms are presented and the effects of their parameters are discussed. The efficiency and accuracy of this new approach are exemplified using three benchmarks. Copyright © 2010 John Wiley & Sons, Ltd.     

An X‐FEM and level set computational approach for image‐based modelling: Application to homogenization

The advances in material characterization by means of imaging techniques require powerful computational methods for numerical analysis. The present contribution focuses on highlighting the advantages of coupling the extended finite elements method and the level sets method, applied to solve microstructures with complex geometries. The process of obtaining the level set data starting from a digital image of a material structure and its input into an extended finite element framework is presented. The coupled method is validated using reference examples and applied to obtain homogenized properties for heterogeneous structures. Although the computational applications presented here are mainly two‐dimensional, the method is equally applicable for three‐dimensional problems. Copyright © 2010 John Wiley & Sons, Ltd.     

X‐FEM in isogeometric analysis for linear fracture mechanics

The extended finite element method (X‐FEM) has proven to be an accurate, robust method for solving problems in fracture mechanics. X‐FEM has typically been used with elements using linear basis functions, although some work has been performed using quadratics. In the current work, the X‐FEM formulation is incorporated into isogeometric analysis to obtain solutions with higher order convergence rates for problems in linear fracture mechanics. In comparison with X‐FEM with conventional finite elements of equal degree, the NURBS‐based isogeometric analysis gives equal asymptotic convergence rates and equal accuracy with fewer degrees of freedom (DOF). Results for linear through quartic NURBS basis functions are presented for a multiplicity of one or a multiplicity equal the degree. Copyright © 2011 John Wiley & Sons, Ltd.     

Crack growth analysis using mesh superposition technique and X‐FEM

In this paper, a new approach to crack propagation analysis for large scale or complicated geometry structures is presented. That is the connection of mesh superposition technique and the extended finite element method. The former is a technique that increases the accuracy of analysis locally by superimposing additional mesh of higher resolution onto the mesh that represents the rough deformation of the structure. In this technique, the boundaries and nodes of both meshes do not have to coincide with each other. In the latter technique, the discontinuity across the crack segment and singularity around the crack tips are represented in the approximation by enriching the nodal degrees using partition of unity condition. This technique does not require meshes to conform to the crack geometry and this enables crack propagation procedure with no re‐meshing process by connecting both techniques, it becomes possible to analyse crack growth models of large scale and complicated shape with highly flexible modelling and no remeshing process. By using this approach, some numerical examples are analysed and appropriate results are obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

Improved implementation and robustness study of the X‐FEM for stress analysis around cracks

Numerical crack propagation schemes were augmented in an elegant manner by the X‐FEM method. The use of special tip enrichment functions, as well as a discontinuous function along the sides of the crack allows one to do a complete crack analysis virtually without modifying the underlying mesh, which is of industrial interest, especially when a numerical model for crack propagation is desired. This paper improves the implementation of the X‐FEM method for stress analysis around cracks in three ways. First, the enrichment strategy is revisited. The conventional approach uses a ‘topological’ enrichment (only the elements touching the front are enriched). We suggest a ‘geometrical’ enrichment in which a given domain size is enriched. The improvements obtained with this enrichment are discussed. Second, the conditioning of the X‐FEM both for topological and geometrical enrichments is studied. A preconditioner is introduced so that ‘off the shelf’ iterative solver packages can be used and perform as well on X‐FEM matrices as on standard FEM matrices. The preconditioner uses a local (nodal) Cholesky based decomposition. Third, the numerical integration scheme to build the X‐FEM stiffness matrix is dramatically improved for tip enrichment functions by the use of an ad hoc integration scheme. A 2D benchmark problem is designed to show the improvements and the robustness. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling arbitrary crack propagation in coupled shell/solid structures with X‐FEM

In this paper, the extended finite element method (X‐FEM) formulation for the modeling of arbitrary crack propagation in coupled shell/solid structures is developed based on the large deformation continuum‐based (CB) shell theory. The main features of the new method are as follows: (1) different kinematic equations are derived for different fibers in CB shell elements, including the fibers enriched by shifted jump function or crack tip functions and the fibers cut into two segments by the crack surface or connecting with solid elements. So the crack tip can locate inside the element, and the crack surface is not necessarily perpendicular to the middle surface. (2) The enhanced CB shell element is developed to realize the seamless transition of crack propagation between shell and solid structures. (3) A revised interaction integral is used to calculate the stress intensity factor (SIF) for shells, which avoids that the auxiliary fields for cracks in Mindlin–Reissner plates cannot satisfy exactly the equilibrium equations. Several numerical examples, including the calculation of SIF for the cracked plate under uniform bending and crack propagation between solid and shell structures are presented to demonstrate the performance of the developed method. Copyright © 2015 John Wiley & Sons, Ltd.     

A stable X‐FEM in cohesive transition from closed to open crack

Ignoring crack tip effects, the stability of the X‐FEM discretizations is trivial for open cracks but remains a challenge if we constrain the crack to be closed (i.e., the bi‐material problem). Here, we develop a formulation for general cohesive interactions between crack faces within the X‐FEM framework. The stability of the new formulation is demonstrated for any cohesive crack stiffness (including the closed crack) and illustrated for a nonlinear cohesive softening law. A benchmark of the new model is carried out with simpler approaches for a closed crack (i.e., Lagrange multipliers) and for a cohesive crack (i.e., penalty approach). Due to the analogies between stable cohesive X‐FEM and Nitsche's methods, the new method simplifies the implementation and is attractive in dynamic explicit codes. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐order generalized FEM for through‐the‐thickness branched cracks

This paper presents high‐order implementations of a generalized finite element method for through‐the‐thickness three‐dimensional branched cracks. This approach can accurately represent discontinuities such as triple joints in polycrystalline materials and branched cracks, independently of the background finite element mesh. Representative problems are investigated to illustrate the accuracy of the method in combination with various discretizations and refinement strategies. The combination of local refinement at crack fronts and high‐order continuous and discontinuous enrichments proves to be an excellent combination which can deliver convergence rates close to that of problems with smooth solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

A 3D beam‐to‐beam contact finite element for coupled electric–mechanical fields

In this paper the formulation of an electric–mechanical beam‐to‐beam contact element is presented. Beams with circular cross‐sections are assumed to get in contact in a point‐wise manner and with clean metallic surfaces. The voltage distribution is influenced by the contact mechanics, since the current flow is constricted to small contacting spots. Therefore, the solution is governed by the contacting areas and hence by the contact forces. As a consequence the problem is semi‐coupled with the mechanical field influencing the electric one. The electric–mechanical contact constraints are enforced with the penalty method within the finite element technique. The virtual work equations for the mechanical and electric fields are written and consistently linearized to achieve a good level of computational efficiency with the finite element method. The set of equations is solved with a monolithic approach. Copyright © 2005 John Wiley & Sons, Ltd.     

Three‐dimensional numerical solution for wear prediction using a mortar contact algorithm

A finite element formulation to compute the wear between three‐dimensional flexible bodies that are in contact with each other is presented. The contact pressure and the bodies displacements are calculated using an augmented Lagrangian approach in combination with a mortar method, which defines the contact kinematics. The objective of this study is to characterize the wear rate coefficients for bimetallic pairs and to numerically predict the wear depths in new component designs. The proposed method is first validated with the classical pin‐on‐disc problem. Then, experimental results of wear for the metallic pairs used in internal combustion engine valves and inserts are presented and are taken as a reference solution. An example is provided that shows agreement of the numerical and experimental solution. Finally, the proposed algorithm is used to predict the wear in an application example: the wear in an internal combustion engine valve. Copyright © 2013 John Wiley & Sons, Ltd.     

Arbitrary discontinuities in space–time finite elements by level sets and X‐FEM

An enriched finite element method with arbitrary discontinuities in space–time is presented. The discontinuities are treated by the extended finite element method (X‐FEM), which uses a local partition of unity enrichment to introduce discontinuities along a moving hyper‐surface which is described by level sets. A space–time weak form for conservation laws is developed where the Rankine–Hugoniot jump conditions are natural conditions of the weak form. The method is illustrated in the solution of first order hyperbolic equations and applied to linear first order wave and non‐linear Burgers' equations. By capturing the discontinuity in time as well as space, results are improved over capturing the discontinuity in space alone and the method is remarkably accurate. Implications to standard semi‐discretization X‐FEM formulations are also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A node‐to‐node scheme with the aid of variable‐node elements for elasto‐plastic contact analysis

A strategy for a two‐dimensional contact analysis involving finite strain plasticity is developed with the aid of variable‐node elements. The variable‐node elements, in which nodes are added freely where they are needed, make it possible to transform the non‐matching meshes into matching meshes directly. They thereby facilitate an efficient analysis, maintaining node‐to‐node contact during the contact deformation. The contact patch test, wherein the contact patch is constructed out of variable‐node elements, is thus passed, and iterations for equilibrium solutions reach convergence faster in this scheme than in the conventional approach based on the node‐to‐surface contact. The effectiveness and accuracy of the proposed scheme are demonstrated through several numerical examples of elasto‐plastic contact analyses. Copyright © 2015 John Wiley & Sons, Ltd.