On a finite-strain viscoplastic law coupled with anisotropic damage: theoretical formulations and numerical applications

Based on a dissipation inequality at finite strains and the effective stress concept, a Chaboche-type infinitesimal viscoplastic theory is extended to finite-strain cases coupled with anisotropic damage. The anisotropic damage is described by a rank-two symmetric tensor. The constitutive law is formulated in the corotational material coordinate system. Thus, the evolution equations of all internal variables can be expressed in terms of their material time derivatives. The numerical algorithm for implementing the material model in a finite element programme is also formulated, and several numerical examples are shown. Comparing the numerical simulations with experimental observations indicates that the present material model can describe well the primary, secondary and tertiary creep. It can also predict the anisotropic damage modes observed in experiments correctly.     

On internal dissipation inequalities and finite strain inelastic constitutive laws: Theoretical and numerical comparisons

Internal dissipation always occurs in irreversible inelastic deformation processes of materials. The internal dissipation inequalities (specific mathematical forms of the second law of thermodynamics) determine the evolution direction of inelastic processes. Based on different internal dissipation inequalities several finite strain inelastic constitutive laws have been formulated for instance by Simo [Simo, J.C., 1992. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering 99, 61–112]; Simo and Miehe [Simo, J.C., Miehe, C., 1992. Associative coupled thermoplasticity at finite strains: formulation, numerical analysis and implementation. Computer Methods in Applied Mechanics and Engineering 98, 41–104]; Lion [Lion, A., 1997. A physically based method to represent the thermo-mechanical behavior of elastomers. Acta Mechanica 123, 1–25]; Reese and Govindjee [Reese, S., Govindjee, S., 1998. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482]; Lin and Schomburg [Lin, R.C., Schomburg, U., 2003. A finite elastic–viscoelastic–elastoplastic material law with damage: theoretical and numerical aspects. Computer Methods in Applied Mechanics and Engineering 192, 1591–1627]; Lin and Brocks [Lin, R.C., Brocks, W., 2004. On a finite strain viscoplastic theory based on a new internal dissipation inequality. International Journal of Plasticity 20, 1281–1311]; and Lin and Brocks [Lin, R.C., Brocks, W., 2005. An extended Chaboche’s viscoplastic law at finite strains: theoretical and numerical aspects. Journal of Materials Science and Technology 21, 145–147]. These constitutive laws are consistent with the second law of thermodynamics. As the internal dissipation inequalities are described in different configurations or coordinate systems, the related constitutive laws are also formulated in the corresponding configurations or coordinate systems. Mathematically, these constitutive laws have very different formulations. Now, a question is whether the constitutive laws provide identical constitutive responses for the same inelastic constitutive problems. In the present work, four types of finite strain viscoelastic and endochronically plastic laws as well as three types of J₂-plasticity laws are formulated based on four types of dissipation inequalities. Then, they are numerically compared for several problems of homogeneous or complex finite deformations. It is demonstrated that for the same inelastic constitutive problem the stress responses are identical for deformation processes without rotations. In the simple shear deformation process with large rotation, the presented viscoelastic and endochronically plastic laws also show almost identical stress responses up to a shear strain of about 100%. The three laws of J₂-plasticity also produce the same shear stresses up to a shear strain of 100%, while different normal stresses are generated even at infinitesimal shear strains. The three J₂-plasticity laws are also compared at three complex finite deformation processes: billet upsetting, cylinder necking and channel forming. For the first two deformation processes similar constitutive responses are obtained, whereas for the third deformation process (with large global rotations) significant differences of constitutive responses can be observed.     

Original ArticleKinematic hardening rules in finite plasticity Part I: A constitutive approach

Plasticity laws exhibiting non-linear kinematichardening are considered within the framework of infinitesimal deformations. The evolution equations governing the response of kinematic hardening are derived as sufficient conditions in order for the intrinsic dissipation inequality to be satisfied in every process. With a view to the extension to finite deformations, two basic possibilities are proposed. In every case, an isotropic elasticity law with respect to the so-called plastic intermediate configuration is assumed to hold. The theory applicable to finite deformations is based on the concept of so-called dual variables and associated time derivatives. Thus, the main difference between the present work and other contributions in this area is the choice of the variables used to formulate the theory. In fact, using dual variables, hardening rules are derived as sufficient conditions for the intrinsic dissipation inequality to be satisfied in every process. This is quite analogous to the case of infinitesimal deformation, but now the hardening rules take a very specific form which is explained in the paper. Received June 14, 1995     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

A coupled temperature and strain rate dependent yield function for dynamic deformations of bcc metals

A coupled temperature and strain rate microstructure physically based yield function is proposed in this work. It is incorporated along with the Clausius–Duhem inequality and an appropriate free energy definition in a general thermodynamic framework for deriving a three-dimensional kinematical model for thermo-viscoplastic deformations of body centered cubic (bcc) metals. The evolution equations are expressed in terms of the material time derivatives of the elastic strain, accumulated plastic strain (isotropic hardening), and the back stress conjugate tensor (kinematic hardening). The viscoplastic multipliers are obtained using both the Consistency and Perzyna viscoplasticity models. The athermal yield function is employed instead of the static yield function in the case of the Perzyna viscoplasticity model. It is found that the static strain rate value, at which the material shows rate-independent behavior, varies with the material deformation temperature. Computational aspects of the proposed model are addressed through the finite element implementation with an implicit stress integration algorithm. Finite element simulations are performed by implementing the proposed viscoplasticity constitutive models in the commercial finite element program ABAQUS/Explicit [ABAQUS, 2003. User Manual, Version 6.3. Habbitt, Karlsson and Sorensen Inc., Providence, RI] via the user material subroutine coded as VUMAT. Numerical implementation for a simple compression problem meshed with one element is used to validate the proposed model implementation with applications to tantalum, niobium, and vanadium at low and high strain rates and temperatures. The analysis of a tensile shear banding is also investigated to show the effectiveness and the performance of the proposed framework in describing the strain localizations at high velocity impact. Results show mesh independency as a result of the viscoplastic regularization used in the proposed formulation.     

A viscoplastic strain gradient analysis of materials with voids or inclusions

A finite strain viscoplastic nonlocal plasticity model is formulated and implemented numerically within a finite element framework. The model is a viscoplastic generalisation of the finite strain generalisation by Niordson and Redanz (2004) [Journal of the Mechanics and Physics of Solids 52, 2431–2454] of the strain gradient plasticity theory proposed by Fleck and Hutchinson (2001) [Journal of the Mechanics and Physics of Solids 49, 2245–2271]. The formulation is based on a viscoplastic potential that enables the formulation of the model so that it reduces to the strain gradient plasticity theory in the absence of viscous effects. The numerical implementation uses increments of the effective plastic strain rate as degrees of freedom in addition to increments of displacement. To illustrate predictions of the model, results are presented for materials containing either voids or rigid inclusions. It is shown how the model predicts increased overall yield strength, as compared to conventional predictions, when voids or inclusions are in the micron range. Furthermore, it is illustrated how the higher order boundary conditions at the interface between inclusions and matrix material are important to the overall yield strength as well as the material hardening.     

A strain localization analysis using a viscoplastic softening model for clay

Strain localization has become an attractive subject in geomechanics during the past decade. Shear bands are well known to develop in clay specimens during the straining process. Strain localization is closely related to plastic instability. In the present paper, a non-linear instability condition for the viscoplastic strain softening model during the creep process is firstly obtained. It is found that the proposed viscoplastic model is capable of describing plastic instability. Secondly, a two-dimensional linear instability analysis is performed and the preferred orientation for the growth of fluctuation and the instability condition are derived. It is worth noting that the two instability conditions are equivalent. Finally, the behavior of the clay is numerically analyzed in undrained plane-strain compression tests by the finite element method, considering a transport of pore water in the material at a quasi-static strain rate. The numerical results show that the model can predict strain localization phenomena, such as shear banding. From the numerical calculations, the effects of strain rate and permeability are discussed.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

Coupled thermoviscoplasticity of glassy polymers in the logarithmic strain space based on the free volume theory

The paper outlines a constitutive model for finite thermo-visco-plastic behavior of amorphous glassy polymers and considers details of its numerical implementation. In contrast to existing kinematical approaches to finite plasticity of glassy polymers, the formulation applies a plastic metric theory based on an additive split of Lagrangian Hencky-type strains into elastic and plastic parts. The analogy between the proposed formulation in the logarithmic strain space and the geometrically linear theory of plasticity, makes this constitutive framework very transparent and attractive with regard to its numerical formulation. The characteristic strain hardening of the model is derived from a polymer network model. We consider the particularly simple eight chain model, but also comment on the recently developed microsphere model. The viscoplastic flow rule in the logarithmic strain space uses structures of the free volume flow theory, which provides a highly predictive modeling capacity at the onset of viscoplastic flow. The integration of this micromechanically motivated approach into a three-dimensional computational model is a key concern of this work. We outline details of the numerical implementation of this model, including elements such as geometric pre- and post-transformations to/from the logarithmic strain space, a thermomechanical operator split algorithm consisting of an isothermal mechanical predictor followed by a heat conduction corrector and finally, the consistent linearization of the local update algorithm for the dissipative variables as well as its relationship to the global tangent operator. The performance of the proposed formulation is demonstrated by means of a spectrum of numerical examples, which we compare with our experimental findings.     

黏塑性本构计算的稳定性分析

提出本构方程计算方法的稳定性问题,针对黏塑性本构计算的显式精确算法的稳定性进行分析,发现该算法并非无条件稳定,使用小扰动方法给出了其计算稳定的必要条件,稳定性条件对数值计算中的时间步长提出限制要求。通过有限元算例验证了分析的正确性,计算结果也表明理论推导得到的稳定性公式能够准确预测满足计算稳定性条件要求的最大时间步长与各参数之间关系。     

Compression of a viscoelastic layer between rough parallel plates

If the maximal friction law is applied, then some generalizations of the Prandtl solution for the compression of a plastic layer between rough plates do not exist. In particular, this pertains to the viscoplastic solutions obtained earlier. In the present paper, we show that these solutions do not exist because of the properties of the model material and introduce a model for which this solution can be constructed. The obtained solution is singular. In particular, the equivalent strain rate tends to infinity as the friction surface is approached, and its asymptotic behavior exactly coincides with that arising in the classical solution. The obtained solution is illustrated by numerical examples, which, in particular, show that an extremely thin boundary layer may arise near the friction surfaces.     

Evolution of plastic zones in dynamically loaded plates using different elastic–viscoplastic laws

In the framework of viscoplastic theory many different laws were developed, accounting for material behaviors like creep, relaxation or evolution of overstresses. Though each model is able to predict in uni-axial material tests the values of stresses depending on plastic strains and plastic strain rates the question is if solutions of simulations are still realistic if the viscoplastic law is applied on structural deformations. In the present study strain rate sensitive metal plates are subjected to shock waves. The purpose is to compare simulation results obtained with different elastic–viscoplastic laws to experiments in order to determine the most appropriate material model. By subjecting circular metal plates experimentally to shock wave loadings realistic deformation histories are measured. The measurements are compared to simulation results obtained with different viscoplastic laws. The aim is to find out the accuracy of each model concerning the predictions of displacements, shape formings, spread of plastic zones and evolutions of inner bending moments.     

Design sensitivity analysis for parameters affecting geometry,elastic–viscoplastic material constant and boundary condition by consistent tangent operator-based boundary element method

This paper presents a design sensitivity analysis method by the consistent tangent operator concept-based boundary element implicit algorithm. The design variables for sensitivity analysis include geometry parameters, elastic–viscoplastic material parameters and boundary condition parameters. Based on small strain theory, Perzyna’s elastic–viscoplastic material constitutive relation with a mixed hardening model and two flow functions is considered in the sensitivity analysis. The related elastic–viscoplastic radial return algorithm and the formula of elastic–viscoplastic consistent tangent operator are derived and discussed. Based on the direct differentiation approach, the incremental boundary integral equations and related algorithms for both geometric and elastic–viscoplastic sensitivity analysis are developed. A 2D boundary element program for geometry sensitivity, elastic–viscoplastic material constant sensitivity and boundary condition sensitivity has been developed. Comparison and discussion with the results of this paper, analytical solution and finite element code ANSYS for four plane strain numerical examples are presented finally.     

采用共旋应变的三维热弹塑性有限变形有限元法

本文采用线性化共旋应变张量和增率型虚功原理，建立了有限变形热力耦合弹塑性有限元法。在该方法中，材料的流动应力取为应变总量、应变速率和温度的函数，推导了包含这种函数关系的本构矩阵。另外在温度场分析中，考虑了塑性功和摩擦功转化的热量。文后给出的算例表明该方法可以很好地模拟热加工过程。     

Generalisation de la theorie du potentiel plastique multiple en viscoplasticite†

J. Mandel has extended the classical theory of the plastic potential to materials in which the plastic strain results from several simultaneous slips that are governed by the shear stress for each slip-mechanism (multiple plastic potential). J.R. Rice has shown that if the slip rate on a particular mechanism is also governed by the same law, then the viscoplastic strain-rate tensor can be derived from a viscoplastic potential. These two workers gave only some phenomenological relations.     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

A novel internal dissipation inequality by isotropy and its implication for inelastic constitutive characterization

In the present work a novel inelastic deformation caused internal dissipation inequality by isotropy is revealed. This inequality has the most concise form among a variety of internal dissipation inequalities, including the one widely used in constitutive characterization of isotropic finite strain elastoplasticity and viscoelasticiy. Further, the evolution term describing the difference between the rate of deformation tensor and the “principal rate” of the elastic logarithmic strain tensor is set, according to the standard practice by isotropy, to equal a rank-two isotropic tensor function of the corresponding branch stress, with the tensor function having an eigenspace identical to the eigenspace of the branch stress tensor. Through that a general form of evolution equation for the elastic logarithmic strain is formulated and some interesting and important results are derived. Namely, by isotropy the evolution of the elastic logarithmic strain tensor is embodied separately by the evolutions of its eigenvalues and eigenprojections, with the evolution of the eigenprojections driven by the rate of deformation tensor and the evolution of the eigenvalues connected to specific material behavior. It can be proved that by isotropy the evolution term in the present dissipation inequality stands for the essential form of the evolution term in the extensively applied dissipation inequality.