A three‐dimensional element‐free Galerkin elastic and elastoplastic formulation

A small strain, three‐dimensional, elastic and elastoplastic Element‐Free Galerkin (EFG) formulation is developed. Singular weight functions are utilized in the Moving‐Least‐Squares (MLS) determination of shape functions and shape function derivatives allowing accurate, direct nodal imposition of essential boundary conditions. A variable domain of influence EFG method is introduced leading to increased efficiency in computing the MLS shape functions and their derivatives. The elastoplastic formulations are based on the consistent tangent operator approach and closely follow the incremental formulations for non‐linear analysis using finite elements. Several linear elastic and small strain elastoplastic numerical examples are presented to verify the accuracy of the numerical formulations. Copyright © 1999 John Wiley & Sons, Ltd.     

A novel non‐linearly explicit second‐order accurate L‐stable methodology for finite deformation: hypoelastic/hypoelasto‐plastic structural dynamics problems

A novel non‐linearly explicit second‐order accurate L‐stable computational methodology for integrating the non‐linear equations of motion without non‐linear iterations during each time step, and the underlying implementation procedure is described. Emphasis is placed on illustrative non‐linear structural dynamics problems employing both total/updated Lagrangian formulations to handle finite deformation hypoelasticity/hypoelasto‐plasticity models in conjunction with a new explicit exact integration procedure for a particular rate form constitutive equation. Illustrative numerical examples are shown to demonstrate the robustness of the overall developments for non‐linear structural dynamics applications. Copyright © 2004 John Wiley & Sons, Ltd.     

Stabilization of mixed tetrahedral elements at large deformations

This paper presents stabilized mixed finite element formulations for tetrahedral elements at large deformations using volume and area bubble functions. To this end, the corresponding weak formulations are derived for the standard two‐field method, the method of incompatible modes and the enhanced strain method. Then, the weak formulations will be linearized. Furthermore, the matrix formulations for the weak formulations and its linearizations are summarized. The numerical results for incompressible rubber‐like materials using a Neo‐Hookean material law show the locking‐free performance and the drastic damping of the stresses for the new stabilized tetrahedral elements in finite deformation problems. This paper is an extension of the works published by the authors regarding small deformation problems for linear elasticity and plasticity. Copyright © 2011 John Wiley & Sons, Ltd.     

Higher‐order accurate time‐step‐integration algorithms by post‐integration techniques

In this paper, a framework to construct higher‐order‐accurate time‐step‐integration algorithms based on the post‐integration techniques is presented. The prescribed initial conditions are naturally incorporated in the formulations and can be strongly or weakly enforced. The algorithmic parameters are chosen such that unconditionally A‐stable higher‐order‐accurate time‐step‐integration algorithms with controllable numerical dissipation can be constructed for linear problems. Besides, it is shown that the order of accuracy for non‐linear problems is maintained through the relationship between the present formulation and the Runge–Kutta method. The second‐order differential equations are also considered. Numerical examples are given to illustrate the validity of the present formulation. Copyright © 2001 John Wiley & Sons, Ltd.     

High‐order Higdon‐like boundary conditions for exterior transient wave problems

Recently developed non‐reflecting boundary conditions are applied for exterior time‐dependent wave problems in unbounded domains. The linear time‐dependent wave equation, with or without a dispersive term, is considered in an infinite domain. The infinite domain is truncated via an artificial boundary ??, and a high‐order non‐reflecting boundary condition (NRBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. The new boundary scheme is based on a reformulation of the sequence of NRBCs proposed by Higdon. We consider here two reformulations: one that involves high‐order derivatives with a special discretization scheme, and another that does not involve any high derivatives beyond second order. The latter formulation is made possible by introducing special auxiliary variables on ??. In both formulations the new NRBCs can easily be used up to any desired order. They can be incorporated in a finite element or a finite difference scheme; in the present paper the latter is used. In contrast to previous papers using similar formulations, here the method is applied to a fully exterior two‐dimensional problem, with a rectangular boundary. Numerical examples in infinite domains are used to demonstrate the performance and advantages of the new method. In the auxiliary‐variable formulation long‐time corner instability is observed, that requires special treatment of the corners (not addressed in this paper). No such difficulties arise in the high‐derivative formulation. Published in 2005 by John Wiley & Sons, Ltd.     

Iterative solvers for 3D linear and nonlinear elasticity problems: Displacement and mixed formulations

We present new iterative solvers for large‐scale linear algebraic systems arising from the finite element discretization of the elasticity equations. We focus on the numerical solution of 3D elasticity problems discretized by quadratic tetrahedral finite elements and we show that second‐order accuracy can be obtained at very small overcost with respect to first‐order (linear) elements. Different Krylov subspace methods are tested on various meshes including elements with small aspect ratio. We first construct a hierarchical preconditioner for the displacement formulation specifically designed for quadratic discretizations. We then develop efficient tools for preconditioning the 2 × 2 block symmetric indefinite linear system arising from mixed (displacement‐pressure) formulations. Finally, we present some numerical results to illustrate the potential of the proposed methods. Copyright © 2010 John Wiley & Sons, Ltd.     

Self-consistent boundary conditions with blurred derivatives

It is shown in this paper that self-consistent boundary conditions for numerical methods based on blurred derivatives can be derived from a suitable change of variables of the fundamental blurred approximation of the differential equation, followed by application of Leibnitz theorem for differentiation of an integral. The simplest scheme obtained in this way resembles the weak Local Petrov-Galerkin approximation, although interpretation of the operators appearing in the final equations is quite different—as is the derivation itself. Subsequent transformation leads to integral equations similar to the starting point for boundary integral methods of solution. In this way, a number of well-known computational methods are shown to be derivable from adequate manipulation of the blurred derivative technique. However, other approximations, which are not derivable with standard methods can also be obtained, hinting at a greater generality of blurred derivatives.     

Dispersion analysis and element‐free Galerkin solutions of second‐ and fourth‐order gradient‐enhanced damage models

Gradient‐dependent damage formulations incorporate higher‐order derivatives of state variables in the constitutive equations. Different formulations have been derived for this gradient enhancement, comparison of which is difficult in a finite element context due to higher‐order continuity requirements for certain formulations. On the other hand, the higher‐order continuity requirements are met naturally by element‐free Galerkin (EFG) shape functions. Thus, the EFG method provides a suitable tool for the assessment of gradient enhanced continuum models. Dispersion analyses have been carried out to compare different gradient enhanced models with the non‐local damage model. The formulation of the additional boundary conditions is addressed. Numerical examples show the objectivity with respect to the discretization and the differences between various gradient formulations with second‐ and fourth‐order derivatives. It is shown that with the same underlying internal length scale, very different results can be obtained. Copyright © 2000 John Wiley & Sons, Ltd.     

Computational mechanics based on the theory of boundary eigensolutions

The theory of boundary eigensolutions for boundary value problems is applied to the development of computational mechanics formulations. The boundary element and finite element methods that result are consistent with the mathematical theory of boundary value problems. Although the approach is quite general, this paper focuses on potential problems. For these problems, both methods employ potential and boundary flux as primary variables. Convergence characteristics of the new flux‐oriented finite element method are also developed. By utilizing suitable boundary weight functions, the formulations are written exclusively in terms of bounded quantities, even for non‐smooth problems involving notches, cracks and mixed boundary conditions. The results of numerical experiments indicate that the algorithms perform in concert with the underlying theory and thus provide an attractive alternative to existing approaches. Beyond this, the approach developed here provides a new perspective from which to view computational mechanics, and can be used to obtain a better understanding of boundary element and finite element methods. Comparisons with closed‐form boundary eigensolutions are also presented in order to provide a means for assessing the numerical methods. Copyright © 2001 John Wiley & Sons, Ltd.     

Linear smoothed extended finite element method

The extended finite element method was introduced in 1999 to treat problems involving discontinuities with no or minimal remeshing through appropriate enrichment functions. This enables elements to be split by a discontinuity, strong or weak, and hence requires the integration of discontinuous functions or functions with discontinuous derivatives over elementary volumes. A variety of approaches have been proposed to facilitate these special types of numerical integration, which have been shown to have a large impact on the accuracy and the convergence of the numerical solution. The smoothed extended finite element method (XFEM), for example, makes numerical integration elegant and simple by transforming volume integrals into surface integrals. However, it was reported in the literature that the strain smoothing is inaccurate when non‐polynomial functions are in the basis. In this paper, we investigate the benefits of a recently developed Linear smoothing procedure which provides better approximation to higher‐order polynomial fields in the basis. Some benchmark problems in the context of linear elastic fracture mechanics are solved and the results are compared with existing approaches. We observe that the stress intensity factors computed through the proposed linear smoothed XFEM is more accurate than that obtained through smoothed XFEM.     

Automatic solver for non‐linear partial differential equations with implicit local laws: Application to unilateral contact

In general, non‐linear continuum mechanics combine global balance equations and local constitutive laws. In this work, frictionless contact between a rigid tool and a thin elastic shell is considered. This class of boundary value problems involves two non‐linear algebraic laws: the first one gives explicitly the stress field as a function of the strain throughout the continuum part, whereas the second one is a non‐linear equation relating the contact forces and the displacement at the boundary.Given the fact that classical computational approaches sometimes require significant effort in implementation of complex non‐linear problems, a computation technique based on automatic differentiation of constitutive laws is presented in this paper. The procedure enables to compute automatically the higher‐order derivatives of these constitutive laws and thereafter to define the Taylor series that are the basis of the continuation technique called asymptotic numerical method. The algorithm is about the same with an explicit or implicit constitutive relation. In the modelling of forming processes, many tool shapes can be encountered. The presented computational technique permits an easy implementation of these complex surfaces, for instance in a finite element code: the user is only required to define the tool geometry and the computer is able to obtain the higher‐order derivatives. Copyright © 2013 John Wiley & Sons, Ltd.     

Exploiting partial or complete geometrical symmetry in 3D symmetric Galerkin indirect BEM formulations

Procedures based on group representation theory, allowing the exploitation of geometrical symmetry in symmetric Galerkin BEM formulations, are investigated. In particular, this investigation is based on the weaker assumption of partial geometrical symmetry, where the boundary has two disconnected components, one of which is symmetric; e.g. this can be very useful for defect identification problems. The main development is expounded in the context of 3D Neumann elastostatic problems, considered as model problems; and then extended to SGBIE formulations for Dirichlet and/or scalar problems. Both Abelian and non‐Abelian finite symmetry groups are considered. The effectiveness of the present approach is demonstrated through numerical examples, where both partial and complete symmetry are considered, in connection with both Abelian and non‐Abelian symmetry groups. Copyright © 2003 John Wiley & Sons, Ltd.     

Hamilton's law of variable mass system and time finite element formulations for time‐varying structures based on the law

Traditional principles of mechanics are primarily conceived for constant mass systems, which are only valid if mass is gained or lost at null velocity with respect to an inertial reference frame for variable mass systems, thus the numerical algorithms for time‐varying structures based on these principles are only suitable for this special case. In this paper, Hamilton's law of variable mass system is derived based on Meshchersky's fundamental equation, and two classes of novel time finite element formulations for linear systems with arbitrary continuous time‐varying parameters are developed based on the previous law. The formulations are verified extensively through numerical examples in which the convergence and effectiveness of algorithms are evaluated. Numerical examples demonstrate that compared with the algorithms for time‐varying structures that developed based on traditional principles of mechanics, the proposed algorithms provide extended capabilities in both time‐varying mass problems that mass is gained or lost at any velocity (such as rocket problem) and moving‐mass problems (such as vehicle‐bridge interaction problem) besides the time‐varying stiffness and damping problems, the proposed algorithms have a wider range of application. In particular, Hamilton's law of variable mass system provides a solid theoretical foundation for further research on the algorithm design for time‐varying structures. Copyright © 2014 John Wiley & Sons, Ltd.     

Discontinuous Galerkin methods for non‐linear elasticity

This paper presents the formulation and a partial analysis of a class of discontinuous Galerkin methods for quasistatic non‐linear elasticity problems. These methods are endowed with several salient features. The equations that define the numerical scheme are the Euler–Lagrange equations of a one‐field variational principle, a trait that provides an elegant and simple derivation of the method. In consonance with general discontinuous Galerkin formulations, it is possible within this framework to choose different numerical fluxes. Numerical evidence suggests the absence of locking at near‐incompressible conditions in the finite deformations regime when piecewise linear elements are adopted. Finally, a conceivable surprising characteristic is that, as demonstrated with numerical examples, these methods provide a given accuracy level for a comparable, and often lower, computational cost than conforming formulations. Stabilization is occasionally needed for discontinuous Galerkin methods in linear elliptic problems. In this paper we propose a sufficient condition for the stability of each linearized non‐linear elastic problem that naturally includes material and geometric parameters; the latter needed to account for buckling. We then prove that when a similar condition is satisfied by the discrete problem, the method provides stable linearized deformed configurations upon the addition of a standard stabilization term. We conclude by discussing the complexity of the implementation, and propose a computationally efficient approach that avoids looping over both elements and element faces. Several numerical examples are then presented in two and three dimensions that illustrate the performance of a selected discontinuous Galerkin method within the class. Copyright © 2006 John Wiley & Sons, Ltd.     

Methods for connecting dissimilar three‐dimensional finite element meshes

Two methods are presented for connecting dissimilar three‐dimensional finite element meshes. The first method combines the concept of master and slave surfaces with the uniform strain approach for finite elements. By modifying the boundaries of elements on a slave surface, corrections are made to element formulations such that first‐order patch tests are passed. The second method is based entirely on constraint equations, but only passes a weaker form of the patch test for non‐planar surfaces. Both methods can be used to connect meshes with different element types. In addition, master and slave surfaces can be designated independently of relative mesh resolutions. Example problems in three‐dimensional linear elasticity are presented. Published in 2000 by John Wiley & Sons, Ltd.     

3D modelling of strong discontinuities in elastoplastic solids: fixed and rotating localization formulations

This paper is concerned with the incorporation of displacement discontinuities into a continuum theory of elastoplasticity for the modelling of localization processes such as cracking in brittle materials. Based on the strong discontinuity approach (SDA) (Computational Mechanics 1993; 12: 277–296) mesh objective 2D and 3D finite element formulations are developed using linear and quadratic 2D elements as well as 8‐noded 3D elements. In the formulation of the finite‐element model proposed in the paper, the analogy with standard formulations is emphasized. The parameter defining the amplitude of the displacement jump within the finite element is condensed out at the material level without employing the standard static condensation technique. This approach results in linearized constitutive equations formally identical to continuum models. Therefore, the standard return mapping algorithm is used to solve the non‐linear equations. In analogy to concepts used in continuum smeared crack models, a rotating formulation of the SDA is proposed in addition to the standard concept of fixed discontinuities. It is shown that the rotating localization approach reduces locking effects observed in analyses based on fixed localization directions. The applicability of the proposed SDA finite‐element model as well as its numerical performance is investigated by means of a three‐dimensional ultimate load analysis of a steel anchor embedded in a concrete block subjected to a shear force. Copyright © 2003 John Wiley & Sons, Ltd.     

Optimal layout design of three‐dimensional geometrically non‐linear structures using the element connectivity parameterization method

The topology design optimization of ‘three‐dimensional geometrically‐non‐linear’ continuum structures is still difficult not only because of the size of the problem but also because of the unstable continuum finite elements that arise during the optimization. To overcome these difficulties, the element connectivity parameterization (ECP) method with two implementation formulations is proposed. In ECP, structural layouts are represented by inter‐element connectivity, which is controlled by the stiffness of element‐connecting zero‐length links. Depending on the link location, ECP may be classified as an external ECP (E‐ECP) or an internal ECP (I‐ECP). In this paper, I‐ECP is newly developed to substantially enhance computational efficiency. The main idea in I‐ECP is to reduce system matrix size by eliminating some internal degrees of freedom associated with the links at voxel level. As for ECP implementation with commercial software, E‐ECP, developed earlier for two‐dimensional problems, is easier to use even for three‐dimensional problems because it requires only numerical analysis results for design sensitivity calculation. The characteristics of the I‐ECP and E‐ECP methods are compared, and these methods are validated with numerical examples. Copyright © 2006 John Wiley & Sons, Ltd.     

Hybrid and enhanced finite element methods for problems of soil consolidation

Hybrid and enhanced finite element methods with bi‐linear interpolations for both the solid displacements and the pore fluid pressures are derived based on mixed variational principles for problems of elastic soil consolidation. Both plane strain and axisymmetric problems are studied. It is found that by choosing appropriate interpolation of enhanced strains in the enhanced method, and by choosing appropriate interpolations of strains, effective stresses and enhanced strains in the hybrid method, the oscillations of nodal pore pressures can be eliminated. Several numerical examples demonstrating the capability and performance of the enhanced and hybrid finite element methods are presented. It is also shown that for some situations, such as problems involving high Poisson's ratio and in other related problems where bending effects are evident, the performance of the enhanced and hybrid methods are superior to that of the conventional displacement‐based method. The results from the hybrid method are better than those from the enhanced method for some situations, such as problems in which soil permeability is variable or discontinuous within elements. Since all the element parameters except the nodal displacements and nodal pore pressures are assumed in the element level and can be eliminated by static condensation, the implementations of the enhanced method and the hybrid method are basically the same as the conventional displacement‐based finite element method. The present enhanced method and hybrid method can be easily extended to non‐linear consolidation problems. Copyright © 2006 John Wiley & Sons, Ltd.