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.     

An evaluation for several kinematic hardening rules on prediction of multiaxial stress-controlled ratchetting

This study evaluates the performance of several non-linear kinematic hardening rules in predicting the various biaxial ratchetting experiments of stainless steel (SS) 304L under various stress-controlled histories performed by Hassan et al. (2008). The non-linear kinematic hardening rules proposed by 9, 32, 33 and 160, 19, 12 and 13 and the different rules of Abdel-Karim (2009) are examined and carefully scrutinized. The considered kinematic hardening rules range from the simple classical ones to more detailed rules, which incorporate additional terms and/or parameters to simulate different factors that affect ratchetting. It is shown that none of the examined kinematic hardening rules is general enough to simulate all of the ratchetting responses for the experiments under consideration.     

Kinematic hardening in large strain plasticity

A finite strain hyper elasto-plastic constitutive model capable to describe non-linear kinematic hardening as well as non-linear isotropic hardening is presented. In addition to the intermediate configuration and in order to model kinematic hardening, an additional configuration is introduced – the center configuration; both configurations are chosen to be isoclinic. The yield condition is formulated in terms of the Mandel stress and a back-stress with a structure similar to the Mandel stress.It is shown that the non-dissipative part of the plastic velocity gradient not governed by the thermodynamical framework and the corresponding quantity associated with the kinematic hardening influence the material behaviour to a large extent when kinematic hardening is present. However, for isotropic elasticity and isotropic hardening plasticity it is shown that the non-dissipative quantities have no influence upon the stress–strain relation.As an example, kinematic hardening von Mises plasticity is considered, which fulfils the plastic incompressibility condition and is independent of the hydrostatic pressure. To evaluate the response and to examine the influence of the non-dissipative quantities, simple shear is considered; no stress oscillations occur.     

Kinematic hardening rules in uncoupled modeling for multiaxial ratcheting simulation

An earlier paper by the authors evaluated the performance of several coupled models in simulating a series of uniaxial and biaxial ratcheting responses. This paper evaluates the performance of various kinematic hardening rules in an uncoupled model for the same set of ratcheting responses. A modified version of the Dafalias–Popov uncoupled model has been demonstrated to perform well for uniaxial ratcheting simulation. However, its performance in multiaxial ratcheting simulation is significantly influenced by the kinematic hardening rules employed in the model. Performances of eight different kinematic hardening rules, when engaged with the modified Dafalias–Popov model, are evaluated against a series of rate-independent multiaxial ratcheting responses of cyclically stabilized carbon steels. The kinematic hardening rules proposed by Armstrong–Frederick, Voyiadjis–Sivakumar, Phillips, Tseng–Lee, Kaneko, Xia–Ellyin, Chaboche and Ohno–Wang are examined. The Armstrong–Frederick rule performs reasonably for one type of the biaxial ratcheting response, but fails in others. The Voyiadjis–Sivakumar rule and its constituents, the Phillips and the Tseng–Lee rules, can not simulate the biaxial ratcheting responses. The Kaneko rule, composed of the Ziegler and the prestress directions, and the Xia–Ellyin rule, composed of the Ziegler and Mroz directions, also fail to simulate the biaxial ratcheting responses. The Chaboche rule, with three decomposed Armstrong–Frederick rules, performs the best for the whole set of ratcheting responses. The Ohno–Wang rule performs well for the data set, except for one biaxial response where it predicts shakedown with subsequent reversal of ratcheting.     

Effect of elastic modulus variation during plastic deformation on uniaxial and multiaxial ratchetting simulations

In order to identify different variables that affect ratchetting simulations, variation of elastic modulus during loading and unloading is considered and discussed based on the experimental observations which pointed out by Morestin and Boivin, 1996, Ishikawa, 1997, Cleveland and Ghosh, 2002, Zhou et al. (2005) and recently by Khan et al., 2009a, Khan et al., 2009b, Khan et al., 2009c. Then the effect of such variation on simulations is scrutinized from the theoretical point of view by considering simulations of ratchetting experiments conducted on stainless steel 304L by Hassan et al. (2008) using the well-known Armstrong–Frederick model. It is shown that, using two different values for the elastic modulus during loading and unloading could have a significant effect on simulations of uniaxial ratchetting. On the other hand, such significant effect hardly occurs in the case of simulations of biaxial ratchetting experiments under consideration. The importance of such findings is that the excessive ratchetting over-prediction resulting from any specific kinematic hardening rule is expected to decrease significantly by taking into consideration this effect. In this case, modeling of kinematic hardening rules could necessitate more attention and reconsideration.     

Cyclic ratchetting of 1070 steel under multiaxial stress states

Cyclic ratchetting behavior of 1070 steel is studied under proportional and nonproportional loading with specific emphasis on the ratchetting rate decay mechanisms for large numbers of loading cycles. Under proportional loading, where the principal stress directions are unchanged, the ratchetting evolves in the mean stress direction. Under nonproportional loading, however, the ratchetting direction is determined by the loading path and can be different from the mean stress direction. The ratchetting rate decreases with increasing loading cycles, displaying a power law relationship with the number of loading cycles. The experimental ratchetting results indicate that under cyclic loading the material exhibits a tendency toward complying with a linear hardening rule with concomitant hysteresis loop closure. Based on the fundamental framework of plasticity theory and detailed evaluation of the stress-strain behaviors, the ratchetting can be classified into two basic types; Type I, which is identifiable with proportional loading where the ratchetting is due to the different values of the plastic modulus function at the symmetric loading points with respect to the mean stress state, and Type II, which represents nonproportional loading where the ratchetting is driven by the noncoincidence of the plastic strain rate vector and the translation direction of the yield surface (backstress rate vector). The Armstrong-Frederick-based plasticity models modified by Chaboche et al. and Bower are ill-suited for describing the experimental results of both types of ratchetting. The Ohno-Wang model, which introduces a threshold concept, can account for the ratchetting rate decay of Type II ratchetting, providing results that agree with experimental observations. Modification may be needed for the Ohno-Wang model so that the model can better describe Type I ratchetting.     

Large strain responses of elastic-perfect plasticity and kinematic hardening plasticity with the logarithmic rate: Swift effect in torsion

A new Eulerian rate type elastic-perfectly plastic model has recently been established by utilizing the newly discovered logarithmic rate. It has been proved that this model is unique among the objective elastic-perfectly plastic models with all objective corotational stress rates and other known objective stress rates by virtue of the self-consistency criterion: the hypoelastic formulation intended for elastic behaviour must be exactly integrable to deliver a hyperelastic relation. The finite simple shear response of this model has been studied and shown to be reasonable for both shear and normal stress components. On the other hand, a kinematic hardening plasticity model may be formulated by adopting the logarithmic rate. The objective of this work is to further study the large deformation responses of the foregoing two kinds of idealized models, in particular the well-known Swift effect, in torsion of thin-walled cylindrical tubes. A complete, rigorous analysis is made for the orders of magnitude of all stress components. A closed-form solution is obtained for the kinematic hardening plastic case, and an analytical perturbation solution is derived for the elastic-perfectly plastic case. It is shown that the simple idealized kinematic hardening model with the logarithmic rate, which uses only two classical material constants, i.e., the initial (tensile) yield stress and the hardening modulus, may arrive at satisfactory explanation for and reasonable accord with salient features of experimental observation.     

Constitutive modelling of cyclic plasticity of metal particulate-reinforced brittle matrix composites under compression-compression loading

The Mori-Tanaka approach is used to modelling metal particulate-reinforced brittle matrix composites under cyclic compressive loading. The J₂-flow theory is considered as the relevant physical law of plastic flow in inclusions. Ratchetting of the composite is prevented by the strong constraint exerted by the matrix on the inclusions, even under the assumption of evanescent kinematic hardening. However, the weakening constraint power of the matrix caused by microfracture damage around inclusions is closely coupled with the plasticity of inclusion and leads to ratchetting even when the plastic deformation of inclusions is described by an isotropic hardening rule. A detailed parametric study has revealed that ratchetting is followed by either plastic or elastic shakedown, depending on the load amplitude, composite parameters and the mean length of microcracks.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

Ratcheting and wrinkling of tubes due to axial cycling under internal pressure: Part II analysis

Part I presented a set of experiments in which pressurized tubes were cycled axially under stress control about a compressive mean stress. This loading history causes biaxial ratcheting involving compressive axial strain and expansion of the diameter of the tube. The compressive strain in turn induces the initiation and growth of axisymmetric wrinkles. Persistent cycling resulted in localization of the wrinkles and collapse. In Part II the problem is first modeled as a shell with initial axisymmetric imperfections while the biaxial ratcheting of the material is modeled using the Dafalias–Popov two-surface nonlinear kinematic hardening model. It is demonstrated that when suitably calibrated this modeling framework reproduces the prevalent ratcheting deformations and the evolution of wrinkling including the conditions at collapse accurately for all experiments. The calibrated model is then used to evaluate the ratcheting behavior of pipes under thermal-pressure cyclic loading histories experienced by axially restrained pipelines.     

Strain gradient plasticity modeling of the cyclic behavior of laminate microstructures

Two recently proposed Helmholtz free energy potentials including the full dislocation density tensor as an argument within the framework of strain gradient plasticity are used to predict the cyclic elastoplastic response of periodic laminate microstructures. First, a rank-one defect energy is considered, allowing for a size-effect on the overall yield strength of micro-heterogeneous materials. As a second candidate, a logarithmic defect energy is investigated, which is motivated by the work of Groma et al. (2003). The properties of the back-stress arising from both energies are investigated in the case of a laminate microstructure for which analytical as well as numerical solutions are derived. In this context, a new regularization technique for the numerical treatment of the rank-one potential is presented based on an incremental potential involving Lagrange multipliers. The results illustrate the effect of the two energies on the macroscopic size-dependent stress–strain response in monotonic and cyclic shear loading, as well as the arising pile-up distributions. Under cyclic loading, stress–strain hysteresis loops with inflections are predicted by both models. The logarithmic potential is shown to provide a continuum formulation of Asaro's type III kinematic hardening model. Experimental evidence in the literature of such loops with inflections in two-phased FFC alloys is provided, showing that the proposed strain gradient models reflect the occurrence of reversible plasticity phenomena under reverse loading.     

Ratchetting under tension–torsion loadings: experiments and modelling

This paper is concerned with the mechanical behaviour of 316 austenitic stainless steel under multiaxial loadings and particular attention is paid to ratchetting under tension–torsion non-proportional loadings. First, a series of uniaxial tests and biaxial tests has been carried out in order to calibrate five different cyclic plasticity models based on an isotropic hardening rule and a non-linear kinematic hardening rule. It is shown that this class of models gives quite good agreement between the experimental and numerical results. Second, another series of ratchetting tests has been carried out under tension–torsion loadings in order to test the prediction capacities of the previous models. It is shown that whereas the models have been calibrated with similar loading paths, four of the five selected models give poor predictions.     

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.     

Ratchetting of viscoplastic material with cyclic softening,part 2: application of constitutive models

Simulation capability on ratchetting of modified 9Cr–1Mo steel at 550 °C was discussed using several constitutive models in the present paper. It was revealed that the authors' previous model, which uses an Armstrong–Frederick kinematic hardening rule, has a strong tendency to overestimate both uniaxial and multiaxial ratchetting of the material. On the contrary, the Ohno–Wang (OW) I model tended to underestimate the uniaxial and multiaxial ratchetting. The OW II and III models predicted the uniaxial and multiaxial ratchetting with better accuracy. Regarding the uniaxial ratchetting under the zero mean stress condition described in part 1 of this study, none of the constitutive models was able to simulate it even qualitatively. On the basis of the OW I model, a constitutive model incorporating a tension–compression asymmetry was proposed to predict the ratchetting behavior under the zero mean stress condition. The simulation capability of the proposed model was discussed in comparison with that of the other constitutive models.     

Anatomy of coupled constitutive models for ratcheting simulation

This paper critically evaluates the performance of five constitutive models in predicting ratcheting responses of carbon steel for a broad set of uniaxial and biaxial loading histories. The models proposed by Prager, Armstrong and Frederick, Chaboche, Ohno-Wang and Guionnet are examined. Reasons for success and failure in simulating ratcheting by these models are elaborated. The bilinear Prager and the nonlinear Armstrong-Frederick models are found to be inadequate in simulating ratcheting responses. The Chaboche and Ohno-Wang models perform quite well in predicting uniaxial ratcheting responses; however, they consistently overpredict the biaxial ratcheting responses. The Guionnet model simulates one set of biaxial ratcheting responses very well, but fails to simulate uniaxial and other biaxial ratcheting responses. Similar to many earlier studies, this study also indicates a strong influence of the kinematic hardening rule or backstress direction on multiaxial ratcheting simulation. Incorporation of parameters dependent on multiaxial ratcheting responses, while dormant for uniaxial responses, into Chaboche-type kinematic hardening rules may be conducive to improve their multiaxial ratcheting simulations. The uncoupling of the kinematic hardening rule from the plastic modulus calculation is another potentially viable alternative. The best option to achieve a robust model for ratcheting simulations seems to be the incorporation of yield surface shape change (formative hardening) in the cyclic plasticity model.     

The appropriate corotational rate, exact formula for the plastic spin and constitutive model for finite elastoplasticity

The exact formulae for the plastic and the elastic spin referred to the deformed configuration are derived, where the plastic spin is a function of the plastic strain rate and the elastic spin a function of the elastic strain rate. With these exact formulae we determine the macroscopic substructure spin that allows us to define the appropriate corotational rate for finite elastoplasticity.Plastic, elastic and substructure spin are considered and simplified for various sub-classes of restricted elastic-plastic strains. It is shown that for the special cases of rigid-plasticity and hypoelasticity the proposed corotational rate is identical with the Green-Naghdi rate, while the ZarembaJaumann rate yields a good approximation for moderately large strains.To compare our exact plastic spin formula with the constitutive assumption for the plastic spin introduced by Dafalias and others, we simplify our result for small elastic-moderate plastic strains and introduce a simplest evolution law for kinematic hardening leading to the Dafalias formula and to an exact determination of its unknown coefficient. It is also shown that contrary to statements in the literature the plastic spin is not zero for vanishing kinematic hardening.For isotropic-elastic material with induced plastic flow undergoing isotropic and kinematic hardening constitutive and evolution laws are proposed. Elastic and plastic Lagrangean and Eulerian logarithmic strain measures are introduced and their material time derivatives and corotational rates, respectively, are considered. Finally, the elastic-plastic tangent operator is derived.The presented theory is implemented in a solution algorithm and numerically applied to the simple shear problem for finite elastic-finite plastic strains as well as for sub-classes of restricted strains. The results are compared with those of the literature and with those obtained by using other corotational rates.     

A plasticity model for pressure-dependent anisotropic cellular solids

The initial and subsequent yield surfaces for an anisotropic and pressure-dependent 2D stochastic cellular material, which represents solid foams, are investigated under biaxial loading using finite element analysis. Scalar measures of stress and strain, namely characteristic stress and characteristic strain, are used to describe the constitutive response of cellular material along various stress paths. The coupling between loading path and strain hardening is then investigated in characteristic stress–strain domain. The nature of the flow rule that best describes the plastic flow of cellular solid is also investigated. An incremental plasticity framework is proposed to describe the pressure-dependent plastic flow of 2D stochastic cellular solids. The proposed plasticity framework adopts the anisotropic and pressure-dependent yield function recently introduced by Alkhader and Vural [Alkhader M., Vural M., 2009a. An energy-based anisotropic yield criterion for cellular solids and validation by biaxial FE simulations. J. Mech. Phys. Solids 57(5), 871–890]. It has been shown that the proposed yield function can be simply calibrated using elastic constants and flow stresses under uniaixal loading. Comparison of stress fields predicted by continuum plasticity model to the ones obtained from FE analysis shows good agreement for the range of loading paths and strains investigated.     

Constitutive modelling of nonproportional cyclic plasticity

A multiplicative hardening function and a unified evolution rule of the hardening factors are proposed. The hardening factorf ₁ is introduced to describe cyclic hardening with respect to the plastic strain range, whilef ₂ andf ₃ describe, respectively, instantaneous and hereditary additional hardening with respect to the nonproportionality of the plastic strain path. Two material dependent memory parametersa ₁ anda ₃ are introduced to keep the memory of the largest cyclic and additional hardening in the previous plastic deformation history. Different hardening mechanisms are then embedded into a thermomechanically consistent constitutive equation through the hardening function. The constitutive response of 304 and 316 stainless steels subjected to biaxial nonproportional cyclic loading is analyzed and the proposed model is critically verified by comparing the results with experimental results obtained by Tanaka et al., and Ohashi et al. The project supported by National Natural Science Foundation of China