Comprehensive Weldability Criterion for Magnetic Pulse Welding of Dissimilar Materials

: Despite its exceptional ability to join dissimilar materials and environmental friendliness, several challenges must be addressed in magnetic pulse welding (MPW). The conventional weldability criterion (i.e., minimum impact velocity) is analytically calculated as a function of material properties without considering the geometry of electromagnetic coil, electrical and physical parameters, making the minimum impact velocity a necessary but not sufﬁcient condition for a sound MPW joint. A new weldability criterion, namely effective impact velocity, is proposed, which overcomes the conventional weldability criterion’s limitations. The effective impact velocity can be inversely modelled to identify shop-ﬂoor relevant process parameters and it eliminates the need to fabricate several coils in the process and product proving stages. The proposed approach is demonstrated by a case study on tubular welding of Aluminium and SS304. The weld’s soundness produced with computed process parameters was corroborated by experimental observations on lap shear tests, hardness measurements, optical and scanning electron microscopy, and surface energy dispersive spectroscopy mapping. This investigation is expected to pave the way for developing the process window for MPW of several material combinations, with high cost and time savings.


Introduction
Developing optimal lightweight structures by joining dissimilar materials has increased remarkably in many industries, such as automotive, aerospace, tooling, power generation, and marine applications.The joining of dissimilar materials by fusion welding causes unfavorable microstructures [1].Chemical interaction between dissimilar materials leads to metallurgical incompatibility and brittle intermetallic compound (IMC) formation.Other issues include differences in the thermal and physical properties (thermal conductivity and coefficient of thermal expansion), the type of suitable heat treatment of the dissimilar joint, and galvanic corrosion [1].When joining dissimilar materials, the product and joining process design must overcome the above-stated challenges [2].Solidstate welding techniques facilitate joints' formation at low temperatures and often very quickly, usually within microseconds [3].Examples include explosive welding (EXW) [4], magnetic pulse welding (MPW) [1], vaporizing foil actuator welding [5], cold welding, and diffusion welding.These processes involve reduced formation of harmful and brittle IMCs, thereby keeping the material properties intact.Among the solid-state processes, MPW is one of the most environmentally friendly methods for joining dissimilar materials, wherein electromagnetic forces impact one metal onto another to form a solid-state cold weld.Am- pere's law governs the process.The force (F) experienced by two parallel current-carrying conductors with the distance of separation (d) is given as [6]: where I 1 and I 2 are the currents flowing in the conductors and µ 0 is the permeability of free space.This force is known as the Lorentz force and is repulsive when the currents flow in opposite directions.Figure 1 shows the setup of a tubular MPW configuration, including the outer flyer tube, the coil, and the field shaper, which drives the flyer tube to cause a high-speed impact with the inner target tube.The flyer and target tubes collide at velocities from 250 to 500 m/s [7], forming a solid-state joint because of severe plastic deformation.
Metals 2022, 12, x FOR PEER REVIEW 2 of 18 other to form a solid-state cold weld.Ampere's law governs the process.The force () experienced by two parallel current-carrying conductors with the distance of separation () is given as [6]: where  1 and  2 are the currents flowing in the conductors and  0 is the permeability of free space.This force is known as the Lorentz force and is repulsive when the currents flow in opposite directions.Figure 1 shows the setup of a tubular MPW configuration, including the outer flyer tube, the coil, and the field shaper, which drives the flyer tube to cause a high-speed impact with the inner target tube.The flyer and target tubes collide at velocities from 250 to 500 m/s [7], forming a solid-state joint because of severe plastic deformation.MPW might offer a joint strength around that of the weaker material without any emission, making the process environmentally friendly.Despite its many advantages, the spread of MPW is limited because of the higher initial investment in electromagnetic coils.Several coils are tested before obtaining a shop-floor applicable coil design.Numerical modeling allows the materials and geometry to be changed without additional expenditure and can be used to determine a process's feasibility with high precision.However, this depends upon how critically numerical results are analyzed and corroborated with experimental results.The literature on combined modeling and experimental work in MPW of dissimilar metals is presented in Table 1.Although considerable effort has been put forward to develop rigorous numerical models combined with experimental work that can accurately capture the MPW process, there is a lack of work that concentrates on developing weldability criteria.

Investigators
Material Pair Work Reported Hisashi et al. [8] Al-Fe Examined magnetic pressure seam welding through numerical simulations Zhang et al. [9] Al-Cu Determination of process parameters through experiments and numerical simulations Kore et al. [10] Al-SS Weld formation criteria for plate joint based on numerical simulations and experiments Shim et al. [11] Al-Cu Electromagnetic force estimation in dissimilar MPW using finite element method (FEM) Kakizaki et al. [12] Al-Cu, Al-Ni Numerical and experimental study on weld interface formation for dissimilar metal lap joints Zhidan et al. [13] Al-Fe Calculation of impact rate in the dissimilar tubular joint through numerical simu- MPW might offer a joint strength around that of the weaker material without any emission, making the process environmentally friendly.Despite its many advantages, the spread of MPW is limited because of the higher initial investment in electromagnetic coils.Several coils are tested before obtaining a shop-floor applicable coil design.Numerical modeling allows the materials and geometry to be changed without additional expenditure and can be used to determine a process's feasibility with high precision.However, this depends upon how critically numerical results are analyzed and corroborated with experimental results.The literature on combined modeling and experimental work in MPW of dissimilar metals is presented in Table 1.Although considerable effort has been put forward to develop rigorous numerical models combined with experimental work that can accurately capture the MPW process, there is a lack of work that concentrates on developing weldability criteria.The investigations listed in Table 1 consider minimum impact velocity, impact angle, and interface morphology as indicators of the success of the joint.Recently, an analytical model was developed to identify the parameter combinations leading to creating a wavy pattern [39].In their review, Ribeiro et al. [40] provided an integrated description of the weldability limits and criterion and its associated equations with the variables and parameters.The minimum impact velocity, also known as the threshold impact velocity criterion, is one-sided and can be considered necessary but not sufficient.The minimum impact velocity is calculated using material properties (as described in a later part).It does not consider the actual scenario, which includes the quality of mating surfaces, the geometry of an electromagnetic coil and electrical and physical parameters, such as air gap and tube/plate thickness.Therefore, the calculated minimum impact velocity might not always ensure a sound weld.The minimum impact velocity criterion does not prescribe any target velocity or corresponding process parameters.The investigation in this article proposes an approach to computationally identify target process parameters used on the shop floor after corroboration with a minimal number of experiments.This study's broader objective is to present an integrated numerical and experimental framework that could overcome technical and economic limitations in product and process development through MPW.This investigation proposes a comprehensive weldability criterion-effective impact velocity-that is sufficient to overcome surface imperfections and incorporates material properties and process parameters in the numerical algorithm used to compute the impact velocity.The investigation also demonstrates inverse modeling of effective impact velocity to obtain coil geometry and operating parameters required for shop-floor applications without fabricating several coils.The following section presents the FEM used in this investigation along with the details of the experimental techniques.Weldability criteria and methods for the selection of process parameters are presented after that.The efficacy of the proposed approach is then evaluated through the proximity of numerical and experimental observations on the interfacial mapping of hardness, plastic strain, and elemental distribution, following the analysis and discussion of the results obtained in this investigation.

Finite Element Model
MPW involves a complex combination of electromagnetism and structural mechanics.Maxwell's equations govern MPW.Substitution of the constitutive equations into Maxwell's equations leads to Equation (2) for the tube region: where magnetic permeability (H/m) is µ, Υ is medium conductivity (S/m), Υ (∂ → A/∂t) is the current density (A/m 2 ) and → A is the magnetic vector potential.The electromagnetic module calculates the transient magnetic forces at each time step.Maxwell's relation gives the magnetic force → f as follows: where → J is the current density in the electromagnetic coil (A/m 2 ) and → B is the magnetic flux density (T).The electromagnetic module outputs Lorentz force, an input for the structural module.The current passing through the coil is the body load for the electromagnetic module.
where charging voltage (kV The material conduct of the flyer and target tubes is described by the Cowper-Symonds constitutive model [41]: where quasi-static flow stress is σ y , plastic strain rate (s −1 ) is .ε, and m and P are constants that depend on the material.
A typical MPW simulation consists of an electromagnetic module (for the coil) and a structural module (for tubes).A nonlinear solver is used in the electromagnetic environment to calculate the transient magnetic forces.The forces are provided to the structural module to obtain deformation.The structural module considers the inertial effects due to timedependent stress.Figure 2 presents a flowchart of the sequentially coupled electromagneticstructural analysis.The FEM model was developed assuming that cracking, the heat generated by friction, deformation, and joule heating is negligible.Thus, temperatureindependent elastoplastic properties are used.The effect of resistance from air compression between the tubes is neglected.An Augmented Lagrangian contact pressure method and a time-dependent multifrontal massively parallel sparse direct solver (MUMPS) were employed.A field shaper has not been employed in this work.As explained next, a 2D axisymmetric problem using FEM was solved in COMSOL for the candidate material pair (Al-SS304).In this study, the geometry of the flyer was not changed during the simulations.Several process parameters that directly or indirectly influence the impact velocity, including input voltage, coil turns, coil length, coil cross-sectional area, capacitance, air gap, and current frequency, were varied for the simulations.Table 2 lists all the process parameters along with the range of variation.The operational details of the FEM are presented in Table 3.For additional details on the FEM employed in this study, the authors' previous work can be referred to [42,43].
Metals 2022, 12, x FOR PEER REVIEW 5 of 18 tural module to obtain deformation.The structural module considers the inertial effects due to time-dependent stress.Figure 2 presents a flowchart of the sequentially coupled electromagnetic-structural analysis.The FEM model was developed assuming that cracking, the heat generated by friction, deformation, and joule heating is negligible.Thus, temperature-independent elastoplastic properties are used.The effect of resistance from air compression between the tubes is neglected.An Augmented Lagrangian contact pressure method and a time-dependent multifrontal massively parallel sparse direct solver (MUMPS) were employed.A field shaper has not been employed in this work.As explained next, a 2D axisymmetric problem using FEM was solved in COMSOL for the candidate material pair (Al-SS304).In this study, the geometry of the flyer was not changed during the simulations.Several process parameters that directly or indirectly influence the impact velocity, including input voltage, coil turns, coil length, coil cross-sectional area, capacitance, air gap, and current frequency, were varied for the simulations.Table 2 lists all the process parameters along with the range of variation.The operational details of the FEM are presented in Table 3.For additional details on the FEM employed in this study, the authors' previous work can be referred to [42,43].

Experimental Method
A tubular bi-metallic joint configuration, as shown in Figure 3, was fabricated.The material composition of the constituent base materials and electromagnetic coil are given in Table 4.The feasible process parameters were obtained by a selection process described in the next section.
Metals 2022, 12, x FOR PEER REVIEW Event tolerance 0.01

Experimental Method
A tubular bi-metallic joint configuration, as shown in Figure 3, was fabricat material composition of the constituent base materials and electromagnetic coil ar in Table 4.The feasible process parameters were obtained by a selection proc scribed in the next section.The welded sample was tested using the lap-shear test, micro-hardness m ments, and metallographic examination using an Olympus STM6-LM optical mic (OM) (OLYMPUS, Tokyo, Japan) and Carl Zeiss AG-Supra 40 field electron sc electron microscope (FESEM) (ZEISS, Oberkochen, Germany).Elemental distr mapping was obtained through energy dispersive spectroscopy (EDS) analysis o and SS304 members' interface.Surface mapping of the interface was also carried check for the diffusion of one metal into the other at the welded interface.The welded sample was tested using the lap-shear test, micro-hardness measurements, and metallographic examination using an Olympus STM6-LM optical microscope (OM) (OLYMPUS, Tokyo, Japan) and Carl Zeiss AG-Supra 40 field electron scanning electron microscope (FESEM) (ZEISS, Oberkochen, Germany).Elemental distribution mapping was obtained through energy dispersive spectroscopy (EDS) analysis of the Al and SS304 members' interface.Surface mapping of the interface was also carried out to check for the diffusion of one metal into the other at the welded interface.Impact velocity is an essential criterion for successful bonding.The impact velocity relations available for EXW are equally applicable for MPW as both processes go through the same interface phenomenon [1].The critical impact pressure (P c ) essential for a successful joint is given as follows: where Z eq represents the equivalent acoustic impedance of the colliding sheets/tubes, Φ is the angle for the jet formation, and V T is the minimum (threshold) impact velocity essential for bond formation.The Z eq is computed by individual acoustic impedances of two members as follows: where Z 1 =ρ 1 S 1 is the flyer sheet/tube acoustic impedance, Z 2 =ρ 2 S 2 is the target sheet/tube acoustic impedance, S 1 and S 2 are the speeds of sound in the flyer and target sheet/tube materials, respectively, and ρ 1 and ρ 2 are the material densities of the two sheets/tubes.For a successful weld, the threshold value of pressure must exceed five times the Hugoniot elastic limit (HEL), which is given as: where K is the bulk modulus, G is the shear modulus, and Y 0 is the tensile yield stress.
The threshold value of impact velocity was calculated for the candidate material pair (Al-SS304) using the corresponding material properties presented in Table 2.The higher impact velocity among the two candidate materials was considered the threshold impact velocity, which in the present case was 161.03 m/s for SS304.
The threshold impact velocity (V T ) is a necessary but not sufficient condition for bond formation.As the interface of the weld is subjected to severe plastic deformation, an increase in impact velocity improves the likelihood of bonding.Moreover, achieving a theoretical threshold impact velocity might not remove the impurities and the oxide layer from the surface, leading to joint failure.On the other hand, excessive impact velocity damages the part itself.Even if the impact velocity remains under the maximum achievable impact velocity without damage (V D ), the possibility of IMC formation is increased at higher impact velocities.Thus, the effective impact velocity (V I ) (i.e., impact velocity expected during the actual experiment) should be higher than the threshold impact velocity but remain lower than V D .A typical jet behavior establishes a steady-state wave pattern at intermediate velocities, and the maximum bond strength is obtained [42].Because of these limitations, the effective impact velocity (V I ) is proposed as a weldability criterion.The desired (i.e., effective) impact velocity is considered to be in the vicinity of the average of the threshold impact velocity (V T ) and maximum achievable impact velocity without damage (V D ) as follows: The maximum achievable impact velocity without damage (V D ) can be estimated by FEM simulation, which is discussed later.However, FEM simulation overpredicts the impact velocity as the effect of resistance from air compression between the tubes is neglected [39].To compensate for this, a velocity correction factor (ε) is introduced, such that where V DS is simulated maximum achievable impact velocity without damage.Using Equation (10) in Equation ( 9) results in the following: The effective velocity can be used to obtain process parameters through inverse modeling.However, as explained earlier, the velocity correction factor (ε) is also applied to the simulated effective impact velocity.Thus, where V IS is simulated effective impact velocity.Applying Equation (12) in Equation ( 11) yields Metals 2022, 12, 1791 9 of 19 The maximum achievable impact velocity without damage (V DS ) is algorithmically obtained and used to calculate the effective impact velocity, which is then used to compute the parameters used on the shop floor.

Algorithm for Selection of Shop-Floor Applicable Parameters
Let V be the impact velocity influenced by m number of process parameters (X: x i ; i = 1 to m), each operating in a range of [x i min , x i max ].The maximum achievable impact velocity by altering x i, while keeping others at mid-range is V i max (i.e., V = V i max ; when x i = x i Vmax ).Let V r: v k be a set of five maximum possible impact velocities in the neighborhood (within 1% difference) that can be reached without damage and obtained through simulation.Once the maximum achievable impact velocity without damage (V DS ) is obtained, V IS (effective impact velocity) is calculated using Equation ( 13).The algorithm for calculating the effective impact velocity is shown in the flowchart in Figure 4.The process parameters corresponding to the effective impact velocity obtained by selective grid inverse modeling vis-à-vis the FEM are run several times by changing process parameters, such that every time the simulated impact velocity nears the effective impact velocity parameter, the matrix is updated.The process converges when no significant change in impact velocity is observed for a certain number of FEM runs.

Process Parameter Identification
The outcome of the process parameter selection exercise is shown in Figure 5.The one-parameter-at-a-time investigation indicates that for a candidate pair of a given geometry, higher voltage values, coil turns, coil cross-section area, capacitance and frequency, and a very narrow air gap can damage the flyer tube because of undesirably high impact velocity.The maximum achievable impact velocity without damage (VDs) and the simulated effective impact velocity (VIS) obtained using the algorithm shown in Figure 4 were 484 m/sand 338 m/s, respectively.The values of process parameters obtained from

Process Parameter Identification
The outcome of the process parameter selection exercise is shown in Figure 5.The oneparameter-at-a-time investigation indicates that for a candidate pair of a given geometry, higher voltage values, coil turns, coil cross-section area, capacitance and frequency, and a very narrow air gap can damage the flyer tube because of undesirably high impact velocity.The maximum achievable impact velocity without damage (V Ds ) and the simulated effective impact velocity (V IS ) obtained using the algorithm shown in Figure 4 were 484 m/s and 338 m/s, respectively.The values of process parameters obtained from the inverse modeling (voltage = 16 kV, air gap = 1.5 mm, coil turns = 6, coil length = 27 mm, coil cross-section area = 9 mm 2 (circular) and frequency = 50,000 rad/s) yielded a simulated impact velocity of 341 m/s, which is close to the calculated value.The welding sample was manufactured with the obtained process parameters and checked for different experimental attributes, as presented below.

Experimental Corroboration of FEM Observations
Experimental measurement of the impact velocity in MPW was largely conducted using optical methods, particularly with the use of a photonic Doppler velocimeter (PDV) [44,45].In their previous work, the authors validated the FEM employed in this study by comparing the simulated impact velocity with the experimental impact velocity measured using PDV for various process conditions.Figure A1a,b (see Appendix A) depicts the comparison of the simulated and experimental impact velocities.In this study, the corroboration of the FEM observations was achieved through a comparison of the interface morphology and mechanical behavior, as presented next.
(a) Interface morphology The pulse-welded Al-SS304 metal pair showed a wavy pattern at the impact zone, typical of MPW (Figure 6).The FEM also predicted a similar wavy morphology.

Experimental Corroboration of FEM Observations
Experimental measurement of the impact velocity in MPW was largely conducted using optical methods, particularly with the use of a photonic Doppler velocimeter (PDV) [44,45].In their previous work, the authors validated the FEM employed in this study by comparing the simulated impact velocity with the experimental impact velocity measured using PDV for various process conditions.Figure A1a,b (see Appendix A) depicts the comparison of the simulated and experimental impact velocities.In this study, the corroboration of the FEM observations was achieved through a comparison of the interface morphology and mechanical behavior, as presented next.
(a) Interface morphology The pulse-welded Al-SS304 metal pair showed a wavy pattern at the impact zone, typical of MPW (Figure 6).The FEM also predicted a similar wavy morphology.
ured using PDV for various process conditions.Figure A1a,b (see Appendix A) depicts the comparison of the simulated and experimental impact velocities.In this study, the corroboration of the FEM observations was achieved through a comparison of the interface morphology and mechanical behavior, as presented next.
(a) Interface morphology The pulse-welded Al-SS304 metal pair showed a wavy pattern at the impact zone, typical of MPW (Figure 6).The FEM also predicted a similar wavy morphology.The flyer tube started moving as soon as the magnetic pressure generated by the electromagnetic coil exceeded the plasticization pressure, attaining a maximum velocity at the time of impact, and decreased as the process progressed.The impact velocity variation across the weld length resulted in different interface phenomena.An MPW tubular joint welded in the longitudinal direction can be divided into three zones, as shown in the cut weld section in Figure 7a.The flyer tube started moving as soon as the magnetic pressure generated by the electromagnetic coil exceeded the plasticization pressure, attaining a maximum velocity at the time of impact, and decreased as the process progressed.The impact velocity variation across the weld length resulted in different interface phenomena.An MPW tubular joint welded in the longitudinal direction can be divided into three zones, as shown in the cut weld section in Figure 7a.The middle part of the mating members is where the joining occurs, and two non-welded zones occur on either side.The non-welded zones on the left and right in Figure 7a are called the run-in and run-out zones, respectively.The flyer tube and the target tube make an angle at the end of the run-out zone.A noticeable deformation occurs in the run-in zone.The target tube's deformation declines gradually from the run-in to the run-out zone (Figure 7a).In the unbonded zones, the flyer tube rebounds and cre- The middle part of the mating members is where the joining occurs, and two nonwelded zones occur on either side.The non-welded zones on the left and right in Figure 7a are called the run-in and run-out zones, respectively.The flyer tube and the target tube make an angle at the end of the run-out zone.A noticeable deformation occurs in the run-in zone.The target tube's deformation declines gradually from the run-in to the run-out zone (Figure 7a).In the unbonded zones, the flyer tube rebounds and creates a gap between the two plates.The distinct bonded and unbonded zones are visible in the numerical results (Figure 7b).The similarity between the experimental and numerical results is an important outcome of the proposed approach.
(b) Mechanical behavior Figure 8 shows the joints after the lap-shear test.For all specimens, fracture occurred outside the welded region and at weaker Al, proving that the weld was sound.This result validates the proposed approach, wherein process parameters obtained with numerical modeling passed the mechanical test.An agreement between experimentally observed hardness mapping and numerically obtained plastic strain also justified the proposed approach.The interface layer showed increased hardness relative to the base material (BM) (Figure 9a).The transition zone had the highest hardness value of 302 HV, much higher than the average BM hardness values of 55 HV for Al and 210 HV for SS304.The micro-hardness promptly increased on both sides of the transition zone and tended to be constant beyond these regions.This behavior was attributable to severe plastic deformation at the interface, observed in the simulated strain distribution in Figure 9b.Higher stress and strain were generated in the bonded center zone compared to the unbonded run-in and run-out zones.
(a) (b) With increasing distance from the interface, the plastic deformation decreased, and with it, the micro-hardness values.The interface layer went through two possible phenomena-solid-solid bonding caused by plastic deformation and liquid-solid bonding due to a very thin layer formed by molten metal.In addition, the high-hardness IMC produced in the transition zone might result in higher hardness than that of the BMs [46].An agreement between experimentally observed hardness mapping and numerically obtained plastic strain also justified the proposed approach.The interface layer showed increased hardness relative to the base material (BM) (Figure 9a).The transition zone had the highest hardness value of 302 HV, much higher than the average BM hardness values of 55 HV for Al and 210 HV for SS304.The micro-hardness promptly increased on both sides of the transition zone and tended to be constant beyond these regions.This behavior was attributable to severe plastic deformation at the interface, observed in the simulated strain distribution in Figure 9b.Higher stress and strain were generated in the bonded center zone compared to the unbonded run-in and run-out zones.An agreement between experimentally observed hardness mapping and numerically obtained plastic strain also justified the proposed approach.The interface layer showed increased hardness relative to the base material (BM) (Figure 9a).The transition zone had the highest hardness value of 302 HV, much higher than the average BM hardness values of 55 HV for Al and 210 HV for SS304.The micro-hardness promptly increased on both sides of the transition zone and tended to be constant beyond these regions.This behavior was attributable to severe plastic deformation at the interface, observed in the simulated strain distribution in Figure 9b.Higher stress and strain were generated in the bonded center zone compared to the unbonded run-in and run-out zones.
(a) (b) With increasing distance from the interface, the plastic deformation decreased, and with it, the micro-hardness values.The interface layer went through two possible phenomena-solid-solid bonding caused by plastic deformation and liquid-solid bonding due to a very thin layer formed by molten metal.In addition, the high-hardness IMC produced in the transition zone might result in higher hardness than that of the BMs [46].The metallurgical investigation in the next section discusses the same phenomenon.With increasing distance from the interface, the plastic deformation decreased, and with it, the micro-hardness values.The interface layer went through two possible phenomena-solid-solid bonding caused by plastic deformation and liquid-solid bonding due to a very thin layer formed by molten metal.In addition, the high-hardness IMC produced in the transition zone might result in higher hardness than that of the BMs [46].The metallurgical investigation in the next section discusses the same phenomenon.

Metallurgical Investigation
The transfer of material between the flyer and target materials resulted in several phenomena.One of the reasons for the increase in hardness (Figure 9a) in the interface layer or transition zone is the material transfer between the mating members due to the highspeed impact.The intermetallic phase formation at the welding interface cannot be avoided during MPW, particularly in the case of dissimilar material joining.Raoelison et al. [47,48] investigated the effect of process conditions on the interface properties and weld features for Al/Al and Al/Cu MPW joints.They found that an Al/Cu pair led to the formation of an intermetallic phase, whereas an Al/Al pair had metal continuity at the bonded interface.The interface wave results were due to flow velocity discontinuities [49].The discontinuities caused waves across the interface.The two fluids with different velocities resulted in instabilities at the interface, which led to mass flow from one material to the other.Thus, element diffusion is inevitable during MPW of dissimilar materials, as shown in the surface EDS map in Figure 10.The intermetallic zone formation could be a direct result of the atomic diffusion of one metal in the other.bonded interface.The interface wave results were due to flow velocity discontinuities [49].The discontinuities caused waves across the interface.The two fluids with different velocities resulted in instabilities at the interface, which led to mass flow from one material to the other.Thus, element diffusion is inevitable during MPW of dissimilar materials, as shown in the surface EDS map in Figure 10.The intermetallic zone formation could be a direct result of the atomic diffusion of one metal in the other.Line EDS is used to determine the distribution of Al and SS304 across the weld interface for the bonded center zone.Line analysis was performed across the intermediate layer (IML) at the crest of a wave, the trough of a wave, and a flat interface, as shown in Figure 11a-c, respectively.The area between the dotted lines corresponds to the IML region.The composition of elements (primarily aluminum) changed sharply in the IML region.Beyond this region, the composition of elements was almost constant.The width of the IML changed along the length of the bonded zone.Figure 11a shows the distribution of elements in the line scan performed along the crest of a wave.The width of the IML, in this case, was approximately 9 μm.The IML zone width increased slightly, by around 12 μm, at the trough of the wavy interface (Figure 11b).The distribution of Al changed gradually by several micrometers and then sharply decreased to zero.Along the flat interface (Figure 11c), there was a considerable decrease in the width of the IML (by around 5 μm).Line EDS is used to determine the distribution of Al and SS304 across the weld interface for the bonded center zone.Line analysis was performed across the intermediate layer (IML) at the crest of a wave, the trough of a wave, and a flat interface, as shown in Figure 11a-c, respectively.The area between the dotted lines corresponds to the IML region.The composition of elements (primarily aluminum) changed sharply in the IML region.Beyond this region, the composition of elements was almost constant.The width of the IML changed along the length of the bonded zone.Figure 11a shows the distribution of elements in the line scan performed along the crest of a wave.The width of the IML, in this case, was approximately 9 µm.The IML zone width increased slightly, by around 12 µm, at the trough of the wavy interface (Figure 11b).The distribution of Al changed gradually by several micrometers and then sharply decreased to zero.Along the flat interface (Figure 11c), there was a considerable decrease in the width of the IML (by around 5 µm).
tion of elements in the line scan performed along the crest of a wave.The width of the IML, in case, was approximately 9 μm.The IML zone width increased slightly, by around 12 μm, at the trough of the wavy interface (Figure 11b).The distribution of Al changed gradually by several micrometers and then sharply decreased to zero.Along the flat interface (Figure 11c), there was a considerable decrease in the width of the IML (by around 5 μm).The chances of intermetallic phase formation increased with the increase in discharge energies.The maximum thickness of the intermetallic phase can increase to more than 25 µ m at higher energies, which is detrimental to the weld quality and strength [50].However, the maximum thickness observed with the parameters obtained from numerical modeling was around 10-12 μm, meaning the interface was free of any detrimental effects that indicate the proposed approach's efficacy.The joint's strength is also a result of interlocking between the two metal surfaces, as shown in Figure 12.A high-speed collision leads to the formation of a discontinuous interface that is mostly wavy, where the vortices of waves take part in mechanical interlocking as a joining mechanism.The in- The chances of intermetallic phase formation increased with the increase in discharge energies.The maximum thickness of the intermetallic phase can increase to more than 25 µm at higher energies, which is detrimental to the weld quality and strength [50].However, the maximum thickness observed with the parameters obtained from numerical modeling was around 10-12 µm, meaning the interface was free of any detrimental effects that indicate the proposed approach's efficacy.The joint's strength is also a result of interlocking between the two metal surfaces, as shown in Figure 12.A high-speed collision leads to the formation of a discontinuous interface that is mostly wavy, where the vortices of waves take part in mechanical interlocking as a joining mechanism.The interlocking is promoted by Al's low strength and high ductility compared to SS304.This is like the combing action in dissimilar friction stir welding, where the low-strength alloy can penetrate high-strength alloy and create mechanical interlocking.However, it is to be noted that the notion of interlocking due to the swirling kinematic flow in MPW may also be a defect site based on the nature of the swirl-affected zone.This is more evident in the case of MPW of dissimilar materials, where an intense swirling motion along the bi-metallic interface can lead to the formation of the intermediate phase, severely degrading the weld properties [47,48].Thus, it is essential to have a proper set of input parameters that satisfy the weldability criterion and allow a significant reduction in the thickness of the intermediate phase.The experimental results of the present study are consistent with the FEM mode The FEM and experimental results show similar weld interface characteristics, whic gives confidence for predicting suitable process parameters for successful welding wit few experimental runs.The proposed approach to identify the effective impact velocit and corresponding process parameters through inverse modeling is the first of its kind.I has immense potential to reduce process development costs using the MPW process.Th developed model can be coupled with soft-computational tools such as artificial neura networks to enable machine learning-based evaluation of the MPW process.

Conclusions
This study investigates MPW and presents a new weldability criterion that over comes existing limitations and offers an economical and timely process developmen approach.The major conclusions from the present investigation are: 1.The conventional weldability criterion (threshold impact velocity) is necessary bu insufficient, as it is calculated using material properties without considering th geometry of the electromagnetic coil, electrical and physical parameters, such as ai gap and plate thickness, and surface imperfections.Such a criterion does not identif suitable process parameters for shop-floor applications.2. A new criterion (effective impact velocity) is proposed, in which the effective impac velocity is the average of the threshold and maximum possible velocity withou damage.The investigation offers a numeric algorithm to compute effective impac The experimental results of the present study are consistent with the FEM model.The FEM and experimental results show similar weld interface characteristics, which gives confidence for predicting suitable process parameters for successful welding with few experimental runs.The proposed approach to identify the effective impact velocity and corresponding process parameters through inverse modeling is the first of its kind.It has immense potential to reduce process development costs using the MPW process.The developed model can be coupled with soft-computational tools such as artificial neural networks to enable machine learning-based evaluation of the MPW process.

Conclusions
This study investigates MPW and presents a new weldability criterion that overcomes existing limitations and offers an economical and timely process development approach.The major conclusions from the present investigation are: 1.
The conventional weldability criterion (threshold impact velocity) is necessary but insufficient, as it is calculated using material properties without considering the geometry of the electromagnetic coil, electrical and physical parameters, such as air gap and plate thickness, and surface imperfections.Such a criterion does not identify suitable process parameters for shop-floor applications.

2.
A new criterion (effective impact velocity) is proposed, in which the effective impact velocity is the average of the threshold and maximum possible velocity without damage.The investigation offers a numeric algorithm to compute effective impact velocity, wherein the maximum possible velocity without damage is obtained through FEM simulation.

3.
The proposed weldability criterion overcomes the existing limitations and can be numerically computed and inversely modeled, and thus can be used to prescribe shop-floor applicable process parameters.4.
The weld samples obtained with numerically computed parameters were consistent with the experiments in terms of interface morphology, the intermediate layer's width, and plastic strain The joints cleared lap shear tests, wherein fracture occurred outside the welded region-increased hardness in and around the interface zone corresponded with the predicted plastic strain in the FEM simulation.Surface energy dispersive spectroscopy also showed that the increased hardness was associated with element transfer at the interface during the severe plastic deformation at the time of impact.

5.
The proposed approach is an experimentally corroborated and time-saving method that will encourage using finite element modeling to obtain shop-floor-applicable process parameters.Moreover, the cost will be reduced as use of the electromagnetic coil will be restricted to that of a few (essential for FEM validation), unlike conventional methods, wherein several coils are to be tested.
Figure A1a,b shows the comparison of experimental and simulated values of workpiece velocity computed by the FEM for dissimilar metal joining.It is seen that there was a close agreement between the simulated and experimental values within a range of ±10% variation.Figure A1a shows the comparison of simulated and experimentally measured impact velocities with variation in input energy and capacitance, and Figure A1b shows the comparison of simulated and experimentally measured impact velocities with variation in input voltage and air gap.[51], and (b) with variation of input voltage and air gap, data from [13,52].

Figure 2 .
Figure 2. Simulation flowchart for the developed FEM.

Figure 2 .
Figure 2. Simulation flowchart for the developed FEM.

Figure 3 .
Figure 3. Arrangement of flyer and target tubes.

Figure 3 .
Figure 3. Arrangement of flyer and target tubes.

Figure 4 .
Figure 4. Flowchart for calculation of effective impact velocity (V IS ).

FOR PEER REVIEW 10 of 18 Figure 5 .
Figure 5. Process parameter selection.

Figure 6 .
Figure 6.Comparison of interface morphology between experimental test and FEM model.

Figure 6 .
Figure 6.Comparison of interface morphology between experimental test and FEM model.

Figure 7 .
Figure 7. (a) Optical micrograph of a pulse-welded Al-SS304 member showing distinct zones at 10× magnification and (b) impact zone predicted by the FEM.

Figure 7 .
Figure 7. (a) Optical micrograph of a pulse-welded Al-SS304 member showing distinct zones at 10× magnification and (b) impact zone predicted by the FEM.

Figure 9 .
Figure 9. (a) Hardness and (b) plastic strain at the weld interface.

Figure 8 .
Figure 8. Lap-shear tests showing failure outside the welded zone.

Figure 9 .
Figure 9. (a) Hardness and (b) plastic strain at the weld interface.

Figure 9 .
Figure 9. (a) Hardness and (b) plastic strain at the weld interface.

Figure 10 .
Figure 10.Surface EDS map of the bi-metallic MPW joint.

Figure 10 .
Figure 10.Surface EDS map of the bi-metallic MPW joint.

Figure 11 .
Figure 11.EDS line scan at (a) crest of a wave, (b) trough of a wave, (c) flat interface.

Figure 11 .
Figure 11.EDS line scan at (a) crest of a wave, (b) trough of a wave, (c) flat interface.

Figure 12 .
Figure 12.Mechanical interlocking at the interface.

Figure A1 .
Figure A1.Comparison of experimental and simulated values of workpiece velocity for model validation: Case of compression joining of dissimilar materials (a) with variation of input energy

Figure A1 .
Figure A1.Comparison of experimental and simulated values of workpiece velocity for model validation: Case of compression joining of dissimilar materials (a) with variation of input energy and capacitance, data from [51], and (b) with variation of input voltage and air gap, data from [13,52].

Table 1 .
Summary of available literature on numerical and experimental studies in MPW of dissimilar metal joints.

Table 1 .
Summary of available literature on numerical and experimental studies in MPW of dissimilar metal joints.

Table 2 .
Process parameters employed in this work.

Table 3 .
Details of the FEM.

Table 2 .
Process parameters employed in this work.

Table 3 .
Details of the FEM.

Table 4 .
Chemical composition and material properties of flyer and target tubes and elec netic coil.

Table 4 .
Chemical composition and material properties of flyer and target tubes and electromagnetic coil.