Cyclic multiaxial and shear finite deformation responses of OFHC Cu. Part II: An extension to the KHL model and simulations

In this part, the Khan–Huang–Liang (KHL) constitutive model was extended to account for kinematic hardening characteristic behavior of materials. The extended model is then generalized and used to simulate experimental response of oxygen free high conductivity (OFHC) copper under cyclic shear straining and biaxial tension–torsion (multiaxial ratchetting) experiments presented in Part I (Khan et al., 2007). In addition, a new modification for the non-linear kinematic hardening rule of Karim–Ohno (Abdel-Karim and Ohno, 2000) is proposed to simulate multiaxial ratchetting behaviors. Although, the kinematic hardening contributes the most to the response, it is shown that, the loading rate effect, and a coupled isotropic and kinematic hardening effect should also be considered while simulating the multiaxial ratchetting behavior of OFHC copper. Furthermore, the newly modified kinematic hardening rules is able to fairly well simulate the multiaxial ratchetting experiments under different loading conditions, irrespective of the value of applied axial tensile stress, shear strain amplitude, pre-cyclic hardening and/or loading sequence.     

爆炸载荷下钢板层裂的二维数值计算

采用二维有限元方法对爆炸载荷下钢板的层裂问题进行数值模拟计算。计算中引入了一种普适的弹塑性本构方程,并采用了损伤积累准则和半解耦的损伤本构模型。结果表明,二维计算的主裂片厚度略小于一维计算,裂片飞散速度要比一维计算低得多,更接近实际结果。     

一般加载规律的弹塑性本构关系

将有关文献给出一般加载规律一维全量理论的简单模型推广到一般加载规律的一维增量理论,进而推广到一般加载规律的多维增量理论,在此基础上,建立了推导一般加载规律的多维增量理论的本构关系的一种途径。应用这种途径,从应力空间的加载函数和应变空间的加载函数出发,推导了等向强化材料和被加热的等向强化材料的一般加载规律的弹塑性本构关系的两种表示形式。理论和实例均表明,这种途径对等向强化材料、随动强化材料和理想弹塑性材料均适用。     

Numerical plane stress elastic–perfectly plastic crack analysis under Tresca yield condition with comparison to Dugdale plastic strip model

A number of plane stress numerical analyses of the mode I elastoplastic fracture mechanics problem have been performed in the past using the Huber–Mises yield criterion. This study employs instead the Tresca yield condition using an incremental theory of plasticity for a stationary crack. A commercial finite element program is used to solve the opening mode of fracture problem (mode I) for a square plate containing a central crack under generalized plane stress loading conditions. A biaxial uniform tensile traction is applied to the edges of a thin plate composed of a linear elastic non-work hardening material under small strain assumptions. The finite element results are compared with the analytical predictions of the Dugdale plastic strip model for a crack in an infinite plate subject to a biaxial uniform load at infinity.     

316L不锈钢随动强化模型改进研究

本文对316L不锈钢进行了单轴与多轴非比例路径下的应力控制棘轮试验,考察了应力幅值、平均应力和加载历程对棘轮特性的影响。同时进行了应变控制循环试验以研究材料的应力松弛特性。试验结果表明轴向棘轮效应在对称剪切荷载下效果明显,同时棘轮应变随应力幅值和平均应力的增加而增加。研究了Chen-Jiao随动强化模型与Jiang-Sehitoglu随动强化模型采用的单轴与多轴参数对背应力分量增量方向的影响,将Chen-Jiao模型中的多轴系数替换为界面饱和率,并在此基础上引入新的参数对塑性模量系数进行修正,计算结果表明修正后的模型能提升应力控制下多轴棘轮的预测精度,并能很好的预测应力松弛现象,表明了新模型的正确性与有效性。     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.     

A two-surface model describing ratchetting behaviors and transient hardening under nonproportional loading

A new superposed rule of Mroz's kinematic hardening rule and Ziegler's kinematic hardening rule based on two-surface model is proposed in the paper. Some experimental results on ratchetting of 2014-T6 aluminum alloy are predicted very well under multiaxial loading. In addition the conformability of the model is discussed for transient cyclic hardening under two kinds of nonproportional cyclic loading paths, i.e. square and rhombic path. The project supported by the National Natural Science Foundation of China     

An advancement in cyclic plasticity modeling for multiaxial ratcheting simulation

In a search for a constitutive model for ratcheting simulations, the models by Chaboche, Ohno–Wang, McDowell, Jiang–Sehitoglu, Voyiadjis–Basuroychowdhury and AbdelKarim–Ohno are evaluated against a set of uniaxial and biaxial ratcheting responses. With the assumption of invariant shape of the yield surface during plastic loading, the ratcheting simulations for uniaxial loading are primarily a function of the plastic modulus calculation, whereas the simulations for multiaxial loading are sensitive to the kinematic hardening rule of a model. This characteristic of the above mentioned models is elaborated in this paper. It is demonstrated that if all parameters of the kinematic hardening rule are determined from uniaxial responses only, these parameters primarily enable a better plastic modulus calculation. However, in this case the role of the kinematic hardening rule in representing the ratcheting responses for multiaxial loading is under-appreciated. This realization motivated many researchers to incorporate multiaxial load dependent terms or parameters into the kinematic hardening rule. This paper evaluates some of these modified rules and finds that none is general enough to simulate the ratcheting responses consistently for the experiments considered. A modified kinematic hardening rule is proposed using the idea of Delobelle and his co-workers in the framework of the Chaboche model. This new rule introduces only one multiaxial load dependent parameter to the Chaboche model, but performs the best in simulating all the ratcheting responses considered.     

On shakedown analysis in hardening plasticity

The extension of classical shakedown theorems for hardening plasticity is interesting from both theoretical and practical aspects of the theory of plasticity. This problem has been much discussed in the literature. In particular, the model of generalized standard materials gives a convenient framework to derive appropriate results for common models of plasticity with strain-hardening. This paper gives a comprehensive presentation of the subject, in particular, on general results which can be obtained in this framework. The extension of the static shakedown theorem to hardening plasticity is presented at first. It leads by min-max duality to the definition of dual static and kinematic safety coefficients in hardening plasticity. Dual static and kinematic approaches are discussed for common models of isotropic hardening of limited or unlimited kinematic hardening. The kinematic approach also suggests for these models the introduction of a relaxed kinematic coefficient following a method due to Koiter. Some models for soils such as the Cam-clay model are discussed in the same spirit for applications in geomechanics. In particular, new appropriate results concerning the variational expressions of the dual kinematic coefficients are obtained.     

A three-dimensional model of anisotropic hardening in metals and its application to the analysis of sheet metal formability

The hardening model proposed by Z. Mróz based on the uniaxial fatigue behavior of many metals is adopted to derive an incremental constitutive equation for general three-dimensional problems. This constitutive law is then employed in the analysis of metal forming problems to assess the influence of loading cycles, of the types involved in standard forming processes, on the ultimate formability of sheet metals. The predicted forming limit curves differ quantitatively from results obtained via an isotropie hardening model and differ qualitatively from those obtained via a kinematic model. Also investigated are the effects of such loading cycles on material response to simple tensile loading, which is often used to characterize a material. Significant differences between the present model and the other two models considered are observed in such characterizers of simple tensile behavior as the stress-strain curve, the anisotropy parameter and the uniform elongation. These differences suggest a rather simple experiment to identify the proper material model to be used in analyses of problems which involve loading cycles. Comparisons with some experimental results reveal that the employment of an anisotropic hardening model, such as the generalized Mróz model derived herein, is indeed crucial in accurately predicting material response to complicated loading histories.     

Constitutive modeling of strain range dependent cyclic hardening

In this paper, a new approach for constitutive modeling of strain range dependent cyclic hardening is proposed by extending the kinematic hardening model based on the critical state of dynamic recovery. It is assumed that isotropic, as well as kinematic, hardening consists of several parts, and that each part of isotropic hardening evolves when the corresponding part of kinematic hardening is in the critical state of dynamic recovery. The extended model is capable of simulating the cyclic hardening behavior in which different characteristics of cyclic hardening appear depending on strain range. The model is verified by simulating the relatively large cyclic straining tests of 304 stainless steel at ambient temperature, in which cyclic hardening does not stabilize before rupture if strain range exceeds a certain value. The model is further verified by predicting the history dependence of cyclic hardening under incremental cyclic loading and the maximum plastic strain dependence of strain hardening in cyclic tension.     

An elasto-plastic damage constitutive theory and its prediction of evolution of subsequent yield surfaces and elastic constants

Based on pair functional potentials, Cauchy-Born rule and slip mechanism, a material model assembling with spring-bundle components, a cubage component and slip components is established to describe the elasto-plastic damage constitutive relation under finite deformation. The expansion/shrink, translation and distortion of yield surfaces can be calculated based on the hardening rule and Bauschinger effect defined on the slip component level. Both kinematic and isotropic hardening are included. Numerical simulations and predictions under tension, torsion, and combined tension-torsion proportional/non-proportional loading are performed to obtain the evolution of subsequent yield surfaces and elastic constants and compare with two sets of experimental data in literature, one for a very low work hardening aluminum alloy Al 6061-T6511, and another for a very high work hardening aluminum alloy annealed 1100 Al. The feature of the yield surface in shape change, which presents a sharp front accompanied by a blunt rear under proportional loading, is described by the latent hardening and Bauschinger effect of slip components. Further, the evolution law of subsequent yield surfaces under different proportional loading paths is investigated in terms of their equivalence. The numerical simulations under non-proportional loading conditions for annealed 1100 Al are performed, and the subsequent yield surfaces exhibit mixed cross effect because the kinematic hardening and isotropic hardening follow different evolution tendency when loading path changes. The results of non-proportional loading demonstrate that the present model has the ability to address the issue of complex loading due to the introduction of state variables on slip components. Moreover, as an elasto-plastic damage constitutive model, the present model can also reflect the variation of elastic constants through damage defined on the spring-bundle components.     

On the performance of kinematic hardening rules in predicting a class of biaxial ratcheting histories

It is well known that strain-symmetric axial cycling of thin-walled metal tubes in the presence of pressure results in a progressive accumulation (ratcheting) of circumferential strain. It was previously demonstrated that the prediction of the rate of ratcheting under constant internal pressure, by nonlinear kinematic hardening models, is very sensitive to the hardening rule adopted. It was shown that the Armstrong-Frederick hardening suitably calibrated and used in a class of models can yield reasonably good predictions of the rate of ratcheting for a range of cycle parameters. In this paper, the subject is revisited in the light of further experimental results involving simultaneous cycling of the internal pressure and the axial strain. Experiments and analyses were performed for a family of five such biaxial loading histories. A similar sensitivity to the kinematic hardening rule used in the models was observed in all the new loading histories. Furthermore the hardening rule calibrated in the constant pressure experiments was found to yield accurate predictions for three of the loading histories considered and poor predictions for the other two. The reasons for this varied performance are analyzed and some recommendations for implementation of such models in structural applications are made.     

Shakedown theorems for elastic–plastic solids in the framework of gradient plasticity

Static and kinematic shakedown theorems are given for a class of generalized standard materials endowed with a hardening saturation surface in the framework of strain gradient plasticity. The so-called residual-based gradient plasticity theory is employed. The hardening law admits a hardening potential, which is a C¹-continuous function of a set of kinematic internal variables and of their spatial gradients, and is required to satisfy a global sign restriction (but not to be necessarily convex). The totally produced, the accumulated and the freely moving dislocations per unit volume, distinguished as statistically stored and geometrically necessary ones, are in this way accounted for. Like for a generalized standard material, the shakedown safety factor is found to depend on the (generally size dependent) yield and saturation limits, but not on the particular differential-type hardening law of the material.     

Effect of the rate dependence of nonlinear kinematic hardening rule on relaxation behavior

Relaxation experiments for metallic materials and solid polymers have exhibited nonlinear dependence of stress relaxation on prior loading rate; the relaxed stress associated with the fastest prior strain rate has the smallest magnitude at the end of the same relaxation periods. Modeling capability for the basic feature of relaxation behavior is qualitatively investigated in the context of unified state variable theory. Unified constitutive models are categorized into three general classes according to the rate dependence of kinematic hardening rule, which defines the evolution of the back (equilibrium) stress and is the major difference among constitutive models. The first class of models adopts the nonlinear kinematic hardening rule proposed by Armstrong and Frederick. In this class, the back stress appears to be rate-independent under loading and subsequent relaxation conditions. In the second class of models, a stress rate term is incorporated into the Armstrong–Frederick rule and the back stress then becomes rate-dependent during relaxation condition even though it remains rate-independent under loading condition. The final class proposed here includes a new nonlinear kinematic hardening rule that causes the back stress to be rate-dependent all the time. It is shown that the apparent rate dependence of the back stress during relaxation enables constitutive models to predict the influence of prior loading rate on relaxation behavior.     

Multiaxial ratcheting with advanced kinematic and directional distortional hardening rules

Ratcheting is defined as the accumulation of plastic strains during cyclic plastic loading. Modeling this behavior is extremely difficult because any small error in plastic strain during a single cycle will add to become a large error after many cycles. As is typical with metals, most constitutive models use the associative flow rule which states that the plastic strain increment is in the direction normal to the yield surface. When the associative flow rule is used, it is important to have the shape of the yield surface modeled accurately because small deviations in shape may result in large deviations in the normal to the yield surface and thus the plastic strain increment in multi-axial loading. During cyclic plastic loading these deviations will accumulate and may result in large errors to predicted strains.This paper compares the bi-axial ratcheting simulations of two classes of plasticity models. The first class of models consists of the classical von Mises model with various kinematic hardening (KH) rules. The second class of models introduce directional distortional hardening (DDH) in addition to these various kinematic hardening rules. Directional distortion describes the formation of a region of high curvature on the yield surface approximately in the direction of loading and a region of flattened curvature approximately in the opposite direction. Results indicate that the addition of directional distortional hardening improves ratcheting predictions, particularly under biaxial stress controlled loading, over kinematic hardening alone.     

Exact integration of constitutive equations of kinematic hardening materials and its extended applications

In this paper, an exact formula for the integration of the constitutive equations of kinematic hardening material is presented. Its algorithms are simple and clear. For isotropic hardening or mixed hardening material, the formula is still an exact solution for the case of radial loading, and it is an approximate solution with reasonable accuracy for the case of non-radial loading. The computation results show that the procedure proposed in this paper improves both accuracy and efficiency of numerical integration schemes adopted widely in elastic-plastic finite element analysis.     

Deformation gradient based kinematic hardening model

A kinematic hardening model applicable to finite strains is presented. The kinematic hardening concept is based on the residual stresses that evolve due to different obstacles that are present in a polycrystalline material, such as grain boundaries, cross slips, etc. Since these residual stresses are a manifestation of the distortion of the crystal lattice a corresponding deformation gradient is introduced to represent this distortion. The residual stresses are interpreted in terms of the form of a back-stress tensor, i.e. the kinematic hardening model is based on a deformation gradient which determines the back-stress tensor. A set of evolution equations is used to describe the evolution of the deformation gradient. Non-dissipative quantities are allowed in the model and the implications of these are discussed. Von Mises plasticity for which the uniaxial stress–strain relation can be obtained in closed form serves as a model problem. For uniaxial loading, this model yields a kinematic hardening identical to the hardening produced by isotropic exponential hardening. The numerical implementation of the model is discussed. Finite element simulations showing the capabilities of the model are presented.