Devising Strain Hardening Models Using Kocks–mecking Plots—a Comparison of Model Development for Titanium Aluminides and Case Hardening Steel

The present study focuses on the development of strain hardening models taking into account the peculiarities of titanium aluminides. In comparison to steels, whose behavior has been studied extensively in the past, titanium aluminides possess a much larger initial work hardening rate, a sharp peak stress and pronounced softening. The work hardening behavior of a TNB-V4 (Ti–44.5Al–6.25Nb–0.8Mo–0.1B) alloy is studied using isothermal hot compression tests conducted on a Gleeble 3500 simulator, and compared to the typical case hardening steel 25MoCrS4. The behavior is analyzed with the help of the Kocks-Mecking plots. In contrast to steel the TNB-V4 alloy shows a non-linear course of θ (i.e., no stage-III hardening) initially and exhibits neither a plateau (stage IV hardening) nor an inflection point at all deformation conditions. The present paper describes the development and application of a methodology for the design of strain hardening models for the TNB-V4 alloy and the 25CrMoS4 steel by taking the course of the Kocks-Mecking plots into account. Both models use different approaches for the hardening and softening mechanisms and accurately predict the flow stress over a wide range of deformation conditions. The methodology may hence assist in further developments of more sophisticated physically-based strain hardening models for TiAl-alloys.


Introduction
Intermetallic titanium aluminides (TiAl) alloys are relatively new class of lightweight materials which can be applied in automotive applications as well as in aviation gas turbines where they have to withstand temperatures up to 800 • C. The high potential of TiAl for these applications stems from the fact that the TiAl-based alloys have remarkable mechanical properties, such as a low specific weight of about 4 g/cm 3 , a high strength at elevated temperatures and a very good oxidation/corrosion resistance [1][2][3].In recent years, forgeable β-solidifying TiAl alloys have been developed with the goal of exploiting the benefits of the hot-worked microstructure, which allows for better performance than the as-cast state.However, as many other intermetallic alloys titanium aluminides are brittle, even at elevated temperatures.
An important task in processing TiAl to products consists in overcoming problems associated with the limited workability of TiAl-based alloys and the robust design of forging processes.With the rapid development of computing techniques, Finite Element (FE) simulation is widely applied nowadays to study metal forming processes.In metal forming, the accuracy of FE simulations depends crucially on the constitutive model.It is hence very important to establish a precise constitutive model for the flow stress which predicts the dependence of flow stress on deformation temperature, strain rate, and strain.In the past, various analytical, phenomenological [4][5][6] and empirical models [7][8][9] have been devised to model the strain hardening behavior of ordinary single phase steels.For steels, hardening is usually modeled using semi-empirical material models.For softening caused by dynamic recrystallization, the Johnson-Avrami-Mehl-Kolmogorov (JMAK) equation is commonly used.
The existing models have been successfully applied to steel and also transferred to TiAl.Extensive studies have been carried out to describe the flow stress as a function of temperature, strain rate and strain of the TiAl-alloys.Cheng et al. [10] adopt the model proposed by Laasraoui and Jonas [11] to TiAl-alloys with a nominal composition of Ti-42Al-8Nb-0.2W-0.1Y(at %).Godor et al. [12] applied two phenomenological constitutive models, the Sellars-McTegart model [13] and the Hensel-Spittel model [14], for the description of the deformation behavior of two γ-TiAl alloys with nominal compositions Ti-41Al-3Mo-0.5Si-0.1Band Ti-45Al-3Mo-0.5Si-0.1B.He et al. [15] and Pu et al. [16] established a constitutive model based on Arrhenius-type equations and the sine hyperbolic description for the Zener-Hollomon parameter for a Ti-45Al-8.5Nb-(W,B,Y)alloy and for a TiAl-Cr-V alloy, respectively.Deng et al. [17] proposed a flow stress model of Ti-47Al-alloy based on a regression analysis using Hensel-Spittel model.All existing descriptions of the flow behavior of TiAl-alloys are based on empirical models, which were derived primarily for steels.Steels are typically hot-worked in the austenitic phase, i.e., the models are designed for microstructure evolution in a single phase material.
In contrast to steel, TiAl-alloys have a complex multi-phase microstructure during hot working.The understanding of the flow behavior of TiAl alloys under hot deformation conditions thus has a great importance for the design of hot working processes.The majority of studies show that the flow curves of TiAl-based intermetallic alloys exhibited a sharp peak stress at a relatively low strain.Then, the flow stress decreased rapidly into a steady-state stress.Previous research [18][19][20][21] reported that dynamic recrystallization contributes to flow softening.Several other influencing factors have been reported, involving most notably flow localization that is partly initiated by dynamic recrystallization and adiabatic heating during deformation [21].Hot compression tests with various TiAl-alloys with different primary structures have shown that the primary structure strongly influences the material behavior during hot forming [22].Semiatin et al. [23,24] and Schaden et al. [25] concluded that the presence of lamellar colonies is responsible for a strong softening during hot forming (i.e., the pronounced maximum in the course of the flow curves).
It is known that the strain hardening rate θ = ∂σ/∂ε as a function of flow stress (this type of plot is often referred to as Kocks-Mecking plot) reveals different stages of microstructure evolution during hot deformation as characteristic stages [26].Ordinary single phase alloys such as steels in the austenitic phase exhibit at first a linearly decreasing strain hardening rate (stage III hardening) followed by a slowly decreasing hardening rate called stage IV, before attaining stage V where the values of θ approach zero [4,27].Steels, which are hot-compression-tested in the austenitic phase, show inflection points in Kocks-Mecking plots that correspond to the onset of DRX [28].In contrast to the ones of steel, the Kocks-Mecking plots of TiAl-alloys initially decrease monotonically, but not with a linear rate.Recent research of the authors has shown that TiAl-alloys initially show a non-linear, concave down course of θ (i.e., no stage-III hardening) and exhibit neither a plateau (stage IV hardening) nor an inflection point under all tested conditions [29].
The aim of this study is to investigate the peculiarities of the work hardening behavior of a Ti-44.5Al-6.25Nb-0.8Mo-0.1B(in at %) TiAl-alloy by direct comparison to a typical 25MoCrS4 steel used for forging by means of isothermal hot compression tests.It is shown how the features of the TiAl-alloy affect the set-up of strain hardening models.Special attention is paid on the accuracy of the strain hardening model in θ-σ-space, which is typically not taken into account in model development.
The paper is structured as follows: Section 2 gives an overview of the TNB-V4 alloy and 25MoCrS4 steel and the performed hot compression tests.Section 3 details the flow curves, Kocks-Mecking plots and metallographic results.The methodology for designing strain hardening models based on the behavior observed in the Kocks-Mecking plots is presented in Section 4, which also details the parameter identification procedure.In Sections 5 and 6, a comparison of the model results to experimental data and final conclusions are presented.

Materials and Hot Compression Testing
We studied a high Nb containing TiAl alloy with a nominal composition of Ti-44.5Al-6.25Nb-0.8Mo-0.1B(in at %) and a 25MoCrS4 steel with the chemical composition of 0.25 wt % C, 0.25% Si, 0.70% Mn, 1.05% Cr, 0.25% Mo, balance Fe.In order to start with chemically homogeneous samples the TiAl alloy was prepared by VAR melting.Residual porosity was removed by hot isostatic pressing at 1260 • C for 4 h with a pressure of 200 bars.A Gleeble 3500 thermo-mechanical simulator (DSI Europe GmbH, Weissenhorn, Germany) was used to simulate the hot forging process by means of hot compression tests.Strain measurement was performed with the standard LVDT (Linear variable differential transformer) of the Gleeble machine (DSI Europe GmbH, Weissenhorn, Germany).The testing temperatures ranged from 1150 • C to 1300 • C and the strain rates ranged from 0.001 s −1 to 0.5 s −1 .Compression tests were conducted in vacuum of 10 −4 mbar.The specimens were heated with a heating rate of 10 K/s and held for 5 min to obtain a homogeneous microstructure before compression testing.
In order to investigate the microstructure evolution under hot working conditions, each specimen was quenched immediately after upsetting.Metallographic examination was carried out on cross sections parallel to the compression direction of the deformed specimens.The sections were also analyzed using scanning electron microscopy (Tescan Orsay Holding, Brno-Kohoutovice, Česká Republika).
The 25MoCrS4 steel was analyzed by compression testing in previous work of the first author of this paper [30].Cylindrical specimens with a diameter of 10 mm and a height of 15 mm were deformed at temperatures ranging from 700 • C to 1200 • C, and strain rates from 0.01 s −1 to 100 s −1 .In all tests, the specimens were first heated from room temperature to 1250 • C and homogenized for 10 min, then placed into the furnace of the compression testing machine (Servotest Testing Systems Ltd., Surrey, UK) and cooled to the specific deformation temperature.The maximum true strain obtained in all tests was ~0.8.All specimens were quenched after testing and examined for grain size and recrystallized volume fractions.

Hot Deformation Behavior
Figure 1 illustrates the true stress-strain curves at deformation temperature of 1200 • C and different strain rates in the isothermal compression of the 25MoCrS4 steel (Figure 1a) and the TNB-V4-alloy (Figure 1b).It can be seen that the hot deformation behavior of the TNB alloy is in strong contrast to the behavior of the 25MoCrS4 steel, which is fully austenitic under the testing conditions.
It is observed that the isothermal compression behavior of the 25MoCrS4 steel is approximately divided into three stages: (i) The flow stress increases to a peak value (σ p ) with increasing strain; (ii) The flow stress decreases to a steady state value.In this stage the softening becomes increasingly predominant due to the onset of DRX; (iii) The flow stress keeps a steady state due to the dynamic balance between work hardening and softening induced by DRX.It is apparent that the flow curves of steel display the typical characteristic of softening by DRX.Typical Kocks-Mecking plots (θ-σ curves) of the 25MoCrS4 steel at 1200 • C and various strain rates are shown in Figure 1c.The steel exhibits a characteristic work-hardening behavior.At first, θ decreases linearly (so called stage-III hardening) until reaching an inflection point, which determines the critical stress (σ c ) and accordingly the onset of DRX [28].Just before the inflection point, the working hardening rate decreases to a lower slope (stage IV hardening).After the onset of DRX, θ decreases rapidly towards the peak stress (σ p ), at which θ = 0.The θ-σ curves can be used to determine the onset of DRX.
hardening) until reaching an inflection point, which determines the critical stress (σc) and accordingly the onset of DRX [28].Just before the inflection point, the working hardening rate decreases to a lower slope (stage IV hardening).After the onset of DRX, θ decreases rapidly towards the peak stress (σp), at which θ = 0.The θ-σ curves can be used to determine the onset of DRX.In contrast to the flow curves of the 25MoCrS4 steel, the TNB-V4 alloy has a high initial hardening with a pronounced flow stress maximum, which is much sharper than for steel in the austenitic range.Then, the flow stress decreases rapidly to a steady-state stress with increasing strain, albeit with a different curvature than steel.The measured flow curves of TNB-V4 show a very strong dependence on strain rate and temperature.The flow stress increases very rapidly with increasing strain rate and decreasing temperature.In contrast to steel, it can be seen (Figure 1d) that the Kocks-Mecking plots of TNB-V4 look absolutely different.The alloy under investigation neither shows a linear initial decrease of θ (i.e., no stage-III hardening) nor does it exhibit a plateau (stage IV hardening) or an inflection point at all forming conditions.The work-hardening behavior of the TNB-V4-alloy can be explained by the complex multiphase microstructure, and many factors that may affect the hot deformation behavior of TiAl-alloys, such as flow localization, additional interior adiabatic heating, rotation or buckling and breakdown of the lamellar colonies during hot forming, cf.[31].
It has to be mentioned here that measurement accuracy, adiabatic heating, inhomogeneous deformation as well as further processing of the recorded data will affect the computation of stress-strain curves and the derivatives of the stress strain curves, which form the basis of the Kocks-Mecking plots used in this paper.In previous work by Lohmar and Bambach [32], the processing of recorded data was treated in detail with respect to the computation of Kocks-Mecking plots.The effect and correction of adiabatic heating was treated in Xiong et al. [33].Data processing in the present paper was performed according to these findings.The stress-strain curves and Kocks-Mecking plots should hence well represent the behavior of the material.In contrast to the flow curves of the 25MoCrS4 steel, the TNB-V4 alloy has a high initial hardening with a pronounced flow stress maximum, which is much sharper than for steel in the austenitic range.Then, the flow stress decreases rapidly to a steady-state stress with increasing strain, albeit with a different curvature than steel.The measured flow curves of TNB-V4 show a very strong dependence on strain rate and temperature.The flow stress increases very rapidly with increasing strain rate and decreasing temperature.In contrast to steel, it can be seen (Figure 1d) that the Kocks-Mecking plots of TNB-V4 look absolutely different.The alloy under investigation neither shows a linear initial decrease of θ (i.e., no stage-III hardening) nor does it exhibit a plateau (stage IV hardening) or an inflection point at all forming conditions.The work-hardening behavior of the TNB-V4-alloy can be explained by the complex multiphase microstructure, and many factors that may affect the hot deformation behavior of TiAl-alloys, such as flow localization, additional interior adiabatic heating, rotation or buckling and breakdown of the lamellar colonies during hot forming, cf.[31].

Microstructure Evolution
It has to be mentioned here that measurement accuracy, adiabatic heating, inhomogeneous deformation as well as further processing of the recorded data will affect the computation of stress-strain curves and the derivatives of the stress strain curves, which form the basis of the Kocks-Mecking plots used in this paper.In previous work by Lohmar and Bambach [32], the processing of recorded data was treated in detail with respect to the computation of Kocks-Mecking plots.The effect and correction of adiabatic heating was treated in Xiong et al. [33].Data processing in the present paper was performed according to these findings.The stress-strain curves and Kocks-Mecking plots should hence well represent the behavior of the material.

Microstructure Evolution
The microstructure of the latest industrially used generation of Titanium aluminides consists of globular γ-grains, α + γ-lamellar colonies and an increased amount of β/β 0 -grains along the boundaries of lamellar colonies.Compared with two phase alloys, the stabilization of the β/β 0 -phase leads to a better strength-ductility ratio.Especially the workability of these TiAl-alloys is better than the one of previous generations of TiAl-alloys [34].The deformation in two-phase α 2 + γalloys is mainly based on dislocation glide and twinning in the γ-phase combined with deformation anisotropy of the lamellar colonies.The anisotropy of the lamellar colonies leads to soft or hard deformation modes depending on the direction of the applied stress.The yield and fracture stresses are low if the deformation occurs parallel to the lamella plane.In the opposite case the stresses are high if the deformation occurs perpendicular to the lamellae [35].Under load hard mode, colonies rotate into flow direction or a kinking of the lamellas occurs.As can be seen in Figure 2b, after kinking both parts of the colonies subsequently rotate into flow direction during forming.It can be assumed that, depending on the bcc structure of the β-phase and its ordered counterpart β 0 , both phases have a sufficient amount of independent slip modes.That leads, combined with the limited deformability of the αand γ-phase, to a high deformation gradient between the different phases.As a consequence, the β/β 0 -phase acts as a lubricant for the whole microstructure and supports the rotation of hard mode lamellar α 2 + γ-colonies [36].It is hence assumed that the rotation of the α 2 + γ-colonies from hard into soft orientations leads to a non-linear course of strain hardening rate in the Kocks-Mecking plots.The microstructure of the latest industrially used generation of Titanium aluminides consists of globular γ-grains, α + γ-lamellar colonies and an increased amount of β/β0-grains along the boundaries of lamellar colonies.Compared with two phase alloys, the stabilization of the β/β0-phase leads to a better strength-ductility ratio.Especially the workability of these TiAl-alloys is better than the one of previous generations of TiAl-alloys [34].The deformation in two-phase α2 + γ-alloys is mainly based on dislocation glide and twinning in the γ-phase combined with deformation anisotropy of the lamellar colonies.The anisotropy of the lamellar colonies leads to soft or hard deformation modes depending on the direction of the applied stress.The yield and fracture stresses are low if the deformation occurs parallel to the lamella plane.In the opposite case the stresses are high if the deformation occurs perpendicular to the lamellae [35].Under load hard mode, colonies rotate into flow direction or a kinking of the lamellas occurs.As can be seen in Figure 2b, after kinking both parts of the colonies subsequently rotate into flow direction during forming.It can be assumed that, depending on the bcc structure of the β-phase and its ordered counterpart β0, both phases have a sufficient amount of independent slip modes.That leads, combined with the limited deformability of the α-and γ-phase, to a high deformation gradient between the different phases.As a consequence, the β/β0-phase acts as a lubricant for the whole microstructure and supports the rotation of hard mode lamellar α2 + γ-colonies [36].It is hence assumed that the rotation of the α2 + γ-colonies from hard into soft orientations leads to a non-linear course of strain hardening rate in the Kocks-Mecking plots.

Derivation of Model Equations
Various microstructure-based strain hardening models have been derived in the past (see, e.g., Bariani et al. [37] for a review).For steel, dynamic recrystallization (DRX) is a dominant softening mechanism which restores the deformed microstructure during forming and allows for large plastic deformations without failure of the workpiece.Most strain hardening models rely on Johnson-Mehl-Avrami-Kolmogorov (JMAK) kinetics for DRX.Poliak and Jonas [28] showed that the onset of DRX corresponds to an inflection point in the course of the strain hardening rate θ as a function of stress σ, i.e., the critical stress is a root of the second derivative of the strain hardening rate with respect to stress: Recently, two problems in modeling the strain hardening behavior under hot working conditions were discovered: (i) inconsistency of JMAK kinetics with the Poliak-Jonas criterion; and (ii) absence of an inflection point in practical experiments despite occurrence of DRX [38].

Derivation of Model Equations
Various microstructure-based strain hardening models have been derived in the past (see, e.g., Bariani et al. [37] for a review).For steel, dynamic recrystallization (DRX) is a dominant softening mechanism which restores the deformed microstructure during forming and allows for large plastic deformations without failure of the workpiece.Most strain hardening models rely on Johnson-Mehl-Avrami-Kolmogorov (JMAK) kinetics for DRX.Poliak and Jonas [28] showed that the onset of DRX corresponds to an inflection point in the course of the strain hardening rate θ as a function of stress σ, i.e., the critical stress is a root of the second derivative of the strain hardening rate with respect to stress: ∂ ∂σ Recently, two problems in modeling the strain hardening behavior under hot working conditions were discovered: (i) inconsistency of JMAK kinetics with the Poliak-Jonas criterion; and (ii) absence of an inflection point in practical experiments despite occurrence of DRX [38].
Since the criterion by Poliak and Jonas was derived from principles of irreversible thermodynamics and has been confirmed with experimental data, a flow stress model that takes DRX into account should adhere to this criterion, not only to be consistent with experimental data but also to make sure that the model is sound in a thermodynamical sense.Hence, the flow stress σ should be three times continuously differentiable with respect to ε.It was recently shown by Bambach [38] that the Poliak-Jonas criterion requires that the Avrami exponent exceeds a value of 3 for the model to be consistent with the criterion.Hence, generalized transformation kinetics was considered to find the root cause of the inconsistency and to alleviate them.As mentioned in the introduction, TiAl alloys do not seem to show an inflection point in θ in spite of the occurrence of DRX [29].
A flow stress model that includes DRX essentially consists of four types of equations: (i) a model for strain hardening and dynamic recovery of the non-recrystallized material; (ii) an initiation criterion for DRX; (iii) DRX kinetics and an equation for the recrystallized grain size; and (iv) a rule of mixture to determine the macroscopic flow stress when recrystallized and non-recrystallized grains coexist.In addition, the evolution of grain size is expressed by an additional equation, which is mostly not coupled to the flow stress model.
Various model choices for the individual parts of the model are available.Due to the large differences in strain hardening behavior for different materials, setting up a model for a certain material requires to take the peculiarities of the strain hardening and softening of the material into account.As a consequence, there does not seem to be a unified model that would be applicable to all kinds of materials.In the following, we describe the features of typical strain hardening models that take DRX into account.We then show how to devise and select the model equations so that a suitable model is attained, at the example of the case hardening steel and the TNB alloy.
Model equations for characteristic points.Semi-empirical material models are the most commonly used models in industrial applications.They are often preferred to physically-based models due to their fast computing times.The basic idea behind semi-empirical models is that characteristic points of the flow curves can be expressed as a function of the Zener-Hollomon parameter Z.A typical semi-empirical model based on the equations proposed by Luton and Sellars [39] and Beynon and Sellars [40] is given in Figure 3 together with an illustration of a flow curve, the strain hardening rate and their characteristic points.
Since the criterion by Poliak and Jonas was derived from principles of irreversible thermodynamics and has been confirmed with experimental data, a flow stress model that takes DRX into account should adhere to this criterion, not only to be consistent with experimental data but also to make sure that the model is sound in a thermodynamical sense.Hence, the flow stress σ should be three times continuously differentiable with respect to ε.It was recently shown by Bambach [38] that the Poliak-Jonas criterion requires that the Avrami exponent exceeds a value of 3 for the model to be consistent with the criterion.Hence, generalized transformation kinetics was considered to find the root cause of the inconsistency and to alleviate them.As mentioned in the introduction, TiAl alloys do not seem to show an inflection point in θ in spite of the occurrence of DRX [29].
A flow stress model that includes DRX essentially consists of four types of equations: (i) a model for strain hardening and dynamic recovery of the non-recrystallized material; (ii) an initiation criterion for DRX; (iii) DRX kinetics and an equation for the recrystallized grain size; and (iv) a rule of mixture to determine the macroscopic flow stress when recrystallized and non-recrystallized grains coexist.In addition, the evolution of grain size is expressed by an additional equation, which is mostly not coupled to the flow stress model.
Various model choices for the individual parts of the model are available.Due to the large differences in strain hardening behavior for different materials, setting up a model for a certain material requires to take the peculiarities of the strain hardening and softening of the material into account.As a consequence, there does not seem to be a unified model that would be applicable to all kinds of materials.In the following, we describe the features of typical strain hardening models that take DRX into account.We then show how to devise and select the model equations so that a suitable model is attained, at the example of the case hardening steel and the TNB alloy.
Model equations for characteristic points.Semi-empirical material models are the most commonly used models in industrial applications.They are often preferred to physically-based models due to their fast computing times.The basic idea behind semi-empirical models is that characteristic points of the flow curves can be expressed as a function of the Zener-Hollomon parameter Z.A typical semi-empirical model based on the equations proposed by Luton and Sellars [39] and Beynon and Sellars [40] is given in Figure 3 together with an illustration of a flow curve, the strain hardening rate and their characteristic points.The strain at which the maximum value of the flow stress is attained is given by Equation ( 2).The critical strain for the initiation of DRX is assumed to be a constant fraction a4 of εp by Equation (3).The strain at which a steady state deformation is reached is given by Equation ( 4), which uses four fit parameters, e1-e4.The peak and steady state flow stress σp and σss are modeled as functions The strain at which the maximum value of the flow stress is attained is given by Equation (2).The critical strain for the initiation of DRX is assumed to be a constant fraction a 4 of ε p by Equation (3).The strain at which a steady state deformation is reached is given by Equation ( 4), which uses four fit parameters, e 1 -e 4 .The peak and steady state flow stress σ p and σ ss are modeled as functions of the Zener-Hollomon parameter Z via Equations ( 5) and ( 6), with material-dependent parameters f 1 , f 2 , f 3 and h 1 , h 2 , h 3 .
Strain hardening.A plenitude of models is available to model strain hardening.Recent research has shown that different kinds of inconsistencies with the Poliak-Jonas criterion can occur depending on the model choice [38].If the strain hardening models do not show an inflection point, they do not reach the critical point.If the strain hardening model shows an inflection point, care must be taken that the inflection point coincides with the critical conditions of DRX nucleation models.In addition, the model might just not capture the course of the strain hardening rate in a Kocks-Mecking plot in the right way, which shows that the underlying physics are not taken into account.
Typical dislocation density models are formulated as ordinary differential equations, which contain terms for the generation and annihilation of dislocations.These models were derived from physical mechanisms of dislocation theory.The Taylor equation correlates dislocation density to the macroscopic flow stress.It can be used to transform the model equations into θ-σ-space (Kocks-Mecking plots), where the course of the strain hardening rate can be analyzed.None of the dislocation-based models detailed in Table 1 are appropriate for the TNB alloy since they cannot reproduce the curvature in a Kocks-Mecking plot.The phenomenological model by Cingara and McQueen, however, possesses the right curvature for C ≥ 1.The constant C determines the curvature of the flow curve up to the peak.The location of the peak is explicitly taken into account by the strain hardening model.Softening will be superimposed after the onset of DRX and will be discussed below.
In contrast, the model by Cingara and McQueen is not the best choice for the 25MoCrS4 steel.The 25MoCrS4 steel shows all hardening stages, i.e., a linear decrease of θ(σ) initially, which is referred to as stage-III hardening, an stage IV hardening with a reduced softening rate and an inflection point, and a subsequently increased softening rate due to DRX.
In physical theories of crystal plasticity, the dislocation density is commonly used as the structural parameter for macroscopic descriptions of plastic flow, with the flow stress governed by Taylor-like hardening following: where σ 0 is the yield stress, at which point ρ = ρ 0 and ε = ε 0 , α is a constant of 0.5-1, G is the shear modulus, and b is the Burger's vector.The Taylor equation can be utilized to transform the model equations into Kocks-Mecking plots.Equally, models can be formulated directly in the θ-σ-space and then transferred into stress-strain curves based on the Taylor equation.The strain hardening model proposed in [38] reproduces the course of the experimentally observed θ-σ-curves using a linear function for the initial stage-III hardening, a transition to a second linear function which represents stage IV and then, starting from the critical point, a cubic function of stress is used to model stage V hardening.Hence, the model is designed to mimic the hardening behavior observed for steel (cf. Figure 1c), but it would not be applicable to TNB-V4.Recrystallization kinetics.In the original work of Luton and Sellars [39], JMAK kinetics was proposed for DRX, which were formulated as function of time.Later, a strain-dependent version of the kinetics was introduced by Sellars [41].Recent research has shown that JMAK kinetics is only consistent with the Poliak-Jonas criterion if the Avrami exponent exceeds a value of 3 [38].As a consequence, DRX kinetics based transformation kinetics that go back on Cahn were proposed, which are detailed in Table 2.In contrast, the softening of TNB is not only determined by DRX but also by the rotation and break-down of lamellar colonies, adiabatic heating etc.Since consistency to the Poliak-Jonas criterion is not an issue for TNB-V4 (it shows no inflection point), ordinary JMAK kinetics may be used, albeit with no physical meaning.
Dynamically recrystallized grain size.The dynamically recrystallized grain size is independent of the initial grain size and only depends on the Zener-Hollomon parameter via This equation is used in almost all semi-empirical models for DRX.In this work, it is applied to the 25MoCrS4 steel and to the phases of the TNB alloy that show recrystallization.
Strain hardening models.Based on the observations made about the strain hardening behavior of TNB-V4 and 25MoCrS4, two strain hardening models have been set up, which are detailed in Table 2.Both models are designed to meet the requirements of the materials Using non-linear regression by least-squares, the models have been fitted to a single flow curve of TNB-V4 and 25MoCrS4, respectively.Figure 4 shows that for a single flow curve, the models reproduce not only the flow curves but also the behavior in the Kocks-Mecking plot with very high accuracy.When the model is applied to a larger range of temperatures and strain rates, parameter identification by non-linear regression may fail.To allow for successful parameter fitting, most of the model equations are first fitted to the experimental data using the characteristic points extracted from the flow curves.Regression is thus performed separately for each model equation.In the model, all equations are coupled.Hence, a non-linear regression of the entire model is performed subsequently, which will improve the initial parameter values obtained by fitting the model equations individually.The procedure is described in detail in the subsequent sections.Dynamically recrystallized grain size.The dynamically recrystallized grain size is independent of the initial grain size and only depends on the Zener-Hollomon parameter via This equation is used in almost all semi-empirical models for DRX.In this work, it is applied to the 25MoCrS4 steel and to the phases of the TNB alloy that show recrystallization.
Strain hardening models.Based on the observations made about the strain hardening behavior of TNB-V4 and 25MoCrS4, two strain hardening models have been set up, which are detailed in Table 2.Both models are designed to meet the requirements of the materials Using non-linear regression by least-squares, the models have been fitted to a single flow curve of TNB-V4 and 25MoCrS4, respectively.Figure 4 shows that for a single flow curve, the models reproduce not only the flow curves but also the behavior in the Kocks-Mecking plot with very high accuracy.When the model is applied to a larger range of temperatures and strain rates, parameter identification by non-linear regression may fail.To allow for successful parameter fitting, most of the model equations are first fitted to the experimental data using the characteristic points extracted from the flow curves.Regression is thus performed separately for each model equation.In the model, all equations are coupled.Hence, a non-linear regression of the entire model is performed subsequently, which will improve the initial parameter values obtained by fitting the model equations individually.The procedure is described in detail in the subsequent sections.

Evaluation of Experimental Data and Parameters Identification
The procedures for parameter identification depend largely on the model equations.However, some steps are common for the two models defined in Table 2 above.As a first step, all characteristic points such as εp, σp, εss, and σss have to be determined from the experimental flow curves.Using σp,

Evaluation of Experimental Data and Parameters Identification
The procedures for parameter identification depend largely on the model equations.However, some steps are common for the two models defined in Table 2 above.As a first step, all characteristic points such as ε p , σ p , ε ss , and σ ss have to be determined from the experimental flow curves.
Using σ p , the activation energy and hence the Zener-Hollomon parameter Z, sometimes referred to as temperature-compensated strain rate, can be determined, where A and α are constants, T is the deformation temperature, σ is the flow stress (MPa), .
ε is the strain rate (s −1 ), n is the stress exponent, R is the gas constant (kJ•mol −1 •K −1 ) and Q w is the activation energy of deformation (kJ•mol −1 ).
For some model parameters (predominantly those with a physical background), dedicated procedures exist for parameter determination.The activation energy of hot working Q w , for instance, can be determined from the following relation: According to Poliak and Jonas [28], the critical strain ε cr for the onset of DRX can be detected from the inflection point in the course of the strain hardening rate θ = ∂σ/∂ε as a function of flow stress σ.To this end, the strain hardening rate has to be computed from flow stress data.All other model parameters of the equations in Figure 3 can be determined for each equation independently by regression analysis.
A specific problem with the TNB alloy is that it shows no inflection point in the Kocks-Mecking plot.Hence, no critical stress or strain can be determined from flow stress data only.
To describe the relation between the dynamically recrystallized volume fraction X and strain an Avrami-type equation is used for TNB-V4 while different kinetics are used for 25MoCrS4 to ensure consistency with the Poliak-Jonas criterion.The model parameters are not determined from the flow curves but from the recrystallized fractions determined by quantitative metallography of quenched specimens.Based on the phase fraction and grain size of the different phases in the not deformed state a histogram for every phase was calculated.These histograms were used as basis to determine the differences of the phase fractions and grain sizes after deformation.
In the same way, the coefficients b 1 and b 2 in Equation ( 8), which describes the dynamically recrystallized grain size, were identified from metallographic evaluation of the specimens that show DRX.

Determination of Model Parameters Common to Both Models
In this section, the experimental data from the cylinder compression tests and the metallographic analysis are used to determine the parameters for the strain hardening models.For this purpose, firstly the characteristic points such as ε p , σ p , ε ss , and σ ss were determined from the flow curves shown in Figure 1.
Based on Equation (10), the slope of the plots of ln[sinh(ασ)] versus ln .
ε and 1/T can be utilized for calculating the value of Q w .The required plots are shown in Figure 5. Values of 481.3 kJ•mol −1 for 25MoCrS4 and 321.2 kJ•mol −1 for TNB-V4 were obtained for the activation energy Q w .
In Figure 6a, a plot of the peak strain as a function of the Zener-Hollomon parameter is shown together with the results of regression analysis.The steady state strain strongly influences the shape of the flow curves.The extraction of values for the steady state strain is complicated for TNB-V4 since only a few flow curves attain a steady state regime.To determine the parameters e 1 , e 2 , and e 3 in Equation ( 4), which defines the dependence of ε ss on the Zener-Hollomon parameter, a regression analysis was performed as shown in Figure 6b.
Since the TNB-V4-alloy shows no inflection point in the Kocks-Mecking plots no critical stress/strain can be determined from flow stress data only and α was set to 0.5 from experience.For 25MoCrS4 previous work showed that a value of 0.7ε p is suitable.The dependence of peak and steady state stress on the Zener-Hollomon parameter was determined by non-linear regression using Equations ( 5) and ( 6).The fitted and measured values are shown in Figure 7.As mentioned above, the equation of Cingara et al. [44] was employed to describe the work-hardening condition of TNB-V4-alloy up to the peak stress.Hence, the fitted peak stress values can be used directly.The dependence of peak and steady state stress on the Zener-Hollomon parameter was determined by non-linear regression using Equations ( 5) and ( 6).The fitted and measured values are shown in Figure 7.As mentioned above, the equation of Cingara et al. [44] was employed to describe the work-hardening condition of TNB-V4-alloy up to the peak stress.Hence, the fitted peak stress values can be used directly.The dependence of peak and steady state stress on the Zener-Hollomon parameter was determined by non-linear regression using Equations ( 5) and (6).The fitted and measured values are shown in Figure 7.As mentioned above, the equation of Cingara et al. [44] was employed to describe the work-hardening condition of TNB-V4-alloy up to the peak stress.Hence, the fitted peak stress values can be used directly.
The dynamically recrystallized grain size of 25MoCrS4 increases with decreasing Zener-Hollomon parameter, i.e., by increase of temperature and decrease of strain rate the resulting grain size increases.
In Figure 8a, the grain sizes determined for various Z values are shown along with the best fit of Equation (8).The kinetics of dynamic recrystallization is displayed in Figure 8b.The dynamically recrystallized grain size of 25MoCrS4 increases with decreasing Zener-Hollomon parameter, i.e., by increase of temperature and decrease of strain rate the resulting grain size increases.In Figure 8a, the grain sizes determined for various Z values are shown along with the best fit of Equation ( 8).The kinetics of dynamic recrystallization is displayed in Figure 8b.The relationship between the dynamically recrystallized grain size of the TNB-V4-alloy and the Zener-Hollomon-Parameter is shown in Figure 9.In the multiphase TiAl-alloys, solely the β/β0-(Figure 9a) and γ-phase (Figure 9b) show recrystallization.The α-phase only forms during hot forming.The typical microstructure of the TNB-V4-alloy after forming can be seen in Figure 9c.The dynamically recrystallized grain size of 25MoCrS4 increases with decreasing Zener-Hollomon parameter, i.e., by increase of temperature and decrease of strain rate the resulting grain size increases.In Figure 8a, the grain sizes determined for various Z values are shown along with the best fit of Equation ( 8).The kinetics of dynamic recrystallization is displayed in Figure 8b.The relationship between the dynamically recrystallized grain size of the TNB-V4-alloy and the Zener-Hollomon-Parameter is shown in Figure 9.In the multiphase TiAl-alloys, solely the β/β0-(Figure 9a) and γ-phase (Figure 9b) show recrystallization.The α-phase only forms during hot forming.The typical microstructure of the TNB-V4-alloy after forming can be seen in Figure 9c.The relationship between the dynamically recrystallized grain size of the TNB-V4-alloy and the Zener-Hollomon-Parameter is shown in Figure 9.In the multiphase TiAl-alloys, solely the β/β 0 -(Figure 9a) and γ-phase (Figure 9b) show recrystallization.The α-phase only forms during hot forming.The typical microstructure of the TNB-V4-alloy after forming can be seen in Figure 9c.The dynamically recrystallized grain size of 25MoCrS4 increases with decreasing Zener-Hollomon parameter, i.e., by increase of temperature and decrease of strain rate the resulting grain size increases.In Figure 8a, the grain sizes determined for various Z values are shown along with the best fit of Equation ( 8).The kinetics of dynamic recrystallization is displayed in Figure 8b.The relationship between the dynamically recrystallized grain size of the TNB-V4-alloy and the Zener-Hollomon-Parameter is shown in Figure 9.In the multiphase TiAl-alloys, solely the β/β0-(Figure 9a) and γ-phase (Figure 9b) show recrystallization.The α-phase only forms during hot forming.The typical microstructure of the TNB-V4-alloy after forming can be seen in Figure 9c.

Final Calibration and Validation of the Model
The direct determination of model parameters detailed in section 4 is not applicable to all model parameters.In addition, since the interplay of the individual parts of the model determines the accuracy of the model, an inverse parameter identification is favorable as final step, i.e., after starting values of most model parameters are available from the direct determination performed in Section 4. Inverse parameter identification is essentially the solution of an optimization problem.The cost function is the sum of squares of the differences between measured and calculated flow stress values for all available values from the experiments: The cost function F(β) is minimized using an optimization algorithm such as the Levenberg-Marquardt method to determine the material parameters, which are assembled in the parameter vector β.
Figure 10 shows a comparison between the predicted and measured flow stress values for the 25MoCrS4 steel for temperatures ranging from 900 °C to 1200 °C and five different strain rates.The model parameters were determined in a similar way as in the previous study of Konovalov et al. [30].However, in the present work, a new model formulation was used.A very good agreement is obtained between the calculation and measurement in view of the shape of the flow curves and the absolute level of the calculated flow stress over the entire range of strain rates and temperatures.

Final Calibration and Validation of the Model
The direct determination of model parameters detailed in section 4 is not applicable to all model parameters.In addition, since the interplay of the individual parts of the model determines the accuracy of the model, an inverse parameter identification is favorable as final step, i.e., after starting values of most model parameters are available from the direct determination performed in Section 4. Inverse parameter identification is essentially the solution of an optimization problem.The cost function is the sum of squares of the differences between measured and calculated flow stress values for all available values from the experiments: The cost function F(β) is minimized using an optimization algorithm such as the Levenberg-Marquardt method to determine the material parameters, which are assembled in the parameter vector β.
Figure 10 shows a comparison between the predicted and measured flow stress values for the 25MoCrS4 steel for temperatures ranging from 900 • C to 1200 • C and five different strain rates.The model parameters were determined in a similar way as in the previous study of Konovalov et al. [30].However, in the present work, a new model formulation was used.A very good agreement is obtained between the calculation and measurement in view of the shape of the flow curves and the absolute level of the calculated flow stress over the entire range of strain rates and temperatures.

Final Calibration and Validation of the Model
The direct determination of model parameters detailed in section 4 is not applicable to all model parameters.In addition, since the interplay of the individual parts of the model determines the accuracy of the model, an inverse parameter identification is favorable as final step, i.e., after starting values of most model parameters are available from the direct determination performed in Section 4. Inverse parameter identification is essentially the solution of an optimization problem.The cost function is the sum of squares of the differences between measured and calculated flow stress values for all available values from the experiments: The cost function F(β) is minimized using an optimization algorithm such as the Levenberg-Marquardt method to determine the material parameters, which are assembled in the parameter vector β.
Figure 10 shows a comparison between the predicted and measured flow stress values for the 25MoCrS4 steel for temperatures ranging from 900 °C to 1200 °C and five different strain rates.The model parameters were determined in a similar way as in the previous study of Konovalov et al. [30].However, in the present work, a new model formulation was used.A very good agreement is obtained between the calculation and measurement in view of the shape of the flow curves and the absolute level of the calculated flow stress over the entire range of strain rates and temperatures.The model shows a very high accuracy for a wide range of temperatures and strain rates.In Figure 11, the comparison between the calculated and measured flow stresses is shown for the TNB-V4-alloy for strain rates varying from 0.001 s −1 to 0.5 s −1 and for a temperature of 1200 °C suitable for isothermal forging.The results indicate that the proposed strain hardening model can properly describe and predict the flow stress of the TiAl-alloy under isothermal forging conditions.The model accurately predicts the flow stress up to the peak stress and to some extent also the transition to the steady state.
The results show that for both 25MoCrS4 and TNB-V4, building a model from the behavior observed in Kocks-Mecking plots has the advantage that the model is consistent with the observed strain hardening stages and shows a high accuracy.If only curve fitting in terms of stress strain The model shows a very high accuracy for a wide range of temperatures and strain rates.In Figure 11, the comparison between the calculated and measured flow stresses is shown for the TNB-V4-alloy for strain rates varying from 0.001 s −1 to 0.5 s −1 and for a temperature of 1200 • C suitable for isothermal forging.The model shows a very high accuracy for a wide range of temperatures and strain rates.In Figure 11, the comparison between the calculated and measured flow stresses is shown for the TNB-V4-alloy for strain rates varying from 0.001 s −1 to 0.5 s −1 and for a temperature of 1200 °C suitable for isothermal forging.The results indicate that the proposed strain hardening model can properly describe and predict the flow stress of the TiAl-alloy under isothermal forging conditions.The model accurately predicts the flow stress up to the peak stress and to some extent also the transition to the steady state.
The results show that for both 25MoCrS4 and TNB-V4, building a model from the behavior observed in Kocks-Mecking plots has the advantage that the model is consistent with the observed strain hardening stages and shows a high accuracy.If only curve fitting in terms of stress strain The results indicate that the proposed strain hardening model can properly describe and predict the flow stress of the TiAl-alloy under isothermal forging conditions.The model accurately predicts the flow stress up to the peak stress and to some extent also the transition to the steady state.
The results show that for both 25MoCrS4 and TNB-V4, building a model from the behavior observed in Kocks-Mecking plots has the advantage that the model is consistent with the observed strain hardening stages and shows a high accuracy.If only curve fitting in terms of stress strain curves is performed and valuable information provided by the Kocks-Mecking plots is neglected, inconsistencies and larger uncertainties of the model performance may be the consequence.For the 25MoCrS4, an intricate model was designed which takes the physics of dynamic recrystallization into account.For TNB-V4, both the strain hardening before the initiation of DRX and the softening kinetics are entirely phenomenological.For TiAl alloys in general, a deeper insight into the deformation mechanism leading to the observed behavior seems necessary.Reflecting the model behavior with the experimental results not only in stress-strain space but also in θ-σ-space, as proposed in this paper, yields important information to guide the model development.

Conclusions
In the present study the peculiarities of the work hardening and softening behavior of a Ti-44.5Al-6.25Nb-0.8Mo-0.1B(in at %) TNB-V4-alloy were investigated and compared to a typical case-hardening steel 25MoCrS4.The flow curves were analyzed using Kocks-Mecking plots in order to illustrate the difference in deformation behavior between a typical steel and the high Nb containing TiAl-alloy.Specific model equations for 25MoCrS4 and Ti-44.5Al-6.25Nb-0.8Mo-0.1Bwere derived based on the course of the Kocks-Mecking plots.The empirical equation proposed by Cingara and McQueen was used to model the work-hardening region for TNB-V4, since this model is able to show a concave down course in the Kocks-Mecking plot, as observed experimentally.A new model formulation capable of showing an inflection point was used for 25MoCrS4.For the DRX kinetics, a new formulation consistent with the Poliak-Jonas criterion was applied to 25MoCrS4.For TNB, a JMAK type kinetics was used, which describes all softening mechanisms without taking into account the physical background.
The following conclusions can be drawn from the presented studies: • Kocks-Mecking plots reveal that in contrast to 25MoCrS4 steel the Ti-44.5Al-6.25Nb-0.8Mo-0.1Balloy neither shows a linear decrease of θ (i.e., no stage-III hardening) nor does it exhibit a plateau (stage IV hardening) or an inflection point (marking the onset of DRX) at all forming conditions.

•
The information obtained from Kocks-Mecking plots should be taken into account in the development of strain hardening models.Otherwise, inconsistencies, e.g., with the Poliak-Jonas criterion, may result.

•
Both models show a high accuracy.They may hence be used in finite element simulations of metal forming processes.For TNB-V4, however, the complex microstructure evolution, i.e., the recrystallization of the individual phases, the reorientation of the lamellar colonies and flow localization effects during deformation need to be taken into account in continuative work.

Figure 2 .
Figure 2. Typical microstructures of TNB-V4: (a) as-cast + HIP (hot isostatic pressing) state and (b) compressed at 1250 °C with a strain rate of 0.25 s −1 to a strain of 60%.

Figure 2 .
Figure 2. Typical microstructures of TNB-V4: (a) as-cast + HIP (hot isostatic pressing) state and (b) compressed at 1250 • C with a strain rate of 0.25 s −1 to a strain of 60%.

Figure 3 .
Figure 3. Flow curve and strain hardening rate with characteristic point and model equations for the Beynon and Sellars model [40].

Figure 3 .
Figure 3. Flow curve and strain hardening rate with characteristic point and model equations for the Beynon and Sellars model [40].

Table 1 .
Comparison of different strain hardening models and their course in Kocks-Mecking plots θ(σ).

Figure 8 .
Figure 8.(a) Relationship between DRX grain size and Zener-Hollomon-Parameter: and (b) volume fraction of DRX obtained at a temperature of 1200 °C and different strain rates.

Figure 8 .
Figure 8.(a) Relationship between DRX grain size and Zener-Hollomon-Parameter: and (b) volume fraction of DRX obtained at a temperature of 1200 °C and different strain rates.

Figure 8 .
Figure 8.(a) Relationship between DRX grain size and Zener-Hollomon-Parameter: and (b) volume fraction of DRX obtained at a temperature of 1200 • C and different strain rates.

Figure 8 .
Figure 8.(a) Relationship between DRX grain size and Zener-Hollomon-Parameter: and (b) volume fraction of DRX obtained at a temperature of 1200 °C and different strain rates.

Figure 10 .
Figure 10.Model prediction vs. experimental data of 25MoCrS4 at various temperatures und strain rates.

Figure 11 .
Figure 11.Model prediction vs. experimental data of TNB-V4-alloy at 1200 °C and various strain rates.

Figure 10 .
Figure 10.Model prediction vs. experimental data of 25MoCrS4 at various temperatures und strain rates.

Figure 10 .
Figure 10.Model prediction vs. experimental data of 25MoCrS4 at various temperatures und strain rates.

Figure 11 .
Figure 11.Model prediction vs. experimental data of TNB-V4-alloy at 1200 °C and various strain rates.

Figure 11 .
Figure 11.Model prediction vs. experimental data of TNB-V4-alloy at 1200 • C and various strain rates.

Table 1 .
Comparison of different strain hardening models and their course in Kocks-Mecking plots θ(σ).