Thermodynamic based model for the evolution equation of the backstress in cyclic plasticity

A nonlinear kinematic hardening rule is developed here within the framework of thermodynamic principles. The derived kinematic hardening evolution equation has three distinct terms: two strain hardening terms and a dynamic recovery term that operates at all times. The proposed hardening rule, which is referred in this paper as the FAPC (Fredrick and Armstrong–Phillips–Chaboche) kinematic hardening rule, shows a combined form of the Frederick and Armstrong backstress evolution equation, Phillips evolution equation, and Chaboche series rule. A new term is incorporated into the Frederick and Armstrong evolution equation that appears to have agreement with the experimental observations that show the motion of the center of the yield surface in the stress space is directed between the gradient to the surface at the stress point and the stress rate direction at that point. The model is further modified in order to simulate nonproportional cyclic hardening by proposing a measure representing the degree of nonproportionality of loading. This measure represents the topology of the incremental stress path. Numerically, it represents the angle between the current stress increment and the previous stress increment, which is interpreted through the material constants of the kinematic hardening evolution equation. This new kinematic hardening rule is incorporated in a material constitutive model based on the von Mises plasticity type and the Chaboche isotropic hardening type. Numerical integration of the incremental elasto-plastic constitutive equations is based on a simple semi-implicit return-mapping algorithm and the full Newton–Raphson iterative method is used to solve the resulting nonlinear equations. Experimental simulations are conducted for proportional and non-proportional cyclic loadings. The model shows good correlation with the experimental results.     

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.     

Time-dependent behavior of passive skeletal muscle

An isotropic three-dimensional nonlinear viscoelastic model is developed to simulate the time-dependent behavior of passive skeletal muscle. The development of the model is stimulated by experimental data that characterize the response during simple uniaxial stress cyclic loading and unloading. Of particular interest is the rate-dependent response, the recovery of muscle properties from the preconditioned to the unconditioned state and stress relaxation at constant stretch during loading and unloading. The model considers the material to be a composite of a nonlinear hyperelastic component in parallel with a nonlinear dissipative component. The strain energy and the corresponding stress measures are separated additively into hyperelastic and dissipative parts. In contrast to standard nonlinear inelastic models, here the dissipative component is modeled using an evolution equation that combines rate-independent and rate-dependent responses smoothly with no finite elastic range. Large deformation evolution equations for the distortional deformations in the elastic and in the dissipative component are presented. A robust, strongly objective numerical integration algorithm is used to model rate-dependent and rate-independent inelastic responses. The constitutive formulation is specialized to simulate the experimental data. The nonlinear viscoelastic model accurately represents the time-dependent passive response of skeletal muscle.     

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.     

A unified constitutive law for cyclic viscoplasticity

A phenomenological constitutive model for cyclic viscoplasticity is presented within the framework of unified state variable theory. The model utilizes three distinct parameters to account for the cyclic (isotropic) hardening: drag stress, isotropic stress and rate sensitivity parameter causing the back stress to be rate-dependent. With the help of a rate-dependent format of the back stress, the constitutive model can reproduce positive, zero and negative strain rate sensitivity of the flow stress in a concise manner. The purpose of the paper is to investigate the influence of the three parameters on stress relaxation behavior and rate-dependent cyclic hardening concerned with the variation in strain rate sensitivity. The applicability of the model to monotonic and cyclic loading is validated by comparing the predictions with experiments of two stainless steels and an aluminum alloy reported in literature.     

Experimental identification and mathematical modeling of viscoplastic material behavior

Uniaxial torsion and biaxial torsion-tension experiments on thin-walled tubes were carried out to investigate the viscoplastic behavior of stainless steel XCrNi18.9. A series of monotonic tests under strain and stress control shows nonlinear rate dependence and suggests the existence of equilibrium states, which are asymptotically approached during relaxation and creep processes. Strain controlled cyclic experiments display various hardening and softening phenomena that depend on strain amplitude and mean strain. All experiments indicate that the equilibrium states within the material depend on the history of the input process, whereas the history-dependence of the relaxation and creep behavior appears less significant. From the experiments the design of a constitutive model of viscoplasticity is motivated: The basic assumption is a decomposition of the total stress into an equilibrium stress and a non-equilibrium overstress: At constant strain, the overstress relaxes to zero, where the relaxation time depends on the overstress in order to account for the nonlinear rate-dependence. The equilibrium stress is assumed to be a rate independent functional of the total strain history. Classical plasticity is utilized with a kinematic hardening rule of the Armstrong-Frederick type. In order to incorporate the amplitude-dependent hardening and softening behavior, a generalized arc length representation is applied [14]. The introduction of an additional kinematic hardening variable facilitates consideration of additional hardening effects resulting from the non-radiality of the input process. Apart from the common yield and loading criterion of classical plasticity, the proposed constitutive model does not contain any further distinction of different cases.The experimental data are sufficient to identify the material parameters of the constitutive model. The results of the identification procedure demonstrate the ability of the model to represent the observed phenomena with satisfactory approximation.     

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.     

Experiments and Material Parameter Identification Using Finite Elements. Uniaxial Tests and Validation Using Instrumented Indentation Tests

In this article, we focus our attention on the relation between instrumented indentation tests and the prediction by means of finite element calculations. To this end, a finite strain viscoplasticity model of Perzyna-type with non-linear isotropic and kinematic hardening is calibrated at experimental data of steel S690QL. A particular concept for conducting uniaxial tensile and compression tests is taken up in order to represent the basic rate-dependent material behavior. In this respect, an algorithmic framework of material parameter identification using finite elements is proposed leading to a two-stage procedure in the case of the underlying rate-dependent constitutive model. On the basis of the termination points of relaxation the rate-independent equilibrium stress state can be identified and all viscous parts of the model are obtained using rate-dependent loading paths. Finally, use is made of finite elements for predicting indentation experiments, which results in a critical view on modeling and parameter identification on the basis of experimental results occurring in instrumented indentation tests.     

Constitutive modeling of cyclic plasticity deformation of a pure polycrystalline copper

Cyclic plasticity experiments were conducted on a pure polycrystalline copper and the material was found to display significant cyclic hardening and nonproportional hardening. An effort was made to describe the cyclic plasticity behavior of the material. The model is based on the framework using a yield surface together with the Armstrong–Frederick type kinematic hardening rule. No isotropic hardening is considered and the yield stress is assumed to be a constant. The backstress is decomposed into additive parts with each part following the Armstrong–Frederick type hardening rule. A memory surface in the plastic strain space is used to account for the strain range effect. The Tanaka fourth order tensor is used to characterize nonproportional loading. A set of material parameters in the hardening rules are related to the strain memory surface size and they are used to capture the strain range effect and the dependence of cyclic hardening and nonproportional hardening on the loading magnitude. The constitutive model can describe well the transient behavior during cyclic hardening and nonproportional hardening of the polycrystalline copper. Modeling of long-term ratcheting deformation is a difficult task and further investigations are required.     

Constitutive relations for isotropic or kinematic hardening at finite elastic–plastic deformations

The rate-type constitutive relations of rate-independent metals with isotropic or kinematic hardening at finite elastic–plastic deformations were presented through a phenomenological approach. This approach includes the decomposition of finite deformation into elastic and plastic parts, which is different from both the elastic–plastic additive decomposition of deformation rate and Lee’s elastic–plastic multiplicative decomposition of deformation gradient. The objectivity of the constitutive relations was dealt with in integrating the constitutive equations. A new objective derivative of back stress was proposed for kinematic hardening. In addition, the loading criteria were discussed. Finally, the stress for simple shear elastic–plastic deformation was worked out.     

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

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

An inelastic material model for filled polytetrafluorethylen

Summary  This paper presents a viscoplastic model for PTFE designed to simulate numerically PTFE shaft seals. A rate-independent elastoplastic model with an endochronic flow rule is coupled in series with a rate-dependent Kelvin model, which has a highly nonlinear damper. In contrast to previous models for PTFE, this unified approach is suitable for numerical simulation of the loading and the stress relaxation behaviour at ambient temperature. Received 30 October 2001; accepted for publication 21 January 2002     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

On the evolution equation for the rest stress in rate-dependent plasticity

Three different structures of the evolution equation for the rest stress are derived. They correspond to different models of viscoplastic response, constructed by three-dimensional generalization of simple one-dimensional rheological models. Isotropic, kinematic and combined isotropic–kinematic hardening, with an evolution equation for the back stress, are incorporated in the constitutive framework. The results may be of interest in the analytical and numerical studies of the rate-dependent inelastic response.     

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.     

Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

Uniaxial and non-proportionally multiaxial ratcheting of SS304 stainless steel at room temperature: experiments and simulations

The uniaxial and non-proportionally multiaxial ratcheting behaviors of SS304 stainless steel at room temperature were initially researched by experiment and then were theoretically described by a cyclic constitutive model in the framework of unified visco-plasticity. The effects of cyclic stress amplitude, mean stress, and their histories on the ratcheting were experimentally investigated under uniaxial and different multiaxial loading paths. The shapes of non-proportional loading paths were linear, circular, elliptical and rhombic, respectively. In the constitutive model, the rate-dependent behavior of the material was reflected by a viscous term; the cyclic flow and cyclic hardening behaviors of the material under asymmetrical stress-controlled cycling were reflected by the evolution rules of kinematic hardening back stress and isotropic deforming resistance, respectively. The effect of loading history on the ratcheting was also considered by introducing two fading memorization functions for maximum inelastic strain amplitude and isotropic deformation resistance, respectively, into the constitutive model. The effect of multiaxial loading path on the ratcheting was reflected by a non-proportional factor defined in this work. The predicting ability of the developed model was proved to be good by comparing the simulations with corresponding experiments.     

Material aspects of dynamic neck retardation

Neck retardation in stretching of ductile materials is promoted by strain hardening, strain-rate hardening and inertia. Retardation is usually beneficial because necking is often the precursor to ductile failure. The interaction of material behavior and inertia in necking retardation is complicated, in part, because necking is highly nonlinear but also because the mathematical character of the response changes in a fundamental way from rate-independent necking to rate-dependent necking, whether due to material constitutive behavior or to inertia. For rate-dependent behavior, neck development requires the introduction of an imperfection, and the rate of neck growth in the early stages is closely tied to the imperfection amplitude. When inertia is important, multiple necks form. In contrast, for rate-independent materials deformed quasi-statically, single necks are preferred and they can emerge in an imperfection-free specimen as a bifurcation at a critical strain. In this paper, the interaction of material properties and inertia in determining neck retardation is unraveled using a variety of analysis methods for thin sheets and plates undergoing plane strain extension. Dimensionless parameters are identified, as are the regimes in which they play an important role.