Effect of martensite volume fraction on cyclic plastic deformation behavior of dual phase steel: micromechanics simulation study

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

Available online at www.sciencedirect.com

www.jmrt.com.br

Original Article

Effect of martensite volume fraction on cyclic plastic deformation behavior of dual phase steel: micromechanics simulation study Amit Kumar Rana a,∗ , Surajit Kumar Paul b , Partha Pratim Dey a a b

Mechanical Engineering Department, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India Department of Mechanical Engineering, Indian Institute of Technology Patna, Bihar 801106, India

a r t i c l e

i n f o

a b s t r a c t

Article history:

The cyclic plastic deformation of an archetypal microstructure dual phase steels has been

Received 1 January 2019

examined via micromechanics modelling based on representative volume elements tech-

Accepted 14 June 2019

nique. The dual phase steel has soft ferrite – matrix with hard martensite – islands, which

Available online 6 July 2019

are distributed in a discrete manner. In the field of automobile industries, the suitable combination of strain hardening, strength and ductility, and their lean combination represent

Keywords:

them as an economically desire option for huge multiple lightweight options. In this numer-

Cyclic plastic deformation

ical approach, different microstructures for micromechanical modelling were constructed

Dual phase steel

to detailed examination and analysis for the tensile and cyclic deformation response of

Martensite volume fraction

dual phase steels. Due to the distinct difference in the stress–strain responses of ferrite and

Plastic strain localization

martensite phases, incompatibility of strain between matrix-softer ferrite and island-harder

Stress triaxiality

martensite phase arises during tensile straining. The effect of strain partitioning in cyclic

Material inhomogeneity

plastic deformation response of dual phase steel has been studied in the present investigation. Apart from this, effect of martensite volume fraction on cyclic stress–strain response of dual phase steels, cyclic plastic deformation distribution in individual phases and cyclic deformation inhomogeneity at microstructural level have been systematically investigated. © 2019 Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1.

Introduction

Ferrite–martensite dual phase (DP) steel have good combination of mechanical properties and which make it suitable to use in engineering structural application such as automotive body section, wheels disc, reinforcement beams in inner bumper etc. Such engineering automobile components are inherently exposed to assorted cyclic loading conditions

∗

during the course of driving and that may cause cyclic plastic deformation of the material. Hence in such structural mechanisms of DP steel, it is significant to recognize the fatigue cracks places. Fatigue cracks can be generally found in places where stress/deformation concentration occurs. Normally deformation concentration takes place because of geometrical variation or alteration of material properties. Severe research has been concluded to comprehend the impact of the microstructural changes on the tensile and fatigue property of DP steels. Present authors Rana et al. [1] reported that the martensite volume fraction has direct important worthy attention effect on partitioning of strain apportioning between

Corresponding author. E-mail: [email protected] (A.K. Rana). https://doi.org/10.1016/j.jmrt.2019.06.022 2238-7854/© 2019 Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

3706

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

soft ferrite and hard martensite phases and hence tensile properties of DP steels. Yang and Chen [2] and Ramazani et al. [3] considered the impact of martensite volume part and martensite morphology on the mechanical collective possessions and failure behaviour of DP steels respectively. Paul et al. [4,5] revealed that the LCF life for apt plastic strain amplitude is fundamentally lessened by an expansion in the martensite volume division in the microstructure [4] and an increase in coarseness in martensite morphology [5]. As reported by Mediratta et al. [6], substantially LCF behaviour of DP steel is affected by the martensite morphology with a martensite volume fraction of 21%. Comparatively superior LCF life may be found in DP steel that has a fine grained matrix-soft ferrite and consistent emplacements of martensite, while considerably inferior LCF life can be seen in DP steel where the ferrite matrix has course martensite particles. All distributed literature on LCF conduct of DP steel reasoned that there is a slight variation in the specific behaviour of LCF life due to change in morphology of martensites. From the review of literature it has been observed that, although researchers unanimously mentioned that the martensite volume percentage of phases has an important effect on LCF life, but the systematic numerical detail investigation has not been done yet. Sherman and Davies [7] demonstrated that work hardening, which after-effects of flow-incited alter in the microstructure, has a considerable supremacy on behavior of fatigue of the dual phase steels. In addition to experimental findings, extensive computational modelling has been vigorously sought after to demonstrate the influence of microstructure on mechanical properties. Sun et al. [8,9] industrialized a microstructurebased modeling technique in which the failure mode and ultimate ductility of dual phase steels are anticipated under various quasi-static tensile loading settings utilizing the local plastic strain restriction hypothesis. For DP980 and DP780 steels Sun et al. [8], for DP 590 and DP 780 steels Paul et al. [10–12], and for DP steel Marvi-Mashhadi et al. [13] demonstrated that inhomogeneous microstructures in DP 980, DP 780 and DP 590 steels have been shown to act as imperfections that initiates the plastic instability causing plastic strain localization during deformation. Sodjit and Uthaisangsuk [14] found that a few short intruded on shear groups created with comparatively low fraction of martensite volume, they also observed that the volume part of martensite increments with expanding continuous pronounced long localize shear band, which is appeared in the DP microstructure prominently. Paul [15] showed that for DP steel ultimate tensile strength increases whereas percentage of elongation decreases with increase in martensite volume fraction. But all investigations are limited for tensile loading condition only. Therefore, computational modeling to study the mechanical behavior of DP steels for cyclic loading is justified.

2. Micromechanical simulation of microstructure with a 2D RVE Dual Phase steels (DP590 and DP780) with their tensile strength and constitutive behavior are used for micromechanical RVE FE simulation in the present work. DP steels macroscopic perceptible mechanical properties are not only

just rely on the stress–strain response of ferrite yet additionally volume portion and morphology of martensite [15–18]. Distributed works depicted that damage originated in DP steels amid tensile loading by the particle cracking of island (martensite) or decohesion of matrix — island interface, trailed by void nucleation and development in the matrix (ferrite) [19–25]. The technique in light of damage mechanics perception can’t be allied without considerable further advancement as the mechanisms of damage evolution are complex in nature of DP steel. Micromechanics representation using Representative Volume Element (RVE) was adopted by number of research groups [8,10–15,26,27] to model tensile deformation response of DP steels in microstructural level. The material input quantities supplied to the RVE are the individual flow response of ferrite and martensite. Detailed interaction between interfaces boundaries of the phases were neglected as it is insignificant (atomic sizes) compared to the RVE. Slight change in numerical value may be observed without noticeable change in trend. The results can be obtained from the micromechanics model without incorporating any separate damage law or failure criteria. Similar two dimensional RVE approach is adopted for the present work. The cyclic deformation responses of DP steels were modelled by constructing four RVEs corresponding to martensite volume fractions of 13%, 23%, 33% and 43% (Fig. 1). The RVEs size for the present investigation is 1000 × 1000 ␮m and the average area of a martensite island is approximately 5 ␮m with aspect ratio of 2. A finite element code Abaqus is used in the present research. To mesh the RVEs, two dimensional plane strain elements are used. To simulate uniaxial tensilecompressive loading condition, all of the nodes along the bottom edge are compelled in the y bearing path and remain free in the x course. In RVEs loading, all of the nodes along the top length edge are given the similar relocations in the y course, whereas they stay free in the x bearing. The engineering strain in the y direction is attained by dividing the displacements of the top length edge by original length of the FE RVE. The engineering stress in the y direction is obtained from load divided by the primary area of cross section.

3.

LCF experimentation

Using cognition from numerical (RVE) analysis and also for further experimentation, 75 mm × 150 mm sheet of 1 mm thickness was used for preparing LCF specimens for DP590 and DP780. The specimen length is considered in the rolling direction. For the fatigue test sample, a 2 mm width gauge and 7.9 mm parallel length was used, as done previously by Paul et al. [4,5]. As per ASTM E606-92 [28], total strain amplitude control LCF tests were conducted under fully reversed (R = −1.0) condition. Since large compressive strains tend to produce buckling, hence anti-buckling guides were used. As per the ASTM, E8 standard tensile test sample has been prepared considering gauge length and width 50 mm and 12.5 mm, respectively. Moreover, it was prepared by machining steel sheets along rolling direction. At room temperature, tests were conducted on the specimen. It was performed using an electro-mechanical testing machine at a crosshead speed of 1 mm min−1 corresponding to a strain rate of 0.333 × 10−3 s−1 .

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

3707

Fig. 1 – Representative volume elements (RVEs) of the intended entire microstructures with various percentage of ferrite and martensite (a) 87% ferrite and 13% martensite (b) 77% ferrite and 23% martensite (c) 67% ferrite and 33% martensite (d) 57% ferrite and 43% martensite.

4. Numerical analysis of kinematic hardening parameter Rana et al. [1] carried out RVE simulation of DP (ferrite–martensite) steels having major applied range of

martensite volume fractions in used of 13%, 23%, 33% and 43% to determine stress–strain response under uniaxial tensile loading. These four grades of DP steel were supplied with individual stress–strain response of the two phases to determine the macroscopic stress–strain relationship of most commonly used range of automotive grade steel. The same

3708

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

Fig. 2 – Monotonic validation of stress–strain relationship of experimental and simulation.

stress–strain relationships of the phases are employed in this study to understand the cyclic deformation characteristic of DP steels. The monotonic simulated stress–strain curve shown in Fig. 2 by 2D RVE model is matching well with the experimental (tensile test) result. By contriving different heat treatment processes and anneals, then it can be generated DP steels with different martensite fractions [4,5]. The tensile stress–strain relationships based on dislocation theory as per Paul and Kumar [10] of the each phases i.e. isle and matrix are employed in this study to understand the cyclic deformation characteristic of DP steels cognitive to Ohno and Wang [29], Abdel-Karim and Ohno [30], and Zhang and Jiang [31]. To investigate the cyclic plastic deformation response of DP (ferrite–martensite) steel, kinematic hardening rule introduced by Chaboche is incorporated means of superposition of several kinematic hardening models to describe the cyclic stress–strain response of matrix and island phases singly in the finite element micromechanics RVE. The commonly accepted equation form of Chaboche kinematic hardening law [32] would be:

d˛ =

3 

d˛j

(1)

j=1

d˛¯ j =

2 p C d¯εp − j ˛¯ j dεeq 3 j

(2)

The material constraints in the above equations are C1 , C2 , C3 ,  1 ,  2 and  3 , where C1 , C2 and C3 are the coefficients of kinematic hardening, and  1 ,  2 and  3 are the exponents of p ¯ d¯εp and dεeq are the increkinematic hardening. Whereas, d˛, ment of back stress vector, increment vector of plastic strain and equivalent plastic strain increment defined respectively in the Chaboche numerical model. In the above law, the result back stress (˛) is the addition of three (1–3) decomposed back stresses (˛ = ˛1 + ˛2 + ˛3 ). The foremost one (˛1 ) is projected to forecast a high plasticity modulus at the onset of yielding point and it stabilizes very rapidly. The next one (˛2 ) is projected to simulate the transient nonlinear segment of the hysteresis loop. The last decomposed back stress (˛3 ) would simulate

Fig. 3 – Prediction of tensile stress–strain curve of phases of DP590 steel by Chaboche model.

the linear hardening (i.e. constant plasticity modulus) in the higher strain range. The equation can be represented for the loading part of the flow curve as 3 

˛j + 0 = 

(3)

j=1

where  0 is the yield stress of cyclic flow. Several earlier studies [32–34] have been made about the details material constant determination procedure. In Fig. 3, it has been shown that the Chaboche model was able to successfully predict monotonic stress–strain curve of ferrite and martensite phases. Stress vs. strain hysteresis loop of both ferrite and martensite for microstructure predicted by Chaboche model for result strain amplitude of 1.0% are shown in Fig. 4.

5.

Results and discussions

The total strain in loading direction distributed in DP steel with 23% martensite fraction for both tensile and compressive peaks are illustrated in Fig. 5. During cyclic plastic deformation, there is inhomogeneous deformation of hard and soft phase. This occurs because the plastic deformation of the softer phase (ferrite matrix) is opposed by the harder phase (martensite particles). Fig. 6 examines the stress triaxiality build up when loading pattern is tension followed by compression for strain amplitude of 1.0% in same DP steel. Nature of stress triaxiality development in tensile and compressive peaks is same while the sign changes. Stress triaxiality builds up in microstructural level of DP steel occurs as the hard martensite particles restrict the plastic deformation of soft ferrite matrix locally and alters the local stress state. Fig. 7 shows plastic strain amplitude in loading direction during low cyclic fatigue for RVE with 13%, 23%, 33% and 43% of martensite volume fractions having strain amplitude of 1.0%. Fig. 7 shows inhomogeneous distribution of plastic strain amplitude in RVE of different grade DP steels with various martensite volume fractions. The plastic strain

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

3709

Fig. 4 – Stress vs. strain hysteresis loop prediction of ±1.0% strain amplitude: (a) for ferrite phase and (b) for martensite phase by Chaboche model.

Fig. 5 – Total strain in loading direction (LE22) during LCF with strain amplitude of 1.0% for RVE with 77% ferrite and 23% martensite (a) tensile peak, and (b) compression peak.

amplitude in each integration points in ferrite phase of different DP steel RVEs having significant grades range martensite volume fractions are computed and plotted in Fig. 8. It can be made clear from Fig. 8 that the distribution become wider and the peak also become lower as the martensite volume fraction gradually ascends from 13 to 43%. The wider distribution indicates higher plastic strain amplitude in more integration points in ferrite phase i.e. higher plastic deformation in some areas of ferrite phase. During straining process,

according to various research [11,15,35,36], DP steels show strain partitioning between the two phases. For a given applied strain level, the deformation in ferrite is far more than the other phase, obviously the extent of deformation in martensite phase depends upon the volume fraction and carbon content in martensite phase. As the martensite volume fraction ascends, it is logical that for a given applied strain level, the deformation in ferrite phase will increase. This explains the wider distribution of plastic strain amplitude

Fig. 6 – Stress triaxiality during LCF with strain amplitude of 1.0% for RVE with 77% ferrite and 23% martensite (a) tensile peak, and (b) compression peak.

3710

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

Fig. 7 – Plastic strain amplitude in loading direction (PEamp) during LCF with strain amplitude of 1.0% for RVE with (a) 87% ferrite + 13% martensite (b) 77% ferrite + 23% martensite (c) 67% ferrite + 33% martensite (d) 57% ferrite + 43%martensite.

in ferrite phase with increasing volume fraction of martensite. Stress–strain hysteresis loop during LCF with total strain amplitude of 1.0% for RVE with 13–43% volume fraction of island phase is plotted in Fig. 9. The strain hardening rate

in the plastic region of the stress–strain hysteresis loop shows an upward trend as the martensite volume fraction ascends. This can be explained in two combine phenomena: stress partitioning and deformation constrain in ferrite by martensite island. Stress partitioning between the dual phases

Fig. 8 – Distribution of plastic strain amplitude in loading direction during LCF with strain amplitude of 1.0% for RVEs with considered ascending martensite volume of 13–43%.

Fig. 9 – Stress–strain hysteresis loop during LCF with total strain amplitude of 1.0% for RVE with significant range of volume fraction of martensite.

3711

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

560

0.78

540

0.76

520

0.74

500

0.72 Total strain amplitude = 1.0%

480

Plastic strain amplitude, %

Stress amplitude, MPa

Stress amplitude, MPa Plastic strain amplitude, %

0.70 10

20

30

40

Martensite voume, %

Fig. 11 – Variation of stress amplitude and plastic strain amplitude with volume fraction of martensite of DP steels for LCF with total strain amplitude of 1.0%.

is maintained constant and as stress amplitude increases with increasing volume fraction of harder phase (martensite), therefore plastic strain amplitude decreases when the volume fraction of harder phase (martensite) increases.

6.

Fig. 10 – Prediction of flow hysteresis loops during LCF at the total strain amplitude for 1.0%: (a) DP 590 steel, (b) DP 780 steel.

during plastic deformation consequences higher stress level in martensite phase in comparison to ferrite phase. Higher martensite volume fraction contributes higher microstructural volume faction of DP steel carries high stress level and hence high strain hardening behaviour. Apart from this, it is well established that martensite island restrict deformation in ferrite phase of DP steel and which results higher strain hardening behaviour. As the volume percentage of martensite ascends, the restriction of deformation in ferrite phase also increases and as a consequence strain hardening behaviour also increases in DP steel. The stress–strain hysteresis loop during LCF with total strain amplitude of 1.0% for distinct DP590 and DP780 steels are plotted in Fig. 10(a) and (b), respectively. It follows a prediction trend obtained from simulation. Variation of stress/plastic strain amplitude of two phase steel with volume fraction of martensite having total strain amplitude of 1.0% during LCF is shown in Fig. 11. Stress amplitude shows enhancement where as the plastic strain amplitude shows an opposite behavior as the volume fraction of martensite ascends in DP steels. Stress partitioning between dual phases and deformation constrain in soft matrix by martensite island are the cause of increase in stress amplitude with ascending martensite. The strain amplitude as a whole

Conclusions

Cyclic plastic deformation responses of dual phase steels with different martensite volume fraction were investigated in the present work with microstructure based micromechanical modelling technique. In this research, the main investigated conclusions are summarized as follows: • All DP steels without reference to ascending martensite volume fraction, the cyclic plastic deformation in microstructural level is inhomogeneous. Extent of inhomogeneous cyclic plastic deformation increases gradually from 13 to 43% of martensite volume fraction in the microstructure. Distributions of plastic strain amplitude in ferrite phase widen and peak becomes lower with increasing martensite volume fraction. • Due to local deformation restriction in the ferrite phase, the local stress triaxiality builds up. Similar type non uniform local stress triaxiality build up is noticed in ferrite phase during tensile and compressive peak. • Shape of the stress vs. strain hysteresis loop of ferrite–martensite dual phase steel alters with ascending martensite volume. The hardening rate becomes higher in the plastic region, stress amplitude enhances and plastic strain amplitude lessens with increasing martensite volume fraction. • The RVE simulated and experimental hysteresis flow loop shows close concurrence. Hence, micromechanics investigation with trivialized experiment is justified and can be considered as a method to predict the cyclic deformation behavior of dual phase steel with volume fractions. However more rigorous investigation is required before engineering application.

3712

j m a t e r r e s t e c h n o l . 2 0 1 9;8(5):3705–3712

Conflicts of interest [18]

The authors declare no conflicts of interest.

references

[19]

[20] [1] Rana AK, Paul SK, Dey PP. Effect of martensite volume faction on strain partitioning behavior of dual phase steel. Phys Mesomech 2018;21(4):333–40. [2] Yang JR, Chen LJ. Dual ferrite-martensite treatments of a high-strength low-alloy ASTM A588 steel. J Mater Sci 1991;26(4):889–98. [3] Ramazani A, Ebrahimi Z, Prahl U. Study the effect of martensite banding on the failure initiation in dual-phase steel. Comp Mater Sci 2014;87:241–7. [4] Paul SK, Stanford N, Hilditch T. Effect of martensite volume fraction on low cycle fatigue behaviour of dual phase steels: experimental and microstructural investigation. Mater Sci Eng A 2015;638:296–304. [5] Paul SK, Stanford N, Hilditch T. Effect of martensite morphology on low cycle fatigue behaviour of dual phase steels: experimental and microstructural investigation. Mater Sci Eng A 2015;644:53–60. [6] Mediratta SR, Ramaswamy V, Rao PR. Influence of ferrite-martensite microstructural morphology on the low cycle fatigue of a dual-phase steel. Int J Fatigue 1985;7(2):107–15. [7] Sherman AM, Davies RG. Fatigue of a dual-phase steel. Metall Trans A 1979;10(7):929–33. [8] Sun X, Choi KS, Liu WN, Khaleel MA. Predicting failure modes and ductility of dual phase steels using plastic strain localization. Int J Plasticity 2009;25(10):1888–909. [9] Sun X, Choi KS, Soulami A, Liu WN, Khaleel MA. On key factors influencing ductile fractures of dual phase (DP) steels. Mater Sci Eng A 2009;526(1):140–9. [10] Paul SK, Kumar A. Micromechanics based modeling to predict flow behavior and plastic strain localization of dual phase steels. Comput Mater Sci 2012;63:66–74. [11] Paul SK, Mukherjee M. Determination of bulk flow properties of a material from the flow properties of its constituent phases. Comp Mater Sci 2014;84:1–12. [12] Paul SK. Real microstructure based micromechanical model to simulate microstructural level deformation behavior and failure initiation in DP 590 steel. Mater Des 2013;44:397–406. [13] Marvi-Mashhadi M, Mazinani M, Rezaee-Bazzaz A. FEM modeling of the flow curves and failure modes of dual phase steels with different martensite volume fractions using actual microstructure as the representative volume. Comput Mater Sci 2012;65:197–202. [14] Sodjit S, Uthaisangsuk V. Microstructure based prediction of strain hardening behavior of dual phase steels. Mater Des 2012;41:370–9. [15] Paul SK. Effect of martensite volume fraction on stress triaxiality and deformation behavior of dual phase steel. Mater Des 2013;50:782–9. [16] Bag A, Ray KK, Dwarakadasa ES. Influence of martensite content and morphology on tensile and impact properties of high-martensite dual-phase steels. Metall Trans A 1999;30(5):1193–202. [17] Erdogan M, Tekeli S. The effect of martensite particle size on tensile fracture of surface-carburised AISI 8620 steel with

[21]

[22]

[23]

[24]

[25]

[26]

[27]

[28]

[29]

[30]

[31]

[32]

[33] [34]

[35]

[36]

dual phase core microstructure. Mater Des 2002;23(7):597–604. Pierman A-P, Bouaziz O, Pardoen T, Jacques PJ, Brassart L. The influence of microstructure and composition on the plastic behaviour of dual-phase steels. Acta Mater 2014;73:298–311. Uthaisangsuk V, Prahl U, Bleck W. Modelling of damage and failure in multiphase high strength DP and TRIP steels. Eng Fract Mech 2011;78(3):469–86. Gou R, Dan W, Zhang W, Yu M, Zhang C, Qiao Y, et al. Microstructure failure in ferrite-martensite dual phase steel under in-situ tensile test. J Iron Steel Res Int 2017;24(3):350–6. Erdogan M. The effect of new ferrite content on the tensile fracture behaviour of dual phase steels. J Mater Sci 2002;37(17):3623–30. Zhang J, Di H, Deng Y, Misra RDK. Effect of martensite morphology and volume fraction on strain hardening and fracture behavior of martensite–ferrite dual phase steel. Mater Sci Eng A 2015;627:230–40. Kadkhodapour J, Butz A, Ziaei-Rad S, Schmauder S. A micro mechanical study on failure initiation of dual phase steels under tension using single crystal plasticity model. Int J Plasticity 2011;27(7):1103–25. Maire E, Bouaziz O, Michiel MD, Verdu C. Initiation and growth of damage in a dual-phase steel observed by X-ray microtomography. Acta Mater 2008;56(18):4954–64. Kim DK, Kim EY, Han J, Woo W, Choi SH. Effect of microstructural factors on void formation by ferrite/martensite interface decohesion in DP980 steel under uniaxial tension. Int J Plasticity 2017;94:3–23. Uthaisangsuk V, Prahl U, Bleck W. Stretch-flangeability characterisation of multiphase steel using a microstructure based failure modeling. Comp Mater Sci 2009;45:617–23. Uthaisangsuk V, Prahl U, Bleck W. Micromechanical modelling of damage behaviour of multiphase steels. Comp Mater Sci 2008;43(1):27–35. ASTM Standard E466. Conducting force controlled constant amplitude axial fatigue tests of metallic materials. West Conshohocken, PA: ASTM International; 2007. Ohno N, Wang J-D. Kinematic hardening rules with critical state of dynamic recovery, part I: formulation and basic features for ratchetting behavior. Int J Plast 1993;9(3):375–90. Abdel-Karim M, Ohno N. Kinematic hardening model suitable for ratchetting with steady-state. Int J Plasty 2000;16(3-4):225–40. Zhang J, Jiang Y. Constitutive modeling of cyclic plasticity deformation of a pure polycrystalline copper. Int J Plast 2008;24(10):1890–915. Chaboche JL. On some modifications of kinematic hardening to improve the description of ratchetting effects. Int J Plast 1991;7(7):661–78. Bari S, Hassan T. Anatomy of coupled constitutive models for ratcheting simulation. Int J Plast 2000;16(3-4):381–409. Paul SK, Sivaprasad S, Dhar S, Tarafder M, Tarafder S. Simulation of cyclic plastic deformation response in SA333 C–Mn steel by a kinematic hardening model. Comput Mater Sci 2010;48(3):662–71. Paul SK. Effect of material inhomogeneity on the cyclic plastic deformation behavior at the microstructural level: micromechanics-based modeling of dual-phase steel. Model Simul Mat Sci Eng 2013;21(5):055001–26. Cong ZH, Jia N, Sun X, Ren Y, Almer J, Wang YD. Stress and strain partitioning of ferrite and martensite during deformation. Metall Mater Trans A 2009;40(6):1383–7.