Prediction of Preheating Temperatures for S690QL High Strength Steel Using FEM-Simulation for High Power Laser Welding

This study investigates a method for predicting the effect of preheating temperatures on the resulting hardness for high power laser welding of high strength steel. An FEM model is introduced containing a hardness calculation based on an existing model. Moreover, the hardness values of experimental results have been measured in order to show the performance of the model. The hardness calculation requires the chemical composition and the t8/5-time at the point of measurement. It is claimed that a calibration of the melt pool width and depth at room temperature only is enough to get reasonable results from the FEM-model for higher preheating temperatures. From the experimental result of a single experiment the width of a weld seam and the depth was deducted. In this study experiments have been done at various preheating temperatures in order to show the correlation between the model and the experimental results at various temperatures. The hardness equation provides suitable results in the verification with the measurements. The prediction of preheating temperature can be done with the resulting t8/5-time of the FEM-model. This method can decrease the amount of time and costs within a production according to testing and analyzing a matrix of process parameters. Moreover it is concluded that this methodology might be used for single item production.


Introduction
The use of high-strength steel (HSS) gets more important today. With the increasing mechanical properties, the quantity of material and the weight of constructional parts can decrease [1]. Due to the reduction of the overall weight, economic and environmental advantages arise in industries like e.g. the automotive industry. In industrial production process chains the takt time determines the amount of time available for specific welding tasks. Increasing productivity and decreasing takt times result in the demand of higher process speed. More often laser deep penetration welding processes are used in industrial applications in order to achieve elevated processing velocities. In addition, it is possible to generate deep welding joints with single passes in comparable small production times [2,3].
However, the laser welding process comes with a high-power density inside the laser spot and comparably small diameters. From this follows a strong temperature gradient with rapid cooling speeds inside the heat affected zone (HAZ). Martensitic phase transformation in the weld seam and the HAZ can be the consequence depending on the material. The original micro structure inside the material will be lost in this region. Moreover, welding distortion and residual stresses are largely affected by phase transformation in the cooling process [4]. The amount of martensite and size of the HAZ are a function of the temperature field and cooling rates respectively. In general the type of welding process and the welding parameters can indirectly control the temperature field and cooling rates. Moreover, the temperature distribution is following the geometry of the joined parts itself. Hence, the unwanted effects are a function of the process conditions. For instance cracks can occur at different process parameters [5] and the formation of pores can be reduced by the process parameters as well [6].
The hardness in laser beam welded high strength steel can be comparably high. Often a post treatment of weld seams is necessary to restore the mechanical properties of the material. This increases the production time and cost of laser beam welded constructions [7]. However, the hardness can be controlled by adapted temperature fields.
Some strategies like preheating can help to reduce the cooling speeds and hence decrease the drawbacks of laser deep penetration welding [8]. An investigation showed that the effect of increasing hardness at the welding seam and HAZ can be reduced by preheating. Correlations between hardness and preheating temperature have been reported in several studies by empirical analysis of measurements and process parameters [9,10]. Nevertheless, finding the right preheating temperature can be a cost intensive experiment. In that method the estimation of hardness of weld seam as a function of the preheating temperature is reported. It was possible to conduct an empirical equation for the calculation of the hardness. However, experiments at different preheating temperatures had to be done.
Another investigation shows the prediction of hardness using a neural network model [11]. However, the amount of requested data is high and the portability to transfer this method to other materials is doubtful. Even though the model is not suitable for the preheating temperature prediction. An investigation showed a method to calculate the hardness in the HAZ with the chemical composition and the cooling time between 800°C and 500°C (t8/5-time) [12]. Measuring the temperature within the weld seam and HAZ in practical conditions is often not possible. Even though the measurements would be affected by the measurement method itself. Therefore, the prediction method is less useful for practical tasks.
Using the finite element method (FEM) can provide a realistic temperature field within three dimensional parts. The increase of requested values also increases the calculation time. Even though the results are heavily affected by the quality of the requested process parameters. Using default setting for complex geometries or process parameters can provide insufficient results. Welding simulations can be subdivided into three parts [13]. A Process simulation can provide weld pool profile and dynamics according to process parameters. The heat source is defined by the interaction of laser beam profile and material parameters [14]. Structure and material simulations can provide insights into residual stresses and distortions within three-dimensional part. The heat source is parametrized previously and remains static during the simulation [15,16].
The target of the following investigation is to create a thermos-mechanical calculation model that provides preheating temperatures according to a defined hardness in the weld seam. Therefore an empirical hardness calculation and a FEM model are used. The model is to be validated by experiments. The resulting methodology can be used in any industry working with high strength steel in order to calculate the correct preheating temperature, e.g. the crane industry.

Methodology of Validation
The basic idea is to use modern FEM-simulation tools for the prediction of properties as a function of preheating temperatures. In practice the validation of the temperature field in the process can be an exhausting process. Hence, in this study it is claimed that a calibration of the melt pool width and depth at room temperature only is enough in order to get reasonable results from the model determining the properties at higher preheating temperatures. From the experimental result of a single experiment the width of a weld seam and the depth was deducted. This was used to model the heat source in the FEM calculation. The methodology of this study is divided into two parts. First the actual method for the prediction of preheat temperatures is made. The second part of the investigation is the validation of the method. Therefore, the results of the method's calculation are compared with the experimental measurements. This validation is not part of the actual prediction method and has only been done in order to show the performance of the model. For the prediction method a welding sample at room temperature (20°C) is made. A metallographic analysis provides the cross-sectional dimensions. Therefore, the depth and wide at the top and bottom of the weld seam are measured. The heat source within the FEM model is calibrated according to this measurement. Furthermore the welding velocity and laser power are integrated into the model. The FEM calculation is done for various preheating temperatures. The resulting t8/5 times for defined elements within the weld seam are exported to the hardness calculation model. Also the chemical composition for the welding sample is provided by an energy dispersive X-ray analysis. The hardness can be calculated according to a specific preheating temperature.
For the validation of the prediction method welding samples are made at various preheating temperatures. The hardness is measured in the weld seam. These measurements are compared with the calculated values for specific preheating temperatures.
The following flow chart describes the method of preheating temperature prediction and the validation of the method itself. See Fig. 1.

Material
Metal sheets made of S690QL were used for the experiments. The metal sheet thickness is 6 mm. The plates have a length of 200 mm and a width of 50 mm. The chemical composition has been generated for the base material and the weld seam by an energy dispersive X-ray analysis (EDXA). This analysis was made with an acceleration voltage of 20 kV. The results are shown in Table 1. Due to the insufficient detectability of carbon the value was set to 0.16 wt% according to [17]. This relates to the material composition in the FEM-model. The absolute error was set to + − 0.04 wt%. The differences between the chemical composition of the base material and the weld seam are neglectable. According to this the following calculations are based on the chemical composition of the base material. The two plates were clamped in the lap joined configuration. The length of the weld seam is 150 mm and 25 mm from each side. The laser source is an IPG YLS-10000. This laser is a fiber laser with a maximum output power of 10 kW. The fiber diameter is 200 μm. A collimator lens of 160 mm and focusing lens of 300 mm were used. The focus was set 4 mm below the surface of the upper metal sheet. The used laser power was 6 kW and a travel speed of 1 m/min. The experimental setup can be seen in Fig. 2.
To investigate the accuracy of the validation the results of the model were compared to the measurements of the experimental welds. An accordance of the validated and measured values is the proof for the model.
For this the welding experiments have been performed with and without preheating. For the preheating the specimen where put into an oven. After the heating process the specimens have been placed into the welding setup. The temperature was measured with a tactile thermometer. The welding process was performed after the plates reached the desired temperature. Experiments have been made at room temperature and at 100°C, 150°C and 200°C respectively.
The specimens have been investigated by metallographic inspection. Hardness measurements were made in the weld seams at 3 mm below the surface. The Vickers hardness was measured with a 10 kg load.

FEM-Model
The FEM welding model has been set up with a software named Simufact welding. The model is shown in Fig. 4. It has already been reported about Simufact in industrial applications, [18]. The temperature dependent thermophysical and mechanical properties of matter have been provided by the software [17,19,20]. The mesh geometry of  The geometry of the heat source was adapted by the laser beam parameters and the cross-sectional area. Width and depth of the weld pool have been taken from the experiment without preheating. Also the process parameters of the experiment are taken into account. See Table 2. The configuration of the heat source is calibrated with the measured depth and width of the weld seam cross section for room temperature. See Table 3. The efficiency was set to a typical value of overall absorption for deep penetration welding processes [21].

Hardness Calculation
The model describes the calculation of maximum hardness in HAZ for weld seams. Adapting this method is used for hardness calculations in weld seams. According to this method the hardness is connected to the t8/5 time by means of the arctangent function. It is necessary to calculate the maximum hardness and t8/5-time of a full martensitic and full bainitic structure. These are based on specific carbon equivalents. The final equation is described as followed, Eqs. 1 and 2. Specific calculations according to [12].
With this equation it is possible to calculate the maximum hardness of a specific t8/5time. These are given by the FEM simulation in this study.

Metallographic Analysis
All weld seams show a sound bonding without any kind of defects like cracks or pores in the metallographic cross sectional investigation. The width of the weld seam at the surface is 3.6 mm, 1.6 mm at the bottom and the depth is 10 mm. In Fig. 5 an exemplary cross section without preheating can be seen. The measured width and depth are used for the calibration of the FEM-simulation. However the depth can be calculated by analytical approaches [22] but it was decided to go for a practical solution and measure it in order to be able to transfer the model to more complex geometries. See Table 3. The geometry of the weld seam geometry and the calculated seam geometry show a sound accordance. The microstructure shows a martensitic structure in the weld seam.

Hardness Measurement
The Vickers hardness of all specimen haven been measured with a 10 kg load. The results are shown in Table 4. The relative error was set to 1.0% of measured value. The maximum measured Vickers hardness occurs at the weld seam without preheating with 407 HV10. The higher the preheating temperature the lower the hardness in the weld seam. At a preheating temperature of 200°C the measured Vickers hardness drops down to 325 HV10. The initial measured Vickers hardness of the specimen is 285 HV10.

FEM-Simulations-Results
The requested t8/5 time is provided by the temperature of specific nodes in the simulation. In Fig. 6 the result from the FEM simulation without preheating is shown. Each node of the FEM-model contains a separate temperature for each time step. The following data contains to the node in the middle of the weld seam at 3 mm below the surface. The results for each preheating simulation are shown in Fig. 7. According to the temperature-time dataset the t8/5 times have been calculated. The accurate times of 800°C and 500°C are calculated with linear interpolation. With the increase of preheating temperature, the t8/5 times also increase. Due to the chosen time increments the total error for the t8/5-time is set to 0.2 s. See Table 5.

Hardness Calculation
The calculation uses the chemical composition of the base material and the t8/5-time from the FEM model. In Table 6 all calculated values can be seen. The hardness decreases with increasing preheating temperature. The error propagation has been used according to Eq. 3 and [23] in order to calculate the error in the hardness.
The error of the calculated value f is given by Δf. The different variables and individual errors are given by x n and Δx n respectively. All values have been integrated in this calculation. The calculated hardness at a preheating temperature of 20°C is 396 HV10. This is the maximum hardness that has been calculated. With the increasing preheating temperature, the calculated hardness decreases. At 200°C preheating temperature the calculated hardness drops down to 342 HV10. This can be seen in Fig. 8 with the arctan function. With the increasing preheating temperature, the absolute error for the hardness calculation decreases.

Conclusion
In this investigation the potential of welding simulation has been shown. The combination of validating the FEM-model with measurements and the use of calculationmodels can provide solid results. In Fig. 10 the calculated and the measured hardness are compared related to the preheating temperature. With this method it is possible to predict preheating temperatures for specific materials. Modern-day computers and userfriendly simulation software make similar calculations possible in reasonable times. The prediction of the proper preheating temperature can be realized with comparably low effort. Moreover, the experimental validation can be done via microstructural investigations of the weld seam geometry and a hardness test. No further analysis is necessary. The essential results of this study are the following: & Preheating can decrease the maximum hardness within the weld seam of S690 high strength steel for high power laser welding & The error of the chemical composition measurement has a major effect on the calculated hardness. An accurate measurement is required to provide suitable results for the hardness calculation & The calibration of FEM-model with experimental dimension is needed to provide realistic results