A comparative study between two smoothing strategies for the simulation of contact with large sliding

The numerical simulation of contact problems is still a delicate matter especially when large transformations are involved. In that case, relative large slidings can occur between contact surfaces and the discretization error induced by usual finite elements may not be satisfactory. In particular, usual elements lead to a facetization of the contact surface, meaning an unavoidable discontinuity of the normal vector to this surface. Uncertainty over the precision of the results, irregularity of the displacement of the contact nodes and even numerical oscillations of contact reaction force may result of such discontinuity. Among the existing methods for tackling such issue, one may consider mortar elements (Fischer and Wriggers, Comput Methods Appl Mech Eng 195:5020–5036, 2006; McDevitt and Laursen, Int J Numer Methods Eng 48:1525–1547, 2000; Puso and Laursen, Comput Methods Appl Mech Eng 93:601–629, 2004), smoothing of the contact surfaces with additional geometrical entity (B-splines or NURBS) (Belytschko et al., Int J Numer Methods Eng 55:101–125, 2002; Kikuchi, Penalty/finite element approximations of a class of unilateral contact problems. Penalty method and finite element method, ASME, New York, 1982; Legrand, Modèles de prediction de l’interaction rotor/stator dans un moteur d’avion Thèse de doctorat. PhD thesis, École Centrale de Nantes, Nantes, 2005; Muñoz, Comput Methods Appl Mech Eng 197:979–993, 2008; Wriggers and Krstulovic-Opara, J Appl Math Mech (ZAMM) 80:77–80, 2000) and, the use of isogeometric analysis (Temizer et al., Comput Methods Appl Mech Eng 200:1100–1112, 2011; Hughes et al., Comput Methods Appl Mech Eng 194:4135–4195, 2005; de Lorenzis et al., Int J Numer Meth Eng, in press, 2011). In the present paper, we focus on these last two methods which are combined with a finite element code using the bi-potential method for contact management (Feng et al., Comput Mech 36:375–383, 2005). A comparative study focusing on the pros and cons of each method regarding geometrical precision and numerical stability for contact solution is proposed. The scope of this study is limited to 2D contact problems for which we consider several types of finite elements. Test cases are given in order to illustrate this comparative study.     

Finite element shakedown analysis of two‐dimensional structures

In the present paper, the formulation proposed by Casciaro and Garcea (Comput. Meth. Appl. Mech. Eng., 2002; 191 :5761–5792) and applied to the shakedown analysis of plane frames, is extended to the analysis of two‐dimensional flat structures in both the cases of plane‐stress and plane‐strain. The discrete formulation is obtained using a mixed finite element in which both stress and displacement fields are interpolated. The material is assumed to be elasto‐plastic and a linearization of the elastic domain is performed. The result is a versatile iterative scheme well suited to implementation in general purpose FEM codes. An extensive series of numerical tests is presented showing the reliability of the proposed formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

The role of continuity in residual-based variational multiscale modeling of turbulence

This paper examines the role of continuity of the basis in the computation of turbulent flows. We compare standard finite elements and non-uniform rational B-splines (NURBS) discretizations that are employed in Isogeometric Analysis (Hughes et al. in Comput Methods Appl Mech Eng, 194:4135–4195, 2005). We make use of quadratic discretizations that are C ⁰-continuous across element boundaries in standard finite elements, and C ¹-continuous in the case of NURBS. The variational multiscale residual-based method (Bazilevs in Isogeometric analysis of turbulence and fluid-structure interaction, PhD thesis, ICES, UT Austin, 2006; Bazilevs et al. in Comput Methods Appl Mech Eng, submitted, 2007; Calo in Residual-based multiscale turbulence modeling: finite volume simulation of bypass transition. PhD thesis, Department of Civil and Environmental Engineering, Stanford University, 2004; Hughes et al. in proceedings of the XXI international congress of theoretical and applied mechanics (IUTAM), Kluwer, 2004; Scovazzi in Multiscale methods in science and engineering, PhD thesis, Department of Mechanical Engineering, Stanford Universty, 2004) is employed as a turbulence modeling technique. We find that C ¹-continuous discretizations outperform their C ⁰-continuous counterparts on a per-degree-of-freedom basis. We also find that the effect of continuity is greater for higher Reynolds number flows.     

Theory and finite element computations of a unified cyclic phase transformation model for monocrystalline materials at small strains

After a survey the refined numerical treatment and verification is presented for a rate-independent macroscopic unified PT material model (including mass conservation with respect to phase fractions and covexified free energy) by Govindjee and Miehe (Comput Methods Appl Mech Eng 191:215–238, 2001) for describing SME and SE effects within a linear kinematic setting. Special attention is given to temperature dependent PTs. The material model was implemented into ABAQUS via the UMAT material interface in 2004. Validation of this PT model is carried out with experimental data supplied by Xiangyang et al. (J Mech Phys Solids 48:2163–2182, 2000), using 3D finite element computations. Experimentally gained material data from different sources are used and numerical results of energy barriers for PTs are given. Another feature is the simulation of suppressed shape memory effects by quasiplastic temperature induced PT. Furthermore, a plane strain problem is treated with comparisons of butterfly shaped expansions of martensitic PT and plastic deformation, correspondingly.     

Boundary element solution of the two-dimensional sine-Gordon equation using continuous linear elements

This article studies the boundary element solution of two-dimensional sine-Gordon (SG) equation using continuous linear elements approximation. Non-linear and in-homogenous terms are converted to the boundary by the dual reciprocity method and a predictor–corrector scheme is employed to eliminate the non-linearity. The procedure developed in this paper, is applied to various problems involving line and ring solitons where considered in references [Argyris J, Haase M, Heinrich JC. Finite element approximation to two-dimensional sine-Gordon solitons. Comput Methods Appl Mech Eng 1991;86:1–26; Bratsos AG. An explicit numerical scheme for the sine-Gordon equation in 2+1 dimensions. Appl Numer Anal Comput Math 2005;2(2):189–211, Bratsos AG. A modified predictor–corrector scheme for the two-dimensional sine-Gordon equation. Numer Algorithms 2006;43:295–308; Bratsos AG. The solution of the two-dimensional sine-Gordon equation using the method of lines. J Comput Appl Math 2007;206:251–77; Bratsos AG. A third order numerical scheme for the two-dimensional sine-Gordon equation. Math Comput Simul 2007;76:271–8; Christiansen PL, Lomdahl PS. Numerical solutions of 2+1 dimensional sine-Gordon solitons. Physica D: Nonlinear Phenom 1981;2(3):482–94; Djidjeli K, Price WG, Twizell EH. Numerical solutions of a damped sine-Gordon equation in two space variables. J Eng Math 1995;29:347–69; Dehghan M, Mirzaei D. The dual reciprocity boundary element method (DRBEM) for two-dimensional sine-Gordon equation. Comput Methods Appl Mech Eng 2008;197:476–86]. Using continuous linear elements approximation produces more accurate results than constant ones. By using this approach all cases associated to SG equation, which exist in literature, are investigated.     

A Multiscale/stabilized Formulation of the Incompressible Navier–Stokes Equations for Moving Boundary Flows and Fluid–structure Interaction

This paper presents a multiscale/stabilized finite element formulation for the incompressible Navier–Stokes equations written in an Arbitrary Lagrangian–Eulerian (ALE) frame to model flow problems that involve moving and deforming meshes. The new formulation is derived based on the variational multiscale method proposed by Hughes (Comput Methods Appl Mech Eng 127:387–401, 1995) and employed in Masud and Khurram in (Comput Methods Appl Mech Eng 193:1997–2018, 2006); Masud and Khurram in (Comput Methods Appl Mech Eng 195:1750–1777, 2006) to study advection dominated transport phenomena. A significant feature of the formulation is that the structure of the stabilization terms and the definition of the stabilization tensor appear naturally via the solution of the sub-grid scale problem. A mesh moving technique is integrated in this formulation to accommodate the motion and deformation of the computational grid, and to map the moving boundaries in a rational way. Some benchmark problems are shown, and simulations of an elastic beam undergoing large amplitude periodic oscillations in a viscous fluid domain are presented.     

On the approximation in the smoothed finite element method (SFEM)

This letter aims at resolving the issues raised in the recent short communication (Int. J. Numer. Meth. Engng 2008; 76 (8):1285–1295. DOI: 10.1002/nme.2460 ) and answered by (Int. J. Numer. Meth. Engng 2009; DOI: 10.1002/nme.2587 ) by proposing a systematic approximation scheme based on non‐mapped shape functions, which both allows to fully exploit the unique advantages of the smoothed finite element method (SFEM) (Comput. Mech. 2007; 39 (6):859–877. DOI: 10.1007/s00466‐006‐0075‐4 ; Commun. Numer. Meth. Engng 2009; 25 (1):19–34. DOI: 10.1002/cnm.1098 ; Int. J. Numer. Meth. Engng 2007; 71 (8):902–930; Comput. Meth. Appl. Mech. Engng 2008; 198 (2):165–177. DOI: 10.1016/j.cma.2008.05.029 ; Comput. Meth. Appl. Mech. Engng 2007; submitted; Int. J. Numer. Meth. Engng 2008; 74 (2):175–208. DOI: 10.1002/nme.2146 ; Comput. Meth. Appl. Mech. Engng 2008; 197 (13–16):1184–1203. DOI: 10.1016/j.cma.2007.10.008 ) and resolve the existence, linearity and positivity deficiencies pointed out in (Int. J. Numer. Meth. Engng 2008; 76 (8):1285–1295). We show that Wachspress interpolants (A Rational Basis for Function Approximation. Academic Press, Inc.: New York, 1975) computed in the physical coordinate system are very well suited to the SFEM, especially when elements are heavily distorted (obtuse interior angles). The proposed approximation leads to results that are almost identical to those of the SFEM initially proposed in (Comput. Mech. 2007; 39 (6):859–877. DOI: 10.1007/s00466‐006‐0075‐4 ). These results suggest that the proposed approximation scheme forms a strong and rigorous basis for the construction of SFEMs. Copyright © 2009 John Wiley & Sons, Ltd.     

Strain-injection and crack-path field techniques for 3D crack-propagation modelling in quasi-brittle materials

This paper presents a finite element approach for modelling three-dimensional crack propagation in quasi-brittle materials, based on the strain injection and the crack-path field techniques. These numerical techniques were already tested and validated by static and dynamic simulations in 2D classical benchmarks [Dias et al., in: Monograph CIMNE No-134. International Center for Numerical Methods in Engineering, Barcelona, (2012); Oliver et al. in Comput Methods Appl Mech Eng 274:289–348, (2014); Lloberas-Valls et al. in Comput Methods Appl Mech Eng 308:499–534, (2016)] and, also, for modelling tensile crack propagation in real concrete structures, like concrete gravity dams [Dias et al. in Eng Fract Mech 154:288–310, (2016)]. The main advantages of the methodology are the low computational cost and the independence of the results on the size and orientation of the finite element mesh. These advantages were highlighted in previous works by the authors and motivate the present extension to 3D cases. The proposed methodology is implemented in the finite element framework using continuum constitutive models equipped with strain softening and consists, essentially, in injecting the elements candidate to capture the cracks with some goal oriented strain modes for improving the performance of the injected elements for simulating propagating displacement discontinuities. The goal-oriented strain modes are introduced by resorting to mixed formulations and to the Continuum Strong Discontinuity Approach (CSDA), while the crack position inside the finite elements is retrieved by resorting to the crack-path field technique. Representative numerical simulations in 3D benchmarks show that the advantages of the methodology already pointed out in 2D are kept in 3D scenarios.     

Continuum phase field modeling of dynamic fracture: variational principles and staggered FE implementation

The modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations of complex crack topologies including branching. This drawback can be overcome by a diffusive crack modeling based on the introduction of a crack phase field as proposed in Miehe et?al. (Comput Methods Appl Mech Eng 19:2765?C2778, 2010a; Int J Numer Meth Eng 83:1273?C1311, 2010b), Hofacker and Miehe (Int J Numer Meth Eng, 2012). In this work, we summarize basic ingredients of a thermodynamically consistent, variational-based model of diffusive crack propagation under quasi-static and dynamic conditions. It is shown that all coupled field equations, in particular the balance of momentum and the gradient-type evolution equation for the crack phase field, follow as the Euler equations of a mixed rate-type variational principle that includes the fracture driving force as the mixed field variable. This principle makes the proposed formulation extremely compact and provides a perfect basis for the finite element implementation. We then introduce a local history field that contains a maximum energetic crack source obtained in the deformation history. It drives the evolution of the crack phase field. This allows for the construction of an extremely robust operator split scheme that updates in a typical time step the history field, the crack phase field and finally the displacement field. We demonstrate the performance of the phase field formulation of fracture by means of representative numerical examples, which show the evolution of complex crack patterns under dynamic loading.     

Study on sub-cycling algorithm for flexible multi-body system: stability analysis and numerical examples

Numerical stability is an important issue for any integral procedure. Since sub-cycling algorithm has been presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18: 259–275, 1979), various kinds of these integral procedures were developed in later 20 years and their stability were widely studied. However, on how to apply the sub-cycling to flexible multi-body dynamics (FMD) is still a lack of investigation up to now. A particular sub-cycling algorithm for the FMD based on the central difference method was introduced in detail in part I (Miao et al. in Comp Mech doi: 10.1007/s00466-007-0183-9) of this paper. Adopting an integral approximation operator method, stability of the presented algorithm is transformed to a generalized eigenvalue problem in the paper and is discussed by solving the problem later. Numerical examples are performed to verify the availability and efficiency of the algorithm further.     

Improved Discontinuity-capturing Finite Element Techniques for Reaction Effects in Turbulence Computation

Recent advances in turbulence modeling brought more and more sophisticated turbulence closures (e.g. k-ɛ, k-ɛ -v ²-f, Second Moment Closures), where the governing equations for the model parameters involve advection, diffusion and reaction terms. Numerical instabilities can be generated by the dominant advection or reaction terms. Classical stabilized formulations such as the Streamline–Upwind/Petrov–Galerkin (SUPG) formulation (Brook and Hughes, comput methods Appl Mech Eng 32:199–255, 1982; Hughes and Tezduyar, comput methods Appl Mech Eng 45: 217–284, 1984) are very well suited for preventing the numerical instabilities generated by the dominant advection terms. A different stabilization however is needed for instabilities due to the dominant reaction terms. An additional stabilization term, called the diffusion for reaction-dominated (DRD) term, was introduced by Tezduyar and Park (comput methods Appl Mech Eng 59:307–325, 1986) for that purpose and improves the SUPG performance. In recent years a new class of variational multi-scale (VMS) stabilization (Hughes, comput methods Appl Mech Eng 127:387–401, 1995) has been introduced, and this approach, in principle, can deal with advection–diffusion–reaction equations. However, it was pointed out in Hanke (comput methods Appl Mech Eng 191:2925–2947) that this class of methods also need some improvement in the presence of high reaction rates. In this work we show the benefits of using the DRD operator to enhance the core stabilization techniques such as the SUPG and VMS formulations. We also propose a new operator called the DRDJ (DRD with the local variation jump) term, targeting the reduction of numerical oscillations in the presence of both high reaction rates and sharp solution gradients. The methods are evaluated in the context of two stabilized methods: the classical SUPG formulation and a recently-developed VMS formulation called the V-SGS (Corsini et al. comput methods Appl Mech Eng 194:4797–4823, 2005). Model problems and industrial test cases are computed to show the potential of the proposed methods in simulation of turbulent flows.     

3D XFEM investigation of the plasticity effect on fatigue propagation under thermo-mechanical loading

The aim of this paper is to propose a computation strategy for fatigue propagation simulation of a crack by taking into account the plasticity. Feulvarch et al. (Comput Methods Appl Mech Eng 361: 112805, 2020) recently proposed a first XFEM formulation capable of overcoming the volumetric locking phenomenon due to plastic incompressibility in 3D. This formulation is here applied to quadratic elements for the mode I propagation of a crack in a valve structure submitted to cyclic thermo-mechanical loading. A simulation strategy is proposed where it is not necessary to compute all the cycles and thus the complete plastic history. This is of great interest because it avoids the treatment of the possible closing of the crack and uses the conventional J-integral. The application reveals the interest of taking plasticity into account for the propagation accuracy.

     

A set of pathological tests to validate new finite elements

The finite element method entails several approximations. Hence it is essential to subject all new finite elements to an adequate set of pathological tests in order to assess their performance. Many such tests have been proposed by researchers from time to time. We present an adequate set of tests, which every new finite element should pass. A thorough account of the patch test is also included in view of its significance in the validation of new elements.     

Hyperelastic modelling of arterial layers with distributed collagen fibre orientations.

Constitutive relations are fundamental to the solution of problems in continuum mechanics, and are required in the study of, for example, mechanically dominated clinical interventions involving soft biological tissues. Structural continuum constitutive models of arterial layers integrate information about the tissue morphology and therefore allow investigation of the interrelation between structure and function in response to mechanical loading. Collagen fibres are key ingredients in the structure of arteries. In the media (the middle layer of the artery wall) they are arranged in two helically distributed families with a small pitch and very little dispersion in their orientation (i.e. they are aligned quite close to the circumferential direction). By contrast, in the adventitial and intimal layers, the orientation of the collagen fibres is dispersed, as shown by polarized light microscopy of stained arterial tissue. As a result, continuum models that do not account for the dispersion are not able to capture accurately the stress-strain response of these layers. The purpose of this paper, therefore, is to develop a structural continuum framework that is able to represent the dispersion of the collagen fibre orientation. This then allows the development of a new hyperelastic free-energy function that is particularly suited for representing the anisotropic elastic properties of adventitial and intimal layers of arterial walls, and is a generalization of the fibre-reinforced structural model introduced by Holzapfel & Gasser (Holzapfel & Gasser 2001 Comput. Meth. Appl. Mech. Eng. 190, 4379-4403) and Holzapfel et al. (Holzapfel et al. 2000 J. Elast. 61, 1-48). The model incorporates an additional scalar structure parameter that characterizes the dispersed collagen orientation. An efficient finite element implementation of the model is then presented and numerical examples show that the dispersion of the orientation of collagen fibres in the adventitia of human iliac arteries has a significant effect on their mechanical response.     

ODDLS: A new unstructured mesh finite element method for the analysis of free surface flow problems

This paper introduces a new stabilized finite element method based on the finite calculus (Comput. Methods Appl. Mech. Eng. 1998; 151 :233–267) and arbitrary Lagrangian–Eulerian techniques (Comput. Methods Appl. Mech. Eng. 1998; 155 :235–249) for the solution to free surface problems. The main innovation of this method is the application of an overlapping domain decomposition concept in the statement of the problem. The aim is to increase the accuracy in the capture of the free surface as well as in the resolution of the governing equations in the interface between the two fluids. Free surface capturing is based on the solution to a level set equation. The Navier–Stokes equations are solved using an iterative monolithic predictor–corrector algorithm (Encyclopedia of Computational Mechanics. Wiley: New York, 2004), where the correction step is based on imposing the divergence‐free condition in the velocity field by means of the solution to a scalar equation for the pressure. Examples of application of the ODDLS formulation (for overlapping domain decomposition level set) to the analysis of different free surface flow problems are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Attachment mode performance of network-modeled ballistic fabric shielding

A central issue in the use of ballistic fabric shielding is the mode of attachment to the structure that it is intended to protect. In order to investigate this issue, a discrete multi-scale yarn-network model is developed for structural fabric undergoing ballistic impact, based on work found in Zohdi and Powell [Zohdi TI, Powell D. Multiscale construction and large-scale simulation of structural fabric undergoing ballistic impact. Comput Meth Appl Mech Eng 2006;195:94–109] and Zohdi [Zohdi TI. Modeling/simulation of progressive penetration of multilayered ballistic fabric shielding. Comput Mech 2002;29:61–7]. The model is comprised of a network of yarn with stochastic properties determined by smaller-scale fibrils, which are randomly misaligned. The effects of stochasticity on the overall response are explored, and the model is compared against macro-scale experiments. The key feature of the model is the fact that it does not depend on phenomenological parameters, and can be calibrated by simply measuring the properties of an individual, smallest-scale, fibril. The properties of a fibril are easily ascertained from a simple tension test. The response of the overall fabric model and ballistic experiments are in excellent agreement. The model indicates that fabric which is attached by being pinned at the corners generally absorbs more energy, relative to fabric clamped along the sides. The basis for this result is discussed at length in the body of this work. Furthermore, it is observed that a uniform-yarn model, one which ignores the stochastic nature of the yarn, over-estimates the amount of energy absorbed.     

A discontinuous Galerkin formulation of non‐linear Kirchhoff–Love shells

Discontinuous Galerkin (DG) methods provide a means of weakly enforcing the continuity of the unknown‐field derivatives and have particular appeal in problems involving high‐order derivatives. This feature has previously been successfully exploited (Comput. Methods Appl. Mech. Eng. 2008; 197 :2901–2929) to develop a formulation of linear Kirchhoff–Love shells considering only the membrane and bending responses. In this proposed one‐field method—the displacements are the only unknowns, while the displacement field is continuous, the continuity in the displacement derivative between two elements is weakly enforced by recourse to a DG formulation. It is the purpose of the present paper to extend this formulation to finite deformations and non‐linear elastic behaviors. While the initial linear formulation was relying on the direct linear computation of the effective membrane stress and effective bending couple‐stress from the displacement field at the mid‐surface of the shell, the non‐linear formulation considered implies the evaluation of the general stress tensor across the shell thickness, leading to a reformulation of the internal forces of the shell. Nevertheless, since the interface terms resulting from the discontinuous Galerkin method involve only the resultant couple‐stress at the edges of the shells, the extension to non‐linear deformations is straightforward. Copyright © 2008 John Wiley & Sons, Ltd.     

A fourth‐order compact scheme for the Helmholtz equation: Alpha‐interpolation of FEM and FDM stencils

We propose a fourth‐order compact scheme on structured meshes for the Helmholtz equation given by R(φ):=f( x )+Δφ+ξ²φ=0. The scheme consists of taking the alpha‐interpolation of the Galerkin finite element method and the classical central finite difference method. In 1D, this scheme is identical to the alpha‐interpolation method (J. Comput. Appl. Math. 1982; 8 (1):15–19) and in 2D making the choice α=0.5 we recover the generalized fourth‐order compact Padé approximation (J. Comput. Phys. 1995; 119 :252–270; Comput. Meth. Appl. Mech. Engrg 1998; 163 :343–358) (therein using the parameter γ=2). We follow (SIAM Rev. 2000; 42 (3):451–484; Comput. Meth. Appl. Mech. Engrg 1995; 128 :325–359) for the analysis of this scheme and its performance on square meshes is compared with that of the quasi‐stabilized FEM (Comput. Meth. Appl. Mech. Engrg 1995; 128 :325–359). In particular, we show that the relative phase error of the numerical solution and the local truncation error of this scheme for plane wave solutions diminish at the rate O((ξ?)⁴), where ξ, ? represent the wavenumber and the mesh size, respectively. An expression for the parameter α is given that minimizes the maximum relative phase error in a sense that will be explained in Section 4.5. Convergence studies of the error in the L² norm, the H¹ semi‐norm and the l^∞ Euclidean norm are done and the pollution effect is found to be small. Copyright © 2010 John Wiley & Sons, Ltd.