Next Article in Journal
Diffusion Mechanism of Cinnamon Essential Oils Release from Calcium Alginate Based Controlled Release Films in Contact with Food Simulating Solvent
Next Article in Special Issue
Hot Deformation and Dynamic Recrystallisation Behaviour of Twin-Roll Cast Mg-6.8Y-2.5Zn-0.4Zr Magnesium Alloy
Previous Article in Journal
Precision Hard Turning of Ti6Al4V Using Polycrystalline Diamond Inserts: Surface Quality, Cutting Temperature and Productivity in Conventional and High-Speed Machining
Previous Article in Special Issue
Microstructure and Texture Evolution of Mg-Gd-Y-Zr Alloy during Reciprocating Upsetting-Extrusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Predictability of Microstructure Evolution during Hot Deformation of Titanium Alloys

by
Ricardo Henrique Buzolin
1,2,*,
Franz Miller Branco Ferraz
1,2,
Michael Lasnik
3,
Alfred Krumphals
3 and
Maria Cecilia Poletti
1,2
1
Christian Doppler Laboratory for Design of High-Performance Alloys by Thermomechanical Processing, Kopernikusgasse 24, 8010 Graz, Austria
2
Institute of Materials Science, Joining and Forming, Graz University of Technology, Kopernikusgasse 24, 8010 Graz, Austria
3
voestalpine BÖHLER Aerospace GmbH & Co KG, Mariazellerstraße 25, 8605 Kapfenberg, Austria
*
Author to whom correspondence should be addressed.
Materials 2020, 13(24), 5678; https://doi.org/10.3390/ma13245678
Submission received: 22 November 2020 / Revised: 7 December 2020 / Accepted: 10 December 2020 / Published: 12 December 2020
(This article belongs to the Special Issue Dynamic Recrystallization Behaviors in Metals and Alloys)

Abstract

:
Two different mesoscale models based on dislocation reactions are developed and applied to predict both the flow stress and the microstructure evolution during the hot deformation of titanium alloys. Three distinct populations of dislocations, named mobile, immobile, and wall dislocations, describe the microstructure, together with the crystal misorientation and the densities of boundaries. A simple model consisting of production and recovery terms for the evolution of dislocations is compared with a comprehensive model that describes the reactions between different type of dislocations. Constitutive equations connect the microstructure evolution with the flow stresses. Both models consider the formation of a high angle grain boundary by continuous dynamic recrystallization due to progressive lattice rotation. The wall dislocation density evolution is calculated as a result of the subgrain size and boundary misorientation distribution evolutions. The developed models are applied to two near-β titanium alloys, Ti-5553 and Ti-17, and validated for use in hot compression experiments. The differences in the predictability between the developed models are discussed for the flow stress, dislocation densities and microstructure evolutions. Only the comprehensive model can predict the different reactions and their contributions to the evolution of mobile and immobile dislocation densities. The comprehensive model also allows for correlating the elastic strain rate with the softening and hardening kinetics. Despite those differences, the selection of the model used has a small influence on the overall prediction of the subgrain size and the fraction of high angle grain boundaries.

Graphical Abstract

1. Introduction

As a high stacking fault energy (SFE) material, the β-phase in titanium and its alloys undergoes a complex microstructure evolution during hot deformation due to dynamic recovery (DRV) and dynamic recrystallization (DRX) [1,2,3,4,5,6,7,8]. Although there are some controversies [1,2,3,4], it is generally accepted that the β-phase undergoes DRV [5] followed by continuous dynamic recrystallization (cDRX) [8] or geometric dynamic recrystallization (gDRX) [6,7] at large strains. cDRX occurs by progressive lattice rotation [6,7] and was also observed in Ti-64 [5], while gDRX occurs at lower strain rates and lower temperatures [6,7].
From discrete dislocation dynamics [9] to phenomenological approaches [10,11], the hot deformation behavior of metallic alloys has been modeled and predicted. Among the most powerful methods to predict the hot deformation of complex parts, the mesoscale models are of great interest [12]. Kocks [13], Bergström [14], Bergström and Roberts [15], Roberts and Bergström [16] and Mecking and Estrin [17] made the initial attempts to develop dislocation density-based models for metal plasticity. Ghoniem et al. [18] proposed a dislocation-based model to describe the dislocation creep behavior of steels, accounting for the several dislocation density reactions that can occur on a mesoscopic scale. A similar dislocation density reaction approach was applied to describe the hot deformation of a Ti-5553 alloy [19] and a Ti-17 alloy [20]. Simpler dislocation density-based models, mostly based on Kock–Mecking formalism [21], were developed to predict the deformation behavior of other Ti alloys. The deformation of a Ti-6Al-4V over a wide temperature range was achieved in Ref. [22], while constitutive equations applied to room temperature deformation describe the complex deformation mechanism, as described in [23]. Picu and Majorell [24] proposed constitutive equations to predict the influence of the chemical composition on the deformation behavior of a Ti-6Al-4V. The Kocks–Mecking formalism was coupled with dDRX to simulate the flow softening observed in the α+β domain [25,26,27]. The approach was extended in Ref. [28] to include Hall–Petch describing the hot deformation behavior of Ti alloys in both α+β and β domains. Babu et al. [29] extended the modeling approach to include the effects of vacancy evolution and globularization of the α-phase. The damage was coupled to dDRX to predict the tensile properties of a TA15 alloy [30]. Although cDRX or gDRX are observed as restoration mechanisms for high stacking fault energy materials, most existing models still consider discontinuous dynamic recrystallization (dDRX) as the deformation mechanism of titanium alloys [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40], and also do not consider the deformation in the α+β domain as a biphase plastic deformation [25,26,27,29,30,31,34,35,37,38,39].
This publication aims to compare and discuss the predictability of the well-established Kock–Mecking formalism [21] with a novel model described in Refs. [19,20] based on the work of Ghoniem et al. [18] where different dislocation reactions can be quantified. Both dislocation rate models are combined with a cDRX model and can predict the evolution of any starting microstructure. Moreover, the plastic deformations of both α and β phases are considered in the α+β domain and a change in flow partitioning is related to the flow softening.

2. Materials and Methods

Cylindrical samples of 10 mm diameters and 15 mm lengths were cut perpendicular to the forging direction and parallel to the cogging direction from disks of 128 mm and 256 mm diameters from commercial Ti-5553 and Ti-17, respectively. The cogged Ti-5553 presented globular primary α grains, and β grains elongated in the cogging direction with a size ~1 mm in length and ~150 µm in width. The as-received microstructure of the Ti-17 consisted of a lamellar microstructure and prior globular-like β-grains of ~500 µm in length. The β-transus temperature is ~860 °C for Ti-5553 in agreement with the literature [41,42] and ~865 °C for Ti-17. These temperatures were estimated using JMatPro® v.10 for the chemical compositions shown in Table 1.
Hot compression tests were performed using a Gleeble® 3800. Experiments were performed in the α+β phase field at 800 °C, 820 °C and 840 °C for the Ti-5553 and at 810 °C, 830 °C, and 850 °C for the Ti-17. Hot compression of the single β-phase was performed at 880 °C, 900 °C and 920 °C for the Ti-5553, and at 930 °C, 950 °C and 970 °C for the Ti-17. Five different constant true strain rates from 0.001 s−1 to 10 s−1 were tested. The specimens were heated at 5 °Cs−1 up to the testing temperature and held for 15 min before deformation, except for the hot compression at the single β-phase-field for Ti-5553, where the holding time was 5 min. The temperature was measured with a thermocouple type K welded at the surface of the sample, and Ar was used as a protective atmosphere to prevent oxidation and α-case formation. Graphite and tantalum foils were used to minimize the friction of the samples with the anvils. The samples were in situ water-quenched immediately after hot deformation to preserve the deformed microstructure. Some samples were water-quenched immediately after annealing to investigate the microstructure before deformation. The stress-strain data were corrected with respect to adiabatic heating according to the procedure described elsewhere [43,44]. The material density was considered as 4300 kg·m3 [45], the specific heat ( C p in (J·K−1·kg−1)) and material conductivity ( K w in (W·m−1·K−1)) of the workpiece are given in Equations (1) and (2), respectively. A thermal conductivity of 10 (W·m−1·K−1) [44] was considered for a SiN die.
C p   =   5.82 · 10 9 · T 2 1.12 · 10 6 · T + 6.18 · 10 3 ,
K w   =   0.0109 · T + 34.22 ,
The samples annealed before deformation and hot compression were cut in half perpendicular to the radial direction, cold mounted, ground using SiC paper, polished using silica oxide polishing suspension (OPS) solution, and vibro-polished before the microstructure analysis. Preparation using a cross-section polisher JEOL SM-09010 (JEOL, Tokyo, Japan) was used for EBSD investigations carried out in Ti-17. A TESCAN Mira3 (TESCAN, Brno, Czech Republic) and a JEOL JSM7001F (JEOL, Tokyo, Japan) electron scanning microscope were used for microstructure investigation. A TSL-OIM Data Collector software package was used for performing the EBSD measurements. The measurements of (500 µm)2 with a step size of 0.5 µm to investigate the β-phase were carried out at an acceleration voltage of 30 kV and a spot size of 80 nm using a TESCAN Mira3 (TESCAN, Brno, Czech Republic) microscope. An acceleration voltage of 15 kV and a spot size of 5 nm was used for the measurements of the α and β phases in Ti-17 using a JEOL JSM7001F (JEOL, Tokyo, Japan) microscope. The OIM Analysis v.8 software was used to treat the data. The minimum grain size of 2 µm was chosen to guarantee that each grain consists of at least four measured pixels. Misorientation angles larger than 12° were assigned to high angle grain boundaries (HAGBs). The confidence index was standardized for each grain. The EBSD data were then cleaned, considering a minimum confidence index of 0.5 correlated to the neighbor grains. A maximum of ~5% of the data were cleaned.

3. Modeling Strategy

The developed model consists of four main elements: microstructure concept, constitutive equations, plastic strain rate partitioning, and rate equations. The model parameters values and expressions are given in Appendix A.

3.1. Microstructure Model

The microstructure was modeled based on its mean properties, i.e., an average subgrain and grain sizes, as well as three different populations of dislocations. The representative microstructure entity named subgrain is considered to be surrounded by low angle grain boundaries (LAGBs) and HAGBs, named subgrain. The LAGBs are considered to be sharp and consisted of n sets of dislocations forming a wall [46]. Two populations of dislocations are distinguished within the subgrain: the mobile and the immobile dislocations. All dislocation densities are related to the whole volume of material. The mobile dislocations slide along the active slip systems [47] and accommodate the plastic deformation. Obstacles, such as other dislocations, boundaries, or phase interfaces [48], block the immobile dislocations. They do not contribute directly to plastic deformation and are responsible for the hardening [49].
The misorientation distribution and the volumetric boundary density define, respectively, the fractions and size of the subgrain. In the case of a substantial fraction of LAGB, a subgrain is surrounded mainly by LAGB. In an ideally fully recrystallized material, the fraction of LAGB is negligible, and the subgrain size is coincident with the grain size. The full description of the microstructure is achieved by including the misorientation distribution. This is modeled as a sum of Rayleigh and Mackenzie [50] distributions [19]. The Rayleigh distribution ( Θ R ( θ ) ) accounts for the misorientation caused by the LAGB formed during deformation. The increase in boundary misorientation leads to the evolution of HAGB. The Mackenzie distribution ( Θ M ( θ ) ) considers a fully recrystallized material with random texture. Texture evolution is not considered in the developed models. The initialization of the microstructure is explained in Appendix B, and the microstructural features values for the Ti-5553 and Ti-17 are defined in Table A1 in Appendix A.

3.2. Constitutive Equations

Constitutive equations are required to couple the microstructure with the flow stress ( σ ). The thermal, σ t h , and athermal, σ a t h , stresses are considered as the constituents of the flow stress for each phase, Equation (3). Any effect of dislocation pile-up along the grain boundary (Hall–Petch effect) is considered negligible due to the small Hall–Petch coefficient of typical Ti alloys at high temperatures [28,51]. However, this effect must be considered for deformation at lower temperatures.
σ x   = σ a t h x + σ t h x ,         x   = α ,   β
The athermal stress is expressed according to Equation (4) as a function of the immobile, mobile, and wall dislocation densities ( ρ i , ρ m , and ρ w , respectively), the Taylor constant (α), the Taylor factor (M), and the shear modulus (G).
σ a t h x   =   α M G b ρ i x + ρ m x + F w x ρ w x ,         x   = α ,   β
The Taylor constant is considered as 0.1 for both alloys, and its low value is attributed to the lattice resistance due to the interaction between parallel dislocations as the main barrier for dislocation movement at high temperature [52]. F w is described in Appendix C. The thermal stress is calculated using the fitted yield stress ( σ Y S ), Equation (5). The procedure to obtain σ Y S is described in Appendix D.
σ t h x   =   σ Y S x σ a t h x 0 ,         x   = α ,   β
σ a t h 0 is calculated considering initial dislocation densities. The σ t h is constant for a constant temperature and strain rate since the σ Y S does not vary with the temperature and the strain rate.

3.3. Load Partitioning

The deformation rate of hard platelets of the α-phase is different than the deformation rate of the soft β-matrix. As explained in [53,54,55], iso-strain, iso-stress, and/or iso-power load conditions have to be considered. The change in load transfer between both phases explains the observed flow softening in the α+β domain. The Ti-5553 and the Ti-17 exhibit an initial globular and lamellar structure, respectively. Thus, it is proposed:
  • Ti-5553—the rotation and accommodation of plastic deformation in the α-particles up to a steady-state condition lead to the change from iso-power to iso-stress regimes;
  • Ti-17—the process of dynamic α-globularization leads to a change from iso-strain to iso-stress regimes.
The eventual formation of texture is not considered in the model. The change in flow partition mechanisms is here used to describe the combined effects that lead to flow softening during deformation in the α+β domain. The equations and assumptions of iso-strain, iso-stress, and/or iso-power load conditions are explained in Appendix E.

3.3.1. Overall Stress

The overall stress is calculated using a rule of mixtures, Equation (6) as a function of the fraction of α-phase F v , and the stress in the α ( σ α ) and β ( σ β ) phases.
σ   =   F v σ α +   ( 1 F v ) σ β

3.3.2. Fraction of Material in Iso-Stress Regime

In the case of the Ti-5553, the volumetric fraction of α and β phases ( f i s o σ ) that evolves from iso-power to iso-stress load partitioning is given empirically by Equation (7) as a function of the overall strain (ε) and strain rate ( ε ˙ ).
f i s o σ   =   1 exp ( A i s o σ ( ε ˙ ε ˙ r e f ) n s ε n s )
A i s o σ and n s are fitting parameters and ε ˙ r e f is a reference strain rate.
In the case of the Ti-17, the f i s o σ is the fraction of α-globularization ( f g l o b ). Finally, the overall strain rate for each phase ( ε ˙ x ) is calculated as a rule of mixtures of both contributions, Equation (8).
ε ˙ x = ( 1 f i s o σ ) ε ˙ x y + f i s o σ ε ˙ x σ ,         x = α ,   β ;   y = W   ( Ti 5553 )   or   ε   ( Ti 17 )

3.4. Dynamic α-Globularization

The dynamic α-globularization consists of the formation of boundaries within the α-platelets and the movement of α/β interfaces leading to the formation of globular α-particles [56,57]. It is proposed that initially the load is transferred to the α-platelets following the iso-strain regime. With the onset of the α-globularization, the strain rate in the α-phase reduces, reaching the iso-stress regime at the flow steady-state condition.
An α-platelet of mean aspect ratio ( A R ) with a mean thickness ( t α ) and mean width ( w α   =   t α · A R ) is considered. The dynamic α-globularization consists of the formation of boundaries within the α-platelet due to plastic deformation followed by the formation of new α/β interfaces by the migration of grooved interfaces with a velocity v g l o b . Thickening of the α-particles can occur, leading particles with an aspect ratio of one. The overall boundary density within the α-phase is a sum of the production of new boundaries modeled similarly to cDRX [46], minus the consumption of existing boundaries within the α-phase due to α/β interface movement through the HAGBs within the α-phase, Equation (9).
d S v α   d t   =   α C D R X b n θ 0 Δ ρ c D R X ε ˙ p α f B G f H A G B α S v α v g l o b t α
ε ˙ p α is the plastic strain rate in the α-phase, f B G is the globularized fraction, S v α is the boundary surface density within the α-platelet, S g l o b is of the α-phase, and v g l o b is the velocity for α/β interface movement. Equation (10) describes the rate of the production of α/β interface density formed during dynamic globularization ( d S g l o b / d t ).
d S g l o b   d t   =   f B G f H A G B α S v α v g l o b t α   | S g l o b   =   0 t ( d S g l o b   d t ) d t
Equation (11) describes the total boundary density formed within the α-platelet.
S T O T A L α   =   S v α + S g l o b
Full globularization of an α-platelet occurs when S v α = 0. Equation (12) gives the globularization fraction.
f g l o b   =   S g l o b S T O T A L α
Appendix F describes the dynamic α-globularization in more detail.

3.5. Rate Equations and Microstructure Evolution

The evolution of the internal variables is expressed as differential equations and is solved using numerical incremental steps for both α and β phases. In this manuscript, two models to describe the rate equations are compared: (1) a more general and simple model based on the work of Kocks and Mecking [21] (herein named “Model KM”); (2) a model with descriptions of dislocation reactions based on the work of Ghoniem et al. [18] (herein named “Model G”).

3.5.1. Mobile Dislocation Density

Model KM
The mobile dislocation density is considered to be constant at a constant strain rate, in agreement with the Orowan’s relationship (Equation (13)). The experimental strain rate ( ε ˙ ) and the plastic strain rate ( ε ˙ p ) are the same.
v g x   =   M x ε ˙ p x   b x ρ m x ,         x   = α ,   β
The mobile dislocation density is assumed invariable with respect to plastic strain rate and strain because:
  • the developed model does not consider any effect of texture and change in the Taylor factor (M);
  • the Burgers vector (b) is constant during deformation;
  • the glide velocity is independent of strain if the plastic strain rate and temperature are kept constant.
Thus, the ratio of the glide velocities for two different strain rates equals the ratio of the strain rates.
Model G
Equation (14) describes the rate for mobile dislocation density. The terms were initially proposed by [18]. The terms related to the production of dislocations at subgrain boundaries and the movement of HAGBs during cDRX are further incorporated into the model as described in Ref. [19].
d ρ m x d t   =   ρ m x 3 / 2 v g + δ S G x ρ m x Φ s g x v g x λ m x 2 ρ m x v g x Φ s g x 8 ρ m x v c m x λ m x δ D R V x ρ m x ( ρ m x + ρ i x ) v g x ρ m x S v x f H A G B x v H A G B x ,         x   = α ,   β
The production of mobile dislocations due to Read’s source and at subgrain boundaries corresponds to the first and second terms, respectively. The third, fourth and fifth terms correspond to the immobilization of mobile dislocations at subgrain boundaries, the static recovery due to climb, and the dynamic recovery, respectively. The last term is related to the consumption of mobile dislocations due to the movement of HAGBs during cDRX. v g is the glide velocity, λ m and λ i are the inter-dislocation distance for the mobile and immobile dislocations, respectively, v c m and v c i are the climb velocities for the mobile and immobile dislocations, respectively, and δ D R V is the critical distance of dislocation annihilation via DRV. If the production of mobile dislocations at subgrain boundaries is neglected, the immobilization at subgrain boundaries can lead to the annihilation of all mobile dislocation. For simplification, δ S G is an internal variable adjusted at each step to maintain a constant glide velocity during plastic deformation at a given strain rate and temperature. f H A G B is the fraction of HAGBs ( f H A G B   =   1 f L A G B ). The velocity of HAGBs ( v H A G B ) is described in Appendix G. The ε ˙ p is determined as a function of the mobile dislocation density by the Orowan’s relationship (Equation (13)), where ε ˙ p is the output. Here, the difference with respect to the “Model KM” is that the glide velocity is calculated phenomenologically according to Equation (15). The expression describes the mechanisms of strengthening [18] as a function of the activation energy for gliding ( W g ), the thermal stress, the atomic volume ( Ω ), the Boltzman’s contant ( k B ) and the temperature ( T ). The pre-factor ( a 1 ) is calculated as described in Ref. [19].
v g x   =   a 1 x exp ( W g x k B T ) σ t h x Ω k B T ,         x   = α ,   β
The experimental strain rate equals the sum of the plastic strain rate and the elastic strain rate ( ε ˙ e ). Consequently, the variation in the flow stress at each iteration is expressed according to the derivative of the Hook’s law, Equation (16), as a function of the Young modulus (E).
σ ˙ x   =   E x ( ε ˙ x ε ˙ p x ) ,         x   = α ,   β

3.5.2. Immobile Dislocation Density

Model KM
Equation (17) describes the immobile dislocation density rate using a modified expression from Kocks and Mecking [21] and Gourdet and Montheillet [46].
d ρ i x d t   =   h 1 x   ε ˙ p x h 2 x   ρ i x   ε ˙ p x ρ i x S v x f H A G B x v H A G B x ,         x   = α ,   β
The rate equation has the following terms:
  • a term with a hardening coefficient (h1) to describe the immobilization of mobile dislocations at forest dislocations, cross-slip around obstacles [47], and the formation of diploes of mobile dislocations with an anti-parallel Burgers vector [49];
  • a term with a recovery coefficient (h2) to describe the consumption of immobile dislocations that occurs due to their rearrangement and annihilation [40] by DRV;
  • a softening term to describe the consumption of immobile dislocations that are swept by the movement of HAGBs during cDRX.
The definitions of the terms h1 and h2 are given in Appendix H.
Model G
The rate equations for immobile dislocation are modified in Ref. [19] from Ref. [18] to include the consumption of immobile dislocations due to the movement of HAGBs during cDRX; Equation (18).
d ρ i x d t   =   ρ m x v g x Φ s g x 8 ρ i x v c i x λ i x δ D R V x ρ m x ρ i x ρ i x S v x f H A G B x v H A G B x ,         x   = α ,   β
The production of immobile dislocations occurs via the immobilization of mobile dislocations (first term), while SRV and DRV are responsible for the annihilation of immobile dislocations (second and third terms, respectively).

3.5.3. Boundary Density and Misorientation Distribution

The model of cDRX [46] assumes that an amount of dislocations ( Δ ρ c D R X ) either are rearranged into LAGBs or migrate into existing LAGBs, increasing their misorientation. The boundary density evolution is given by Equation (19) [46]. α C D R X is the fraction of dislocations that form new LAGBs, while (1 − α C D R X ) is the fraction that increases the misorientation of existing LAGBs. θ n e w is the mean misorientation angle of a newly formed boundary. The second term corresponds to the boundary surface density that is swept by the movement of HAGBs during cDRX.
d S v x d t = α C D R X x b x n x θ n e w x Δ ρ c D R X x     S v x 2 f H A G B x v H A G B x , x = α ,   β
The evolution of the misorientation distribution is given by the increase in the mean boundary misorientation angle of the Rayleigh distribution, and is given in Equation (20) [46].
d θ R x ¯ d t   =   ( 1 α C D R X x ) b n x S v x Δ ρ c D R X x , x = α ,   β
The increase in the average boundary misorientation of the Rayleigh distribution leads to a change in f R (Equation (A12) in Appendix B). The boundary misorientation angle distribution is calculated using Equation (A13) for the new values of f R and A 1 .
Model KM
The amount of dislocations that contribute to the boundary density and misorientation distribution evolution is the amount of recovered immobile dislocations, Equation (21).
Δ ρ c D R X   =   h 2   ρ i   ε ˙ p Δ t
Model G
The amount of dislocations that contribute to the boundary density and misorientation distribution evolution is a fraction ( f C D R X , see Appendix A) of the mobile and immobile dislocations that are annihilated by SRV and DRV, and is given in Equation (22).
Δ ρ c D R X = f C D R X ( 8 ρ m v c m λ m S R V + δ D R V ρ m ( ρ m + ρ i ) v g D R V + 8 ρ i v c i λ i S R V + δ D R V ρ m ρ i v g D R V ) Δ t

3.5.4. Wall Dislocation Density, Grain and Subgrain Size

The evolution of the boundary dislocations density is related to the evolution of the subgrain size; Equation (9). The grain size ( G s ) and the subgrain size ( Φ s g ) of the β-phase are calculated according to Equations (23) and (24), respectively.
G s   =   2 f H A G B S v
Φ s g   =   2 S v
Finally, the wall dislocation density ( ρ w ) is updated according to Equation (A14) at each iteration for both phases.

4. Results and Discussion

The measured and simulated flow stresses are compared for both Ti-5553 and Ti-17 alloys. Dynamic α-globularization and cDRX are interpreted from the experimental results. The differences in the predictions of the two developed models are discussed and compared.

4.1. Measured Microstructure

DRV and SRV occur by the annihilation and reorganization of dislocations during plastic deformation. The formed boundaries:
  • in the platelets of the α-phase promote α/β interface migration and the progressive globularization of the α-phase;
  • within the β-phase as the consequence of the progressive increase in boundary misorientation angle.
Figure 1, Figure 2 and Figure 3 show typical microstructures before and after deformation for the two near-β alloys. The microstructures of the α- and β-phases evolve as follows.

4.1.1. α-Phase

  • Ti-5553 exhibits an initial globular α-phase microstructure (Figure 2a), and slightly elongated α-particles due to plastic deformation after hot deformation at 820 °C and 0.001 s−1; Figure 2b. Ti-17 has an initial lamellar α-phase; Figure 2c. The α-phase becomes partially dynamically globularized after hot deformation at 810 °C and 0.001 s−1; Figure 2d.

4.1.2. β-Phase

  • Ti-5553: Figure 1a shows elongated β-grains perpendicular to the forging direction before deformation at 820 °C. LAGBs and HAGBs were formed in the β-phase after deformation at 820 °C and 0.001 s−1; Figure 3a. The dissolution of the α-phase beyond 920 °C and the higher HAGBs mobility promotes the formation of a fully recrystallized microstructure after annealing, Figure 1b;
  • Ti-17: Figure 1c,d show partially statically recrystallized β-grains before deformation at 820 °C and 970 °C, respectively. Figure 3c shows a well-established β-subgrain structure with newly formed HAGBs between the partially globularized α-particles after deformation at 810 °C and 0.001 s−1.
  • Figure 3b,d show larger recovered subgrains with very few newly formed HAGBs for the Ti-5553 and Ti-17 deformed in the β-domain.

4.2. Identification of the Restoration Phenomena

The intricate α-lamellar structure in the Ti-17 leads to the rotation, kinking or bending of the lamellas during plastic deformation. Dislocations are formed and rearranged into LAGBs as indicated in the KAM map in Figure 4a by the red arrows. If boundaries are formed within the α-phase, migration of the α/β interface occurs, leading to the formation of globular α-grains, as shown in Figure 4b (yellow arrows).
The high SFE of the β-phase facilitates the reorganization of dislocations into subgrain boundaries via DRV. The microstructure is comparable for the Ti-5553 deformed at 920 °C (Figure 3b) and the Ti-17 deformed at 970 °C (Figure 3d). The low fraction of newly formed HAGBs indicates the early stage of cDRX for those deformation conditions. On the other hand, the behavior of the β-phase deformed in the α+β domain notably differs between the Ti-5553 (Figure 3a) and the Ti-17 (Figure 3d). Figure 4 shows that the high values of KAM are only associated with LAGBs (white arrows) and not accumulated at α/β interfaces. Thus, the rotation, bending or kinking of the α-phase not only leads to α-globularization but also promotes the formation of boundaries within the β-phase [20,57]. The globular and lower fraction of the α-phase leads to an earlier stage of cDRX of the β-phase for the Ti-5553 (Figure 3a) in comparison to the larger α-fraction with lamellar morphology in Ti-17 (Figure 3c). Figure 4d shows a grain reference orientation deviation axis map. The different axis of rotation of each subgrain shows that progressive lattice rotation leads to an increment in boundary misorientation, i.e., cDRX.

4.3. Flow Curves

Figure 5 shows the flow curves for Ti-5553 (Figure 5a,b) and Ti-17 (Figure 5c,d). The initial load partitioning regime of the initial α-globular Ti-5553 microstructure corresponds to the iso-power regime, while the lamellar Ti-17 is associated with an initial iso-strain load partitioning. The steady-state load partition is the iso-stress regime. The influence of the hard α-phase in the overall flow stress is stronger for the iso-strain regime compared to the iso-power one. Thus, Ti-17 exhibits a more pronounced hardening as well as softening in the α+β-domain compared to the Ti-5553. For the β-domain, both alloys exhibited a fast work-hardening reaching a plateau-like steady-state.
The overall differences in the flow stress evolution between the “Model KM” and the “Model G” are minor. Figure 6a,b shows the evolution of the stress components predicted for the Ti-17 deformed at 970 °C, and 0.001 s−1 and 10 s−1, respectively, up to large strains. A slight softening is observed at 0.001 s−1 and corresponds to the decrement of the factor of strengthening ( F w , see Appendix C) due to the progressive increase in the average boundary misorientation angle. The differences between “Model KM” and “Model G” for the thermal stress are related to the different initial dislocation densities assumed as initial values for each model. However, if both models are initialized with the same dislocation densities, there is no difference in thermal stress between “Model KM” and “Model G”. “Model G” exhibits a faster increase in the athermal stress at the beginning of deformation compared to “Model KM”.
The derivative of the flow stress with respect to the strain is shown in Figure 7 for the Ti-17 deformed at 970 °C. “Model G” predicts a notably faster work-hardening at the beginning of deformation compared to “Model KM” for a given strain rate. The measured values are closer to “Model KM” than to “Model G”. The lower the strain rate, the smaller the strain required to achieve the same value of the derivative at the beginning of the deformation.

4.4. Microstructure Evolution Predictions

“Model G” can predict dislocation reactions. Moreover, the existing fitting parameters δ S G and δ D R V can be associated with a physical concept. The δ S G is the fraction of the mobile dislocations that contribute to the production of new mobile dislocations at boundaries, and δ D R V is a critical distance for the annihilation of dislocations via DRV, thus is not smaller than the Burgers vector, but is in the µm or nm range. “Model G” can predict the evolution of internal variables, such as climb velocity, that are not included in “Model KM” due to simplified dislocation evolution rates. Moreover, every term of the dislocation rate represents a different phenomenon, of which the contribution can be identified and correlated to the deformation conditions (temperature, strain and strain rate). The following chapters discuss the differences between both models in terms of the prediction of the evolution of the α- and β-phases.

4.4.1. α-Phase

Figure 8 shows the influence of the different models on the predictions of the α-evolution for Ti-17 deformed at 810 °C. A fast increase in the strain in the α-phase occurs at the beginning of deformation (Figure 8a) during the iso-strain regime. With the increase in the fraction of globularization (Figure 8e), the strain rate in the α-phase decreases drastically and the iso-stress regime becomes active. The globularization fraction of 100% does not necessarily correspond to a steady-state aspect ratio of 1 (Figure 8c) or that the formation of boundaries within the α-phase is finished (Figure 8d). The fraction of Rayleigh boundary misorientation distribution within the α-phase (Figure 8f) requires larger strains to reach negligible values compared to the globularization fraction.
The differences between “Model KM” and “Model G” for the evolution of the microstructure features for the α-phase are relatively large—a faster α-globularization is predicted by “Model KM” because a faster increase in average boundary misorientation is calculated. The differences between both models are associated with the difference in the amount of dislocations that contribute to the cDRX predicted by both models. They are adjusted by the constant f C D R X for “Model G”, while the behavior and kinetics are given by α C D R X and θ 0 , A g l o b and B g l o b . The difference in the glide velocity predictions between “Model KM” and “Model G” leads to the difference in the kinetics of the α/β interface movement related to the constant A g l o b . The slight differences in the flows stress evolutions between “Model KM” and the “Model G” in the α+β-domain in Figure 5 are related to the different α-globularization kinetics predicted by each model.

4.4.2. β-Phase

Figure 9 shows the grain and subgrain size evolutions for Ti-17 deformed at 970 °C and different strain rates. A fast decrease in subgrain size occurs at the beginning of deformation due to the formation of LAGBs. After a specific strain, the average grain size starts to decrease due to the formation of new HAGBs. The grain size finally decreases with slower kinetics at larger strains, reaching a steady-state subgrain size. The difference in kinetics between the “Model KM” and the “Model G” is associated with the difference in the amounts of dislocations that contribute to the cDRX predicted by both models. They are adjusted by the constant f C D R X for “Model G”, while the behavior and kinetics are given by α C D R X and θ 0 . However, negligible differences can be observed in the behavior since the same cDRX rate equations are used for both models.
Figure 10 shows the evolutions of the immobile and wall dislocation densities, as well as the fraction of HAGBs. Wall dislocations (Figure 10b,e) are rapidly produced due to the reorganization of immobile dislocations in the case of “Model KM” and due to the reorganization of mobile and immobile dislocations via SRV and DRV in case of “Model G”. Higher strain rates lead to faster kinetics initially for the formation of wall dislocation due to the higher rates of DRV. A maximum wall dislocation density follows the rapid increase. It corresponds to the strain whereat the consumption counterbalances the production of new LAGBs. After this maximum, a progressive decrease occurs. Higher strain rates form smaller subgrains sizes (Figure 9). Consequently, a higher rate of dislocation is required to evolve the formed LAGBs. Thus, the higher the strain rate, the slower the pace at which the wall dislocation density reaches the steady state. The fraction of HAGBs shows a fast initial decrease related to the formation of new LAGBs, followed by a valley that corresponds to the strain whereat the boundary misorientation distribution yields the maximum fraction of LAGBs.
The differences are more marked in the evolution of the immobile dislocation density between “Model KM” and “Model G”. “Model KM” predicts a rapid increase in immobile dislocation followed by saturation. Thus, the fast production of immobile dislocations at the early stage of deformation is balanced by their annihilation via DRV. The “Model G” predicts a more complex behavior for the immobile dislocation density evolution. An initial increase occurs due to the immobilization of mobile dislocations at LAGBs and HAGBs. The immobilization competes with the consumption of immobile dislocations via SRV due to the climb, DRV and movement of HAGBs. The immobile dislocation density does not reach an apparent plateau because the subgrain size varies during deformation (Figure 9), and so does also the immobilization rate. It is also notable that the predicted immobile dislocation density using “Model G” is smaller compared to “Model KM”. The main difference is that the Kocks–Mecking formalism [21] adopted in “Model KM” considers the production of immobile dislocations to be independent of the microstructure features. In contrast, in “Model G” the production rate of immobile dislocation density is dependent on the boundary density.
The mobile dislocation density is considered constant and independent of temperature and strain rate for the hot deformation modeled using “Model KM”. Thus, the increase in strain rate leads to the same proportional increment in glide velocity in “Model KM” following Orowan’s relationship (Equation (13)). “Model G” predicts an initial near-elastic regime, where the plastic strain rate is negligible (Figure 11a,b). With an increase in stress, the mobile dislocation starts to glide once the thermal stress becomes positive. A peak in the glide velocity and a sharp increase in the plastic strain rate are predicted once the yield stress is reached. The fast multiplication of mobile dislocations follows this peak in the glide velocity. Finally, the reduction in glide velocity and stabilization of the plastic strain rate occurs. A steady-state value of mobile dislocation density, plastic strain rate and mobile dislocation density is reached at very low strains. The sharp increase in mobile dislocation density at the beginning of deformation is responsible for the fast work-hardening predicted by “Model G” (Figure 7) since the predicted mobile dislocation density is higher than the immobile one in “Model G” (Figure 10a,d). Thus, the evolution of the flow stress for “Model G” is related to the mobile dislocation density evolution, with the minor influence of the wall and immobile dislocation density. The differences in mobile dislocation density between 930 °C and 970 °C, and between the different strain rates, are substantially smaller (from ~2 × 1013 for 0.001 s−1 to ~1 × 1014 for 10 s−1; Figure 11c) in comparison to the differences in the immobile and wall dislocation densities for the same strain rate difference (Figure 10).

4.5. “Model G”: Contributions to Deformation

One notable difference between “Model KM” and “Model G” is the capability of predicting the contribution of the different deformation phenomena and their impact on the mobile and immobile dislocation density evolutions. Figure 12 shows the evolution of the different dislocation reactions for the mobile (Figure 12a,b) and immobile dislocation density (Figure 12c,d) for Ti-17 deformed at 970 °C. The production rate of mobile dislocations by the Frank–Read source ( ρ m 3 / 2 v g ) follows the same kinetics as the annihilation, mainly DRV ( δ D R V ρ m ( ρ m + ρ i ) v g ). The small contributions of the immobilization rate ( ρ m v g / Φ s g ), SRV due to climb ( 8 ρ m v c m / λ m ) and consumption of mobile dislocations due to the movement of HAGBs ( ρ m S v f H A G B v H A G B ) are also predicted.
Annihilation via DRV ( δ D R V ρ m ρ i ) balances the production of mobile dislocations via the immobilization of mobile ones ( ρ m v g / Φ s g ) (Figure 12c,d). A minor contribution of the SRV due to climb ( 8 ρ i v c i / λ i ) and the consumption of mobile dislocations due to the movement of HAGBs ( ρ i S v f H A G B v H A G B ) are also predicted. A higher strain rate (Figure 12b,d) yields a lower annihilation rate of dislocations via SRV compared to a lower strain rate (Figure 12a,c) for both mobile and immobile dislocation densities. This occurs due to the decrease in time for SRV with the increase in strain rate.

4.6. “Model G”: Strain Rates

“Model G” predicts elastic and plastic strain rates. The elastic strain rate causes the variation in the stress in each phase following Hook’s law (Equation (16)). Thus, the evolution of the elastic strain rate is an indicator of the hardening/softening behavior; negative elastic strain rates correspond to flow softening, while positive ones to work-hardening. Figure 13a shows the influence of the thickness of the initial α-platelets on the hardening/softening behavior of the Ti-17. An initial hardening that occurs due to the multiplication of dislocations within the α- and β-phases is followed by a softening up to the end of the α-globularization. A more pronounced softening is observed for an initial thickness of 0.3 µm compared to 5 µm. the blue arrow in Figure 13a indicates a change in the kinetics of α-globularization at strain ~1 due to a change in the evolution of the α-globularization once the aspect ratio reaches the value of 1. Figure 13b shows the influence of the strain rate on the hardening kinetics of in Ti-17 deformed at the β-domain. An initial pronounced hardening due to the fast production of mobile dislocations (Figure 11a,b) is followed by a progressive decrease in the hardening. Finally, a slight softening occurs due to the reduction in the factor of strengthening of the wall dislocations when the average boundary misorientation increases. The higher the strain rate, the larger the strain required to achieve this flow softening (Figure 13b).

4.7. Predictions for Different Initial Microstructure

The developed models are robust in predicting the evolution of any starting microstructure during deformation. The influence of the previous β-grain size on the immobile dislocation density (Figure 14a,b) and the fraction of HAGB (Figure 14c,d) evolutions for Ti-17 deformed at 830 °C and 0.01 s−1 is shown in Figure 14 for both “Model KM” (Figure 14a,c) and “Model G” (Figure 14b,d). The smaller the initial β-grain size, the higher the predicted immobile dislocation density for a given temperature, strain and strain rate (Figure 14b). This occurs due to the higher immobilization rate (Equation (18)) for smaller grain sizes. On the other hand, “Model KM” does not have any term that describes the effect of the boundary density on the production or annihilation of the immobile dislocations. Thus, the previous β-grain size does not influence the immobile dislocation density evolution.
Smaller initial β-grain sizes evolve slower than larger ones (Figure 14c,d). The larger the initial grain size, the faster the increment in boundary average misorientation (Equation (20)). Moreover, smaller grain sizes promote the consumption of the formed low angle grain boundaries due to the movement of high angle grain boundaries (Equation (19)). Therefore, larger initial β-grain sizes lead to a fast evolution of the fraction of high angle grain boundary for a given temperature and strain rate.

4.8. Model and Measurements

Figure 15 shows the comparison between measured and simulated subgrain size and the fraction of HAGBs for Ti-5553 and Ti-17. The differences between “Model KM” and “Model G” are higher for the Ti-5553, mainly due to the adopted values of α C D R X , θ 0 , f C D R X . The simulated and measured values exhibit non-linear behaviors, i.e., an intricate dependency on the strain, strain rate and the temperature is attained. Figure 9 and Figure 10 show that the evolution of the microstructure is a complex function of the input parameters. However, as can also be seen in Figure 14, the developed models are robust in predicting the evolution of any initial microstructure up to any strain, temperature and strain rate.
The differences that are observed can be related to both model assumption and measurements issues.
  • Model assumptions: the mobile dislocation density is considered constant in “Model KM”, and no term to describe its influence on the cDRX is given. This issue is overcome in “Model G”, but with the assumption that the mobile dislocation is adjusted to yield a constant value once a maximum value is achieved. Additionally, its calculation is performed, assuming a constant thermal stress for a given temperature and strain rate. The limitations on the prediction of mobile dislocation density affect directly both the stress and microstructure evolutions. The effect of texture formation on the thermal stress, or the Taylor factor, is not considered, affecting the overall prediction of the stress evolution, and indirectly the microstructure evolution. Another limitation is the phenomenological description of the HAGB velocity and its power–law dependency uniquely via the strain rate (see Appendix G). This strongly affects the evolution of the boundary density and the average boundary misorientation. In the case of “Model G”, the dependence of the boundary evolution on the applied strain rate is greatly influenced by the phenomenological description of the δ D R V (see Appendix A), as well as by the constant f C D R X (see Appendix A).
  • Measurements issues: the limited measurement area size impacts the statistical significance of the measured values. The large previous β-grain size limits the accuracy of the determination of the fraction of HAGBs, since more/less initial β-HAGBs can be considered depending on the measurement site. The formation of the substructure occurs heterogeneously within the β-grain since recovery occurs preferentially in the vicinity of the previous β-HAGBs [58]. Moreover, DRV and consequently cDRX occur differently in each grain, since the Schmid factor of each grain in a polycrystal is different. Thus, the limited measurement area size leads to local information on subgrain size and the fraction of HAGBs that is strongly influenced by the area wherein the EBSD measurement was performed.

5. Summary and Conclusions

Two physical models were developed to predict the flow stress and microstructure evolutions of two near-β titanium alloys. A simple dislocation density reaction approach (“Model KM”) is compared to a comprehensive one (“Model G”). cDRX and DRV are the predominant restoration mechanisms. On the other hand, SRV and other phenomena play a negligible role in the restoration of dislocations, as described by “Model G”. The evolution of the mean grain size, the mean subgrain size, the boundary misorientation distribution, and the dislocation densities of all phases can be predicted during deformation. The models can predict the evolution of any initial microstructure. The change in load partitioning mechanism in the α+β-domain is considered to be the reason for the flow softening, and was correlated with the dynamic globularization of the α-phase. The measurements and simulations allow the following conclusions:
  • The mobile dislocation density can be estimated by “Model G” for any deformation condition, while it needs to be assumed for the “Model KM”;
  • The kinetics of immobile dislocation density evolution depend on the kinetics subgrain size for “Model G”. “Model KM” predicts a fast increase in immobile dislocation density followed by saturation;
  • The applied and plastic strain rate are considered the same for “Model KM”, and the calculations start at the yield point. The separation between elastic and plastic strains allows the modeling of the elastic part using “Model G”, as well as the interpretation of the kinetics of hardening and softening in the material;
  • Both models predict the dynamic globularization of the α-phase for any initial thickness and aspect ratio, and its evolution is not only described in terms of the fraction of globularization but also the evolution of the microstructural features (thickness, width, the fraction of HAGBs within the α-phase);
  • The load partition model describes the overall flow softening in the α+β-domain. The predicted change in the mechanism describes the phenomena related to a decrease in dislocation density, strain rate, and consequently, the stress in the α-phase;
  • cDRX is coupled with the dislocation reactions in the same manner for the “Model KM” and “Model G”—recovered dislocations are the source of the formation of new boundaries and increases in their misorientation. However, in “Model G”, only a fraction of the recovered dislocation participates in cDRX;
  • The influence of the dislocation reaction approach on cDRX-related variables (subgrain and grain sizes, misorientation distribution, the fraction of HAGBs) is small since the amount of dislocations that contribute to cDRX are comparable between both models;
  • “Model G” can be used to predict any dislocation-based phenomena, since the dislocation reactions are defined, and the strain is separated between elastic and plastic.

Author Contributions

Conceptualization, R.H.B., F.M.B.F., M.L., A.K. and M.C.P.; data curation, M.C.P.; investigation, R.H.B.; methodology, R.H.B., F.M.B.F. and M.C.P.; project administration, M.L., A.K. and M.C.P.; resources, M.L., A.K. and M.C.P.; software, R.H.B.; supervision, M.C.P.; validation, R.H.B.; visualization, R.H.B.; writing—original draft, R.H.B.; writing—review and editing, F.M.B.F., M.L., A.K. and M.C.P. All authors have read and agreed to the published version of the manuscript.

Funding

The authors carried out this work in the framework of the D-1303000107/CD-Laboratory for Design of High-Performance Alloys by Thermomechanical Processing with the support of the Christian Doppler Forschungsgesellschaft. Fonds zur Förderung der wissenschaftlichen Forschung, FWF, is fostered for the project F37729.

Acknowledgments

Shoichi Nambu is kindly acknowledged for his support and allowance to use the cross-section polisher (JEOL SM-09010) and the scanning electron microscope (JEOL JSM7001F). Authors accept Open Access Funding by the Graz University of Technology.

Conflicts of Interest

The authors declare no conflict of interest.

Data Availability

The raw/processed data that support the findings of this study are available from the corresponding author, RHB, upon request.

Appendix A

The parameters of the initial microstructure were obtained by the EBSD and BSE measurements of the microstructure after annealing and quenching, Table A1.
Table A1. Measured microstructural parameters corresponding to the initial microstructure used in the developed model.
Table A1. Measured microstructural parameters corresponding to the initial microstructure used in the developed model.
α-Phaseβ-Phase
Ti-5553Ti-17Ti-5553Ti-17
Initial microstructure G s 0 in α+β domain (µm)3-300500
G s 0 in β domain (µm)--200500
Φ s g 0 in α+β domain (µm)3- 250
Φ s g 0 in β domain (µm)--200250
t α 0 (µm)-3--
A R (-)-15--
θ R ¯ in α+β domain (°)0.10.440.5
θ R ¯ in β domain (°)--0.10.01
Other parameters are directly obtained from values established in the literature, Table A2.
Table A2. Parameters obtained from the literature and used in the developed model.
Table A2. Parameters obtained from the literature and used in the developed model.
α-phaseβ-phase
Ti-5553Ti-17Ti-5553Ti-17
Parameters obtained from literatureb (nm)0.2950.2950.2860.286
M (-)553.053.05
Ω (m3)1.7668 × 10−29
Other constants are adjusted according to values established in the literature, Table A3. They are fitted by simultaneously adjusting the flow stress and microstructural results with respect to the measured ones.
Table A3. Parameters adjusted according to a range established in the literature and used in the developed model.
Table A3. Parameters adjusted according to a range established in the literature and used in the developed model.
α-Phaseβ-Phase
Ti-5553Ti-17Ti-5553Ti-17
Parameters fitted from a range established in the literature α C D R X (-)Model KM0.40.150.20.19
Model G0.40.150.20.19
Q H A G B (Jmol−1)Model KM1.97 × 1051.97 × 1051.97 × 1051.97 × 105
Model G1.97 × 1051.97 × 1051.97 × 1051.97 × 105
n (-)Model KM30153030
Model G30153030
θ 0 (°)Model KM10.20.80.7
Model G10.20.80.7
θ C (°)Model KM12121212
Model G12121212
ε ˙ r e f (s−1)Model KM1 × 1061 × 1061 × 1061 × 106
Model G1 × 1061 × 1061 × 1061 × 106
Other parameters were fitted by simultaneously adjusting the simulated microstructural features with respect to the measured ones, Table A4.
Table A4. Fitting parameters used in the developed model.
Table A4. Fitting parameters used in the developed model.
α-Phaseβ-Phase
Ti-5553Ti-17Ti-5553Ti-17
Model KM Q 1 (Jmol−1)1.2 × 1041.2 × 1041.2 × 1041.2 × 104
Q 2 (Jmol−1)1.2 × 1041.2 × 1041.2 × 1041.2 × 104
h 1 0 (m−2)Equation (A1)
h 1 K M (m−2)2 × 10162.2 × 10164.80 × 10164 × 1014
A h 1 S v (m)01.25 × 10−23.5 × 10−31.25 × 10−2
η h 1 K M (-)-6.6 × 10−116.6 × 10−1
h 2 0 (-)Equation (A3)
h 2 K M (-)3 × 10−13 × 10−31.31.1 × 101
A h 2 S v (m)01.25 × 10−23.8 × 10−31.25 × 10−2
η h 2 K M (-)-6.6 × 10−116.6 × 10−1
E M O B 0 (Wm−1)3.5 × 10206 × 1012 × 104
F w 0 (-)3 × 10−21 × 10−23 × 10−21 × 10−2
n v (-)4 × 10−1-4 × 10−16.5 × 10−1
m 1 (-)3.3 × 10−18 × 10−23.3 × 10−18 × 10−2
m 2 (-)3.3 × 1016 × 10−13.1 × 10−18 × 10−2
A g l o b (-)-4 × 102--
B g l o b (-)-3 × 10−3--
A i s o s t r e s s (-)6 × 10−2 (only Ti-5553)
n s (-)3 (only Ti-5553)
Model G W g (Jmol−1)2.1 × 1052.1 × 1052.1 × 1052.1 × 105
δ D R V (m−2)Equation (A4)
δ D R V 0 (m−2)6.75 × 10−73.20 × 10−75.62 × 10−73.50 × 10−7
η D R V (-)−1 × 25 × 10−1−6.5 × 10−2−1 × 25 × 10−1−7.75 × 10−2
η (-)2 × 10−32 × 10−32 × 10−32 × 10−3
E M O B 0   (Wm−1)--2.4 × 1042.4 × 104
η v (-)--6.5 × 10−16.5 × 10−1
f C D R X (-)--Equation (A5)Equation (A6)
f C D R X 0 (-)2 × 10−22 × 10−21.25 × 10−17.5 × 10−3
  f g l o u l a r i n t e r f a c e (m2)---1.2 × 10−10
  f p l a t e s i n t e r f a c e (m)---2.1 × 10−5
n f C D R X (-)---6.7 × 10−1
A g l o b (-)-4 × 102--
B g l o b (-)-3 × 10−2--
Some parameters are empirically adjusted according to expressions. The h 1 0 is given by Equation (A1) to consider the role of the α/β interface density ( S i n t e r f a c e , Equation (A2)) in the production of immobile dislocations.
h 1 0 = ( 1 + ( A h 1 S v S i n t e r f a c e ) η h 1 K M ) h 1 K M
S i n t e r f a c e = { 2 F v t α 0 , T i 17 3 F v r α 0 , T i 5553
r α 0 is the mean radius of the α-particle in Ti-5553 and is equal to G s 0 . The h 2 0 is also given as a function of the α/β interface density; Equation (A3).
h 2 0 = ( 1 + ( A h 2 S v S i n t e r f a c e ) η h 2 K M ) h 2 K M
The δ D R V is given as a power–law function of the plastic strain rate; Equation (A4).
δ D R V = δ D R V 0 ( ε ˙ p ε ˙ r e f ) η D R V
The f C D R X is given as a function of the S i n t e r f a c e and is given by Equations (A5) and (A6) for the Ti-5553 and Ti-17, respectively.
f C D R X   =   f C D R X 0 +     f g l o u l a r i n t e r f a c e S i n t e r f a c e λ i n t e r s p a c e ,         T i 5553
f C D R X   =   f C D R X 0 +   (   f p l a t e s i n t e r f a c e S i n t e r f a c e   ) n f C D R X ,         T i 17
λ i n t e r s p a c e is the particle interspace and is calculated according to Equation (A7) [58] assuming a random distribution of spherical α-phase for the Ti-5553.
λ i n t e r s p a c e   = ln ( 3 ) 2 π n d R α 0 + 8 R α 0 2 3 8 3 R α 0 ,         n d =   3 f α 4 π ( R α 0 ) 3

Appendix B

The initial subgrain ( Φ s g 0 ) and grain ( G s 0 ) sizes, average boundary misorientation angle and dislocation densities define the initial microstructure. Equation (A8) [46] describes the maximum boundary surface density that is given by the total surface of boundaries (HAGB + LAGB) over the total volume.
S v 0 = 2 Φ s g 0
Equation (A9) describes the HAGB surface density corresponding to the initial grain size.
S v H A G B 0 = 2 G s 0
Equation (A10) describes the initial boundary surface density prone to evolve during the deformation.
S v C D R X 0 = S v 0 S v H A G B 0
For an initial fully recrystallized microstructure, S v 0 S v H A G B 0 . Finally, a multiplicative factor A 1 0 that represents the fraction of boundaries prone to evolve in misorientation during the deformation is calculated using Equation (A11).
A 1 0 = S v C D R X 0 S v 0
The misorientation ranges given by the Rayleigh distribution, and the Mackenzie distributions range between 0 and 60.7° ( 1.06 rad) [50]. Equation (A12) represents the initial fraction of Rayleigh distribution, f R 0 .
f R 0 = 0 60.7 Θ R 0 ( θ ) d θ
The fraction ( 1 A 1 0 f R 0 ) describes the boundary misorientation angle distribution with a Mackenzie distribution. Equation (A13) describes the initial boundary misorientation distribution.
Θ 0 ( θ ) = A 1 0 f R 0 Θ R 0 ( θ ) + ( 1 A 1 0 f R 0 ) Θ M 0 ( θ )
If the initial microstructure is fully recrystallized, A 1 0 = f R 0 = 0, and the boundary misorientation follows a Mackenzie distribution, while the boundary misorientation distribution of an initial microstructure with a high fraction of LAGB is mainly represented by a Rayleigh distribution ( A 1 0 f R 0 1).
The initial wall dislocation density ( ρ w 0 ) is calculated according to Equation (9).
ρ w 0 = n S v 0 θ ¯ L A G B 0 f L A G B 0 b
θ ¯ L A G B 0 is the mean misorientation of the boundary distribution, and f L A G B 0 is the initial fraction of LAGBs and is given according to Equation (A15).
f L A G B 0 = 0 12 Θ 0 ( θ ) d θ
The value 12° ( 0.21 rad) is the transition angle between an LAGB and an HAGB.

Appendix C

Non-sharp LAGBs contributes to a higher internal athermal stress in comparison to sharp LAGBs (considered negligible [18]). An empirical factor ( F w ) gives the strain field size caused by a non-sharp LAGB (Equation (A16)).
F w x = F w x 0 0 3 Θ x ( θ ) d θ ,         x =   α ,   β
The transition angle of 3° ( 0.05 rad) describes the transition between non-sharp and sharp LAGBs. When the substructure is forming, θ ¯ < 3°, and the contribution of the non-sharp LAGBs is higher. An increase in boundary misorientations decreases the F w value asymptotically to zero.

Appendix D

A simple approach is considered to calculate the yield stress assuming an Arrhenius-type behavior (Equation (A17)) [59].
σ Y S = 1 α Y S ln { ( ε ˙ exp ( Q Y S R T ) A Y S ) 1 n Y S + [ ( ε ˙ exp ( Q Y S R T ) A Y S ) 1 2 n Y S + 1 ] 1 2 }
where ε ˙ is the strain rate, R is the universal constant of gases (8.314 J.mol−1.K−1) and T the temperature. The values of the activation energy ( Q Y S ), α Y S , A Y S , and n Y S that define the yield stress of the α- and β-phases are listed in Table A5. A detailed procedure for the parameter calculation was published elsewhere [59]. For both alloys, the parameters of the yield stress of the β-phase are obtained by fitting the measured yield stress in the β-domain with Equation (A17). The yield stress of the β-phase is then extrapolated from the values obtained for the α+β-domain. The yield stress of the α-phase is estimated based on the yield stress values in the α+β-domain using a simple law of mixtures for the fractions of α- and β-phases. Finally, the parameters of the yield stress of the α-phase are obtained by fitting with Equation (A17) the estimated yield stress of the α-phase in the α+β-domain.
Table A5. Constants used for the calculation of the yield stress.
Table A5. Constants used for the calculation of the yield stress.
α-Phaseβ-Phase
Ti-5553Ti-17Ti-5553Ti-17
Q Y S (Jmol−1)1.58 × 1051.58 × 1052.38 × 1052.38 × 105
α Y S (MPa−1)2.85 × 10−32.80 × 10−31.16 × 10−21.85 × 10−2
A Y S (s−1)2.22 × 1061.78 × 1062.17 × 1091.80 × 106
n Y S (-)3.334.243.232.94

Appendix E

In the iso-strain regime, Equation (A18) is fulfilled:
ε ˙ = ε ˙ α ε = ε ˙ β ε ,
ε ˙ , ε ˙ α ε and ε ˙ β ε are the applied strain rate, the strain rate in the α-phase, and the strain rate in the β-phase, respectively.
Equations (A19) and (A20) describe the iso-power model [55].
ε ˙ = F v ε ˙ α W + ( 1 F v ) ε ˙ β W
σ α ε ˙ α W = σ β ε ˙ β W
ε ˙ α W and ε ˙ β W are strain rates in the iso-power condition of the α- and β-phases, respectively. F v is the fraction of the α-phase.
For the iso-stress regime, the flow stress of the β-phase is equal to the flow stress for the α-phase. The strain rates in the α- and β-phases are calculated by Equation (A21) [59] and Equation (A22), respectively.
ε ˙ α σ = exp ( n Y S α ln ( sinh ( α Y S α σ β ) ) Q Y S α R T + ln ( A Y S α ) )
ε ˙ β σ = F v ε ˙ α σ ε ˙ ( 1 F v )

Appendix F

Figure A1 shows the process of dynamic α-globularization schematically. The migration of the formed and existing interfaces promoting the thickening of the α-platelet accompanies the formation of new α/β interfaces. The movement of the α/β interfaces consuming the existing HAGBs within the α-phase occurs with a velocity, v g l o b , that is empirically correlated with the glide velocity in the β-phase by a constant and fitting parameter of the model A g l o b (Equation (A23)).
v g l o b = A g l o b v g β
Semiatin et al. [60] proved that the growth of the spherical particles during dynamic globularization is a faster Ostwald ripening process driven by pipe diffusion. During the Ostwald ripening of non-spherical particles [61], complex solutions are necessary. Here, the growth of the α-particles is simplified assuming an empirical velocity of growth, v g r o w t h , related to the climb velocity of the mobile dislocation density in the β-phase by a constant and fitting parameter of the model B g l o b (Equation (A24)).
v g r o w t h =   B g l o b v c l i m b β
Figure A1. Schematic representation of the dynamic globularization of the α-phase for thin (ad) and thicker (eg) α-platelets: (a,e) formation of HAGBs within the α-phase; (b,f) migration of α/β interfaces though formed boundaries within the α-phase; (c) thickening of the formed particles; (d,g) formation of dispersed and globular α-phase.
Figure A1. Schematic representation of the dynamic globularization of the α-phase for thin (ad) and thicker (eg) α-platelets: (a,e) formation of HAGBs within the α-phase; (b,f) migration of α/β interfaces though formed boundaries within the α-phase; (c) thickening of the formed particles; (d,g) formation of dispersed and globular α-phase.
Materials 13 05678 g0a1
The evolution of the morphological features in the α-phase depends on the aspect ratio of the formed particles:
  • Aspect ratio larger than one:
Equations (A25)–(A27)describe the thickness evolution rate, the width evolution rate, and the aspect ratio, respectively.
d t α d t = 2 f g l o b v g r o w t h
d w α d t = A R 2 t α 3 d t α d t
A R = 2 w α 2 t α ( S g l o b w α 2 + 2 w α )
Equations F3 to F5 are used until the aspect ratio is larger than one. Even if f g l o b = 1, the aspect ratio is not necessarily one. A small α-thickness can lead to the consumption of the formed boundaries within the α-phase quickly, and the α-globularization ends with an aspect ratio larger than one (Figure A1b). This is followed by thickening (Figure A1c) until reaching the steady state (Figure A1d).
  • Aspect ratio equals to one:
Equation (A28) describes the thickness update for the iteration i + 1 ( t α i + 1 ) as a function of the density of the α/β interface formed during dynamic globularization of the α-phase ( S g l o b ) and the thickness from the iteration i ( t α i ).
t α i + 1   = 4 t α i S g l o b t α i   +   4 t α i
A large initial thickness of the α-platelet leads to their division into particles of smaller thickness than the initial one (Figure A1f). This occurs when the aspect ratio reaches one while α-globularization fraction is still not finished.

Appendix G

Equation (A29) [46] describes phenomenologically the velocity of HAGBs, v H A G B .
v H A G B x = v 0 x ( ε ˙ p x ε ˙ r e f ) n v x exp ( Q H A G B x R T ) ,         x =   α ,   β
v 0 is a coefficient, n v is an exponent and Q H A G B is the activation energy for the movement of HAGBs.

Appendix H

Equations (A30) and (A31) [46] describe the hardening and recovery coefficients, respectively.
h 1 x = h 1 x 0 ( ε ˙ x ε ˙ r e f ) m 1 x exp ( Q 1 x R T ) ,         x =   α ,   β
h 2 x = h 2 0 x ( ε ˙ x ε ˙ r e f ) m 2 x exp ( Q 2 x R T ) ,         x =   α ,   β
where h 1 0 and h 2 0 are the constants for the production of dislocations and dynamic recovery coefficients, respectively; m 1 and   m 2 are exponents related to the strain hardening coefficient; Q 1 and Q 2 are activation energies for the production of dislocations and DRV, respectively.

References

  1. Fan, X.G.; Zhang, Y.; Gao, P.F.; Lei, Z.N.; Zhan, M. Deformation behavior and microstructure evolution during hot working of a coarse-grained Ti-5Al-5Mo-5V-3Cr-1Zr titanium alloy in beta phase field. Mater. Sci. Eng. A 2017, 694, 24–32. [Google Scholar] [CrossRef]
  2. Chuan, W.; Liang, H. Hot deformation and dynamic recrystallization of a near-beta titanium alloy in the β single phase region. Vacuum 2018, 156, 384–401. [Google Scholar] [CrossRef]
  3. Qu, F.S.; Zhou, Y.H.; Zhang, L.Y.; Wang, Z.H.; Zhou, J. Research on hot deformation behavior of Ti–5Al–5Mo–5V–1Cr–1Fe alloy. Mater. Des. 2015, 69, 153–162. [Google Scholar] [CrossRef]
  4. Tan, K.; Li, J.; Guan, Z.; Yang, J.; Shu, J. The identification of dynamic recrystallization and constitutive modeling during hot deformation of Ti55511 titanium alloy. Mater. Des. 2015, 84, 204–211. [Google Scholar] [CrossRef]
  5. Poletti, C.; Germain, L.; Warchomicka, F.; Dikovits, M.; Mitsche, S. Unified description of the softening behavior of beta-metastable and alpha+beta titanium alloys during hot deformation. Mater. Sci. Eng. A 2016, 651, 280–290. [Google Scholar] [CrossRef]
  6. Warchomicka, F.; Poletti, C.; Stockinger, M. Study of the hot deformation behavior in Ti–5Al–5Mo–5V–3Cr–1Zr. Mater. Sci. Eng. A 2011, 528, 8277–8285. [Google Scholar] [CrossRef]
  7. Dikovits, M.; Poletti, C.; Warchomicka, F. Deformation Mechanisms in the Near-β Titanium Alloy Ti-55531. Metall. Mater. Trans. A 2014, 45, 1586–1596. [Google Scholar] [CrossRef]
  8. Matsumoto, H.; Kitamura, M.; Li, Y.; Koizumi, Y.; Chiba, A. Hot forging characteristic of Ti–5Al–5V–5Mo–3Cr alloy with single metastable β microstructure. Mater. Sci. Eng. A 2014, 611, 337–344. [Google Scholar] [CrossRef]
  9. Hiratani, M.; Zbib, H.M.; Khaleel, M.A. Modeling of thermally activated dislocation glide and plastic flow through local obstacles. Int. J. Plast. 2003, 19, 1271–1296. [Google Scholar] [CrossRef]
  10. Sun, Z.C.; Yang, H.; Han, G.J.; Fan, X.G. A numerical model based on internal-state-variable method for the microstructure evolution during hot-working process of TA15 titanium alloy. Mater. Sci. Eng. A 2010, 527, 3464–3471. [Google Scholar] [CrossRef]
  11. Peng, W.; Zeng, W.; Wang, Q.; Yu, H. Comparative study on constitutive relationship of as-cast Ti60 titanium alloy during hot deformation based on Arrhenius-type and artificial neural network models. Mater. Des. 2013, 51, 95–104. [Google Scholar] [CrossRef]
  12. Matsumoto, H.; Naito, D.; Miyoshi, K.; Yamanaka, K.; Chiba, A.; Yamabe-Mitarai, Y. Forging property, processing map, and mesoscale microstructural evolution modeling of a Ti-17 alloy with a lamellar (α+β) starting microstructure. Sci. Technol. Adv. Mater. 2017, 18, 893–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kocks, U.F. A statistical theory of flow stress and work-hardening. Philos. Mag. 1966, 13, 541–566. [Google Scholar] [CrossRef]
  14. Bergström, Y. A dislocation model for the stress-strain behavior of polycrystalline α-Fe with special emphasis on the variation of the densities of mobile and immobile dislocations. Mater. Sci. Eng. 1970, 5, 193–200. [Google Scholar] [CrossRef]
  15. Bergstrom, Y.; Roberts, W. The dynamical strain ageing of α-iron: Effects of strain rate and nitrogen content in the jerky-flow region. Acta Metall. 1973, 21, 741–745. [Google Scholar] [CrossRef]
  16. Roberts, W.; Bergström, Y. The stress-strain behavior of single crystals and polycrystals of face-centered cubic metals-a new dislocation treatment. Acta Metall. 1973, 21, 457–469. [Google Scholar] [CrossRef]
  17. Mecking, H.; Estrin, Y. The effect of vacancy generation on plastic deformation. Scr. Metall. 1980, 14, 815–819. [Google Scholar] [CrossRef]
  18. Ghoniem, N.M.; Matthews, J.R.; Amodeo, R.J. Dislocation model for creep in engineering materials. Res. Mech. 1990, 29, 197–219. [Google Scholar]
  19. Buzolin, R.H.; Lasnik, M.; Krumphals, A.; Poletti, M.C. A dislocation-based model for the microstructure evolution and the flow stress of a Ti5553 alloy. Int. J. Plast. 2021, 136. [Google Scholar] [CrossRef]
  20. Buzolin, R.H.; Lasnik, M.; Krumphals, A.; Poletti, M.C. Hot deformation and dynamic α-globularization of a Ti-17 alloy: Consistent physical model. Mater. Des. 2021, 197, 109266. [Google Scholar] [CrossRef]
  21. Mecking, H.; Kocks, U.F. Kinetics of flow and strain-hardening. Acta Metall. 1981, 29, 1865–1875. [Google Scholar] [CrossRef]
  22. Gao, C.Y.; Zhang, L.C.; Yan, H.X. A new constitutive model for HCP metals. Mater. Sci. Eng. A 2011, 528, 4445–4452. [Google Scholar] [CrossRef]
  23. Nemat-Nasser, S.; Guo, W.G.; Cheng, J.Y. Mechanical properties and deformation mechanisms of a commercially pure titanium. Acta Mater. 1999, 47, 3705–3720. [Google Scholar] [CrossRef] [Green Version]
  24. Picu, R.C.; Majorell, A. Mechanical behavior of Ti–6Al–4V at high and moderate temperatures—Part II: Constitutive modeling. Mater. Sci. Eng. A 2002, 326, 306–316. [Google Scholar] [CrossRef]
  25. Ning, Y.Q.; Luo, X.; Liang, H.Q.; Guo, H.Z.; Zhang, J.L.; Tan, K. Competition between dynamic recovery and recrystallization during hot deformation for TC18 titanium alloy. Mater. Sci. Eng. A 2015, 635, 77–85. [Google Scholar] [CrossRef]
  26. Xiao, Y.-W.; Lin, Y.C.; Jiang, Y.-Q.; Zhang, X.-Y.; Pang, G.-D.; Wang, D.; Zhou, K.-C. A dislocation density-based model and processing maps of Ti-55511 alloy with bimodal microstructures during hot compression in α+β region. Mater. Sci. Eng. A 2020, 790, 139692. [Google Scholar] [CrossRef]
  27. Souza, P.M.; Beladi, H.; Singh, R.; Rolfe, B.; Hodgson, P.D. Constitutive analysis of hot deformation behavior of a Ti6Al4V alloy using physical based model. Mater. Sci. Eng. A 2015, 648, 265–273. [Google Scholar] [CrossRef]
  28. Fan, X.G.; Yang, H. Internal-state-variable based self-consistent constitutive modeling for hot working of two-phase titanium alloys coupling microstructure evolution. Int. J. Plast. 2011, 27, 1833–1852. [Google Scholar] [CrossRef]
  29. Babu, B.; Lindgren, L.-E.E. Dislocation density based model for plastic deformation and globularization of Ti-6Al-4V. Int. J. Plast. 2013, 50, 94–108. [Google Scholar] [CrossRef]
  30. Gao, P.F.; Guo, J.; Zhan, M.; Lei, Z.N.; Fu, M.W. Microstructure and damage based constitutive modelling of hot deformation of titanium alloys. J. Alloys Compd. 2020, 831, 154851. [Google Scholar] [CrossRef]
  31. Gao, P.; Yang, H.; Fan, X.; Zhu, S. Unified modeling of flow softening and globularization for hot working of two-phase titanium alloy with a lamellar colony microstructure. J. Alloys Compd. 2014, 600, 78–83. [Google Scholar] [CrossRef]
  32. OuYang, D.; Fu, M.W.; Lu, S.Q. Study on the dynamic recrystallization behavior of Ti-alloy Ti–10V–2Fe–3V in β processing via experiment and simulation. Mater. Sci. Eng. A 2014, 619, 26–34. [Google Scholar] [CrossRef]
  33. Chuan, W.; He, Y.; Wei, L.H. Modeling of discontinuous dynamic recrystallization of a near-α titanium alloy IMI834 during isothermal hot compression by combining a cellular automaton model with a crystal plasticity finite element method. Comput. Mater. Sci. 2013, 79, 944–959. [Google Scholar] [CrossRef]
  34. Ning, Y.Q.; Xie, B.C.; Liang, H.Q.; Li, H.; Yang, X.M.; Guo, H.Z. Dynamic softening behavior of TC18 titanium alloy during hot deformation. Mater. Des. 2015, 71, 68–77. [Google Scholar] [CrossRef]
  35. Liang, H.; Guo, H.; Ning, Y.; Peng, X.; Qin, C.; Shi, Z.; Nan, Y. Dynamic recrystallization behavior of Ti–5Al–5Mo–5V–1Cr–1Fe alloy. Mater. Des. 2014, 63, 798–804. [Google Scholar] [CrossRef]
  36. Bobbili, R.; Venkata Ramudu, B.; Madhu, V. A physically-based constitutive model for hot deformation of Ti-10-2-3 alloy. J. Alloys Compd. 2017, 696, 295–303. [Google Scholar] [CrossRef]
  37. Quan, G.; Luo, G.; Liang, J.; Wu, D.; Mao, A.; Liu, Q. Modelling for the dynamic recrystallization evolution of Ti–6Al–4V alloy in two-phase temperature range and a wide strain rate range. Comput. Mater. Sci. 2015, 97, 136–147. [Google Scholar] [CrossRef]
  38. Li, H.; Sun, X.; Yang, H. A three-dimensional cellular automata-crystal plasticity finite element model for predicting the multiscale interaction among heterogeneous deformation, DRX microstructural evolution and mechanical responses in titanium alloys. Int. J. Plast. 2016, 87, 154–180. [Google Scholar] [CrossRef]
  39. Sun, X.; Li, H.; Zhan, M. Full-stage prediction of discontinuous dynamic recrystallization of a titanium alloy through a sub-mesh internal state variables method. Model. Simul. Mater. Sci. Eng. 2019, 27, 015004. [Google Scholar] [CrossRef]
  40. Li, H.; Wu, C.; Yang, H. Crystal plasticity modeling of the dynamic recrystallization of two-phase titanium alloys during isothermal processing. Int. J. Plast. 2013, 51, 271–291. [Google Scholar] [CrossRef]
  41. Jones, N.G.; Dashwood, R.J.; Dye, D.; Jackson, M. Thermomechanical processing of Ti–5Al–5Mo–5V–3Cr. Mater. Sci. Eng. A 2008, 490, 369–377. [Google Scholar] [CrossRef]
  42. Jones, N.G.; Dashwood, R.J.; Dye, D.; Jackson, M. The Flow Behavior and Microstructural Evolution of Ti-5Al-5Mo-5V-3Cr during Subtransus Isothermal Forging. Metall. Mater. Trans. A 2009, 40, 1944–1954. [Google Scholar] [CrossRef]
  43. Khan, A.S.; Liang, R. Behaviors of three BCC metal over a wide range of strain rates and temperatures: Experiments and modeling. Int. J. Plast. 1999, 15, 1089–1109. [Google Scholar] [CrossRef]
  44. Goetz, R.L.; Semiatin, S.L. The adiabatic correction factor for deformation heating during the uniaxial compression test. J. Mater. Eng. Perform. 2001, 10, 710–717. [Google Scholar] [CrossRef]
  45. Song, B.; Dong, S.; Zhang, B.; Liao, H.; Coddet, C. Effects of processing parameters on microstructure and mechanical property of selective laser melted Ti6Al4V. Mater. Des. 2012, 35, 120–125. [Google Scholar] [CrossRef]
  46. Gourdet, S.; Montheillet, F. A model of continuous dynamic recrystallization. Acta Mater. 2003, 51, 2685–2699. [Google Scholar] [CrossRef]
  47. Shanthraj, P.; Zikry, M.A. Dislocation density evolution and interactions in crystalline materials. Acta Mater. 2011, 59, 7695–7702. [Google Scholar] [CrossRef]
  48. Hansen, B.L.; Beyerlein, I.J.; Bronkhorst, C.A.; Cerreta, E.K.; Dennis-Koller, D. A dislocation-based multi-rate single crystal plasticity model. Int. J. Plast. 2013, 44, 129–146. [Google Scholar] [CrossRef]
  49. Ma, A.; Roters, F.; Raabe, D. A dislocation density based constitutive model for crystal plasticity FEM including geometrically necessary dislocations. Acta Mater. 2006, 54, 2169–2179. [Google Scholar] [CrossRef]
  50. Mason, J.K.; Schuh, C.A. The generalized Mackenzie distribution: Disorientation angle distributions for arbitrary textures. Acta Mater. 2009, 57, 4186–4197. [Google Scholar] [CrossRef]
  51. Semiatin, S.; Bieler, T. The effect of alpha platelet thickness on plastic flow during hot working of TI–6Al–4V with a transformed microstructure. Acta Mater. 2001, 49, 3565–3573. [Google Scholar] [CrossRef]
  52. Mughrabi, H. The α-factor in the Taylor flow-stress law in monotonic, cyclic and quasi-stationary deformations: Dependence on slip mode, dislocation arrangement and density. Curr. Opin. Solid State Mater. Sci. 2016, 20, 411–420. [Google Scholar] [CrossRef]
  53. Taylor, G.I. Plastic strain in metals. J. Inst. Met. 1938, 62, 307–324. [Google Scholar]
  54. Sachs, G. Plasticity problems in metals. Trans. Faraday Soc. 1928, 24, 84–92. [Google Scholar] [CrossRef]
  55. Montheillet, F.; Damamme, G. Simple Flow Rules for Modeling the Behavior of Inhomogeneous Viscoplastic Materials. Adv. Eng. Mater. 2005, 7, 852–858. [Google Scholar] [CrossRef]
  56. Sun, J.-Z.; Li, M.-Q.; Li, H. Fragmentation of α Grains Accelerated by the Growth of β Phase in Ti–5Al–2Sn–2Zr–4Mo–4Cr during Hot Deformation. Adv. Eng. Mater. 2018, 20, 1700200. [Google Scholar] [CrossRef]
  57. Sun, J.Z.; Li, M.Q.; Li, H. Interaction effect between alpha and beta phases based on dynamic recrystallization of isothermally compressed Ti-5Al-2Sn-2Zr-4Mo-4Cr with basketweave microstructure. J. Alloys Compd. 2017, 692, 403–412. [Google Scholar] [CrossRef]
  58. Sonderegger, B.; Kozeschnik, E. Particle strengthening in fcc crystals with prolate and oblate precipitates. Scr. Mater. 2012, 66, 52–55. [Google Scholar] [CrossRef]
  59. Souza, P.M.; Beladi, H.; Singh, R.P.; Hodgson, P.D.; Rolfe, B. An Analysis on the Constitutive Models for Forging of Ti6Al4V Alloy Considering the Softening Behavior. J. Mater. Eng. Perform. 2018, 27, 3545–3558. [Google Scholar] [CrossRef]
  60. Semiatin, S.L.; Corbett, M.W.; Fagin, P.N.; Salishchev, G.A.; Lee, C.S. Dynamic-coarsening behavior of an α/β titanium alloy. Metall. Mater. Trans. A 2006, 37, 1125–1136. [Google Scholar] [CrossRef]
  61. Marsh, S.P.; Glicksman, M.E. Ostwald ripening in non-spherical morphologies. Mater. Sci. Eng. A 1997, 238, 140–147. [Google Scholar] [CrossRef]
Figure 1. Microstructures before deformation of the (a,b) Ti-5553 and (c,d) Ti-17 after heat treatment at (a) 820 °C for 15 min, (b) 920 °C for 5 min, (c) 810 °C for 15 min, and (d) 970 °C for 15 min. Only the β-phase is indexed in the electron backscattered difraction (EBSD) inverse pole figure maps, and the black lines are HAGBs.
Figure 1. Microstructures before deformation of the (a,b) Ti-5553 and (c,d) Ti-17 after heat treatment at (a) 820 °C for 15 min, (b) 920 °C for 5 min, (c) 810 °C for 15 min, and (d) 970 °C for 15 min. Only the β-phase is indexed in the electron backscattered difraction (EBSD) inverse pole figure maps, and the black lines are HAGBs.
Materials 13 05678 g001
Figure 2. Microstructures of the (a,b) Ti-5553 and (c,d) Ti-17. Microstructures after heat treatment for 15 min at (a) 820 °C, (c) 810 °C. Microstructures after deformation at 0.001 s−1 and (b) 820 °C, (d) 810 °C. The α- and β-phases are the darker and lighter phases, respectively.
Figure 2. Microstructures of the (a,b) Ti-5553 and (c,d) Ti-17. Microstructures after heat treatment for 15 min at (a) 820 °C, (c) 810 °C. Microstructures after deformation at 0.001 s−1 and (b) 820 °C, (d) 810 °C. The α- and β-phases are the darker and lighter phases, respectively.
Materials 13 05678 g002
Figure 3. Microstructures of the (a,b) Ti-5553 and (c,d) Ti-17 after deformation at 0.001 s−1 and (a) 820 °C, (b) 920 °C, (c) 810 °C, (d) 970 °C. Only the β-phase is indexed in the EBSD inverse pole figure maps, and the black lines are HAGBs.
Figure 3. Microstructures of the (a,b) Ti-5553 and (c,d) Ti-17 after deformation at 0.001 s−1 and (a) 820 °C, (b) 920 °C, (c) 810 °C, (d) 970 °C. Only the β-phase is indexed in the EBSD inverse pole figure maps, and the black lines are HAGBs.
Materials 13 05678 g003
Figure 4. Typical microstructures of the Ti-17 deformed at 0.001 s−1 and (ac) 810 °C, (d) 970 °C; representing the process of (a,b) dynamic α-globularization; (c,d) cDRX.
Figure 4. Typical microstructures of the Ti-17 deformed at 0.001 s−1 and (ac) 810 °C, (d) 970 °C; representing the process of (a,b) dynamic α-globularization; (c,d) cDRX.
Materials 13 05678 g004
Figure 5. Simulated flow curves using the “Model KM” and the “Model G” compared to the measured ones for (a,b) Ti-5553 and (c,d) Ti-17 at (a) 800 °C, (b) 920 °C, (c) 810 °C, and (d) 970 °C.
Figure 5. Simulated flow curves using the “Model KM” and the “Model G” compared to the measured ones for (a,b) Ti-5553 and (c,d) Ti-17 at (a) 800 °C, (b) 920 °C, (c) 810 °C, and (d) 970 °C.
Materials 13 05678 g005
Figure 6. Simulated stress components ( σ a t h and σ t h ) and total stress ( σ ) for the Ti-17 deformed at 970 °C and (a) 0.001 s−1, (b) 10 s−1 up to large strains.
Figure 6. Simulated stress components ( σ a t h and σ t h ) and total stress ( σ ) for the Ti-17 deformed at 970 °C and (a) 0.001 s−1, (b) 10 s−1 up to large strains.
Materials 13 05678 g006
Figure 7. Derivative of the total stress with respect to the strain for the Ti-17 deformed at 970 °C with the strain rates of 0.001 s−1 and 10 s−1 and for ”Model KM” and “Model G”.
Figure 7. Derivative of the total stress with respect to the strain for the Ti-17 deformed at 970 °C with the strain rates of 0.001 s−1 and 10 s−1 and for ”Model KM” and “Model G”.
Materials 13 05678 g007
Figure 8. Evolutions of the α-phase for the Ti-17 deformed at 810 °C and 0.001 s−1 predicted up to large strains using “Model KM” and “Model G”: (a) plastic deformation, (b) wall dislocation density, (c) α-aspect ratio, (d) fraction of HAGB, (e) α-globularized fraction, (f) fraction of Rayleigh distribution.
Figure 8. Evolutions of the α-phase for the Ti-17 deformed at 810 °C and 0.001 s−1 predicted up to large strains using “Model KM” and “Model G”: (a) plastic deformation, (b) wall dislocation density, (c) α-aspect ratio, (d) fraction of HAGB, (e) α-globularized fraction, (f) fraction of Rayleigh distribution.
Materials 13 05678 g008
Figure 9. Simulated β-grain and β-subgrain sizes for the Ti-17 deformed at 970 °C predicted using “Model KM” and “Model G”.
Figure 9. Simulated β-grain and β-subgrain sizes for the Ti-17 deformed at 970 °C predicted using “Model KM” and “Model G”.
Materials 13 05678 g009
Figure 10. Simulated evolutions of (a,d) immobile dislocation density; (b,e) wall dislocation density; (c,f) fraction of HAGBs for the (ac) Ti-5553 deformed at 920 °C and (df) Ti-17 deformed at 970 °C.
Figure 10. Simulated evolutions of (a,d) immobile dislocation density; (b,e) wall dislocation density; (c,f) fraction of HAGBs for the (ac) Ti-5553 deformed at 920 °C and (df) Ti-17 deformed at 970 °C.
Materials 13 05678 g010
Figure 11. Evolutions of the mobile dislocation density, glide velocity and plastic strain rate predicted by “Model G” for Ti-17 deformed at 970 °C and (a) 0.001 s−1, (b) 10 s−1; (c) steady-state mobile dislocation density predicted for the of the Ti-17 using “Model G” for different strain rates and temperatures and the comparison with the assumed invariable mobile dislocation density of 1013 m−2 for the “Model KM”.
Figure 11. Evolutions of the mobile dislocation density, glide velocity and plastic strain rate predicted by “Model G” for Ti-17 deformed at 970 °C and (a) 0.001 s−1, (b) 10 s−1; (c) steady-state mobile dislocation density predicted for the of the Ti-17 using “Model G” for different strain rates and temperatures and the comparison with the assumed invariable mobile dislocation density of 1013 m−2 for the “Model KM”.
Materials 13 05678 g011
Figure 12. Evolution of the dislocation density reactions for “Model G” for the Ti-17 deformed at 970 °C and (a,c) 0.001 s−1; and (b,d) 10 s−1. (a,b) mobile dislocation density reactions; (c,d) immobile dislocation density reactions. The production terms are indicated with a “+” while the deduction ones with a “−“.
Figure 12. Evolution of the dislocation density reactions for “Model G” for the Ti-17 deformed at 970 °C and (a,c) 0.001 s−1; and (b,d) 10 s−1. (a,b) mobile dislocation density reactions; (c,d) immobile dislocation density reactions. The production terms are indicated with a “+” while the deduction ones with a “−“.
Materials 13 05678 g012
Figure 13. Evolution of the global elastic strain rate predicted by “Model G” for the Ti-17 deformed at (a) 810 °C and 0.001 s−1, and (b) 970 °C. The evolution of the aspect ratio and α-thickness are shown in (a) for the initial α-thicknesses of 0.3 µm and 5 µm.
Figure 13. Evolution of the global elastic strain rate predicted by “Model G” for the Ti-17 deformed at (a) 810 °C and 0.001 s−1, and (b) 970 °C. The evolution of the aspect ratio and α-thickness are shown in (a) for the initial α-thicknesses of 0.3 µm and 5 µm.
Materials 13 05678 g013
Figure 14. Contour maps of Ti-17 deformed at 830 °C and 0.01 s−1 showing the influence of the previous β-grain size on (a,b) the immobile dislocation density and the (c,d) fraction of HAGB, predicted by (a,c) “Model KM” and (b,d) “Model G”.
Figure 14. Contour maps of Ti-17 deformed at 830 °C and 0.01 s−1 showing the influence of the previous β-grain size on (a,b) the immobile dislocation density and the (c,d) fraction of HAGB, predicted by (a,c) “Model KM” and (b,d) “Model G”.
Materials 13 05678 g014
Figure 15. Measured and simulated subgrain size (a,c) and fraction of HAGBs (b,d) for the Ti-5553 (a,b) and Ti-17 (c,d).
Figure 15. Measured and simulated subgrain size (a,c) and fraction of HAGBs (b,d) for the Ti-5553 (a,b) and Ti-17 (c,d).
Materials 13 05678 g015
Table 1. Chemical composition of the investigated Ti-5553 and Ti-17 alloys (wt. %).
Table 1. Chemical composition of the investigated Ti-5553 and Ti-17 alloys (wt. %).
AlloyAlVMoCrSnZrFeONC
Ti-55535.625.034.932.99--0.350.1170.00940.0075
Ti-174.960.013.913.681.951.910.080.1120.050.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Buzolin, R.H.; Miller Branco Ferraz, F.; Lasnik, M.; Krumphals, A.; Poletti, M.C. Improved Predictability of Microstructure Evolution during Hot Deformation of Titanium Alloys. Materials 2020, 13, 5678. https://doi.org/10.3390/ma13245678

AMA Style

Buzolin RH, Miller Branco Ferraz F, Lasnik M, Krumphals A, Poletti MC. Improved Predictability of Microstructure Evolution during Hot Deformation of Titanium Alloys. Materials. 2020; 13(24):5678. https://doi.org/10.3390/ma13245678

Chicago/Turabian Style

Buzolin, Ricardo Henrique, Franz Miller Branco Ferraz, Michael Lasnik, Alfred Krumphals, and Maria Cecilia Poletti. 2020. "Improved Predictability of Microstructure Evolution during Hot Deformation of Titanium Alloys" Materials 13, no. 24: 5678. https://doi.org/10.3390/ma13245678

APA Style

Buzolin, R. H., Miller Branco Ferraz, F., Lasnik, M., Krumphals, A., & Poletti, M. C. (2020). Improved Predictability of Microstructure Evolution during Hot Deformation of Titanium Alloys. Materials, 13(24), 5678. https://doi.org/10.3390/ma13245678

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

Article Metrics

Back to TopTop