Next Article in Journal
Surface Characterization and Corrosion Resistance of 36Cr-Ni-Mo4 Steel Coated by WC-Co Cermet Electrode Using Micro-Electro Welding
Next Article in Special Issue
Control of Porosity and Spatter in Laser Welding of Thick AlMg5 Parts Using High-Speed Imaging and Optical Microscopy
Previous Article in Journal
Microstructural, Mechanical, Texture and Residual Stress Characterizations of X52 Pipeline Steel
Previous Article in Special Issue
Effects of Pulsed Nd:YAG Laser Welding Parameters on Penetration and Microstructure Characterization of a DP1000 Steel Butt Joint
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FEM Simulation of Dissimilar Aluminum Titanium Fiber Laser Welding Using 2D and 3D Gaussian Heat Sources

1
Innovation Engineering Department, University of Salento, Via per Arnesano s.n., 73100 Lecce, Italy
2
DMMM Politecnico di Bari, Viale Japigia, 182, 70126 Bari, Italy
*
Author to whom correspondence should be addressed.
Metals 2017, 7(8), 307; https://doi.org/10.3390/met7080307
Submission received: 1 June 2017 / Revised: 2 August 2017 / Accepted: 8 August 2017 / Published: 10 August 2017
(This article belongs to the Special Issue Laser Welding)

Abstract

:
For a dissimilar laser weld, the model of the heat source is a paramount boundary condition for the prediction of the thermal phenomena, which occur during the welding cycle. In this paper, both two-dimensional (2D) and three-dimensional (3D) Gaussian heat sources were studied for the thermal analysis of the fiber laser welding of titanium and aluminum dissimilar butt joint. The models were calibrated comparing the fusion zone of the experiment with that of the numerical model. The actual temperature during the welding cycle was registered by a thermocouple and used for validation of the numerical model. When it came to calculate the fusion zone dimensions in the transversal section, the 2D heat source showed more accurate results. The 3D heat source provided better results for the simulated weld pool and cooling rate.

Graphical Abstract

1. Introduction

Laser welding is recognized as an effective process to weld metals with a laser beam of high-power, high-energy density. In fact, the power density of a laser beam is much higher than that of arc or plasma. Consequently, a deep narrow penetration weld can be effectively produced. These properties have made laser welding a suitable technology for weldments that are made from metals of different compositions and properties [1,2]. A dissimilar joint is as strong as the weaker of the two metals being joined, i.e., possesses sufficient tensile strength and ductility so that the joint will not fail in the weld, has good fatigue behavior [3].
Among them, Al/Ti dissimilar joints are of major interest in aeronautics and automotive applications, where weight reduction, coupled with high mechanical strength and corrosion resistance, are paramount. Between the different Al/Ti welding processes laser welding offers numerous advantages. Especially in aluminum alloys when used in keyhole mode improves the absorption of the beam due to the multiple reflections in the cavity [4]. Moreover, high energy density, high cooling and heating rate allow for reducing the importance of mixing and diffusion phenomena, and thus reduce the formation of intermetallic compounds in the case of dissimilar joints. An Al/Ti joint has a remarkably lower elongation due to the high residual stresses, which facilitate the crack ignition and propagation. Therefore, the quality depends heavily on the process parameters, which determine the magnitude of thermal stresses [5].
The selection of the welding parameters is crucial for obtaining a satisfactory quality weld. Residual stresses and temperature field in laser welding joints can be predicted by numerical analysis such as a finite element one. The Finite Element Method has been one of the performing techniques to predict the joint properties in the welding process, which involves thermal, metallurgical and mechanical phenomena.
The computation of thermal field relies strongly on the heat source model. Rosenthal was the first researcher who proposed a model for the heat source in welding [6]. He proposed an analytical solution considering a punctual or a line heat source. Since then, other more realistic models have been proposed. For arc welding, several heat source configurations have also been proposed. Two and three-dimensional approaches can be used. Zeng et al. described the thermal elastic-plastic analysis using finite element techniques to analyze the thermos-mechanical behavior and evaluate the residual stresses and welding distortion on the AZ31B magnesium alloy and 304L steel butt joint in laser-TIG hybrid welding [7]. A modified three-dimensional conical heat source was used for performing the simulation in arc welding [8,9]. In certain cases, the finite element model is integrated with other computational techniques like artificial intelligence to establish an automated and iterative optimization algorithm [10]. Zeng et al. calculated the thermal cycles and temperature distribution of MIG welding of 5A06 aluminum alloy structure during discontinuous welding. The finite element method transient heat transfer analysis was used to save computing time and improve calculation accuracy [11].
For the fiber laser, Casalino et al. [12] developed and applied a stationary process with a surface heat source model based on thermal load through several specific elements next to the welding line. Other researchers [13] proposed the combination of Gaussian distribution on the surface and distribution along the thickness to consider 3D distribution by applying the conical Gaussian heat source model. They found that 3D conical Gaussian heat distribution can obtain better results with high depth to width ratio (defined as the ratio between the weld penetration evaluated on the axis of fused zone, and the width of the welded seam in the horizontal direction on the sample surface). Nagel et al. [14] proposed some strategies for the optimization of the laser welding of high alloy steel sheets using two different heat sources.
In this paper, two-dimensional and three-dimensional heat distribution was used for welding laser simulation of dissimilar Al/Ti button joint. The objective of this study is to compare the two approaches to establish the best model. The simulated fusion zone was compared with the macrograph obtained from the experiment to calibrate the model. Then, the validation is based on the comparison of the temperature profile measured with thermocouples during the welding cycle.

2. Experimental Setup

A laser butt joint has been produced from two plates of aluminum and titanium (3 and 2 mm, respectively), according to the scheme of Figure 1. Table 1 and Table 2 report the chemical composition of the two alloys and the mechanical and thermo-physical properties at room temperature. Particularly, the thermal conductivity has been considered as a function of temperature (Table 3 and Table 4).
The experimental trials were carried out using an Ytterbium Fiber Laser System (IPG YLS-4000), with a maximum output power equal to 4 kW (IPG Laser, Barbuch, Germany). The laser beam was delivered through a 200 µm optical fiber with a Beam Parameter Product (BPP) equal to 6.3 mm·mrad, the product of a laser beam’s divergence angle (half-angle) and the radius of the beam at its narrowest point. The laser beam, whose wavelength was 1070.6 nm, has been focused continuously through a lens with focal distance of 250 mm producing a spot diameter of 0.4 mm on the workpiece surface. Argon and helium were employed as shielding gas with 10 L/min volumetric flow rate, particularly Argon has been employed on the upper surface and helium on the bottom surface. The laser beam axis was placed on the titanium side, some 1 mm far from the interface (laser offset). The welding parameters have been 1200 W power at 1000 m/min welding rate.
The microstructure of the weld was analyzed by optical microscopy (OM; Nikon Epiphot 200, Nikon, Tokyo, Japan) and Zeiss EVO scanning electron microscope (SEM, Zeiss, Oberkochen, Germany). For optical microscopy observation, the transversal sections of the samples were cut and prepared using the standard metallographic grinding and polishing techniques and attached using Keller reagent (95 mL H2O, 2.5 mL HNO3, 1.5 mL HCl, HF 1 mL). The dimension of the fusion zone (FZ) was evaluated using NIS-Element software (version 4.5, Nikon, Tokio, Japan, 2016)) for the image analysis. NIS-Elements is a Nikon software supplied with Epiphot 200 OM. The software is tailored to facilitate image capture, object measurement and counting. Vickers microhardness profile (0.3/15) was collected using a Vickers Affri Wiky 200JS2 microhardness tester (Affri, Wood Dale, IL, USA) at half of weld cross section thickness. The distance between indentations was equal to 300 μm. The hardness at the Al/Ti fusion zones interface, where an intermetallic compound layer was observed, had been identified with nanoindentation (0.01/15) using Leica VMHT (Leica Microsystems Wetzlar GmbH, Wetzlar Germany) due to the reduced size of the layer.
Two thermocouple recording systems were placed in the middle of the plates, 2 mm distant from the weld centre line in both titanium and aluminum side. The calibration of the two models was carried out by comparing the size and shape of the fusion zone of the numerical model and the experimental one. The validation of the model was made by comparing the experimental and numerical thermal cycle.

3. Numerical Model

3.1. Model for the Plates

The plates of titanium and aluminum (200 mm × 50 mm) were joined along the long side. The plates have different thickness, i.e., 2 mm for titanium one and 3 mm for Aluminum one. The adopted mesh is the same for the two models. To obtain accurate results, a fine mesh was adopted close to the welding line. Mesh size was determined rapidly by trial-and-error; likewise, it is generally done in the literature. Mesh size is equal to 1 × 0.5 × 0.5 mm3 inside of 10 mm from the interface and the number of the nodes is equal to 20 along the x-axis. At a distance from the interface higher than 10 mm, the mesh size increases along the x-axis. Particularly from 10 to 50 mm far from the interface (the limit of the Al sheet), the number of the nodes is equal to 20, but the ratio in size between the last element and first element in the distribution is 0.1 (arithmetic sequence). The mapped mesh has 40,000 elements and 49,692 nodes. Figure 2 shows the results of the meshing procedure.
The numerical simulation was performed using the finite element code COMSOL Multiphysics. Comsol Multiphysics is a Multiphysics modeling tool that solves all types of problems based on Finite Element [15].
During the welding process, the transfer of heat is governed by the general equation heat flow (1a):
ρ C p δ T δ t   +   ρ C p v T   +   q   =   h f
where
ρ   [ g / mm 2 ] is the density of the metal as a function of temperature,
C p [J/(g·K)] is specific heat of the metal as a function of temperature,
v [m/s] is the velocity field,
q   [ W / mm 2 ] is the heat lost to the surroundings by combination of radiation and convection and conduction,
h f   [ W / mm 2 ] is the heat source.
The heat lost is given by Equation (1b):
q = ε σ ( T r 4 T 4 ) + h ( T r T )
where ε is the emissivity of the surface and is taken as 0.5 for titanium and 0.3 for aluminum. σ is the Stefan–Boltzmann constant and is taken as 5.67 × 10−8 W/(m2·K4), Tr is the room temperature and was taken at 293 K. h is the heat transfer coefficient assumed equal to 20 W/(m2·K) for air and 200 W/(m2·K) for the bottom surface in contact with the workspace.
When the contact is formed by pressing two similar or dissimilar metallic materials together, only a small fraction of the nominal surface area is in contact because of the roughness and irregularities of the contacting surfaces. When a heat flux is imposed across the junction, there are only a limited number and size of the contact spots that results in an actual contact area. The actual contact area is significantly smaller than the apparent contact area and causes a thermal contact resistance. There are several analytic expressions for predicting the contact conductance, and several values in literature [16]. Based on these values, by acting on a contact resistance value, the calibration has been carried out by comparing the fusion zones of the numerical model with the experimental results. Finally, the value of contact conductance is assumed to be equal to 30,000 W/(m2·K), for both of the two models [17]. Mechanical constraints were imposed for the simulation of the clamping system. Constraints were imposed to the four exterior nodes of the plates so degrees of freedom were zero (no displacement is permitted).

3.2. 2D Heat Source

Since the pioneering work of Rosenthal [6] that proposed punctual and linear heat sources, several more realistic sources have been proposed. When the distribution along the thickness is not important like in thin plates, the surface Gaussian heat source model is a good proposal for bed-on-plate cases when both TIG and conductive laser welding must be simulated [18]. For the surface laser heat source, the radius ( R ) should be calculated first by the formula written as [19]:
R   =   2 M 0 2 λ f π D 0 ,
where M 0 2 is the beam quality equal to 1.1 for fiber laser, f is the focal length of the focusing lens and D 0 is the diameter of the lens.
For laser welding, the heat source of the laser beam was simulated by a traveling two-dimensional distribution of heat source (Figure 3a). Particularly, the heat surface distribution was built by combining two Gaussian distributions one in the plane XZ and one moving with welding rate v in plane YZ. Those two Gaussian distributions were obtained, and, with an energy of about 99.7% of the total laser energy, a fusion of the titanium alloys was obtained. The heat source radius ( R ) was assumed equal to three standard deviations of the Gaussian pulse (Figure 3b):
g p ( s ) = 1 σ 2 π e t 2 2 σ 2 = 3 R 2 π e 9 s 2 2 R 2
The 2D total heat source (W) (Figure 4) is expressed in Equation (4):
h s 2 D ( x , y ) 2 D = P laser [ g p ( x ) × g p ( y v t ) ]

3.3. 3D Heat Source

For the three-dimensions laser source, the heat flux was described in the numerical model as a volume and surface distributed heat flux:
h f 3 D = ( 1 φ ) [ h f sur ( x , y ) ] load + φ [ h f vol ( z ) ] load
where the coefficient φ is the energy fraction that will be introduced through the cylinder, and the remaining energy will be introduced via the surface heat source (the value being used is 0.9).
For the surface, a Gaussian distribution was used, the heat source radius was equal to two standard deviations of the Gaussian pulse under the assumption that 95.44% of the total fusion energy of the titanium alloys was applied. In this case, considering the surface heat distribution proposed by Goldak and Akhlagi [20], the surface heat flux h f sur [W/mm2] can be expressed as shown in Equation (6):
h f sur ( x , y ) = 2 η P laser π R 2 e 2 r 2 R 2
where η is the process efficiency (the value being used is 1), P laser is the laser power in W and R is the heat source radius.
For volume heat source, a constant heat distribution cylinder was considering, with a radius R FZ equal to those of the molten cylinder (Figure 5). The numeric value for volumetric heat source is given by the ratio between the laser power and the volume of cylinder with a radius R FZ .
Thus, finally, the keyhole is modelled in the software using the heat source radius, the radius of the molten cylinder, measured by experimental results, and the energy fraction.

4. Results and Discussion

4.1. Metallurgical Characterization of Weld

The appearance of Ti6Al4V/AA5754 laser welded joint after chemical etching is shown in Figure 6. Both titanium and aluminum alloy melted at the joint interface and separate fusion zones were observed (Al side and Ti side) as well as two heat affected zone (HAZ) between FZ and the base material. Good weld appearance of the cross section was obtained with full penetration and low level of porosity. It was demonstrated that in aluminum-titanium dissimilar weld, porosity tended to be produced in the fusion line. It happens that gases in the seam are hard to escape and concentrated in the middle of fusion line during the solidification process [21]. In Figure 6, some gas trapping occurred on both sides of the intermetallic layer but not in the intermetallic layer.
Figure 7 shows the Ti microstructure of the different zones showed in Figure 6, i.e., base material, heat affected zone and fusion zone. The basic mental (BM) is composed of dark β phase in the dominating bright α matrix. Particularly, the β phase is distributed at the boundary of the α grains. This is a typical microstructure for α-β titanium alloys in mill-annealed conditions [22]. The microstructure within the joint depends on the heat received from the laser beam, and varies according to the distance from it. In FZ, a predominantly martensitic microstructure of acicular type (α’) is present. In fact, as reported in literature [22,23,24], the microstructure of the laser fusion zone of Ti-6Al-4V alloy is completely martensitic due to the high cooling speed from β field. The heat affected zone is a mixture of martensitic and primary α grain. During the heat thermal cycle due to the laser heating, an increase of β phase results. However, due to the rapid cooling, that transformation is never complete, so the microstructure of the HAZ is mixed, formed from martensitic grains and grain α.
The base material aluminum alloy AA5754 was supplied in the annealed condition (Figure 7). In the aluminum matrix, there are some second phases. These phases have been identified in literature as (Fe,Mn)Al6, (Fe,Mn)3SiAl12, Mg2Si and Mg2Al3 [25,26,27].
The microstructure of the fusion zone exhibits a very fine dendritic microstructure due to the high cooling rate. The dendrites grow in direction parallel to the direction of heat flow giving to columnar grains [28,29]. In the HAZ microstructure, no detectable changes are evident to the OM; however, in this zone, solubilization of magnesium based has been reported in the literature [30]. Figure 8 shows the micro hardness profile (0.3/15) in the transverse section of the weld. The microhardness was very high in the titanium FZ where the microstructure was martensitic. In the HAZ, the value diminished with the lower amount of martensitic microstructure. The rise in the microhardness of the FZ of aluminum in the weld was caused by the rapid cooling that produced a very fine solidification structure and solid solution strength. The increment in the HAZ of aluminum was due to the dissolution of magnesium compounds during the welding cycle [30].
A layer of intermetallic compounds (IMC) formed between the two fusion zones (Figure 9), whose stoichiometry was clarified in some paper [26,31,32]. This layer is due to the reaction in the temperature of the two alloys, and the size is variable as a function of the process parameters [32]. Particularly for the welding process parameters concerning this study, the average thickness of the IMC layer was equal 50 ± 5 µm to and the average nanohardness value (0.01/15) was equal to 2485 ± 232.

4.2. Calibration of the Model

After the simulation, the fusion zone was compared with the experimental one for calibration purposes. Figure 10 shows the comparison between the simulated and experimental fusion zone profile for aluminum and titanium. For both 2D and 3D heat source modeling, a complete penetration of keyhole was obtained, and, particularly for the 2D heat source, the boundary of the fusion zone obtained numerically had similar shapes in comparison with the experimental ones.
In Table 5 and Table 6, the dimensions of the fusion zone profile, respectively, in the top, middle and bottom, are reported both for 2D and 3D numerical simulation and the experimental measure. A better matching between the numerical and experimental results for the size of fusion zone has been obtained using a 2D source.
In fact, by correctly adjusting keyhole parameter in 3D heat source, only the top width of the weld was close to the experimental result, while the 2D heat source could almost precisely calculate the top, middle and bottom width. The observed performance can be explained as follows. It is well known that the behavior of heat flow depends on welding conditions [14,33]. For higher welding powers and thinner plates, the heat flow is predominantly 2D, whereas for lower welding powers and thicker plates, the heat flow was 3D. Therefore, there must be a transition thickness in which the heat flow has a behavior, which is a mix of the 2D and 3D ones. Probably for the thin plates examined in this study, the heat flow was predominantly 2D, and that is why 2D heat flow modelling is more appropriate.
Figure 11 and Figure 12 show the temperature distribution on the top when the weld bead profile reached the quasi-steady state. The thermal conductivity of the two materials influenced the position of the maximum temperature in the melt pool and away from the melt pool. For both numeric simulation, the maximum temperature was recorded within the titanium side, where the beam laser was directed. The width of the heat affected zone in the aluminum side was greater due to the high thermal conductivity of the aluminum alloy. On the contrary, the low thermal conductivity of the titanium alloy leads to accumulating the heat over the metal by reducing the area interested from metallurgical transformation.
The shape of the weld pool longitudinal plane (Figure 13) showed a teardrop shape for the 2D heat source. Thus, this means that probably the 2D heat source cannot perfectly simulate the laser heat source. However, better results have been obtained using the 3D heat source; in this case, the resulting molten puddle had an elliptical shape, which is more similar to the experimental pool shape as shown by microstructure evolution of the grains (Figure 14) in the fusion zone. This shape of the grains derives from the elliptical molten pool [24,34].

4.3. Validation of the Model

The temperature cycle was measured in two different points that were placed 2 mm away from the weld centerline on both sides of the weld (Figure 15 and Figure 16). The temperature profiles were compared with those obtained from the thermocouples. High thermal gradients, and fast cooling rate were present, due to the laser process. When it concerns the temperature peak, the numerical results matched the experimental one for both of the heat sources. Nonetheless, regarding the cooling rate, the 3D heat source was more precise than the 2D one.

5. Conclusions

In this paper, the numerical models to simulate the laser welding process of butt dissimilar Al/TI joints were developed. Two different 2D and 3D heat source modelling processes have been utilized to simulate the proper heat flux during the welding. The numeric results were compared with the experimental ones to calibrate and to validate the two models. The metallurgical analyses showed that the titanium fusion zone was principally martensitic, and the heat affected zone was a mixture of martensitic and primary α grain. In the aluminum fusion zone, a dendritic structure was present and the heat affected zone was characterized by the solubilization of the magnesium compounds.
The FEM simulation of the thermal cycle of fiber offset welding was satisfactory. The following points were demonstrated:
(1)
The calculations for the fusion zone dimensions were accurate both for the 2D and the 3D heat source. By using that 2D heat source, a better matching of numeric and experimental results was obtained at the three levels at which the molten zone sizes were taken.
(2)
In the longitudinal section, the numerical results were not as accurate for both of the heat sources. For the 2D one, a teardrop shape of the molten weld pool formed while the 3D heat source produced an elliptical one. It is possible to conclude that the 3D heat source can better approximate the heat flux during laser welding and the maximum temperature gradients, which determined the change in the grain growth direction in the titanium side.
(3)
The overall thermal cycle accuracy was good for 2D and 3D heat sources. However, the 3D heat source provided better results for the cooling rate simulation.

Author Contributions

The welds were fabricated in the TISMA ((Innovative Welding for Advanced Materials) laboratory of Bari, which is led by Giuseppe Casalino. Paola Leo looked after the metallographic preparation and analysis of the microstructure and Sonia D’Ostuni built the numerical model. Discussion and conclusions were written with the contribution of all authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oliveira, J.P.; Zeng, Z.; Andrei, C.; Braz Fernandes, F.M.; Miranda, R.M.; Ramirez, A.J.; Omori, T.; Zhou, N. Dissimilar laser welding of superelastic NiTi and CuAlMn shape memory alloys. Mater. Des. 2017, 128, 166–175. [Google Scholar] [CrossRef]
  2. Casalino, G.; Guglielmi, P.; Lorusso, V.D.; Mortello, M.; Peyre, P.; Sorgente, D. Laser offset welding of AZ31B magnesium alloy to 316 stainless steel. J. Mater. Process. Technol. 2017, 242, 49–59. [Google Scholar] [CrossRef]
  3. Zeng, X.; Oliveira, J.P.; Yang, M.; Song, D.; Peng, B. Functional fatigue behavior of NiTi-Cu dissimilar laser welds. Mater. Des. 2017, 114, 282–287. [Google Scholar] [CrossRef]
  4. Katayama, S. Handbook of Laser Welding Technologies; Woodhead Publishing Limited: Sawston, UK, 2013. [Google Scholar]
  5. Casalino, G.; Mortello, M.; Peyre, P. Yb–YAG laser offset welding of AA5754 and T40 butt joint. J. Mater. Process. Technol. 2015, 223, 139–149. [Google Scholar] [CrossRef]
  6. Rosenthal, D. Mathematical theory of heat distribution during welding and cutting. Weld. J. 1941, 20, 220–234. [Google Scholar]
  7. Zeng, Z.; Li, X.; Miao, Y.; Wu, G.; Zhao, Z. Numerical and experiment analysis of residual stress on magnesium alloy and steel butt joint by hybrid laser-TIG welding. Comput. Mater. Sci. 2011, 50, 1763–1769. [Google Scholar] [CrossRef]
  8. Pham, S.M.; Tran, V.P. Study on the Structure Deformation in the Process of Gas Metal Arc Welding (GMAW). Am. J. Mech. Eng. 2014, 2, 120–124. [Google Scholar]
  9. Dhinakaran, V.; Suraj Khope, N.S.S.; Sankaranarayanasamy, K. Numerical Prediction of Weld Bead Geometry in Plasma Arc Welding of Titanium Sheets Using COMSOL. In Proceedings of the 2014 COMSOL Conference in Bangalore, Bangalore, India, 13–14 October 2014. [Google Scholar]
  10. Casalino, G.; Hu, S.J.; Hou, W. Deformation prediction and quality evaluation of the gas metal arc welding butt weld. Proc. Inst. Mech. Eng. Part B 2003, 217, 1615–1622. [Google Scholar] [CrossRef]
  11. Zeng, Z.; Wang, L.; Wang, Y.; Zhang, H. Numerical and experimental investigation on temperature distribution of the discontinuous welding. Comput. Mater. Sci. 2009, 44, 1153–1162. [Google Scholar] [CrossRef]
  12. Casalino, G.; Michelangelo, M. A FEM model to study the fiber laser welding of Ti6Al4V thin sheets. Int. J. Adv. Manuf. Technol. 2016, 86, 1339–1346. [Google Scholar] [CrossRef]
  13. Azizpour, M.; Ghoreishi, M.; Khorram, A. Numerical simulation of laser beam welding of Ti6Al4V sheet. J. Comput. Appl. Res. Mech. Eng. 2015, 4, 145–154. [Google Scholar]
  14. Falk, N.; Flaviu, S.; Benjamin, K.; Jean, P.B.; Jorg, H. Optimization strategies for laser welding high alloy steel sheets. Phys. Proced. 2014, 56, 1242–1251. [Google Scholar]
  15. Comsol Multiphysics, version 5.2a; Software for Multiphysics Modeling; COMSOL: Stockholm, Sweden, 2016.
  16. Contuzzi, N.; Campanelli, S.L.; Casalino, G.; Ludovico, A.D. On the role of the Thermal Contact Conductance during the Friction Stir Welding of an AA5754-H111 butt joint. Appl. Therm. Eng. 2016, 104, 263–273. [Google Scholar] [CrossRef]
  17. Casalino, G.; Mortello, M.; Peyre, P. FEM Analysis of Fiber Laser Welding of Titanium and Aluminum. Procedia CIRP 2016, 41, 992–997. [Google Scholar] [CrossRef]
  18. Teixeira, P.R.D.F.; Araújo, D.B.D.; Cunda, L.A.B.D. Study of the gaussian distribution heat source model applied to numerical thermal simulations of tig welding processes. Ciênc. Eng. 2014, 23, 115–122. [Google Scholar] [CrossRef]
  19. Steen, W.M.; Mazumder, J. Laser Material Processing; Springer: Berlin, Germany, 2010. [Google Scholar]
  20. Goldak, J.; Akhlagi, M. Computational Welding Mechanics; Springer: Ottawa, ON, Canada, 2005; pp. 26–27. [Google Scholar]
  21. Lv, S.X.; Jing, X.J.; Huang, Y.X.; Xu, Y.Q.; Zheng, C.Q.; Yang, S.Q. Investigation on TIG arc welding–brazing of Ti/Al dissimilar alloys with Al based fillers. Sci. Technol. Weld. Join. 2012, 17, 519–524. [Google Scholar] [CrossRef]
  22. Lutjering, G.; Williams, J.C. Titanium (Engineering Materials and Processes), 2nd ed.; Springer: Berlin, Germany; New York, NY, USA, 2007. [Google Scholar]
  23. Akman, E.; Demir, A.; Canel, T.; Sınmazcelik, T. Laser welding of Ti6Al4V titanium alloys. J. Mater. Process. Technol. 2009, 209, 3705–3713. [Google Scholar] [CrossRef]
  24. Xu, P.; Li, L.; Zhang, C. Microstructure characterization of laser welded Ti-6Al-4V fusion zones. Mater. Charact. 2014, 87, 179–185. [Google Scholar] [CrossRef]
  25. Kumar, S.; Nadendla, H.B.; Scamans, G.M.; Eskin, D.G.; Fan, Z. Solidification behavior of an AA5754 alloy ingot cast with high impurity content. Int. J. Mater. Res. 2012, 103, E1–E7. [Google Scholar] [CrossRef]
  26. Raghavan, V. Phase Diagram Evaluations: Section II. J. Ph. Equilib. Diffus. 2005, 26, 171–172. [Google Scholar] [CrossRef]
  27. ASM Specialty Handbook: Aluminum and Aluminum Alloys; Davis, J.R. (Ed.) ASM International: Novelty, OH, USA, 1993. [Google Scholar]
  28. Porter, D.A.; Easterling, K.E. Phase Transformations in Metals and Alloys; Chapman & Hall: London, UK, 1992. [Google Scholar]
  29. Verhoeven, J.D. Fundamentals of Physical Metallurgy; Wiley: Hoboken, NJ, USA, 1975. [Google Scholar]
  30. Leo, P.; D’Ostuni, S.; Casalino, G. Hybrid welding of AA5754 annealed alloy: Role of post weld heat treatment on microstructure and mechanical properties. Mater. Des. 2016, 90, 777–786. [Google Scholar] [CrossRef]
  31. Bailey, N. Welding Dissimilar Metals; The Welding Institute: Cambridge, UK, 1986. [Google Scholar]
  32. Gao, M.; Chen, C.; Gu, Y.; Zeng, X. Microstructure and tensile behavior of laser arc hybrid welded dissimilar Al and Ti alloys. Materials 2014, 7, 1590–1602. [Google Scholar] [CrossRef]
  33. Sorensen, M.B. Simulation of Welding Distortion in Ship Section. Ph.D. Thesis, University of Denmark, Lyngby, Denmark, 1999. [Google Scholar]
  34. Messler, R.W. Principles of Welding; Wiley-VCH, John Wiley & Sons, Inc.: Singapore, 2004. [Google Scholar]
Figure 1. Scheme of laser welding on AA5754/Ti6Al4V.
Figure 1. Scheme of laser welding on AA5754/Ti6Al4V.
Metals 07 00307 g001
Figure 2. Mesh outlook.
Figure 2. Mesh outlook.
Metals 07 00307 g002
Figure 3. (a) Combining of the two Gaussian distribution in plane XZ and XY, and (b) Gaussian distribution.
Figure 3. (a) Combining of the two Gaussian distribution in plane XZ and XY, and (b) Gaussian distribution.
Metals 07 00307 g003
Figure 4. Heat flux in 2D Gaussian heat distribution.
Figure 4. Heat flux in 2D Gaussian heat distribution.
Metals 07 00307 g004
Figure 5. Total heat flux in 3D heat distribution.
Figure 5. Total heat flux in 3D heat distribution.
Metals 07 00307 g005
Figure 6. Appearance of Ti6Al4V/AA5754 laser welded joint after chemical etching.
Figure 6. Appearance of Ti6Al4V/AA5754 laser welded joint after chemical etching.
Metals 07 00307 g006
Figure 7. Microstructure of the various zones in Ti6Al4V and AA5754.
Figure 7. Microstructure of the various zones in Ti6Al4V and AA5754.
Metals 07 00307 g007
Figure 8. Microhardness profile at half thickness of the weld cross section.
Figure 8. Microhardness profile at half thickness of the weld cross section.
Metals 07 00307 g008
Figure 9. SEM Micrograph of intermetallic compounds (IMC) layer at the Al/Ti joint interface.
Figure 9. SEM Micrograph of intermetallic compounds (IMC) layer at the Al/Ti joint interface.
Metals 07 00307 g009
Figure 10. Calibration of (a) 3D and (b) 2D heat source models by the cross section of the weld.
Figure 10. Calibration of (a) 3D and (b) 2D heat source models by the cross section of the weld.
Metals 07 00307 g010
Figure 11. Isometric view of temperature distributions (in Kelvin) using 2D heat source modelling.
Figure 11. Isometric view of temperature distributions (in Kelvin) using 2D heat source modelling.
Metals 07 00307 g011
Figure 12. Isometric view of temperature distributions (in Kelvin) using 3D heat source modelling.
Figure 12. Isometric view of temperature distributions (in Kelvin) using 3D heat source modelling.
Metals 07 00307 g012
Figure 13. Weld pool numerical results for the numerical simulations.
Figure 13. Weld pool numerical results for the numerical simulations.
Metals 07 00307 g013
Figure 14. Interface zoom-up in the weld cross section ns.
Figure 14. Interface zoom-up in the weld cross section ns.
Metals 07 00307 g014
Figure 15. Experimental and numerical thermal cycle using 2D heat source modelling.
Figure 15. Experimental and numerical thermal cycle using 2D heat source modelling.
Metals 07 00307 g015
Figure 16. Experimental and numerical thermal cycle using 3D heat source modelling.
Figure 16. Experimental and numerical thermal cycle using 3D heat source modelling.
Metals 07 00307 g016
Table 1. Chemical composition of AA5754 aluminum and Ti6Al4V titanium alloys (wt %).
Table 1. Chemical composition of AA5754 aluminum and Ti6Al4V titanium alloys (wt %).
AA5754
SiFeCuMnMgCrZnTiAl
0.400.400.100.502.6–3.60.300.20<0.15balance
Ti6Al4V
CFeN2O2AlVH2TiC
<0.08<0.25<0.05<0.25.53.5<0.0375balance<0.08
Table 2. Mechanical and thermo-physical properties of the two alloys.
Table 2. Mechanical and thermo-physical properties of the two alloys.
PropertyAA5754Ti6Al4V
Young modulus [GPa]70114
Poisson ratio0.30.3
Density [g/cm3]2.74.4
Liquidus Temperature [K]8701923
Solidus Temperature [K]8561880
Table 3. Temperature dependent thermal conductivity of AA5754.
Table 3. Temperature dependent thermal conductivity of AA5754.
AA5754
TemperatureThermal Conductivity [W/mK]
293138
373147.2
473152.7
573162.7
673152.7
773158.75
873138
1773138
Table 4. Temperature dependent thermal conductivity of Ti6Al4V.
Table 4. Temperature dependent thermal conductivity of Ti6Al4V.
Ti6Al4V
TemperatureThermal Conductivity [W/mK]
2936.01
773.1514.78
793.1515
823.1515.15
953.1517.20
993.1517.80
1013.1518.30
1053.1518.80
1093.1519.50
1113.1520
1133.1520.50
1153.1521
1173.1521.60
1273.1523.91
1933.1534.3
Table 5. Dimension (mm) of the fusion zone profile in the aluminum side.
Table 5. Dimension (mm) of the fusion zone profile in the aluminum side.
Aluminum Fusion Zone2D Heat Source3D Heat SourceExperimental Data
Top118136116
Middle120135112
Bottom114135108
Table 6. Dimension (mm) of the fusion zone profile in the titanium side.
Table 6. Dimension (mm) of the fusion zone profile in the titanium side.
Titanium Fusion Zone2D Heat Source3D Heat SourceExperimental Data
Top232204225
Middle207196198
Bottom226187196

Share and Cite

MDPI and ACS Style

D’Ostuni, S.; Leo, P.; Casalino, G. FEM Simulation of Dissimilar Aluminum Titanium Fiber Laser Welding Using 2D and 3D Gaussian Heat Sources. Metals 2017, 7, 307. https://doi.org/10.3390/met7080307

AMA Style

D’Ostuni S, Leo P, Casalino G. FEM Simulation of Dissimilar Aluminum Titanium Fiber Laser Welding Using 2D and 3D Gaussian Heat Sources. Metals. 2017; 7(8):307. https://doi.org/10.3390/met7080307

Chicago/Turabian Style

D’Ostuni, Sonia, Paola Leo, and Giuseppe Casalino. 2017. "FEM Simulation of Dissimilar Aluminum Titanium Fiber Laser Welding Using 2D and 3D Gaussian Heat Sources" Metals 7, no. 8: 307. https://doi.org/10.3390/met7080307

APA Style

D’Ostuni, S., Leo, P., & Casalino, G. (2017). FEM Simulation of Dissimilar Aluminum Titanium Fiber Laser Welding Using 2D and 3D Gaussian Heat Sources. Metals, 7(8), 307. https://doi.org/10.3390/met7080307

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop