# A Characterization for the Deformation Behavior of As-Cast P91 Alloy Steel and Utilization in Hot Extrusion Process

Authors: Lu Jia; Yongtang Li; Yang Zhang
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

The constitutive models of as-cast P91 alloy steel were investigated by unidirectional thermal compression tests, which were conducted over a wide range of temperature and strain rate on Gleeble-3500 thermal simulation machine. The constitutive models were coupled with finite element method to simulate the hot extrusion process. Furthermore, the industrial hot extrusion experiment on actual pipe was conducted on 500 MN vertical extrusion machine and the load-displacement data during the extrusion were recorded. Comparisons between the numerical simulation and industrial experimental results were carried out and confirmed that the developed constitutive models can be used in the simulation and industrial research of compact hot extrusion process.


---

## Body

## 1. Introduction

Heavy caliber thick-wall pipes are widely used in the fields such as large supercritical thermal power, nuclear power, and petrochemical engineering. The current production process of heavy caliber thick-wall pipe is the punching and hot extrusion (PHE) process shown in Figure1. This PHE process has several defects including long production cycle, serious waste in material, and investments of large scale equipment. In order to solve these disadvantages, a new compact hot extrusion (CHE) process shown in Figure 2 was proposed by Li et al. [1]. During the CHE process, the as-cast billet is demolding under hot extrusion temperature and directly extruded on the extrusion machine. The short forming process realizes the integration on shape forming and performance improving. Thus, the hot deformation behaviors of as-cast billet during CHE process are required to be further investigated.Figure 1
The punching and hot extrusion process.Figure 2
The compact hot extrusion process.P91 martensite alloy steel (10Cr9Mo1VNb) is one of the leading materials of heavy caliber thick-wall pipe due to its high thermal conductivity, high coefficient of linear expansion, corrosion-resistant ability, and moderate price [2–6]. Investigations and great progress on the behaviors of as-forged P91 alloy steel have been done based on the experimental research and simulation results [7–11]. However, the reports of hot deformation performance on as-cast P91 alloy steel have not been seen yet.The object of this study is to establish an innovative constitutive model which can accurately be applied to the simulation and experimental research on the CHE process. The hot deformation behaviors of as-cast P91 alloy steel were investigated by hot compression tests. In addition, the constitutive models of as-cast P91 alloy steel were coupled with the finite element analysis of hot extrusion on heavy caliber thick-wall pipe. The comparisons between the simulation and the industrial hot extrusion experiment were then discussed with the load-displacement data.


## 2. Experimental

### 2.1. Materials and Methods

The chemical composition of as-cast P91 alloy steel was (wt.%) 0.084C -0.28Mn -0.38P -0.005S -0.28Si -8.63Cr -0.96Mo -0.209V -0.084Nb -0.21Ni -0.031N -(bal) Fe. The hot compression tests were carried out on a Gleeble-3500 thermomechanical simulator. Cylindrical specimens were machined with 8 mm in diameter and 12 mm in height.The specimens were deformed at the reduction of 60% in height. Different deformation temperatures (from 1173 to 1473 K) and different strain rates (from 0.01 to 5 s−1) were used in the hot compression tests. The specimens were firstly preheated to 1473 K at a rate of 5 K/s and held for 5 min at isothermal conditions for fully austenitizing. Then the specimens were cooling to the deformation temperature at the rate of 5 K/s and held for 30 s for the homogenization of temperature and compressed at the selected conditions. Finally, the compressed specimens were rapidly quenched in water.


### 2.2. Experimental Results and Analysis

The true strain-stress curves of as-cast P91 alloy steel under different deformation conditions are depicted in Figure3. It can be obtained that the flow stress curves are sensitively dependent on the strain rate and the temperature. It increases with the increase of strain rate and decreases with the increase of temperature. The dynamic recrystallization (DRX) occurs under higher temperature and lower strain rate.Figure 3
The true stress-strain curves of different deformation temperatures at strain rate of 0.1 s−1 (a); different strain rates at deformation temperature of 1473 K (b).
(a)
(b)During the hot deformation process, the flow stress is sensitive to the work hardening and the dynamic softening [10]. At the beginning of the deformation, the flow stress increases with the increasing strain. The material is under an unstable status due to the increasing dislocation density. The stored energy turns into a driving force for the dislocation migration. When the strain reaches to the critical strain, the dynamic recovery occurs and the flow stress increases owing to the confliction between softening and hardening. Moreover, with the strain continuing to increase, the work hardening and softening reach to a certain balance and the flow stress reaches to a plateau as well. Thus, the critical strain and work hardening are the significant parameters for the flow stress and should be considered in the constitutive models.


## 2.1. Materials and Methods

The chemical composition of as-cast P91 alloy steel was (wt.%) 0.084C -0.28Mn -0.38P -0.005S -0.28Si -8.63Cr -0.96Mo -0.209V -0.084Nb -0.21Ni -0.031N -(bal) Fe. The hot compression tests were carried out on a Gleeble-3500 thermomechanical simulator. Cylindrical specimens were machined with 8 mm in diameter and 12 mm in height.The specimens were deformed at the reduction of 60% in height. Different deformation temperatures (from 1173 to 1473 K) and different strain rates (from 0.01 to 5 s−1) were used in the hot compression tests. The specimens were firstly preheated to 1473 K at a rate of 5 K/s and held for 5 min at isothermal conditions for fully austenitizing. Then the specimens were cooling to the deformation temperature at the rate of 5 K/s and held for 30 s for the homogenization of temperature and compressed at the selected conditions. Finally, the compressed specimens were rapidly quenched in water.


## 2.2. Experimental Results and Analysis

The true strain-stress curves of as-cast P91 alloy steel under different deformation conditions are depicted in Figure3. It can be obtained that the flow stress curves are sensitively dependent on the strain rate and the temperature. It increases with the increase of strain rate and decreases with the increase of temperature. The dynamic recrystallization (DRX) occurs under higher temperature and lower strain rate.Figure 3
The true stress-strain curves of different deformation temperatures at strain rate of 0.1 s−1 (a); different strain rates at deformation temperature of 1473 K (b).
(a)
(b)During the hot deformation process, the flow stress is sensitive to the work hardening and the dynamic softening [10]. At the beginning of the deformation, the flow stress increases with the increasing strain. The material is under an unstable status due to the increasing dislocation density. The stored energy turns into a driving force for the dislocation migration. When the strain reaches to the critical strain, the dynamic recovery occurs and the flow stress increases owing to the confliction between softening and hardening. Moreover, with the strain continuing to increase, the work hardening and softening reach to a certain balance and the flow stress reaches to a plateau as well. Thus, the critical strain and work hardening are the significant parameters for the flow stress and should be considered in the constitutive models.


## 3. Constitutive Model for As-Cast P91 Alloy Steel

The experimental true stress-strain curves were fitted with the logarithmic from Hollomon equation (σ=s·εp), where σ is the flow stress, s is the strength coefficient, ε is the plastic strain, and p is the work hardening exponent. The work hardening rate θ (dσ/dε) in this investigation can be expressed from the experimental researches as [12, 13](1)θi=dσdε=σi+1-σi-1εi+1-εi-1=lσ+k,where l is a parameter which effected by the deformation temperature and strain rate and k is a constant. The relationships between work hardening rate and flow stress can be shown in Figure 4 (the deformation temperature was 1473 K), σc is the critical stress, σp is the peak stress, and σs is the saturation stress.Figure 4
The relationships between work hardening rate and flow stress.It can be obtained in Figure4, by extending the imaginary line with X-axis, and the intersection (where θ=0, σ=σs) was gotten. From (1), constant k can be found as k=-lσs. In Figure 4, the flow stress has a linear relationship between work hardening rate and flow stress before critical stress, and l is the slope of the line, so the point of critical strain and stress is the inflection point of the curves. From (1), the flow stress can be expressed as(2)σ=exp⁡l·ε·C-kl.With the condition ofσ=σc and ε=εc, the constant C can be described as(3)C=σc-σsexp⁡lεc.Substituting (3) into (2) the constitutive model of flow stress can be obtained as(4)σ=exp⁡l·ε-εcσc-σs+σs.The parameters in this model will be established and determined in the following parts. The values ofσc, σp, and σs under different deformation conditions can be obtained and illustrated in Table 1.Table 1
The values ofσc, εc, σp, and σs under different deformation conditions.
Strain rate (s−1)
T (K)
σ
c (MPa)
ε
c
σ
p (MPa)
σ
s (MPa)
0.01
1473
34.65
0.15
34.87
41.30
1423
49.35
0.16
50.14
55.01
1373
64.63
0.22
65.33
68.57
1323
88.2
0.37
88.79
104.84
0.05
1473
57.01
0.21
57.13
65.23
1423
70.2
0.25
71.10
76.54
1373
85.45
0.30
86.80
92.87
1323
114.14
0.34
115.15
127.31
0.1
1473
62.44
0.22
62.97
78.30
1423
76.01
0.26
77.38
81.74
1373
96.05
0.31
96.70
107.84
0.5
1473
83.28
0.27
84.07
98.14
1423
89.05
0.37
90.45
113.35
1
1473
87.25
0.31
88.84
104.57
1423
104.31
0.41
105.61
120.31


### 3.1. Calculations of the Activation Energy andZ Parameter

Arrhenius equations were represented [14] which are used in describing the relationship by hyperbolic sine, and it can be expressed with the Zener-Hollomom (Z) parameter as follows:(5)Z=ε˙exp⁡QRT=Asinh⁡ασn,where Q is the activation energy (kJ/mol), R is the constant (8.314 J·mol−1·K−1), T is deformation temperature (K), and A and α are material constants. According to (5), the strain rate can be expressed as(6)ε˙=A·Fσexp⁡-QRT,where F(σ) is the function of the stress which can be described in different conditions [15]; the functional relations between flow stress and strain rate can be described as follows:(7)Fσ=σn1ασ<0.8Fσ=exp⁡βσασ<1.2,where β and n1 are the materials dependent constants. Substituting the value of σp into (5)~(7) and taking natural logarithm of (7), the relationships between ln⁡ε˙ and ln⁡σ, ln⁡ε˙ and σ are illustrated in Figures 5 and 6, respectively. n1 (n1=8.78215) and β (β=0.06965) are determined from the slopes of straight line in Figures 5 and 6, where α, β, and n1 are related by α=β/n1, and α can be obtained as 0.00793.Figure 5
Relationship betweenln⁡σ and ln⁡ε˙.Figure 6
Relationship betweenσ and ln⁡ε˙.Taking natural logarithm of (6), activation energy Q can be expressed as(8)Q=Rn∂ln⁡sinh⁡ασ∂1/T.Similarly, taking natural logarithm on both sides of (5), the relationship between ln⁡sinh⁡ασ and ln⁡ε˙ can be obtained. Using the linear plot of ln⁡sinh⁡ασ versus ln⁡ε˙ (Figure 7), the value of n can be calculated (n=6.2154). The value of Q can be obtained by using the same method (in Figure 8, the linear fit between ln⁡sinh⁡ασ and 1000/T); the average value of Q is 328.1 kJ/mol.Figure 7
Relationship betweenln⁡ε˙ and ln⁡[sinh⁡(ασ)].Figure 8
Relationship between1000/T and ln⁡[sinh⁡(ασ)].


### 3.2. Calculation of Parameterl

According to (4), the flow stress can be expressed as(9)ln⁡σ-σsσc-σs=l·ε-εc.The determination of parameterl can be obtained from (9). The linear plot of ln⁡[(σ-σs)/(σc-σs)] versus (ε-εc) can be seen in Figure 9, and l is determined as the slopes of the lines. As mentioned above,l is sensitive to strain rate and deformation temperature which is similar with parameter Z. So establish a relationship between Z and l as(10)l=a·exp⁡bZ+c,where a, b, and c are the material parameters. The values of Z can be obtained in (6), and plotting Z and l in Figure 10. There has an exponential relationship between Z and l in Figure 10. Using the nonlinear regression in (10), all the three parameters can be obtained. Parameters a, b, and c are calculated to be 16.019, −0.107, and 2.692, respectively. Thus, (10) becomes(11)l=16.019·exp⁡-0.107Z+2.692.Figure 9
Relationship between(ε-εc) and ln⁡[(σ-σs)/(σc-σs)] at 1473 K (a); 1423 K (b).
(a)
(b)Figure 10
Relationship betweenl and Z.


### 3.3. Calculations ofεc and σc

The critical stressσc and critical strain εc cannot be directly obtained from the experiment data. Kim et al. [16–18] established the relationships between εc, σc, and Z under hot deformation, and it can be expressed as(12)εc=a1ln⁡Z+b1σc=a2ln⁡Z+b2.The values ofa1, b1, a2, and b2 can be obtained by using linear regression as shown in Figure 11. The values of a1, b1, a2, and b2 are 0.035, −1.03, 10.55, and −312.61, respectively. Hence, the critical strain and stress can be represented as(13)εc=0.035ln⁡Z-1.03σc=10.55ln⁡Z-312.61.Figure 11
The relationships betweenεc and ln⁡Z (a); σc and ln⁡Z (b).
(a)
(b)


### 3.4. Calculation ofσs

Saturated stressσs cannot be directly obtained from the strain-stress curves as critical stress. It can be obtained according to Figure 4. From (5), the mathematical model of saturated stress is shown in (14):(14)Z=A2sinh⁡ασsn2.Taking the nature logarithm on both sides of (14) gives(15)ln⁡Z=ln⁡A2+n2ln⁡sinh⁡ασs.The relationship ofσs and Z can be expressed as(16)σs=1αsinh-1ZA21/n2.The value of constantsn2 and A2 can be obtained using linear regression as presented in Figure 12. In summary, the model for σs can be represented as follows(17)σs=126.103·sinh-15.63×10-16Z6.25.Figure 12
Relationship betweenln⁡Z and ln⁡[sinh⁡(ασs)].


### 3.5. Verification of the Proposed Constitutive Model

The mathematical model for predicting flow stress of as-cast P91 alloy steel during the hot deformation can be summarized as(18)σ=exp⁡l·ε-εcσc-σs+σsl=16.019·exp⁡-0.107Z+2.692εc=0.035ln⁡Z-1.03σc=10.55ln⁡Z-312.61σs=126.103·sinh-15.63×10-16Z6.25Z=ε˙exp⁡328100RT.Figure13 shows the comparisons of the experimental and predicted flow stress by the proposed constitutive model. It is undeniable that there have deviations between the experimental data and the calculations of the value fitting. So taking the correlation coefficient R to study the accuracy of the calculated data, R can be obtained by the following equation:(19)R=∑i=1i=Nσexpi-σ-expσPi-σ-P∑i=1i=Nσexpi-σ-exp2∑i=1i=NσPi-σ-P2,where σexpi is the experimental data, σPi is the predicted data, σ-exp is the experimental average data, and σ-P is the average predicted data. The correlation coefficient R under different deformation conditions is illustrated in Table 2. The smaller the distinctness between R and 1, the lesser the error of the constitutive model. From Table 2, the result indicated that the proposed model can give an accurate and precise estimate of flow stress for as-cast P91 alloy steel.Table 2
The correlation coefficients under different deformation conditions.
T (K)
Strain rate (s−1)
0.01
0.05
0.1
0.5
1
5
1173
0.813
0.865
0.981
0.946
0.815
0.823
1273
0.836
0.872
0.869
0.968
0.965
0.857
1323
0.865
0.886
0.846
0.864
0.945
0.978
1370
0.908
0.928
0.947
0.885
0.907
0.958
1423
0.904
0.921
0.927
0.912
0.934
0.879
1473
0.876
0.875
0.978
0.897
0.982
0.887Figure 13
Comparisons between predicted and experiment flow stress by the proposed constitutive model under different deformation temperatures at strain rate of 0.01 s−1 (a); 0.05 s−1 (b); 0.1 s−1 (c); 1 s−1 (d).
(a)
(b)
(c)
(d)


## 3.1. Calculations of the Activation Energy andZ Parameter

Arrhenius equations were represented [14] which are used in describing the relationship by hyperbolic sine, and it can be expressed with the Zener-Hollomom (Z) parameter as follows:(5)Z=ε˙exp⁡QRT=Asinh⁡ασn,where Q is the activation energy (kJ/mol), R is the constant (8.314 J·mol−1·K−1), T is deformation temperature (K), and A and α are material constants. According to (5), the strain rate can be expressed as(6)ε˙=A·Fσexp⁡-QRT,where F(σ) is the function of the stress which can be described in different conditions [15]; the functional relations between flow stress and strain rate can be described as follows:(7)Fσ=σn1ασ<0.8Fσ=exp⁡βσασ<1.2,where β and n1 are the materials dependent constants. Substituting the value of σp into (5)~(7) and taking natural logarithm of (7), the relationships between ln⁡ε˙ and ln⁡σ, ln⁡ε˙ and σ are illustrated in Figures 5 and 6, respectively. n1 (n1=8.78215) and β (β=0.06965) are determined from the slopes of straight line in Figures 5 and 6, where α, β, and n1 are related by α=β/n1, and α can be obtained as 0.00793.Figure 5
Relationship betweenln⁡σ and ln⁡ε˙.Figure 6
Relationship betweenσ and ln⁡ε˙.Taking natural logarithm of (6), activation energy Q can be expressed as(8)Q=Rn∂ln⁡sinh⁡ασ∂1/T.Similarly, taking natural logarithm on both sides of (5), the relationship between ln⁡sinh⁡ασ and ln⁡ε˙ can be obtained. Using the linear plot of ln⁡sinh⁡ασ versus ln⁡ε˙ (Figure 7), the value of n can be calculated (n=6.2154). The value of Q can be obtained by using the same method (in Figure 8, the linear fit between ln⁡sinh⁡ασ and 1000/T); the average value of Q is 328.1 kJ/mol.Figure 7
Relationship betweenln⁡ε˙ and ln⁡[sinh⁡(ασ)].Figure 8
Relationship between1000/T and ln⁡[sinh⁡(ασ)].


## 3.2. Calculation of Parameterl

According to (4), the flow stress can be expressed as(9)ln⁡σ-σsσc-σs=l·ε-εc.The determination of parameterl can be obtained from (9). The linear plot of ln⁡[(σ-σs)/(σc-σs)] versus (ε-εc) can be seen in Figure 9, and l is determined as the slopes of the lines. As mentioned above,l is sensitive to strain rate and deformation temperature which is similar with parameter Z. So establish a relationship between Z and l as(10)l=a·exp⁡bZ+c,where a, b, and c are the material parameters. The values of Z can be obtained in (6), and plotting Z and l in Figure 10. There has an exponential relationship between Z and l in Figure 10. Using the nonlinear regression in (10), all the three parameters can be obtained. Parameters a, b, and c are calculated to be 16.019, −0.107, and 2.692, respectively. Thus, (10) becomes(11)l=16.019·exp⁡-0.107Z+2.692.Figure 9
Relationship between(ε-εc) and ln⁡[(σ-σs)/(σc-σs)] at 1473 K (a); 1423 K (b).
(a)
(b)Figure 10
Relationship betweenl and Z.


## 3.3. Calculations ofεc and σc

The critical stressσc and critical strain εc cannot be directly obtained from the experiment data. Kim et al. [16–18] established the relationships between εc, σc, and Z under hot deformation, and it can be expressed as(12)εc=a1ln⁡Z+b1σc=a2ln⁡Z+b2.The values ofa1, b1, a2, and b2 can be obtained by using linear regression as shown in Figure 11. The values of a1, b1, a2, and b2 are 0.035, −1.03, 10.55, and −312.61, respectively. Hence, the critical strain and stress can be represented as(13)εc=0.035ln⁡Z-1.03σc=10.55ln⁡Z-312.61.Figure 11
The relationships betweenεc and ln⁡Z (a); σc and ln⁡Z (b).
(a)
(b)


## 3.4. Calculation ofσs

Saturated stressσs cannot be directly obtained from the strain-stress curves as critical stress. It can be obtained according to Figure 4. From (5), the mathematical model of saturated stress is shown in (14):(14)Z=A2sinh⁡ασsn2.Taking the nature logarithm on both sides of (14) gives(15)ln⁡Z=ln⁡A2+n2ln⁡sinh⁡ασs.The relationship ofσs and Z can be expressed as(16)σs=1αsinh-1ZA21/n2.The value of constantsn2 and A2 can be obtained using linear regression as presented in Figure 12. In summary, the model for σs can be represented as follows(17)σs=126.103·sinh-15.63×10-16Z6.25.Figure 12
Relationship betweenln⁡Z and ln⁡[sinh⁡(ασs)].


## 3.5. Verification of the Proposed Constitutive Model

The mathematical model for predicting flow stress of as-cast P91 alloy steel during the hot deformation can be summarized as(18)σ=exp⁡l·ε-εcσc-σs+σsl=16.019·exp⁡-0.107Z+2.692εc=0.035ln⁡Z-1.03σc=10.55ln⁡Z-312.61σs=126.103·sinh-15.63×10-16Z6.25Z=ε˙exp⁡328100RT.Figure13 shows the comparisons of the experimental and predicted flow stress by the proposed constitutive model. It is undeniable that there have deviations between the experimental data and the calculations of the value fitting. So taking the correlation coefficient R to study the accuracy of the calculated data, R can be obtained by the following equation:(19)R=∑i=1i=Nσexpi-σ-expσPi-σ-P∑i=1i=Nσexpi-σ-exp2∑i=1i=NσPi-σ-P2,where σexpi is the experimental data, σPi is the predicted data, σ-exp is the experimental average data, and σ-P is the average predicted data. The correlation coefficient R under different deformation conditions is illustrated in Table 2. The smaller the distinctness between R and 1, the lesser the error of the constitutive model. From Table 2, the result indicated that the proposed model can give an accurate and precise estimate of flow stress for as-cast P91 alloy steel.Table 2
The correlation coefficients under different deformation conditions.
T (K)
Strain rate (s−1)
0.01
0.05
0.1
0.5
1
5
1173
0.813
0.865
0.981
0.946
0.815
0.823
1273
0.836
0.872
0.869
0.968
0.965
0.857
1323
0.865
0.886
0.846
0.864
0.945
0.978
1370
0.908
0.928
0.947
0.885
0.907
0.958
1423
0.904
0.921
0.927
0.912
0.934
0.879
1473
0.876
0.875
0.978
0.897
0.982
0.887Figure 13
Comparisons between predicted and experiment flow stress by the proposed constitutive model under different deformation temperatures at strain rate of 0.01 s−1 (a); 0.05 s−1 (b); 0.1 s−1 (c); 1 s−1 (d).
(a)
(b)
(c)
(d)


## 4. Utilization in the Research on the Hot Extrusion of CHE Process

The hot extrusion is a process with the characteristics of high strain rate and large deformation degree, and the deformation conditions are more complicated than the hot compression tests. It is vital to study the flow stress during the hot extrusion process. Thus, the application of the as-cast P91 alloy steel constitutive model on the CHE hot extrusion process should be further studied and explored. The finite element (FE) method is the efficient and economical alternative method for studying the flow stress behavior during hot deformation [19–21]. Hence, the constitutive model of the as-cast p91 alloy steel was encoded into the DEFORM platform by the FORTRAN language. The hot extrusion numerical simulation results were compared with the industrial experimental research.


### 4.1. The Thermal-Mechanical FE Model

According to the geometric symmetry of the billet and extrusion pipe, the 3D axisymmetric simplified FE models were established in Figure14. The thermal-mechanical FE model combined with the constitutive model of as-cast P91 alloy steel was encoded into DEFORM platform. During the simulation, the as-cast billet moves with the movement of terrace die while the die and needle keep static. The material of the billet was as-cast P91 steel, and the die was H-13 die steel. The required simulation conditions for the reliability evaluation of the thermal-mechanical coupled FE model are shown in Table 3.Table 3
Hot extrusion parameters setting for the thermal-mechanical coupled FE model.
Parameter
Content
Billet size (mm)
ϕ
1253
×
ϕ
500
×
1000
Initial billet temperature (K)
1533
Initial die temperature (K)
573
Extrusion ratio
9
Extrusion speed (mm/s)
50
Die taper angle (°)
42
Coulomb friction coefficient
0.04
Thermal emissivity of blank
0.7
Environment temperature (K)
293
Heat convection coefficient between billet and extrusion tolls (W·m−2·K−1)
1100
Heat convection coefficient between billet and air (W·m−2·K−1)
20Figure 14
The geometric model of the hot extrusion process. 1: die; 2: as-cast billet; 3: extrusion needle; 4: terrace die.


### 4.2. Industrial Hot Extrusion Experiment on Heavy Caliber Thick-Wall Pipe

The hot extrusion experiment was carried out by the 500 MN vertical hot extrusion equipment (Figure15) in the factory. The hot extrusion process was investigated by using the same extrusion parameters (Table 3) as FE analysis. The billet was presprayed with glass lubricant. The 80 mm in height glass lubrication gasket was placed between the billet and the terrace die. During the extrusion process, the load-displacement data of the terrace die were measured.Figure 15
The structural diagram (a); 500 MN vertical hot extrusion machine (b).
(a)
(b)


### 4.3. The Comparisons between FE Analysis and Hot Extrusion Experiment

Figure16(a) shows the experimental billet; Figure 16(b) is the simulated initial temperature distribution of the billet. It can be seen that the billet temperature decreases due to the heat convection and heat emissivity. The temperature of the surfaces is lower than the inner. Figure 17 shows the extrusion production pipe of this experiment, and the dimensional accuracy of the pipe meets the requirements of the industrial.Figure 16
Experimental billet (a); billet temperature distribution before extrusion (b).
(a)
(b)Figure 17
The extruded pipe under the extrusion process (a); after the extrusion (b).
(a)
(b)The distributions of temperature field, effective strain, and stress during the hot extrusion process are illustrated in Figure18. It can be found that the heat is generated during the extrusion process, and the temperature increases in the large deformation area (the conical surface and the die orifice). The temperature increases due to the transfer of the plastic deformation heat and the dynamic friction heat. The effective strain and stress field during the hot extrusion process indicated that the extrusion force increased during the process. When the material flows out of the die orifice, the effective strain remains stable while the effective stress reaches to the peak value. Furthermore, it also can be seen in Figure 18(c) that the values of the effective stress are larger than zero; this illustrates that the stress state during the hot extrusion is the compressive stresses rather than the tensile stresses. It shows that the hot cracking will not appear on the extruded pipe which is verified by the industrial experiment shown in Figure 17.Figure 18
Simulation results of hot extrusion at temperature field (a); effective strain field (b); effective stress field (c).
(a)
(b)
(c)The comparisons of the load-displacement curves between the measured and simulated results are illustrated in Figure19. The results showed that the maximum load value of the experiment is 430 MN, while the simulation value is 410 MN. The difference was due to the simplification of the FE analysis; the actual industrial extrusion conditions are more complicated. Thus the error is in a reasonable deviation range. The computed results present a good estimate of the relationship between load and displacement during the hot extrusion.Figure 19
The comparisons between measured and simulated displacement—load curves.The simulation and experimental results indicated that the proposed constitutive model of as-cast P91 alloy steel can give an accurate and precise estimate of the flow stress and it can be used to the analyze of the CHE process. In other words, the constitutive models established from the compressive tests can be applied to the industrial research and simulation.


## 4.1. The Thermal-Mechanical FE Model

According to the geometric symmetry of the billet and extrusion pipe, the 3D axisymmetric simplified FE models were established in Figure14. The thermal-mechanical FE model combined with the constitutive model of as-cast P91 alloy steel was encoded into DEFORM platform. During the simulation, the as-cast billet moves with the movement of terrace die while the die and needle keep static. The material of the billet was as-cast P91 steel, and the die was H-13 die steel. The required simulation conditions for the reliability evaluation of the thermal-mechanical coupled FE model are shown in Table 3.Table 3
Hot extrusion parameters setting for the thermal-mechanical coupled FE model.
Parameter
Content
Billet size (mm)
ϕ
1253
×
ϕ
500
×
1000
Initial billet temperature (K)
1533
Initial die temperature (K)
573
Extrusion ratio
9
Extrusion speed (mm/s)
50
Die taper angle (°)
42
Coulomb friction coefficient
0.04
Thermal emissivity of blank
0.7
Environment temperature (K)
293
Heat convection coefficient between billet and extrusion tolls (W·m−2·K−1)
1100
Heat convection coefficient between billet and air (W·m−2·K−1)
20Figure 14
The geometric model of the hot extrusion process. 1: die; 2: as-cast billet; 3: extrusion needle; 4: terrace die.


## 4.2. Industrial Hot Extrusion Experiment on Heavy Caliber Thick-Wall Pipe

The hot extrusion experiment was carried out by the 500 MN vertical hot extrusion equipment (Figure15) in the factory. The hot extrusion process was investigated by using the same extrusion parameters (Table 3) as FE analysis. The billet was presprayed with glass lubricant. The 80 mm in height glass lubrication gasket was placed between the billet and the terrace die. During the extrusion process, the load-displacement data of the terrace die were measured.Figure 15
The structural diagram (a); 500 MN vertical hot extrusion machine (b).
(a)
(b)


## 4.3. The Comparisons between FE Analysis and Hot Extrusion Experiment

Figure16(a) shows the experimental billet; Figure 16(b) is the simulated initial temperature distribution of the billet. It can be seen that the billet temperature decreases due to the heat convection and heat emissivity. The temperature of the surfaces is lower than the inner. Figure 17 shows the extrusion production pipe of this experiment, and the dimensional accuracy of the pipe meets the requirements of the industrial.Figure 16
Experimental billet (a); billet temperature distribution before extrusion (b).
(a)
(b)Figure 17
The extruded pipe under the extrusion process (a); after the extrusion (b).
(a)
(b)The distributions of temperature field, effective strain, and stress during the hot extrusion process are illustrated in Figure18. It can be found that the heat is generated during the extrusion process, and the temperature increases in the large deformation area (the conical surface and the die orifice). The temperature increases due to the transfer of the plastic deformation heat and the dynamic friction heat. The effective strain and stress field during the hot extrusion process indicated that the extrusion force increased during the process. When the material flows out of the die orifice, the effective strain remains stable while the effective stress reaches to the peak value. Furthermore, it also can be seen in Figure 18(c) that the values of the effective stress are larger than zero; this illustrates that the stress state during the hot extrusion is the compressive stresses rather than the tensile stresses. It shows that the hot cracking will not appear on the extruded pipe which is verified by the industrial experiment shown in Figure 17.Figure 18
Simulation results of hot extrusion at temperature field (a); effective strain field (b); effective stress field (c).
(a)
(b)
(c)The comparisons of the load-displacement curves between the measured and simulated results are illustrated in Figure19. The results showed that the maximum load value of the experiment is 430 MN, while the simulation value is 410 MN. The difference was due to the simplification of the FE analysis; the actual industrial extrusion conditions are more complicated. Thus the error is in a reasonable deviation range. The computed results present a good estimate of the relationship between load and displacement during the hot extrusion.Figure 19
The comparisons between measured and simulated displacement—load curves.The simulation and experimental results indicated that the proposed constitutive model of as-cast P91 alloy steel can give an accurate and precise estimate of the flow stress and it can be used to the analyze of the CHE process. In other words, the constitutive models established from the compressive tests can be applied to the industrial research and simulation.


## 5. Conclusions

The deformation characteristics of as-cast P91 alloy steel were investigated by means of the unidirectional hot compression experiments over a practical range of temperatures and strain rates. According to the results, the flow stress is sensitive to the work hardening and the critical strain. A parameterl was proposed to express the constitutive models. The thermal-mechanical finite element model was further developed by encoding the constitutive models into DEFORM platform. According to the actual thick-wall pipe, the hot extrusion process was simulated and experimented, respectively. Comparisons between the simulated and experimental results were carried out and confirmed that the developed constitutive models can be used in the industrial and simulation research of compact hot extrusion process.



---

*Source: 6582739-2017-07-13.xml*