Next Article in Journal
Advancements in Nickel-Phosphate/Boron Based Electroless Composite Coatings: A Comprehensive Review of Mechanical Properties and Recent Developments
Next Article in Special Issue
Multi-Objective Optimization of The Low-Pressure Casting of Large-Size Aluminum Alloy Wheels through a Systematic Optimization Idea
Previous Article in Journal
A Green Synthesis Strategy for Cobalt Phosphide Deposited on N, P Co-Doped Graphene for Efficient Hydrogen Evolution
Previous Article in Special Issue
Investigation of Pre-Aged Hardening Single-Point Incremental Forming Process and Mechanical Properties of AA6061 Aluminum Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formation and Evolution Mechanism of Intermetallic Compounds of Friction Stir Lap Welded Steel/Aluminum Joints

1
Hubei Longzhong Laboratory, Wuhan University of Technology, Xiangyang 441000, China
2
Hubei Key Laboratory of Advanced Technology for Automotive Components, Wuhan University of Technology, Wuhan 430070, China
3
Hubei Engineering Research Center for Green & Precision Material Forming, Wuhan University of Technology, Wuhan 430070, China
4
Department of Mechanical and Electrical Engineering, Wuhan Donghu University, Wuhan 430212, China
5
Shanghai Aerospace Equipments Manufacturer Co., Ltd., Shanghai 200245, China
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(18), 6118; https://doi.org/10.3390/ma16186118
Submission received: 1 August 2023 / Revised: 25 August 2023 / Accepted: 31 August 2023 / Published: 7 September 2023

Abstract

:
Interfacial layers with brittle intermetallic compounds (IMC) greatly influence the performance of steel–aluminum friction stir lap welding (FSLW) joints. Thus, the formation and evolution of IMC between 7075-T6 aluminum alloy and galvanized DP590 steel in steel–aluminum FSLW joints were investigated. An FSLW numerical model was developed using the computational fluid dynamics method to analyze the interface temperature between the aluminum alloy and steel. Scanning electron microscopy (SEM) was conducted to observe the microstructure characterization and measure the IMC thickness. Phases among different joint zones were analyzed by X-ray diffraction (XRD) and energy dispersive spectroscopy (EDS). IMC layer formation was predicted by the effective Gibbs free energy model presented in this paper according to thermodynamic principles. The Monte Carlo method was utilized to predict the thickness of IMC layers. It was found that the IMC layer at the interface of the welded joint is composed of Fe2Al5, FeAl3, and Al-Zn eutectic. The IMC thickness decreased from 4.3 μm to 0.8 μm with the increasing welding speed, which was consistent with the Monte Carlo simulation results.

1. Introduction

The automobile lightweight is one of the hot spots in the automobile manufacturing industry. The steel–aluminum integrated body frame structure realizes the lightweight and high strength of the body frame structure by using aluminum alloy material to replace steel components in the traditional body frame steel structure [1]. The steel–aluminum composite structure can not only consider the lightweight and cost requirements but also give play to the strength advantages of high-strength steel and the weight reduction and energy absorption advantages of aluminum alloy [2]. Although there are many methods used to connect aluminum alloy to steel, including MIG arc brazing [3], laser welding [4], self-stamping riveting [5], and other technologies, to further improve the reliability of aluminum steel joints, the solid-state welding method of friction stir welding (FSW) is preferred [6].
During the welding process, solid phase diffusion of Fe and Al atoms occurs between base materials under the welding thermal cycle and severe plastic deformation. An intermetallic compound (IMC) is formed at the joint overlap interface, achieving metallurgical bonding between steel and aluminum alloy. It is clear that IMC has a significant role in achieving steel/aluminum metallurgical bonding and affecting joint performance [7,8,9]. The intervention of the Zn coat promoted the formation of an Al–Zn eutectic structure with a low melting point at the interface, which significantly improved the weldability of aluminum alloy and steel [10]. Furthermore, Zn coatings on steel were found to be beneficial for the formation of intermetallic reaction zones in solid/liquid joining procedures [11].
Both IMC composition and thickness have a significant impact on metallurgical bonding strength, which in turn determines joint performance [6]. According to the content of the element, IMC can be divided into Fe-rich IMC (FeAl, Fe3Al) and Al-rich IMC (FeAl2, Fe2Al5, and FeAl3). There are significant differences in performance between the two different types of IMCs, with Fe-rich IMC exhibiting better toughness, while Al-rich IMC is a brittle hard phase. Therefore, Fe-rich IMCs are more conducive to improving the metallurgical bonding strength of joints [8]. Due to significant differences in thermal cycling and plastic deformation experienced by different regions of the overlap interface during the welding process, there are usually two or more different components of IMC at the interface [6,8], and the IMC composition is influenced by welding temperature and element concentration. At present, some studies [12,13] find that the hard and brittle nature of IMC has an adverse impact on joint performance, and the thickness of IMC in steel/aluminum joints should be minimized as much as possible. However, some studies [14,15,16] find that there is a critical value for the thickness of IMC in steel/aluminum joints. When the thickness of IMC is thick, due to its hard and brittle properties, cracks are prone to occur under residual welding stress, which has a negative impact on the performance of the joint. However, when the thickness of IMC is below the critical value, the presence of IMC will not harm the strength of the joint and can even improve the performance of the joint. However, there are relatively few studies on the formation and evolution mechanism of IMC at the interface of steel/aluminum FSLW joints.
This article focuses on the unclear formation and evolution mechanism of IMC in the FSLW process of the 7075-T6 aluminum alloy and galvanized DP590 steel. The FSLW numerical simulation model based on computational fluid dynamics was established to clarify the temperature distribution of the joint. The microstructure of the joint interface was analyzed through OM and SEM. XRD and EDS were utilized to figure out the composition of IMC. Based on the temperature distribution of the joint during the welding process and the changes in Gibbs free energy, the formation order and evolution of IMCs were studied to clarify the formation mechanism of intermetallic compounds, while the IMC thickness was simulated by the Monte Carlo method.

2. Materials and Methods

2.1. Materials and Experiment

The base materials used in the present study were a 3 mm thick 7075-T6 aluminum alloy sheet and a 1.2 mm thick galvanized DP590 steel sheet. Figure 1 shows a schematic illustration of the plate arrangement in the lap configuration before welding. A welding tool (made of H13 die steel) with a smooth pin was adopted in the experiment. The shoulder diameter of the tool is 12 mm, and the conical pin has a length of 2.9 mm. The root and tip diameter of the pin is 5 mm and 3.6 mm. The tilt angle of the tool was set at 2.5 degrees. As shown in Figure 2, the zinc coating in steel sheet is uniform and continuous with average thickness of 5.5 μm.
After welding, the joint was cross-sectioned perpendicular to the welding direction. The microstructure of the cross-sectioned joint was analyzed using optical microscopy (Zeiss Scope A1, Zeiss, Oberkochen, Germany) and scanning electron microscopy (SEM, JEM-7500F, JEOL, Tokyo, Japan). Energy dispersive spectrometer (EDS) analysis was carried out to determine the element distribution at the joint interfaces. The phase composition of the specimen was analyzed by X-ray diffraction (XRD, Empyrean). The cross-sections of the metallographic specimens were polished and etched with Keller’s reagent (2.5 mL HF, 1 mL HCl, 1.5 mL HNO3, and 95 mL H2O) and 4% nitric acid alcohol for the aluminum alloy and steel, respectively.

2.2. Simulations

2.2.1. Friction Stir Lap Welding Simulation Model

An FSLW numerical model was developed using the computational fluid dynamics method in ANSYS/FLUENT to analyze the temperature distribution under the steady-state welding of the 7075-T6 aluminum alloy and DP590 steel under the conditions listed in Table 1. Only the overlap area of workpieces was set as computational domain, which is illustrated in Figure 3. The swirling process of stirring pin was not taken into consideration in the model.
To model the dissimilar FSLW of 7075 aluminum alloy and DP590 steel, the material is treated as an incompressible non-Newtonian fluid, with a viscoplastic flow [7]. The material flow constitutive equation is expressed by the hyperbolic sinusoidal relation proposed by Sellar and Tegart [17].
ε ˙ ¯ = A ( sinh α σ ) n exp ( Q R T ) ,
where ε ˙ ¯ is the effective strain rate, σ is the flow stress, Q is the thermal activation energy of deformation, R is the gas constant, T is the temperature, n is the stress index, A and α are constants related to the material. The material constitutive parameters are listed in Table 2 [18,19,20].
During the FSLW process, the heat generation of the contact surface between the workpiece and the tool involves friction heat generation and plastic deformation heat generation [18]. In this model, both heat inputs from friction and plastic deformation were taken into consideration. The heat generation rates at the contact surface between the shoulder/pin and the workpiece are defined by Equations (2)–(4) [7,21], respectively.
q s b ( r ) = [ η ( 1 δ ) τ c + δ μ f P N ] ( ω r u w e l d sin θ ) ,
q p b ( r ) = [ η ( 1 δ ) τ c + δ μ f P N ] ( ω r u w e l d sin θ ) ,
q p s ( r ) = [ η ( 1 δ ) τ c + δ μ f σ y ] ( ω r u w e l d sin θ ) ,
where q s b ( r ) , q p b ( r ) , and q p s ( r ) represent the heat generation rates of the shoulder bottom, the pin bottom, and the pin side surface; η is the conversion efficiency of plastic deformation; δ is the slip coefficient of the interface; τ c is the shear yield stress of plastic material;   μ f is the friction coefficient of the stirring pin and the workpiece contact surface;     P N is the axial pressure of the stirring pin, which is equal to 45 MPa; σ y is the yield strength of the workpiece material; ω   is the welding speed; r represents the distance from the elemental area to the tool axis; u w e l d is the material flow speed; θ is the angle between the welding direction and the r radial vector direction. δ A l   = 0.4, δ S t e e l = 0.7, μ f A l   = 0.35, μ f s t e e l = 0.3.
In addition, the plastic deformation heat source in the shear zone was also taken into consideration, which can be expressed by [7,21]
S v = f m μ Φ ,
where f m is an arbitrary constant that represents the extent of atomic mixing in the system, and Φ is given by [7,18]
Φ = 2 [ ( u x ) 2 + ( v y ) 2 + ( w z ) 2 ] + ( u y + v x ) 2 + ( u z + w x ) 2 + ( w y + v z ) 2 ,
where u, v, and w are the velocity components of material flow in the x-, y-, and z-directions, respectively. The relative velocity at the contact area between the shoulder/pin and the workpiece can be expressed by Equations (7)–(9) [22].
u = ( 1 δ ) ( ω r sin θ u w e l d ) ,
v = ( 1 δ ) ω r cos θ ,
w = ψ ω 2 π ,
where ψ is the thread pitch.
The convection heat transfer coefficient on the side and bottom surfaces of the workpiece that contacts the backing plate and clamping equipment was defined as 500.0 W/(m2·K) [7], the top surface of the workpiece was taken as a free surface contacting with air, and the convection heat transfer coefficient on the top surface was defined as 30.0 W/(m2·K) [7].

2.2.2. Monte Carlo Simulation Model

The Monte Carlo simulation model was developed in MATLAB code. The microstructure is generated based on a two-dimensional N × N array of a uniform-sized square cell (lattice site). N is the number of the discretized lattice points in two directions. To describe each lattice site, parameter q, which represents the crystallographic orientation of the lattice site (ranging from 1 to Q), is used. The initial lattice properties are set so that the lattice site belongs to Al (q = 0) and otherwise belongs to Fe (q = 1).
In this paper, the chemical-free energy item and the boundary energy term are introduced into Hamiltonian of the system in the q-state Potts model as follows:
H = 1 2 i n j m G S i S j b + i n G i c ,
where G i c is the chemical-free energy of cell i with its orientation Si; G S i S j b is the boundary energy between MC cell i and cell j with their corresponding orientations.
The boundary energy of the grain boundary between two neighboring cells (cell i and cell j) can be expressed with the Read–Shockley model as follows [23]:
γ = γ m ( θ θ m ) ( 1 ln θ θ m ) ,
where γ is the interface energy, θ is the orientation angle of grain boundary dislocation; γ m = 0.56 J m 2 , θ m = 15 [24]. When the grain boundary becomes a high-angle grain boundary, the energy of the grain boundary has nothing to do with the direction and angle of the dislocation.
In each Monte Carlo step (MCS), cell i is selected at random, and then neighbors of the selected lattice site and their crystallographic orientations were determined. One neighbor site was chosen at random out of the eight neighbors. If there is a certain proportion of Fe and Al in the nine lattices centered on i, there is a certain probability that a phase will form. If the selected lattice site (i) and chosen neighboring site (j) belong to two different phases, the chosen neighboring site orientation (qj) will be exchanged with the first selected lattice site orientation (qi). Calculate the energies of the selected lattice and chosen neighboring sites (H1) at the initial state, using Equation (10). Calculate the energies of the selected lattice and chosen neighboring sites (H2) after reorientation or orientation exchanges, using Equation (10). Calculate the net energy change (ΔH = H2H1) before and after reorientation or orientation exchanges. The new orientation assignment or exchange is accepted with a transition probability, Px, of [25,26]:
P x = { 1 Δ H < 0 exp ( Δ H k B T ) Δ H > 0 ,
where T is the simulation temperature, kB is the Boltzmann’s constant. Generate a random probability R between 0 and 1. The reorientation is accepted if RPx, otherwise the lattice site orientation is unchanged.

3. Results and Discussion

3.1. Temperature Distribution and Microstructure Morphology

Figure 4 shows the simulation results of temperature distribution at transverse cross-sections with different welding speeds. The maximum temperature point appears at the contact surface between the stirring pin and the steel. The flow stress of steel is greater than aluminum alloy, which accounts for higher heat production from plastic deformation in the welding process.
Figure 5 shows the macroscopic morphology of a typical steel–aluminum FSLW joint under the parameter of 1200 r/min–100 mm/min. The joint cross-section can be divided into the stirring zone (SZ), interface zone (IZ), thermo-mechanical affected zone (TMAZ), heat affected zone (HAZ), and base material zone (BZ). Under the mechanical and thermal effects, severe plastic deformation is produced in the SZ. Meanwhile, the material in the SZ reaches a good plastic flow state, and the material within a certain range near the stirring pin is stirred, resulting in a significantly wider SZ than the diameter of the stirring pin. During the welding process, the TMAZ does not come into direct contact with the welding tool, resulting in a low degree of plastic deformation. The TMAZ shows a significant difference between the advancing side and the retreating side of the joint, and the boundary between the advancing side and the SZ is more obvious. This is because, during the welding process, the material shear rate on the advancing side is relatively high. The structure of the HAZ only undergoes welding thermal cycles without any mechanical effect, and the area is also not affected by material plastic flow.

3.2. Microstructure Characterization

According to the XRD results in Figure 6, the IMCs of the galvanized steel–Al joint were determined to be Fe2Al5, FeAl3, and the Al-Zn eutectics. Figure 7 shows the IMC in different areas of a typical joint. During the steel/aluminum FSLW process, the Al and Fe atoms in the bonding area of the lap interface undergo solid-phase diffusion, and a continuum Al-Fe IMC layer was formed due to welding heat and plastic deformation. The results of the EDS point analysis are shown in Table 3. During the welding process, the Zn coating at the interface of the IZ is extruded by the stirring pin. The interface IMC in the IZ is generated by the solid phase diffusion of Al and Fe atoms. Its composition is only Fe2Al5 with a thickness of 0.8 µm (Figure 7a). The welding pressure on the interface of the TMAZ is relatively small, and there are some remaining Al-Zn eutectic components at the interface of this area, and the presence of Zn can promote the growth of Al-Fe IMC, so the IMC thickness of the TMAZ is 30~40 μm (Figure 7b). From the results of EDS point scanning, there are three components of Fe2Al5, FeAl3, and Al-Zn eutectic at the same time in IMC. The IMC distribution at the interface of the HAZ of the joint is shown in Figure 7c. Compared with TMAZ, the IMC thickness is reduced to only 15 μm. The results show that the IMC in the HAZ is divided into two layers, the upper layer is Al-Zn eutectic composition, and the lower layer is mainly the mixture of FeAl3 and Al-Zn eutectic composition.

3.3. Interfacial Layer Formation and Evolution

According to the effective heat of formation (EHF) model proposed by Pretorius et al. [25], the effective Gibbs free energy of formation can be obtained by analogy. The formation of primary phases of the Al-Fe binary system on the reaction and the order of formation of each phase can be predicted.
Δ G ( T ) = c e f f c 0 Δ G m ( T ) ,
where ∆G′(T) is the change in effective Gibbs free energy of formation at different temperatures, ∆Gm (T) is Moore Gibbs free energy change at different temperatures, ceff is the effective concentration of the limiting element on the reaction interface, c0 is the concentration of the limiting element in the compound. Since the effective concentration of elements at the reaction interface is difficult to obtain, the eutectic component at the lowest eutectic temperature in the reaction system is usually selected as the effective concentration [26]. For the Al-Fe binary system, according to the Al-Fe binary phase diagram, the effective concentration of Al Fe at the lowest eutectic temperature (928 K) is 99.1% and 0.9%, respectively.
The molar Gibbs free energy of the Al-Fe system can be expressed by [27]:
G m = G m s r f + G m c f g + G m E = i x i G i ( T ) + R T i x i ln ( x i ) + x i x j v L i j ( x i x j ) v
where G m s r f , G m c f g and G m E are the surface of reference, configurational, and the excess contributions, respectively; xi is the molar fraction of component i ; Gi( T ) represents the Gibbs energy of component i, which is a function of the temperature; R is the gas constant; T is the absolute temperature; Lij is the interaction parameter between components i   and j ; and the Lij and Gi( T ) values of Fe and Al are listed in Table 4 [7,27].
For the reaction xFe + yAl = Fe x Al y in the Al–Fe system, the Gibbs free energy change can be defined as:
Δ G = G m F e x A l y x G F e y G A l ,
According to the Al–Fe phase diagram [27], there are five possible components in the IMC layers: FeAl, Fe3Al, FeAl2, FeAl3, and Fe2Al5. Therefore, the Gibbs free energy change in the Al-Fe intermetallic compounds can be described by Equations (16)–(20), according to Equations (13)–(15):
Δ G m F e A l ( T ) = 33490 . 691 345 . 968 T + 55 . 262 TlnT 0 . 0141 × T 2 151451 × T 1
Δ G m F e 3 A l ( T ) = 59480.674 461 . 42 T + 76 . 72 TlnT 0 . 004 × T 2 229616.75 × T 1
Δ G m F e A l 2 ( T ) = 19813.351 380 . 803 T + 58 . 01 TlnT 0 . 0128 × T 2 150362 × T 1
Δ G m F e A l 3 ( T ) = 7552.39 593 . 22 T + 89 . 07 TlnT 0 . 0384 × T 2 224726.25 × T 1
Δ G m F e 2 A l 5 ( T ) = 6242.03 1167 . 51 T + 176 . 37 TlnT 0 . 072 × T 2 450152.57 × T 1
Figure 8 shows the ∆G′ of each Al-Fe IMC at the corresponding concentration of the lowest eutectic point at different temperatures. In the welding temperature range of 600 K~950 K, the ∆G′ of Fe2Al5 is the smallest, so Fe2Al5 is the first to be formed at the joint interface, while FeAl3 is formed after Fe2Al5. In the welding process, Fe2Al5 first nucleates and grows at the lap interface. When the continuous Fe2Al5 layer is formed, the diffusion of the Al and Fe atom is hindered, resulting in the decrease in the Fe atom concentration on the side of the Fe2Al5 layer away from the steel, and the growth of Fe2Al5 is inhibited. Meanwhile, the concentration of the Al atom on this side is relatively high, so FeAl3 is gradually formed on this side, and excessive Al can also react with Fe2Al5 to form FeAl3 [28].
In the steel/aluminum FSLW process, the IMC formation and evolution mechanism of the joint interface are shown in Figure 9. The formation process includes IMC nucleation and growth. First, as shown in Figure 9a, the temperature in the area in front of the stirring pin gradually rises under the effect of heat conduction, and the diffusion of various elements in this area gradually increases. When the temperature reaches 654 K, the zinc-rich coating on the surface of the steel plate reacts with Al atoms diffused to the interface to form a low melting point Al-Zn eutectic composition. The melted Al-Zn eutectic component is extruded by the stirring pin from the IZ and flows to other areas of the joint. While the welding pressure and temperature of the TMAZ and HAZ are lower than that of IZ, there remain some Al-Zn eutectic components, as shown in Figure 9b. As the welding process continues (Figure 9c), a hook-like structure will form at the interface. The aluminum alloy and the steel at the interface of the IZ are in direct contact with each other in a high-temperature and high-strain environment. When Al and Fe atoms start solid-phase diffusion, the element concentration at the joint interface gradually increases. With the distance increasing from the IZ, the plastic deformation and the temperature of the material gradually decrease, and the element diffusion effect is weakened. Therefore, the Al and Fe elements in the TMAZ and HAZ are not as concentrated as those in the IZ.
As the diffusion of atoms progresses (Figure 9d), a solid solution structure is formed at the interface of the joint. When the diffusion reaches a certain level, the solid solution at the interface transforms into a supersaturated state. In this state, a stable phase compound will be formed at the interface. In the Al-Fe system, the reaction always proceeds in the direction of the decrease in Gibbs free energy, that is, the lower the Gibbs free energy of the compound, the easier it is to form. Therefore, Fe2Al5 first nucleates at the interface between the IZ and TMAZ. At the interface of the HAZ, FeAl3 is the first to nucleate due to the limitation of element concentration. After the formation of the IMC crystal nucleus, the solute atoms in the solid solution at the interface precipitate out of the supersaturation state. When Al and Fe atoms continue to diffuse to the interface of the joint, the solid solution transforms into a supersaturated state again, causing the IMC nucleus to grow, and the IMC nucleus continues to be generated at the higher grain boundary energy, and finally, IMC grows along with the interface and forms a continuous IMC layer (Figure 9e).
The formation of a continuous IMC layer at the interface of the IZ blocks the mutual contact between the steel and aluminum alloy (Figure 9f), and the interface changes from the Al-Fe interface to the Al-Fe2Al5, Fe-Fe2Al5 interfaces. At this time, the growth of IMC is mainly affected by the influence of the diffusion rate of Al and Fe atoms. Since the diffusion rate of Fe atoms in the Fe2Al5 layer is much lower than that of Al atoms, the growth of Fe2Al5 mainly occurs at the Fe-Fe2Al5 interface. In the TMAZ, the presence of the Fe2Al5 layer hinders the diffusion of Fe atoms to the aluminum side, resulting in a lower concentration of Fe atoms in the Fe2Al5 layer near the aluminum side, so that the Al-rich phase FeAl3 is formed on this side and is mixed with the Al-Zn eutectic component. Finally, a dense Fe2Al5 layer is formed at the interface of the IZ of the joint, and the IMC of the TMAZ is divided into two layers, namely the Fe2Al5 layer and a mixed layer of FeAl3 + Al-Zn eutectic composition, and the IMC of the HAZ is divided into Al -Zn eutectic and FeAl3 + Al-Zn eutectic composition mixed layer.

3.4. IMC Thickness

In the FSLW process, the growth of IMC at the lap interface is closely related to the welding parameters. The formation and growth of IMC are determined by the solid phase diffusion of Fe and Al atoms. The change in welding parameters can cause changes in the welding heat input, and cooling rate, which will cause changes in the diffusion of Al and Fe atoms, and ultimately have an impact on the growth of IMC. Therefore, combining the atomic diffusion process and IMC growth kinetics, the growth of the IMC layer at the joint interface is analyzed. Although there are different regions in the welding joint, the complexity of IMC composition at the interface of HAZ and TMAZ leads to uncertainty in the growth of IMC, and the IZ plays an important role in mechanical properties such as tensile strength and shear strength. Therefore, only the growth process of the IMC layer in the bonding interface area of the IZ is simulated.
According to the mechanism of IMC formation, the formation of IMC includes nucleation and growth. In the nucleation stage, the reaction rate between Al and Fe atoms is the main factor limiting the nucleation of IMC. In the planar growth process of IMC, the diffusion of Al and Fe atoms is inhibited by the Fe2Al5 layer. The growth of IMC is limited by the diffusion rate of Al and Fe atoms.
Figure 10 shows the interface morphology of numerical simulation and the experimental and simulated IMC thickness changes under different welding parameters. The experimental results can be viewed in our previous work [2]. With the increase in welding speed, the welding heat input decreases, the diffusion between atoms weakens, and the growth of IMC is inhibited, resulting in a decrease in IMC thickness. There is a certain difference between the simulated values and the experimental values. In this model, only the grain boundary energy and Gibbs free energy are considered, while the deformation energy is not taken into consideration. Meanwhile, the atomic diffusion rate is set as the function of the distance from the initial steel–aluminum interface fails to accurately consider the actual physical diffusion effect.

4. Conclusions

(1)
In the process of FSLW of steel and aluminum alloy, due to the friction heat and plastic deformation, the composition of the IMC layer in the lap interface of the joint is complex, including Fe2Al5, FeAl3, and Al-Zn eutectic.
(2)
The effective Gibbs free energy change indicates that Fe2Al5 is the first phase formed in the Al-Fe system in the welding temperature range of the joint interface between 600 K and 950 K.
(3)
During the welding process, when the continuous Fe2Al5 layer is formed at the IZ, the diffusion of Al and Fe atoms is hindered, so the growth of IMC in the IZ is inhibited, and the thickness of IMC is relatively thin. However, due to the weak plastic deformation of the matrix in the TMAZ and the existence of the Fe2Al5 layer, the diffusion of the Fe element is greatly inhibited, resulting in the formation and growth of the FeAl3 layer near the side of the Fe2Al5 layer to the aluminum material. In addition, the presence of Zn elements promotes the growth of the IMC in the TMAZ. Due to the low welding temperature and no plastic deformation in the HAZ, resulting in weak element diffusion, the overlapping interface only forms FeAl3 and a hybrid layer with Al-Zn eutectic composition.
(4)
The simulation results of the IMC thickness in the IZ are relatively close to the experimental results. The Monte Carlo model can predict the IMC thickness in the IZ within a certain parameter range.

Author Contributions

Y.L.: software; data curation; formal analysis; investigation; methodology; visualization; writing—original draft. Z.Z.: formal analysis; investigation. L.H.: validation; resources. Q.P.: conceptualization; methodology; funding acquisition; project administration; resources; supervision; validation; writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (52075400, 52275368); Independent Innovation Projects of the Hubei Longzhong Laboratory (2022ZZ-04); the Key Research and Development Program of Hubei Province (2022AAA001); the Science and Technology Innovation Talent Program of Hubei Province (2022EJD012). This study is also sponsored by the Shanghai Rising-Star Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Long, J.Q.; Lan, F.C.; Chen, J.Q. New technology of lightweight and steel-aluminum hybrid structure car body. J. Mech. Eng. 2008, 44, 27–35. [Google Scholar] [CrossRef]
  2. Yu, H.Y.; Hu, Z.L. Investigation of interfacial microstructure and mechanical properties for DP590 steel/7075 aluminum dissimilar materials friction stir lap welding joints. J. Mech. Eng. 2020, 56, 65–72. [Google Scholar] [CrossRef]
  3. Kumar, T.; Kiran, D.V.; Cheon, J.; Arora, N. Probing the Influence of Cold Wire Gas Tungsten Arc Welding Current Waveforms on the Aluminum-Steel Joining. J. Manuf. Process. 2020, 59, 378–388. [Google Scholar] [CrossRef]
  4. Wallerstein, D.; Salminen, A.; Lusquiños, F.; Comesaña, R.; García, J.D.V.; Rodríguez, A.R.; Badaoui, A.; Pou, J. Recent Developments in Laser Welding of Aluminum Alloys to Steel. Metals 2021, 11, 622. [Google Scholar] [CrossRef]
  5. Abe, Y.; Mori, K. Mechanical Clinching and Self-Pierce Riveting for Sheet Combination of 780-MPa High-Strength Steel and Aluminium Alloy A5052 Sheets and Durability on Salt Spray Test of Joints. Int. J. Adv. Manuf. Technol. 2021, 113, 59–72. [Google Scholar] [CrossRef]
  6. Wan, L.; Huang, Y. Friction Stir Welding of Dissimilar Aluminum Alloys and Steels: A Review. Int. J. Adv. Manuf. Technol. 2018, 99, 1781–1811. [Google Scholar] [CrossRef]
  7. Hu, Z.; Yu, H.; Pang, Q. Investigation of Interfacial Layer for Friction Stir Welded AA7075-T6 Aluminum to DP1180 Steel Joints. J. Manuf. Sci. Eng. 2020, 142, 091002. [Google Scholar] [CrossRef]
  8. Pourali, M.; Abdollah-zadeh, A.; Saeid, T.; Kargar, F. Influence of Welding Parameters on Intermetallic Compounds Formation in Dissimilar Steel/Aluminum Friction Stir Welds. J. Alloys Compd. 2017, 715, 1–8. [Google Scholar] [CrossRef]
  9. Xiang, Q.; Lv, N.C.; Xue, P.; Ma, Z.Y. Research status and prospect on friction stir welded dissimilar Al-Steel joints. J. Mech. Eng. 2017, 5, 28–37. [Google Scholar] [CrossRef]
  10. Chen, Y.C.; Komazaki, T.; Tsumura, T.; Nakata, K. Role of Zinc Coat in Friction Stir Lap Welding Al and Zinc Coated Steel. Mater. Sci. Technol. 2008, 24, 33–39. [Google Scholar] [CrossRef]
  11. Springer, H.; Szczepaniak, A.; Raabe, D. On the Role of Zinc on the Formation and Growth of Intermetallic Phases during Interdiffusion between Steel and Aluminium Alloys. Acta Mater. 2015, 96, 203–211. [Google Scholar] [CrossRef]
  12. Tanaka, T.; Morishige, T.; Hirata, T. Comprehensive Analysis of Joint Strength for Dissimilar Friction Stir Welds of Mild Steel to Aluminum Alloys. Scr. Mater. 2009, 61, 756–759. [Google Scholar] [CrossRef]
  13. Elrefaey, A.; Gouda, M.; Takahashi, M.; Ikeuchi, K. Characterization of Aluminum/Steel Lap Joint by Friction Stir Welding. J. Mater. Eng. Perform. 2005, 14, 10–17. [Google Scholar] [CrossRef]
  14. Movahedi, M.; Kokabi, A.H.; Reihani, S.M.S.; Najafi, H. Effect of Tool Travel and Rotation Speeds on Weld Zone Defects and Joint Strength of Aluminium Steel Lap Joints Made by Friction Stir Welding. Sci. Technol. Weld. Join. 2012, 17, 162–167. [Google Scholar] [CrossRef]
  15. Liu, X.; Lan, S.; Ni, J. Analysis of Process Parameters Effects on Friction Stir Welding of Dissimilar Aluminum Alloy to Advanced High Strength Steel. Mater. Des. 2014, 59, 50–62. [Google Scholar] [CrossRef]
  16. Bozzi, S.; Helbert-Etter, A.L.; Baudin, T.; Criqui, B.; Kerbiguet, J.G. Intermetallic Compounds in Al 6016/IF-Steel Friction Stir Spot Welds. Mater. Sci. Eng. A 2010, 527, 4505–4509. [Google Scholar] [CrossRef]
  17. Shi, L.; Wu, C.S.; Gao, S.; Padhy, G.K. Modified Constitutive Equation for Use in Modeling the Ultrasonic Vibration Enhanced Friction Stir Welding Process. Scr. Mater. 2016, 119, 21–26. [Google Scholar] [CrossRef]
  18. Sheppard, T.; Jackson, A. Constitutive Equations for Use in Prediction of Flow Stress during Extrusion of Aluminium Alloys. Mater. Sci. Technol. 1997, 13, 203–209. [Google Scholar] [CrossRef]
  19. Mahmoudiniya, M.; Kokabi, A.H.; Kheirandish, S.; Kestens, L.A.I. Microstructure and Mechanical Properties of Friction Stir Welded Ferrite-Martensite DP700 Steel. Mater. Sci. Eng. A 2018, 737, 213–222. [Google Scholar] [CrossRef]
  20. Yang, Y.; Qian, H.; Su, Y. Effect of Mn Addition on Deformation Behaviour of 23% Cr Low Nickel Duplex Stainless Steel. Mater. Charact. 2018, 145, 606–618. [Google Scholar] [CrossRef]
  21. Shi, L.; Wu, C.S.; Padhy, G.K.; Gao, S. Numerical Simulation of Ultrasonic Field and Its Acoustoplastic Influence on Friction Stir Welding. Mater. Des. 2016, 104, 102–115. [Google Scholar] [CrossRef]
  22. Sharghi, E.; Farzadi, A. Simulation of Strain Rate, Material Flow, and Nugget Shape during Dissimilar Friction Stir Welding of AA6061 Aluminum Alloy and Al-Mg2Si Composite. J. Alloys Compd. 2018, 748, 953–960. [Google Scholar] [CrossRef]
  23. Read, W.T.; Shockley, W. Dislocation Models of Crystal Grain Boundaries. Phys. Rev. 1950, 78, 275–289. [Google Scholar] [CrossRef]
  24. Xiao, N.; Tong, M.; Lan, Y.; Li, D.; Li, Y. Coupled Simulation of the Influence of Austenite Deformation on the Subsequent Isothermal Austenite–Ferrite Transformation. Acta Mater. 2006, 54, 1265–1278. [Google Scholar] [CrossRef]
  25. Mohammadi, H.; Eivani, A.R.; Seyedein, S.H.; Ghosh, M. Modified Monte Carlo Approach for Simulation of Grain Growth and Ostwald Ripening in Two-Phase Zn–22Al Alloy. J. Mater. Res. Technol. 2020, 9, 9620–9631. [Google Scholar] [CrossRef]
  26. Tan, Y.; Maniatty, A.M.; Zheng, C.; Wen, J.T. Monte Carlo Grain Growth Modeling with Local Temperature Gradients. Model. Simul. Mater. Sci. Eng. 2017, 25, 065003. [Google Scholar] [CrossRef]
  27. Sundman, B.; Ohnuma, I.; Dupin, N.; Kattner, U.R.; Fries, S.G. An Assessment of the Entire Al–Fe System Including D03 Ordering. Acta Mater. 2009, 57, 2896–2908. [Google Scholar] [CrossRef]
  28. Bouayad, A.; Gerometta, C.; Belkebir, A.; Ambari, A. Kinetic Interactions between Solid Iron and Molten Aluminium. Mater. Sci. Eng. A 2003, 363, 53–61. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the FSLW for dissimilar joining between the Al alloy and steel.
Figure 1. Schematic representation of the FSLW for dissimilar joining between the Al alloy and steel.
Materials 16 06118 g001
Figure 2. The cross-section view of DP590 galvanized steel.
Figure 2. The cross-section view of DP590 galvanized steel.
Materials 16 06118 g002
Figure 3. The computational domain and its dimensions.
Figure 3. The computational domain and its dimensions.
Materials 16 06118 g003
Figure 4. Temperature distribution at transverse cross−section with different welding speeds.
Figure 4. Temperature distribution at transverse cross−section with different welding speeds.
Materials 16 06118 g004
Figure 5. Morphology of steel–aluminum FSLW joint.
Figure 5. Morphology of steel–aluminum FSLW joint.
Materials 16 06118 g005
Figure 6. XRD results of the joint interface.
Figure 6. XRD results of the joint interface.
Materials 16 06118 g006
Figure 7. IMC distribution in different regions (a) IZ, (b) TMAZ, (c) HAZ.
Figure 7. IMC distribution in different regions (a) IZ, (b) TMAZ, (c) HAZ.
Materials 16 06118 g007
Figure 8. Values of ∆G′ at different temperatures.
Figure 8. Values of ∆G′ at different temperatures.
Materials 16 06118 g008
Figure 9. IMC formation and evolution mechanism. (af) describes the IMC formation and evolution order through the FSLW process.
Figure 9. IMC formation and evolution mechanism. (af) describes the IMC formation and evolution order through the FSLW process.
Materials 16 06118 g009
Figure 10. Simulation results of IMC under different welding parameters (ac) and IMC thickness between experiment and simulation (d).
Figure 10. Simulation results of IMC under different welding parameters (ac) and IMC thickness between experiment and simulation (d).
Materials 16 06118 g010
Table 1. The welding parameters for FSLW.
Table 1. The welding parameters for FSLW.
Rotational Speed (rpm)Welding Speed (mm/min)Depth of Penetration (mm)
120025/50/1001.2
Table 2. Material constitutive parameters [18,19,20].
Table 2. Material constitutive parameters [18,19,20].
MaterialsαAnQ
7075-T6 1.03   × 1090.01415.41129,000
DP590 3.612   × 10130.007556.897329,000
Table 3. EDS results of each point in Figure 7 (at%).
Table 3. EDS results of each point in Figure 7 (at%).
PointsAlFeZnValue
171.5028.170.33Fe2Al5
272.9526.110.94Fe2Al5
370.3628.531.11Fe2Al5
468.5823.068.36FeAl3 + Al-Zn
575.2924.550.16FeAl3
685.320.9113.77Al-Zn
767.4524.867.69FeAl3 + Al-Zn
Table 4. Computational parameters in the Gibbs free energy change [7,27].
Table 4. Computational parameters in the Gibbs free energy change [7,27].
SymbolFormula/Value
G A l 12061.144 + 219.68 T 41.75 T ln T + 0.01853 × T 2 + 74092 T 1
G F e 22557.71 + 131.817 T 23.5143 T ln T 0.0043975 × T 2 + 77359 T 1
L F e , A l 0 91976.5 + 22.1314 T
L F e , A l 1 5672.58 + 4.8728 T
L F e , A l 2 121.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Pang, Q.; Zhang, Z.; Hu, L. Formation and Evolution Mechanism of Intermetallic Compounds of Friction Stir Lap Welded Steel/Aluminum Joints. Materials 2023, 16, 6118. https://doi.org/10.3390/ma16186118

AMA Style

Liu Y, Pang Q, Zhang Z, Hu L. Formation and Evolution Mechanism of Intermetallic Compounds of Friction Stir Lap Welded Steel/Aluminum Joints. Materials. 2023; 16(18):6118. https://doi.org/10.3390/ma16186118

Chicago/Turabian Style

Liu, Yongzhi, Qiu Pang, Zhichao Zhang, and Lan Hu. 2023. "Formation and Evolution Mechanism of Intermetallic Compounds of Friction Stir Lap Welded Steel/Aluminum Joints" Materials 16, no. 18: 6118. https://doi.org/10.3390/ma16186118

APA Style

Liu, Y., Pang, Q., Zhang, Z., & Hu, L. (2023). Formation and Evolution Mechanism of Intermetallic Compounds of Friction Stir Lap Welded Steel/Aluminum Joints. Materials, 16(18), 6118. https://doi.org/10.3390/ma16186118

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