Velocity‐based formulations for standard and quasi‐incompressible hypoelastic‐plastic solids

We present three velocity‐based updated Lagrangian formulations for standard and quasi‐incompressible hypoelastic‐plastic solids. Three low‐order finite elements are derived and tested for non‐linear solid mechanics problems. The so‐called V‐element is based on a standard velocity approach, while a mixed velocity–pressure formulation is used for the VP and the VPS elements. The two‐field problem is solved via a two‐step Gauss–Seidel partitioned iterative scheme. First, the momentum equations are solved in terms of velocity increments, as for the V‐element. Then, the constitutive relation for the pressure is solved using the updated velocities obtained at the previous step. For the VPS‐element, the formulation is stabilized using the finite calculus method in order to solve problems involving quasi‐incompressible materials. All the solid elements are validated by solving two‐dimensional and three‐dimensional benchmark problems in statics as in dynamics. Copyright © 2016 John Wiley & Sons, Ltd.     

P1‐nonconforming quadrilateral finite element for topology optimization

This investigation focuses on an alternative approach to topology optimization problems involving incompressible materials using the P1‐nonconforming finite element. Instead of using the mixed displacement‐pressure formulation, a pure displacement‐based approach can be employed for finite element formulation owing to the Poisson locking‐free property of the P1‐nonconforming element. Moreover, because the P1‐nonconforming element has linear shape functions that are defined at element vertices, it has considerably fewer degrees of freedom than other quadrilateral nonconforming elements and its implementation is as simple as that of the conforming bilinear element. Various problems dealing with incompressible materials and pressure‐loaded structures found in published works are solved to verify the applicability of the proposed method. The application of the method is extended to the optimal design of fluid channels in the Stokes flow. This is done by expressing pressure in terms of volumetric strain rates and developing a velocity‐field‐only finite element formulation. The optimization results obtained from all the problems considered in this study are in close agreement with those found in the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

A new enhanced assumed strain quadrilateral membrane element with drilling degree of freedom and modified shape functions

A new four‐node quadrilateral membrane finite element with drilling rotational degree of freedom based on the enhanced assumed strain formulation is presented. A simple formulation is achieved by five incompatible modes that are added to the Allman‐type interpolation. Furthermore, modified shape functions are used to improve the behaviour of distorted elements. Numerical results show that the proposed new element exhibits good numerical accuracy and improved performance, and in many cases, superior to existing elements. In particular, Poisson's locking in nearly incompressible elasticity fades and the element performs well when it becomes considerably distorted even when it takes almost triangular shape. Copyright © 2016 John Wiley & Sons, Ltd.     

Patch‐averaged assumed strain finite elements for stress analysis

A finite element model for linear‐elastic small deformation problems is presented. The formulation is based on a weighted residual that requires a priori the satisfaction of the kinematic equation. In this approach, an averaged strain‐displacement matrix is constructed for each node of the mesh by defining an appropriate patch of elements, yielding a smooth representation of strain and stress fields. Connections with traditional and similar procedure are explored. Linear quadrilateral four‐node and linear hexahedral eight‐node elements are derived. Various numerical tests show the accuracy and convergence properties of the proposed elements in comparison with extant finite elements and analytic solutions. Specific examples are also included to illustrate the ability to resist numerical locking in the incompressible limit and insensitive response in the presence of shape distortion. Furthermore, the numerical inf‐sup test is applied to a selection of problems to show the stability of the present formulation. Copyright © 2012 John Wiley & Sons, Ltd.     

On the enhancement of low‐order mixed finite element methods for the large deformation analysis of diffusion in solids

A stabilized scheme is developed for mixed finite element methods for strongly coupled diffusion problems in solids capable of large deformations. Enhanced assumed strain techniques are employed to cure spurious oscillation patterns of low‐order displacement/pressure mixed formulations in the incompressible limit for quadrilateral elements and brick elements. A study is presented that shows how hourglass instabilities resulting from geometrically nonlinear enhanced assumed strain methods have to be distinguished from pressure oscillation patterns due to the violation of the inf‐sup condition. Moreover, an element formulation is proposed that provides stable results with respect to both types of instabilities. Comparisons are drawn between material models for incompressible solids of Mooney–Rivlin type and models for standard diffusion in solids with incompressible matrices such as polymeric gels. Representative numerical examples underline the ability of the proposed element formulation to cure instabilities of low‐order mixed formulations. Copyright © 2015 John Wiley & Sons, Ltd.     

On using linear elements in incompressible plane strain problems: a simple edge based approach for triangles

In this paper a simple iterative method is presented for finite element solution of incompressible plane strain problems using linear elements. Instead of using a mixed formulation approach, we use an equivalent displacement/velocity approach in an iterative manner. Control volumes are taken for regions which are to exhibit incompressible behaviour. For triangular elements the control volume is chosen as the area built on the parts of each pair of elements at the sides of an edge. In this case, elements are let to exchange volume. It is shown that the proposed edge based approach removes the deficiency of the linear triangular elements i.e. locking effect. Similar edge based approach is applied to the linear quadrilateral elements. However, if the control volume is chosen as the element volume the formulation gives similar results as the discontinuous mixed formulation using one pressure point without exhibiting instability behaviour. The formulation is based on decomposition of the displacement/velocity field into deviatoric and volumetric parts. The volumetric part is iteratively eliminated without confronting locking or instability phenomenon. The iterative procedure is very cheap and simple to be implemented in any FEM code. Several examples are given to demonstrate the performance of the procedure. Copyright © 2004 John Wiley & Sons, Ltd.     

The development of hybrid axisymmetric elements based on the Hellinger–Reissner variational principle

We present a general procedure for the development of hybrid axisymmetric elements based on the Hellinger–Reissner principle within the context of linear elasticity. Similar to planar elements, the stress interpolation is obtained by an identification of the zero‐energy modes. We illustrate our procedure by designing a lower‐order (four‐node) and a higher‐order (nine‐node) element. Both elements are of correct rank, and moreover use the minimum number of stress parameters, namely seven and 17. Several examples are presented to show the excellent performance of both elements under various demanding situations such as when the material is almost incompressible, when the thickness to radius ratio is very small, etc. When the variation of the field variables is along the radial direction alone, when the mesh is uniform, and the loading is of pressure type, the developed elements are superconvergent, i.e. they yield the exact nodal displacement values. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite calculus formulation for incompressible solids using linear triangles and tetrahedra

Many finite elements exhibit the so‐called ‘volumetric locking’ in the analysis of incompressible or quasi‐incompressible problems.In this paper, a new approach is taken to overcome this undesirable effect. The starting point is a new setting of the governing differential equations using a finite calculus (FIC) formulation. The basis of the FIC method is the satisfaction of the standard equations for balance of momentum (equilibrium of forces) and mass conservation in a domain of finite size and retaining higher order terms in the Taylor expansions used to express the different terms of the differential equations over the balance domain. The modified differential equations contain additional terms which introduce the necessary stability in the equations to overcome the volumetric locking problem. The FIC approach has been successfully used for deriving stabilized finite element and meshless methods for a wide range of advective–diffusive and fluid flow problems. The same ideas are applied in this paper to derive a stabilized formulation for static and dynamic finite element analysis of incompressible solids using linear triangles and tetrahedra. Examples of application of the new stabilized formulation to linear static problems as well as to the semi‐implicit and explicit 2D and 3D non‐linear transient dynamic analysis of an impact problem and a bulk forming process are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

A stabilized nodally integrated tetrahedral

A stabilized, nodally integrated linear tetrahedral is formulated and analysed. It is well known that linear tetrahedral elements perform poorly in problems with plasticity, nearly incompressible materials, and acute bending. For a variety of reasons, low‐order tetrahedral elements are preferable to quadratic tetrahedral elements; particularly for nonlinear problems. But the severe locking problems of tetrahedrals have forced analysts to employ hexahedral formulations for most nonlinear problems. On the other hand, automatic mesh generation is often not feasible for building many 3D hexahedral meshes. A stabilized, nodally integrated linear tetrahedral is developed and shown to perform very well in problems with plasticity, nearly incompressible materials and acute bending. The formulation is analytically and numerically shown to be stable and optimally convergent for the compressible case provided sufficient smoothness of the exact solution u ∈ C² ∩ (H¹)³. Future work may extend the formulation to the incompressible regime and relax the regularity requirements; nonetheless, the results demonstrate that the method is not susceptible to locking and performs quite well in several standard linear and nonlinear benchmarks. Published in 2006 by John Wiley & Sons, Ltd.     

A meshfree‐enriched finite element method for compressible and near‐incompressible elasticity

In this paper, a two‐dimensional displacement‐based meshfree‐enriched FEM (ME‐FEM) is presented for the linear analysis of compressible and near‐incompressible planar elasticity. The ME‐FEM element is established by injecting a first‐order convex meshfree approximation into a low‐order finite element with an additional node. The convex meshfree approximation is constructed using the generalized meshfree approximation method and it possesses the Kronecker‐delta property on the element boundaries. The gradient matrix of ME‐FEM element satisfies the integration constraint for nodal integration and the resultant ME‐FEM formulation is shown to pass the constant stress test for the compressible media. The ME‐FEM interpolation is an element‐wise meshfree interpolation and is proven to be discrete divergence‐free in the incompressible limit. To prevent possible pressure oscillation in the near‐incompressible problems, an area‐weighted strain smoothing scheme incorporated with the divergence‐free ME‐FEM interpolation is introduced to provide the smoothing on strains and pressure. With this smoothed strain field, the discrete equations are derived based on a modified Hu–Washizu variational principle. Several numerical examples are presented to demonstrate the effectiveness of the proposed method for the compressible and near‐incompressible problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Displacement/pressure mixed interpolation in the method of finite spheres

The displacement‐based formulation of the method of finite spheres is observed to exhibit volumetric ‘locking’ when incompressible or nearly incompressible deformations are encountered. In this paper, we present a displacement/pressure mixed formulation as a solution to this problem. We analyse the stability and optimality of the formulation for several discretization schemes using numerical inf–sup tests. Issues concerning computational efficiency are also discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Consistent coupling of beam and shell models for thermo‐elastic analysis

In this paper, the finite element formulation of a transition element for consistent coupling between shell and beam finite element models of thin‐walled beam‐like structures in thermo‐elastic problems is presented. Thin‐walled beam‐like structures modelled only with beam elements cannot be used to study local stress concentrations or to provide local mechanical or thermal boundary conditions. For this purpose, the structure has to be modelled using shell elements. However, computations using shell elements are a lot more expensive as compared to beam elements. The finite element model can be more efficient when the shell elements are used only in regions where the local effects are to be studied or local boundary conditions have to be provided. The remaining part of the structure can be modelled with beam elements. To couple these two models (i.e. shell and beam models) at transitional cross‐sections, transition elements are derived here for thermo‐elastic problems. The formulation encloses large displacement and rotational behaviour, which is important in case of thin‐walled beam‐like structures. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical and spectral investigations of Trefftz infinite elements

The numerical and spectral performance of novel infinite elements for exterior problems of time‐harmonic acoustics are examined. The formulation is based on a functional which provides a general framework for domain‐based computation of exterior problems. Two prominent features simplify the task of discretization: the infinite elements mesh the interface only and need not match the finite elements on the interface. Various infinite element approximations for two‐dimensional configurations with circular interfaces are reviewed. Numerical results demonstrate the good performance of these schemes. A simple study points to the proper interpretation of spectral results for the formulation. The spectral properties of these infinite elements are examined with a view to the representation of physics and efficient numerical solution. Copyright © 1999 John Wiley & Sons, Ltd.     

Remeshing strategies for large deformation problems with frictional contact and nearly incompressible materials

We present a framework to efficiently solve large deformation contact problems with nearly incompressible materials by implementing adaptive remeshing. Specifically, nodally integrated elements are employed to avoid mesh locking when linear triangular or tetrahedral elements are used to facilitate mesh re‐generation. Solution variables in the bulk and on contact surfaces are transferred between meshes such that accuracy is maintained and re‐equilibration on the new mesh is possible. In particular, the displacement transfer in the bulk is accomplished through a constrained least squares problem based on nodal integration, while the transfer of contact tractions relies on parallel transport. Finally, a residual‐based error indicator is chosen to drive adaptive mesh refinement. The proposed strategies are applicable to both two‐dimensional or three‐dimensional analysis and are demonstrated to be robust by a number of numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Boundary element formulation for 3D transversely isotropic cracked bodies

The boundary traction integral representation is obtained in elasticity when the classical displacement representation is differentiated and combined according to Hooke's law. The use of both traction and displacement integral representations leads to a mixed (or dual) formulation of the BEM where the discretization effort for crack problems is much smaller than in the classical formulation. A boundary element analysis of three‐dimensional fracture mechanics problems of transversely isotropic solids based on the mixed formulation is presented in this paper. The hypersingular and strongly singular kernels appearing in the formulation are regularized by using two terms of the displacement series expansion and one term of the traction expansion, at the collocation point. All the remaining integrals are analytically evaluated or transformed by means of Stokes' theorem into regular or weakly singular integrals, which are numerically computed. The method is general and can be used for elements of any shape including quarter‐point crack front elements. No change of co‐ordinates is required for the integration. The formulation as presented in this paper is something as clear, general and easy to handle as the classical BE formulation. It is used in combination with three‐dimensional quadratic and quarter‐point elements to obtain accurate results for several different crack problems. Cracks in boundless and finite transversely isotropic domains are studied. The meshes are simple and include only discretization of the crack and the external boundary. The obtained results are in good agreement with those existing in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Computational electro‐elasticity and magneto‐elasticity for quasi‐incompressible media immersed in free space

In this work, a mixed variational formulation to simulate quasi‐incompressible electro‐active or magneto‐active polymers immersed in the surrounding free space is presented. A novel domain decomposition is used to disconnect the primary coupled problem and the arbitrary free‐space mesh update problem. Exploiting this decomposition, we describe a block‐iterative approach to solving the linearised multiphysics problem, and a physically and geometrically based, three‐parameter method to update the free space mesh. Several application‐driven example problems are implemented to demonstrate the robustness of the mixed formulation for both electro‐elastic and magneto‐elastic problems involving both finite deformations and quasi‐incompressible media. Copyright © 2016 John Wiley & Sons, Ltd.     

A particle method for two‐phase flows with compressible air pocket

When using particle methods to simulate water–air flows with compressible air pockets, a major challenge is to deal with the large differences in physical properties (e.g., density and viscosity) between water and air. In addition, the accurate modeling of air compressibility is essential. To this end, a new two‐phase strategy is proposed to simulate incompressible and compressible fluids simultaneously without iterations between the solvers for incompressible and compressible flows. Water is modeled by the recently developed 2‐phase Consistent Particle Method for incompressible flows. For air modeling, a new compressible solver is proposed based on the ideal gas law and thermodynamics. The formulation avoids the problem of determining the actual sound speed that is dependent on the temperature and is therefore not necessarily constant. In addition, the compressible air solver is seamlessly integrated with the incompressible solver 2‐phase Consistent Particle Method because they both use the same predictor–corrector scheme to solve the governing equations. The performance of the proposed method is demonstrated by three benchmark problems as well as an experimental study of sloshing impact with entrapped air pockets in an oscillating tank. Copyright © 2016 John Wiley & Sons, Ltd.     

Enrichment of enhanced assumed strain approximations for representing strong discontinuities: addressing volumetric incompressibility and the discontinuous patch test

We present a geometrically non‐linear assumed strain method that allows for the presence of arbitrary, intra‐finite element discontinuities in the deformation map. Special attention is placed on the coarse‐mesh accuracy of these methods and their ability to avoid mesh locking in the incompressible limit. Given an underlying mesh and an arbitrary failure surface, we first construct an enriched approximation for the deformation map with the non‐linear analogue of the extended finite element method (X‐FEM). With regard to the richer space of functions spanned by the gradient of the enriched approximation, we then adopt a broader interpretation of variational consistency for the construction of the enhanced strain. In particular, in those elements intersected by the failure surface, we construct enhanced strain approximations which are orthogonal to piecewise‐constant stress fields. Contrast is drawn with existing strong discontinuity approaches where the enhanced strain variations in localized elements were constructed to be orthogonal to constant nominal stress fields. Importantly, the present formulation gives rise to a symmetric tangent stiffness matrix, even in localized elements. The present modification also allows for the satisfaction of a discontinuous patch test, wherein two different constant stress fields (on each side of the failure surface) lie in the solution space. We demonstrate how the proposed modifications eliminate spurious stress oscillations along the failure surface, particularly for nearly incompressible material response. Additional numerical examples are provided to illustrate the efficacy of the modified method for problems in hyperelastic fracture mechanics. Copyright © 2003 John Wiley & Sons, Ltd.