Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Effective pressure-sensitive elastoplastic behavior of particle-reinforced composites and porous media under isotropic loading

The composite under investigation consists of an elastoplastic matrix reinforced by elastic particles or weakened by pores. The material forming the matrix is pressure-sensitive. The Drucker–Prager yield criterion and a one-parameter non-associated flow rule are employed to formulate the yield behavior of the matrix. The objective of this work is to estimate the effective elastoplastic behavior of the composite under isotropic tensile and compressive loadings. To achieve this objective, the composite sphere assemblage model of Hashin [Z. Hashin, The elastic moduli of heterogeneous materials, ASME J. Appl. Mech. 29 (1962) 143–150] is used. Exact solutions are thus derived as estimations for the effective secant and tangent bulk moduli of the composite. The effects of the loading modes and phase properties on the effective elastoplastic behavior of the composite are analytically and numerically evaluated.     

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.     

Elastic–plastic FE analysis of a notched shaft under multiaxial nonproportional synchronous cyclic loading

An elastic–plastic finite element analysis is presented for a notched shaft subjected to multiaxial nonproportional synchronous cyclic tension/torsion loading. The elastic–plastic material property is described by the von Mises yield criterion and the kinematic hardening rule of Prager/Ziegler. The finite element program system ABAQUS is used to solve the boundary value problem. Special emphasis is given to explore the effects of the stress amplitude, the mean-stress, and the mutual interactions on the local stress–strain responses at the notch root.     

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.     

Reconcile the perfectly elastoplastic model to simulate the cyclic behavior and ratcheting

Perfectly elastoplastic constitutive model is modified through a smoothing factor introduced by Liu [Liu, C.-S., 2003. Smoothing elastoplastic stress–strain curves obtained by a critical modification of conventional models. Int. J. Solids Struct. 40, 2121–2145]. The new model allows plasticity to happen in a non-zero-measure yield volume in stress space, rather than that of conventional zero-measure yield surface, and within the yield volume the plastic modulus is varying continuously. It endows a specific strain-hardening rule of flow stress and is able to describe the phenomena of strain hardening, cyclic hardening, the Bauschinger effect, mean-stress relaxation, strain ratcheting, out-of-phase hardening, as well as erasure-of-memory. In order to suppress the over prediction of ratcheting we consider a scalar function of smoothing factor, which can simulate the saturation behavior of uniaxial/multiaxial strain ratcheting. These effects are demonstrated through numerical examples. The existence of stress equilibrium point and limiting surface is a natural result without requiring an extra design. Moreover, the non-linear constitutive equations can be converted into a linear system for augmented stress in the Minkowski space, of which the symmetry group is a proper orthochronous Lorentz group SO_o(5, 1). The augmented stress is a time-like vector, moving on hyperboloids inside the cone. When taking the Prager kinematic hardening rule into account we can simulate some cyclic behaviors of SAE 4340 and grade 60 steels within a certain accuracy through the use of only three material constants and a fixed smoothing factor. To simulate the ratcheting behaviors of SS304 stainless steel we allow the smoothing factor to be an exponential decaying function of λ.     

Experimental investigation of cyclic and time-dependent deformation of titanium alloy at elevated temperature

A novel cyclic deformation test program was undertaken to characterize macroscopic time dependent deformation of a titanium alloy for use in viscoplastic model development. All tests were conducted at a high homologous temperature, 650 °C, where there are large time dependent and loading rate dependent effects. Uninterrupted constant amplitude tests having zero mean stress or a tensile mean stress were conducted using three different control modes: strain amplitude and strain rate, stress amplitude and stress rate, and a hybrid stress amplitude and strain rate. Strain ratcheting occurred for all cyclic tests having a tensile mean stress and no plastic shakedown was observed. The shape of the strain ratcheting curve as a function of time is analogous to a creep curve having primary, steady state and tertiary regions, but the magnitude of the ratchet strains are higher than creep strains would be for a constant stress equal to the mean stress. Strain cycles interrupted with up to eight 2-h stress relaxation periods around the hysteresis loop, including hold times in each quadrant of the stress–strain diagram, were also conducted. Stress relaxation was path-dependent and in some cases the stress relaxed to zero. The cyclic behavior of these interrupted tests was similar even though each cycle was very complex. These results support constitutive model development by providing exploratory, characterization and validation data.     

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.     

Fatigue life predictions by integrating EVICD fatigue damage model and an advanced cyclic plasticity theory

An event independent cumulative damage (EVICD) fatigue prediction model was previously developed for the fatigue damage prediction under general multiaxial stress state and loading conditions. The model takes the plastic strain energy as the major contributor to the fatigue damage. The application of the EVICD model does not require a cycle counting method for general random loading. In the current effort, derivations were made to explicitly and directly relate the material constants in the fatigue model to the parameters in the Manson–Coffin equations and the cyclic stress–strain curve of the material. In addition, an advanced cyclic plasticity theory was implemented for the determination of the detailed stress–strain response that was required as the input for the EVICD fatigue model. Three metallic materials were used to demonstrate the capability of the modified fatigue model for the predictions of fatigue lives under different loading conditions. The results show that the fatigue model can provide fatigue life predictions in close agreement with the experimental observations.     

Texture development and plastic anisotropy of B.C.C. strain hardening sheet metals

A Taylor-like polycrystal model is adopted here to investigate the plastic behavior of body centered cubic (b.c.c.) sheet metals under plane-strain compression and the subsequent in-plane biaxial stretching conditions. The <111> pencil glide system is chosen for the slip mechanism for b.c.c. sheet metals. The {110} <111> and {112} <111> slip systems are also considered. Plane-strain compression is used to simulate the cold rolling processes of a low-carbon steel sheet. Based on the polycrystal model, pole figures for the sheet metal after plane-strain compression are obtained and compared with the corresponding experimental results. Also, the simulated plane-strain stress—strain relations are compared with the corresponding experimental results. For the sheet metal subjected to the subsequent in-plane biaxial stretching and shear, plastic potential surfaces are determined at a given small amount of plastic work. With the assumption of the equivalence of the plastic potential and the yield function with normality flow, the yield surfaces based on the simulations for the sheet metal are compared with those based on several phenomenological planar anisotropic yield criteria. The effects of the slip system and the magnitude of plastic work on the shape and size of the yield surfaces are shown. The plastic anisotropy of the sheet metal is investigated in terms of the uniaxial yield stresses in different planar orientations and the corresponding values of the anisotropy parameter R, defined as the ratio of the width plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial tensile loading. The uniaxial yield stresses and the values of R at different planar orientations from the polycrystal model can be fitted well by a yield function recently proposed by Barlat et al. (1997b).     

Evaluation of identification methods for YLD2004-18p

Four calibration methods have been evaluated for the linear transformation-based anisotropic yield function YLD2004-18p (Barlat, F., Aretz, H., Yoon, J.W., Karabin, M.E., Brem, J.C., Dick, R.E., 2005. Linear transformation-based anisotropic yield functions. Int. J. Plasticity 21, 1009–1039) and the aluminium alloy AA5083-H116. The different parameter identifications are based on least squares fits to combinations of uniaxial tensile tests in seven directions with respect to the rolling direction, compression (upsetting) tests in the normal direction and stress states found using the full-constraint (FC) Taylor model for 690 evenly distributed strain paths. An elastic–plastic constitutive model based on YLD2004-18p has been implemented in a non-linear finite element code and used in finite element simulations of plane-strain tension tests, shear tests and upsetting tests. The experimental results as well as the Taylor model predictions can be satisfactorily reproduced by the considered yield function. However, the lacking ability of the Taylor model to quantitatively reproduce the experiments calls for more advanced crystal plasticity models.     

考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

Mechanism-based strain gradient Drucker–Prager elastoplasticity for pressure-sensitive materials

Modeling strain gradient plasticity effects has achieved considerable success in recent years. However, incorporating the full mechanisms of the pressure-sensitive yielding and the size dependence of plastic deformation still remains an open challenge. In this work, a mechanism-based stain gradient (MSG) plasticity theory for pressure sensitive materials with a variable material length-scale parameter is presented. The flow theory of MSG plasticity based on the Drucker–Prager yield function is established following the same hierarchical framework of MSG plasticity proposed by Gao et al., 1999, Huang et al., 2000 and Qiu et al. (2003) in order to link the strain gradient plasticity theory on the mesoscale to the Taylor dislocation model on the micro-scale. The incremental constitutive relation based on the associated flow rule is derived for the Drucker–Prager yield function on the micro-scale, including the higher-order stress introduced as the thermodynamic conjugate of strain gradient at the mesoscale. The proposed theory successfully predicts the experimental results. The numerical results show that when the pressure-sensitivity index defined by the Drucker–Prager yield function takes different values, the material response curves are different and the material strength increases with the increase of pressure-sensitivity index. It proves that this procedure is able to represent the material behavior of pressure-sensitive materials such as geomaterials, polymeric materials, metallic foams and metallic glass.     

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.     

Evaluation of cyclic plasticity models in ratcheting simulation of straight pipes under cyclic bending and steady internal pressure

This paper evaluates seven cyclic plasticity models for structural ratcheting response simulations. The models evaluated are bilinear (Prager), multilinear (Besseling), Chaboche, Ohno–Wang, Abdel Karim–Ohno, modified Chaboche (Bari and Hassan) and modified Ohno–Wang (Chen and Jiao). The first three models are already available in the ANSYS finite element package, whereas the last four were implemented into ANSYS for this study. Experimental responses of straight steel pipes under cyclic bending with symmetric end rotation history and steady internal pressure were recorded for the model evaluation study. It is demonstrated that when the model parameters are determined from the material response data, none of the models evaluated perform satisfactorily in simulating the straight pipe diameter change and circumferential strain ratcheting responses. A detailed parameter sensitivity study with the modified Chaboche model was conducted to identify the parameters that influence the ratcheting simulations and to determine the ranges of the parameter values over which a genetic algorithm can search for refinement of these values. The refined parameter values improved the simulations of straight pipe ratcheting responses, but the simulations still are not acceptable. Further, improvement in cyclic plasticity modeling and incorporation of structural features, like residual stresses and anisotropy of materials in the analysis will be essential for advancement of low-cycle fatigue response simulations of structures.     

非比例循环塑性内时本构关系研究

在Valanis的内时本构理论的框架中,引入内结构张量以反映由于非比例加载而引起金属材料的附加等向强化及异向强化效应,同时提出材料强化程度的度量采用沿路径法线方向的塑性应变分量来描述.这些考虑的有效性已经通过用所建模型对304不锈钢材料在一些典型非比例循环加载路径下的响应进行的理论预测得到了验证;将该模型应用于U71Mn材料室温单轴棘轮行为描述中,结果显示内结构张量的引入不仅能较好地反映应变控制下的非比例附加效应,而且也能较好地反映应力控制下塑性应变的累积及变化率.     

Orthotropic yield criteria for description of the anisotropy in tension and compression of sheet metals

In this paper, yield functions describing the anisotropic behavior of textured metals are proposed. These yield functions are extensions to orthotropy of the isotropic yield function proposed by Cazacu et al. (Cazacu, O., Plunkett, B., Barlat, F., 2006. Orthotropic yield criterion for hexagonal close packed metals. Int. J. Plasticity 22, 1171–1194). Anisotropy is introduced using linear transformations of the stress deviator. It is shown that the proposed anisotropic yield functions represent with great accuracy both the tensile and compressive anisotropy in yield stresses and r-values of materials with hcp crystal structure and of metal sheets with cubic crystal structure. Furthermore, it is demonstrated that the proposed formulations can describe very accurately the anisotropic behavior of metal sheets whose tensile and compressive stresses are equal.     

Multiscale experimental investigations about the cyclic behavior of the 304L SS

This work follows a series of experiments carried out earlier at INSA of Rouen (Hassan, T., Taleb, L., Krishna, S., 2008. Influence of non-proportional loading paths on ratcheting responses and simulations by two recent cyclic plasticity models. Int. J. Plast. 24, 1863–1889). It investigates the elastoplastic cyclic behavior of a 304L stainless steel at room temperature. In a first step the cross path effect on ratcheting is confirmed, as well as the crucial role of the loading path non-proportionality. Strain controlled tests are also conducted for different strain amplitudes and loading paths. Cross-hardening effect appears more important when the shearing sequence is followed by the axial one. Moreover for alternating axial and shearing cycles, this phenomenon occurs after each crossing sequence leading to a very significant strain hardening, at least of the same order as the one obtained after a circular strain path. Yet, the magnitude of the observed over hardening does not necessarily seem a function of the cumulated plastic strain.