You are currently viewing a beta version of our website. If you spot anything unusual, kindly let us know.

Preprint
Article

Dynamic Mechanical Properties and Modified Johnson-Cook Model Considering Recrystallization Softening for Nickel Based Powder Metallurgy Superalloys

Altmetrics

Downloads

110

Views

31

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

05 December 2023

Posted:

07 December 2023

You are already at the latest version

Alerts
Abstract
The material undergoes high temperature and high strain rate deformation process during cutting process, may induce the dynamic recrystallization behavior and result the evolution of dynamic mechanical properties of the material to be machined. In this paper, the modified Johnson-Cook(J-C) model for Nickel based powder metallurgy superalloy considering dynamic recrystallization behavior in high strain rate and temperature is proposed. The dynamic mechanical properties of the material under different strain rates and temperature conditions are obtained by quasi-static compression test and split Hopkinson pressure bar (SHPB) test. The coefficients of the modified J-C model are obtained by linear regression method. The modified model is verified by comparison with experimental and model prediction results. The results show that the modified J-C model proposed in this paper can accurately describe the mechanical properties of nickel-based powder metallurgy superalloys at high temperature and high strain rate. Which provides help for studying the cutting mechanism and finite element simulation of nickel-based powder metallurgy superalloy.
Keywords: 
Subject: Engineering  -   Industrial and Manufacturing Engineering

1. Introduction

Comparing with cast and wrought superalloy, the powder metallurgy (PM) superalloy has widely used in advanced aero-engine turbine disks [1,2] due to its homogeneous microstructure, fine grains, higher yield strength and no mace segregation [3,4], etc. Although PM is a near-forming process, machining is still a necessary process to meet the final accuracy requirements. However, because of low thermal conductivity and low elastic modulus, PM superalloy are deemed too difficult to machining materials [5,6].
The external power is converted into deformation energy in the cutting process. At the same time, the heat is not easy to release, causing the tool and workpiece material temperature rise, resulting in material softening [7]. On the other hand, with the increase of strain, the material work hardening occurs. The interaction between the two makes the material strength not only a function of strain, but also a function of strain rate and deformation temperature, which is called the constitutive relation equation [8,9]. Establishing the accurate constitutive relationship of the workpiece material in the cutting deformation zone is the key to study the cutting deformation mechanism.
The constitutive equation for the cutting process is mainly divided into two categories, namely phenomenological models and physical models [10]. Compared with physical models, phenomenological models have become the widely used one due to its simple form and wide applicability [11]. As the most widely used phenomenological model, the Johnson-Cook (J-C) model [12] is given in the form of three-part product, representing strain hardening, strain rate hardening and thermal softening, respectively. However, the above three items of the J-C model are considered to be independent of each other, and the effects of the microstructure of the material itself and the adiabatic temperature rise on the flow stress are not considered [13,14,15].
The occurrence of dynamic recrystallization in the cutting process has been confirmed by researchers [16,17]. Dynamic recrystallization not only refines the grain size of the workpiece material, but also produces a softening phenomenon that cannot be ignored in the flow stress of the material during the deformation process.
A constitutive equation that can accurately describe the microscopic deformation process is still to be realized. Calamaz M et al. [18] established a JC-TANH modified model to describe the phenomenon of soft flow caused by dynamic recrystallization of cutting titanium alloy materials, and applied it to the finite element software. This model was verified by Ulutan [18] in software Deform. Other similar corrections to the J-C constitutive equation can be found in the literature [19,20].
However, the above modified models only use the temperature (Tm/2) condition as the critical condition for recrystallization. According to the modified J-C model of DENGUIR et al. [21], the recrystallization strain threshold was proposed. This critical condition takes into account the combined effects of temperature and strain rate.
The cutting process is a strong coupling process of transient high strain, high strain rate and temperature. The occurrence of recrystallization is very important to the dynamic mechanical properties of PM superalloy. However, the current research on the modified constitutive equation considering recrystallization for difficult-to-machine materials focuses on titanium alloys. The research on the dynamic mechanical properties of PM superalloy at high temperature and high strain rate, especially the influence of dynamic recrystallization, is still very scarce.
In this paper, the dynamic mechanical properties of powder metallurgy superalloy were analyzed by quasi-static test and SHPB test. A modified constitutive model considering the dynamic recrystallization behavior of the cutting process is established. The parameters of the modified model are obtained by linear fitting of the experimental data. Finally, it is verified that the model has good accuracy.

2. Experimentals

2.1. Specimen preparation

The FGH96 superalloy used in this study is prepared by hot isostatic pressing (HIP) and its composition is shown in Table 1. The experiment method is consisted with two parts, which are quasi-static state compression tests and SHPB tests.
The samples sizes for quasi-static state compression tests and SHPB tests are φ3mm×3mm andφ2mm×2mm, respectively. By grinding the end face and outer surface of the sample obtained by wire cutting, the sample has higher dimensional accuracy and lower surface roughness (Ra 0.8mm). In addition, it is necessary to ensure that the parallelism of the two ends is less than 0.01 mm, so as to reduce the experimental error caused by the sample.

2.2. Quasi-Static State Compression Tests

The quasi-static compression test of FGH96 superalloy is carried out in the electric universal testing machine. All the experiments are carried out at room temperature. The strain rate is set at 0.001s-1 by controlling the loading speed (0.2mm/min). The force and displacement in the experiment are output by the control software, and then processed by Origin software to obtain the stress-strain curve. In order to ensure the reliability of the data, three repeated tests were carried out for each specimen, and the average value was obtained.

2.3. SHPB Tests

The SHPB equipment system are schematically illustrated in Figure 1. In order to obtain higher strain rate, the diameter of the incident bar and the transmitted bar is reduced to 5mm.
The strains of sample could be recorded by strain gauges, which are attachment on the incident bar and transmitted bar. According to the one-dimensional stress wave theory, the engineering strain εeng, engineering stress σeng and strain rate ε ˙ eng sample can be calculated by Equation (1)-(3). The temperature is controlled by resistance-heated furnace. To ensure the accuracy of the data, three experiments are repeated for each sample.
Preprints 92339 i001
Preprints 92339 i002
Preprints 92339 i003
Where, c0 is the wave speed in the incident bar, E and A are the Young’s modulus and the cross-sectional area of the bars. A0 and l0 are the cross-sectional area and the length of the cylindrical specimen. εr and εt are the reflected and transmitted strains, respectively, when the reflected and transmitted waves propagate independently.
Under the condition that the material is incompressible, the true strain εT, true stress σT and true strain rate of ε ˙ T the specimen can be obtained, by Equation(4)-(6).
Preprints 92339 i004
Preprints 92339 i005
Preprints 92339 i006
The experimental conditions of SHPB tests are shown in Table 2.

3. Modified Constitutive Model

The modified model proposed in this paper considers the coupling effect of strain, strain rate and temperature and the influence of microstructure transformation on the basis of JC model, as shown in Equation (7).
Preprints 92339 i007
Where, A represents the yield stress, B and n represent the material strain hardening coefficients, ε represents the true strain, C represents the viscosity coefficient, ε ˙  represents the strain rate, 𝜀0 represents the reference strain rate, T represents the temperature, Tm represents the material melting temperature, T0 represents the room temperature, and m represents the thermal softening coefficient. The effect of dynamic recrystallization on flow stress is characterized by adding a correction item H(ε, ε ˙ , T) to J-C model. The mathematical expressions of H(ε, ε ˙ , T) are shown in Equation (8)-(10). Considering that the recrystallization threshold is affected by both temperature and strain rate, it is expressed by Equation (11).
Preprints 92339 i008
Preprints 92339 i009
Preprints 92339 i010
Preprints 92339 i011
Where, hi (i=0,1,2) and ri(i =0, 1…5) are material constants, εr is the critical strain for the occurrence of dynamic recrystallization, which is expressed by Equation (11).

4. Results and Discussion

4.1. Dynamic Mechanical Properties

4.1.1. Characteristics of Stress-Strain Curve

The stress-strain curve of FGH96 superalloy obtained by quasi-static compression test is shown in Figure 2. According to Figure 2, in the first stage of compression process (OA), the material undergoes elastic deformation with the stress increases linearly with the increase of strain. When the strain reaches point A, the growth rate of stress slows down with the increase of strain, indicating that the material enters the plastic deformation stage (AB), showing obvious work hardening phenomenon. When the strain exceeds the strain at point B, the material enters the damage deformation stage which is marked by the decrease of stress with the increase of strain.
According to the above analysis, there is no obvious yield phenomenon in the quasi-static compression process of FGH96 superalloy. The yield stress of bulk samples was measured by 0.2% strain shift method, that is, the stress value when 0.2% plastic deformation occurs is taken as its yield strength (point A in Figure 2). The yield strength of FGH96 superalloy is σ0.2=773 MPa.
The effect of strain rate on the stress-strain curve of FGH96 superalloy obtained by SHPB tests at different temperatures (25, 200, 400, 600, 700, 800) is shown in Figure 3. According to Figure 3a, when the temperature is 25℃, the maximum strains in the plastic deformation stage corresponding to the strain rates of 4000s-1,6000s-1,10000s-1 and 12000s-1 are 0.3, 0.4, 0.6 and 0.7, respectively. That is to say, the plastic deformation stage at room temperature increases significantly with the increase of strain rate, and the increase is as high as 200 %. According to Figure 3b–f, under high temperature deformation conditions (200, 400℃, 600, 700, 800), compared with the lower strain rate level, the plastic deformation stage of the material under high strain rate conditions also increases significantly. In summary, the PM superalloy has obvious plasticizing effect in the deformation process over a wide range of strain rate.
The effect of temperature on the stress-strain curve of FGH96 superalloy at different strain rates (4000s-1, 6000s-1, 10000s-1, 12000s-1) is shown in Figure 4. According to Figure 4a–d, at the same strain rate, the stress in the plastic deformation stage decreases with the increase of temperature, showing obvious temperature softening effect.
Table 3 shows the yield stress at different temperatures and strain rates under SPHB test conditions. Compared with the yield stress (773MPa) obtained by quasi-static compression test at room temperature, the yield stress obtained by SPHB tests at high temperature and high strain rate increases significantly. When the strain rate is 12000s-1, the yield strength reaches 1616MPa. According to Table 3, when the temperature is 200, 400, 600, 700 and 800, the yield strength increases with the increase of strain rate. It is worth noting that the flow stress of the material also increases with the increase of strain rate (Figure 3), which indicates that the PM superalloy material exhibits a significant strain rate hardening effect.
In addition, under the same strain rate loading conditions, the yield strength of the PM superalloy decreases with the increase of temperature. Combined with Figure 4, the flow stress and yield strength decrease with the increase of temperature, which further confirms the temperature softening effect of the PM superalloy.
Table 3. Yield strength at different temperatures and strain rates (MPa).
Table 3. Yield strength at different temperatures and strain rates (MPa).
Strain rate (s-1) Temperature, T ()
25 200 400 600 700 800
4000 1315 1305 1234 1097 903 883
6000 1359 1318 1280 1142 1001 926
10000 1544 1544 1322 1227 1182 994
12000 1616 1616 1355 1271 1205 1043

4.1.2. Strain-Hardening

In order to quantitatively analyze the strain hardening phenomenon of PM superalloy, the strain hardening rate Q is calculated by Equation (12).
Preprints 92339 i012
Where, εi and σi represent the strain and the stress of the ith experiment, respectively. εi−1 and σi−1 represent the strain and the stress of the (i-1)th experiment.
Based on the quasi-static compression experimental data, the strain hardening rate-strain curve of the PM superalloy obtained according to the Equation (12) is shown in Figure 5. According to Figure 5, in the elastic deformation stage, the strain hardening rate Q of the material decreases sharply with the increase of strain. In the plastic deformation stage (AB), the decrease rate of strain hardening rate slows down, and at point B, it decreases to 0 (it is no longer in the hardening state). In the damage deformation stage, the strain hardening rate is less than 0, and the material is damaged and deformed. Therefore, it can be concluded that during the quasi-static compression test, the strain hardening phenomenon of FGH96 superalloy is significant, and with the continuous increase of strain, the strain hardening effect gradually weakens, and reaches equilibrium at point B. This is because in the process of quasi-static compression deformation, with the increase of strain, thermal softening occurs inside the material[22].
When the temperature is 700°C, the strain hardening rate-strain curve of the PM superalloy at different strain rates is shown in Figure 6. According to Figure 6, the PM superalloy still exhibits obvious strain hardening under high temperature and high strain rate conditions. However, it is worth noting that in the plastic deformation stage, the strain hardening rate Q of the material at strain rate of 4000s-1 always remains positive. In contrast, the Q of the material at the strain rate of 6000s-1 continues to decrease and reaches a negative state after reaching zero. This is because the material not only has the thermal softening effect, but also superimposes the recrystallization softening effect, so that the hardening effect in the plastic deformation stage is greatly weakened, and the flow softening phenomenon appears [23].

4.1.3. Adiabatic induced Increase of Temperature

According to Section 4.1.1, the PM superalloy has obvious plasticizing effect during high temperature and high strain rate deformation. Researchers have found that the adiabatic temperature rise during material deformation is one of the important factors that cannot be ignored in the plasticizing effect [24]. In addition, the adiabatic temperature rise makes the actual temperature inside the material higher than the experimental loading temperature, which reduces the dislocation slip resistance and strengthens the internal softening of the material. Under the experimental conditions in this paper, the adiabatic temperature rise ΔT in the process of material deformation can be obtained by Equation (13).
Preprints 92339 i013
Where, η is the plastic work-heat conversion coefficient, for the experimental deformation conditions in this paper, it take 0.9 [25]; ρ is the material density; Cp is the specific heat capacity at atmospheric pressure. The specific heat capacity Cp of the PM superalloy FGH96 superalloy at different temperatures is shown in Table 4.
Table 4. Specific heat capacity of FGH96 superalloy at different temperatures.
Table 4. Specific heat capacity of FGH96 superalloy at different temperatures.
Temperature, T (℃) 25 200 400 600 700 800
Cp (kJ/ (kg‧K)) 0.391 0.422 0.455 0.487 0.503 0.525
Figure 7 shows the influence of loading temperature on adiabatic temperature rise of FGH96 superalloy under different strain rates. According to Figure 7, at the same strain rate, the adiabatic temperature rise decreases with the increase of temperature. For example, at the strain rate of 12000s-1, as the loading temperature increases from 25 to 800, the adiabatic temperature rise decreases from 100 to 30℃ (reduced by 70%). At the same time, when the strain rate is 4000s-1, the adiabatic temperature rise decreases from 35 to 15 with the increase of loading temperature (reduced by 60%). It can be seen that with the decrease of strain rate, the rate of adiabatic temperature rise decreasing with the increase of loading temperature is obviously reduced. From the above, with the decrease of strain rate, the reduction rate of adiabatic temperature rise decreases significantly with the increase of loading temperature. At the same loading temperature, the higher the strain rate, the higher the adiabatic temperature rise. For example, at 200, the corresponding adiabatic temperature rise at strain rates of 4000 s-1, 6000 s-1, 10000 s-1 and 12000 s-1 are 30, 40, 65 and 70, respectively. At the same time, compared with the temperature of 25, 200 and 400, the adiabatic temperature rises at 600, 700 and 800 decreases significantly with the increase of strain rate.

4.1.4. The Strain Rate Sensitivity

According to Section 4.1.1, FGH96 superalloy exhibits strain rate strengthening effect in the plastic deformation stage. In order to describe the degree of strain rate strengthening of the material, the strain rate sensitivity coefficient q is introduced, as shown in Equation (14). The larger the strain rate sensitivity coefficient q, the stronger the strain rate sensitivity of the material.
Preprints 92339 i014
Figure 8 shows the influence of strain on the strain rate sensitivity coefficient of FGH96 superalloy at different temperatures. According to Figure 8, the strain rate sensitivity coefficient decreases with the increase of strain under the experimental temperature conditions. For example, when the temperature is 600, the strain rate sensitivity coefficient decreases from 0.11 to 0.09 with the increase of strain (from 0.1 to 0.4). It is worth noting that the influence of temperature on the strain rate sensitivity coefficient is not obvious. For example, for the temperature conditions of 200 and 400, the strain rate sensitivity coefficient is lower than the strain rate sensitivity coefficient of 25. However, when the temperature is greater than 600, the strain rate sensitivity coefficient becomes larger and increases with the increase of temperature. At the same time, when the temperature is 800 and the strain is 0.1, the biggest strain rate sensitivity coefficient is 0.155 (less than 0.2), which shows that FGH96 superalloy exhibits weak strain rate sensitivity.

4.1.5. The Temperature Sensitivity

According to section 4.1.1, FGH96 superalloy shows a strong temperature softening effect in the plastic deformation stage. In order to quantitatively describe the temperature sensitivity of the material, the temperature sensitivity coefficient s is introduced, as shown in Equation (15). The greater the temperature sensitivity coefficient s, the stronger the temperature sensitivity of the material.
Preprints 92339 i015
The variation of the temperature sensitivity coefficient of FGH96 superalloy with temperature at different strain rates calculated by Equation (15) is shown in Figure 9. According to Figure 9, the temperature sensitivity coefficient increases significantly with the increase of temperature under the same strain rate. When the strain rate is 4000s-1, the temperature sensitivity coefficient is 0.1 at 200, and reaches 1.2 at 800℃, which is increased by 1100%. At the same time, the temperature sensitivity coefficient decreases with the increase of strain rate. When the temperature is 800, the temperature sensitivity coefficient is 1.3 at 4000s-1, and decreases to 0.9 at 12000-1, with a decrease of 30%. From the above, the temperature sensitivity of the PM superalloy increases significantly with the increase of temperature and decreases with the decrease of strain rate. In summary, FGH96 superalloy exhibits strong temperature sensitivity.

4.2. Construction of Constitutive Model

4.2.1. Recrystallization Critical Condition

The effect of dynamic recrystallization on the stress-strain curve of the material is shown in Figure 10. If there is no dynamic recrystallization during the deformation process, the flow stress of the material increases slowly with the increase of strain, as shown in the black solid line in Figure 10. When the strain reaches the critical strain of dynamic recrystallization, the flow stress decreases, which is manifested as the recrystallization softening effect as shown in the red line in Figure 10. However, when the strain reaches a certain value (critical strain εr), the material undergoes dynamic recrystallization. At this time, the flow stress shows a significant downward trend, as shown in the red solid line in Figure 10, which is the flow softening phenomenon. Therefore, the determination of the critical condition of dynamic recrystallization, that is, the critical strain, is the key to the study of dynamic recrystallization flow softening.
The dynamic recrystallization critical strain of the PM superalloy obtained under the experimental conditions in this paper is shown in Table 5. It can be seen from the data in Table 5 that the critical strain of dynamic recrystallization of PM superalloy is not only related to the deformation temperature, but also related to the strain rate, which verifies the conclusion of Denguir [21]. The critical strain decreases with the increase of temperature and decreases with the increase of strain rate.
Table 5. Critical strain εr under different temperature and strain rate conditions.
Table 5. Critical strain εr under different temperature and strain rate conditions.
strain rate (s-1) Temperature, T (℃)
200 400 600 700 800
6000
10000
12000
0.4233 0.3621 0.3211
0.4033 0.3640 0.3199 0.3144 0.2911
0.3999 0.3443 0.3091 0.3051 0.2698
Based on the data in Table 3, the fitting surface of the critical strain is obtained by polynomial fitting (Equation (11)), is shown in Figure 11.
The equation of calculating the critical strain of dynamic recrystallization of PM superalloy under experimental conditions in this paper is shown in Equation (16).
Preprints 92339 i016

4.2.2. Identification of Constitutive Model’s Coefficients

(1)
Linear regression method
The modified J-C constitutive model (Equation (7)) represents strain hardening effect, strain rate strengthening effect, thermal softening effect and recrystallization softening effect from left to right. According to the linear regression parameter solving method, A, B, C, n, m and Hi are the parameters to be fitted, ε ˙ 0 , T0 and Tm are 0.001s-1, 25°C and 1350°C, respectively.
The strain hardening coefficient can be obtained by processing the quasi-static compression test data at room temperature. The proposed constitutive model is simplified as shown in Equation (17). The quasi-static compression tests permitted to determine the yield stress, which is represented by the coefficient A (Figure 12a). The Equation (18) was obtained by taking the logarithm on both sides of the Equation (17). The solution of the coefficients n and B can be obtained by a slope and an intercept of fitting straight line (Figure 12b). As shown in Figure 12, A = σ0.2 = 773 MPa, n = 0.667, B = 1271MPa.
Preprints 92339 i017
Preprints 92339 i018
Figure 12. Solution of coefficients A, B and n. a) Solution of coefficient A; b) Solution of coefficients B and n.
Figure 12. Solution of coefficients A, B and n. a) Solution of coefficient A; b) Solution of coefficients B and n.
Preprints 92339 g012
The strain rate sensitivity coefficient C, the thermal softening index m and the recrystallization softening correction coefficient Hi in the modified J-C constitutive model were obtained from the result of SHPB experiments.
According to SHPB tests at room temperature for different strain-rates, the constitutive equation simplified to Equation (19). The stress-strain curves of PM materials at different strain rates at room temperature in this paper are shown in Figure 13.
Preprints 92339 i019
The thermal softening coefficient m is determined according to Equation (20). The data of the SHPB tests under high temperature conditions were used. Finally, the fitting relationship between m and strain-rate is obtained as shown in Figure 14.
Preprints 92339 i020
The coefficients hi (i=0, 1, 2) related to the dynamic recrystallization can be obtained according to Equation (21).
Preprints 92339 i021
In summary, the modified constitutive model obtained by the linear regression method is shown in Table 6.
(2)
Function iteration method
According to Function iteration method [26,27], the proposed model is shown in Equation (22). The iterative function method determines the coefficients of the prediction model by continuously iterating the function through the optimization algorithm.
Preprints 92339 i022
The stress values corresponding to the room temperature T0 and the reference strain rate ε ˙ 0 conditions are selected as the initial values. The quasi-static compression test data at room temperature are selected for polynomial fitting to obtain the results of f (ε), which are shown in Figure 15a. Then, according to the relationship between measured stress in SHPB experiment at room temperature and f (ε), the relationship between stress and strain rate is obtained by iteration, which is f ( ε ˙ ), which are shown in Figure 15b. Finally, the results of f (T) and f (Hi) can be obtained according to the high temperature SHPB experimental data by iteration process, as shown in Figure 15c,d, respectively.
In the iterative process, the error R2 is used to judge the accuracy of the results. When the result meet Equation (23), the result is considered to be accepted.
Preprints 92339 i023
Where, R k 2 is the error of determination of the ith iterated function, R k 1 2 is the error of determination of the (i-1)th iterated function.
Figure 15. Function iteration method model solution results. a) Optimization results of f (ε); b) Optimization results of f ( ε ˙ ); c) Optimization results of f (T); d) Recrystallization fitting f (Hi).
Figure 15. Function iteration method model solution results. a) Optimization results of f (ε); b) Optimization results of f ( ε ˙ ); c) Optimization results of f (T); d) Recrystallization fitting f (Hi).
Preprints 92339 g015
Finally, the constitutive model constructed by the function iteration method is obtained, which is shown in Equation (24).
Preprints 92339 i024
(3)
Comparison of different methods
In order to optimize the solution method of the parameters, the stress-strain curves obtained by the constitutive equations solving by the two methods described in the previous section are compared with the experimental results. Under the condition of strain rate is 10000s-1 and the temperature is 25℃, 200, 400, 600, 700, 800, the comparison results of stress-strain curves are shown in Figure 16. Figure 16a,b are the experimental comparison results of linear regression solving method and functional iteration solving method, respectively. According to Figure 16, both of the equations obtained by linear regression solving method and functional iteration solving method can predict the trend of flow stress.
In order to quantitatively evaluate the overall error of the two methods, the scatter plot is used to calculate the correlation value. The results of the correlation between the calculated stress and the experimental stress obtained by the two methods are shown in Figure 17. Figure 17a,b are the results of linear regression solving method and functional iteration solving method, respectively. Compared with Figure 17b (R2 = 0.889), the data correlation index in Figure 17a is higher (R2 = 0.985), that is, the data concentration is higher.
The maximum relative error θ between the measured stress and the predicted stress is also calculated to evaluate prediction accuracy, as shown in Equation (25).
Preprints 92339 i025
Where, σp is the predicted stress, σm is the measured stress. The maximum relative error of the constitutive model obtained by the linear fitting method and the function iteration method are shown in Table 7. According to Table 7, the accuracy of the model obtained by the linear regression method (11.21%) is much higher than that obtained by the function iteration method (4.74%) in the non-dynamic recrystallization stage. Correspondingly, in the recrystallization stage, the accuracy of the model obtained by the function iteration method is improved (4.11%), which is very small compared with the accuracy of the model obtained by the linear regression method (5.11%). By calculating the average value of the maximum error θ - before and after dynamic recrystallization, the model accuracy obtained by the linear regression method is greater than the model obtained by the functional regression method.
In summary, the J-C constitutive equation of PM superalloy proposed in this paper solved by linear regression method is shown in Equation (26).
Preprints 92339 i026

4.3. Validation of Modified J-C constitutive Model

Temperature-dependent J-C model [28,29] is also widely used to predict the flow stress behavior of materials in cutting process, which is shown in Equation (27).
Preprints 92339 i027
Where, a, b are the coefficients related to flow softening caused by dynamic recrystallization. According to the stress-strain data of the PM superalloy under the experimental conditions in this paper, the linear regression method is used to solve Equation (24), and the established constitutive model is shown in Equation (28).
Preprints 92339 i028
The comparison between the predicted stress and the experimental stress of the two modified models at 10000s-1 strain rate are shown in Figure 18. Figure 18a,b are the comparison results of the model proposed in this paper and the Temperature-dependent J-C model, respectively.
The results of the correlation analysis between the calculated stress and the experimental stress obtained by the two models are shown in Figure 19. Figure 19a,b are the correlation analysis results of model proposed in this paper and Temperature-dependent J-C model, respectively. Compared with Figure 19b (R2=0.956), the data correlation index in Figure 19a is higher (R2=0.985), that is, the data concentration is higher.
According to Equation (25), the average values of the maximum relative error of the two models were calculated. The average value of the maximum relative errors of flow stress predicted by the modified J-C model in this paper and Temperature-dependent J-C model are 5.11% and 8.14%, respectively. Therefore, the constitutive model proposed in this paper can predict the flow stress of PM superalloy under the influence of recrystallization softening more accurately.

5. Conclusion

In this works, the dynamic mechanical properties of Nickel based PM superalloy are analyzed by quasi-static compression test and SHPB test at strain rate of 4000s-1 ~ 12000s-1 and temperature of 25°C ~ 800°C. A J-C constitutive model considering dynamic recrystallization softening effect is proposed. The following conclusions can be drawn in this paper.
(1) Dynamic mechanical properties of Nickel based PM superalloy was obtained by experiments. The PM superalloy has obvious plasticizing effect in the deformation process over a wide range of strain rate. The PM superalloy material exhibits a significant strain rate hardening effect and temperature softening effect. During the compression test, the strain hardening phenomenon of PM superalloy is significant and the strain hardening effect gradually weakens with the continuous increase of strain. FGH96 superalloy exhibits weak strain rate sensitivity and strong temperature sensitivity.
(2) Considering the influence of temperature and strain rate, the formula for the flow stress softening term H(ε, ε ˙ , T) and critical strain εr of dynamic recrystallization were obtained. On this basis, the modified J-C constitutive model of PM superalloy considering dynamic recrystallization behavior in a wide range of strain rate and temperature is proposed.
(3) The coefficients of the modified constitutive equation established in this paper were obtained by linear fitting and function iteration respectively. Through correlation analysis and maximum error analysis, the modified constitutive equation solved by the linear fitting method has higher accuracy.
(4) Compared with the temperature-dependent constitutive model, it is found that the modified constitutive model established in this paper is more suitable for the description of the stress-strain relationship of PM superalloy at high temperature and high strain rate.

Data Availability Statement

In this section, please provide details regarding where data supporting reported results can be found, including links to publicly archived datasets analyzed or generated during the study. Please refer to suggested Data Availability Statements in section “MDPI Research Data Policies” at https://www.mdpi.com/ethics. You might choose to exclude this statement if the study did not report any data.

Acknowledgments

The project is supported by National Natural Science Foundation of China (No. 52105460). This work is also supported by grants from China Postdoctoral Science Foundation (2022M71190), Shandong Provincial Key Research and Development Program (Major Scientific and Technological Innovation Project) (No.2020CXGC010204), Natural Science Foundation of Shandong Province (ZR202111150191), Key Laboratory of High-efficiency and Clean Mechanical Manufacture at Shandong University and Taishan Scholar Foundation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boyer R, Cotton J, Mohaghegh M, Schafrik R J M B. Materials considerations for aerospace applications[J]. 2015, 40, 1055–1066. [CrossRef]
  2. Yang L, Ren X, Ge C, Yan Q J I J o M R. Status and development of powder metallurgy nickel-based disk superalloys[J]. 2019, 110, 901–910. [CrossRef]
  3. Yang L, Ren X, Ge C, Yan Q. Status and development of powder metallurgy nickel-based disk superalloys[J]. International Journal of Materials Research. 2019, 110, 901–910. [CrossRef]
  4. Selvaraj S K, Sundaramali G, Jithin Dev S, et al. Recent Advancements in the Field of Ni-Based Superalloys[J]. Advances in Materials Science and Engineering. 2021, 2021, 1–60. [CrossRef]
  5. Zhang X, Chen Y, Hu J. Recent advances in the development of aerospace materials[J]. Progress in Aerospace Sciences, 2018, 97: 22-34. [CrossRef]
  6. Guo N, Leu M C. Additive manufacturing: technology, applications and research needs[J]. Frontiers of mechanical engineering, 2013, 8: 215-243. [CrossRef]
  7. Zhang H, Zeng H, Yan R, Wang W, Peng F J T I J o A M T. Analytical modeling of cutting forces considering material softening effect in laser-assisted milling of AerMet100 steel[J]. 2021, 113: 247-260. [CrossRef]
  8. Melkote S N, Grzesik W, Outeiro J, et al. Advances in material and friction data for modelling of metal machining[J]. 2017, 66, 731–754. [CrossRef]
  9. Ucun İ, Aslantas K J T I J o A M T. Numerical simulation of orthogonal machining process using multilayer and single-layer coated tools[J]. 2011, 54: 899-910. [CrossRef]
  10. Hou X, Liu Z, Wang B, et al. Stress-strain curves and modified material constitutive model for Ti-6Al-4V over the wide ranges of strain rate and temperature[J]. Materials. 2018, 11, 938. [CrossRef]
  11. Liu H, Xu X, Zhang J, et al. The state-of-the-art on numerical simulation of the effect of microstructure and its evolution in metal cutting processes[J]. International Journal of Machine Tools and Manufacture, 2022: 103890. [CrossRef]
  12. Johnson G, R. A constitutive model and data for materials subjected to large strains, high strain rates, and high temperatures[J]. Proc 7th Inf Sympo Ballistics, 1983: 541-547.
  13. He A, Xie G, Zhang H, Wang X. A comparative study on Johnson–Cook, modified Johnson–Cook and Arrhenius-type constitutive models to predict the high temperature flow stress in 20CrMo alloy steel[J]. Materials & Design (1980-2015), 2013, 52: 677-685. [CrossRef]
  14. Shokry A, Gowid S, Mulki H, Kharmanda G. On the Prediction of the Flow Behavior of Metals and Alloys at a Wide Range of Temperatures and Strain Rates Using Johnson–Cook and Modified Johnson–Cook-Based Models: A Review[J]. Materials. 2023, 16, 1574. [CrossRef]
  15. Su X, Wang G, Li J, Rong Y J S. Dynamic mechanical response and a constitutive model of Fe-based high temperature alloy at high temperatures and strain rates[J]. 2016, 5, 1–14. [CrossRef]
  16. Sakai T, Belyakov A, Kaibyshev R, Miura H, Jonas J J. Dynamic and post-dynamic recrystallization under hot, cold and severe plastic deformation conditions[J]. Progress in materials science, 2014, 60: 130-207. [CrossRef]
  17. Shin H, Ju Y, Choi M K, Ha D H. Flow stress description characteristics of some constitutive models at wide strain rates and temperatures[J]. . Technologies 2022, 10, 52. [CrossRef]
  18. Calamaz M, Coupard D, Girot F J I J o M T, Manufacture. A new material model for 2D numerical simulation of serrated chip formation when machining titanium alloy Ti–6Al–4V[J]. 2008, 48(3-4): 275-288. [CrossRef]
  19. Ulutan D, Sima M, Özel T J A M R. Prediction of machining induced surface integrity using elastic-viscoplastic simulations and temperature-dependent flow softening material models in titanium and nickel-based alloys[J]. 2011, 223: 401-410. [CrossRef]
  20. Rotella G, Umbrello D. Finite element modeling of microstructural changes in dry and cryogenic cutting of Ti6Al4V alloy[J]. Cirp Annals 2014, 63, 69–72. [CrossRef]
  21. Liang X, Liu Z, Wang B, Hou X J I J o M S. Modeling of plastic deformation induced by thermo-mechanical stresses considering tool flank wear in high-speed machining Ti-6Al-4V[J]. 2018, 140: 1-12. [CrossRef]
  22. Denguir L, Outeiro J, Fromentin G, Vignal V, Besnard R J J o M P T. A physical-based constitutive model for surface integrity prediction in machining of OFHC copper[J]. 2017, 248: 143-160. [CrossRef]
  23. Pan Z, Wu Z, Xiong J J P T. High-speed infrared imaging and mesostructural analysis of localized temperature rise in damage and failure behavior of 3-D braided carbon/epoxy composite subjected to high strain-rate compression[J]. 2019, 80: 106151. [CrossRef]
  24. Liang H, Guo H, Nan Y, et al. The construction of constitutive model and identification of dynamic softening mechanism of high-temperature deformation of Ti–5Al–5Mo–5V–1Cr–1Fe alloy[J]. 2014, 615: 42-50. [CrossRef]
  25. Guo Y, Ruan Q, Zhu S, et al. Temperature rise associated with adiabatic shear band: causality clarified[J]. 2019, 122, 015503. [CrossRef]
  26. Su X, Wang G, Li J, Rong Y J S. Dynamic mechanical response and a constitutive model of Fe-based high temperature alloy at high temperatures and strain rates[J]. 2016, 5, 1–14. [CrossRef]
  27. Lin Y, Chen X-M, Liu G J M S, A E. A modified Johnson–Cook model for tensile behaviors of typical high-strength alloy steel[J]. 2010, 527, 6980–6986. [CrossRef]
  28. Murugesan M, Jung D W J M. Johnson Cook material and failure model parameters estimation of AISI-1045 medium carbon steel for metal forming applications[J]. 2019, 12, 609. [CrossRef]
  29. Andrade U, Meyers M, Vecchio K, Chokshi A H J A m e m. Dynamic recrystallization in high-strain, high-strain-rate plastic deformation of copper[J]. 1994, 42, 3183–3195. [CrossRef]
  30. Ashrafian M M, Hosseini Kordkheili S A J P o t I o M E, Part G: Journal of Aerospace Engineering. A novel phenomenological constitutive model for Ti-6Al-4V at high temperature conditions and quasi-static strain rates[J]. 2021, 235, 1831–1842. [CrossRef]
Figure 1. Schematic diagram of SHPB equipment system.
Figure 1. Schematic diagram of SHPB equipment system.
Preprints 92339 g001
Figure 2. Stress-strain curve of FGH96 superalloy in the quasi-static compression test.
Figure 2. Stress-strain curve of FGH96 superalloy in the quasi-static compression test.
Preprints 92339 g002
Figure 3. Stress-strain curves of FGH96 superalloy over the wide range of temperature in the SHPB tests. a)25; b)200; c)400; d)600; e)700; f)800.
Figure 3. Stress-strain curves of FGH96 superalloy over the wide range of temperature in the SHPB tests. a)25; b)200; c)400; d)600; e)700; f)800.
Preprints 92339 g003
Figure 4. Stress-strain curves of FGH96 superalloy over the wide range of strain rate in the SHPB tests. a) 4000s-1; b) 6000s-1; c) 10000s-1; d) 12000s-1.
Figure 4. Stress-strain curves of FGH96 superalloy over the wide range of strain rate in the SHPB tests. a) 4000s-1; b) 6000s-1; c) 10000s-1; d) 12000s-1.
Preprints 92339 g004
Figure 5. Strain hardening rate-strain curve under quasi-static compression test conditions.
Figure 5. Strain hardening rate-strain curve under quasi-static compression test conditions.
Preprints 92339 g005
Figure 6. The strain hardening rate-strain curve at 700°C obtained by SHPB tests. a) 4000s-1; b) 6000s-1.
Figure 6. The strain hardening rate-strain curve at 700°C obtained by SHPB tests. a) 4000s-1; b) 6000s-1.
Preprints 92339 g006
Figure 7. The influence of loading temperature on adiabatic temperature rise during deformation under different strain rates.
Figure 7. The influence of loading temperature on adiabatic temperature rise during deformation under different strain rates.
Preprints 92339 g007
Figure 8. Strain rate sensitivity coefficient-strain curves of FGH96 superalloy at different temperatures.
Figure 8. Strain rate sensitivity coefficient-strain curves of FGH96 superalloy at different temperatures.
Preprints 92339 g008
Figure 9. Temperature sensitivity coefficient-temperature curves of FGH96 superalloy at different strain rate.
Figure 9. Temperature sensitivity coefficient-temperature curves of FGH96 superalloy at different strain rate.
Preprints 92339 g009
Figure 10. Flow softening of stress-strain curve induced by dynamic recrystallization.
Figure 10. Flow softening of stress-strain curve induced by dynamic recrystallization.
Preprints 92339 g010
Figure 11. Solving method and results of dynamic recrystallization critical strain.
Figure 11. Solving method and results of dynamic recrystallization critical strain.
Preprints 92339 g011
Figure 13. The results of strain rate sensitivity coefficient C. a) 4000s-1; b)6000s-1; c)10000s-1; d)12000s-1.
Figure 13. The results of strain rate sensitivity coefficient C. a) 4000s-1; b)6000s-1; c)10000s-1; d)12000s-1.
Preprints 92339 g013
Figure 14. Solution of thermal softening coefficient m.
Figure 14. Solution of thermal softening coefficient m.
Preprints 92339 g014
Figure 16. Comparison of experimental results and predicted results ( ε ˙ = 10000 s−1). a) Linear regression method; b) Function iteration method.
Figure 16. Comparison of experimental results and predicted results ( ε ˙ = 10000 s−1). a) Linear regression method; b) Function iteration method.
Preprints 92339 g016
Figure 17. Comparison of correlation between predicted stress and measured stress by different solving methods. a) Linear regression method; b) Function iteration method.
Figure 17. Comparison of correlation between predicted stress and measured stress by different solving methods. a) Linear regression method; b) Function iteration method.
Preprints 92339 g017
Figure 18. Comparison of correlation between predicted stress and measured stress by two models. (a) The modified J-C model proposed in this paper; (b) Temperature-dependent J-C model.
Figure 18. Comparison of correlation between predicted stress and measured stress by two models. (a) The modified J-C model proposed in this paper; (b) Temperature-dependent J-C model.
Preprints 92339 g018
Figure 19. Comparison of correlation between predicted stress and measured stress by different modified models. (a) The modified J-C model proposed in this paper; (b) Temperature-dependent J-C model.
Figure 19. Comparison of correlation between predicted stress and measured stress by different modified models. (a) The modified J-C model proposed in this paper; (b) Temperature-dependent J-C model.
Preprints 92339 g019
Table 1. Chemical composition of FGH96 superalloy.
Table 1. Chemical composition of FGH96 superalloy.
Element C Al Ti Cr Co Nb Mo W Ni
Wt. (%) 0.03 2.2 3.7 16 13 0.8 4 4 Bal
Table 2. Experimental conditions of SHPB tests.
Table 2. Experimental conditions of SHPB tests.
Variable Values
Strain rate ε ˙ (s-1) 4000, 6000, 10000, 12000
Temperature T (°C) 25, 200, 400, 600, 700, 800
Table 6. Constitutive model coefficients of linear regression method.
Table 6. Constitutive model coefficients of linear regression method.
Coefficients Value
A (MPa) 773
B (MPa) 1271
C 0.031
n 0.667
m Preprints 92339 i029
h0
h1
h2
-0.015
-0.015
0.046
Preprints 92339 i030 Preprints 92339 i031
Table 7. Comparison of the maximum relative errors of two different methods.
Table 7. Comparison of the maximum relative errors of two different methods.
Maximum relative errors Linear regression method Function iteration method
Preprints 92339 i032 4.74% 11.21%
Preprints 92339 i033 5.11% 4.11%
θ - 4.93% 7.66%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated