Numerical integration and FEM-implementation of a viscoplastic Chaboche-model with static recovery

This paper is concerned with the implementation of a viscoplastic material model of the Chaboche type in the framework of the finite element method (FEM). The equations of the used constitutive law, that incorporates isotropic hardening, back stress evolution with static recovery terms and drag stress evolution, are introduced. A representation of their numerical integration using the implicit backward Euler method under the assumption of small deformations and an isothermal formulation follows. The use of the backward Euler method leads to a nonlinear algebraic system of three equations, which is solved by a combination of the Pegasus method and a fixed-point iteration. After considering the accuracy of the presented integration algorithm in form of iso-error maps, the derivation of the consistent viscoplastic tangent operator is shown. The integration scheme and the calculation of the consistent viscoplastic tangent operator are implemented in the commercial finite element code ABAQUS, using the possibility of the user-defined material subroutine (UMAT). Finally a numerical example in form of a notched bar under tension is presented.     

Degenerated shell element with composite implicit time integration scheme for geometric nonlinear analysis

A degenerated shell element with composite implicit time integration scheme is developed in the present paper to solve the geometric nonlinear large deformation and dynamics problems of shell structures. The degenerated shell element is established based on the eight‐node solid element, where the nodal forces, mass matrices, and stiffness matrices are firstly obtained upon virtual velocity principle and then translated to the shell element. The strain field is modified based on the mixed interpolation of tensorial components method to eliminate the shear locking, and the constitutive relation is modified to satisfy the shell assumptions. A simple and practical computational method for nonlinear dynamic response is developed by embedding the composite implicit time integration scheme into the degenerated shell element, where the composite scheme combines the trapezoidal rule with the three‐point backward Euler method. The developed approach can not only keep the momentum and energy conservation and decay the high frequency modes but also lead to a symmetrical stiffness matrix. Numerical results show that the developed degenerated shell element with the composite implicit time integration scheme is capable of solving the geometric nonlinear large deformation and dynamics problems of the shell structures with momentum and energy conservation and/or decay. Copyright © 2015 John Wiley & Sons, Ltd.     

A computational model for the finite element analysis of thermoplasticity coupled with ductile damage at finite strains

An updated Lagrangian implicit FEM model for the analysis of large thermo‐mechanically coupled hyperelastic‐viscoplastic deformations of isotropic porous materials is considered. An appropriate framework for constitutive modelling is introduced that includes a stress‐free thermally expanded configuration and a plastically deformed unstressed damaged configuration. A two‐level iterative scheme is employed at each time increment to solve the field equations governing the conservation of momentum (mechanical step) and the conservation of energy (thermal step) for the coupled thermo‐mechanical problem. Exact linearizations for the calculation of the tangent stiffness are performed in each of these solution steps. A fully implicit, thermo‐mechanically coupled and incrementally objective Euler‐backward radial return based map is developed for the time integration of the constitutive equations. The present model is used to analyse a number of benchmark examples including metal forming processes wherein temperature and the accumulated damage play an important role in influencing the deformation mechanism and the nature of the deformed workpiece. Copyright © 1999 John Wiley & Sons, Ltd.     

Dynamic elasto-plastic response of shells in an acoustic medium—EPSA code

A nonlinear, large deflection, elasto-plastic finite element code (EPSA) has been developed for the analysis of shells in an acoustic medium subjected to dynamic loadings. The nonlinear equations of shells are discretized with the aid of a finite difference/finite element method based upon the principle of virtual work. The resulting system of equations contains the nodal displacements as the generalized co-ordinates of the problem. The integration in time of the equations of motion is done explicitly via a central difference scheme. Shell strain-displacement relations are established by a two-dimensional finite difference scheme. The shell constitutive equations are formulated in terms of the shell stress resultants and the shell strains and curvatures. The fluid-structure interaction is accounted for by means of the doubly asymptotic approximation (DAA) expressed in terms of orthogonal fluid expansion functions. The analytically produced results satisfactorily reproduce available experimental data for dynamically loaded shells.     

Fully implicit integration and consistent tangent modulus in elasto-plasticity

This paper deals with the numerical integration of a class of rate-independent elasto-plastic models. The backward Euler scheme is used to integrate the rate constitutive relations. The non-linear equations obtained are solved by the Newton method. The consistent tangent modulus is obtained by exact linearization of the algorithm. In the case of J₂ elasto-plasticity with non-linear isotropic hardening and non-linear kinematic hardening (Chaboche-Marquis model), explicit formulas are derived, without any approximations.     

An implicit time integration scheme for inelastic constitutive equations with internal state variables

In a broad class of inelastic constitutive models for the deformation of metals the inelastic strain rates are functions of the current state of stress and internal state variables only. All known models are in some regions of application mathematically stiff and therefore difficult to integrate. The unconditionally stable implicit Euler rule is used for integration. It leads to a system of highly nonlinear algebraic equations which have to be solved by an iterative process. The general Newton-Raphson method, which converges under very broad conditions, requires repeated solution of the finite element system and is infeasible for large inelastic problems. But for the inelastic strains and internal state variables the Jacobian can be computed analytically and therefore the NRI can be used. For the stresses the Jacobian cannot be computed analytically and therefore the accelerated Jacobi iteration is used. A new method for computing the relaxation parameter is introduced which increases the rate of convergence significantly. The new algorithm is applied on Hart's model. A comparison with prior computations using an approximation is made.     

A return mapping algorithm for isotropic elastoplasticity

A return mapping algorithm is presented for the numerical time integration of the constitutive equations for elastoplasticity with isotropic yield surfaces, constructed from all three invariants of the stress tensor. Based on the first-order backward Euler difference formula (BDF), the governing equations for the stresses are solved in the space of the invariants and the discretized persistence parameter. The stresses are recovered afterwards. The solution concept is applied to a pressure-independent yield function, expressed in terms of the second and third invariant of the stress tensor. The numerical performance of the method is demonstrated with two examples.     

循环稳定材料的棘轮行为:II.隐式应力积分算法和有限元实现

针对第一部分发展的、能够合理描述循环稳定材料棘轮行为的粘塑性本构模型,详细讨论该模型的数值计算方法和有限元实现。在径向回退(RadialReturn)和向后欧拉积分方法的基础上,结合连续迭代(SuccessiveSubstitution)方法,推导并建立了针对循环粘塑性本构模型的、新的隐式应力积分算法。为了本构模型在大型有限元分析程序(如ABAQUS等)中的实现,针对有限元的整体节点迭代计算,推导和确立了一个新的、考虑率相关塑性的一致切线刚度矩阵(ConsistentTangentModulus)表达式。通过对一些算例的有限元分析,讨论了建立的隐式应力积分算法的优越性,同时对特定构件的棘轮行为进行了数值模拟,进而检验了有限元实现的合理性和必要性。     

A semi‐implicit integration scheme for a combined viscoelastic‐damage model of plastic bonded explosives

This paper presents a new implementation of a constitutive model commonly used to represent plastic bonded explosives in finite element simulations of thermomechanical response. The constitutive model, viscoSCRAM, combines linear viscoelasticity with isotropic damage evolution. The original implementation was focused on short duration transient events; thus, an explicit update scheme was used. For longer duration simulations that employ significantly larger time step sizes, the explicit update scheme is inadequate. This work presents a new semi‐implicit update scheme suitable for simulations using relatively large time steps. The algorithm solves a nonlinear system of equations to ensure that the stress, damaged state, and internal stresses are in agreement with implicit update equations at the end of each increment. The crack growth is advanced in time using a sub‐incremental explicit scheme; thus, the entire implementation is semi‐implicit. The theory is briefly discussed along with previous explicit integration schemes. The new integration algorithm and its implementation into the finite element code, Abaqus, are detailed. Finally, the new and old algorithms are compared via simulations of uniaxial compression and beam bending. The semi‐implicit scheme has been demonstrated to provide higher accuracy for a given allocated computational time for the quasistatic cases considered here. Published 2014. This article is a US Government work and is in the public domain in the USA.     

Explicit numerical integration algorithm for a class of non-linear kinematic hardening model

 An explicit updating algorithm has been developed for the Armstrong–Frederick family of non-linear kinematic hardening model, based on the trapezoidal and the backward Euler integration method. The algorithm provides a computationally efficient method for implementing the non-linear kinematic hardening model in finite element codes. It is shown that the trapezoidal method performs better with the original Armstrong–Frederick rule, while the backward Euler rule provides an improved accuracy to the modified multiple back-stress model that incorporates a weight function for dynamic recovery. Numerical examples are presented to illustrate the performance of the algorithm developed, and a comparison with the experimental observation shows that the modified constitutive model indeed provides a more accurate prediction to the long term mean stress relaxation.     

A comprehensive implicit substepping integration scheme for multisurface plasticity

A complex elastoplastic model requires a robust integration procedure of the evolution equations. The performance of the finite element solution is directly affected by the convergence characteristics of the state-update procedure. Thereby, this study proposes a comprehensive numerical integration scheme to deal with generic multisurface plasticity models. This algorithm is based on the backward Euler method aiming at accuracy and stability, and on the Newton–Raphson method to solve the unconstrained optimization problem. In this scenario, a line search strategy is adopted to improve the convergence characteristics of the algorithm. The golden section method, an exact line search, is considered. Also, a substepping scheme is implemented to provide additional robustness to the state-update procedure. Therefore, this work contributes to computational plasticity presenting an adaptive substep size scheme and a consistent tangent modulus according to the substepping technique. Finally, some numerical problems are evaluated using the proposed algorithm. Single-surface and novel multisurface plasticity models are employed in these analyses. The results testify how the line search and substepping strategies can improve the robustness of the nonlinear analysis.     

Numerical simulation of martensitic phase transitions in shape memory alloys using an improved integration algorithm

The numerical simulation of structures made of shape memory materials is of increasing interest in different fields. Among others, the computation of pipe connectors or medical devices like endoscopic instruments and stents is a challenge. In such practical applications the pseudoelastic effect as well as the one‐way and two‐way shape memory effects are utilized. These material properties are caused by martensitic phase transitions between austenite and martensite. In the present contribution, a recently proposed constitutive theory is numerically treated in the context of the finite element method. This constitutive theory is formulated in the framework of continuum thermomechanics for geometrically linear problems and is able to represent the occurring martensitic phase transitions in shape memory alloys. For the numerical integration of the evolution equations, the backward Euler method is applied. In spite of the complexity of the constitutive theory, it is shown that an improved integration procedure can be formulated, which merely involves the solution of three non‐linear equations for three scalar‐valued unknown variables. Numerical examples show the capability of the proposed model and the improved integration algorithm. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical implementation of constitutive integration for rate-independent elastoplasticity

In this paper, constitutive integration for rate-independent, small deformation elastoplasticity is studied. Smooth yield surfaces and work/strain hardening are assumed. Both associative or non-associative flow rules are considered. An Euler backward algorithm is applied for constitutive integration. Tangent moduli that are consistent with the Euler backward algorithm, i.e. a so-called consistent tangent operator, are derived. Emphasis is placed on numerical implementation of the Eular backward algorithm into finite element codes using such a consistent tangent operator. In particular, a commercial code ANSYS is considered. Numerical examples, including materials sensitive and insensitive to hydrostatic stress, are used for the verification of the implementation. A comparison of the algorithmic performance to an explicit Euler forward algorithm is given and the superiority of the Euler backward algorithm is demonstrated.The work described in the present paper has been sponsored by The Research Council of Norway, The North Calotte Education and Research Council, Statoil and Norsk Hydro.     

Efficient time stepping for the multiplicative Maxwell fluid including the Mooney‐Rivlin hyperelasticity

A popular version of the finite‐strain Maxwell fluid is considered, which is based on the multiplicative decomposition of the deformation gradient tensor. The model combines Newtonian viscosity with hyperelasticity of the Mooney‐Rivlin type; it is a special case of the viscoplasticity model proposed by Simo and Miehe in 1992. A simple, efficient, and robust implicit time‐stepping procedure is suggested. Lagrangian and Eulerian versions of the algorithm are available, with equivalent properties. The numerical scheme is iteration free, unconditionally stable, and first order accurate. It exactly preserves the inelastic incompressibility, symmetry, and positive definiteness of the internal variables and w‐invariance. The accuracy of the stress computations is tested using a series of numerical simulations involving a nonproportional loading and large strain increments. In terms of accuracy, the proposed algorithm is equivalent to the modified Euler backward method with exact inelastic incompressibility; the proposed method is also equivalent to the classical integration method based on exponential mapping. Since the new method is iteration free, it is more robust and computationally efficient. The algorithm is implemented into MSC.MARC, and a series of initial boundary value problems is solved to demonstrate the usability of the numerical procedures.     

Thermomechanical finite element simulations of selective electron beam melting processes: performance considerations

The present contribution is concerned with the macroscopic modelling of the selective electron beam melting process by using the finite element method. The modelling and simulation of the selective electron beam melting process involves various challenges: complex material behaviour, phase changes, thermomechanical coupling, high temperature gradients, different time and length scales etc. The present contribution focuses on performance considerations of solution approaches for thermomechanically coupled problems, i.e. the monolithic and the adiabatic split approach. The material model is restricted to nonlinear thermoelasticity with temperature-dependent material parameters. As a numerical example a straight scanning path is simulated, the predicted temperatures and stresses are analysed and the performance of the two algorithms is compared. The adiabatic split approach turned out to be much more efficient for linear thermomechanical problems, i.e. the solution time is three times less than with the monolithic approach. For nonlinear problems, stability issues necessitated the use of the Euler backward integration scheme, and therefore, the adiabatic split approach required small time steps for reasonable accuracy. Thus, for nonlinear problems and in combination with the Euler backward integration scheme, the monolithic solver turned out to be more efficient.     

完全隐式返回映射算法对岩土地基问题的求解

针对岩土工程中的复杂力学问题,在弹塑性力学理论框架和非线性有限元理论基础上,采用非关联等向硬化Drucker-Prager模型的完全隐式积分算法—返回映射算法(Return Mapping Algorithm)编制了有限元求解程序。该算法可以避免预测应力漂移屈服面的现象,对准静态变形条件下的本构方程可以获得准确解,在迭代中使用Newton-Raphson法获得近似平方的收敛速率,具有较高的精确性和稳定性。对岩土工程中的地基问题进行求解,计算得出位移、应力等结果,模拟了塑性区随载荷步增加的演化过程,对地基极限承载力进行了解析解和数值解的对比。结果表明了算法的优越性、程序的正确性和实用性。     

Computational modeling of cardiac electrophysiology: A novel finite element approach

The key objective of this work is the design of an unconditionally stable, robust, efficient, modular, and easily expandable finite element‐based simulation tool for cardiac electrophysiology. In contrast to existing formulations, we propose a global–local split of the system of equations in which the global variable is the fast action potential that is introduced as a nodal degree of freedom, whereas the local variable is the slow recovery variable introduced as an internal variable on the integration point level. Cell‐specific excitation characteristics are thus strictly local and only affect the constitutive level. We illustrate the modular character of the model in terms of the FitzHugh–Nagumo model for oscillatory pacemaker cells and the Aliev–Panfilov model for non‐oscillatory ventricular muscle cells. We apply an implicit Euler backward finite difference scheme for the temporal discretization and a finite element scheme for the spatial discretization. The resulting non‐linear system of equations is solved with an incremental iterative Newton–Raphson solution procedure. Since this framework only introduces one single scalar‐valued variable on the node level, it is extremely efficient, remarkably stable, and highly robust. The features of the general framework will be demonstrated by selected benchmark problems for cardiac physiology and a two‐dimensional patient‐specific cardiac excitation problem. Copyright © 2009 John Wiley & Sons, Ltd.     

Stabilized finite element methods for vertically averaged multiphase flow for carbon sequestration

A computationally efficient numerical model that describes carbon sequestration in deep saline aquifers is presented. The model is based on the multiphase flow and vertically averaged mass balance equations, requiring the solution of two partial differential equations – a pressure equation and a saturation equation. The saturation equation is a nonlinear advective equation for which the application of Galerkin finite element method (FEM) can lead to non‐physical oscillations in the solution. In this article, we extend three stabilized FEM formulations, which were developed for uncoupled systems, to the governing nonlinear coupled PDEs. The methods developed are based on the streamline upwind, the streamline upwind/Petrov–Galerkin and the least squares FEM. Two sequential solution schemes are developed: a single step and a predictor–corrector. The range of Courant numbers yielding smooth and oscillation‐free solutions is investigated for each method. The useful range of Courant numbers found depends upon both the sequential scheme (single step vs predictor–corrector) and also the time integration method used (forward Euler, backward Euler or Crank–Nicolson). For complex problems such as when two plumes meet, only the SU stabilization with an amplified stabilization parameter gives satisfactory results when large time steps are used. Copyright © 2016 John Wiley & Sons, Ltd.     

Adaptive backward Euler time stepping with truncation error control for numerical modelling of unsaturated fluid flow

An automatic time stepping scheme with embedded error control is developed and applied to the moisture‐based Richards equation. The algorithm is based on the first‐order backward Euler scheme, and uses a numerical estimate of the local truncation error and an efficient time step selector to control the temporal accuracy of the integration. Local extrapolation, equivalent to the use of an unconditionally stable Thomas–Gladwell algorithm, achieves second‐order temporal accuracy at minimal additional costs. The time stepping algorithm also provides accurate initial estimates for the iterative non‐linear solver. Numerical tests confirm the ability of the scheme to automatically optimize the time step size to match a user prescribed temporal error tolerance. An important merit of the proposed method is its conceptual and computational simplicity. It can be directly incorporated into existing or new software based on the backward Euler scheme (currently prevalent in subsurface hydrologic modelling), and markedly improves their performance compared with simple fixed or heuristic time step selection. The generality of the approach also makes possible its use for solving PDEs in other engineering applications, where strong non‐linearity, stability or implementation considerations favour a simple and robust low‐order method, or where there is a legacy of backward Euler codes in current use. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical implementation and analysis of a class of metal plasticity models coupled with ductile damage

This paper deals with a class of rate-independent metal plasticity models which exhibit non-linear isotropic hardening, non-linear kinematic hardening (Chaboche-Marquis model) and ductile damage (Lemaitre-Chaboche model). The backward Euler scheme is used to integrate the rate constitutive relations. The non-linear equations obtained are solved by the Newton method. The consistent tangent operator is obtained by exact linearization of the algorithm. Despite the complexity of the constitutive equations, closed-form expressions are derived, without any approximations. Analytical, numerical and experimental results are presented and discussed.