Local and non‐local thermomechanical modeling of elastic‐plastic materials undergoing large strains

This paper deals with the numerical analysis of instabilities for elastic‐plastic materials undergoing large deformations in non‐isothermal conditions. The considered isotropic model is fully thermomechanically coupled and includes temperature‐induced softening, which is another source of strain localization next to geometrical effects. Due to complexity of the model, a symbolic‐numerical tool Ace is used for the preparation of user‐supplied subroutines for the finite element method. The computational verification is performed using two benchmark tests: necking of circular bar in tension and shear banding of elongated rectangular plate in plain strain conditions. The attention is focused on mesh dependence of the numerical results and the regularizing effect of heat conduction. The research reveals that the conductivity influences the shear band width and ductility of the material response; however, for the adiabatic case, the results are discretization sensitive, and another regularization is needed. A new gradient‐enhanced thermomechanical model is developed that introduces an internal length parameter governing the size of the shear band caused by thermal softening. The numerical verification of the non‐local model is performed for the adiabatic case. Subsequently, the simultaneous application of the gradient enhancement and heat conduction in the model is analyzed, which reproduces an evolving shear band. Copyright © 2016 John Wiley & Sons, Ltd.     

Mixed finite element method for coupled thermo‐hydro‐mechanical process in poro‐elasto‐plastic media at large strains

A mixed finite element for coupled thermo‐hydro‐mechanical (THM) analysis in unsaturated porous media is proposed. Displacements, strains, the net stresses for the solid phase; pressures, pressure gradients, Darcy velocities for pore water and pore air phases; temperature, temperature gradients, the total heat flux are interpolated as independent variables. The weak form of the governing equations of coupled THM problems in porous media within the element is given on the basis of the Hu–Washizu three‐filed variational principle. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of the THM constitutive model for unsaturated porous media based on the CAP model. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elasto‐plastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is utilized. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization and the softening behaviours caused by thermal and chemical effects. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐local boundary integral formulation for softening damage

A strongly non‐local boundary element method (BEM) for structures with strain‐softening damage treated by an integral‐type operator is developed. A plasticity model with yield limit degradation is implemented in a boundary element program using the initial‐stress boundary element method with iterations in each load increment. Regularized integral representations and boundary integral equations are used to avoid the difficulties associated with numerical computation of singular integrals. A numerical example is solved to verify the physical correctness and efficiency of the proposed formulation. The example consists of a softening strip perforated by a circular hole, subjected to tension. The strain‐softening damage is described by a plasticity model with a negative hardening parameter. The local formulation is shown to exhibit spurious sensitivity to cell mesh refinements, localization of softening damage into a band of single‐cell width, and excessive dependence of energy dissipation on the cell size. By contrast, the results for the non‐local theory are shown to be free of these physically incorrect features. Compared to the classical non‐local finite element approach, an additional advantage is that the internal cells need to be introduced only within the small zone (or band) in which the strain‐softening damage tends to localize within the structure. Copyright © 2003 John Wiley & Sons, Ltd.     

Dynamic inflation of non‐linear elastic and viscoelastic rubber‐like membranes

The present paper deals with the dynamic inflation of rubber‐like membranes.The material is assumed to obey the hyperelastic Mooney's model or the non‐linear viscoelastic Christensen's model. The governing equations of free inflation are solved by a total Lagrangian finite element method for the spatial discretization and an explicit finite‐difference algorithm for the time‐integration scheme. The numerical implementation of constitutive equations is highlighted and the special case of integral viscoelastic models is examined in detail. The external force consists in a gas flow rate, which is more realistic than a pressure time history. Then, an original method is used to calculate the pressure evolution inside the bubble depending on the deformation state. Our numerical procedure is illustrated through different examples and compared with both analytical and experimental results. These comparisons yield good agreement and show the ability of our approach to simulate both stable and unstable large strain inflations of rubber‐like membranes. Copyright © 2001 John Wiley & Sons, Ltd.     

Computational simulation of chemical dissolution‐front instability in fluid‐saturated porous media under non‐isothermal conditions

This paper primarily deals with the computational aspects of chemical dissolution‐front instability problems in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. After the dimensionless governing partial differential equations of the non‐isothermal chemical dissolution‐front instability problem are briefly described, the formulation of a computational procedure, which contains a combination of using the finite difference and finite element method, is derived for simulating the morphological evolution of chemical dissolution fronts in the non‐isothermal chemical dissolution system within two‐dimensional fluid‐saturated porous media. To ensure the correctness and accuracy of the numerical solutions, the proposed computational procedure is verified through comparing the numerical solutions with the analytical solutions for a benchmark problem. As an application example, the verified computational procedure is then used to simulate the morphological evolution of chemical dissolution fronts in the supercritical non‐isothermal chemical dissolution system. The related numerical results have demonstrated the following: (1) the proposed computational procedure can produce accurate numerical solutions for the planar chemical dissolution‐front propagation problem in the non‐isothermal chemical dissolution system consisting of a fluid‐saturated porous medium; (2) the Zhao number has a significant effect not only on the dimensionless propagation speed of the chemical dissolution front but also on the distribution patterns of the dimensionless temperature, dimensionless pore‐fluid pressure, and dimensionless chemical‐species concentration in a non‐isothermal chemical dissolution system; (3) once the finger penetrates the whole computational domain, the dimensionless pore‐fluid pressure decreases drastically in the non‐isothermal chemical dissolution system.     

Non‐linear analysis of thin‐walled members of open cross‐section

The paper presents a means of determining the non‐linear stiffness matrices from expressions for the first and second variation of the Total Potential of a thin‐walled open section finite element that lead to non‐linear stiffness equations. These non‐linear equations can be solved for moderate to large displacements. The variations of the Total Potential have been developed elsewhere by the authors, and their contribution to the various non‐linear matrices is stated herein. It is shown that the method of solution of the non‐linear stiffness matrices is problem dependent. The finite element procedure is used to study non‐linear torsion that illustrates torsional hardening, and the Newton–Raphson method is deployed for this study. However, it is shown that this solution strategy is unsuitable for the second example, namely that of the post‐buckling response of a cantilever, and a direct iteration method is described. The good agreement for both of these problems with the work of independent researchers validates the non‐linear finite element method of analysis. Copyright © 1999 John Wiley & Sons, Ltd.     

A stabilized finite element formulation for monolithic thermo‐hydro‐mechanical simulations at finite strain

An adaptively stabilized monolithic finite element model is proposed to simulate the fully coupled thermo‐hydro‐mechanical behavior of porous media undergoing large deformation. We first formulate a finite‐deformation thermo‐hydro‐mechanics field theory for non‐isothermal porous media. Projection‐based stabilization procedure is derived to eliminate spurious pore pressure and temperature modes due to the lack of the two‐fold inf‐sup condition of the equal‐order finite element. To avoid volumetric locking due to the incompressibility of solid skeleton, we introduce a modified assumed deformation gradient in the formulation for non‐isothermal porous solids. Finally, numerical examples are given to demonstrate the versatility and efficiency of this thermo‐hydro‐mechanical model. Copyright © 2015 John Wiley & Sons, Ltd.     

A simple finite element‐based framework for the analysis of elastic pseudo‐rigid bodies

This article advocates a general procedure for the numerical investigation of pseudo‐rigid bodies. The equations of motion for pseudo‐rigid bodies are shown to be mathematically equivalent to those corresponding to certain constant‐strain finite element approximations for general deformable continua. A straightforward algorithmic implementation is achieved in a classical finite element framework. Also, a penalty formulation is suggested for modelling contact between pseudo‐rigid bodies. Representative planar simulations using a non‐linear elastic model demonstrate the predictive capacity of the pseudo‐rigid theory, as well as the robustness of the proposed computational procedure. Copyright © 1999 John Wiley & Sons, Ltd.     

Mixed finite element method for saturated poroelastoplastic media at large strains

A mixed finite element for hydro‐dynamic analysis in saturated porous media in the frame of the Biot theory is proposed. Displacements, effective stresses, strains for the solid phase and pressure, pressure gradients, and Darcy velocities for the fluid phase are interpolated as independent variables. The weak form of the governing equations of coupled hydro‐dynamic problems in saturated porous media within the element are given on the basis of the Hu–Washizu three‐field variational principle. In light of the stabilized one point quadrature super‐convergent element developed in solid continuum, the interpolation approximation modes for the primary unknowns and their spatial derivatives of the solid and the fluid phases within the element are assumed independently. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of pressure‐dependent non‐associated plasticity. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elastoplastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is used. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization due to strain softening in poroelastoplastic media subjected to dynamic loading at large strain. Copyright © 2003 John Wiley & Sons, Ltd.     

Energy and momentum conserving elastodynamics of a non‐linear brick‐type mixed finite shell element

The paper presents aspects of the finite element formulation of momentum and energy conserving algorithms for the non‐linear dynamic analysis of shell‐like structures. The key contribution is a detailed analysis of the implementation of a Simó–Tarnow‐type conservation scheme in a recently developed new mixed finite shell element. This continuum‐based shell element provides a well‐defined interface to strain‐driven constitutive stress updates algorithms. It is based on the classic brick‐type trilinear displacement element and is equipped with specific gradient‐type enhanced strain modes and shell‐typical assumed strain modifications. The excellent performance of the proposed dynamic shell formulation with respect to conservation properties and numerical stability behaviour is demonstrated by means of three representative numerical examples of elastodynamics which exhibit complex free motions of flexible structures undergoing large strains and large rigid‐body motions. Copyright © 2001 John Wiley & Sons, Ltd.     

A particle‐in‐cell solution to the silo discharging problem

The problem of flow of a granular material during the process of discharging a silo is considered in the present paper. The mechanical behaviour of the material is described by the use of the model of the elastic–plastic solid with the Drucker–Prager yield condition and the non‐associative flow rule. The phenomenon of friction between the stored material and the silo walls is taken into account—the Coulomb model of friction is used in the analysis. The problem is analysed by means of the particle‐in‐cell method—a variant of the finite element method which enables to solve the pertinent equations of motion on an arbitrary computational mesh and trace state variables at points of the body chosen independently of the mesh. The method can be regarded as an arbitrary Lagrangian–Eulerian formulation of the finite element method, and overcomes the main drawback of the updated Lagrangian formulation of FEM related to mesh distortion. The entire process of discharging a silo can be analysed by this approach. The dynamic problem is solved by the use of the explicit time‐integration scheme. Several numerical examples are included. The plane strain and axisymmetric problems are solved for silos with flat bottoms and conical hoppers. Some results are compared with experimental ones. Copyright © 1999 John Wiley & Sons, Ltd.     

Moving‐least‐squares‐particle hydrodynamics—I. Consistency and stability

The Smooth‐Particle‐Hydrodynamics (SPH) method is derived in a novel manner by means of a Galerkin approximation applied to the Lagrangian equations of continuum mechanics as in the finite‐element method. This derivation is modified to replace the SPH interpolant with the Moving‐Least‐Squares (MLS) interpolant of Lancaster and Saulkaskas, and define a new particle volume which ensures thermodynamic compatibility. A variable‐rank modification of the MLS interpolants which retains their desirable summation properties is introduced to remove the singularities that occur when divergent flow reduces the number of neighbours of a particle to less than the minimum required. A surprise benefit of the Galerkin SPH derivation is a theoretical justification of a common ad hoc technique for variable‐h SPH. The new MLSPH method is conservative if an anti‐symmetric quadrature rule for the stiffness matrix elements can be supplied. In this paper, a simple one‐point collocation rule is used to retain similarity with SPH, leading to a non‐conservative method. Several examples document how MLSPH renders dramatic improvements due to the linear consistency of its gradients on three canonical difficulties of the SPH method: spurious boundary effects, erroneous rates of strain and rotation and tension instability. Two of these examples are non‐linear Lagrangian patch tests with analytic solutions with which MLSPH agrees almost exactly. The examples also show that MLSPH is not absolutely stable if the problems are run to very long times. A linear stability analysis explains both why it is more stable than SPH and not yet absolutely stable and an argument is made that for realistic dynamic problems MLSPH is stable enough. The notion of coherent particles, for which the numerical stability is identical to the physical stability, is introduced. The new method is easily retrofitted into a generic SPH code and some observations on performance are made. Copyright © 1999 John Wiley & Sons, Ltd.     

Explicit finite element artificial boundary scheme for transient scalar waves in two‐dimensional unbounded waveguide

To simulate the transient scalar wave propagation in a two‐dimensional unbounded waveguide, an explicit finite element artificial boundary scheme is proposed, which couples the standard dynamic finite element method for complex near field and a high‐order accurate artificial boundary condition (ABC) for simple far field. An exact dynamic‐stiffness ABC that is global in space and time is constructed. A temporal localization method is developed, which consists of the rational function approximation in the frequency domain and the auxiliary variable realization into time domain. This method is applied to the dynamic‐stiffness ABC to result in a high‐order accurate ABC that is local in time but global in space. By discretizing the high‐order accurate ABC along artificial boundary and coupling the result with the standard lumped‐mass finite element equation of near field, a coupled dynamic equation is obtained, which is a symmetric system of purely second‐order ordinary differential equations in time with the diagonal mass and non‐diagonal damping matrices. A new explicit time integration algorithm in structural dynamics is used to solve this equation. Numerical examples are given to demonstrate the effectiveness of the proposed scheme. Copyright © 2011 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.     

Three‐dimensional coupled heat,moisture, and air transfer in a deformable unsaturated soil

A three‐dimensional numerical model is presented for three‐phase flow (moisture, air, and heat) in a deformable partly saturated soil with deformation calculated via a non‐linear elastic theory. The present work is an extension of a two‐dimensional analysis presented by Thomas and He. The objective of this work is the solution of problems of greater geometric complexity. The mathematical formulation of this coupled problem consists of four governing equations, developed from the principles of mass and energy conservations as well as the stress equilibrium equation. Darcy's flow law is used to describe the motion of liquid and air in the porous medium, and a Philip and de Vries type vapour flow approach is employed in the formulation. A Galerkin finite element method coupled with a finite difference recurrence relationship is used to obtain simultaneous solutions to the governing equations where pore liquid, pore air pressures, temperature and displacements are the primary variables. The method allows the non‐linear nature of the soil parameters to be modelled. Three‐dimensional 20‐noded isoparametric elements are used to simulate different types of cases for the verification of the work. Results are presented of the application of the new model to four problems, two of which are isothermal and two heating simulations. The three‐dimensional nature of the results achieved is highlighted. Copyright © 1999 John Wiley & Sons, Ltd.     

A hierarchical finite element for geometrically non‐linear vibration of doubly curved,moderately thick isotropic shallow shells

A p‐version, hierarchical finite element for doubly curved, moderately thick, isotropic shallow shells is derived and geometrically non‐linear free vibrations of panels with rectangular planform are investigated. The geometrical non‐linearity is due to large displacements, and the effects of the rotatory inertia and transverse shear are considered. The time domain equations of motion are obtained by applying the principle of virtual work and the d'Alembert's principle. These equations are mapped to the frequency domain by the harmonic balance method, and are finally solved by a predictor–corrector method. The convergence properties of the element proposed and the influence of several parameters on the dynamic response are studied. These parameters are the shell's thickness, the width‐to‐length ratio, the curvature‐to‐width ratio and the ratio between curvature radii. The first and higher order modes are analysed. Some results are compared with results published or calculated using a commercial finite element package. It is demonstrated that with the proposed element low‐dimensional, accurate models are obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part II: 3‐D non‐linear behaviour

In a companion paper, the effects of approximations in the flexural‐torsional stability analysis of beams was studied, and it was shown that a second‐order rotation matrix was sufficiently accurate for a flexural‐torsional stability analysis. However, the second‐order rotation matrix is not necessarily accurate in formulating finite element model for a 3‐D non‐linear analysis of thin‐walled beams of open cross‐section. The approximations in the second‐order rotation matrix may introduce ‘self‐straining’ due to superimposed rigid‐body motions, which may lead to physically incorrect predictions of the 3‐D non‐linear behaviour of beams. In a 3‐D non‐linear elastic–plastic analysis, numerical integration over the cross‐section is usually used to check the yield criterion and to calculate the stress increments, the stress resultants, the elastic–plastic stress–strain matrix and the tangent modulus matrix. A scheme of the arrangement of sampling points over the cross‐section that is not consistent with the strain distributions may lead to incorrect predictions of the 3‐D non‐linear elastic–plastic behaviour of beams. This paper investigates the effects of approximations on the 3‐D non‐linear analysis of beams. It is found that a finite element model for 3‐D non‐linear analysis based on the second‐order rotation matrix leads to over‐stiff predictions of the flexural‐torsional buckling and postbuckling response and to an overestimate of the maximum load‐carrying capacities of beams in some cases. To perform a correct 3‐D non‐linear analysis of beams, an accurate model of the rotations must be used. A scheme of the arrangement of sampling points over the cross‐section that is consistent with both the longitudinal normal and shear strain distributions is needed to predict the correct 3‐D non‐linear elastic–plastic behaviour of beams. Copyright © 2001 John Wiley & Sons, Ltd.     

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

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

Finite element analysis of time‐dependent semi‐infinite wave‐guides with high‐order boundary treatment

A new finite element (FE) scheme is proposed for the solution of time‐dependent semi‐infinite wave‐guide problems, in dispersive or non‐dispersive media. The semi‐infinite domain is truncated via an artificial boundary ??, and a high‐order non‐reflecting boundary condition (NRBC), based on the Higdon non‐reflecting operators, is developed and applied on ??. The new NRBC does not involve any high derivatives beyond second order, but its order of accuracy is as high as one desires. It involves some parameters which are chosen automatically as a pre‐process. A C⁰ semi‐discrete FE formulation incorporating this NRBC is constructed for the problem in the finite domain bounded by ??. Augmented and split versions of this FE formulation are proposed. The semi‐discrete system of equations is solved by the Newmark time‐integration scheme. Numerical examples concerning dispersive waves in a semi‐infinite wave guide are used to demonstrate the performance of the new method. Copyright © 2003 John Wiley & Sons, Ltd.