Finite Element Analysis on Welding-Induced Distortion of Automotive Rear Chassis Component

Welding-induced distortion is a major concern for the industrial joining practice. The welding-induced distortion at the weldment between the coupled torsion beam axle (CTBA) of the automotive rear chassis parts and trailing arm connected to the wheel hub axle module seriously affects the tow angle, camber angle, and caster of the wheel axle. In this paper, the welding process between CTBA and trailing arm was numerically analyzed via SYSWELD (i.e., a finite element analysis code), using the material properties predicted via J-MatPro SW, material properties software, considering the thermo-mechanical-metallurgical properties of materials. From the numerical study for the weldment of CTBA and trailing arm, we predicted the welding and thermal distortions, temperature variation, microstructure, and residual stress at the concerned area. As a result, the temperature of the welded area was predicted between 102 ◦C and 840 ◦C at the end of weldment and converged to room temperature after 1000 s. The maximum portion of the martensite structure at the welded area was expected to be 55%. The expected distortions of the trailing arm after the weldment were 0.52 mm, −1.47 mm, and 0.44 mm in the x, y, and z-directions, respectively. Finally, the limitations of this research and recommendations are presented.


Introduction
Among the several methods employed in the automobile manufacturing process wherein tens of thousands of parts are assembled, welding is an immensely efficient methods for joining parts. Welding is widely used in various engineering applications such as automotive, ship building, aerospace, construction, thermal power plants, and pressure vessels, etc. [1][2][3][4].
The welding processes in the automobile manufacturing industry have quickly been replaced by robot technology due to a few workloads involved in the latter. Therefore, process variable optimization is one of the important factors in the welding process with robot technology. The welding methods used in the automobile manufacturing process include spot welding, connecting the vehicle exterior panel and body, ND: YAG (Yttrium Aluminum Garnet) laser, metal inert gas (MIG) welding, etc. Spot welding involves inducing a temperature increase in the material using electrical resistance. The YAG welding method involves a laser using a YAG rod as the laser medium. MIG welding involves generating an arc between the wire (uncoated bare wire) and the base material, with subsequent melting of the wire. SYSWELD [27], a welding software developed by ESI, introduced Gaussian-type heat source distribution and coupled thermal-mechanical-metallic (phase transformation) characteristics. This program has been widely used to control the welding deformation of structures, evaluate the effects of phase transformation on residual stress [16,17], and analyze the thermal-mechanical-metallurgical phenomena of the heat-affected zone (HAZ), etc. in the welding process. Bate et al. [28] suggested that SYSWELD predicts the thermal deformation of the welding process well, but to evaluate the residual stress, a kinematic hardening model is preferred over the isotropic hardening model.
The main purpose of performing computational simulation analysis for the welding process is to predict the welding deformation and residual stress through analysis in advance, and to design the optimal welding process to reduce their occurrence. Jang and Kang et al. [29,30] conducted an experiment on the welding process of automobile bracket products and an analysis study by applying the finite element method. They showed that the method of changing the jig clamping conditions and welding sequence or forced deformation method of the welded deformed product into a shape is effective in suppressing the warpage of the product within the assembly tolerance. Park et al. [31] obtained the residual welding stress by coupling the heat transfer analysis results to the thermal stress analysis for the welding process of automobile frames and compared the analysis and experimental results for welding deformation. Additionally, Lee et al. [32] predicted the welding-induced deformation and residual stress in the MIG welding of an automotive muffler by using the SYSWELD software. They showed that the simulation results were similar to the experimental results, and optimization of the welding procedure by SYSWELD improved the muffler productivity.
Although most research studies have proven that FEM simulation is a good alternative approach, obtaining the material properties related to high temperature such as strain hardening, heat coefficient, TTT curve, etc., remain the research topic.
In this study, the coupled torsion beam axle (CTBA), which is mainly connected to the rear wheel of a small and medium-sized vehicle directly and plays an important role in the vehicle's behavior, was selected because of the difficulty in obtaining the target shape after the welding process. Figure 1 shows the CTBA, reflecting a relatively stable turning behavior compared to the combination method of the lateral rod and the trailing arm because the vehicle wheel is directly connected to the CTBA. In the lateral rod method, since the overall length of the lateral rod does not change when inducing the up and down bounce of the vehicle body, it is difficult to obtain a stable movement of the vehicle body because the wheels move slightly from side to side. Additionally, when passing through unevenness concentrated on one wheel, the camber angle of the opposite wheel may change to the outside. However, in CTBA, since the up and down bounces are offset by the torsion of the CTBA itself, it shows a relatively stable body behavior compared to the lateral road method. Furthermore, when one-wheel passes through the unevenness, the force is distributed or transmitted to the opposite wheel, and the rolling of the car body is relatively reduced [23][24][25].
The last and most important process in the CTBA manufacturing process is the joining process with the trailing arm using MIG welding. In most cases, CTBA is manufactured differently from the initial design value during the manufacturing process due to the thermal deformation in the welding process above-mentioned. In this case, since it affects the toe and camber values of the wheels related to the driving performance of the vehicle, thermal deformation in the welding process must be strictly managed to secure the driving stability of the vehicle. To minimize thermal deformation in the welding process of the CTBA and the trailing arm, a method of restraining the deformation by fastening a jig at each important point of the CTBA and the trailing arm has been used. However, the restraint position of the jig and the time (cooling time) for biting the jig were intuitively determined in most cases, with some time wasted by not being analytical and systematic. The last and most important process in the CTBA manufacturing process is the joining process with the trailing arm using MIG welding. In most cases, CTBA is manufactured differently from the initial design value during the manufacturing process due to the thermal deformation in the welding process above-mentioned. In this case, since it affects the toe and camber values of the wheels related to the driving performance of the vehicle, thermal deformation in the welding process must be strictly managed to secure the driving stability of the vehicle. To minimize thermal deformation in the welding process of the CTBA and the trailing arm, a method of restraining the deformation by fastening a jig at each important point of the CTBA and the trailing arm has been used. However, the restraint position of the jig and the time (cooling time) for biting the jig were intuitively determined in most cases, with some time wasted by not being analytical and systematic.
In this study, for the purpose of analyzing the welding process of CTBA and trailing arm, the thermal-mechanical-metal characteristics of the material of CTBA and trailing arm were first identified. For this purpose, J-MatPro (Java-based Material Properties) [33], based on the CALPHAD method, was used. In addition, using SYSWELD, a finite element analysis program for welding quantified the thermal deformation at critical points during the welding process, temperature change over time, phase transformation, and residual stress. Through this welding process simulation analysis, the welding deformation was identified, the texture distribution of the material was predicted, and the welding residual stress related to the durability of the chassis parts was investigated.

Specification and Modeling of CTBA and Trailing Arms
In general, CTBA is ~1200 mm in length, ~100 mm in diameter, and ~2.5 mm in thickness. The central part of the CTBA forms a V-shape, and both ends form a square shape for fastening (welding) to the trailing arm. In this study, CTBA processed with 22MnB5, a hot forming material, was used. Hot forming generates automotive structure parts with a high strength of about 1500 MPa by heating a steel material, containing numerous B and Mn atoms, to a high temperature to form an austenite structure, and then forming a martensitic structure through a rapid cooling process such as die quenching.
Conversely, the trailing arm used a round or square tube with a tensile strength of 590 MPa. Tables 1-4 show the chemical composition and mechanical properties of the 22MnB5 material and the SAPH590 material. In Tables 2 and 4, E, YS, TS, El, n-value, and R-value represent the Young's modulus, initial yield stress, yield strength, elongation, strain hardening exponent, and plastic anisotropy ratio, respectively. In this study, for the purpose of analyzing the welding process of CTBA and trailing arm, the thermal-mechanical-metal characteristics of the material of CTBA and trailing arm were first identified. For this purpose, J-MatPro (Java-based Material Properties) [33], based on the CALPHAD method, was used. In addition, using SYSWELD, a finite element analysis program for welding quantified the thermal deformation at critical points during the welding process, temperature change over time, phase transformation, and residual stress. Through this welding process simulation analysis, the welding deformation was identified, the texture distribution of the material was predicted, and the welding residual stress related to the durability of the chassis parts was investigated.

Specification and Modeling of CTBA and Trailing Arms
In general, CTBA is~1200 mm in length,~100 mm in diameter, and~2.5 mm in thickness. The central part of the CTBA forms a V-shape, and both ends form a square shape for fastening (welding) to the trailing arm. In this study, CTBA processed with 22MnB5, a hot forming material, was used. Hot forming generates automotive structure parts with a high strength of about 1500 MPa by heating a steel material, containing numerous B and Mn atoms, to a high temperature to form an austenite structure, and then forming a martensitic structure through a rapid cooling process such as die quenching.
Conversely, the trailing arm used a round or square tube with a tensile strength of 590 MPa. Tables 1-4 show the chemical composition and mechanical properties of the 22MnB5 material and the SAPH590 material. In Tables 2 and 4, E, YS, TS, El, n-value, and R-value represent the Young's modulus, initial yield stress, yield strength, elongation, strain hardening exponent, and plastic anisotropy ratio, respectively.   Figure 2 shows the actual welding process between the CTBA and the trailing arm, and the numbered part shows the position of the jig fastened during the welding processes.   Figure 2 shows the actual welding process between the CTBA and the trailing arm, and the numbered part shows the position of the jig fastened during the welding processes. To analyze the welding process joining CTBA and the trailing arm, FE modeling was performed using shell elements, as shown in Figure 3a,b. For accurate finite element analysis, only the middle surface was taken, and the welding part was modeled. In Figure 3b, part A is fastened to the car body using a bush pipe, part B is connected to the wheel through the hub bearing via the spindle bracket, and part C is where the spring of the suspension is coupled. Therefore, parts A and B affect the toe, camber, and caster angles, thus, the dimensional accuracy of these two parts is particularly important. To analyze the welding process joining CTBA and the trailing arm, FE modeling was performed using shell elements, as shown in Figure 3a,b. For accurate finite element analysis, only the middle surface was taken, and the welding part was modeled. In Figure 3b, part A is fastened to the car body using a bush pipe, part B is connected to the wheel through the hub bearing via the spindle bracket, and part C is where the spring of the suspension is coupled. Therefore, parts A and B affect the toe, camber, and caster angles, thus, the dimensional accuracy of these two parts is particularly important. The upper figure of Figure 4a is an enlarged view of the welding part. This part, which is the tip of the CTBA, is joined to the trailing arm, and the welding process is performed along the tip periphery. Figure 4b shows the welding start and end points, line to be welded, welding reference line, and welding start and end positions. The positions of the welding start and the end are the same, but here, they are shown as separate for the The upper figure of Figure 4a is an enlarged view of the welding part. This part, which is the tip of the CTBA, is joined to the trailing arm, and the welding process is performed along the tip periphery. Figure 4b shows the welding start and end points, line to be welded, welding reference line, and welding start and end positions. The positions of the welding start and the end are the same, but here, they are shown as separate for the convenience of indicating the welding path.  The upper figure of Figure 4a is an enlarged view of the welding part. This which is the tip of the CTBA, is joined to the trailing arm, and the welding proce performed along the tip periphery. Figure 4b shows the welding start and end points to be welded, welding reference line, and welding start and end positions. The posit of the welding start and the end are the same, but here, they are shown as separate fo convenience of indicating the welding path.

Setting Parameters for Welding Simulation
The welding process in the actual production line consists of the following steps: (1) Fasten the CTBA and the trailing arm in place; (2) As shown in Figure 3, restrain the bush pipe and the connecting part A, the hub bearing and the connecting part B, and the center of the CTBA; (3) In this state, weld for 40 s along the welding line in Figure 4 with a MIG robot welding machine; and (4) After welding is completed, cool in the air for 1000 s with all restraints dismantled. Table 5 shows the detailed welding conditions and parameters. The temperature of HAZ used for the analysis was set to 867 °C, an average of the temperatures measured at the locations of 20 actual welds. The efficiency coefficient was assumed to be 0.6. Table 5. Parameters of the welding process.

Setting Parameters for Welding Simulation
The welding process in the actual production line consists of the following steps: (1) Fasten the CTBA and the trailing arm in place; (2) As shown in Figure 3, restrain the bush pipe and the connecting part A, the hub bearing and the connecting part B, and the center of the CTBA; (3) In this state, weld for 40 s along the welding line in Figure 4 with a MIG robot welding machine; and (4) After welding is completed, cool in the air for 1000 s with all restraints dismantled. Table 5 shows the detailed welding conditions and parameters. The temperature of HAZ used for the analysis was set to 867 • C, an average of the temperatures measured at the locations of 20 actual welds. The efficiency coefficient was assumed to be 0.6. Among the welding types, overlapping welding is a method of welding overlapping parts by adding two plates, as shown in Figure 6a. Figure 6b shows the type of heat source used for the MIG welding adopted in this study.
Fasten the CTBA and the trailing arm in place; (2) As shown in Figure 3, restrain the bush pipe and the connecting part A, the hub bearing and the connecting part B, and the center of the CTBA; (3) In this state, weld for 40 s along the welding line in Figure 4 with a MIG robot welding machine; and (4) After welding is completed, cool in the air for 1000 s with all restraints dismantled. Table 5 shows the detailed welding conditions and parameters. The temperature of HAZ used for the analysis was set to 867 °C, an average of the temperatures measured at the locations of 20 actual welds. The efficiency coefficient was assumed to be 0.6. Among the welding types, overlapping welding is a method of welding overlapping parts by adding two plates, as shown in Figure 6a. Figure 6b shows the type of heat source used for the MIG welding adopted in this study.

Material Properties for High Temperature
For accurate welding process analysis, the database for thermal-physical properties of 22MnB5 and SAPH590 materials was obtained using the J-MatPro program.
J-MatPro provides the thermo-physical properties of a material such as thermal conductivity, specific heat, Young's modulus, flow curve, phase transform characteristics, continuous cooling transform curve (CCT), and time-temperature transform (TTT) by inputting the chemical composition of the material. The obtained database of the thermalphysical properties is an important variable for the analysis of the welding process. Figure 7a-e shows the TTT diagram, CCT diagram, thermal conductivity, yield stress, and Young's modulus, which are the thermal-physical properties of 22MnB5 obtained using J-MatPro, respectively. Detailed information on the thermo-physical properties of the SAPH590 material is described in [27], omitted here for reasons of space.
putting the chemical composition of the material. The obtained database of the thermalphysical properties is an important variable for the analysis of the welding process. Figure 7a-e shows the TTT diagram, CCT diagram, thermal conductivity, yield stress, and Young's modulus, which are the thermal-physical properties of 22MnB5 obtained using J-MatPro, respectively. Detailed information on the thermo-physical properties of the SAPH590 material is described in [27], omitted here for reasons of space.

Temperature Analysis
As in the previous section, after modeling the welding process between the CTBA and the trailing arm, the welding process was analyzed using SYSWELD, and the results for important properties were derived. Figure 8 shows the results of calculating the temperature change according to the process time based on the position of the welding start point. As shown in the figure, at

Temperature Analysis
As in the previous section, after modeling the welding process between the CTBA and the trailing arm, the welding process was analyzed using SYSWELD, and the results for important properties were derived. Figure 8 shows the results of calculating the temperature change according to the process time based on the position of the welding start point. As shown in the figure, at the beginning of the welding process, the joint between CTBA and trailing arm is heated to 1500 • C or higher, that is, to the melting point or higher, and when the welding time passes 40 s, it drops to 800 • C or lower. After welding is completed, the welding jig is dismantled and cooled to room temperature for up to 1000 s.

Temperature Analysis
As in the previous section, after modeling the welding process between the CTBA and the trailing arm, the welding process was analyzed using SYSWELD, and the results for important properties were derived. Figure 8 shows the results of calculating the temperature change according to the process time based on the position of the welding start point. As shown in the figure, at the beginning of the welding process, the joint between CTBA and trailing arm is heated to 1500 °C or higher, that is, to the melting point or higher, and when the welding time passes 40 s, it drops to 800 °C or lower. After welding is completed, the welding jig is dismantled and cooled to room temperature for up to 1000 s. After welding completion, if the welding jig is dismantled, the welded part is cooled to RT after 1000 s.
As shown in Figure 8, the cooling rate for 40 s while restraining the thermal deformation of the material with the welding jig was about 17.5 °C/s. The cooling rate was about 3.6 °C/s between 40 and 100 s after dismantling the welding jig. Therefore, as predicted from the CCT curve in Figure 8, a mixed structure of martensite, ferrite, and bainite must be obtained at the welding starting point.
At 100 s after the welding completion and dismantling of the welding jig, CTBA showed a linear decrease in temperature by convective heat transfer (cooling rate: ~0.64 °C/s) After welding completion, if the welding jig is dismantled, the welded part is cooled to RT after 1000 s.
As shown in Figure 8, the cooling rate for 40 s while restraining the thermal deformation of the material with the welding jig was about 17.5 • C/s. The cooling rate was about 3.6 • C/s between 40 and 100 s after dismantling the welding jig. Therefore, as predicted from the CCT curve in Figure 8, a mixed structure of martensite, ferrite, and bainite must be obtained at the welding starting point.
At 100 s after the welding completion and dismantling of the welding jig, CTBA showed a linear decrease in temperature by convective heat transfer (cooling rate:~0.64 • C/s).
The temperature at~100 s from the starting point of welding rises to the melting temperature of the material and then rapidly decreases due to conduction and convective heat transfer from the heat source moving along the welding path. Therefore, the first 100 s of the cooling dominated by conduction and convection is a relatively short time compared to the total cooling time of 1000 s, although rapid. Therefore, the cooling behavior up to the first 100 s after the welding process is important in controlling the welding heat deformation of the CTBA products. Figure 9a,b shows the temperature of each part after the welding process. Figure 9a shows the temperature distribution around the weld area after 40 s immediately post welding, with the highest and lowest temperatures of 839.6 • C and 101.9 • C, respectively. Figure 9b shows the temperature distribution around the welding part after cooling for 1000 s immediately post welding with the highest and lowest temperatures of 23.8 • C and 20.3 • C. When performing the cooling process by convective heat transfer for about 1000 s, the welded member is cooled to the room temperature in all parts. Figure 9a,b shows the temperature of each part after the welding process. Figure 9a shows the temperature distribution around the weld area after 40 s immediately post welding, with the highest and lowest temperatures of 839.6 °C and 101.9 °C, respectively. Figure 9b shows the temperature distribution around the welding part after cooling for 1000 s immediately post welding with the highest and lowest temperatures of 23.8 °C and 20.3 °C. When performing the cooling process by convective heat transfer for about 1000 s, the welded member is cooled to the room temperature in all parts.  Figure 10 shows the fraction of the martensite structure obtained in the HAZ area where welding started post-welding analysis by SYSWELD. The part showing the highest fraction of martensitic structure among welded parts is a very small area where the maximum heat source near the welding start point is applied. The region had a martensite fraction of about 55% and other HAZ regions were about 35% or less.   Figure 10 shows the fraction of the martensite structure obtained in the HAZ area where welding started post-welding analysis by SYSWELD. The part showing the highest fraction of martensitic structure among welded parts is a very small area where the maximum heat source near the welding start point is applied. The region had a martensite fraction of about 55% and other HAZ regions were about 35% or less. Figure 9a,b shows the temperature of each part after the welding process. Figure 9a shows the temperature distribution around the weld area after 40 s immediately post welding, with the highest and lowest temperatures of 839.6 °C and 101.9 °C, respectively. Figure 9b shows the temperature distribution around the welding part after cooling for 1000 s immediately post welding with the highest and lowest temperatures of 23.8 °C and 20.3 °C. When performing the cooling process by convective heat transfer for about 1000 s, the welded member is cooled to the room temperature in all parts.  Figure 10 shows the fraction of the martensite structure obtained in the HAZ area where welding started post-welding analysis by SYSWELD. The part showing the highest fraction of martensitic structure among welded parts is a very small area where the maximum heat source near the welding start point is applied. The region had a martensite fraction of about 55% and other HAZ regions were about 35% or less.  In the air cooling of austenitic steel, the material is cooled by convective heat transfer with a low cooling rate, and the fraction of martensite structure does not exceed 50%.

Phase Transformation
In this study, cooling is accelerated by the convective heat transfer and heat conduction in most of the HAZ region during the cooling process after welding, but the cooling rate is slow, thus, the fraction of martensite structure is not high.
The cooling rate at the HAZ position is shown on the continuous cooling transformation curve in Figure 8. As predicted from the CCT diagram, the region contained a mixed structure of martensite, ferrite, and bainite.
To confirm the phase located at the HAZ position, scanning electron microscopy was performed. Figure 11 shows a picture of the structure at the central part of the HAZ area post-welding completion, which is composed of a mixed phase of martensite, ferrite, and bainite.
The cooling rate at the HAZ position is shown on the continuous cooling transformation curve in Figure 8. As predicted from the CCT diagram, the region contained a mixed structure of martensite, ferrite, and bainite.
To confirm the phase located at the HAZ position, scanning electron microscopy was performed. Figure 11 shows a picture of the structure at the central part of the HAZ area post-welding completion, which is composed of a mixed phase of martensite, ferrite, and bainite. Figure 11. Microstructure in the weld starting point.

Thermal Deformation
As mentioned before, the welding process is accompanied by deformation due to the difference in the thermal strain and transformation due to the local heat generated during the heating and cooling cycles. Figure 12a,b shows the thermal deformation at 40 s immediately post welding and 1000 s after dismantling the welding jig, respectively. As shown in Figure 12a, 40 s immediately after the welding process, the maximum thermal strain was 0.41 mm at the hub bearing connection, and the minimum thermal strain at the bush pipe connection was 0.04 mm. However, as shown in Figure 12b, 1000 s post-cooling, the maximum thermal deformation was 1.33 mm in the suspension spring-mounted part and the minimum thermal deformation was 0.13 mm in the middle of the trailing arm. The predicted deformation in the suspension spring-mounted part was within the maximum allowable thermal deformation amount in the actual production line.

Thermal Deformation
As mentioned before, the welding process is accompanied by deformation due to the difference in the thermal strain and transformation due to the local heat generated during the heating and cooling cycles. Figure 12a,b shows the thermal deformation at 40 s immediately post welding and 1000 s after dismantling the welding jig, respectively. As shown in Figure 12a, 40 s immediately after the welding process, the maximum thermal strain was 0.41 mm at the hub bearing connection, and the minimum thermal strain at the bush pipe connection was 0.04 mm. However, as shown in Figure 12b, 1000 s post-cooling, the maximum thermal deformation was 1.33 mm in the suspension spring-mounted part and the minimum thermal deformation was 0.13 mm in the middle of the trailing arm. The predicted deformation in the suspension spring-mounted part was within the maximum allowable thermal deformation amount in the actual production line. mixed structure of martensite, ferrite, and bainite.
To confirm the phase located at the HAZ position, scanning electron microscopy was performed. Figure 11 shows a picture of the structure at the central part of the HAZ area post-welding completion, which is composed of a mixed phase of martensite, ferrite, and bainite. Figure 11. Microstructure in the weld starting point.

Thermal Deformation
As mentioned before, the welding process is accompanied by deformation due to the difference in the thermal strain and transformation due to the local heat generated during the heating and cooling cycles. Figure 12a,b shows the thermal deformation at 40 s immediately post welding and 1000 s after dismantling the welding jig, respectively. As shown in Figure 12a, 40 s immediately after the welding process, the maximum thermal strain was 0.41 mm at the hub bearing connection, and the minimum thermal strain at the bush pipe connection was 0.04 mm. However, as shown in Figure 12b, 1000 s post-cooling, the maximum thermal deformation was 1.33 mm in the suspension spring-mounted part and the minimum thermal deformation was 0.13 mm in the middle of the trailing arm. The predicted deformation in the suspension spring-mounted part was within the maximum allowable thermal deformation amount in the actual production line.  Figure 13a-c also shows the amount of deformation in each direction (X, Y, Z) at 1000 s post-cooling. As shown in Figure 13, the amount of thermal deformation after the welding process was up to 0.53 mm at the C position in the X direction, −1.47 mm at the maximum at B in the Y direction, and up to 0.44 mm at C in the Z direction. The amount of thermal deformation increased along with the distance from the welding site.
Figure 13a-c also shows the amount of deformation in each direction (X, Y, Z) at 1000 s post-cooling. As shown in Figure 13, the amount of thermal deformation after the welding process was up to 0.53 mm at the C position in the X direction, -1.47 mm at the maximum at B in the Y direction, and up to 0.44 mm at C in the Z direction. The amount of thermal deformation increased along with the distance from the welding site.  Figure 14 shows the residual stress (Von mises stress) distribution after welding. In this figure, the residual stress is formed near HAZ due to its two welded parts. The stress remaining in the product is distributed between the maximum of 501 MPa and the minimum of 50 MPa 1000 s post welding. The residual stress can affect the deformation of the part after welding.  Figure 14 shows the residual stress (Von mises stress) distribution after welding. In this figure, the residual stress is formed near HAZ due to its two welded parts. The stress remaining in the product is distributed between the maximum of 501 MPa and the minimum of 50 MPa 1000 s post welding. The residual stress can affect the deformation of the part after welding. To confirm the effect of residual stress on the deformation of the part, the deformation contours are depicted in Figure 15. This figure compares the shape after thermal deformation during welding and the shape before the welding process (design shape) after 1000 s (i.e., corresponding to residual stress contour) have passed after the welding. To confirm the effect of residual stress on the deformation of the part, the deformation contours are depicted in Figure 15. This figure compares the shape after thermal deformation during welding and the shape before the welding process (design shape) after 1000 s (i.e., corresponding to residual stress contour) have passed after the welding. Additionally, the solid line in Figure 15 shows the shape before the welding process. The deformation of part A connected to the bush pipe was 0.8 mm, and the deformation of part B connected to the hub bearing via the spindle bracket was 1.2 mm. As in the previous results, the greater the distance from the welding part, the greater the amount of thermal deformation due to welding. To confirm the effect of residual stress on the deformation of the part, the deformation contours are depicted in Figure 15. This figure compares the shape after thermal deformation during welding and the shape before the welding process (design shape) after 1000 s (i.e., corresponding to residual stress contour) have passed after the welding. Additionally, the solid line in Figure 15 shows the shape before the welding process. The deformation of part A connected to the bush pipe was 0.8 mm, and the deformation of part B connected to the hub bearing via the spindle bracket was 1.2 mm. As in the previous results, the greater the distance from the welding part, the greater the amount of thermal deformation due to welding. If the deformation at the spindle bracket B position of the trailing arm is predicted through the finite element analysis of the welding process, the deformation generated in the welding must be post-processed and removed for coupling with the hub bearing. In this study, there must be no problem in securing the driving stability of the designed vehicle only after welding and removing the spindle bracket protruding 1.20 mm from the design shape after welding. In the actual production process, all welded CTBAs are inspected, and post-processing is usually performed to a diameter of about 1.5 mm.
At the manufacturing site of automobile parts, to reduce the number of machining man-hours in the manufacturing process and quality problems caused by machining errors, it is desirable to design a process that suppresses welding deformation within an allowable range, if possible. For this purpose, finite element analysis is meaningful in predicting the deformation in the welding process and to review the process design to reduce possible welding deformation.
In this study, we did not consider the effect of the deposition bodies considering the filler material deposition in the real welding process as discussed in [34], where each welding seam is sliced into several deposition bodies to predict the welding distortion realistically. This approach will be used in future studies. If the deformation at the spindle bracket B position of the trailing arm is predicted through the finite element analysis of the welding process, the deformation generated in the welding must be post-processed and removed for coupling with the hub bearing. In this study, there must be no problem in securing the driving stability of the designed vehicle only after welding and removing the spindle bracket protruding 1.20 mm from the design shape after welding. In the actual production process, all welded CTBAs are inspected, and post-processing is usually performed to a diameter of about 1.5 mm.
At the manufacturing site of automobile parts, to reduce the number of machining man-hours in the manufacturing process and quality problems caused by machining errors, it is desirable to design a process that suppresses welding deformation within an allowable range, if possible. For this purpose, finite element analysis is meaningful in predicting the deformation in the welding process and to review the process design to reduce possible welding deformation.
In this study, we did not consider the effect of the deposition bodies considering the filler material deposition in the real welding process as discussed in [34], where each welding seam is sliced into several deposition bodies to predict the welding distortion realistically. This approach will be used in future studies.

Conclusions
In this study, the following conclusions were drawn by analyzing the welding deformation behavior during the welding process between the vehicle rear wheel chassis component CTBA and the trailing arm using SYSWELD, a finite element analysis program that considers the physical properties of the thermo-mechanical-metallurgical properties.
(1) The welding process between the vehicle rear wheel chassis component CTBA and the trailing arm was simulated and analyzed using SYSWELD. (2) For the accurate analysis of the welding process, the thermo-mechanical-metallurgical properties of 22MnB5 were obtained using J-MatPro. (3) As a result of the finite element analysis, the welding site temperature immediately post welding was predicted to be a maximum of 840 • C and a minimum of 102 • C, while that after cooling for 1000 s was predicted to be a maximum of 24 • C and a minimum of 20 • C. (4) As a result of the finite element analysis, the martensite fraction of the welded region is a very small region where the maximum heat source near the welding start point is applied, but in most HAZ regions, the martensite tissue fraction is judged less thañ 35%. (5) As a result of the finite element analysis, the maximum thermal deformation of the trailing arm was 1.33 mm with the minimum 0.13 mm, and by direction, the maximum was 0.52 mm in the X direction, 1.47 mm in the Y direction, and 0.44 mm in the Z direction. (6) The finite element analysis of the welding process in this study can be used to design the optimal welding process or welding schedule to reduce the welding deformation of the trailing arm combined with the hub bearing.
In this study, computational simulation analysis of the CTBA welding process was performed to identify thermal deformation and residual stress. However, the simulation result was not validated due to the absence of the experimental results. Moreover, to secure the accuracy of the analysis and its application value to industry, it is necessary not only to measure the temperature and hardness of the welded part and compare it with the results of this analysis, but also to compare it with many simplified analysis techniques that do not consider the phase transformation characteristics. We intend to leave this study as a follow-up study by the authors.