Experimental and Numerical Investigations of the Development of Residual Stresses in Thermo-Mechanically Processed Cr-Alloyed Steel 1 . 3505

Residual stresses in components are a central issue in almost every manufacturing process, as they influence the performance of the final part. Regarding hot forming processes, there is a great potential for defining a targeted residual stress state, as many adjustment parameters, such as deformation state or temperature profile, are available that influence residual stresses. To ensure appropriate numerical modeling of residual stresses in hot forming processes, comprehensive material characterization and suitable multiscale Finite Element (FE) simulations are required. In this paper, experimental and numerical investigations of thermo-mechanically processed steel alloy 1.3505 (DIN 100Cr6) are presented that serve as a basis for further optimization of numerically modeled residual stresses. For this purpose, cylindrical upsetting tests at high temperature with subsequently cooling of the parts in the media air or water are carried out. Additionally, the process is simulated on the macroscale and compared to the results based on the experimental investigations. Therefore, the experimentally processed specimens are examined regarding the resulting microstructure, distortions, and residual stresses. For the investigation on a smaller scale, a numerical model is set up based on the state-data of the macroscopic simulation and experiments, simulating the transformation of the microstructure using phase-field theory and FE analysis on microand meso-scopic level.


Introduction
During the manufacturing of components, various residual stress states may occur in the material because of the process.In the past, there have been several investigations of the relationship between forming parameters such as deformation state or temperature profile and resulting residual stresses by experiments and simulations; see for instance [1][2][3].The work of [4] investigated process parameters, deformation state, and temperature time profile with special regard to the spring-back behavior after a hot stamping process.It is noteworthy that the residual stresses resulting from metal-forming processes are generally known to have negative effects, regarding e.g., lifetime, so that avoidance or minimization during the manufacturing process is usually targeted instead of specifically influencing component properties.For example, positive effects on the fatigue life of components due to residual compressive stresses have already been achieved by the surface treatment process of shot peening applied as an additional process step [5].The overall objective pursued in this paper is to analyze and adjust hot forming processes to achieve a favorable residual stress profile directly from the forming stage.To achieve defined material properties in the formed component, targeted cooling or heat treatment is necessary [6].For example, quenching from the forging temperature is an energy-efficient way of tempering components [7].Hot forming as a thermo-mechanically coupled process offers the potential for many adjustment parameters to influence residual stresses.As a consequence, complex material phenomena arising during the individual process steps interact and emphasize each other, decisively affecting the final component properties and making the virtual process design rather complicated, cf. Figure 1.
In general, according to [10], there are three different types of residual stresses which are explained in the following.Residual stresses of the first type refer to the complete component and can be defined as macroscopic.The difference between the mean stress of a grain and the macroscopic residual stress is regarded as residual stress of the second type, while the location-dependent deviations of the stresses within the grain from the sum of first and second type are called residual stresses of the third type, cf.[11][12][13].
The resulting mechanical properties of steels strongly depend on the physical processes, e.g., hardening and phase transformations, occurring during the material production and the manufacturing.Thus, in numerical modeling, the consideration of these aspects is crucial to obtain good reliability and prediction.
The diffusionless phase transformation, i.e., austenite to martensite, is a key aspect in the improvement of metal properties.In [14], as well as [15], macroscopic constitutive models for its consideration in a large strain plasticity setting is proposed.In the investigations carried out by [16,17], a macroscopic approach to simulate the aforementioned transformation process as well as the effect of transformation induced plasticity (TRIP) is introduced.In their research [8,18] are contributing to the Finite Element (FE) supported calculation of phase transformations regarding the anisotropic transformation strains during martensitic transformation by using a user-defined FE code based on the mathematical network of [19,20].In general, for fully coupled thermo-mechanical-metallurgical process simulations on the macroscale following [21,22], FE codes based on the theory of additive strain decomposition (ASD) are used.The total strain increment is decomposed into an elastic-, plastic-, thermal-, transformation-and a TRIP-strain component, respectively.In the work of [16], first steps towards a multiscale approach within a homogenization scheme are presented.
Based on these investigations in [23,24], the multiscale approach is extended to further phases, e.g., bainite, pearlite, and steel types, e.g., low alloy ones.All the previously mentioned works use the FE Method for the numerical simulation.A further technique for the description of the physical aspects is the Multi-Phase-Field (MPF) model, which is introduced in [25][26][27].This model is based on the consistent formulation of diffusional and surface driven phase transitions and allows the calculation of interfacial tension, multiphase thermodynamics and elastic stress balance in multiple junctions between an arbitrary number of grains and phases.Numerous publications have picked up on this model for the description of a variety of multi-physics effects, starting from single grain simulation up to the evolution of dendrites.In context of this publication, it is referred to [28][29][30], wherein the transformation process inside metals during production and manufacturing are taken into account.An implementation of the MPF model is available in the software library OpenPhase, which is distributed under the GNU General Public License (GPL) and is free of charge; see [31].Actual studies concerning martensite transformation using the phase-field theory can be found in [32][33][34][35].
In the scope of this work, an insight into the experimental and numerical analysis of the material behavior in thermo-mechanical forming processes as well as simulation approaches of residual stresses on different scales is provided.The paper is organized as follows: Section 2 describes the experimental treatment of the forming process of cylindrical specimen.The macroscopic FE simulation model is presented in Section 3. In Section 4, the macroscopic results of the experiment and the simulation are compared.A numerical process of the investigation of residual stresses on mesoscale and microscale is described in Section 5. Section 6 provides a conclusion.

Experimental Procedure
To investigate the development of residual stresses, the steel alloy 1.3505 (DIN designation according to [36]: 100Cr6) was thermo-mechanically processed by upsetting at 1000 • C with subsequent cooling.Following [20], cylindrical specimen with eccentric holes were used, as inhomogeneous residual stresses are expected over the circumference of the specimen.The specimen were prepared with six thermocouples as shown in Figure 2. Three thermocouples are located on the thin-walled side (positions are indicated as TC1 to TC3) and three thermocouples on the thick-walled side (TC4 to TC6), respectively.The prepared specimens were placed in a thermobox together with the tools as shown in Figure 3a which enables the upsetting test to be carried out isothermally.The temperature time profile is shown schematically in Figure 4, illustrating the process procedure.The prepared specimens were heated to 1000 • C inside the thermobox including the upper and lower tools in an oven, which has been manufactured by Nabertherm GmbH.After the desired temperature of 1000 • C had been reached, a ten-minute holding time followed to achieve a homogeneous austenitic structure.Since the austenitization conditions have an influence on the grain size and the grain size in turn has an impact on the phase-transformation behavior, all experiments were carried out using the identical austenitization conditions (1000 • C, 10 min), which result in a grain with the ASTM grain size number 8 according to [37,38].Subsequently, the thermobox including the tools and the specimen was transferred to the servohydraulic forming simulator.The specimens were upset from initial height of 50 mm to a final height of about 28 mm at a punch speed of 200 mm/s.In the next step, the specimen were removed from the thermobox and either quenched in cold water or cooled down to room temperature in still air.During the entire process chain, the temperature was measured via the thermocouples of type K (NiCr-Ni, manufacturer Günther GmbH).The bare thermocouple wires with a diameter of 0.5 mm have been welded directly onto the specimen.Using the measuring amplifier MX1609 (max.600 Hz, manufacturer HBM GmbH), the temperatures have been successfully recorded at a sampling rate of 100 Hz for quenching in water and 5 Hz for air cooling.
Due to the two cooling media of water and air, different temperature time profiles and consequently different transformation kinematics (see Figure 4) have been observed in the steel alloys.During rapid cooling in water, a diffusionless transformation occurs, where austenite is transformed into the martensitic body-centered tetragonal (bct) lattice structure.When cooled by air, a body-centered cubic (bcc) lattice with a ferritic, pearlitic, or bainitic microstructure depending on the cooling path is formed, as a consequence of a diffusion-controlled transformation.As a result, various residual stress conditions and distortions are expected in the differently cooled specimen.After thermo-mechanical processing, both axial and tangential residual stresses were measured on the surface of the specimen at the mid-height.The surface of the specimen at the locus of the measuring point was electrolytically polished, in order to avoid the influence of scale layers on the measurement.Residual stresses were determined by X-ray diffraction employing the sin 2 Ψ method in accordance with the corresponding DIN norm [39].For this purpose a Ψ-diffractometer type XSTRESS X3000 G2 from Stresstech GmbH using Cr Kff radiation was applied.The measuring spot was determined by a 2 mm point collimator, whereas the tilting area of Ψ was varied with a total of nine tilt positions.The evaluation of the residual stress measurements was carried out with the software XTronic (Stresstech GmbH).The X-ray elasticity constants (s 1 = −1.27• 10 −6 mm 2 /N, 1 2 s 2 = 5.81 • 10 −6 mm 2 /N) and the reference values for the unstressed material (d 211 = 1.1703Å, 2θ Cr = 156.084• ) were taken from the tabular data collection of [40].These data of X-ray elasticity constants were calculated as mean values from the so-called Voigt [41] and Reuß [42] approaches for the pure alpha iron lattice.
Furthermore, flow curves for the 1.3505 steel alloy were determined by means of standard cylindrical upsetting tests, in order to characterize the flow behavior in the hot forming process.The characterization tests were carried out on the servohydraulic forming simulator Gleeble 3800-GTC at IFUM for the relevant temperatures 900 • C, 1000 • C, 1100 • C and 1200 • C with the equivalent plastic strain rate 1 s −1 .In addition, the tests were carried out at a temperature of 1000 • C with the equivalent plastic strain rates 10 s −1 and 50 s −1 .Temperature control by means of thermocouples allowed the test to be carried out isothermally.

Numerical Procedure for Macroscopic FE Simulation
The previously described thermo-mechanical cylinder upsetting tests were numerically modeled with the commercial FE software Simufact.forming14 using MSC.marc solver.A time-temperature-transformation (ttt) diagram, accounting for polymorphic behavior of the investigated steel alloy 1.3505 was implemented in the FE-software; see [43].To simulate the viscoplastic material flow during hot forming, the flow stress of the material was implemented in the model as a function of equivalent plastic strain ε eq , equivalent plastic strain rate εeq and temperature T. To calibrate this multidimensional function, flow curves were experimentally determined as described in Section 3. The results of the characterization experiments as well as the calibration of the regression coefficients are presented in Section 4.
To numerically reproduce the thermo-mechanical forming process including the cooling step, phase-specific material data, such as flow curves, coefficients of TRIP, and specific heat capacity are required.In most cases, it is difficult to generate experimental material data for pure structural phases [44].For this reason, material data has been generated with the thermodynamic simulation software [45] as an initial estimation.With help of this software, it is possible to generate mechanical and thermal material data by means of thermodynamic calculation approaches from [46,47] based on the chemical composition of steel alloys.The calculated phase-dependent flow curves for austenite and martensite are exemplarily shown in Figure 5.In addition, more detailed information on the calculated phase-specific material data is presented in tabular form in the Appendix A. The respective chemical composition employed for the material data computation is shown in Table 1 and was taken from the inspection certificate of the investigated steel grade from the steel supplier.With this data in hand, simulations of the thermo-mechanical process are performed, in the manner described below.As mentioned before, it can be assumed that almost no residual stresses from previous process steps are present in the material due to austenitization.For this reason, it was sufficient to start the process simulation under the initial condition of a 100% austenitic structure and a thermally expanded workpiece at 1000 • C. Furthermore, the symmetry of the specimen was exploited by using a half model for the simulations.
Figure 6 shows the FE model in Simufact.formingand the considered boundary conditions.The tools were treated as heat-conducting rigid bodies, as it is usual in hot forming, and provided with a mesh of 1 mm edge length.For the tools made of 2.4668 (DIN NiCr19NbMo), constant values for specific heat capacity c p = 435 J/kgK and heat conductibility λ = 11.4W/mK were taken from the supplier's data [48].For the workpiece, an adaptive hex-mesh with a minimal element edge length of 0.5 mm has been employed, which was determined within a mesh sensitivity analysis.In total, 56,500 hex elements were used for the workpiece, whereby each tool was modeled with 24,000 hex elements.

Comparison of Experimental and Numerical Results for Macroscopic Model Validation
The determined flow curves at the investigated temperatures and equivalent plastic strain rates are shown in Figure 7a.The parameters (A, m 1 , m 2 , m 3 , m 4 ) of the following approach according to [49], cf.Equation ( 1) and Table 2, were determined using the method of the smallest error squares to continuously describe the flow stress k f as a function of equivalent plastic strain ε eq , equivalent plastic strain rate εeq and temperature T, eq e m 4 /ε eq εm 3 eq .
Table 2. Parameters determined for the calculation of the flow stress k f of the considered steel alloy 1.3505 according to [49].
As shown in Figure 7b, the plasticity model for the hot forming process step was validated by force-displacement-curves from experiment and simulation.The temperature time profiles are subject to several influencing parameters.Simufact.formingtherefore considers thermal data, such as specific heat capacity, heat conduction, and latent heat, which were calculated with JMatPro; see Appendix A. Since the modeling of the correct temperature profile is particularly important for the simulation of the phase transformations and finally the residual stresses, heat transfer coefficients (HTC) were also optimized.The HTC between steel and the cooling media air or water were determined in iterative numerical optimizations to accurately represent the cooling of the specimen in the simulation.Thus, they can be regarded as empirically fitted for the experiment under consideration.Initial values for temperature-dependent HTC for air and water were calculated based on literature data [50].In an iterative optimization procedure, the HTCs were adjusted using the experimental data from thermocouple 3 (TC).The calculated HTC values are presented in Table 3. Afterwards, the temperature profiles of the different measuring points from experiment and simulation, based on the calculated HTC data, were compared.As shown exemplarily for TC 2 in Figure 8, the temperature profiles could be successfully simulated.The upset specimen from experiment were optically investigated to evaluate the resulting distortions from thermo-mechanical processing.The CAD-model of a specimen from experiment could be created by 3D scanning, applying the optical measurement system ATOS (GOM GmbH).Subsequently, the dimensional accuracy of the FE simulation was visualized using the software INSPECT (GOM GmbH).Figure 9 shows the geometry of the specimen from experiment.The contour plot of the deviation e on this surface indicates the three-dimensional distance in perpendicular direction to the surface of the geometry from the simulation.A deviation of 0 mm shows a perfect match between the geometry from the simulation and the experiment.The algebraic sign of the deviation value indicates if the corresponding area of the geometry from the simulation lies outside (positive sign) or inside (negative sign) the workpiece shape from the experiment.As can be seen in Figure 9, a maximum deviation of about e = 0.4 mm was determined.
The hardness of the material depends on the resulting phase structure which is related to the transformation mechanisms and the cooling time.Hence, the hardness measurement according to Vickers [51] provides an important reference point for assessing the quality of the microstructure transformation simulation.For these metallographic investigations, the specimen were cut into two parts by wire erosion as shown in Figure 10a.Afterwards, Vickers hardness measurements were carried out at six measuring points (MP1 to MP6).The comparison of hardness values from experiment and simulation shows that the evolution of microstructure during cooling could be predicted well by the numerical model, c.f. Figure 10b.In addition, the corresponding microstructure in MP2 of the differently treated specimen is shown exemplarily in Figure 10c.As expected, a martensitic structure was formed by cooling in water, recognizable by the needle-shaped structure.A fine striped pearlitic structure was found on the specimen cooled by air, respectively.Consequently, the higher hardness values on the martensitic specimen after cooling in water and the lower hardness values on the pearlitic specimen after cooling in air are plausible.Here, the arithmetic mean with the absolute deviations of the tangential σ t as well as the axial σ ax residual stresses on the specimen surface at the locations of MP1 and MP6 on the mid-height are considered.In the X-ray measurements, negative stresses between 50 MPa and 100 MPa were detected on the specimen cooled by air.The simulation also predicts low residual stress values in the range of −5 MPa to −7 MPa for this cooling route.These minor residual stresses can be attributed to the diffusion-controlled phase transformation during cooling in air.For the specimen cooled by water, different values were measured on the surface of the thick-walled specimen side (MP1) and the thin-walled specimen side (MP6).The axial residual stress at MP1 is σ ax = 307 MPa and at MP6 σ ax = 100 MPa respectively.In the tangential direction, residual stresses of σ t = 183 MPa at MP1 and σ t = 138 MPa at MP6 were measured.As expected, the geometry of the specimen with an eccentric hole causes an inhomogeneous residual stress distribution.The developed FE model was able to predict the occurrence of tension or compression stresses for all the considered cases.In addition, the qualitative amount, especially the differences of the residual stresses between MP1 and MP6, could be predicted.However, the simulation overestimates the axial residual stresses at MP1 and MP6 as well as the tangential residual stress at MP1.The tangential residual stress at MP6 could be predicted with a good accuracy.
As shown in Figure 1, thermal, metallurgical and mechanical interactions in the material must be taken into account for the numerical description of the residual stresses.Due to these numerous interacting parameters, an exact calculation of the residual stresses is challenging.In particular, the TRIP strains take a significant impact on the resulting residual stresses [8].Therefore, the calculation of TRIP strains used in the simulation model by MSC.Marc solver could be extended.On the one hand, the load-dependent and reversible back-stress effects resulting from a stress relief during martensitic transformation can be included in the simulations.On the other hand, the work in [18] shows that the TRIP behavior has different effects on the deformation at tensile and compressive loads.For a further improvement of the model, these influences could be considered separately.With the simulation model presented in this paper, it was possible to make an appropriate qualitative prediction of the macroscopic residual stresses resulting from different cooling routes.
However, it is assumed that local stress gradients have also been developed, for example due to dislocations within the grains resulting from plastic deformations or phase transformations at the microstructural scale.These local deviations represented by residual stresses of second and third type can have a particular influence on the macroscopic stresses.Therefore, in future work in-depth numerical investigations on the microscale are to be performed in order to further analyze the interactions of the residual stresses of first, second, and third type and to quantitatively improve the prediction accuracy of the residual stresses; see Section 5.

Numerical Investigation of Residual Stresses on Micro-and Mesoscale
In the previous section, the thermo-mechanical process was numerically modeled on a macroscopic scale, taking into account the boundary conditions from the experiment.Since complex thermal, metallurgical, and mechanical interactions occur on the microstructural level of the material, a consideration of the residual stress development on these minor scales is motivated.By investigating the residual stress development on the microscale, further knowledge about the influence of the residual stresses of the second and third type on the measurable macroscopic residual stresses can be obtained, in order to further improve the prediction accuracy of the model.Therefore, based on data obtained from the experimental and simulative investigations of the previous sections, a numerical model is developed in the following, which accounts for the phase transformation from austenite to martensite (A-M) during rapid cooling and investigates the resulting residual stress distribution on a mesoscopic level.Herein, the mesoscopic level describes a polycrystal with multiple grains.Therefore, MPF simulations in combination with FE simulation are used and applied to the investigated material 1.3505.At this point, the complex process of phase transformation is simulated with a rather simple approach, which shall be understood as a first attempt towards a realistic model.
Since the A-M phase transformation changes the crystalline structure inside the grain but not the intergranular boundary, the latter is assumed to be fixed throughout the simulation.Starting from the microscopic scale on the level of one grain, the evolution of different martensite variants is simulated.The homogenized eigenstrain of the resulting laminate structures is then transferred to the mesoscopic level and thereby accounts for the strain jump during the evolution of martensite.A polycrystalline structure consisting of multiple grains is considered, which is obtained from MPF simulation applying normal grain growth.In this first attempt, the A-M transformation is considered by prescribing a homogenized eigenstrain to each austenitic grain.Using an energy minimization principle, the final eigenstrain distribution on the mesoscale is obtained from the MPF model.Since only elastic strains are considered, an FE model is used to incorporate the elasto-plastic behavior of the material.In a two-scale FE simulation associated with the experimental setting in Section 2, the eigenstrain distribution of the mesoscale is applied and the elasto-plastic response is computed.At the end, the resulting distribution of residual stresses is obtained.The following section elucidates the details on the briefly explained procedure above, which is additionally illustrated in Figure 12.

Multi-Phase-Field Simulation
The used MPF model, introduced in [25], can consider an arbitrary number of phases in the simulation, which differ from each other by e.g., orientation.An implementation is available as open-source software [31].The underlying governing equation is represented by the free energy density F in a domain Ω as where f intf represents the interfacial energy, f chem the chemical energy and f elast the elastic energy density, respectively.They can be calculated by where φ [•] denotes the phase-field parameter of phase-field [•].
The evolution of the phase-field parameter is controlled by the diffusion equation reading Microscale (FE-Simulation) Further quantities used in Equations ( 3) and ( 4) are described in Table 4.  2) and (3).i, j components α, β identifier of phase-fields N total number of phase-fields σ αβ energy of the interface between phase-fields α and β η αβ interface width between phase-fields α and β f α (c i α ) chemical potential for phase-field α c i α concentration of component i in phase-field α μi generalized chemical potential or diffusion potential as Lagrange multiplier for satisfying the balance of mass for each component (c i = ∑ α=1,....,N φ α c i α ) ε α total strain of phase-field α ε * α eigenstrain of phase-field α ε i α strain of component i in phase-field α C α elasticity tensor of phase-field α h(φ α ) polynomial function fulfilling h(1) = 1 and h(0) = 0 with h(φ α ) = φ α µ αβ effective interfacial mobility between phase-fields α and β

Homogenization of Eigenstrains on Microscale
To describe the transformation from austenite to martensite (A-M), different martensitic variants must be taken into account.This diffusionless transformation is caused by fast cooling of heated austenite in water as described in Section 2. In this context, the three so-called Bain-groups are considered (see [52]) which describe possible types of tetragonal transformation of the cubic unit cell, cf. Figure 13a.Shearing and rotation has to be considered additionally but will be neglected due to simplicity.In a first step, in order to simulate the martensite microstructure evolution inside one grain, austenite with three martensitic nuclei is considered in a MPF simulation, where each nucleus represents a possible martensite variant.A transversally isotropic eigenstrain state given by ε * I , ε * II or ε * III with diag ε * I = [−0.20.12 0.12] , diag ε * II = [ 0.12 − 0.2 0.12] , diag ε * III = [ 0.12 0.12 − 0.2] , represents the volume change during tetragonal transformation, I, II and III, respectively.Although these tensors do not fall into the range of small strains anymore, the linear approach is suitable due to high temperatures during application of the transformation scheme.This results in three possible laminate microstructures consisting of two variants, (i.e., I, II or II, III or I, III), following the principle of minimizing the total energy in the system; see Figure 13b.
To obtain the mechanical energy and the mechanical driving force, a homogenization scheme is applied to calculate the individual phase-field strain distribution, which is presented in [53][54][55] and applied to OpenPhase in [56].Therefore, a partial rank-one homogenization is used, which satisfies the static equilibrium and the strains compatibility condition at the interfaces.
The strain jump between phase-fields α and β for linear elastic materials in a small strain framework is given by cf. [54], where n is the normal vector of the interface and a characterizes the jump of the displacement gradient All other parameters are explained in Table 4.The effective (homogenized) eigenstrain tensor representing the internal laminate structure of one martensite grain can be computed by volume averaging, whereas the laminate consists of equal amounts of the two variants.This describes the transition between the microscopic single grain level and the mesoscopic polycrystal.Hence, the three effective eigenstrain tensors, which describe the eigenstrains homogenized over the domain of one martensitic grain on the mesoscale, read diag ε * I,II = [−0.04− 0.04 0.12] .

Mesoscopic Phase Transformation
For the description of phase transformation on the mesoscale, a polycrystalline structure with multiple grains is considered.At the beginning of the simulation, all grains are considered to be austenite, whereas the grains differ by the orientation of the crystallographic lattice.This granular structure is obtained from a MPF simulation involving normal grain growth (see Equation (3)1-2) where a cubic grid with 121 points in each direction is defined.It consists of one austenitic phase with 24 randomly placed nuclei with different orientations.The material properties for the considered material 1.3505 for nuclei and matrix are given by an isotropic elasticity tensor, which is derived by interpolation of Young's modulus E and Poisson's ratio ν generated using JMatPro, cf.Appendix A, for austenite and martensite at the martensitic start temperature taken as 375 • C, i.e., In the MPF simulation, the 24 grains start to grow from the nuclei and the matrix vanishes.The simulation is continued until the matrix has completely vanished.In order to account for the retained austenite, three grains of the evolved granular structure are considered to be austenite in the subsequent calculation of austenite-to-martensite (A-M) transformation.The amount of austenite in the mesostructure was observed in a conceptional study related to the experiments; see Section 3. Thus, the other 21 grains are considered for the simulation of martensitic phase transformation.Next, the deformation of the structure during compression must be accounted for.By assuming a linear mapping of the geometry of the cubic grid under the assumption that no volume change occurs, the bulging of the cylindrical specimen, resulting in a reduction in height by 50%, is described on the mesostructure.Subsequently, the linear mapping is performed with a factor of 0.5 for the reduction of height and the respective elongation in the two remaining directions considering volume preservation.The deformed grid is now defined by 61 grid points in height and 171 in width and length, respectively; see Figure 14b Hence, the structure represents a polycrystal consisting yet of austenitic grains.The A-M transformation is modeled by applying one homogenized eigenstrain state for each grain.Since it is not known which of the three possible homogenized eigenstrain states is dominant in each grain, a simple minimization scheme based on the total elastic energy f elast is used to identify a possible set of eigenstrain states.Therefore, the set of all eigenstrain states is defined by Initially, a random effective eigenstrain tensor ε * i is defined for each of the 21 grains.Next, it is checked whether one of the other two effective eigenstrain states in {ε * II,III , ε * I,III , ε * I,II } leads to a smaller energy level, whereas only one grain can be changed at a time.This iterative procedure is repeated several times until no change of eigenstrain state leads to a smaller energy level.The procedure can be described by the minimization scheme whereas Ṽ defines the set of eigenstrain states in the grains which leads to a minimization of the elastic energy.According to the application of this final set of eigenstrains Ṽ , a mesoscopic strain distribution ε MPF (x) results from the A-M phase transformation.Figure 14c shows the strain distribution on the mesoscopic structure, which serves as input, i.e., driving force, for the computation of residual stresses in the following section.

FE Simulation of Mesoscopic Residual Stress
For the mesoscopic modeling of the stress distribution because of the austenite-martensite transformation a two-scale FE simulation is considered.In particular, the calculation of the elasto-plastic mechanical behavior due to the evolving eigenstrains should be enabled.Based on that, the total strains ε MPF (x) resulting from the MPF model for the optimal set of applied eigenstrains (see Equation (10)) are taken into account as the driving force in the FE simulation.The mechanical behavior of each phase is represented by the application of an elasto-plastic material law with a von Mises yield criterion with exponential hardening; see e.g., [57].The set of necessary equations for the formulation is given in Table 5.
Table 5. Set of equations for von Mises plasticity with exponential hardening considering additive split of strain tensor ε into elastic ε e and plastic part ε p following [58]. kinematic The material parameters describing the mechanical behavior in the elastic domain, i.e., bulk modulus κ and Lamé constant µ, are recalculated from the values used in Section 5.1.2.The plastic behavior of the material model is fitted to the flow curves generated using JMatPro for both phases at 375 • C, cf. Figure 5a,b.The resulting material parameters are listed in Table 6.To incorporate the strain distribution ε MPF (x) resulting from the austenite-to-martensite (A-M) transformation evaluated using the MPF model into the FE simulation, the strain field on element level is incrementally modified over pseudo-time t, i.e., Therefore, the last and actual time step is indicated by t n and t n+1 , respectively.For the strain distribution ε MPF , which is computed for each grid point in the MPF model, a one on one mapping from grid point to FE is considered, i.e., ε MPF is constant over each element domain.To enable a simplified mapping, a voxel-like representation of the structure is achieved using hexahedral FEs with 20 nodes and quadratic shape functions.The same mapping is applied for the phase-field parameter to reconstruct the mesoscopic grain structure in the FE model.In order to reduce the computational costs, the grid of the MPF model from Section 5.1 is coarsed to a voxel-based representation of 35 × 35 × 13; see Figure 15.Based on the afore mentioned discretization of a granular mesoscopic structure, a two-scale boundary value problem is solved; see Figure 16.On the macroscale, a cylindrical section motivated by the specimen used in the experimental setup shown in Section 3, is considered.Therefore, the boundary conditions are enforced as shown in Figure 16 due to the assumption of an axis symmetry, where r, ϕ and z denote the radial and tangential and height component, respectively.On the mesoscale, periodic boundary conditions are applied.The classical micro-macro transition approach, also known as FE 2 -method, solves a mesoscopic boundary value problem at each macroscopic FE integration point.The macroscopic quantities are computed by suitable averages and the quantities associated with the macroscale are denoted by an overline.For further details see e.g., [59,60].In Figure 17a-c, the resulting mesoscopic stress distributions due to the application of the A-M transformation are shown.Since the equilibrium state on the macroscale is considered in the quasi-static simulation and no external forces are applied to the cylindrical specimen, the resultant stress distribution can be interpreted as residual stress, rather than mechanical stress.Stresses tangential to the outer surface of the cylindrical specimen (see Figure 17b) are in the main focus of further investigation.These quantities are most relevant with respect to the durability of the manufactured component.Additionally, it has to be noticed that the radial stresses vanish macroscopically, since σr = 0.875 • 10 −12 ≈ 0 holds, which is enforced by the boundary conditions, which state that σr = σz = 0 and σϕ = 0 as shown in Figure 16.As stated earlier, such stress states can be interpreted as residual stress due to the chosen boundary conditions.

Conclusion and Outlook
In this work, the material behavior of the steel alloy 1.3505 during the thermo-mechanical forming process was investigated by experiments and simulation approaches on different scales regarding microstructural transformation behavior and the resulting residual stresses.Upsetting tests of cylindrical specimen with eccentric holes were performed at 1000 • C with subsequent cooling in the media water or air to achieve different inhomogeneous residual stress states.Following this, the properties of the specimen were experimentally analyzed.Using an optical 3D scanning method of the experimental specimen, the dimensional accuracy of the geometries from the simulations could be evaluated.In addition, the resulting microstructure of the material was examined by light microscopy and the hardness was measured using the Vickers testing method.In addition, the residual stresses at certain points on the surface of the specimen were measured by X-ray diffraction.A macroscopic simulation model was created with the commercial software Simufact.formingusing experimental determined ttt-diagrams and flow curves.Furthermore, HTC were optimized for the considered process.The macroscopic FE model was able to qualitatively predict the residual stresses evolving from the considered thermo-mechanical forming process.
For the description of the evolution of mechanical fields due to microstructural effects on the lower scales, a multiscale approach was proposed based on a MPF model and a two-scale FE simulation.This method is used to describe the residual stresses on the microscale resulting from the austenite-to-martensite transformation.
In a future step, by analyzing the interactions of the residual stresses from micro-and macro-simulations, the prediction accuracy of the arising residual stresses in thermo-mechanical forming processes shall be further improved.From that, the influence of the residual stresses of second and third type on the measurable residual stresses of first type can be derived.Therefore, the improvement of the cooling process in the experiment and the validation of the macroscopic simulation to achieve a stable and useful residual stress state can be enabled.

Figure 2 .
Figure 2. Shape of the investigated specimen with dimensions in mm and location of the thermocouples on the specimen surface.

Figure 3 .
Figure 3. (a) Schematic representation of the specimen prepared with thermocouples in the thermobox and (b) servohydraulic forming simulator of the company Instron with the glowing thermobox inserted.

Figure 4 .
Figure 4. Principle representation of the temperature time profile in the experiments with respect to phase-transformation kinematics.

Figure 6 .
Figure 6.Schematics of the considered FE model with the positions of thermocouples (TC1 to TC6) and the boundary conditions.

Figure 7 .
Figure 7. (a) Determined flow curves for hot forming applications from servohydraulic forming simulator Gleeble 3800-GTC and (b) comparison of force-displacement-curves from experiments and simulations.

Figure 8 .
Figure 8. Time-temperature-profiles from experiments and optimized simulations from TC 2.

Figure 9 .
Figure 9. Visualization of the dimensional accuracy between the simulated geometry and the real specimen, analyzed with the optical measurement system ATOS (GOM GmbH).

Figure 10 .
Figure 10.(a) Positions of the measuring points (MP1 to MP6) along the horizontal axis S on the differently treated specimen, (b) results of Vickers hardness and (c) resulting microstructure after thermo-mechanical processing at MP2.

Figure 11
Figure 11 presents a comparison of simulation results with XRD residual stress measurements.Here, the arithmetic mean with the absolute deviations of the tangential σ t as well as the axial σ ax residual stresses on the specimen surface at the locations of MP1 and MP6 on the mid-height are considered.In the X-ray measurements, negative stresses between 50 MPa and 100 MPa were detected on the specimen cooled by air.The simulation also predicts low residual stress values in the range of −5 MPa to −7 MPa for this cooling route.These minor residual stresses can be attributed to the diffusion-controlled phase transformation during cooling in air.For the specimen cooled by water, different values were measured on the surface of the thick-walled specimen side (MP1) and the thin-walled specimen side (MP6).The axial residual stress at MP1 is σ ax = 307 MPa and at MP6 σ ax = 100 MPa respectively.In the tangential direction, residual stresses of σ t = 183 MPa at MP1 and σ t = 138 MPa at MP6 were measured.As expected, the geometry of the specimen with an eccentric hole causes an inhomogeneous residual stress distribution.

Figure 11 .
Figure 11.Arising residual stresses in (a) axial σ ax and (b) tangential σ t direction on the surface of differently treated specimen at MP1 (thick-walled side) and MP6 (thin-walled side).

Figure 12 .
Figure12.Simulation procedure of macro-, meso-, and microscale using Multi-Phase-Field Model and FE Simulation.Note that ε, ε p , ε MPF , ε * denotes the total strains, the plastic strains, the mesoscopic strain distribution resulting from the A-M phase transformation and the eigenstrain distribution per grain, respectively.

Figure 13 .
Figure 13.(a) Variants of martensite, cf.[52], with a 0 as edge length of an austenitic unit cell and a, b as edge length of a tetragonally transformed martensitic unit cell and (b) laminate consisting of two variants of martensite.

Figure 14 .
Figure 14.(a) Austenitic grain growth with different orientations for the grains, which are indicated by coloring, (b) austenitic cube with 50 % deformation regarding height and (c) first component of mesoscopic strain distribution from the MPF model ε MPF .11 .

Figure 15 .
Figure 15.Reconstruction of the periodic mesoscopic grain structure containing 21 martensitic grains and 3 austenitic grains: FE model with 15,925 cubic elements with 136,107 nodes.

Figure 16 .
Figure 16.Two-scale boundary value problem for the incorporation of the strains ε MPF (x) in the FE model.

Author Contributions:Funding:
Conceptualization, funding acquisition as well as supervision of this research project was done by B.-A.B., J.S. and D.B.; The experimental investigations as well as the numerical FE-analyses on the macroscale were carried out by C.K. with the support of A.C.; S.U.adopted the Multi-Phase-Field model based on work of M.S. and R.N. and performed the numerical simulations with help of L.S. and D.B.; C.K. and S.U.took the lead in writing the manuscript; All authors provided critical feedback and helped shape the research, analysis and manuscript.Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)-374871564 (BE 1691/223-1, BR 5278/3-1, SCHR 570/33-1) within the priority program SPP 2013.

Table 1 .
Chemical composition of the investigated steel alloy 1.3505 used for material data generation with JMatPro.

Table 3 .
Numerically identified heat transfer coefficients between specimen and the cooling media air or water for material 1.3505.

Table 4 .
Quantities governing the energy densities in Equations (

Table 6 .
Material parameters for austenite and martensite for 1.3505 at temperature T = 375 • C.