Next Article in Journal
An Overview of Progressive Collapse Behavior of Steel Beam-to-Column Connections
Previous Article in Journal
CPW-Fed Transparent Antenna for Vehicle Communications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cyclic Plasticity Model with Martensite Transformation for S30408 and Its Finite Element Implementation

1
School of Chemical Engineering and Technology, Hebei University of Technology, Tianjin 300130, China
2
School of Control Engineering, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China
3
School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(17), 6002; https://doi.org/10.3390/app10176002
Submission received: 9 July 2020 / Revised: 26 August 2020 / Accepted: 27 August 2020 / Published: 29 August 2020
(This article belongs to the Section Chemical and Molecular Sciences)

Abstract

:
The ability of the constitutive model to simulate the ratcheting behavior of metastable austenitic stainless steel S30408 is significant to ensure the safety of the liquefied natural gas (LNG) semi-trailer tanks in the lightweight process of the inner containers. This is because the lightweight inner vessels often encounter cyclic stresses due to the road inertia loads together with high mean stresses due to internal pressures. In this study, we performed cryogenic uniaxial tension experiments and a series of ratcheting experiments to investigate the cyclic plasticity behavior of the metastable austenitic stainless steel S30408. Based on the Ohno-Wang II model, we proposed a new cyclic plasticity constitutive model with martensitic transformation, which relates the content of deformation-induced martensite with isotropic hardening and kinematic hardening. The ratcheting behaviors of S30408 were first simulated by the proposed model with the incremental loading method using MATLAB. The results showed that the model could reasonably predict the ratcheting behavior of S30408, and the evolution law of martensite content could well predict the content of deformation-induced martensite. Under the assumption of the von Mises yield criterion and normal plasticity flow rule, we developed a numerical algorithm of plastic strain with the proposed model to implement the finite element calculation of the model. Internal iteration in the numerical algorithm was implemented with the Euler backward method, which calculated the trial strain for each equilibrium iteration using the consistent tangent matrix. With a user subroutine, the proposed model was programmed into ANSYS for a user - executable version. By simulating the uniaxial ratcheting of a S30408 round bar, we found that the calculated results were in good agreement with the experimental results, which promises further applications in the design of structures, such as LNG semi-trailer tanks.

1. Introduction

With excellent plasticity and toughness, S30408 austenitic stainless steel (type 304 stainless steel in American Iron and Steel Institute (AISI) materials designation) is often employed in the inner containers of LNG semi-trailer tanks [1]. In practical engineering applications, due to the complex road spectrum load and inertial forces, such as braking and steering, alternating stress is generated in the supporting regions of the inner container. The alternating stress together with the mean stress due to internal pressure becomes sufficiently high to make the material in this area encounter cyclic plasticity and produce gradual plastic strain accumulation, i.e., ratcheting. In the design of practical engineering structures, the cumulative damage caused by the ratcheting strain is an important issue to consider for the structural integrity [2].
Many scholars have produced extensive research on the ratcheting behavior of various materials and established a series of constitutive models. The classic Armstrong and Frederick (A-F) [3] model employed nonlinear characteristics due to the introduction of dynamic recovery of backstress and overpredicted ratcheting. Chaboche [4,5] described ratcheting with several A-F models. Ohno and Wang [6,7] introduced a critical state of dynamic recovery; however, multiaxial ratcheting prediction is yet be successfully demonstrated. Many scholars have modified and developed the constitutive model by associating different parameters [8,9,10,11,12,13].
Cyclic softening and hardening can have a great influence on ratcheting in unsymmetrical stress cycles [14,15]. The cyclic softening and hardening behaviors of materials are typically reflected in the isotropic hardening law of the model [16,17,18]. In order to improve the prediction of ratcheting strain of materials in the constitutive model, Halama [19] modified the model proposed by Abdel-Karim and Ohno, and embedded the cyclic hardening/softening characteristics of materials and additional hardening caused by non-proportional load into the Abdel-Karim model. Chen et al. [20] introduced the evolution rule of function ϕ(p) and the parameter μi of material cyclic hardening or softening characteristics into the superposition model of the Ohno-Wang II model and A-F model, and improved the prediction effect of strain cyclic and uniaxial ratcheting strain. Wang et al. [21] found that both isotropic hardening and kinematic hardening contributed to the cyclic softening of materials, and modified the Abdel-Karim and Ohno (AKO) model by quantifying the effects of isotropic hardening and kinematic hardening on the cyclic softening of materials.
These studies demonstrated that different materials are of different ratcheting behaviors. S30408 is a metastable austenitic stainless steel. Under the influence of mechanical stress, plastic deformation and a cryogenic temperature service environment, martensitic transformation will be induced [22,23,24], resulting in significant changes in the mechanical behavior, thus, affecting the fatigue life of the materials. The high strength and toughness of S30408 stainless steel results from the transformation strengthening and plastic growth mechanism. Therefore, the influence of martensitic transformation should be considered in the establishment of a constitutive model.
Sitko et al. [25] studied the plastic behavior of stainless steel at cryogenic temperatures, and put forward a constitutive equation including martensite volume fraction. Yi et al. [26] found that cryogenic temperatures can improve the martensitic transformation rate of austenitic stainless steel AISI 316LN. Naghizadeh et al. [27] modified the expression of the martensite content by discussing different kinetic models of martensite transformation. Cryogenic temperature experiments of metastable austenitic stainless steel found that the martensite content in S30408 showed the rule of rapid growth first and then stability [28]. To correctly describe the effect of the martensite content on the material model, the parameters related to martensite can be associated with kinematic hardening and isotropic hardening in the constitutive model [29,30].
The constitutive models used in commercial software, such as ANSYS and ABAQUS, to predict ratcheting strain are relatively simple; however, these models are not ideal for predicting ratcheting strain. Many scholars [31,32,33,34] have written user subroutines by using the secondary development function provided by ANSYS and ABAQUS, and embedded some new constitutive models into the software to predict ratchet deformation, and achieved good prediction results.
In this paper, we tested the ratcheting of this material under different loading conditions. Based on the Ohno-Wang II model, the content of martensite is related to the isotropic hardening and kinematic hardening of materials, and we propose a constitutive model that can reasonably describe the martensitic transformation induced by S30408 at cryogenic temperatures. Then, in the numerical algorithm, we utilized the radial return method to realize the internal equilibrium iteration and determined the tentative strain after each equilibrium iteration by the consistent tangent matrix. Finally, we applied the constitutive model with martensitic transformation to the finite element analysis to simulate the ratcheting strain of S30408 round bar. We examined the capability of the proposed constitutive model to predict the uniaxial ratcheting strain of the specimen.

2. Ratcheting Experiments

2.1. Specimen

The material used in the study was S30408 stainless steel, which was formed as a round bar specimen with a gauge length of 10 mm and a diameter of 3 mm. The chemical composition of the material is (wt.%): C 0.02, Si 0.45, Mn 1.16, P 0.026, S 0.001, Cr 18.13, Ni 8.03, and N 0.05. The specimen dimensions used in the experiments are shown in Figure 1.
In the experiments, liquid nitrogen was used as the coolant, and the round bar specimen was completely immersed in the cryogenic temperatures box. When the temperature was reduced to 110 K, this was maintained for 15 mins to obtain a uniform temperature of the specimen. During the cryogenic temperature experiments, the temperature deviation was ±2 K.

2.2. Uniaxial Experiments

2.2.1. Uniaxial Tension Experiments

The uniaxial tension experiments were conducted with the electric servo fatigue multifunctional machine (Xi’an Li Chuang company, Xi’an, China), and the loading rate was 1.5 × 10−3/s. Figure 2 shows the uniaxial tension curve of the material. Under cryogenic temperatures, a plastic plateau occurs due to the formation of martensite, and secondary hardening can be observed after the plastic plateau due to the deformation-induced martensite and the strong interactions of dislocations of the austenitic with the hardening phase martensite. Although the nominal yield strength σ0.2 was 330 MPa, the initial yield stress σ0 corresponding to the size of yield surface was only 120 MPa, which is used in the later plasticity analysis. It can be seen from Figure 2 that the strain corresponding to the yield platform is 8.3%. In the actual engineering structure, the cyclic load will not exceed 8.3%; thus, the segment of plasticity including the yield platform but without the second hardening section is considered in the subsequent analysis. The mechanical property parameters are shown in Table 1.

2.2.2. Uniaxial Ratcheting Experiments

We performed the uniaxial ratcheting experiments with a EUM - 25K20 tension - torsion fatigue testing machine (CARE Measure & Control Ltd. Co., Tianjin, China). The triangular fluctuating load was applied to the round bar specimen. The load was controlled with a load rate of 312 MPa/s. The ratcheting strain was measured by an Epsilon cryogenic temperature extensometer (Epsilon Tech. Comp., Jackson, WY 83001, USA). The volume fraction of martensite was determined by a ferrite measuring instrument, FMP30 (HELMUT FISCHER GMBH Co., Hunenberg, Germany). The loading conditions of the ratcheting experiments are shown in Table 2.
Figure 3 shows the stress-strain response of materials under different loading conditions. The strains accumulated in the direction of the mean stress and the accumulating rates decreased continuously with the increasing cyclic number. In order to quantitatively study the law of plastic strain accumulation under cyclic loading, in this paper, ratcheting strain is defined as:
ε r = ε max + ε min 2
where εmax and εmin are the maximum and minimum plastic engineering strain in a cycle.
The ratcheting strains for various loading conditions are shown in Figure 4. The ratcheting strain of materials accumulated rapidly in the first few cycles, and with the increase of cycles, the ratcheting deformation rate gradually showed a slow trend. As can be seen from Figure 4, when the mean stress remained unchanged, the ratcheting strain increased with the increase of amplitude stress.

3. Constitutive Model and Simulation Results with MATLAB

3.1. Several Nonlinear Constitutive Models

Five advanced constitutive models were tentatively used to simulate the ratcheting strain of cryogenic temperatures for S30408 austenitic stainless steel.
  • Ohno-Wang model
To better describe ratcheting behavior, Ohno and Wang [6] proposed an equation representing the critical state of the dynamic recovery term, and assumed that each component of backstress α i only works when its value reaches the critical value (OW-I). The OW-I model is proposed as follows
α = i = 1 M α i , d α i = 2 3 r i γ i d ε p γ i H ( f i ) d ε p : α i α i ¯ α i .
where f i = α i ¯ 2 r i 2 . However, the OW-I model produces closed hysteresis loops and hence cannot simulate uniaxial ratcheting. To eliminate this limitation, Ohno and Wang [7] proposed a slight nonlinearity for each rule by introducing an exponential relation, and before reaching its critical state, the dynamic recovery term is partially operative (OW-II). The formula is proposed as follows
α = i = 1 M α i , d α i = 2 3 r i γ i d ε p γ i ( α i ¯ r i ) m i d ε p : α i α i ¯ α i
2.
AKO model and its modification
To improve the description of the uniaxial ratcheting behavior of the OW-I model, Abdel-Karim [8] combined the OW-I model with the A-F model (AKO). The AKO model is proposed in the following form
α = i = 1 M α i , d α i = 2 3 r i γ i d ε p μ i γ i α i d p γ i H ( f i ) d λ i α i .
where d λ i = d ε p : α i r i μ i d p .
Hamala [19] revised the μi in the AKO model (AKO IV), where the μi is expressed as follows
μ i = η f σ : α i α i ¯ χ , d η = d η 1 + d η 2 .
d η 1 = ω 1 ( η 1 η 1 ) d p , d η 2 = ω 2 ( η 2 η 2 ) d p .
where μi is the ratcheting parameter, η is a parameter affecting the ratchet deformation, η∞1, η∞2, η1, η2, ω1 and ω2 are material constants, dp is the magnitude of the plastic strain increment tensor, and χ is the multiaxial parameter. The uniaxial ratcheting studied in this study is χ = 0.
3.
OW II-AF and its modified model
Juan Zhang [10] superimposed the A-F model and the OW-II model (OW II-AF) to improve the prediction of the uniaxial ratcheting effect of the OW-II model. The formula is proposed as follows
α = i = 1 M α i , d α i = 2 3 r i γ i d ε p μ i γ i α i d p γ i ( α i ¯ r i ) m i d λ i α i .
Xiaohui Chen [20] proposed the modified OW II-AF model (M OW II-AF) by superimposing the OW-II model and the A-F model. The kinematic hardening is given by
α = i = 1 M α i , d α i = γ i ( 2 3 r i d ε p μ i ϕ ( p ) α i d p ϕ ( p ) ( α i ¯ r i ) m i d ε p : α i r i / ϕ ( p ) α i ) .
where ϕ(p) is cyclic hardening/softening function of materials ϕ ( p ) = ϕ + ( 1 ϕ ) e ω ϕ p , and ϕ is the stable value of Marquis.

3.2. Parameter Determination and Ratchet Simulation of Several Typical Kinematic Hardening Models

The initial yield stress σ0 of the material was 120 MPa, and the elastic modulus was the slope of the elastic section. By fitting the linear stress-strain curve, the slope was the elastic modulus E = 194 GPa. The material parameters γi and ri are determined by the following formula through the uniaxial tension curve.
r i = ( σ i σ i 1 ε i p ε i 1 p σ i + 1 σ i ε i + 1 p ε i p ) ε i p , i 1 , γ i = 1 ε i p .
Finally, ri is obtained by
i = 1 M r i + σ 0 = σ m a x .
where σi and εip are the stress and plastic strain values corresponding to the endpoints of each section. σ0 is the stress value when the plastic strain value on the curve is 0. The number of components that are needed for evaluation of the accuracy of the backstress are set as M = 8.
The parameters of the models are shown in Table 3.
With the incremental loading method in MATLAB, the ratchetting strains were simulated with the above five models for the uniaxial ratcheting experiments discussed in Section 2.2.2. The simulation effects of these models on ratcheting strain of S30408 stainless steel are shown in Figure 5. The simulation of ratcheting strain by these models had low prediction at the early stage of the cycle and over prediction at the later stage of the cycle. The yield platform stress of the material was 580 MPa. Therefore, when simulating the ratcheting with a maximum stress of 580 MPa, the ratcheting of the material could no longer be accurately obtained through the stress control experimental method. Thus, the ratcheting of the mean stress of 290 MPa and the amplitude stress of 290 MPa will not be discussed further.

3.3. The Constitutive Model with Martensite Transformation

The Garion model [30] used the linear mixed hardening model to calculate the backstress. The content of deformed martensite affects the linear kinematic hardening and dynamic recovery term, which can better simulate the secondary hardening phenomenon under uniaxial loading. However, it is difficult to describe the primary hardening phenomenon of materials after they enter plasticity. The accurate simulation of the ratcheting strain depends on the simulation effect of the model on the hysteresis loop. As the Garion model cannot describe hysteresis well, it cannot accurately describe the ratcheting strain. However, the description of the content of deformed martensite and its influence on the hardening of the material by this model is worth referencing.
The volume fraction of martensite is the most popular macroscopic internal variable specifying the growth of the martensitic phase [35,36]. Bain strain is added to the strain calculation due to the influence of the volume strain caused by the martensitic transformation. The martensite content and cumulative plastic strain values are obtained by experiments, and the evolution law of martensitic content is obtained by fitting the curve. In this paper, the martensite content dependent variable ξ is correlated with the Ohno Wang II model and isotropic hardening model, which can simulate the ratcheting of S30408. We propose the following constitutive model.
  • According to the study of Garion, we assumed that strain ε is divided into elastic part εe inelastic part εp and phase transformation part εbs
    ε = ε e + ε p + ε bs .
  • The elastic part obeys Hooke’s law
    σ = E : ( ε ε p ε bs ) .
  • The plastic part can be stated as
    d ε p = 3 2 d s · n h s n .
  • The martensitic transformation part can be expressed in terms of relative volume change Δv, due to the phase transformation [29], as
    ε bs = 1 3 Δ v ξ I = 1 3 ( V m V a ) V a ξ I .
    where Va and Vm represent the unstresses specific volumes occupied by the austenite and martensite, respectively. The value of Δv is approximately 0.02–0.05, depending on the alloy composition [24].
  • We assumed the material follow von Mises yield criterion, which can be given by
    f ( σ , α , R ) = J 2 ( σ α ) σ y R .
  • The backstress α is altered by the induced martensite, and the evolution law is postulated in the following form:
    α = i = 1 M α i .
    d α i = d α i a + d α i a ξ = γ i ( 2 3 r i d ε p ( α i ¯ r i ) m i d ε p : α i α i ¯ α i ) + h ξ ε p .
    where h is a parameter that depends on the material and can be determined by trial. Due to the existence of martensite, the backstress increment can be regarded as the sum of pure austensite d α i a and the interaction between the dislocations into the austenitic matrix and the martensite d α i a ξ [29].
  • The presence of martensite also affects the cyclic hardening law of materials, and the evolution law is postulated in the following general form:
    d R = F ( ξ ) d p = ( R ( ξ ) - R ) d p .
    According to [28], austenitic stainless steel S30408 is characterized by rapid hardening due to martensitic transformation. The martensite content of the material increases rapidly with the increase of cumulative plastic strain, and then tends to be stable. The occurrence of martensite transformation has a certain influence on the cyclic hardening of the material, which can be attributed to isotropic hardening.
    d R = 2 3 ( 1 ξ ) b ( Q ξ R ) d p .
    where R is the isotropic hardening parameter, and p is the cumulative plastic strain. The term (1 − ξ) is added to compensate for the assumption that the martensite is elastic [29]. In reality, the martensite shall rather be considered elasto-plastic. Therefore, the contribution to the hardening of the material is slightly smaller.
  • According to [28], the evolution law of the martensite content of S30408 austenitic stainless steel at 110 K presents a trend of rapid increase at first and then stability, which can be expressed as follows.
    d ξ = ω m ( ξ ξ ) d p .
By changing the rate of reaching the saturation value ωm, the value of deformation induced martensite due to the cumulative plastic strain can be controlled. ωm is mainly influenced by the temperature, the loading rate, etc., which requires further investigations.

3.4. Determination of Constitutive Model Parameters

In addition to the parameters in the Ohno-Wang model determined above, the model parameters related to the martensitic transformation can be determined as follows.
Cryogenic temperature torsion experiments were performed on specimens with a torque amplitude of 2.5 Nm and a cycle of 3 s [37]. It can be seen from Figure 6 that S30408 austenitic stainless steel showed rapid cyclic hardening in the first 10 cycles, and then the hardening rate decreased until 200 cycles. The maximum angular amplitude was approximately twice that of the stable angle. Therefore, S30408 is a cyclic hardening material at 110 K, and tends to be stable when the number of cycles exceeds 200.
In this study, the cyclic hardening properties of S30408 austenitic stainless steel are all attributed to isotropic hardening. According to the isotropic hardening rule, the peak stress and accumulated plastic strain of the symmetric strain cycle meet the requirement
σ p e a k = σ Q ξ [ 1 exp ( b p ) ] .
where σpeak is the peak stress. As shown in Figure 7a, the isotropic parameters C and Q can be obtained by fitting the peak stress and cumulative plastic strain. The parameters related to the martensitic transformation are obtained by fitting the martensite content and cumulative plastic strain curve, as shown in Figure 7b, and the saturation value of the martensite content ξ = 0.9, which is influenced by the temperature. Parameter ωm is determined by trial at present, as shown in Figure 7c. The material parameters of the model are presented in Table 4.

3.5. Model Prediction Results

With the incremental loading method by MATLAB, the uniaxial tension and uniaxial ratchetting behavior were simulated. With the material parameters determined above, it can be seen from Figure 8 that the model has a good fit for the uniaxial tension curve of the material.
From Figure 9, the predicted stress-strain response of this constitutive model is not greatly different from the experimental value, and the stress-strain response of different loading conditions can be reasonably predicted.
It can be seen from Figure 10 that the ratcheting strain predictions are quite effective for the load case with a mean stress of 290 MPa and amplitude stress of 240 MPa and for the load case with a mean stress of 290 MPa and amplitude stress of 215 MPa. The predictions are a little higher for the other two load cases, which are further from the load case used for the model parameter determination. However, this does overcome the disadvantages of early lower prediction and later over prediction by other models. Although there is little difference between the Ohno-Wang II model and the proposed model in predicting uniaxial ratcheting strain, the Ohno-Wang II model always overpredicts ratcheting strain in late cycles, and is not as effective as the proposed model. This is because the body centered tetragonal martensite is much harder than the face centered cubic austenite, which can improve the hardness of materials. It is necessary to consider the martensitic transformation.
As the experiments were conducted at cryogenic temperatures, to observe any change of martensite content at any time was not convenient. Therefore, by measuring the content of the deformation-induced martensite after the experiments and comparing it with the simulated values, the deviation between the predicted value of martensite content and the experimental value can be observed in Figure 11, within the acceptable range. Therefore, we considered that the model could reasonably predict the content of the deformation induced martensite at cryogenic temperatures.

4. Finite Element Implementations

In this section, we detail the finite element implementation of the proposed model with martensitic transformation. In general, the nonlinear equation is solved by Newton-Raphson (N-R) iteration (Figure 12) for every load step, such as σ1 and σ2, in which the trial strain increment is determined by the consistent tangent matrix. However, in the local N-R iteration, namely from each tangential prediction point back to the yield surface, we used the substitution method, namely the radial backward substitution (Figure 13).

4.1. The Basic Equation of Cyclic Plasticity Constitutive Model with Martensite Transformation

  • If the material obeys von Mises yield criterion, the yield function of the material is
    f = 3 2 ( s α ) T [ M 1 ] ( s α ) .
  • Plasticity flow rule
    d ε p = 3 2 d s n h s n .
  • Kinematic hardening law
    α = i = 1 M α i .
    d α i = γ i ( 2 3 r i d ε p ( α i ¯ r i ) m i d ε p : α i α i ¯ α i ) + h ξ ε p .
  • Isotropic hardening rule
    σ y = σ y 0 + R .
    d R = 2 3 ( 1 ξ ) b ( Q ξ R ) d p .
  • Evolution of the martensite content
    d ξ = ω m ( ξ ξ ) d p .

4.2. Solution of Plastic Strain by Euler Backward Method

Given all constitutive variables at t = t n , the convergent solutions contain σ n , α n , ε n p , ε n b s , ξ n , and R n .
The elastic tentative stress:
σ n + 1 T = D e : ε n + 1 T .
The tentative strain:
ε n + 1 T = ε n + 1 ε n p ε n b s .
Due to the elastic isotropy and plastic incompressibility, taking the deviatoric part of tentative stress σ n + 1 T
S n + 1 T = σ n + 1 T 1 3 t r ( σ n + 1 T ) I s .
According to the radial return algorithm [38], the deviatoric stress on the yield surface corresponds to the tentative deviatoric stress
S n + 1 = S n + 1 T 2 G [ M 2 ] ( Δ ε n + 1 p + Δ ε n + 1 b s ) .
Bain strain increment Δ ε n + 1 b s
Δ ε n + 1 b s = 1 3 Δ v Δ ξ n + 1 I s .
Martensite content Δ ξ n + 1
Δ ξ n + 1 = ω m ( ξ ξ n ) Δ p n + 1 .
Isotropic hardening variable Δ R n + 1
Δ R n + 1 = 2 3 ( 1 Δ ξ n + 1 ) b ( Q Δ ξ n + 1 R n ) Δ p n + 1 .
Normal direction of the yield surface
n n + 1 = S n + 1 n n + 1 S n + 1 n n + 1 = 3 2 S n + 1 n n + 1 σ y ( n + 1 ) Δ p n + 1 .
If the yield surface function is f n + 1 = 0 , the plastic strain increment can be calculated by the following formula
Δ ε n + 1 p = 3 2 Δ p [ M 1 ] n n + 1 = 3 2 S n + 1 n n + 1 σ y ( n + 1 ) Δ p n + 1 .
Plastic multiplier in the martensitic transformation model can be derived as follows:
α n + 1 = α n + i = 1 M Δ α n + 1 ( i ) .
Δ α n + 1 ( i ) = 2 3 r i γ i [ M 2 ] Δ ε n + 1 p + 2 3 r i γ i h Δ ξ n + 1 [ M 2 ] Δ ε n + 1 p γ i ( α n + 1 ( i ) ¯ r i ) m i Δ ε n + 1 p : α n + 1 ( i ) α n + 1 ( i ) ¯ α n + 1 ( i ) .
α n + 1 ( i ) = α n ( i ) + 2 3 r i γ i [ M 2 ] Δ ε n + 1 p + 2 3 r i γ i h Δ ξ n + 1 [ M 2 ] Δ ε n + 1 p γ i ( α n + 1 ( i ) ¯ r i ) m i Δ ε n + 1 p : α n + 1 ( i ) α n + 1 ( i ) ¯ α n + 1 ( i ) .
α n + 1 ( i ) = θ n + 1 ( i ) ( α n ( i ) + 2 3 r i γ i Δ p n + 1 [ M 1 ] n n + 1 + 2 3 r i γ i Δ p n + 1 h Δ ξ n + 1 [ M 1 ] n n + 1 ) .
θ n + 1 ( i ) = 1 1 + γ i ( α n + 1 ( i ) ¯ r i ) m i n n + 1 : α n + 1 ( i ) α n + 1 ( i ) ¯ Δ p n + 1 .
The yield function is:
f y t r i a l = 3 2 ( s α ) T [ M 1 ] ( s α ) ( σ y ( n + 1 ) - R n + 1 ) .
If f y t r i a l 0 , the material does not yield at this strain increment, σ n + 1 = σ n + 1 T .
If f y t r i a l > 0 , the material yields and enters into the plastic state. Under this strain increment condition, there is:
σ n + 1 = σ n + 1 T D e : Δ ε n + 1 p .
Equation (41) subtracted from Equation (32) gives
S n + 1 α n + 1 = S n + 1 T 2 G [ M 2 ] ( Δ ε n + 1 p + Δ ε n + 1 b s ) i = 1 M θ n + 1 ( i ) ( α n ( i ) + 2 3 r i γ i ( h Δ ξ n + 1 + 1 ) [ M 2 ] Δ ε n + 1 p ) .
Substitution of Equation (37) into Equation (44) provides
S n + 1 α n + 1 = ( S n + 1 T 2 G [ M 2 ] Δ ε n + 1 b s i = 1 M θ n + 1 ( i ) α n ( i ) ) ( 1 + 3 G Δ p n + 1 / σ y ( n + 1 ) i = 1 M r i γ i ( h Δ ξ n + 1 + 1 ) θ n + 1 ( i ) Δ p n + 1 / σ y ( n + 1 ) ) .
S n + 1 - α n + 1 satisfies yield function, thus the following formula can be obtained:
Δ p = [ 3 2 ( S n + 1 T 2 G [ M 2 ] Δ ε n + 1 b s i = 1 M θ n + 1 ( i ) α n ( i ) ) T [ M 1 ] ( S n + 1 T 2 G [ M 2 ] Δ ε n + 1 b s i = 1 M θ n + 1 ( i ) α n ( i ) ) ] 1 2 σ y ( n + 1 ) 3 G + i = 1 M r i γ i ( h Δ ξ n + 1 + 1 ) θ n + 1 ( i ) .
When plastic multipliers meet the following limitations, the results converge
| Δ p ( k + 1 ) Δ p ( k ) | e p s .
The following formulas can be obtained when the iterative calculation converges
Δ ε n + 1 p = Δ p ( k + 1 ) 3 2 [ M 1 ] n n + 1 .
Δ ε n + 1 b s = 1 3 Δ v Δ ξ n + 1 I s .
Δ ξ n + 1 = ω m ( ξ ξ n ) Δ p n + 1 .
σ n + 1 = D e : ( ε n + 1 tr Δ ε p ) ε b s .
ε n + 1 p = ε n p + Δ ε n + 1 p .
ε n + 1 = ε n + 1 e + ε n + 1 p + ε n + 1 b s .
Δ R n + 1 = 2 3 ( 1 Δ ξ n + 1 ) b ( Q Δ ξ n + 1 R n ) Δ p n + 1 .
R n + 1 = R n + Δ R n + 1 .
σ y ( n + 1 ) = σ y 0 + R n + 1 .
The above process is represented by a flow chart in Figure 14.

4.3. Consistent Tangent Modulus and Calculation of Strain Increment

The consistent tangent modulus [39], d Δ σ n + 1 / d Δ ε n + 1 , is calculated for the martensitic transformation constitutive model.
d Δ σ n + 1 d Δ ε n + 1 = D e 4 G 2 L n + 1 1 : I d .
L n + 1 = 2 G I + i = 1 M H n + 1 ( i ) + 2 3 ( d Y d p ) n + 1 n n + 1 n n + 1 + 2 3 Y n + 1 Δ p n + 1 ( I n n + 1 n n + 1 ) .
d α n + 1 ( i ) = 2 3 r i γ i ( 1 + h ξ ) θ n + 1 ( i ) [ M 2 ] d Δ ε n + 1 p + d θ n + 1 ( i ) θ n + 1 ( i ) α n + 1 ( i ) .
To obtain H n + 1 ( i ) , we use the differentiating Equation (41).
Differentiation Equation (42) and combined with Equation (37), d θ n + 1 ( i ) , can be expressed by
d θ n + 1 ( i ) = 2 3 θ n + 1 ( i ) 2 γ i ( α n + 1 ( i ) ¯ r i ) m i n n + 1 T : α n + 1 ( i ) α n + 1 ( i ) ¯ : d Δ ε n + 1 p .
The substitution of Equation (61) into Equation (60) provides
d α n + 1 ( i ) = 2 3 r i γ i θ n + 1 ( i ) ( [ M 2 ] ( 1 + h ξ ) m n + 1 ( i ) ( α n + 1 ( i ) ¯ r i ) m i ( m n + 1 ( i ) ) T ) : d Δ ε n + 1 p .
where
m n + 1 ( i ) = 3 2 α n + 1 ( i ) r i
or
m n + 1 ( i ) = 3 2 α n + 1 ( i ) α n + 1 ( i ) ¯ , α n + 1 ( i ) ¯ = 3 2 { α n + 1 ( i ) } T [ M 1 ] { α n + 1 ( i ) } .
The constitutive model, H n + 1 ( i ) , with martensitic transformation can be calculated by the following formula.
H n + 1 ( i ) = d α n + 1 ( i ) d Δ ε n + 1 p = 2 3 r i γ i θ n + 1 ( i ) ( [ M 2 ] ( 1 + h ξ ) ( α n + 1 ( i ) ¯ r i ) m i n n + 1 T : α n + 1 ( i ) α n + 1 ( i ) ¯ ) . .
The substitution of Equation (23) into Equation (39) provides
d α n + 1 ( i ) = ( 1 + h ξ ) r i γ i d s n n + 1 h s n n + 1 γ i ( α n + 1 ( i ) ¯ r i ) m i n n + 1 T : α n + 1 ( i ) α n + 1 ( i ) ¯ α n + 1 ( i ) d p .
The expression of hs, in the martensitic transformation model, can be obtained by the method provided in [40].
Y = σ y = σ y 0 + R .
h s = i = 1 M r i γ i ( 1 + h ξ ) 3 2 n n + 1 T [ M 1 ] i = 1 M γ i α n + 1 ( i ) n n + 1 T : α n + 1 ( i ) α n + 1 ( i ) ¯ + σ y p .
d R = 2 3 ( 1 ξ ) b ( Q ξ R ) d p .
d ξ = ω m ( ξ ξ ) d p .
Thus, the nonlinear numerical calculation of the constitutive model with martensitic transformation can be realized.

4.4. Embedment of Constitutive Model in the ANSYS Program and Ratchetting Prediction

The constitutive model with the martensitic transformation was embedded in ANSYS to form the user’s ANSYS executive file by compiling and connecting the subroutine userspl.f of ANSYS. The programming was verified by the uniaxial tension and used to predict the ratchetting strain for uniaxial loading.

4.4.1. Uniaxial Tension Verification

Based on the proposed model, the uniaxial tensile behavior of the material was calculated using the user version of ANSYS. For simplicity, only one element PLANE42 with a side length of 1 was used for the calculation, as shown in Figure 15. The x-direction displacement constraint was applied to the left side of the model, the y-direction displacement constraint was applied to the lower side, and the -direction load was applied to the right side.
The predictions of the uniaxial tension are plotted in Figure 16 together with those predicted by the incremental loading method of MATLAB. The comparison results show that the results of the finite element method with a single element were consistent with the experimental results and MATLAB predicted results.

4.4.2. Uniaxial Ratcheting Prediction

We simulated the ratcheting strains of the specimen under four different uniaxial loading conditions, first taking the mean stress of 290 MPa and the amplitude stress of 240 MPa as an example. The finite element model of the specimen was established with element SOLID45. In the finite element model, the bottom surface was fully constrained, and the axial load of 353–3746 N was cyclic and applied on the master node of the rigid region for the top surface. Figure 17 indicates the boundary condition of the finite element model. The cyclic loading will generate a uniform cyclic stress of 50–530 MPa in the measuring cross section of the specimen.
As shown in Figure 18, the maximum plastic strain of the specimen appeared in its middle section. Therefore, the ratcheting strain at the middle of the specimen was selected for comparison.
Figure 19 shows the uniaxial ratcheting curve for the ratcheting strain of the specimen under different loading conditions by ANSYS. It can be seen from Figure 19 that the simulation of the uniaxial ratcheting strain with ANSYS was consistent with the experimental results and the MATLAB predictions.

5. Conclusions

In this paper, we experimentally studied the uniaxial ratcheting behavior of S30408 austenitic stainless steel at 110 K. Based on the Ohno-Wang II model, we proposed a new cyclic plasticity constitutive model with martensitic transformation. This was embedded into ANSYS with verification. The main results of the present work can be summarized as follows:
  • The proposed constitutive model with martensitic transformation reasonably predicted the ratcheting strain of S30408 at cryogenic temperatures. In the model, the influences of deformation-induced martensite on both isotropic hardening and kinematic hardening are considered.
  • By comparing the simulated and the experimental values of the deformation-induced martensite content under different loading conditions, we found that the proposed model could reasonably predict the martensite content at cryogenic temperatures.
  • The determination method of the model parameter ωm to control the rate of reaching the martensite saturation value requires further investigations, and may be influenced by the temperature, the loading rate, etc.

Author Contributions

Carried out the experiment: Y.C. and K.Z.; performed the numerical simulations: Y.C. and X.C. (Xiaohui Chen); methodology: Y.C. and K.Z.; conceptualization: B.G. and X.C. (Xu Chen); writing—original draft preparation: Y.C.; writing—review and editing: X.C. and C.Y.; project administration: B.G; and funding acquisition: B.G. All authors provided critical feedback and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key Research and Development Program of China, FC08086002017YFC0805600, and the Natural Science Foundation of Hebei Province of China, No. E2017202078.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

ε total strain tensor
ε e elastic strain tensor
ε p inelastic strain tensor
ε b s Bain strain tensor
ξ martensite content
σ stress tensor (MPa)
E Young’s modulus for elasticity (GPa)
υ Poisson’s ratio
σ 0 size of yield surface (MPa)
σ 0.2 yield strength (MPa)
σ b tensile strength(MPa)
δ elongation (%)
f yield surface function
J 2 second deviatoric stress tensor invariant
R isotropic hardening parameter
α deviatoric backstress tensor
α i the ith component of deviatoric backstress α
α i ¯ the magnitude of αi
γ i , r i material constants
m i material constants
s deviatoric stress
H Heaviside step function, x ≥ 0, H(x) = 1; x < 0, H(x) = 0
Macauley’s bracket, x ≥ 0, < x > = x; x < 0, < x > = 0
b rate of change of R
Q saturation value of R (MPa)
ξ 0 initial value of martensite evolution
ξ limit value of martensite evolution
ω m evolution rate of martensite
[ M 1 ] non-zero auxiliary diagonal matrix, [M1] = diag [1,1,1,2,2,2]
[ M 2 ] non-zero auxiliary diagonal matrix, [M2] = diag [1,1,1,1/2,1/2,1/2]
D e isotropic elastic stiffness tensor
h s plastic modulus
I s second-rank unit tensor
n normal direction of yield surface
H n + 1 ( i ) fourth-rank tensor

References

  1. Shi, X.P. Gas and LNG pricing and trading hub in East Asia: An introduction. Nat. Gas Ind. 2016, 3, 352–356. [Google Scholar] [CrossRef]
  2. Gao, B.J.; Wang, B.; Zhai, L.H.; Shang, C.W. Random vibration analysis of the strip-supported cryogenic storage tank in transporting. J. Hebei Univ. Technol. 2018, 47, 48–52. [Google Scholar]
  3. Armstrong, P.J.; Frederick, C.-O. A mathematical representation of the multiaxial bauschinger effect. Mater. High. Temp. 1966, 24, 1–26. [Google Scholar]
  4. Chaboche, J.L. Time-independent constitutive theories for cyclic plasticity. Int. J. Plast. 1986, 2, 149–188. [Google Scholar] [CrossRef]
  5. Chaboche, J.L. On some modifications of kinematic hardening to improve the description of ratcheting effects. Int. J. Plast. 1991, 7, 661–678. [Google Scholar] [CrossRef]
  6. 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, 375–390. [Google Scholar] [CrossRef]
  7. Ohno, N.; Wang, J.D. Kinematic hardening rules with critical state of dynamic recovery, Part II: Application to experiments of ratchetting bahavior. Int. J. Plast. 1993, 9, 391–403. [Google Scholar] [CrossRef]
  8. Abdel-Karim, M.; Ohno, N. Kinematic hardening model suitable for ratcheting with steady-state. Int. J. Plast. 2000, 16, 225–240. [Google Scholar] [CrossRef]
  9. Chen, X.; Jiao, R. Modified kinematic hardening rule for multiaxial ratcheting prediction. Int. J. Plast. 2004, 20, 871–898. [Google Scholar] [CrossRef]
  10. Zhang, J. Constitutive Description for Non-Proportionally Ratcheting of Cyclically Hardening Material and Its Finite Element Implementation at High Temperatures. Ph.D. Thesis, Southwest Jiaotong University, Chengdu, China, 2002. [Google Scholar]
  11. Chen, G.; Zhao, X. Constitutive modelling on the whole-life uniaxial ratcheting behavior of sintered nano-scale silver paste at room and high temperatures. Microelectron. Reliab. 2018, 80, 47–54. [Google Scholar] [CrossRef]
  12. Liu, X.Y.; Dong, Y.W.; Xiao, X. Study on constitutive model of cyclic ratcheting considering memory recovery of plastic strain. Chin. J. Sol. Mech. 2018, 39, 429–437. [Google Scholar]
  13. Varvani-Farahani, A. A comparative study in descriptions of coupled kinematic hardening rules and ratcheting assessment over asymmetric stress cycles. Fatigue Fract. Eng. Mater. 2017, 40, 882–893. [Google Scholar] [CrossRef]
  14. Rahman, S.M. Finite Element Analysis and Related Numerical Schemes for Ratcheting Simulation. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, 2006. [Google Scholar]
  15. Rahman, S.M.; Hassan, T.; Corona, E. Evaluation of cyclic plasticity models in ratcheting simulation of straight pipes under cyclic bending and steady internal pressure. Int. J. Plast. 2008, 24, 1756–1791. [Google Scholar] [CrossRef]
  16. Hassan, T.; Kyriakides, S. Ratcheting of cyclically hardening and softening materials. Int. J. Plast. 1994, 10, 149–212. [Google Scholar] [CrossRef]
  17. Kang, G.; Ohno, N.; Nebu, A. Constitutive modeling of strain range dependent cyclic hardening. Int. J. Plast. 2003, 19, 1801–1819. [Google Scholar] [CrossRef]
  18. Kalnins, A.; Rudolph, J.; Willuweit, A. Using the nonlinear kinematic hardening material model of Chaboche for elastic-plastic ratcheting analysis. J. Press. Vessel Technol. 2015, 137, 031006. [Google Scholar] [CrossRef]
  19. Halama, R. A modification of abdelkarim-ohno model for ratcheting simulations. Tech. Gaz. 2008, 15, 3–9. [Google Scholar]
  20. Chen, X.H. A Study on the Ratcheting and Shakedown Behavior of Pressurized Piping Components under Cyclic Loading. Ph.D. Thesis, Tianjin University, Tianjin, China, 2013. [Google Scholar]
  21. Wang, Q.; Liu, X. Non-saturated cyclic softening and uniaxial ratcheting of a high-strength steel: Experiments and viscoplastic constitutive modeling. Mech. Mater. 2017, 113, 112–125. [Google Scholar] [CrossRef]
  22. Olson, G.B.; Cohen, M. Kinetics of strain-induced martensitic nucleation. Met. Trans. 1975, 6A, 791–795. [Google Scholar] [CrossRef]
  23. Stringfellow, R.G.; Parks, D.M.; Olson, G.B. A constitutive model for transformation plasticity accompanying strain-induced martensitic transformations in metastable austenitic steels. Acta. Metall. Mater. 1992, 40, 1703–1716. [Google Scholar] [CrossRef]
  24. Tomita, Y.; Iwamoto, T. Constitutive modeling of trip steel and its application to the improvement of mechanical properties. Int. J. Mech. Sci. 1995, 37, 1295–1305. [Google Scholar] [CrossRef]
  25. Sitko, M.; Skoczeń, B. Effect of γ-α’ phase transformation on plastic adaptation to cyclic loads at cryogenic temperatures. Int. J. Solids Struct. 2012, 49, 613–634. [Google Scholar] [CrossRef] [Green Version]
  26. Xiong, Y.; Yue, Y.; He, T.T.; Fan, M.X.; Ren, F.Z.; Cao, W. Effect of rolling temperature on microstructure evolution and mechanical properties of AISI316LN austenitic stainless steel. MDPI Mater. 2018, 11, 1557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Naghizadeh, M.; Mirzadeh, H. Modeling the kinetics of deformation-induced martensitic transformation in AISI 316 metastable austenitic stainless steel. Vacuum 2018, 157, 243–248. [Google Scholar] [CrossRef]
  28. Ji, S.Y. A Study on Cyclic Plasticity and Fatigue Mechanism of Strain Hardening S30408 Austenitic Stainless Steel. Master’s Thesis, Tianjin University, Tianjin, China, 2019. [Google Scholar]
  29. Garion, C.; Skoczen, B. Modeling of plastic strain-induced martensitic transformation for cryogenic applications. Int. J. Appl. Mech. 2002, 69, 755–762. [Google Scholar] [CrossRef]
  30. Garion, C.; Skoczen, B.; Sgobba, S. Constitutive modelling and identification of parameters of the plastic strain-induced martensitic transformation in 316L stainless steel at cryogenic temperatures. Int. J. Plast. 2006, 22, 1234–1264. [Google Scholar] [CrossRef]
  31. Gao, B.J. Modeling of Material Multiaxial Ratcheting and Ratcheting Prediction of Pressure Piping. Ph.D. Thesis, Tianjin University, Tianjin, China, 2005. [Google Scholar]
  32. Hassan, T.; Zhu, Y.; Matzen, V.-C. Improved ratcheting analysis of piping components. Int. J. Press. Vessel. Pip. 1998, 75, 643–652. [Google Scholar] [CrossRef]
  33. Kolasangiani, K.; Varvani-Farahani, A. Local ratcheting assessment of steel samples undergoing various step and block loading conditions. Theor. Appl. Fract. Mech. 2020, 107, 102533. [Google Scholar] [CrossRef]
  34. Kang, G. A visco-plastic constitutive model for ratcheting of cyclically stable materials and its finite element implementation. Mech. Mater. 2004, 36, 299–312. [Google Scholar] [CrossRef]
  35. Fischer, F.D. Transformation induced plasticity in triaxially loaded steel specimens subjected to a martensitic transformation. Eur. J. A Solids. 1992, 11, 233–244. [Google Scholar]
  36. Cherkaoui, M.; Berveiller, M.; Sabar, H. Micromechanical modeling of martensitic transformation induced plasticity (TRIP) in austenitic single crystals. Int. J. Plast. 1998, 14, 597–626. [Google Scholar] [CrossRef]
  37. Li, L.L. Study on Low Temperature Ratcheting Effect of Pre-Strain S30408. Master’s Thesis, Hebei University of Technology, Tianjin, China, 2019. [Google Scholar]
  38. Simo, J.C.; Taylor, R.L. Consistent tangent operators for rate-independent elastoplasticity. Comput. Method. Appl. Mech. Eng. 1985, 48, 101–118. [Google Scholar] [CrossRef]
  39. Kobayashi, M.; Ohno, N. Implementation of cyclic plasticity models based on a general form of kinematic hardening. Int. J. Numer. Meth. Eng. 2002, 53, 2217–2238. [Google Scholar] [CrossRef]
  40. Chen, X.H.; Chen, X.; Yu, D.J.; Gao, B.J. Recent progresses in experimental investigation and finite element analysis of ratcheting in pressurized piping. Int. J. Press. Vessel. Pip. 2013, 101, 113–142. [Google Scholar] [CrossRef]
Figure 1. Specimen geometry (mm).
Figure 1. Specimen geometry (mm).
Applsci 10 06002 g001
Figure 2. Uniaxial tension curve.
Figure 2. Uniaxial tension curve.
Applsci 10 06002 g002
Figure 3. The stress-strain response in the first 50 cycles under different loading conditions: (a) σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240 MPa (d) σm = 240 MPa σa = 240 MPa (e) σm = 290 MPa σa = 290 MPa.
Figure 3. The stress-strain response in the first 50 cycles under different loading conditions: (a) σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240 MPa (d) σm = 240 MPa σa = 240 MPa (e) σm = 290 MPa σa = 290 MPa.
Applsci 10 06002 g003aApplsci 10 06002 g003b
Figure 4. Experimental values of ratcheting.
Figure 4. Experimental values of ratcheting.
Applsci 10 06002 g004
Figure 5. Ratcheting strain prediction of different loading conditions: (a)σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240 MPa (d) σm = 240 MPa σa = 240 MPa (e) σm = 290 MPa σa = 290 MPa.
Figure 5. Ratcheting strain prediction of different loading conditions: (a)σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240 MPa (d) σm = 240 MPa σa = 240 MPa (e) σm = 290 MPa σa = 290 MPa.
Applsci 10 06002 g005
Figure 6. Angular amplitude vs. cycles.
Figure 6. Angular amplitude vs. cycles.
Applsci 10 06002 g006
Figure 7. Parameter determination of the constitutive model: (a) isotropic hardening, (b) martensite transformation, (c) determination of parameter ωm.
Figure 7. Parameter determination of the constitutive model: (a) isotropic hardening, (b) martensite transformation, (c) determination of parameter ωm.
Applsci 10 06002 g007
Figure 8. Simulation of uniaxial tension.
Figure 8. Simulation of uniaxial tension.
Applsci 10 06002 g008
Figure 9. Simulation and experiment values of stress-strain response: (a) σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240MPa (d) σm = 240MPa σa = 240 MPa.
Figure 9. Simulation and experiment values of stress-strain response: (a) σm = 290 MPa σa = 190 MPa (b) σm = 290 MPa σa = 215 MPa (c) σm = 290 MPa σa = 240MPa (d) σm = 240MPa σa = 240 MPa.
Applsci 10 06002 g009
Figure 10. The ratcheting strain under different working conditions by the proposed model.
Figure 10. The ratcheting strain under different working conditions by the proposed model.
Applsci 10 06002 g010
Figure 11. Martensite content prediction under different loading conditions by the proposed model.
Figure 11. Martensite content prediction under different loading conditions by the proposed model.
Applsci 10 06002 g011
Figure 12. Sketch of the Newton-Raphson (N-R) iteration.
Figure 12. Sketch of the Newton-Raphson (N-R) iteration.
Applsci 10 06002 g012
Figure 13. Sketch of the radial return method.
Figure 13. Sketch of the radial return method.
Applsci 10 06002 g013
Figure 14. Flow diagram of the proposed model.
Figure 14. Flow diagram of the proposed model.
Applsci 10 06002 g014
Figure 15. Uniaxial tension of the FEA model.
Figure 15. Uniaxial tension of the FEA model.
Applsci 10 06002 g015
Figure 16. Uniaxial tensile verification curve.
Figure 16. Uniaxial tensile verification curve.
Applsci 10 06002 g016
Figure 17. Boundary conditions applied to the specimen finite element model.
Figure 17. Boundary conditions applied to the specimen finite element model.
Applsci 10 06002 g017
Figure 18. The total plastic strain of the finite element model.
Figure 18. The total plastic strain of the finite element model.
Applsci 10 06002 g018
Figure 19. The simulation results of uniaxial ratcheting for the specimen.
Figure 19. The simulation results of uniaxial ratcheting for the specimen.
Applsci 10 06002 g019
Table 1. The mechanical properties.
Table 1. The mechanical properties.
Mechanical Properties
E/GPaυσ0/MPaσ0.2/MPaσb/MPaδ/%
1940.33120330132475.4
Table 2. Ratcheting experiments scheme.
Table 2. Ratcheting experiments scheme.
Mean Stress
σm/MPa
Amplitude Stressσa/MPa
190215240290
240
290
Table 3. Parameter determination of different models.
Table 3. Parameter determination of different models.
ModelParameter
OW-IIr(1-8) 80, 56, 69, 72, 80, 61, 31, 15 γ(1-8) 8000, 1500, 350, 180, 150, 72, 65, 12  mi = 7
AKOr(1-8) 90, 72, 78, 60, 61, 58, 34, 25  γ(1-8) 8000, 1500, 350, 180, 150, 72, 65, 12  μi = 0.2
AKO IVr(1-8) 80, 56, 69, 72, 80, 61, 31, 15  γ(1-8) 8000, 1500, 350, 180, 150, 72, 65, 12
η01 = 0.1, ω1 = 5, η∞2 = 0.05, χ = 0  η02 = 0.04, ω2 = 3, η∞2 = 0.03
OWII-AFr(1-8) 80, 56, 69, 72, 80, 61, 31, 15  γ(1-8) 8000, 1500, 350, 180, 150, 72, 65, 12
μi = 0.04, mi = 7
M OWII-AFr(1-8) 80, 56, 69, 72, 80, 61, 31, 15 γ(1-8) 8000, 1500, 350, 180, 150, 72, 65,
12 mi = 7 η0 = 0.04, ω = 8, η = 0.01, ϕ = 0.4, ωϕ = 1 χ = 0
Table 4. Material parameters of the constitutive model.
Table 4. Material parameters of the constitutive model.
StyleParameter
Uniaxial Tensionr(1-8) 80, 56, 69, 72, 80, 61, 31, 15  γ(1-8) 8000, 1500, 350, 180, 150, 72, 65, 12
Uniaxial Ratchetmi=7
Martensite ContentΔv = 0.05  h = 0.1  ωm = 1  ξ0 = 0.03 ξ= 0.9
Isotropic HardeningQ = 560  b = 10

Share and Cite

MDPI and ACS Style

Chen, Y.; Chen, X.; Gao, B.; Chen, X.; Zhang, K.; Yu, C. A Cyclic Plasticity Model with Martensite Transformation for S30408 and Its Finite Element Implementation. Appl. Sci. 2020, 10, 6002. https://doi.org/10.3390/app10176002

AMA Style

Chen Y, Chen X, Gao B, Chen X, Zhang K, Yu C. A Cyclic Plasticity Model with Martensite Transformation for S30408 and Its Finite Element Implementation. Applied Sciences. 2020; 10(17):6002. https://doi.org/10.3390/app10176002

Chicago/Turabian Style

Chen, Yanan, Xiaohui Chen, Bingjun Gao, Xu Chen, Kai Zhang, and Chulin Yu. 2020. "A Cyclic Plasticity Model with Martensite Transformation for S30408 and Its Finite Element Implementation" Applied Sciences 10, no. 17: 6002. https://doi.org/10.3390/app10176002

APA Style

Chen, Y., Chen, X., Gao, B., Chen, X., Zhang, K., & Yu, C. (2020). A Cyclic Plasticity Model with Martensite Transformation for S30408 and Its Finite Element Implementation. Applied Sciences, 10(17), 6002. https://doi.org/10.3390/app10176002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop