Parameter Identiﬁcation for Thermo-Mechanical Constitutive Modeling to Describe Process-Induced Residual Stresses and Phase Transformations in Low-Carbon Steels

Featured Application: Fatigue-resistant reinforcing bars in buildings and highway bridges where dynamic loading causes premature fatigue fractures. Abstract: Reinforcing steel bars (rebars) are widely manufactured using the Tempcore™ process. Several studies have been conducted analyzing the effect of the heat treatment route on the strength and corrosion resistance of rebars, but knowledge of its effects on the residual stresses of the ﬁnished product are largely lacking. This paper presents experimental investigations to identify the material parameters necessary to simulate the Tempcore™ process using thermo-elasto-plastic constitutive modeling in order to study the generation of residual stresses during the manufacturing process. Mechanical parameters such as yield strength at elevated temperatures and elastic constants were determined experimentally. A continuous cooling transformation diagram needed to model the phase transformations was also identiﬁed and is presented here. Residual stress distributions in the surface region of the rebar were characterized using X-ray diffraction. Further characterizations of microstructure, chemical composition, and hardness were carried out. The constitutive modeling approach for the numerical simulation is brieﬂy described for which the experimentally determined parameters are required as input.


Introduction
Tempcore™ is the trademark for the heat treatment route that was developed in the 1970s by CRM Group and has since been widely used for the manufacturing of reinforcing steel bars (rebars) [1,2]. It is a quench and self-tempering thermo-mechanical process wherein the hot rolled rebars come out of the last hot-rolling mill at a temperature above A c3 at which the rebars are still in their austenitic state and then are rapidly quenched with water sprayed onto their surface at high volumetric flow rate. The outer layer of the rebar transforms from austenite to martensite up to a certain depth. The time of quenching and the volume of water sprayed are controlled to achieve the desired thickness of the outer martensitic layer for a given diameter of the rebar, and below this depth the core is still austenitic. The rebar is then left to cool in the cooling bed, where the core transforms to a ferrite and pearlite mixture while the heat dissipating from the core tempers the martensitic layer, hence the name Tempcore. The process is schematically presented in Figure 1. This composite nature of the rebar with an outer hard layer and a ductile inner core is achieved with a relatively simple heat treatment process that produces high-strength rebars that are bendable and weldable without the need for micro-alloying with V or Nb [1][2][3]. The Tempcore™ process often has been investigated from the point of view of strength and corrosion resistance [4][5][6][7], fatigue behavior during service [8], behavior at high temperatures [9], and effects of their geometry [10]. Very few studies could be found in the literature that investigate the residual stresses that are generated during the manufacturing process itself. The research project from which this report is extracted aims at understanding the links between fatigue performance and residual stresses in the rebars, with the goal to enhance their fatigue life. It is well established that macroscopic compressive surface residual stresses enhance the fatigue performance of the material by closing fatigue cracks that usually initiate near the surface. On the other hand, tensile residual stresses lead to premature failure of parts by aiding crack growth and could also lead to stress corrosion cracking [11,12]. In service, the residual stresses are superimposed with the external load stresses and can be either beneficial or sometimes catastrophically detrimental by lowering the stresses at which cracks initiate to dangerously low levels [13].

Origin of Residual Stresses during Quenching
The residual stresses in the finished product are a consequence of the asynchronous thermal shrinkage and phase transformation of the core and the surface regions. During quenching all materials are exerted to shrinkage stresses since the surface always cools faster than the core. Additionally, if phase transformations occur during quenching, transformation stresses ensue as well and superimpose on the aforementioned shrinkage stresses. In order to understand the origin of the residual stresses during quenching and thereby the factors that affect it, a closer look at the generation of these two types of stress is necessary. The top part of Figure 2 shows the cooling curves of the surface and the core of a cylinder during quenching (T s , T c ) together with the yield stress, R y of these regions, which depends on the temperature and hence on process time and is assumed to be the same for compression and tension. For simplicity, only the longitudinal stress component is shown, although criteria for plastic flow also consider the tangential and radial components [14][15][16].
Concentrating only on the shrinkage stresses on the left side, the surface cools faster than the core at the initial stage of the quench. This results in longitudinal tensile stresses at the surface, whereas the core experiences longitudinal compressive stresses. As long as there is no plastic deformation, these stresses will be balanced at the end of the process with no residual stresses, as can be seen in the middle part of Figure 2a. Usually, however, these shrinkage stresses exceed the yield stress of the material shown as a gray shaded area in the figure. This is promoted by the low yield stress at high temperatures where the temperature gradients are very high [10,15]. Therefore, the material flows plastically leaving residual stresses at the end as shown in the bottom part of Figure 2a. Upon further cooling the core temperature falls rapidly while the surface gradually reaches the temperature of the cooling medium, thus reducing the magnitude of the shrinkage stresses and eventually reversing their signs. The plastic stretching of the surface causes compression at the surface and vice versa in the core. Once the temperature balance is established, the cylinder is left with compressive stresses at the surface and tensile stresses in the core due to the shrinkage stresses.  If the phase transformation in the material is to be considered, the same line of reasoning could be laid out for pure transformation stresses ignoring the thermal shrinkage stresses and transformation plasticity. The top part of Figure 2b shows the cooling curves of the surface and the core of the cylinder, the martensite start and finish temperatures, M s and M f , and the temperature-dependent yield stress of the material. Upon reaching M s it is assumed that the whole material transforms into martensite, which occurs at times t 1 at the surface and t 2 in the core. The volume expansion associated with the austenite to martensite transformation, which is ca. 3.4% for 0.19 % C steel [17], leads to compressive stresses arising near the surface compensated by tensile stresses in the core as shown in the middle part of Figure 2b. Once the core reaches M s and starts to transform and expand, the stresses are reduced in both regions. Due to the transformation into much harder martensite, the yield stress curves undergo a much larger increase in magnitude than in the case of mere thermal shrinkage as shown in Figure 2a, left. The plastic flow occurring due to the transformation stresses exceeding the yield point (shown in gray) results in volume incompatibilities between the core and the surface at the end of the cycle. This leads to residual stresses, which are shown at the bottom of Figure 2b by arrows, which are reversed in sign compared to the left image.
Combining these two types of stress by simple superposition is not accurate, since the material has deformed plastically. The fundamental concept from this discussion is, however, helpful in understanding the factors involved in solving such problems and the material parameters required for modeling such a quenching process successfully. It is clear that the yield stress of the material at high temperatures determines plastic flow and has to be quantified. The relative position of the time in which the transformation initiates and ends at the surface and in the core significantly influences the magnitude as well as the sign of the final residual stresses [14].
To the best of the authors' knowledge these parameters are not available in literature for the B500B alloy, which constitutes a large fraction of commercially available rebars and, hence, have to be determined experimentally. Some attempts have been made to predict the strength and transformation temperatures with models and theoretical calculations [5]. The thermal history or the cooling curve across the cross section of the bar during quenching and tempering is, however, impractical and prohibitively expensive to obtain directly from the manufacturing facility, if for this purpose the production of large tonnages has to be stopped. Therefore, several studies have been conducted to predict the transient temperature profile across the cross-section of the rebar with computational fluid dynamics simulations of the coolant coupled to the heat transfer problem in the rebar [5,18,19].
In summary, the following facts are important to consider. The plastic deformation determines the magnitude of the residual stresses. The temperature difference between the surface and the core is dependent on the following factors: the starting temperature of the quenching process, the diameter, the heat transfer coefficient, and the thermal properties of the material. The residual stresses are larger in magnitude if the yield strength at high temperature is very low, due to larger plastic deformation. Finally, the residual stresses greatly depend on the relative initiation time of the transformation in the surface and core regions.
The phenomena that were neglected in the above discussion for simplicity need to be considered for solving real problems. The material as a whole will not transform into martensite, but rather a mixture of all the phases including ferrite, pearlite, bainite, and martensite is present in the rebar as will be seen in Section 3.1. Only longitudinal stresses were considered, but the criteria for plastic flow usually consider all the components of stress [16]. Furthermore, transformation plasticity effects that were neglected may have a significant influence on the residual stress distribution and should be examined [20]. In addition to the material parameters, the geometry of the ribs on the rebar and its surface finish also significantly influence the fatigue behavior [8]. Experimental determination of residual stresses, e.g., using X-ray diffraction, will reach their limitations for such complex surface geometry. Furthermore, the resolution is not high enough to measure the residual stresses at the foot of the ribs, which usually act as the origin of cracks in these bars. Hence, it is only through simulating the quenching process with the correct material parameters, material constitutive model, and rebar geometry that a better understanding of the process can be gained.

Constitutive Modeling for Numerical Investigations of the Tempcore™ Process
Since it is not feasible to inspect the transient stresses shown in Figure 2 by in-process measurements, numerical simulations play an important role in understanding the generation of residual stresses during the heat treatment and thereby potentially attain further optimization of the process parameters. Given the requirement to resolve the macro-geometry of the considered component within the spatial discretization, constitutive modeling for heat treatment simulations is often based on component-scale continuum formulations where the microstructure is accounted for only be means of volume proportions of the existing phases [5,6,19,20]. In the simplest scenario, this approach typically involves solving the following initial boundary value problems (IBVPs): • The transient non-homogeneous isotropic heat equation governing heat conduction in the spatial domain Ω ∈ R n spanned by the solid and the temporal domain [t 0 , t e ] ⊂ R with prescribed surface convection on ∂Ω and a known initial temperature field T 0 (x) [21]: Here, ρ, c p , and λ denote the density, specific heat capacity, and thermal conductivity of the compound, respectively, T ∞ the coolant temperature, h the heat transfer coefficient, andq the internal heat generation, e.g., due to latent heat associated with phase transformations or plastic dissipation. • The mechanical compatibility and (quasi-)static equilibrium equations [21] where geometric linearity and negligence of body forces are typically suitable assumptions for heat treatment simulations [20]: In this system, σ denotes the Cauchy stress tensor, the linear strain tensor, and u the displacement field. Equations (4) and (5) are usually complemented with displacement-or traction-free boundary conditions on a partition ∂Ω u∪ ∂Ω σ = ∂Ω of the solid's boundary (the outward normal of which is denoted as n) as well as stressand displacement-free initial conditions: To complete the description of the mechanical IBVP, a constitutive model that particularizes Equation (5) must be specified, where the assumption of isotropic thermoelasto-plastic behavior according to the classical flow theory of plasticity is almost universal in literature [20]. This amounts to computing the stresses σ from the elastic strains e as where th is the eigenstrain due to thermal expansion (and potentially phase transformations) and p is the plastic strain computed by integration of the (associated) flow rule˙ subject to the constraints γ ≥ 0, f (σ) = 3 2 σ : σ − σ Y ≤ 0 (von Mises criterion) and γ f (σ) = 0 [16]. In these equations, I and I denote the identity tensors of rank 2 and 4, respectively, σ the stress deviator, γ the plastic multiplier, σ Y the yield strength, and E and ν Young's modulus and the Poisson's ratio, respectively.
In order to apply the above IBVPs to the numerical simulation of heat treatment of a solid compound of more than one phase, the material parameters encountered in above equations are most commonly computed via a linear volume-fraction averaged mixture rule from the respective parameters of the individual phases [20]. For example, the yield strength σ Y of a mixture of P phases, each constituting a volume fraction z i = V i V −1 , i = 1...P, of the microstructure at the considered point x ∈ Ω of the continuum, may be assumed to adhere to where σ y i denotes the yield strength of phase i. Furthermore, evolution equations describing the phase transformation kinetics in terms of the volume fractions z i of each phase are required to evaluate such mixture rules.
It is apparent from the above presentation that any numerical computation of processinduced residual stresses based on a thermo-mechanical constitutive model of comparable complexity as above relies on the accurate identification of at least the following parameters for each phase encountered during the process: density, heat capacity, and conductivity for the heat transfer problem and Young's modulus, Poisson's ratio, and the yield strength for the mechanical problem. The aim of this paper is therefore to enable numerical heat-treatment simulations for rebars by experimental determination of the mechanical parameters for the alloy B500B.

Materials and Methods
Rebars are made from low-carbon steels with very little amounts of other alloying elements. The specimens for the studies were extracted from B500B rebar samples [22] with a nominal diameter of 28 mm obtained from industry partners. The chemical composition of the samples was measured by spectrometry on the cross-sections of the rebars. The carbon composition of the alloy varied between 0.17% and 0.23% between various samples. The mean chemical content from four rebar samples analyzed and the mean is listed in the last line of Table 1. Trace amounts of other elements less than 0.01% are not displayed.

Tensile Tests at Elevated Temperatures
The yield stress of individual phases are essential in modeling the residual stresses as described in Section 1.1. The tensile tests were carried out in a quenching and deformation dilatometer, Dil805 A/D/T from Baehr Thermoanalyse (TA Instruments). The flat specimens used for the tensile tests had a geometry as shown in Figure 3. The temperatures at which the phases were tested were chosen according to the temperatures the individual phases experience during the Tempcore™ process. As shown in Figure 1, the martensite that is generated during the initial quenching step sustains from 200°C up to 600°C in the following tempering step. The austenite in the core experiences slow cooling and would start transforming at A c3 ; hence, it was tested in the temperature regime above and slightly below A c3 [2,19,23].
The reason for using the dilatometer instead of a conventional tensile testing machine is the fast inductive heating in the dilatometer, which allows the specimens to be heated quickly to the testing temperature. This is important because the specimens have to be tested at some temperatures where the phases are unstable and will start transforming. In the conventional tensile testing machine with an oven, the specimen has to be held at the testing temperature for approx. 30 min to reach a homogeneous temperature distribution. This would lead to substantial tempering of the martensite [24]. In the dilatometer the specimens were inductively heated at ca. 10 K s −1 to the testing temperature, and the strain rate was 0.02 s −1 . This allowed the specimens to be tested within two minutes. In order to obtain the stress-strain curve of austenite, a thermocouple was tack welded in the middle of the specimens, where the temperature reaches its maximum. After mounting the specimens in the dilatometer, they were heated above A c3 to 850°C for 60 s to ensure homogeneous austenitization. Thereafter they were quickly cooled with a cooling rate of ca. 100 K s −1 to the testing temperature and tested.
To investigate the martensitic phase, specimens were first austenitized by soaking in a heat treatment oven at 850°C for 30 min and afterwards water-quenched, since the quenching rate in the dilatometer was insufficient for in situ transformation to a fully martensitic state. The Vickers hardness of these quenched specimens exceeded 500 HV30, which corresponds to a fully martensitic structure [25]. Thereafter, isothermal tensile tests were conducted at various temperatures reached with a heating rate of 10 K s −1 . The tests lasted ca. 2 min, including heating the specimen from room temperature to the testing temperature. In practice, the tempering time is directly proportional to the mass of the workpiece [26]. Since the specimens are thin sheets of metal, they are likely to be tempered to some extent even within the duration of testing [27,28]. Tempering, however, also occurs during the Tempcore™ process when the martensite in the surface region is heated up; therefore, the results of the experiments are still relevant to the process.

Elastic Constants
Determining the Young's modulus E using the stress-strain curve is erroneous given the maximum length of the tensile specimen in the dilatometer is only 10 mm. Impulse Excitation Technique (IET) is a precise and reliable measuring method for the determination of elastic properties. The sound waves emitted from a free vibrating specimen are recorded by a sensitive microphone, allowing to evaluate the fundamental resonant frequencies [29]. In case of metallic alloys the resonance frequencies depend mainly on the geometry of the specimen, the elastic properties, and the density of the material. By analysis of simple geometries, such as cylinders or rectangular bars with respect to the fundamental flexural eigenfrequency, it is possible to calculate the Young's modulus in case of homogeneous and isotropic elastic properties by analytical approaches provided by the ASTM-standard [30]. The simplicity of the experimental set-up allows to perform the measurements at room temperature as well as at elevated temperatures.
The specimens were rectangular prisms with dimensions of 60 × 15 × 4 mm 3 . A standard specimen with the ferritic-pearlitic microstructure and a hardened specimen with martensitic microstructure were measured using a single specimen. To achieve this, the measurements were done while heating a hardened specimen from room temperature to 900°C and also measured while slowly cooling it down after 5 min of soaking to let it austenitize completely. For hardening the specimen, it was austenitized in a heat treatment oven at 860°C for 30 min and then quenched in water. Then Young's modulus was measured while heating from room temperature at a rate of 16 K min −1 and while slowly cooling down at 5 K min −1 .
Following the ASTM-standard, the Young's modulus is calculated by where l is the length, b the width, and t the thickness of the specimen. The specimen mass is denoted by m and the measured fundamental flexural resonance frequency by f f . The formula is based on Euler-Bernoulli beam theory, whereby the correction factor T 1 accounts for the finite thickness of the specimen. A sensitivity analysis of Equation (11) shows that the calculated Young's modulus is highly sensitive to the specimen thickness. It is challenging to evaluate the thickness exactly due to the distortion caused during the quenching. Hence, the specimen was subsequently ground with sand paper by hand in order to remove the slight curvature. Finally, the dimensions of the specimen were measured by use of a high precision micrometer screw at six different positions for determining its thickness and width and at three positions for the determination of its length.

Residual Stresses
Residual stress measurement techniques fall mainly into either of the two categories of destructive or non-destructive methods. Destructive techniques include sectioning, hole drilling, nanoindentation, cut compliance methods, etc. [31][32][33][34]. The fundamental concept behind these techniques is to mechanically remove a section of the workpiece by slitting or sectioning and measure the resulting deformation, usually using strain gauges, and then reconstruct the initial residual stress field. Non-destructive methods include diffraction using X-ray or thermal neutrons from reactors. The principle is to use the lattice structure of crystalline materials as an intrinsic strain gauge. X-rays generated in the lab using standard X-ray sources have an attenuation length of a few micrometers, whereas X-rays from synchrotron sources and neutrons can penetrate the material up to a 7-8 mm [35,36].
The near-surface residual stresses in this work were measured by X-ray diffraction (XRD) using the device Xstress3000 with a G2R goniometer from the manufacturer stresstech (Vaajakoski, Finland). The experimental parameters for the stress measurements with XRD are shown in Table 2. The residual stresses were measured in between the non-uniform ribs of the rebars using the sin 2 ψ method [36,37]. The penetration depth of the Cr K α radiation was calculated to be 11 µm in this alloy using the software xraylib [38]. Thus, the effective depth of measurement is about half of that, i.e., 5 µm when the radiation is directed perpendicular to the surface. The in-depth stress profile was determined by removing thin layers of material using electrochemical etching and measuring the stresses on the new surface using XRD. The depth increment was measured using a high-precision dial gauge with a resolution of 1 µm. Due to the roughness on the etched surface and the need to remove the specimen between etching and measuring steps, the accuracy reduced so that the realistic accuracy was ca. 10 µm. The stress relaxation due to the removal of layers was corrected by using the stresses measured before the removal and the amount of material removed. The foundations for these corrections were laid by Moore and Evans [39] and are integrated into the software of the Xstress device.

Micro-Hardness and Metallography
The micro-hardness profile along the radius of the rebar shown in Figure 4 gives a good picture of the effect of the heat treatment process that the material has undergone. It is further useful to check if the product meets the tolerances specified and to compare the specimens with specific microstructures produced by heat treating in the lab [40]. The micro-hardness was measured using a LECO Micro-Hardness tester with a load of 500 g and a dwell time of 15 s. The microstructure was investigated in the transverse as well as in the longitudinal cross-section of the rebar, which after polishing were etched with Nital for optical microscopy. Figure 5 shows the microstructures observed at different positions on the transverse cross-section of the rebar. On the top of Figure 5 where the positions of the individual microstructures are marked, a clear demarcation can be seen about 3 mm beneath the surface of the rebar. The core shows a mixture of ferrite and pearlite whereas near the surface tempered martensite can be seen inside prior austenite grains. In the transition zone where the cooling rate was not fast enough for martensite formation, a bainitic microstructure is observed. The microstructures obtained from the longitudinal cross-section were similar to those of the transverse cross-section and are not shown here. The longitudinal cross-section was further used to investigate any difference in microstructure on either sides of the transverse ribs that might have arisen due to the water being sprayed from one direction only inside the quenching chamber of the manufacturing unit. No differences were observed in the grain sizes, hardness, as well as in the thickness of the tempered martensite layer.

Tensile Tests at Elevated Temperatures
The stress-strain curves as a function of temperature are shown in Figure 6. The martensite specimens were deformed up to fracture (Figure 6a), whereas for the austenite samples the elongation at fracture could not be determined since it is beyond the limit of admissible elongations of the displacement sensor of the dilatometer. Testing austenite at 700°C did not result in reliable data, since it was not possible to avoid its transformation into ferrite and pearlite. A higher strain rate of 0.05 s −1 during the tensile test was also tried to shorten the testing time, but it had no effect as the recorded stress-strain curve was very similar to the one shown in Figure 6b; therefore, the decomposition of austenite could not be avoided since it occurs within 5 s. The yield stress R p0.2 and the ultimate tensile strength R m from each curve is listed in Table 3. The temperature-dependent Young's moduli determined using the resonance frequency measurements, shown in Figure 7, were utilized in order to deduce the R p0.2 from the stress-strain curves.
The yield stress of austenite at temperatures above A c3 ranges from 37 to 77 MPa. This is very low compared to the shrinkage stresses that are generated during quenching. Referring back to the discussion in Section 1.1 and Figure 2a, it can be seen that the austenite in the surface region and in the core would undergo a large amount of plastic deformation during quenching.
The hardness of the martensitic specimens after the tensile test was around 200 HV30, which corresponds to tempered martensite. The same values were also observed in the near-surface tempered martensite region of the rebar in Figure 4. Thus, it was not possible to avoid the tempering of the martensite during the tensile test, although the duration of the test was kept as short as possible. At high temperatures, the migration of carbon atoms started within seconds, followed by precipitation of transient iron carbides in the first few minutes [28]. Consequently, the martensitic specimen gained substantial ductility at the higher temperatures, as seen in Figure 6a.
The tensile test at room temperature was conducted in the universal tensile testing machine on a standard cylindrical specimen. The extensometer was removed after reaching R p0.2 , and the strain until fracture was evaluated on the basis of the crosshead displacement. The specimen tested at 300°C in the dilatometer also had a slightly smaller cross-section and shorter length than the other flat specimens. Due to these different specimen geometries, the elongation values from these tests are not directly comparable to the other tests when strain localizes beyond uniform elongation.

Elastic Constants
The Young's modulus, bulk modulus, and Poisson's ratio were determined for the base material (ferritic+pearlitic microstructure) at room temperature and were 206.6 GPa, 81.26 GPa, and 0.273, respectively, which corresponds to literature data [41]. The measured Young's moduli over the temperature range of 25-900°C for the specimen with initial martensitic and ferritic-pearlitic microstructures are shown in Figure 7. Both the specimens undergo phase transformation between A c1 and A c3 , which can be seen as a pronounced decrease in Young's modulus above about 700°C. Young's modulus of austenite remains almost constant up to 900°C. The furnace was purged with argon gas while heating up to reduce the oxidation of the samples at high temperatures. Practically the oxidation cannot be completely prevented, which leads to changes in the specimen dimensions and their mass during the experiment, which causes uncertainties in the calculated Young's modulus. Thickness and mass were also measured after the experiment to determine the thickness of the oxide layer. Overall, this effect can lead to an error of up to ±4 GPa in E. The same experiment was performed for a specimen with initial ferritic-pearlitic microstructure in order to check for differences resulting from the two variants of the initial microstructure. It turned out that these differences were in range of the estimated uncertainty.

Continuous Cooling Transformation Diagram
The continuous cooling transformation (CCT) diagram for the particular alloy is indispensable for studying its transformation characteristics. The diagram shown in Figure 8 is taken from literature [42] and published here with the authors' permission. The chemical composition of the alloy of the diagram is compared to the rebar composition in Table 4. It should be noted that both the carbon content of the tested rebars (0.17% to 0.23% wt% (cf. Table 1)) and the content of Mn and Si matched the composition of the CCT-steel quite well. Hence, it is reasonable to assure that depicted CCT diagram is relevant for the rebar steel studied in this work. Furthermore, three heating and cooling schedules for the rebar steel were realized with different cooling rates and were compared with the CCT diagram as shown in Table 5. The hardness of the specimens from this heat treatment corresponding to the ones in the CCT diagram are shown in Figure 8 as a, b, and c. During heating, the temperatures A c1 and A c3 agreed with those noted in the diagram 715°C and 835°C, respectively. During cooling, the starting and finishing temperatures of the transformations, T 1 and T 2 , and the hardness values at the end of the transformations showed good agreement with the values given in the CCT diagram.

Residual Stresses
The near-surface residual stresses in the rebar were measured at the location between two transverse non-parallel ribs and are shown in Figure 9. Both the longitudinal (axial) and tangential stresses are compressive, and they vary in magnitude almost in a linear fashion from 10 MPa on the surface to 90 MPa up to a depth of 100 µm. Beyond this depth the magnitude remains the same for up to a depth of about 300 µm. Compressive stresses of about 50 MPa near the surface of the rebar were also observed by Rocha et al. [10], but were documented only to a depth of 50 µm. The stresses shown in Figure 9 are compressive and much lower than the tensile yield strength of tempered martensite at room temperature, which is the predominant microstructure constituent near the surface of the rebar. No stress reversal of compressive to tensile stresses was observed in the studied depth range. The measurement uncertainty depicted is evaluated in the course of the linear regression analysis of the lattice spacing to the sin 2 ψ values according to the sin 2 ψ-method [36,43].

Summary
Conclusive studies on the generation and distribution of residual stresses in rebars manufactured through the Tempcore™ process are still lacking in literature. To help understand the phenomena contributing to the residual stress generation by means of numerical simulations, the parameters of the underlying thermo-mechanical constitutive model have to be determined. In this contribution, the parameters needed for the rebar alloy B500B were characterized as follows: 1.
The source of residual stresses was clarified, and the necessary material parameters for modeling stress distribution were identified.

2.
Metallographic analysis and micro-hardness measurements of the rebar across the cross-section were performed. Tensile tests were conducted to characterize the mechanical parameters of individual phases of the rebar, viz. austenite and (tempered) martensite, as a function of temperature.

3.
The elastic constants were determined using the resonance frequency method to obtain more accurate values than from the tensile tests. The results were also presented as a function of temperature. It was seen that there is a negligible difference in Young's modulus between different phases below A c1 . 4.
The CCT diagram was identified and presented for the specific alloy of the rebar that is necessary for further investigations in which phase transformation kinetics are to be considered. The relevance of the diagram matched well with the chemical composition of the alloy as well as results from dilatometric investigations.

5.
Residual stress distribution at the surface and in the subsurface region of the rebar were characterized using X-ray diffraction. The stresses were compressive with a magnitude of about 40 MPa near the surface and about 100 MPa beyond a depth of 100 µm. Stresses did not reverse their sign in the studied range of 300 µm below the surface.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.