A conformal decomposition finite element method for arbitrary discontinuities on moving interfaces

Interface capturing methods using enriched finite element formulations are well suited for solving multimaterial transport problems that contain weak or strong discontinuities. The conformal decomposition FEM decomposes multimaterial elements of a non‐conforming background mesh into sub‐elements that conform to material interfaces captured using a level set method. As the interface evolves, interfacial nodes move, and background nodes may change material. The present work describes approaches for handling moving interfaces in the context of the conformal decomposition FEM for both weakly and strongly discontinuous fields. Dynamic discretization methods using extrapolation and moving mesh approaches are considered and developed with first‐order and second‐order time integration methods. The moving mesh approach is demonstrated to be a stable method that preserves both weak and strong discontinuities on a variety of one‐dimensional and two‐dimensional test problems, while achieving the expected second‐order error convergence rate in space and time. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite cover method with mortar elements for elastoplasticity problems

Finite cover method (FCM) is extended to elastoplasticity problems. The FCM, which was originally developed under the name of manifold method, has recently been recognized as one of the generalized versions of finite element methods (FEM). Since the mesh for the FCM can be regular and squared regardless of the geometry of structures to be analyzed, structural analysts are released from a burdensome task of generating meshes conforming to physical boundaries. Numerical experiments are carried out to assess the performance of the FCM with such discretization in elastoplasticity problems. Particularly to achieve this accurately, the so-called mortar elements are introduced to impose displacement boundary conditions on the essential boundaries, and displacement compatibility conditions on material interfaces of two-phase materials or on joint surfaces between mutually incompatible meshes. The validity of the mortar approximation is also demonstrated in the elastic-plastic FCM.     

A combined space–time extended finite element method

The Newmark method for the numerical integration of second order equations has been extensively used and studied along the past fifty years for structural dynamics and various fields of mechanical engineering. Easy implementation and nice properties of this method and its derivatives for linear problems are appreciated but the main drawback is the treatment of discontinuities. Zienkiewicz proposed an approach using finite element concept in time, which allows a new look at the Newmark method. The idea of this paper is to propose, thanks to this approach, the use of a time partition of the unity method denoted Time Extended Finite Element Method (TX‐FEM) for improved numerical simulations of time discontinuities. An enriched basis of shape functions in time is used to capture with a good accuracy the non‐polynomial part of the solution. This formulation allows a suitable form of the time‐stepping formulae to study stability and energy conservation. The case of an enrichment with the Heaviside function is developed and can be seen as an alternative approach to time discontinuous Galerkin method (T‐DGM), stability and accuracy properties of which can be derived from those of the TX‐FEM. Then Space and Time X‐FEM (STX‐FEM) are combined to obtain a unified space–time discretization. This combined STX‐FEM appears to be a suitable technique for space–time discontinuous problems like dynamic crack propagation or other applications involving moving discontinuities. Copyright © 2005 John Wiley & Sons, Ltd.     

New development in extended finite element modeling of large elasto‐plastic deformations

This paper presents new achievements in the extended finite element modeling of large elasto‐plastic deformation in solid problems. The computational technique is presented based on the extended finite element method (X‐FEM) coupled with the Lagrangian formulation in order to model arbitrary interfaces in large deformations. In X‐FEM, the material interfaces are represented independently of element boundaries, and the process is accomplished by partitioning the domain with some triangular sub‐elements whose Gauss points are used for integration of the domain of elements. The large elasto‐plastic deformation formulation is employed within the X‐FEM framework to simulate the non‐linear behavior of materials. The interface between two bodies is modeled by using the X‐FEM technique and applying the Heaviside‐ and level‐set‐based enrichment functions. Finally, several numerical examples are analyzed, including arbitrary material interfaces, to demonstrate the efficiency of the X‐FEM technique in large plasticity deformations. Copyright © 2008 John Wiley & Sons, Ltd.     

Multi‐time‐step explicit–implicit method for non‐linear structural dynamics

We present a method with domain decomposition to solve time‐dependent non‐linear problems. This method enables arbitrary numeric schemes of the Newmark family to be coupled with different time steps in each subdomain: this coupling is achieved by prescribing continuity of velocities at the interface. We are more specifically interested in the coupling of implicit/explicit numeric schemes taking into account material and geometric non‐linearities. The interfaces are modelled using a dual Schur formulation where the Lagrange multipliers represent the interfacial forces. Unlike the continuous formulation, the discretized formulation of the dynamic problem is unable to verify simultaneously the continuity of displacements, velocities and accelerations at the interfaces. We show that, within the framework of the Newmark family of numeric schemes, continuity of velocities at the interfaces enables the definition of an algorithm which is stable for all cases envisaged. To prove this stability, we use an energy method, i.e. a global method over the whole time interval, in order to verify the algorithms properties. Then, we propose to extend this to non‐linear situations in the following cases: implicit linear/explicit non‐linear, explicit non‐linear/explicit non‐linear and implicit non‐linear/explicit non‐linear. Finally, we present some examples showing the feasibility of the method. Copyright © 2001 John Wiley & Sons, Ltd.     

A face‐based smoothed finite element method (FS‐FEM) for 3D linear and geometrically non‐linear solid mechanics problems using 4‐node tetrahedral elements

This paper presents a novel face‐based smoothed finite element method (FS‐FEM) to improve the accuracy of the finite element method (FEM) for three‐dimensional (3D) problems. The FS‐FEM uses 4‐node tetrahedral elements that can be generated automatically for complicated domains. In the FS‐FEM, the system stiffness matrix is computed using strains smoothed over the smoothing domains associated with the faces of the tetrahedral elements. The results demonstrated that the FS‐FEM is significantly more accurate than the FEM using tetrahedral elements for both linear and geometrically non‐linear solid mechanics problems. In addition, a novel domain‐based selective scheme is proposed leading to a combined FS/NS‐FEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The implementation of the FS‐FEM is straightforward and no penalty parameters or additional degrees of freedom are used. The computational efficiency of the FS‐FEM is found better than that of the FEM. Copyright © 2008 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 2‐D meshless model for jointed rock structures

According to the characteristic structural features of jointed rock structures, a meshless model is proposed for the mechanics analysis of jointed rock structures based on the moving least‐squares interpolants. In this model, a jointed rock structure is regarded as a system of relatively intact rock blocks connected by joints or planes of discontinuity; these rock blocks are modelled by general shaped anisotropic blocks while these joints and planes of discontinuity are modelled by interfaces. The displacement field of each block is constructed by the moving least‐squares interpolants with an array of points distributed in the block. To deal with the discontinuities of rock structures, the displacement fields are constructed to be discontinuous between blocks. The displacement fields and their gradients are continuous in each block, hence no post processing is required for the output of strains and stresses. The finite element mesh is totally unnecessary, so the time‐consuming mesh generation is avoided. The rate of convergence can exceed that of finite elements significantly, and a high resolution of localized steep gradients can be achieved. Furthermore, the discontinuities of rock structures are also fully taken into consideration. The present method is developed for two‐dimensional linear elastic analysis of jointed rock structures, and can be extended to three‐dimensional and non‐linear analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

Design sensitivity analysis in large deformation elasto‐plastic and elasto‐viscoplastic problems

Implicit time integration algorithm derived by Simo for his large‐deformation elasto‐plastic constitutive model is generalized, for the case of isotropy and associative flow rule, towards viscoplastic material behaviour and consistently differentiated with respect to its input parameters. Combining it with the general formulation of design sensitivity analysis (DSA) for non‐linear finite element transient equilibrium problem, we come at a numerically efficient, closed‐form finite element formulation of DSA for large deformation elasto‐plastic and elasto‐viscoplastic problems, with various types of design variables (material constants, shape parameters). The paper handles several specific issues, like the use of a non‐algorithmic coefficient matrix or sensitivity discontinuities at points of instantaneous structural stiffness change. Computational examples demonstrate abilities of the formulation and quality of results. Copyright © 2005 John Wiley & Sons, Ltd.     

Higher‐order extended FEM for weak discontinuities – level set representation,quadrature and application to magneto‐mechanical problems

In this article, we present the application of bilinear and biquadratic extended FEM (XFEM) formulations to model weak discontinuities in magnetic and coupled magneto‐mechanical boundary value problems. For properly resolving the location of curved interfaces and the discontinuous physical behaviour, the major part of the contribution is devoted to review and develop methods for level set representation of curved interfaces and numerical integration of the weak form in higher‐order XFEM formulations. In order to reduce the complexity of the representation of curved interfaces, an element local approach that allows for an automated computation of the level set values and also improves the compatibility between the level set representation and the integration subdomains is proposed. Integration rules for polygons and strain smoothing are applied in conjunction with biquadratic elements and compared with curved integration subdomains. Eventually, a coupled magneto‐mechanical demonstration problem is modelled and solved by XFEM. For demonstration purposes, a magneto‐mechanical coupling due to magnetic stresses is considered. Errors and convergence rates are analysed for the different level set representations and numerical integration procedures as well as their dependence on the ratio of material parameters at an interface. Copyright © 2012 John Wiley & Sons, Ltd.     

Remarks on the interpretation of current non‐linear finite element analyses as differential–algebraic equations

For the numerical solution of materially non‐linear problems like in computational plasticity or viscoplasticity the finite element discretization in space is usually coupled with point‐wise defined evolution equations characterizing the material behaviour. The interpretation of such systems as differential–algebraic equations (DAE) allows modern‐day integration algorithms from Numerical Mathematics to be efficiently applied. Especially, the application of diagonally implicit Runge–Kutta methods (DIRK) together with a Multilevel‐Newton method preserves the algorithmic structure of current finite element implementations which are based on the principle of virtual displacements and on backward Euler schemes for the local time integration. Moreover, the notion of the consistent tangent operator becomes more obvious in this context. The quadratical order of convergence of the Multilevel‐Newton algorithm is usually validated by numerical studies. However, an analytical proof of this second order convergence has already been given by authors in the field of non‐linear electrical networks. We show that this proof can be applied in the current context based on the DAE interpretation mentioned above. We finally compare the proposed procedure to several well‐known stress algorithms and show that the inclusion of a step‐size control based on local error estimations merely requires a small extra time‐investment. Copyright © 2001 John Wiley & Sons, Ltd.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.     

A boundary cohesive grain element formulation for modelling intergranular microfracture in polycrystalline brittle materials

In this paper, a cohesive grain boundary integral formulation is proposed, for simulating intergranular microfracture evolution in polycrystalline brittle materials. Artificially generated polycrystalline microstructures are discretized using the proposed anisotropic boundary element method, considering the random location, morphology and material orientation of each grain. Each grain is assumed as a single crystal with general elastic orthotropic mechanical behaviour. Crack initiation and propagation along the grain boundaries interfaces are modelled using a linear cohesive law, considering mixed mode failure conditions. Furthermore, a non‐linear frictional contact analysis is performed over cracked grain interfaces to encounter cases where crack surfaces come into contact, slide or separate. The effect of randomly located pre‐existing flaws on the overall behaviour and microcracking evolution of a polycrystalline material is also investigated for different Weibull moduli. The stochastic effects of each grain morphology‐orientation, internal friction and randomly distributed pre‐existing flaws, under different loading conditions, are studied probabilistically by simulating various randomly generated microstructures. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical finite element formulation of the Schapery non‐linear viscoelastic material model

This study presents a numerical integration method for the non‐linear viscoelastic behaviour of isotropic materials and structures. The Schapery's three‐dimensional (3D) non‐linear viscoelastic material model is integrated within a displacement‐based finite element (FE) environment. The deviatoric and volumetric responses are decoupled and the strain vector is decomposed into instantaneous and hereditary parts. The hereditary strains are updated at the end of each time increment using a recursive formulation. The constitutive equations are expressed in an incremental form for each time step, assuming a constant incremental strain rate. A new iterative procedure with predictor–corrector type steps is combined with the recursive integration method. A general polynomial form for the parameters of the non‐linear Schapery model is proposed. The consistent algorithmic tangent stiffness matrix is realized and used to enhance convergence and help achieve a correct convergent state. Verifications of the proposed numerical formulation are performed and compared with a previous work using experimental data for a glassy amorphous polymer PMMA. Copyright © 2003 John Wiley & Sons, Ltd.     

A general discontinuous Galerkin method for finite hyperelasticity. Formulation and numerical applications

A discontinuous Galerkin formulation of the boundary value problem of finite‐deformation elasticity is presented. The primary purpose is to establish a discontinuous Galerkin framework for large deformations of solids in the context of statics and simple material behaviour with a view toward further developments involving behaviour or models where the DG concept can show its superiority compared to the continuous formulation. The method is based on a general Hu–Washizu–de Veubeke functional allowing for displacement and stress discontinuities in the domain interior. It is shown that this approach naturally leads to the formulation of average stress fluxes at interelement boundaries in a finite element implementation. The consistency and linearized stability of the method in the non‐linear range as well as its convergence rate are proven. An implementation in three dimensions is developed, showing that the proposed method can be integrated into conventional finite element codes in a straightforward manner. In order to demonstrate the versatility, accuracy and robustness of the method examples of application and convergence studies in three dimensions are provided. Copyright © 2006 John Wiley & Sons, Ltd.     

On non‐linear behaviour of spherical shallow shells bonded with piezoelectric actuators by the differential quadrature element method (DQEM)

The static behaviour of spherical shallow shells bonded with piezoelectric actuators and subjected to electrical loading are studied in this paper by using the differential quadrature element method (DQEM). Geometrical non‐linear effects are considered. Detailed formulations for the DQ circular spherical shallow shell element and the DQ annular spherical shallow shell element are given for the first time. Numerical studies are performed to evaluate the effects of actuator size, thickness and boundary conditions. Very accurate results are obtained by the DQEM. Based on the results reported in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions for smart materials and structures exhibiting geometric non‐linear behaviours. Thickness effects cannot be neglected when the actuator thickness is comparable to that of the base material. Snap‐through may occur when the applied voltage reaches a critical value even without mechanical loading for certain geometric configurations. Copyright © 2001 John Wiley & Sons, Ltd.     

Experimental characterisation of a CuAg alloy for thermo‐mechanical applications. Part 1: Identifying parameters of non‐linear plasticity models

Despite the wide use of copper alloys in thermo‐mechanical applications, there is little data on their cyclic plasticity behaviour, particularly for CuAg alloys. This prevents the behaviour of the materials from being correctly described in numerical simulations for design purposes. In this work CuAg0.1 alloy used for thermo‐mechanical applications was tested by strain‐controlled cyclic loading at 3 different temperatures (room temperature, 250°C, 300°C). In each test, stress‐strain cycles were recorded until the alloy had completely stabilised. These cycles were then used to identify material parameters of non‐linear kinematic and isotropic models. The focus was on plasticity models (Armstrong‐Frederick, Chaboche, Voce) that are usually implemented in commercial finite element codes. Simulated cyclic responses with the identified material models were compared with experiments and showed a good agreement. The identified material parameters for the CuAg alloy under investigation can be used directly in finite element models for cyclic plasticity simulations, thus enabling a durability analysis of components under thermo‐mechanical loads to be performed, particularly in the field of steel‐making plants.     

Experimental evaluation of the elastic limit of carbon‐fibre reinforced epoxy composites under a large range of strain rate and temperature conditions

The mechanical behaviour of organic matrix composite materials such as T700GC/M21 carbon fibre reinforced polymer (CFRP) is generally considered by the industry as being orthotropic elastic for the sizing of aeronautical structures under normal isothermal “static” flight loads. During the aircraft lifetime, it may be exposed to severe loading conditions at various temperatures. However, the mechanical behaviour of CFRP is known to exhibit a linear behaviour or a non‐linear behaviour according to the types of loads that are considered creep or extreme conditions. The observed non‐linearity can be commonly attributed to several physical phenomena such as non‐linear viscosity, plasticity, or damage. In the literature, different models can be found that are based on three components: a first elastic reversible behaviour, a second non‐linear behaviour, and a failure criterion. An important issue is to understand and characterize the transition between the elastic reversible behaviour and the non‐linear behaviour. To answer this question, the present paper describes an experimental methodology that permits to evaluate this transition thanks to raw experimental data, and its application to a range of constant but different strain rate and temperature tests performed on the T700GC/M21 CFRP material.     

A geometrically non‐linear piezoelectric solid shell element based on a mixed multi‐field variational formulation

This paper is concerned with a geometrically non‐linear solid shell element to analyse piezoelectric structures. The finite element formulation is based on a variational principle of the Hu–Washizu type and includes six independent fields: displacements, electric potential, strains, electric field, mechanical stresses and dielectric displacements. The element has eight nodes with four nodal degrees of freedoms, three displacements and the electric potential. A bilinear distribution through the thickness of the independent electric field is assumed to fulfill the electric charge conservation law in bending dominated situations exactly. The presented finite shell element is able to model arbitrary curved shell structures and incorporates a 3D‐material law. A geometrically non‐linear theory allows large deformations and includes stability problems. Linear and non‐linear numerical examples demonstrate the ability of the proposed model to analyse piezoelectric devices. Copyright © 2005 John Wiley & Sons, Ltd.     

Hybrid finite element/volume method for shallow water equations

A hybrid numerical scheme based on finite element and finite volume methods is developed to solve shallow water equations. In the recent past, we introduced a series of hybrid methods to solve incompressible low and high Reynolds number flows for single and two‐fluid flow problems. The present work extends the application of hybrid method to shallow water equations. In our hybrid shallow water flow solver, we write the governing equations in non‐conservation form and solve the non‐linear wave equation using finite element method with linear interpolation functions in space. On the other hand, the momentum equation is solved with highly accurate cell‐center finite volume method. Our hybrid numerical scheme is truly a segregated method with primitive variables stored and solved at both node and element centers. To enhance the stability of the hybrid method around discontinuities, we introduce a new shock capturing which will act only around sharp interfaces without sacrificing the accuracy elsewhere. Matrix‐free GMRES iterative solvers are used to solve both the wave and momentum equations in finite element and finite volume schemes. Several test problems are presented to demonstrate the robustness and applicability of the numerical method. Copyright © 2010 John Wiley & Sons, Ltd.