Thermal analysis of Wire Arc Additive Manufacturing process

Online since 30 March 2021

Article

Abstract

Wire arc additive manufacturing process (WAAM) is an innovative technology that offers freedom in terms of designing functional parts, due to its ability to manufacture large and complex workpieces with a high rate of deposition. This technology is a metal AM process using an electric arc heat source. The parts manufactured are affected by thermal residual stresses due to high-energy input between wire and workpiece despite numerous advantages with this technology. It could cause severe deformation and change the global mechanical response. A 3D transient thermal model was created to evaluate the thermal gradients and fields during metal deposition. The material used in this study is a steel alloy (S355JR-AR). This numerical model takes into account the heat dissipation through the external environment and the heat loss through the cooling system under the base plate. Birth-element activation strategy was used to generate warm solid part following the movement of the heat source. The metal deposition is defined with constant welding speed. Goldak model was used to simulate the heat source in order to have a realistic heat flow distribution. Results were in concordance for thermal cycles at different points comparing with experimental results issued from bibliography in terms of: (1) Temperature maximum, (2) Thermal cycles and (3) Cooling gradient phase. This study enabled to check the numerical model and used as a predictive tool

Keywords

Additive Manufacturing, Wire Arc Additive Manufacturing, Heat transfer, Finite Element Modelling

Table of contents

Text

1 Introduction

Additive manufacturing (AM) was devoted initially to prototyping. However, this process has become a serious competitor for the production of metal parts with technological evolution [1].

Wire arc additive manufacturing (WAAM) has become a promising alternative to conventional manufacturing processes due to its high deposition rate (1-4 kg/h) [2], environmentally friendly processes [3], productivity [4], cost reduction [1] and reduction in the buyto-fly ratio (from around 10-20 for a conventionally machined component to 1 for WAAM) [5].

There are a variety of WAAM technologies that can be classified in 4 groups depending on the welding generator used: gas tungsten arc welding (GTAW), gas metal arc welding (GMAW), cold metal transfer (CMT) and plasma arc welding (PAW) [6][7]. WAAM is a process that allows parts to be produced by successive addition of material in the form of layers composed of several weld-lines. The manufacturing of a part layer by layer presents a freedom of design. Therefore, this process allows the realization of parts with complex shapes [8] but heat transfer need to be evaluated before completion.

This process has distortion problems such as strain and crack formation [8][9][10] due to high gradient temperature field and high rate cooling. Numerical simulation of this process is hence necessary to evaluate temperature field. It can be link to strain and mechanical properties. The objective could be to optimize parameters on the deposition process using the numerical predictive tool. Modeling of WAAM process is a complex step. It is necessary to study the couplings between physical domains such as metallurgy, mechanics and thermic [11].

In this study, a model of the thermal behavior is developed. The process studied is the CMT (cold metal transfer) process. The main aim of this study is to determine the evolution of the temperature field during the deposition of three weld-lines (wall) and to analyze the thermal cycles of the weld-line.

2 Material and geometry specifications of the WAAM wall

The geometry of the base plate is shown in Fig.1 and the dimensions are 200x60x12 mm3.The wall is built in the middle of the base plate and in the direction of its length. The geometry of these weld-lines is modeled by a rectangular section. It is defined by the width and height of deposit to simplify the problem. The assumption is based on experimental results [12]. In this project, metallurgical transformations are not taken into account. To simplify the problem, the thermo-mechanical coupling can be reduced to a weak coupling, which takes into account only the influence of thermal loading on the mechanical behavior of the part to be manufactured. Hence the mechanical and thermal analysis can be dissociated[13][11]. The thermal analysis can be used as input data for the mechanical analysis (residual stresses and strains).

Image 1000000000000642000001F35708F23358C84A04.png

Fig.1. Geometry used to simulate WAAM process

In this study, three weld-lines are simulated. In the experimental work illustrated in literature studies[12], the first weld-lines have a width less than 5 mm. It is due to quick heat flow (base plate conduction and cooling system). The out heat flow becomes weaker and stable after few layers. However, according to J.Ding et al.[12], this change in width does not have a significant influence on the results. Therefore, a constant width and a constant height is used for all weld-lines.

The base plate is made of S355AR-JR steel. Its chemical composition is summarized in Table 1 and the chemical composition of the wire to be used is summarized in Table 2. This steel is widely used in AM. The high cooling rate does not influence strongly the mechanical properties [12].

Table 1. Chemical composition of the S355JR baseplate (in wt. %)

Table 1. Chemical composition of the S355JR baseplate (in wt. %)

Table 2. Chemical composition of the consumable electrode (in wt. %)

Table 2. Chemical composition of the consumable electrode (in wt. %)

2.1 Thermo-physical properties

In order to solve the equations for solid and liquid phase, it is necessary to define the thermo-physical properties (density, thermal conductivity, latent heat and specific heat). The same properties are used for the wall and base plate materials.

Thermal conductivity and specific heat depend on temperature. Work of Michaleris and Debiccari [14] is used as a reference to determine these thermo-physical properties. As shown in Fig.2, the specific heat decreases when the temperature is above 723°C due to phase transformation. A high thermal conductivity is used for temperatures above 1500°C.

A constant value of density is used in the numerical simulation (equal to 7860 Kg.m-3). The latent heat of fusion is equal to 270 KJ.Kg-1 between the solidus temperature (1450°C) and the liquidus temperature (1500°C) [14].

Image 1000000000000611000001CA55F83D7DE29A62ED.png

Fig. 2. Thermal properties of S355[14]: a) Thermal conductivity; b) Specific heat

2.2 Thermocouple positions

In the work of Ding et al.[12], temperature evolution is carried out from thermocouples placed during the experiment. The coordinates are shown Fig.3. Temperature variable were placed at the same position for the comparison between experimental and numerical studies.

Image 1000000000000681000001D59F4FB1CEADEB377D.png

Fig. 3. The coordinates of the thermocouples[12]

3 Thermal finite element model

The numerical model is a transient FE analysis. The heat source is mobile in the x-axis direction and the workpiece remains fixed. It is based on the modelling of physical phenomena involved in the manufacturing process, such as heat transfer.

The finite element (FE) software ABAQUS is used to develop the thermal model. The objective was to have a predictive tool for temperature distribution which could be coupled to metallurgy phase or mechanical behavior (residual stress). Mesh is structured with linear hexagonal elements with 8 nodes DC3D8 (8-node linear heat transfer brick) [18]. A refined meshes are made in three areas since high thermal gradient:

  • Along the x-direction, three elements were used to cover the length of the melt pool. The total number of elements along the direction of the deposit is equal to 80 elements with the size of 2.5 mm length.

  • Along the y-direction, six elements were used to describe the melt pool with the size of 0.833 mm.

  • Along the z-direction, three elements were used to describe the temperature variation with the size of 0.666 mm.

For the base plate, the constructed mesh is coarse since thermal gradients are low in this region. A representation is shown in Fig.4. The mesh dimension is optimized checking variance stabilization.

Image 100000000000055700000196C78A40D3BACFFF3D.png

Fig. 4. The mesh of the transient model

The Finite-element activation strategy is used to simulate liquid material input. Initially, the wall is completely constructed of inactive elements, i.e. the mechanical and thermal properties are zero. An element is activated if it crosses the electric arc by attributing thermal and mechanical properties to it, otherwise it remains inactive [15]. This strategy has a disadvantage in terms of computational time but it allows a very realistic description of the process.

The liquid material deposition by electric arc is simulated as a volumetric heat source that moves along x-axis, based on Goldak's model (1), (2) [16]. This source is integrated with ABAQUS DFLUX subroutines.

Image 10000000000005C9000000ED0EE8033917AE30D1.png

Where a is the width of the ellipsoid, b is the height of the ellipsoid, cf is the length of the front ellipsoid and cr is the length of the rear ellipsoid. f is a coefficient that allows to impose higher thermal gradients at the front of the torch than at the rear, f=ff at the front

(x > xc), f=fr at the rear (x < xc) with ff+ff =2. For reasons of continuity, a condition (3) should be imposed on the coefficient C which also takes 2 different values at the front and at the rear of the heat source [12].

The term Q0=ŋ.I.U represents the maximum power where ŋ is the process efficiency, I is the current and U is the voltage, all parameters are summarized in Table 3.

Table. 3. Parameters used in the Goldak model[12]

Table. 3. Parameters used in the Goldak model[12]

During the deposition phase, the weld-line is introduced into the model using the element activation method with the 'MODEL CHANGE' option of ABAQUS. Initially, all steps are inactive (phase pre-step) and then the steps will be activated, step-by-step (weld line-i), by using an activation time, which depends on the welding speed. Smaller the size of the step is, closer the simulation results are to the experimental results but longer the computational time is. The size of the Step is equal to 10 mm and the activation time of each Step is equal to 1.2 second (For the inter-layer phase, the duration of the step is equal to 400 seconds, which is the inter-layer time required for a maximum temperature of 50°C. It is representative of a water-cooled aluminum machine plate. An organigram of the complete process is shown Fig.5.

Image 10000000000006230000035CD3D7BFBC68CA2A8F.png

Fig. 5. Organigram of the numerical simulation

4 Results

A comparison between the predicted data obtained by this study and the experiment results of Ding et al.[12] is presented in Fig.6. It shows numerical transient and the experimental thermal histories at one of the measurement positions on the base plate (P1) that are shown in Fig.3. The predicted temperatures were extracted from the nodal points. The gap between transient model and experimental result is equal to 5.6 % for the first peak, 1.3 % for the second peak and 8.6 % for the third peak. Predicted heating and cooling rates are correlated with the experimental results. A comparison of the five positions will be discussed.

Image 100000000000062F0000019595283F0352EAB5D5.png

Fig. 6. Temperature verification on the measuring positions: a) Thermal cycle found by the transient model (P1); b) Thermal cycle detected by the thermocouple (P1)[12].

Transient thermal cycles at five points (Fig.7.a) located on a cross section (x=100; z=0) are illustrated in Fig.7.b. All point are equidistant, however the thermal gap decrease along the y-axis from ≈ 250°C between C1 and C2 to less than 30°C between C4 and C5. The thermal gradient decrease until it reaches a few degrees after 30 seconds.

Image 100000000000063E000001D2CEF57AB27C50A419.png

Fig. 7. Transient thermal cycles at different points located on a cross section perpendicular to the x-axis: a) Positions of measurement points; b) Transient thermal cycles.

The temperature cycles of eleven points (Fig.8.a) located along the weld-axis (y=0; z=0) are shown in Fig.8.b. The maximum temperature reached remains constant (≃2200°C) for all points after about 20 mm (distance travelled after the first 2.4 second of deposition) except the point at the end of the weld-line (x=200 mm), which corresponds to the extinction zone of the heat source. It reaches a slightly higher maximum temperature (≃200°C).

Image 100000000000065A000001EEDBF616C61E6336FE.png

Fig. 8. Thermal cycles during the deposition along the x-axis: a) Positions of measurement points; b) Transient thermal cycles.

The thermal cycles at different heights of the deposited wall were acquired from the mid-length cross section (Fig.9.a) of the model, and shown in Fig.9.b. The maximum temperature increases after each deposition, it increases by 6 % between weld-line 1 and weld line 2, and by 3 % between weld-line 2 and weld-line 3. The peak of the first weld-line is lower due to the influence of the cooling system under the base plate.

Image 100000000000062700000387DAA7E07EC7910032.png

Fig. 9. Thermal cycles at different points located in the middle of the three weld-lines: a) Positions of measurement points; b) Thermal cycle in the middle of the first weld-line; c) Thermal cycle in the middle of the second weld-line; d) Thermal cycle in the middle of the third weld-line.

5 Discussion

The process starts at time t=0 s with an initial temperature equal to 20°C. Then the deposition phase begins, the heat source is applied, the temperature gradually increases at the front and the rear of the melt pool where the temperature is higher than the melting temperature of the material used (1500°C) as shown in Fig.10. A high thermal gradient is generated from the beginning of the deposition phase, in fact the maximum temperature varies from 20°C to 2300°C at the time t=0.298 second.

Image 10000000000005540000010DDEEF33DD8CA2DEE3.png

Fig. 10. Goldak double ellipsoidal heat source model

By analyzing the thermal cycles illustrated in Fig.6:

  • A good concordance in the temperature rise during the deposition phase in the different measurement points.

  • A similar evolution between the simulated temperatures and the measured temperatures.

  • A slight difference between the maximum simulated and measured temperatures.

The comparison shows that the numerical model developed can predict the temperature evolution during the WAAM process with reasonable accuracy. Indeed, the Table 4 shows the temperature gaps |𝑇𝑡𝑟𝑎𝑛𝑠𝑖𝑒𝑛𝑡 − 𝑇𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙| between the transient model and the experimental results[12] in each point of measurement during the deposition of three layers. All gaps do not exceed 25°C except in the case of point P5, in which the gap reaches 150°C (21 %) during the deposition of the 3rd weld-line. This difference may be due to fluid mechanics phenomena which is not taken into account in this study[17].

Table 4. Thermal gaps

Table 4. Thermal gaps

Fig.7 shows the temperature evolution as a function of time for five points located at different distances from the weld-line. It is noticeable that the thermal gradient is high in a cross section. This gradient decreases deeply after a few seconds, it is only a few degrees as shown in Fig.7. The farther the point is along the y-axis, the lower the maximum temperature reached. The maximum is reached after a higher time when the point is further away from the source. A similar thermal behavior has been shown by the experimental determination of the temperature evolution during welding presented by A.Ortega [5]. This delay is due to the variation of the heating speed along the y-axis.

Thermal cycles along the x-axis show that a quasi-stationary regime is established after a few seconds at a speed of 8.33 mm/s, where the maximum temperature remains constant, which explains the poor wettability at the beginning of the deposition [5] as shown in Fig.11. The maximum temperature reached for the point at the end of the first layer is higher than the maximum temperature reached in the quasi-stationary regime due to edge effect.

Image 10000000000005CC000000EF4E5FA19A0B4E8E14.png

Fig. 11. An example of a WAAM sample fabricated using mild steel wire

Fig.9 illustrates the thermal cycles recorded in the middle of the three deposited weld-lines. It is noticeable that the maximum temperature increases progressively until it reaches about 2480°C during the 3rd weld-lines. According to previous research [2] the maximum temperature stabilizes after 10 - 15 layers. Indeed, heat transfers by convection and radiation balance the heat input generated by the electric arc.

6 Conclusion

A thermal model was developed and compared with experimental results. The objective was to validate the model as a predictive tool. The results present the temperature evolution during the deposition phase, the inter-layer phase and the cooling phase. The evolution of the thermal gradient layer by layer can be analyzed. The simulation showed that:

  1. Transient models can accurately predict the heating and cooling cycles compared to experimental results.

  2. A quasi-stationary regime is established during each layer after a few seconds (3-5s), during which the temperature field is stable. Deposited layer keeps hence the same height and width.

  3. The instability of the melt pool at the beginning of each layer leads to an accumulation of material or a ‘macro-drop’ which causes a poor wettability.

  4. The superposition of layers generates a temperature increased layer by layer.

  5. At a cross section, the high thermal gradient decreases quickly after few seconds.

Finally, this thermal model can be carried out for thermo-mechanical prediction, in order to analyze the evolution of the residual stresses and strain.

Bibliography

[1] T. Artaza and A. Suárez, “Wire Arc Additive Manufacturing of Mn4Ni2CrMo Conference Steel : Comparison of Mechanical and Metallographic Properties of PAW Wire Arc Additive Manufacturing of Mn4Ni2CrMo Steel : and GMAW Comparison of Mechanica,” Procedia Manuf., vol. 41, pp. 1071–1078, 2020, doi: 10.1016/j.promfg.2019.10.035.

[2] V. Querard, “Réalisation de pièces aéronautiques de grandes dimentions par fabrication additive WAAM,” 2019.

[3] R. Tangestani, G. H. Farrahi, M. Shishegar, B. P. Aghchehkandi, S. Ganguly, and A. Mehmanparast, “Effects of Vertical and Pinch Rolling on Residual Stress Distributions in Wire and Arc Additively Manufactured Components,” J. Mater. Eng. Perform., vol. 29, no. 4, pp. 2073–2084, 2020, doi: 10.1007/s11665-020-04767-0.

[4] W. Liu, C. Jia, M. Guo, J. Gao, and C. Wu, “Compulsively constricted WAAM with arc plasma and droplets ejected from a narrow space,” Addit. Manuf., vol. 27, no. June 2018, pp. 109–117, 2019, doi: 10.1016/j.addma.2019.03.003.

[5] A. G. Ortega, “Prototypage rapide de pièces en alliage d’aluminium : étude du dépôt de matière et d’énergie lors de la fusion à l’arc d’un fil par le procédé MIG-CMT,” 2018.

[6] A. M. Waam and J. P. Oliveira, “Current Status and Perspectives on Wire and Arc Additive Manufacturing (WAAM),” 2019, doi: 10.3390/ma12071121.

[7] B. Wu et al., “A review of the wire arc additive manufacturing of metals: properties, defects and quality improvement,” J. Manuf. Process., vol. 35, no. February, pp. 127–139, 2018, doi: 10.1016/j.jmapro.2018.08.001.

[8] F. Montevecchi, G. Venturini, A. Scippa, and G. Campatelli, “Finite element modelling of Wire-Arc-Additive-Manufacturing process,” Procedia CIRP, vol. 55, pp. 109–114, 2016, doi: 10.1016/j.procir.2016.08.024.

[9] T. Shao, “Welding Thermal Simulation and Metallurgical Characteristics Analysis in WAAM for 5CrNiMo Hot Forging Die Remanufacturing,” Procedia Eng., vol. 207, pp. 2203–2208, 2017, doi: 10.1016/j.proeng.2017.10.982.

[10] S. Leuders et al., “On the mechanical behaviour of titanium alloy TiAl6V4 manufactured by selective laser melting: Fatigue resistance and crack growth performance,” Int. J. Fatigue, vol. 48, pp. 300–307, 2013, doi: 10.1016/j.ijfatigue.2012.11.011.

[11] J. Bergheau, “Modélisation numérique des procédés de soudage,” pp. 1–15, 1990.

[12] J. Ding et al., “Thermo-mechanical analysis of Wire and Arc Additive Layer Manufacturing process on large multi-layer parts,” Comput. Mater. Sci., vol. 50, no. 12, pp. 3315–3322, 2011, doi: 10.1016/j.commatsci.2011.06.023.

[13] C. Cambon, I. Bendaoud, C. Bordreuil, S. Rouquette, and F. Soulie, “Modélisation thermomécanique multicouche du procédé WAAM,” 2019.

[14] P. Michaleris and A. Debiccari, “Prediction of welding distortion,” Weld. J. (Miami, Fla), vol. 76, no. 12, pp. 172-s, 1997.

[15] M. Chiumenti et al., “Numerical modelling and experimental validation in Selective Laser Melting,” Addit. Manuf., vol. 18, pp. 171–185, 2017, doi:10.1016/j.addma.2017.09.002.

[16] J. Goldak, A. Chakravarti, and M. Bibby, “A New Finite Element Model for Welding Heat Sources,” Scand. J. Urol., vol. 52, no. 1, pp. 1–7, 2018, doi: 10.1080/21681805.2017.1363816.

[17] S. Cadiou, “Modélisation magnéto-thermohydraulique de procédés de fabrication additive arc-fil ( WAAM ), ” HAL Id : tel-02496906, 2020.

[18] ABAQUS Documentation, V 6.7, http://130.149.89.49:2080/v6.7/books/usb/default.htm?startat=pt06ch22s01ael03.html.

Illustrations

Table 1. Chemical composition of the S355JR baseplate (in wt. %)

Table 1. Chemical composition of the S355JR baseplate (in wt. %)

Table 2. Chemical composition of the consumable electrode (in wt. %)

Table 2. Chemical composition of the consumable electrode (in wt. %)

Table. 3. Parameters used in the Goldak model[12]

Table. 3. Parameters used in the Goldak model[12]

Table 4. Thermal gaps

Table 4. Thermal gaps

References

Electronic reference

Mohamed Belhadj, Sana Werda, Asma Belhadj, Robin Kromer and Philippe Darnis, « Thermal analysis of Wire Arc Additive Manufacturing process  », ESAFORM 2021 [Online], Online since 30 March 2021, connection on 06 December 2024. URL : https://popups.uliege.be/esaform21/index.php?id=4095

Authors

Mohamed Belhadj

University of Bordeaux, CNRS, ENSAM, I2M Bordeaux, 351 cours de la Libération, Talence, France.

Corresponding author: Mohamed Belhadj. E-mail address: mohamed.belhadj@u-bordeaux.fr

Sana Werda

By this author

Asma Belhadj

Robin Kromer

Philippe Darnis

By this author

Details

Title
Thermal analysis of Wire Arc Additive Manufacturing process
Language
en
Author, co-author
Mohamed Belhadj, Sana Werda, Asma Belhadj, Robin Kromer and Philippe Darnis,
Publication date
14 April 2021
Journal title
ESAFORM 2021
Copyright
CC-BY
DOI
10.25518/esaform21.4095
Permanent URL
https://popups.uliege.be/esaform21/index.php?id=4095

Statistics

Views : 1938 (3 ULiège)
Downloads : 0 (0 ULiège)

Cite

Belhadj, M., Werda, S., Belhadj, A., Kromer, R., & Darnis, P. (2021). Thermal analysis of Wire Arc Additive Manufacturing process . Paper presented at ESAFORM 2021. 24th International Conference on Material Forming, Liège, Belgique. doi: 10.25518/esaform21.4095