Floating stress‐point integration meshfree method for large deformation analysis of elastic and elastoplastic materials

A new meshfree formulation of stress‐point integration, called the floating stress‐point integration meshfree method, is proposed for the large deformation analysis of elastic and elastoplastic materials. This method is a Galerkin meshfree method with an updated Lagrangian procedure and a quasi‐implicit time‐advancing scheme without any background cell for domain integration. Its new formulation is based on incremental equilibrium equations derived from the incremental virtual work equation, which is not generally used in meshfree formulations. Hence, this technique allows the temporal continuity of the mechanical equilibrium to be naturally achieved. The details of the new formulation and several examples of the large deformation analysis of elastic and elastoplastic materials are presented to show the validity and accuracy of the proposed method in comparison with those of the finite element method. Copyright © 2012 John Wiley & Sons, Ltd.     

A comparative study of stress update algorithms for rate‐independent and rate‐dependent crystal plasticity

The paper presents a comparative discussion of stress update algorithms for single‐crystal plasticity at small strains. The key result is a new unified fully implicit multisurface‐type return algorithm for both the rate‐independent and the rate‐dependent setting, endowed with three alternative approaches to the regularization of possible redundant slip activities. The fundamental problem of the rate‐independent theory is the possible ill condition due to linear‐dependent active slip systems. We discuss three possible algorithmic approaches to deal with this problem. This includes the use of alternative generalized inverses of the Jacobian of the currently active yield criterion functions as well as a new diagonal shift regularization technique, motivated by a limit of the rate‐dependent theory. Analytical investigations and numerical experiments show that all three approaches result in similar physically acceptable predictions of the active slip of rate‐independent single‐crystal plasticity, while the new proposed diagonal shift method is the most simple and efficient concept. Copyright © 2001 John Wiley & Sons, Ltd.     

Implicit integration and consistent tangent modulus of a time‐dependent non‐unified constitutive model

This paper describes the implicit integration and consistent tangent modulus of an inelastic constitutive model with transient and steady strain rates, both of which are time‐ and temperature‐dependent; the transient rate is influenced by the evolution of back stress decomposed into parts, while the steady rate depends only on applied stress and temperature. Such a non‐unified model is useful for high‐temperature structural analysis and is practical owing to the ease in determining material constants. The implicit integration is shown to result in two scalar‐valued coupled equations, and the consistent tangent modulus is derived in a quite versatile form by introducing a set of fourth‐rank constitutive parameters into the discretized evolution rule of back stress. The constitutive model is, then, implemented in a finite element program and applied to a lead‐free solder joint analysis. It is demonstrated that the implicit integration is very accurate if the multilinear kinematic hardening model of Ohno and Wang is employed, and that the consistent tangent modulus certainly affords quadratic convergence to the Newton–Raphson iteration in solving nodal force equilibrium equations. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical integration of a pressure‐dependent,non‐linear kinematic hardening constitutive model for large strain cyclic plasticity of metals

A Gurson‐based constitutive model is presented, which includes non‐linear mixed isotropic–kinematic hardening and creep, and allows the analysis of problems involving arbitrarily large plastic strains. This model was developed with the main objective of allowing, on the basis of a single set of material parameters, the numerical simulation of all the main features of cold metal forming processes, which usually imply severe loading–unloading cycles with very large plastic strains, difficult to be correctly reproduced numerically. A suitable integration scheme of the rate equations is described and implemented into a finite element code. The results obtained are compared with some reference experimental ones; an application of the model for the simulation of wire drawing processes is also presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Stress integration and mesh refinement for large deformation in geomechanics

This paper first discusses alternative stress integration schemes in numerical solutions to large‐ deformation problems in hardening materials. Three common numerical methods, i.e. the total‐Lagrangian (TL), the updated‐Lagrangian (UL) and the arbitrary Lagrangian–Eulerian (ALE) methods, are discussed. The UL and the ALE methods are further complicated with three different stress integration schemes. The objectivity of these schemes is discussed. The ALE method presented in this paper is based on the operator‐split technique where the analysis is carried out in two steps; an UL step followed by an Eulerian step. This paper also introduces a new method for mesh refinement in the ALE method. Using the known displacements at domain boundaries and material interfaces as prescribed displacements, the problem is re‐analysed by assuming linear elasticity and the deformed mesh resulting from such an analysis is then used as the new mesh in the second step of the ALE method. It is shown that this repeated elastic analysis is actually more efficient than mesh generation and it can be used for general cases regardless of problem dimension and problem topology. The relative performance of the TL, UL and ALE methods is investigated through the analyses of some classic geotechnical problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Assumed‐deformation gradient finite elements with nodal integration for nearly incompressible large deformation analysis

An assumed‐strain finite element technique for non‐linear finite deformation is presented. The weighted‐residual method enforces weakly the balance equation with the natural boundary condition and also the kinematic equation that links the elementwise and the assumed‐deformation gradient. Assumed gradient operators are derived via nodal integration from the kinematic‐weighted residual. A variety of finite element shapes fits the derived framework: four‐node tetrahedra, eight‐, 27‐, and 64‐node hexahedra are presented here. Since the assumed‐deformation gradients are expressed entirely in terms of the nodal displacements, the degrees of freedom are only the primitive variables (displacements at the nodes). The formulation allows for general anisotropic materials and no volumetric/deviatoric split is required. The consistent tangent operator is inexpensive and symmetric. Furthermore, the material update and the tangent moduli computation are carried out exactly as for classical displacement‐based models; the only deviation is the consistent use of the assumed‐deformation gradient in place of the displacement‐derived deformation gradient. Examples illustrate the performance with respect to the ability of the present technique to resist volumetric locking. A constraint count can partially explain the insensitivity of the resulting finite element models to locking in the incompressible limit. Copyright © 2008 John Wiley & Sons, Ltd.     

An energy–momentum conserving algorithm for non‐linear hypoelastic constitutive models

This paper presents an extension of the energy momentum conserving algorithm, usually developed for hyperelastic constitutive models, to the hypoelastic constitutive models. For such a material no potential can be defined, and thus the conservation of the energy is ensured only if the elastic work of the deformation can be restored by the scheme. We propose a new expression of internal forces at the element level which is shown to verify this property. We also demonstrate that the work of plastic deformation is positive and consistent with the material model. Finally several numerical applications are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

A large deformation solid‐shell concept based on reduced integration with hourglass stabilization

In this paper a new eight‐node (brick) solid‐shell finite element formulation based on the concept of reduced integration with hourglass stabilization is presented. The work focuses on static problems. The starting point of the derivation is the three‐field variational functional upon which meanwhile established 3D enhanced strain concepts are based. Important additional assumptions are made to transfer the approach into a powerful solid‐shell. First of all, a Taylor expansion of the first Piola–Kirchhoff stress tensor with respect to the normal through the centre of the element is carried out. In this way the stress becomes a linear function of the shell surface co‐ordinates whereas the dependence on the thickness co‐ordinate remains non‐linear. Secondly, the Jacobian matrix is replaced by its value in the centre of the element. These two assumptions lead to a computationally efficient shell element which requires only two Gauss points in the thickness direction (and one Gauss point in the plane of the shell element). Additionally three internal element degrees‐of‐freedom have to be determined to avoid thickness locking. One important advantage of the element is the fact that a fully three‐dimensional stress state can be modelled without any modification of the constitutive law. The formulation has only displacement degrees‐of‐freedom and the geometry in the thickness direction is correctly displayed. Copyright © 2006 John Wiley & Sons, Ltd.     

A rate‐consistent and intrinsically path‐dependent incremental method for plastic buckling problems

This paper presents an incremental predictor–corrector method which is able to handle the continuous spreading of elastic unloading and, therefore, is particularly well suited to solve plastic bucking problems. The method, which deals explicitly with rate variables and equations, is (i) rate consistent, because it leads to the ‘true’ tangent matrix, and (ii) intrinsically path‐dependent, because it enables an adequate identification and characterization of the onset of elastic unloading within each incremental step. Copyright © 2000 John Wiley & Sons, Ltd.     

Gradient and dilatational stabilizations for stress‐point integration in the element‐free Galerkin method

Stabilized stress‐point integration schemes based on gradient stabilization and dilatational stabilization methods are presented for linear elastostaticity problems in the framework of element‐free Galerkin (EFG) method. The instability in stress fields associated with the stress‐point integration is treated by the addition to the Galerkin weak form of stabilization terms which contain product of the gradient of the residual or the trace of the gradient of the residual; the latter is called dilatational stabilization. Numerical results show that the oscillations in the stress fields are successfully removed by the presented stabilization methods, and that the convergence and stability properties of direct stress‐point integration are greatly improved. These stabilization methods are particularly suitable for the solution of non‐linear continua with explicit time integration methods. Copyright © 2008 John Wiley & Sons, Ltd.     

An implicit unitless error and step‐size control method in integrating unified viscoplastic/creep ODE‐type constitutive equations

A series of numerical analyses are carried out to investigate the difficulties in numerical integration of unified viscoplastic/creep constitutive equations, which are normally represented as a system of ordinary differential equations (ODEs). The problems of numerically integrating the constitutive equations are identified and analysed. To overcome the stiffness problems, implicit methods are used for the numerical integration and a generic technique is introduced to calculate the Jacobian matrix. A normalization technique is introduced in the paper to convert the integration errors for each time increment to unitless errors. Thus, a single tolerance can be used to control the accuracy and step size in integrating a set of unified viscoplastic/creep constitutive equations. In addition, an implicit step‐size control method is proposed and used in the integrations. This method reduces the possibility of rejection of an integration increment due to poor accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

Finite elements in space and time for the analysis of generalised visco‐elastic materials

Numerical analysis of linear visco‐elastic materials requires robust and stable methods to integrate partial differential equations in both space and time. In this paper, symmetric space–time finite element operators are derived for the first time for elementary linear elastic spring and linear viscous dashpot. These can thereafter be assembled in parallel and in series to simulate an arbitrarily complex linear visco‐elastic behaviour. The flexibility of the proposed method allows the formulation of the behaviour, which closely reflects physical processes. An efficient algorithm is proposed to use the generated elementary matrices in a way that is comparable with finite difference schemes, in terms of both processor and memory costs. This unconditionally stable and convergent procedure is equally valid for space domains in which geometry or material properties evolve with time. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of different forms of strain energy functions in hyperelasticity‐based crystal plasticity models on texture evolution and mechanical response of face‐centered cubic crystals

Micro‐mechanical and macro‐mechanical behavior of face‐centered cubic (FCC) crystals is investigated by using different forms of strain energy functions in hyperelastic material models in crystal plasticity finite element framework. A quadratic strain energy function with anisotropic elastic constants, a polyconvex strain energy function with invariants associated with the cubic symmetry, and a strain energy function from an inter‐atomic potential are considered in hyperelastic material models to describe the elastic deformation of FCC crystals. In our numerical experiments, the trajectories of {111} poles in the pole figure and the accumulated plastic slips of FCC coppers under uniaxial tension and simple shear depend on the choice of strain energy functions when the slip resistance of the slip systems is high. The ability of strain energy functions in this study to represent elastic lattice distortions in crystals varies with the amount of elastic deformation and the shape of deformed lattice. However, numerical results show that the change of macroscopic mechanical behavior of FCC coppers is not significant for the choice of strain energy functions, compared with the change of crystallographic texture evolution. Copyright © 2014 John Wiley & Sons, Ltd.     

Frequency‐domain analysis of time‐integration methods for semidiscrete finite element equations—part II: Hyperbolic and parabolic–hyperbolic problems

Time‐integration methods for semidiscrete finite element equations of hyperbolic and parabolic– hyperbolic types are analysed in the frequency domain. The discrete‐time transfer functions of six popular methods are derived, and subsequently the forced response characteristics of single modes are studied in the frequency domain. Three characteristic numbers are derived which eliminate the parameter dependence of the frequency responses. Frequency responses and L₂‐norms of the phase and magnitude errors are calculated, and comparisons are given of the methods. As shown; the frequency‐domain analysis explains all time‐domain properties of the methods, and gives more insight into the subject. Copyright © 2001 John Wiley & Sons, Ltd.     

A fast and robust hybrid method for block‐structured mesh deformation with emphasis on FSI‐LES applications

The present work introduces an efficient technique for the deformation of block‐structured grids occurring in simulations of fluid–structure interaction (FSI) problems relying on large‐eddy simulation (LES). The proposed hybrid approach combines the advantages of the inverse distance weighting (IDW) interpolation with the simplicity and low computational effort of transfinite interpolation (TFI), while preserving the mesh quality in boundary layers. It is an improvement over the state‐of‐the‐art currently in use. To reach this objective, in a first step, three elementary mesh deformation methods (TFI, IDW, and radial basis functions) are investigated based on several test cases of different complexities analyzing not only their capabilities but also their computational costs. That not only allows to point out the advantages of each method but also demonstrates their drawbacks. Based on these specific properties of the different methods, a hybrid methodology is suggested that splits the entire grid deformation into two steps: first, the movement of the block‐boundaries of the block‐structured grid and second, the deformation of each block of the grid. Both steps rely on different methodologies, which allows to work out the most appropriate method for each step leading to a reasonable compromise between the grid quality achieved and the computational effort required. Finally, a hybrid IDW‐TFI methodology is suggested that best fits to the specific requirements of coupled FSI‐LES applications. This hybrid procedure is then applied to a real‐life FSI‐LES case. Copyright © 2016 John Wiley & Sons, Ltd.     

A practical large‐strain solid finite element for sheet forming

An alternative approach for developing practical large‐strain finite elements has been introduced and used to create a three‐dimensional solid element that exhibits no locking or hourglassing, but which is more easily and reliably derived and implemented than typical reduced‐integration schemes with hourglassing control. Typical large‐strain elements for forming applications rely on reduced integration to remove locking modes that occur with the coarse meshes that are necessary for practical use. This procedure introduces spurious zero‐energy deformation modes that lead to hourglassing, which in turn is controlled by complex implementations that involve lengthy derivations, knowledge of the material model, and/or undetermined parameters. Thus, for a new material or new computer program, implementation of such elements is a daunting task. Wang–Wagoner‐3‐dimensions (WW3D), a mixed, hexahedral, three‐dimensional solid element, was derived from the standard linear brick element by ignoring the strain components corresponding to locking modes while maintaining full integration (8 Gauss points). Thus, WW3D is easily implemented for any material law, with little chance of programming error, starting from programming for a readily available linear brick element. Surprisingly, this approach and resulting element perform similarly or better than standard solid elements in a series of numerical tests appearing in the literature. The element was also tested successfully for an applied sheet‐forming analysis problem. Many variations on the scheme are also possible for deriving special‐purpose elements. Copyright © 2005 John Wiley & Sons, Ltd.