Hot Deformation Behaviors and Intrinsic Hot Workability Map of Ti-12Mo-4Zr-5Sn Alloy Based on Physical Model and Polar Reciprocity Model

: The hot compression tests of Ti-12Mo-4Zr-5Sn alloy were tested on the thermo-mechanical simulator of Gleeble-3500 under isothermal and constant strain rate. We studied this alloy’s behavior during thermal deformation at the conditions of T = 670~820 °C,   = 0.001~10 s −1 , and deformation degree 70%. The rheological stress curves of the alloy were modified, the characteristics of the rheological stress curves were analyzed, and the activation energy map of hot deformation was established. A physical constitutive model of the alloy based on strain compensation was established, which has taken the relationship between Young’s modulus and self-diffusion coefficient and temperature into account. Moreover, the intrinsic hot workability ξ map of the alloy was established based on the polar reciprocity model. The results show that this alloy’s rheological stress will descend when the temperature of deformation rises and grow when the strain rate increases, and has negative sensitivity of temperature and positive sensitivity of strain rate. According to the error calculation, the physical constitutive model’s correlation coefficient is 0.9910 and the average relative error is 3.97%, which has good accuracy. Through the analysis of the microstructures of the instability zone and the stability zone, it was found that the instability mode of the instability zone was dominated by the local flow, and the deformation mechanism of the stability zone was dominated by the dynamic recrystallization. The optimum processing parameters of the alloy known from ξ map and metallographic structure are the following: T = 790~820 °C and   = 0.001~0.01 s −1 .


Introduction
Titanium alloys are widely used as biomedical metal materials for artificial joints, fracture fixation, spine implantations, soft tissue reconstruction, and other situations because of their good biocompatibility, stable physical and chemical properties, and elastic modulus similar to human bone [1][2][3][4]. Many scholars have reported the hot deformation behaviors of titanium alloys. Gao et al. [5] studied the hot deformation flow behavior and microstructure evolution of TA15 titanium alloy by using isothermal compression test and electron backscatter diffraction technique. They found that the flow stress of the inhomogeneous structure showed an obvious flow softening at low strain after reaching the peak stress, which is similar to that of the layered structure. Mandal et al. [6] carried out the thermal compression experiments on Ti-5Al-5Mo-5V-3Cr alloy at different temperatures and strain rates. The viscoplastic self-consistent model was used to simulate the mechanical response of the material. Ti-Mo-Zr-Sn alloys have good mechanical properties and belong to the third-generation biomedical titanium alloy. Among them, Ti-12Mo-4Zr-5Sn alloy has a strength that greater than 1000 MPa and an elastic modulus of 60~80 GPa. It is a metastable, low modulus, and high-strength medical titanium alloy material, which has a good practical application. For Ti-12Mo-4Zr-5Sn alloy, thermal processing is the most important processing method before obtaining the product, and also an important way of improving the internal microstructure. Therefore, it is necessary to study the thermal processing of this alloy.
The constitutive relation model of metal material is the premise of the process design and control of metal plastic forming. It can express the influence of deformation conditions (temperature, strain rate, and strain) on the rheological behavior, and is the basis of the finite element numerical simulation technology and the formulation and optimization of thermal deformation process parameters. In recent years, there have been many reports on the application of the constitutive model to exploring the thermal deformation behavior of titanium alloy [7,8], but most of the constitutive models are based on Arrhenius constitutive equation, which belongs to the macroscopic phenomenological constitutive model. Those cannot effectively reveal the physical properties of materials. The physical-based constitutive model effectively integrates the macro-variables with the physical variables to predict the rheological stress and can also reflect the physical properties of materials. However, there are many equations involved in the complex physical constitutive model, it is difficult to identify the material parameters, and the amount of computation of parameters is large [9,10]. A physical constitutive model based on creep theory has good applicability. The model takes into account the influences of temperature on Young's modulus and self-diffusion coefficient of materials and has achieved the prediction of the hot deformation behavior of materials by adopting the hyperbolic sine function. The plastic deformation of Ti-12Mo-4Zr-5Sn alloy at higher temperatures can not only prevent the cracking caused by deformation at high strain rate, but also obtain the required product shape. However, metal plastic deformation is more important to meet the requirements of microstructure and performance, so as to meet the actual production and application. Therefore, it is very important to select appropriate technological parameters to avoid various defects in the material forming process (such as instability, torsion, adiabatic shear, etc.). Murty et al. [11] discussed the relationship of intrinsic hot workability parameter between dynamic material model (DMM) and polar reciprocity model (PRM). However, most of the current researches use DMM to draw the processing map and get some appropriate process parameters, while there are relatively few reports about PRM. Compared with DMM, PRM not only considers the effect of strain rate on flow stress but also the effect of strain history on flow stress. It is suitable for optimizing the hot working parameters of the alloy. Although the effect of strain history on flow stress is not as significant as that of strain rate, it cannot be ignored. For example, the flow stress of material with insensitive strain rate is mainly affected by the strain history [12][13][14].
In this paper, the rheological properties of Ti-12Mo-4Zr-5Sn alloy at high temperatures were analyzed through hot compression experiments, and the rheological stress was modified by friction and temperature. Based on the modified data, the activation energies of hot deformation were calculated under different temperatures with different strains and strain rates, and the activation energy maps were constructed. A physical constitutive model based on strain compensation was constructed by introducing the strain into the physical constitutive equation. Finally, the intrinsic hot workability ξ map of the alloy was established by using the PRM theory. The microstructure of the corresponding region of the intrinsic hot workability ξ map was verified, and the processing parameters of the alloy were optimized. The results can provide theoretical guidance and data support for the actual hot forming of the alloy.
The hot compression of small cylindrical specimens with Φ 8 mm × 12 mm obtained from raw ingot was tested on the thermo-mechanical simulator of Gleeble-3500 (Dynamic Systems Inc., Portland, OR, USA). To minimize the friction between contact surfaces, thin graphite sheets were attached as a lubricant material to both ends of the specimen. Samples were heated to 670, 700, 730, 760, 790 and 820 °C at a speed of 5 °C /s, respectively, and held for 300 s. Thereafter, the compression tests were conducted at different strain rates (0.001, 0.01, 0.1, 1 and 10 s −1 ) with a compression rate of 70%. When the reduction in the height reached 70%, to retain the deformation microstructure, the specimen should be quenched immediately. The flow diagram of hot compression deformation is shown in Figure 1, which shows four stages: heating, holding, compression, and quenching. Finally, specimens for the compression testing were sectioned parallel to the deformation axis and the cut sides were used for metallographic examination after polishing. Microstructure characterizations were observed through FL7500 optical microscope (Jinan Metallographic Instrument Equipment Co., Jinan, China) after etching by the solution of 6%HNO3 + 4%HF + 90%H2O (vol.%).

Discussion of Flow Behaviors
In thermal compression experiments, the existence of friction will cause an obvious bulging phenomenon of the compressed sample, and the accumulation of deformation heat effects will increase the experiment temperature, which is the main source of stress error. However, these effects cannot be eliminated by the existing thermomechanical simulator, it is necessary to modify the experimental data by friction and temperature. The revised data are closest to the real hot deformation behavior of materials.
The equation derived from the finite element method can be used to modify the rheological stress measured in the experiment by friction: where pave is the measured stress, MPa; σf is the friction-corrected stress, MPa; Hi is the height of the specimen at any time, mm; Ri is the radius at any time, mm, it can be calculated by ; m and b are the friction factor and bulging parameters, respectively, and these equations can be obtained by , mm; RT is the radius of the contact surface between the end face and the indenter after compression, mm, it can be calculated by . ΔH is the height change of the specimen before and after compression, mm; R and H are the radius and height of the specimen after compression, mm.
The temperature rise caused by deformation heat effect can be expressed as Equation (4): where σ is the measured stress, MPa; ΔT is the temperature rise caused by material deformation, °C; p is the heat transfer coefficient related to strain rate; C is the heat capacity coefficient of the material, N·mm −2 ·°C −1 , and the heat capacity coefficient of titanium alloy is generally C = 4 [15].
Others are: The rheological stress change Δσ caused by temperature rise ΔT can be determined by the Equation (6): where Tn is the deformation temperature in the experiment, K.
According to Equations (1)~(6), the friction and temperature corrections were made for the rheological stresses of Ti-12Mo-4Zr-5Sn alloy measured at deformation temperatures of 670~820 °C and strain rates of 0.001~10 s −1 . The corrected rheological stress curves are shown in the dotted line of Figure 2. The effects of friction and temperature rise on the flow stress were interactive when the specimen was compressed. That is to say, the friction makes the experimental value higher than the true value, while the rise of temperature makes the experimental value lower than the true value.
According to the comparison of the rheological stress curves before and after friction and temperature correction in Figure 2, we found that the effects of friction and temperature rise on the flow stress become more and more obvious with the decrease in deformation temperature or the increase in strain rate. With the increase in strain, the trends between the revised curve and the measured curve present two situations: one is that the revised curve is below the original curve during the whole period, as shown in some curves in Figure 2a,b. The reason is that the true stress corrected by friction is less than the measured stress, and at a low strain rate, the main factor affecting the true stress is friction, while the effect of temperature rise has little effect on the true stress. Another case is that the revised curve lies below the original curve at first; with the accumulation of deformation, the dotted line will gradually approach the solid line first, and then deviate from the solid line gradually, shown in the curve at low temperature-high strain rate. This is because the effect of temperature rise on the true stress is less than that of friction when the degree of deformation is small. When the deformation accumulates to a certain level, the temperature gathered in the material will produce a temperature rise effect; the effect of temperature rise on the true stress will be higher than that of friction at this situation, and the corresponding correction curve will gradually deviate upward.  Figure 2 depicts a series of experimental rheological stress curves of Ti-12Mo-4Zr-5Sn alloy at different conditions, which have been revised through the friction and temperature corrections. It can be seen that the effects of strain rate and temperature of deformation on the rheological stress are more significant, the machining hardening phenomenon is very obvious at a high strain rate, and the steady state flow characteristics are presented at a low strain rate. When the temperature of deformation and strain rate are constant, as the strain increases, the curve has the characteristic of a typical dynamic recrystallization, that is, the rheological stress appears as a significant rheological softening after the peak, which is because the alloy has a low stacking fault energy and is prone to dynamic recrystallization under the influence of thermal deformation. Figure 3 shows the rheological stress of the alloy at strain 0.4, which can be seen from the figure that when the temperature or strain rate is constant, the stress increases with the increase in strain rate and decreases with the increase in deformation temperature, indicating that the alloy has negative temperature sensitivity and positive strain rate sensitivity. This is because when the temperature is constant, with the increase in strain rate, larger strain will be produced in a short period of time, and a large number of dislocations will proliferate. At the same time, the rheological stress increases with the increase in strain rate because there is not enough time to produce softening effect. When the strain rate is constant, the increase in temperature intensifies the movement of atoms, which makes the dislocation movement easier, thus promoting the nucleation and growth of dynamic recrystallization. At the same time, because the critical shear stress value of slip is reduced, dislocation density and rheological stress are reduced.

Determination of Hot Deformation Activation Energy
The metal material's high temperature deformation may include the following processes: the slip or climb of dislocations and the deformation caused by offset of opposite sign dislocations; the deformation caused by relative sliding between grain boundaries; the deformation caused by relative movement of interfaces; and the deformation caused by directional diffusion of interstitial atoms or vacancies. Each deformation process needs to overcome a certain energy barrier which is called the deformation activation energy of the process [16][17][18]. Taking the modified peak stress data of flow stress curve as an example, the hot deformation activation energy of Ti-12Mo-4Zr-5Sn alloy is calculated based on the Arrhenius hyperbolic sine function equation, as shown in Equation (7): where  is strain rate, s −1 ; σ is stress, MPa; Q is the activation energy of hot deformation, J/mol; R is gas constant, and R = 8.314 J/(mol·K); A, n, and α are material constants; T is temperature of deformation, K; the parameter α can be determined by the constants β and n1 in Arrhenius exponential equation and power function equation [19,20], as shown in Equation (8): The natural logarithms on both sides of Equation (7) can be written into the following two forms: And the equation of activation energy of hot deformation Q: Figure 4a According to Equations (9) and (10) Through linear fitting, we can get k = 6.038624, and n = 5.605, respectively. The average hot deformation activation energy obtained by introducing the calculated n and k values into Equation (11) is Q = 281.40 kJ/mol. )] [sinh( ln ln In the experiment, the material deformation degree was 70%. Similarly, the material parameters α, k, and n, with the modified strain between 0.1~1.2 (the interval is 0.1), were also calculated. The deformation activation energy Q at different temperatures and different strains was obtained, as shown in Table 1. According to the calculation results, the deformation activation energy is higher than the self-diffusion activation energies of pure α and pure β titanium alloys, 204 and 161 kJ/mol [21], respectively. It shows that the alloy's hot deformations under these conditions may correspond to different deformation mechanisms, which are controlled by processes other than high-temperature diffusion. To further analyze the physical mechanism of hot deformation of the alloy, the activation energy Q at different strain rates and different temperatures was calculated, as shown in Table 2. The activation energy maps of deformation were constructed according to Tables 1 and 2, respectively, as shown in Figure 6a  It can be seen from Figure 6 that the activation energy of deformation decreases with the increase in strain. It may be due to the increase in strain, the enhancement of slip or climb of dislocations, and the enhancement of dynamic softening effect in alloys, which decreases the activation energy required for material deformation [22]. The activation energy of the alloy at 670~730 °C is obviously higher than the self-diffusion activation energies of pure α and pure β titanium alloys. In this situation, the deformation may be mainly caused by dislocation and pinning mechanism, requiring larger energy to initiate dislocations. In the range of strain rates from 0.01 to 10s −1 , the required activation energy of deformation decreases with the increase in strain rate, which may be that the defect densities of dislocations and vacancies among grains increase accordingly with the increases in strain rate and the driving force required for dynamic softening. Thus, the energy consumed decreases for extricating the attraction from atoms during partial recrystallization of grains, resulting in the decrease in recrystallization activation energy.

Establishment and Analysis of Constitutive Model Based on Physics
The physical constitutive model studied in this paper is based on the hyperbolic sine function equation, which is similar to the Arrhenius hyperbolic sine constitutive equation. They were first proposed in creep applications [23] and later in hot working [10,[24][25][26]. The equations are as follows: where D(T) and E(T) represent the functional relationship between self-diffusion coefficient and Young's modulus and temperature, respectively. Among them, D0 is diffusion constant; Qsd is self-diffusion activation energy, J/mol; G0 and E0 represent shear modulus and Young's modulus at 300 K, GPa; Tm is melting point, K; α´, B´, and n´ are the material parameter to be determined.
The primary problem of constructing the above constitutive model with a physical basis is used to solve the material parameters in the equations of D(T) and E(T). Based on the literature [23], the material parameters of Ti-12Mo-4Zr-5Sn alloy are shown in Table 3.  By substituting the related material parameters in Table 3 into the Equations (12) and (13), the functions between the self-diffusion coefficient and Young's modulus and temperature of the alloy can be obtained: Peak stress reflects the maximum load-bearing capacity of materials under experimental conditions. Therefore, taking the peak stress of the modified rheological stress curve as an example, the process of solving material constants α´, n´, and B´ in the physical constitutive equation is Taking natural logarithms on both sides of Equations (17) and (18), and taking derivatives after that, it can be concluded that: Taking the natural logarithm on both sides of Equation (14), the Equation (21) is obtained: According to Equation (21) The physical constitutive equation of the alloy with peak stress data as substitution parameters can be obtained by combining the above-mentioned material constants n´, α´, and lnB´ and Equations (14)~ (16). The equation is shown in Equation (22) To sum up, the material constants in the physical constitutive model can be calculated by the steps in the flow chart, as shown in Figure 9. However, the physical constitutive model mentioned above only considers the peak stress but does not consider the effect of strain on the flow stress. Therefore, the strain factor was introduced into the physical-based constitutive model constructed in this paper for optimization. The rheological stress values under the strain conditions of 0.1~1.2 (the interval is 0.1, 12 strain variables in total) were selected to compute the material constants n´, α´, and lnB´ under different strains according to the above method. The results are shown in Table 4. The results were fitted by linear regression with a six-degree polynomial, and the functions between material constants and strain can be obtained, which can be used to construct the physical constitutive model with strain compensation.  Figure 10 shows the fitting results of six-degree polynomial between material constants n´, α´, and lnB´ and strain. It was found that the graphs that the linear correlation coefficient R of the fitting results are all above 0.99, which indicates that the data points in the graphs can be described by the functions obtained by the fitting of the six-degree polynomial. Medina et al. [27] discussed the general expression of the Zener-Hollomon parameter as a function of the chemical composition of low alloy and microalloyed steels through torsion tests. An expression between activation energy Q and the content of each alloying element is given and the equation is completed by determining the optimum values of α and n for all the steels. They demonstrated that A is a function of the activation energy. Hernandez et al. [28] developed a model to predict the austenite flow curves of low alloy and microalloyed steels. The parameters in model are a function of the dimensionless parameter Z/A of the equation of Sellars and Tegart, which has been modeled not only as a function of temperature and strain rate but also as a function of the chemical composition. The previous two methods are mainly based on the traditional Arrhenius-type hyperbolic sine function . Compared with them, the method used in this paper has the biggest difference that the influence of strain variables on material parameters in the physical constitutive model is considered. The results showed that these parameters (α, n, activation energy Q, and A) are material constants and are not dependent on deformation. The functions determined by Figure 10 are shown in Equation (23):  The relationship between material constants n´, α´, and lnB´and strain in Equation (23) is embedded into Equation (14), and the constitutive model with physical basis of the alloy optimized under compression deformation conditions can be obtained by combining Equations (15) and (16). This equation is shown in Equation (24): To verify the accuracy of the physical-based constitutive model for predicting flow stress, the predicted and experimental values are compared, and the results are shown in Figure 11a,b. It can be seen that the rheological stress predicted by the constitutive model with physical basis basically coincides with the experimental value and has certain reliability. Besides, this model's prediction accuracy was further evaluated by correlation coefficient (R) and average absolute relative error (AARE). The equations are as follows: (26) where Ci is the rheological stress's experimental data; Fi is the rheological stress's prediction data; C is the experimental data's average values; F is the prediction data's average values; n is the sample's total number. All the prediction data points and the corresponding experimental values were contrasted, as shown in Figure 12. We found that AARE between prediction stress and the experimental stress is 3.97%, R is 0.9910, and the data points of the constitutive model of which predicted value deviation is less than 10% account for 97.73%. Moreover, it can be obtained from Figure 12b that values of absolute error between 0 and 15 MPa account for approximately 92.42% of the total data. This indicates that the physical-based constitutive equation established by strain compensation is relatively good in accuracy and can accurately predict the high-temperature rheological stress of this alloy. The solving process of material parameters based on the physical constitutive model is not only simple and effective but also has a certain physical basis. It can be used as a new prediction model for the hot deformation rheological stress.

Construction and Analysis of Intrinsic Hot Workability ξ Map
The polar reciprocity model (PRM) is based on Hill correlation plastic flow rule. Based on the PRM, the intrinsic hot workability map can accurately reflect the evolution mechanism of the microstructure under different deformation conditions after the observation and verification of the microstructure. Thus, the safety zone and instability zone of the material in the thermal processing process are obtained, and the technological parameters (temperature, strain rate, strain variable, etc.) in the thermal processing process are optimized. Finally, the goal of controlling the evolution of microstructure and avoiding defects is achieved, and the microstructure and performance are obtained. In the polar reciprocity model, the total instantaneous power is divided into two parts, and then the intrinsic hot working parameters are defined. Figure 13 shows a convex potential surface Φ = K in stress space with origin O1. At point P on this surface, the radius vector O1P with arbitrary diameter and normal PN of Φ are made. The magnitude of normal PN is equal to KK´/Φ and KK´ is the dissipation power. A vector O1P´ that is parallel to the normal PN and equals to its size is made through O1, a point P on Φ can have a corresponding point P´ on Φ´. The trajectory of point P´ defines a potential surface Φ´ in strain rate space. K and K´ are the scalar functions of temperature and deformation history, which are affected by strain and deformation history. Under the given temperature and strain history conditions, Equation (27) can be obtained by the one-to-one correspondence between the strain rate space and the stress space: Two integral expressions on the left side of the equation represent the hardening power H W  and the dissipative power D W  , respectively. In the earlier dynamic material model (DMM), the total dissipative power is also divided into two parts, but the theories used in the two models to divide the total dissipative power are different. Thus, the intrinsic hot workability parameter is proposed. The equation is shown in Equation (28) where σ denotes stress, MPa; ε is strain;  is strain rate, s −1 ; H(ε P ) and C are two functions related to strain history. The equation of H(ε P ) is shown in Equation (30), which records the strain history from the beginning of plastic flow to strain ε1. σmin is the minimum flow stress from the beginning of plastic flow to strain ε = ε1. When the material presents ideal plastic flow, H(ε P ) = 0.
When there is a polar reciprocity between potential functions Φ and Φ´, the Equation (31) can be obtained: where m´ is a modified strain rate sensitivity coefficient. Equation (32) can be obtained from Equations (28)-(31): Raj and Kutumbarao [12] concluded the criterion of the instability of the polar reciprocity model by studying the hot working behavior of metals and basing on Hill correlation plastic flow rule: the material will lose its stability when ξ approaches 1.
To reduce the experimental error, Matlab was used to interpolate the flow stress data by the cubic spline function. The data can be substituted into Equation (30) to calculate the strain function value H(ε P ), Equation (31) to obtain the modified strain rate sensitivity coefficient m´ under different deformation conditions, and Equation (32) to calculate the value of ξ. The intrinsic hot workability ξ maps with true strains of 0.3, 0.6, 0.9 and 1.2 are drawn on the two-dimensional plane, as shown in Figure 14. According to the PRM criterion mentioned above, the region where ξ approaches 1 is the unstable region. To avoid structural instability, it is necessary to consider the instability of the material under different strains to ensure that the material has qualified properties after processing and forming. Figure 14 shows that the instability region of ξ map varies under different strains. The area of the rheological instability region at T = 670~820 °C and  =0.316~10 s −1 increases with the increase in strain and gradually moves toward the direction of low strain rate. In the deformation process, the processing window of the rheological instability zone of Ti-12Mo-4Zr-5Sn alloy can be characterized by the following parameters: T = 670~820 °C and  = 0.316~10 s −1 . In the region, the corresponding instability parameters of the material are relatively high, and instability phenomena such as local plastic flow or adiabatic shear band are prone to occur.
Among the four types of true strain ξ maps, there is mainly one ξ minimum region. The region is characterized by the following parameters: T = 790~820 °C and  = 0.001~0.01 s −1 . When the true strains are 0.3, 0.6, 0.9, and 1.2, the minimum values are 0.42, 0.44, 0.44, and 0.47, respectively, which is less affected by the strain. The basin area in ξ map corresponds to favorable deformation mechanisms such as dynamic recovery, dynamic recrystallization, and superplasticity; while the peak area corresponds to disadvantageous deformation mechanisms such as fracture, local flow, and adiabatic shear; and the grain growth leads to the increase in ξ value [13]. When determining the parameters of the hot working process, the range of parameters of favorable deformation should be given priority, and whether the range of parameters of the hot working process in a stable region should be considered. The ranges of favorable deformation mechanism in Figure 14 are the following: T = 790~820 °C and = 0.001~0.01 s −1 , which are the suitable processing parameters.

Verification of Microstructure
The validity of the PRM ξ map is verified by observing the microstructure. Figure 15 shows the microstructures of the deformation instability zone of the alloy. Figure 15a corresponds to deformation parameters: T = 670 °C and  = 1 s −1 ; Figure 15b corresponds to deformation parameters: T = 760 °C and = 10 s −1 ; Figure 15c corresponds to deformation parameters: T = 820 °C and  = 10 s −1 . Figure 15a-c all have local plastic flow phenomena, which is due to the low thermal conductivity of the alloy and the uneven deformation of different parts of the sample, resulting in obvious local temperature rise. That means that, during plastic deformation, the alloy produces a lot of heat, which cannot diffuse to the surrounding area, resulting in local rheological instability. Local plastic flow generally occurs at a high strain rate, forming a zigzag-shaped micro-band with an angle of 35~40° from the direction of the principal stress. Local plastic flow occurs in the deformed structure, and the uniformity of deformation and structure will be poor, which will affect the properties of the material after forming.  Figure 16 shows the microstructures of the deformation stable zone of the alloy. Figure 16a corresponds to deformation parameters: T = 790 °C and  = 0.001 s −1 ; Figure 16b corresponds to deformation parameters: T = 820 °C and  = 0.001 s −1 . Figure 16a's corresponding ξ value is 0.65, and its corresponding intrinsic hot workability ξ map is located near the basin of the PRM ξ map. It can be seen from Figure 16a that the dynamic recrystallization is more sufficient, the structure is uniform, the β grain is smaller, the size difference between grains is small, and the deformation is stable. Figure 16b's corresponding ξ value is 0.62, which is close to the minimum of ξ. The dynamic 100μm 100μm 100μm recrystallization under this deformation parameter is more sufficient and the structure is finer. The volume fraction of recrystallization increases. The average grain size is about 20 μm, and the structure is in a rheological stable state. It can also be seen from Figure 16 that the temperature has an important effect on the dynamic recrystallization of the alloy at a strain rate of 0.001 s −1 . Dynamic recrystallization is generally a continuous process in which the phase difference increases gradually through the accumulation of dislocations, and the small-angle grain boundary transforms gradually to the large-angle grain boundary, thus producing new grains [29]. In the process of deformation, the grain boundaries are serrated due to the local intergranular migration. At a low deformation rate, the small-angle grain boundaries have enough time to transform to the large-angle grain boundaries, and dynamic recrystallization easily occurs. Deformation temperature also has an effect on grain boundary migration, so when the temperature increases, atomic diffusion speed increases, grain boundary mobility increases, and dynamic recrystallization is more likely to occur.
Through the prediction of the intrinsic hot workability ξ map and the observation of the microstructure, the validation was carried out. The optimum processing parameters of the alloy are the following: T = 790~820 °C and  = 0.001~0.01 s −1 . In addition, dynamic recrystallization is the main deformation mechanism in the deformation stabilization region. Fine recrystallization is conducive to the subsequent processing.

Conclusions
When the temperature or strain rate is constant, the peak stress will rise when the strain rate grows and descend when the deformation temperature rises; Ti-12Mo-4Zr-5Sn alloy has negative sensitivity of temperature and positive sensitivity of strain rate.
1) The hot deformation activation energy of Ti-12Mo-4Zr-5Sn alloy is higher than that of pure α and pure β titanium alloys. The results show that the hot deformation of the alloy is controlled by a process other than high-temperature diffusion and may correspond to different deformation mechanisms.
2) A physical-based constitutive model considered the effect of deformation temperature on the self-diffusion coefficient and Young's modulus of the alloy was established, and the model was optimized. Strain compensation was also considered in the optimized model. After error calculation, this model's correlation coefficient is 0.9910, the average relative error is 3.97%, and the data points within the deviation of 10% accounting for 97.73%, respectively, which have good accuracy. The model is not only simple and effective but also has a certain physical basis. It can be used as a new prediction model for the hot deformation rheological stress.

3)
Based on the polar reciprocity model, the intrinsic hot workability ξ maps of Ti-12Mo-4Zr-5Sn alloy were drawn. The range of rheological instability zone of the alloy was obtained by observing the microstructures: T = 670~820 °C and  = 0.316~10 s −1 . The optimum processing parameters of the alloy are as follows: T = 790~820 °C and  = 0.001~0.01 s −1 . Local flow is the main instability form in the instability region, while dynamic recrystallization is the main deformation mechanism in the stability region.