On a new integration scheme for von‐Mises plasticity with linear hardening

Limiting the discussion to an associative von‐Mises plasticity model with linear kinematic and isotropic hardening, we compare the performance of the classical radial return map algorithm with a new integration scheme based on the computation of an integration factor. The numerical examples clearly show the improved accuracy of the new method. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

An integration procedure based on multi‐surface plasticity for plane stress plasticity: Theoretical formulation and application

An integration procedure designed to satisfy plane stress conditions for any constitutive law initially described in 3D and based on classical plasticity theory is presented herein. This method relies on multi‐surface plasticity, which allows associating in series various mechanisms. Three mechanisms have ultimately been used and added to the first one to satisfy the plane stress conditions. They are chosen to generate a plastic flow in the 3 out‐of‐plane directions, whose stresses must be canceled (σ₃₃,σ₁₃, and σ₂₃). The advantage of this method lies in its ease of use for every plastic constitutive law (in the general case of the non‐associated flow rule and with both nonlinear kinematic and isotropic hardening). Method implementation using a cutting plane algorithm is presented in its general framework and then illustrated by the example of a J2‐plasticity material model considering linear kinematic and isotropic hardening. The approach is compared with the same J2‐plasticity model that has been directly derived from a projection of its equations onto the plane stress subspace. The performance of the multi‐surface plasticity method is shown through the comparison of iso‐error and iso‐step contours in both formulations, and lastly with a case study considering a hollow plate subjected to tension. Copyright © 2014 John Wiley & Sons, Ltd.     

Limit analysis of plates using the EFG method and second‐order cone programming

The meshless element‐free Galerkin (EFG) method is extended to allow computation of the limit load of plates. A kinematic formulation that involves approximating the displacement field using the moving least‐squares technique is developed. Only one displacement variable is required for each EFG node, ensuring that the total number of variables in the resulting optimization problem is kept to a minimum, with far fewer variables being required compared with finite element formulations using compatible elements. A stabilized conforming nodal integration scheme is extended to plastic plate bending problems. The evaluation of integrals at nodal points using curvature smoothing stabilization both keeps the size of the optimization problem small and also results in stable and accurate solutions. Difficulties imposing essential boundary conditions are overcome by enforcing displacements at the nodes directly. The formulation can be expressed as the problem of minimizing a sum of Euclidean norms subject to a set of equality constraints. This non‐smooth minimization problem can be transformed into a form suitable for solution using second‐order cone programming. The procedure is applied to several benchmark beam and plate problems and is found in practice to generate good upper‐bound solutions for benchmark problems. Copyright © 2009 John Wiley & Sons, Ltd.     

A novel ‘optimal’ exponential‐based integration algorithm for von‐Mises plasticity with linear hardening: Theoretical analysis on yield consistency,accuracy, convergence and numerical investigations

In this communication we propose a new exponential‐based integration algorithm for associative von‐Mises plasticity with linear isotropic and kinematic hardening, which follows the ones presented by the authors in previous papers. In the first part of the work we develop a theoretical analysis on the numerical properties of the developed exponential‐based schemes and, in particular, we address the yield consistency, exactness under proportional loading, accuracy and stability of the methods. In the second part of the contribution, we show a detailed numerical comparison between the new exponential‐based method and two classical radial return map methods, based on backward Euler and midpoint integration rules, respectively. The developed tests include pointwise stress–strain loading histories, iso‐error maps and global boundary value problems. The theoretical and numerical results reveal the optimal properties of the proposed scheme. Copyright © 2006 John Wiley & Sons, Ltd.     

Implicit time integration for the material point method: Quantitative and algorithmic comparisons with the finite element method

An implicit integration strategy was developed and implemented for use with the material point method (MPM). An incremental‐iterative solution strategy was developed around Newton's method to solve the equations of motion with Newmark integration to update the kinematic variables. Test problems directly compared the implicit MPM solutions with those obtained using an explicit MPM code and implicit finite element (FE) code. Results demonstrated very good agreement with FE predictions and also illustrated several advantages in comparison to calculations using the explicit MPM code. In particular, the accuracy of the implicit solution was superior to the explicit MPM when compared to validated FE solutions, and by definition the implicit time integration is unconditionally stable. Similarities between the assembly of the implicit MPM equations and those of the FE method were identified and should allow further algorithmic improvement. Copyright © 2003 John Wiley & Sons, Ltd.     

Micropolar hyper‐elastoplasticity: constitutive model,consistent linearization,and simulation of 3D scale effects

A computational model for micropolar hyperelastic‐based finite elastoplasticity that incorporates isotropic hardening is developed. The basic concepts of the non‐linear micropolar kinematic framework are reviewed, and a thermodynamically consistent constitutive model that features Neo‐Hooke‐type elasticity and generalized von Mises plasticity is described. The integration of the constitutive initial value problem is carried out by means of an elastic‐predictor/plastic‐corrector algorithm, which retains plastic incompressibility. The solution procedure is developed carefully and described in detail. The consistent material tangent is derived. The micropolar constitutive model is implemented in an implicit finite element framework. The numerical example of a notched cylindrical bar subjected to large axial displacements and large twist angles is presented. The results of the finite element simulations demonstrate (i) that the methodology is capable of capturing the size effect in three‐dimensional elastoplastic solids in the finite strain regime, (ii) that the formulation possesses a regularizing effect in the presence of strain localization, and (iii) that asymptotically quadratic convergence rates of the Newton–Raphson procedure are achieved. Throughout this paper, effort is made to present the developments as a direct extension of standard finite deformation computational plasticity. Copyright © 2012 John Wiley & Sons, Ltd.     

An implicit algorithm using explicit correctors for the kinematic hardening model with multiple back stresses

This paper presents an efficient mathematical algorithm for a class of non‐linear kinematic hardening models with multiple back stresses, as an extension of the implicit integration algorithm for a single back stress hardening model. Explicit formulations for general three‐dimensional stress states as well as plane stress and plane strain are given. The new formulation is implemented in a general‐purpose finite element code, ABAQUS, and is verified by comparison with the existing formulation for the single back‐stress constitutive model. Comparison is also made with the experimental results obtained from a plate containing a circular hole subjected to cyclic loading, demonstrating the validity of new method. Copyright © 2001 John Wiley & Sons, Ltd.     

Smooth finite element methods: Convergence,accuracy and properties

A stabilized conforming nodal integration finite element method based on strain smoothing stabilization is presented. The integration of the stiffness matrix is performed on the boundaries of the finite elements. A rigorous variational framework based on the Hu–Washizu assumed strain variational form is developed. We prove that solutions yielded by the proposed method are in a space bounded by the standard, finite element solution (infinite number of subcells) and a quasi‐equilibrium finite element solution (a single subcell). We show elsewhere the equivalence of the one‐subcell element with a quasi‐equilibrium finite element, leading to a global a posteriori error estimate. We apply the method to compressible and incompressible linear elasticity problems. The method can always achieve higher accuracy and convergence rates than the standard finite element method, especially in the presence of incompressibility, singularities or distorted meshes, for a slightly smaller computational cost. It is shown numerically that the one‐cell smoothed four‐noded quadrilateral finite element has a convergence rate of 2.0 in the energy norm for problems with smooth solutions, which is remarkable. For problems with rough solutions, this element always converges faster than the standard finite element and is free of volumetric locking without any modification of integration scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

Integration schemes for von‐Mises plasticity models based on exponential maps: numerical investigations and theoretical considerations

We consider three different exponential map algorithms for associative von‐Mises plasticity with linear isotropic and kinematic hardening. The first scheme is based on a different formulation of the time continuous plasticity model, which automatically grants the yield consistency of the method in the numerical solution. The second one is the quadratically accurate but non‐yield consistent method already proposed in Auricchio and Beirão da Veiga (Int. J. Numer. Meth. Engng 2003; 56 : 1375–1396). The third method is an improved version of the second one, in which the yield consistency condition is enforced a posteriori. We also compare the performance of the three methods with the classical radial return map algorithm. We develop extensive numerical tests which clearly show the main advantages and disadvantages of the three methods. Copyright © 2005 John Wiley & Sons, Ltd.     

Accelerating cyclic plasticity simulations using an adaptive wavelet transformation based multitime scaling method

Cyclic finite element simulations of complex materials, for example, polycrystalline metals, are widely used to study fatigue failure due to plasticity and damage. Typically, this requires the simulation of a large number of cycles to failure for accurate determination of evolving deformation variables. Modeling cyclic deformation using conventional methods of time integration in semidiscretization techniques can however be computationally challenging. Single time scale integration methods typically follow the high frequency characteristics and discretize each cycle into a number of time steps over which integration is performed. To overcome this computational challenge, the wavelet transformation‐based multitime scale (WATMUS) method proposed in an earlier work by the authors is advanced and validated in this paper to perform accelerated finite element simulations of materials undergoing rate‐dependent plasticity for large number of cycles. Specifically, the WATMUS algorithm is integrated with crystal plasticity finite element method to perform accelerated simulations of polycrystalline alloys. The WATMUS methodology is also endowed with adaptive capabilities to optimally construct the wavelet basis functions and determine coarse‐scale cycle steps. Accuracy and efficiency of the WATMUS methodology is conclusively demonstrated by comparing the results with cyclic single‐time scale crystal plasticity finite element simulations performed on image‐based microstructure of titanium alloys. Copyright © 2012 John Wiley & Sons, Ltd.     

Remarks on the interpretation of current non‐linear finite element analyses as differential–algebraic equations

For the numerical solution of materially non‐linear problems like in computational plasticity or viscoplasticity the finite element discretization in space is usually coupled with point‐wise defined evolution equations characterizing the material behaviour. The interpretation of such systems as differential–algebraic equations (DAE) allows modern‐day integration algorithms from Numerical Mathematics to be efficiently applied. Especially, the application of diagonally implicit Runge–Kutta methods (DIRK) together with a Multilevel‐Newton method preserves the algorithmic structure of current finite element implementations which are based on the principle of virtual displacements and on backward Euler schemes for the local time integration. Moreover, the notion of the consistent tangent operator becomes more obvious in this context. The quadratical order of convergence of the Multilevel‐Newton algorithm is usually validated by numerical studies. However, an analytical proof of this second order convergence has already been given by authors in the field of non‐linear electrical networks. We show that this proof can be applied in the current context based on the DAE interpretation mentioned above. We finally compare the proposed procedure to several well‐known stress algorithms and show that the inclusion of a step‐size control based on local error estimations merely requires a small extra time‐investment. Copyright © 2001 John Wiley & Sons, Ltd.     

A convergent adaptive finite element method for the primal problem of elastoplasticity

The boundary value problem representing one time step of the primal formulation of elastoplasticity with positive hardening leads to a variational inequality of the second kind with some nondifferentiable functional. This paper establishes an adaptive finite element algorithm for the solution of this variational inequality that yields the energy reduction and, up to higher order terms, the R‐linear convergence of the stresses with respect to the number of loops. Applications include several plasticity models: linear isotropic‐kinematic hardening, linear kinematic hardening, and multisurface plasticity as model for nonlinear hardening laws. For perfect plasticity, the adaptive algorithm yields strong convergence of the stresses. Numerical examples confirm an improved linear convergence rate and study the performance of the algorithm in comparison with the more frequently applied maximum refinement rule. Copyright © 2006 John Wiley & Sons, Ltd.     

Stress update algorithm for the combined viscoplastic and plastic behaviours of single‐crystal superalloys

A crystallographic constitutive model is developed, which accounts for both rate‐sensitive and rate‐insensitive flow. Single‐crystal plasticity and viscoplasticity are the limiting cases of the model, so that it properly reflects the material response over a wide temperature range. A non‐linear dynamic recovery is included to properly describe ratchetting. We provide a robust integration scheme based on generalization of the return‐mapping algorithm and of the procedure for active set search. The implicit integration and consistent tangent are implemented through the UMAT subroutine in the ABAQUS finite element program. The capability of the model to account for both high and low strain rates is demonstrated in numerical examples. Finally, the stability of integration scheme and quadratic convergence of the global Newton–Raphson equilibrium iterations are demonstrated on the example of a notched bar under tension. Copyright © 2011 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.     

A finite element model for thermomechanical analysis of sheet metal forming

A thermal model based on explicit time integration is developed and implemented into the explicit finite element code DYNA3D to model simultaneous forming and quenching of thin‐walled structures. A staggered approach is used for coupling the thermal and mechanical analysis, wherein each analysis is performed with different time step sizes. The implementation includes a thermal shell element with linear temperature approximation in the plane and quadratic in the thickness direction, and contact heat transfer. The material behaviour is described by a temperature‐dependent elastic–plastic model with a non‐linear isotropic hardening law. Transformation plasticity is included in the model. Examples are presented to validate and evaluate the proposed model. The model is evaluated by comparison with a one‐sided forming and quenching experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

Implicit plane stress algorithm for multilayer J2‐plasticity using the Prager–Ziegler translation rule

A totally implicit algorithm for plane stress multilayer plasticity is presented. The algorithm presents the plane stress version of the 3D/plane strain model for multilayer plasticity presented in a previous work. As in the 3D/plane strain case, the model is consistent with the principle of maximum dissipation and it may be considered as an extension of classical J₂‐plasticity for anisotropic non‐linear kinematic behaviour preserving Masing's rules. In order to obtain the asymptotic second‐order convergence of the Newton algorithm, the consistent tangent moduli are also given. Copyright © 2003 John Wiley & Sons, Ltd.     

基于参数二次规划与精细积分方法的动力弹塑性问题分析

给出了将参数二次规划方法与精细积分方法相结合进行结构弹塑性动力响应分析的一条新途径。基于参变量变分原理与有限元参数二次规划方法建立了动力弹塑性问题的求解方程，方法对于关联与非关联问题的求解在算法上是完全一致的。对于动力非线性方程求解则进一步采用了被线性问题分析所广泛采用的精细积分方法，推导了方法在动力弹塑性问题求解上的算法列式。所给出的数值算例在验证本文理论与算法的同时，进一步证实了精细积分方法在动力学分析中所具有的各种良好性态。     

A robust and efficient substepping scheme for the explicit numerical integration of a rate‐dependent crystal plasticity model

This paper describes the development of efficient and robust numerical integration schemes for rate‐dependent crystal plasticity models. A forward Euler integration algorithm is first formulated. An integration algorithm based on the modified Euler method with an adaptive substepping scheme is then proposed, where the substepping is mainly controlled by the local error of the stress predictions within the time step. Both integration algorithms are implemented in a stand‐alone code with the Taylor aggregate assumption and in an explicit finite element code. The robustness, accuracy and efficiency of the substepping scheme are extensively evaluated for large time steps, extremely low strain‐rate sensitivity, high deformation rates and strain‐path changes using the stand‐alone code. The results show that the substepping scheme is robust and in some cases one order of magnitude faster than the forward Euler algorithm. The use of mass scaling to reduce computation time in crystal plasticity finite element simulations for quasi‐static problems is also discussed. Finally, simulation of Taylor bar impact test is carried out to show the applicability and robustness of the proposed integration algorithm for the modelling of dynamic problems with contact. Copyright © 2014 John Wiley & Sons, Ltd.