Next Article in Journal
Numerical and Experimental Research on the Brushing Aluminium Alloy Mechanism Using an Abrasive Filament Brush
Next Article in Special Issue
FEM-Based Thermogram Correction for Inconel 625 Joint Hardness Clustering
Previous Article in Journal
Effect of Nanoparticles Surface Bonding and Aspect Ratio on Mechanical Properties of Highly Cross-Linked Epoxy Nanocomposites: Mesoscopic Simulations
Previous Article in Special Issue
Evolution of Microstructure in Friction Stir Processed Dissimilar CuZn37/AA5056 Stir Zone
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapid Calculation of Residual Stresses in Dissimilar S355–AA6082 Butt Welds

1
Department of Mechanical and Industrial Engineering, Norwegian University of Science and Technology, Richard Birkelands vei 2b, 7491 Trondheim, Norway
2
Department of Computational Materials Processing, Institute for Energy Technology, 2027 Kjeller, Norway
3
Department of Management and Engineering, University of Padova, Stradella S. Nicola, 2, I-36100 Vicenza, Italy
*
Author to whom correspondence should be addressed.
Materials 2021, 14(21), 6644; https://doi.org/10.3390/ma14216644
Submission received: 22 September 2021 / Revised: 26 October 2021 / Accepted: 28 October 2021 / Published: 4 November 2021

Abstract

:
An analytical model is proposed to rapidly capture the thermal and residual stresses values induced by the hybrid metal extrusion and bonding (HYB) process on dissimilar-metal butt-welded joints. The power input for two welding velocities is first assessed using a thermal–mechanical model solved by a heat generation routine written in MATLAB code. Subsequently, the obtained temperature history is used as input to solve the equilibrium and compatibility equations formulated to calculate the thermal and residual stresses. To verify the soundness of the analytical approach, a Finite Element numerical model of the entire process is carried out and results are compared with those coming from the proposed rapid method. It is found that the degree of accuracy reached by the analytical model is excellent, especially considering the tremendous time reduction when compared to that characterizing the standard numerical approach.

1. Introduction

Metal structures subjected to welding processes are characterized by residual stresses that are a consequence of heating and cooling cycles occurring in the joint and in the parent material. The localized heating induced by welding results in non-uniform temperature fields, which are associated with material thermal expansions and contractions. In the case of butt-welded plates, during the heating phase, compression stresses arise near the joint along the weld line. This happens because the thermal expansion of the hot area is constrained by the mass of the parent material that remains cold. Localized high plastic deformations in the metal, softened by the high temperature, are therefore promoted. If the metal deforms plastically in a non-uniform manner, the elastic recovery during the subsequent cooling phase does not fully occur, and the material remains in a state of residual tension. It is known that residual tensile stresses can contribute to brittle fractures or cause premature failures in components subjected to fatigue, while compression stresses are responsible for deformations and buckling in thin large panels, such as those used in the automotive industry. For these reasons, an assessment of the magnitude and distribution of residual stresses can be extremely helpful in the design phase of a welded structure [1,2].
In the past, welding residual stress investigation employed only experimental measurement methods. These methods can be divided in two categories: destructive and non-destructive. The first category comprises, among the others, hole drilling and strain gauge techniques [3,4,5,6], while non-destructive methods are X-ray diffraction (XRD), neutron diffraction, and ultrasonic waves [6,7,8]. However, all these experimental techniques are time-consuming and unable to capture the complete stress distributions over the entire component geometry. For instance, although the XRD method is highly accurate, it provides limited information about only a relatively thin surface layer [9,10].
In recent years, the development of numerical Finite Element codes made it possible to simulate the whole welding process, considering all the related effects such as the precipitates’ dissolution, solid-state phase transformation, and elastic-plastic and creep behavior of the material under investigation. For example, in his work, Deng [11] used a coupled thermal, metallurgical, and mechanical 3D finite element model to investigate the effects of solid-state phase transformation on welding residual stress and distortion in low carbon and medium carbon steels, revealing that the final residual stress and distortion in low carbon steel did not seem influenced by the solid-state phase transformation, while for the medium carbon steel, the final residual stresses and distortions appeared to be affected by the martensitic transformation. In addition, Chen et al. [12] developed a three-dimensional finite element model to investigate the thermal history and thermomechanical process of friction stir (FS) butt-welded AA6061-T6 plates. They investigated the evolution of longitudinal, transversal, and through-thickness stress and used X-ray diffraction to validate their results. D’Ostuni et al. [13] developed a FE model to simulate dissimilar aluminum–titanium fiber-lased welds employing two-dimensional and three-dimensional Gaussian heat sources. Finally, Obeid et al. [14] developed a 3D FE model to simulate a circumferential single-pass weld overlay and two-pass girth welding and validated the results coming from the model with experimental measurements. In general, numerical methods provide useful and accurate information, but they need a great computational effort and can be relatively time-consuming, depending on the dimension of the joint and the number of welding runs to be simulated.
In the case of thin plates butt-welded, some simplifications are possible, and this allows one to face the problem via analytical approaches. The system of equations governing the phenomenon can be solved by an iterative procedure easily implemented on current computers with programs, allowing drastic reductions in computational effort compared with more sophisticated non-linear FE analyses. The use of analytical models for the study of residual stresses has already been proposed by different authors in the past. Goff [15] simplified the problem by assuming constant thermal properties of the material and a linear temperature distribution in a direction perpendicular to the weld line. Tall [16] used the classical Rosenthal’s solution [17] for the temperature assessment and suggested an iterative method in which for each temperature increase, the thermal stresses needed to be in equilibrium; the results obtained for each temperature increment were then added to the stresses calculated in the previous time interval. Agapakis and Masubuchi [18] extended the Tall’s previous work by solving, with an iterative procedure, the constitutive and equilibrium equations that govern the problem. J. Cañas et al. [19,20] proposed a model in which the plates were divided into bars, and the governing equations were written in matrix form, allowing them to obtain a more efficient calculation algorithm. In all of these works, the welded plates did not present any constraint. Finally, Ferro et al. [21] proposed a more general model for the calculation of residual stresses induced by fusion welding, starting from the model proposed by Cañas et al. [19,20] and extending the model also to the case of fully constrained plates.
Recently, a new solid-state joining method known as the hybrid metal extrusion and bonding (HYB) process has appeared on the horizon [22]. So far, various experimental studies of the HYB process applied to similar and dissimilar metals welds were carried out in order to better understand and predict the effect of process parameters on welded materials; however, only recently was a great effort spent in the development of numerical models [23,24]. Recently, good experimental results were found for dissimilar aluminum–steel butt welds obtained by the HYB technique [25,26,27]. Numerical models were widely used in the literature to study the thermal field and residual stress arising during similar and dissimilar metal welds production; however, analytical investigations on dissimilar welded joints produced through this novel welding technique are still missing in the literature.
The present work is aimed at filling this gap by extending the analytical model proposed by Ferro [21] to dissimilar aluminum–steel butt-welded joints produced via HYB. The equations that solve the problem were written in matrix form and solved by an iterative procedure using MATLAB® code. A 3D finite element model was also developed that uses a double ellipsoid volume function with Gaussian power distribution to simulate the moving heat source. Results coming from the two different approaches, numerical and analytical, were finally compared with each other with the goal of verifying the soundness of the analytical method.

2. Materials and Methods

2.1. Geometry and Materials

A 4 mm thick aluminum alloy AA6082 plate butt-welded to a 4 mm thick structural steel S355 plate was considered. All dimensions are shown in Figure 1.
Materials properties, as a function of temperature, are shown in Figure 2 and Figure 3. More in detail, thermal properties were kept constant with temperature since previous studies have shown that this simplifying hypothesis gives acceptable results when considering the residual stress computation. The thermal conductivity, thermal diffusivity, and melting temperature of the aluminum alloy were set equal to 170 W/(mK), 62 mm2/s, and 638 °C, respectively. On the other side, with reference to same thermal properties, their values for the steel were 44 W/(mK), 11.5 mm2/s, and 1510 °C. The Poisson’s ratio was also kept constant with temperature and equal to 0.34 and 0.28 for aluminum and steel, respectively.
Material properties as a function of temperature can be also written in closed form as in Appendix A.
For the sake of simplicity, microstructural changes, with particular reference to secondary phases evolution in AA6082 aluminum alloy, were not considered herein. However, their effects on thermal and residual stresses were implicitly taken into account by temperature-dependent mechanical material properties.

2.2. The HYB Process

The HYB process is based on the principles of continuous extrusion and allows joining to be performed using aluminum filler metal (FM), similar to gas metal arc welding (GMAW) but without any melting involved. Figure 4 highlights the core parts forming the HYB PinPoint extruder. In the situation of butt-welding, the two plates to be joined are divided from each other by a fixed spacing so that an I-groove forms between them. During the welding operation, the PinPoint extruder shown in Figure 4 moves along the joint line at a constant travel speed, while the rotating pin with its moving dies is submerged between the plates. This enables the extrudate to flow downwards in the axial direction and mix with the base metal (BM) in the groove under conditions of high pressure and severe plastic deformation. Metallic bonding between the FM and the BM occurs as a result of these two contributions, along with oxide dispersion [22]. By taking into consideration the mass balance involved, the extrudate flow rate can be adjusted (by using the rotational speed of the drive spindle as the main process variable) in such a way that the groove between the plates is filled with solid aluminum in a continuous manner. More information about the HYB welding technique can be found elsewhere [25,29,30].

3. Modelling

Figure 5 shows a flowchart of the numerical/analytical work proposed. As a starting point, the welding conditions such as welding speed, rotational speed, and geometry of the tool were fixed. These input parameters were used to calculate the corresponding heat input by using a semi-analytical heat generation model already proposed by Leoni et al. [24]. The obtained heat input was then used both in FE and analytical model to calculate the thermal and residual stresses.
Table 1 shows the welding conditions investigated.

3.1. Thermal FE Model

The numerical model was developed using WELDSIM Finite Element code. The geometry was modeled with 47,032 and 17,000 hexahedral finite elements for the plates and the backing plate, respectively. Moreover, a finer mesh was created in the vicinity of the weld line in order to capture the high temperature gradients that characterize that region (Figure 6a) with the smallest element size of 0.75 mm. The effect of the filler material was also taken into account by element activation; finally, the moving heat source was modeled as a double ellipsoid power density distribution function, as proposed by Goldak [31]. The power density of the double-ellipsoid heat source describing the heat flux distribution in the front and rear of the heat source can be expressed as follows:
p f ( x , y , z , t ) = 6 3 p 0 k 1 k 2 k 3 , f π π e 3 x 2 / k 1 2 e 3 z 2 / k 2 2 e 3 ( y + v ( τ t ) ) 2 / k 3 , f 2
p r ( x , y , z , t ) = 6 3 p 0 k 1 k 2 k 3 , r π π e 3 x 2 / k 1 2 e 3 z 2 / k 2 2 e 3 ( y + v ( τ t ) ) 2 / k 3 , r 2
where x, y, and z (mm) are the coordinates as in Figure 6a, k 1 , k 2 , k 3 , f , k 3 , r (mm) are the half-width, the depth, and the front and rear length of the double ellipsoid heat source, respectively, τ (s) is the lag time, p 0 (W) is the net power input, and v (mm/s) the welding speed.
The heat transfer coefficient with the backing plate was considered to be 200 W/m2 °C, while convective heat loss was modeled using a convective heat transfer coefficient equal to 20 W/m2 °C, in agreement with a high air-flow assumption [32].

3.2. Thermal Analytical Model

The Rosenthal’s equation was used for the rapid thermal field calculation in the analytical model:
T = T 0 + p 0 / d 2 π λ exp [ v y 2 α ] K 0 ( v r 2 α )
In Equation (3), T 0 is the reference temperature set equal to 20 °C, y is the coordinate of a point of the plate in the same direction of the weld line, r is the distance from the heat source center, v is the welding speed, K 0 is the modified Bessel function of the second type and zero order, λ is the thermal conductivity, and α is the diffusivity. The Rosenthal’s Equation (3) is based on the following assumptions [33]:
  • the plates are considered semi-infinite and of thin thickness;
  • the equation describes the temperature range induced by a linear source: temperature gradient along the thickness of the plates is neglected;
  • the thermal field refers to a quasi-stationary condition of the welding process.
It is worth noting that p 0 , the actual thermal power of the pin, is calculated employing the model described by Leoni et al. [24] and summarized with the following equation:
p 0 = p 0 ( ω , τ m , T int , v , r t i p , r s h , d , h 1 , h 2 , m · , T F M )
where ω (rad/s) is the angular velocity of the tool, τ m (MPa) is the shear stress at the tool-matrix interface, T int (°C) is the temperature at the tool-matrix interface, v (mm/s) is the welding speed, r t i p (mm) is the radius of the tip of the tool, r s h (mm) is the radius of the shoulder of the tool, d (mm) is the thickness of the plates, h 1 (mm) is the difference between the radius of the tip and half width of the groove, h 2 (mm) is the difference between the radius of the shoulder and half width of the groove, m · (kg/s) is the flow of the hot filler material being extruded, and T F M (°C) is the temperature of the hot extrudate. Equation (4) describes the physical framework of the heat generation during the hybrid metal extrusion and bonding process (i.e., the extrusion of filler material and the friction between the tool and the base material). It is noted that the coupling between the mechanical (shear stress) and thermal field (T) (Equation (3)) is solved by an iterative routine written in MATLAB code and described elsewhere [24].

4. Residual Stress Assessment

4.1. Finite Element Analysis

Following the computational welding mechanical (CWM) approach, the temperature history at each node of the model was used as load for the mechanical computation. The mechanical boundary conditions are summarized in Figure 7 (Ux = 0 for x = 0 and for x = 2B, where Ux is the displacement in x direction, and B is the plate width). The plates were left free to deform in the y direction. In addition to the in-plane boundary constraints, in the FE model, the out of plane movement of the plates was constrained (Figure 7). The interface between the two plates was modelled as perfectly bonded as the weld proceeded. A sequential and uncoupled thermal and mechanical simulation was used since the thermo-mechanical coupling was already considered in the previous power input calculation (Equation (4)). It is worth noting that the materials constitutive laws were kept the same for both the numerical and analytical model. In particular, for the sake of simplicity, an elastic perfectly plastic law was used for the two materials to be welded.

4.2. Analytical Model

With reference to Figure 8, in order to write the equations that govern the problem in matrix form (constitutive, equilibrium, and compatibility equations), the welded plates were divided into a discrete number (n) of bars. For the sake of simplicity, all bars had the same width (b), while n was chosen sufficiently high to account for the temperature gradient near the weld line with a good approximation. It was assumed that at each time step the temperature of each bar was uniform, depending only on its x coordinate (i.e., T = T (x, t)). In other words, the temperature gradient in the y direction was neglected, a hypothesis that tends to be valid more and more as the welding speed increases. As for the FE model, boundary conditions were applied at the plate ends according to the HYB process (Ux = 0 for x = 0 and for x = 2B, where Ux is the displacement in x direction, and B is the plate width, while the boundary conditions in y direction are described by the compatibility Equation (11)). The Von Mises yield criterion was used, and materials viscous-plastic effects were neglected; the yield stress, the elastic modulus, and the thermal expansion coefficient of both materials were considered functions of temperature, as described by Equations (A1)–(A14).
If N y , i and N x , i are the forces on cross sections in bar i-th originating from the internal stresses and normal to the x and y axes, respectively, the equilibrium equations of forces and momentum can be written in matrix form (Figure 8) as follows [21]:
C ¯ T N ¯ y = 0 ¯
N x i = N x   f o r   i = 1 n
where
C ¯ = [ 1 d 1 1 d i 1 d n ]
N ¯ y = [ N y 1 N y i N y n ]
where n is the number of discrete bars, and d i is the x-distance of bar i-th from the origin (Figure 8).
Under the hypothesis of materials elastic perfectly plastic behavior, the constitutive equations are:
q ¯ y = q ¯ e y + q ¯ t y + q ¯ p y + Δ q ¯ p y
q i x = q e i x + q t i x + q p i x + Δ q p i x f o r i = 1 n
where q ¯ y is the vector containing the total elongations of each discrete bar for the longitudinal direction; q e i y is the elastic contribution to the elongation of the i-th bar governed by Hooke’s law; q t i y is the elongation of the i-th bar due to the thermal expansion; q p i y is the plastic deformation accumulated in the previous time steps for the i-th bar; Δ q p i y is the plastic deformation of the i-th bar at current time step; q i x is the total elongation of the i-th bar in the x direction, q e i x is the elastic contribution to the elongation in x direction of the i-th bar, and it is governed by Hooke’s law; q t i x is the elongation in x of the i-th bar due to the thermal expansion; q p i x is the plastic deformation in x direction for the i-th bar accumulated in the previous time steps; Δ q p i x is the plastic deformation of the i-th bar in x direction calculated at the current time step.
Finally, the compatibility equations are:
q ¯ y = C ¯ u ¯
i = 1 n q i x = 0
where u ¯ is the displacement vector associated with the degrees of freedom δ and θ , as shown in Figure 9.
Using the Von Mises criterion, the material yields when:
σ V M = σ y i e l d
where σ y i e l d is the yield stress of the alloy (i.e., S355 or AA6082) at the corresponding temperature.
Equation (9) can be written as follows:
q ¯ y = A ¯ N ¯ y A ¯ x + A ¯ t ( T ¯ T ¯ 0 ) + q ¯ p y + Δ q ¯ p y
where A ¯ is a diagonal matrix whose terms are ( L / ( E ( T i ) d b ) ) (where T i is the temperature of the i-th bar at the time step considered), A x ¯ is a vector in which each element is represented by the term ( ν N x / ( E ( T i ) d ) ), A t ¯ is another diagonal matrix whose terms are L α ( T i ) , and T ¯ is the vector containing the temperatures for each bar.
Using Equations (11) and (14), the vector N y ¯ can be expressed by the equation:
N ¯ y = A ¯ 1 ( C ¯ u ¯ + A ¯ x A ¯ t ( T ¯ T ¯ 0 ) q ¯ p y Δ q ¯ p y )
And from the combination of (15) and (5):
u ¯ = [ C ¯ T A ¯ 1 C ¯ ] 1 C ¯ T A ¯ 1 ( A ¯ x + A ¯ t ( T ¯ T ¯ 0 ) + q ¯ p y + Δ q ¯ p y )
Finally, from Equations (10) and (12), one can obtain the following equation:
N x = 1 b d L 1 E ( T i ) ( v b N i y E ( T i ) b α ( T i ) ( T i T 0 ) ( q p i x + Δ q p i x ) )
The unknown terms were Δ q p i y and Δ q p i x . The calculation routine started by assuming the increase in plastic deformation ( Δ q p i y ) as well as the force in the transverse (x) direction to be initially null. Using Equation (3), the temperature vector (temperature of each bar) was calculated, and therefore the material properties (which depend on temperature), and matrixes A ¯ , A x ¯ , and A t ¯ . Under such conditions, a first value of N y ¯ , using Equations (15) and (16), was determined. The following condition was then imposed:
| N i y | N p i f o r i = 1 n
where N p i indicates the yield force of the i-th bar (temperature dependent). It was then possible to calculate a first value of the plastic elongation ( Δ q p i y ) by using the previous value of the vector u ¯ and Equations (11) and (14).
This first approximation of plastic elongation was then used in (15), (16), and (18) for a second approximation of the vector N y ¯ , and such procedure was repeated until convergence was achieved for the time step considered.
Subsequently, a first approximation of Nx was calculated with Equation (17), using the values of N i y from the previous step and assuming to be null the increase of plastic deformation ( Δ q p i x ).
After having imposed the following condition (19) (derived from (13)),
[ ( N x d L ) 2 + ( N i y d b ) 2 N x N i y d 2 L b ] 0.5 σ y i e l d , i
and using Equation (17), a first approximation of the term Δ q p i x was obtained, and this was used for a second approximation of Nx; this was repeated until convergence was achieved. The value of Nx was then used to determine the new value of N y ¯ , and if the last value was not equal to the previous one (within the tolerance imposed), Nx and Δ q p i x were recalculated using the last obtained value of N y ¯ ; this iteration was repeated until convergence was achieved.
Ferro et al. [21] showed that the term coupling the longitudinal with the transversal stress (i.e., ν Nx/L) is much smaller than the main term, N i y / b , and that the transversal stress is no longer negligible only below 300 °C, where the yield strength is not significantly affected by the temperature (Figure 3). This means that, at each time step, the first calculated value of Nx does not have a significant influence on the vector N y ¯ obtained by assuming initially Nx = 0, and therefore N y ¯ is limited only by the yield stress of the material.
The algorithm described (see also Figure A1 in Appendix A) was written in MATLAB® code. The input data required by the program were the geometry of the tool (radius of the pin and shoulder), the rotational speed of the tool ( ω ), the welding speed ( v ), the plate size ( L , 2 B , d ), the groove width, and the total number of bars ( n ) (in this work, b = 1 mm, while the time step ( Δ t ) was set equal to 1 s).

5. Results

5.1. Temperature Results

In Figure 10, the analytical solution of the temperature evolution along the direction perpendicular to the weld line and at different distances from the heat source is compared with that coming from FE model. The results of the two models showed good agreement despite some small differences in the weld line where the Rosenthal’s solution (asymptotic) overestimates the temperature field. However, it is known that no melting occurs during the process [26]; therefore, the latent heat was not considered in the simulations. Moreover, since the analytical solution was considered uncoupled for aluminum and steel, a discontinuity at the interface is observed. In Figure 11 the temperature evolutions at points taken at the same distance from the weld line for both sides are shown. As expected, the analytical solution approximates better the temperature at the aluminum alloy side since it is already known that the Rosenthal’s equation is more suitable for materials that have high thermal conductivity and thermal diffusion coefficients [33]. The differences between analytical and FE results may be due to the natural simplifications required by the Rosenthal’s equation solution, such as the semi-infinite plate geometry deriving from the quasi-stationary thermal field distribution hypothesis. However, despite these apparent discrepancies between the temperatures of the two approaches, the results are in overall agreement, especially if one considers that the most important effect comes from the distribution of the temperature, which finds good agreement between the two methods.

5.2. Thermal and Residual Stress Results

Figure 12 compares the evolution of the longitudinal stress obtained with the numerical and the analytical method proposed. For the comparison with the analytical model, FEM results were taken from the mid-thickness of a cross section at the mid-length of the plates. Because the plates were thin, the residual stress was assumed constant throughout the thickness. A good match was found, and the validity of the procedures adopted, as well as the hypotheses formulated, were confirmed.
In more detail, the curve showing the longitudinal residual stress value as a function of the transversal direction (x) coming from the analytical model overlaps with that obtained using the FE model (Figure 13). At the interface between the two materials, a discontinuity is observed, which is due to the different yield stresses characterizing the two materials. Again, good agreement was observed for both the chosen welding speeds. The slightly wider plasticized zone found with the analytical solution, steel side, can be attributed to the temperature value in that zone, which tends to be overestimated by the Rosenthal’s equation with respect to the more accurate FE calculation (Figure 13). Except for this minor difference, the solutions were in good agreement.
Finally, a comparison between the computational time required by the two methods is presented in Table 2 (RAM: 32 GB, Processor: Intel® Xeon® E-2174G CPU 3.8 GHz)

6. Conclusions

A rapid analytical model for the prediction of residual stresses was developed. The hypotheses on which the analytical model formulation is based were validated by comparing the results obtained with a more accurate FE analysis. The numerical analysis was carried out using a 3D model and a double ellipsoid heat source with Gaussian power distribution, while the Rosenthal’s equation was used to calculate the thermal field in the analytical model. In the numerical analysis a thermal source was applied along the welding line, having a constant travel speed, while in the analytical model it was assumed that at each time step each section parallel to the weld line was subjected to the same temperature across its entire length (T = T (x, t)). Despite this substantial difference, a good agreement was found between the thermal and mechanical results coming from the two approaches. The simplified model was found to be more advantageous in terms of computational time; it allows one to carry out parametric studies, and it can conceivably be further developed to analyze the effects of different welding setup or stress relief heat treatments on the residual stresses induced by the welding process. In conclusion, although the FE analysis allows one to avoid drastic simplifying hypotheses and provides a series of more detailed information, the proposed simplified method becomes more attractive when a rapid estimation of the thermal and residual stress arising in dissimilar metal welds is required.

Author Contributions

Conceptualization, F.L. and P.F.; methodology, F.L.; software, F.L. and H.G.F.; validation, F.L.; formal analysis, F.L.; investigation, F.L.; data curation, F.L.; writing—original draft preparation, F.L.; writing—review and editing, F.L., P.F., H.G.F., and F.B.; visualization, F.L.; supervision, P.F. and F.B.; project administration, F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

SymbolValueDescription
A ¯ E from Equations (A5)–(A8), Equation (A13),
L = 500 mm
d = 4 mm
b = 1 mm
Diagonal matrix containing ( L / ( E ( T i ) d b ) ) terms
A t ¯ L = 500 mm
α from Equations (A9) and (A14)
T from Equation (3)
Diagonal matrix containing ( L α ( T i ) ) terms
A x ¯ υsteel = 0.28
υalu = 0.34
E from Equations (A5)–(A8) and (A13)
d = 4 mm
Diagonal matrix containing ( ν N x / ( E ( T i ) d ) ) terms
B 100 mmWidth of the plate, mm
b 1 mmWidth of the discrete bars used in the analytical model, mm
C ¯ Equation (7) Vector containing the values of the distance of each bar from the origin
d 4 mmPlates thickness, mm
E A l u Equation (A13)Elastic modulus of aluminum, MPa
E S t e e l Equations (A5)–(A8)Elastic modulus of steel, MPa
h 1 0.75 mmDifference between the radius of the tip of the tool and half width of the groove, mm
h 2 3.5 mmDifference between the radius of the shoulder of the tool and half width of the groove, mm
K 0 FunctionBessel function of the second kind and zero order
k 1 4 mmHalf-width of welding pool for the Goldak’s equation, mm
k 2 4 mmDepth of welding pool for the Goldak’s equation, mm
k 3 , f 5 mmHalf-length of welding pool for the Goldak’s equation, mm
k 3 , r 5 mmHalf-length of welding pool for the Goldak’s equation, mm
L 500 mmLength of the plates, mm
m · 5.2 × 10−4 kg/s
6.9 × 10−4 kg/s
Mass flow of the extrudate, kg/s
n 200Number of discrete bars
N x Equation (17)Force in x direction, N
N y ¯ Equation (15)Vector containing the forces in y direction for each bar, N
p 0 1000 and 1100 WPower input, W
p f Equation (1) Power density of the double-ellipsoid heat source, front, W/mm3
p r Equation (2) Power density of the double-ellipsoid heat source, rear, W/mm3
q e , i x CalculatedElastic deformation of i-th bar in the x direction, mm
q e y ¯ CalculatedVector containing the elastic deformations in the y direction, mm
q i x CalculatedDeformation of i-th bar in the x direction, mm
q p , i x CalculatedPlastic deformation of the i-th bar in the x direction, mm
q p y ¯ CalculatedVector containing plastic deformations in the y direction, mm
q t , i x CalculatedThermal deformation of the i-th bar in the x direction, mm
q t y ¯ CalculatedVector containing the thermal deformations in the y direction, mm
q y ¯ CalculatedVector containing the deformations in the y direction, mm
r Generic pointDistance of a generic point from the heat source center, mm
r s h 5.5 mmPin shoulder radius, m
r t i p 2.75 mmPin tip radius, m
T Equation (3)Temperature, °C
T ¯ From Equation (3)Vector containing the temperature of each bar, °C
T 0 20 °CReference temperature, °C
T ¯ 0 Vector of T0Vector containing the reference temperature of each bar, °C
T F M 450 °CTemperature of the hot extrudate, °C
t VariableTime, s
u ¯ Vector with δ and θ Vector of the degrees of freedom in the compatibility equation
U x 0 mmDisplacement in x direction
v 12 and 16 mm/sWelding speed, mm/s
x -x-coordinate, transverse direction, mm
X -coordinate in the local heat source reference system, mm
y -y-coordinate, welding direction, mm
Y -coordinate in the local heat source reference system, mm
z -z-coordinate, mm
α Equations (A9) and (A14)Thermal diffusivity, mm2/s
α s t e e l Equation (A9)Thermal diffusivity of steel, mm2/s
α A l u Equation (A14)Thermal diffusivity of aluminum, mm2/s
δ CalculatedFirst degree of freedom of the compatibility equation, mm
Δ t 1 sTime step, s
Δ q p , i x CalculatedPlastic deformation in x direction at the current time step, mm
Δ q p y ¯ CalculatedVector of plastic deformation in y direction at the current time step, mm
λ See belowThermal conductivity, W/mm°C
λ A l u 170 W/mm°CThermal conductivity of aluminum, W/mm°C
λ S t e e l 44 W/mm°CThermal conductivity of steel, W/mm°C
ω 37 rad/sAngular rotation speed, rad/s
σ y i e l d , A l u Equations (A10)–(A12)Yield stress of aluminum, MPa
σ y i e l d , S t e e l Equations (A1)–(A4)Yield stress of steel, MPa
σ V M CalculatedVon Mises stress, MPa
τ 20.8 and 15.6 sLag factor in Goldak’s equation, s
τ m CalculatedShear stress at tool-matrix interface, MPa
θ CalculatedSecond degree of freedom of the compatibility equation, rad
υ υ s t e e l = 0.28
υ A l u = 0.34
Poisson ratio

Appendix A

σ y i e l d , S t e e l = 350 f o r T 420
σ y i e l d , S t e e l = 1.127 E 5 T 3 0.01942 T 2 + 10 T 1254.3 f o r 420 < T 700
σ y i e l d , S t e e l = 165.6 0.1037 T f o r 700 < T 1500
σ y i e l d , S t e e l = 0.01 f o r 1500 < T
E S t e e l = 206000 44.44 T f o r T 200
E S t e e l = 26300 325 T f o r 200 < T 700
E S t e e l = 1.229 E 5 159.9 T + 0.05094 T 2 f o r 700 < T 1200
E S t e e l = 5000 f o r 1200 < T
α S t e e l = 1.29 E 5 + 2.218 E 9 T
σ y i e l d , A l u = 120 0.03392 T f o r T 300
σ y i e l d , A l u = 283 0.58 T f o r 300 < T 460
σ y i e l d , A l u = 69 0.107 T f o r 460 < T 638
E A l u = 0.1155 T 2 41.06 T + 71000
α A l u = 2.188 e 8 T + 2.25 e 5
Figure A1. Flow chart of the analytical model.
Figure A1. Flow chart of the analytical model.
Materials 14 06644 g0a1

References

  1. Leggatt, R.H. Residual stresses in welded structures. Int. J. Press. Vessel. Pip. 2008, 85, 144–151. [Google Scholar] [CrossRef]
  2. Masubuchi, K. Analysis of Welded Structures: Residual Stresses, Distortion, and Their Consequences; Elsevier: Amsterdam, The Netherlands, 2013; Volume 33, ISBN 1483188434. [Google Scholar]
  3. Pang, H.L.; Pukas, S.R. Residual stress measurements in a cruciform welded joint using hole drilling and strain gauges. Strain 1989, 25, 7–14. [Google Scholar] [CrossRef]
  4. Niku-Lari, A.; Lu, J.; Flavenot, J.F. Measurement of residual-stress distribution by the incremental hole-drilling method. J. Mech. Work. Technol. 1985, 11, 167–188. [Google Scholar] [CrossRef]
  5. Schajer, G.S. Advances in hole-drilling residual stress measurements. Exp. Mech. 2010, 50, 159–168. [Google Scholar] [CrossRef]
  6. Withers, P.J.; Bhadeshia, H. Residual stress. Part 1–measurement techniques. Mater. Sci. Technol. 2001, 17, 355–365. [Google Scholar] [CrossRef]
  7. Chu, S.L.; Peukrt, H.; Schnider, E. Residual stress in a welded steel plate and their measurement using ultrasonic techniques. MRL Bull. Res. Dev. 1987, 1, 45–50. [Google Scholar]
  8. Turan, M.E.; Aydin, F.; Sun, Y.; Cetin, M. Residual stress measurement by strain gauge and X-ray diffraction method in different shaped rails. Eng. Fail. Anal. 2019, 96, 525–529. [Google Scholar] [CrossRef]
  9. Leoni, M.; Scardi, P.; Rossi, S.; Fedrizzi, L.; Massiani, Y. (Ti, Cr) N and Ti/TiN PVD coatings on 304 stainless steel substrates: Texture and residual stress. Thin Solid Films 1999, 345, 263–269. [Google Scholar] [CrossRef]
  10. Suh, C.-M.; Hwang, B.-W.; Murakami, R.-I. Behaviors of residual stress and high-temperature fatigue life in ceramic coatings produced by PVD. Mater. Sci. Eng. A 2003, 343, 1–7. [Google Scholar] [CrossRef]
  11. Deng, D. FEM prediction of welding residual stress and distortion in carbon steel considering phase transformation effects. Mater. Des. 2009, 30, 359–366. [Google Scholar] [CrossRef]
  12. Chen, C.M.; Kovacevic, R. Finite element modeling of friction stir welding—thermal and thermomechanical analysis. Int. J. Mach. Tools Manuf. 2003, 43, 1319–1326. [Google Scholar] [CrossRef]
  13. 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. [Google Scholar] [CrossRef] [Green Version]
  14. Obeid, O.; Alfano, G.; Bahai, H.; Jouhara, H. Experimental and numerical thermo-mechanical analysis of welding in a lined pipe. J. Manuf. Process. 2018, 32, 857–872. [Google Scholar] [CrossRef]
  15. Goff, R.F.D.P. A simplified analysis of the residual longitudinal stresses and strains due to the gas-cutting and welding of thin steel plate. Int. J. Mech. Sci. 1979, 21, 287–300. [Google Scholar] [CrossRef]
  16. Tall, L. Residual Stresses in Welded Plates: A Theoretical Study; Fritz Engineering Laboratory, Department of Civil Engineering, Lehigh University: Bethlehem, PA, USA, 1961. [Google Scholar]
  17. Rosenthal, D. Mathematical theory of heat distribution during welding and cutting. Weld. J. 1941, 20, 220–234. [Google Scholar]
  18. Agapakis, J.E.; Masubuchi, K. Analytical Modeling of Thermal Stress Relieving in Stainless and High-Strength Steel Weldments. Weld. J. 1984, 63, 187. [Google Scholar]
  19. Cañas, J.; Picòn, R.; París, F.; Del Río, J.I. A one-dimensional model for the prediction of residual stress and its relief in welded plates. Int. J. Mech. Sci. 1996, 38, 735–751. [Google Scholar] [CrossRef]
  20. Canas, J.; Picon, R.; Pariis, F.; Blazquez, A.; Marin, J.C. A simplified numerical analysis of residual stresses in aluminum welded plates. Comput. Struct. 1996, 58, 59–69. [Google Scholar] [CrossRef]
  21. Ferro, P.; Tiziani, A.; Bonollo, F. Analisi numerica e teorica del campo di tensioni residue indotte dal processo di saldatura in piastre in lega leggera. In Proceedings of the AIAS, Parma, Italy, 18–21 September 2002. [Google Scholar]
  22. Sandnes, L.; Grong, Ø.; Torgersen, J.; Welo, T.; Berto, F. Exploring the hybrid metal extrusion and bonding process for butt welding of Al–Mg–Si alloys. Int. J. Adv. Manuf. Technol. 2018, 98, 1059–1065. [Google Scholar] [CrossRef] [Green Version]
  23. Leoni, F.; Grong, Ø.; Fjær, H.G.; Ferro, P.; Berto, F. A First Approach on Modelling the Thermal and Microstructure Fields During Aluminium Butt Welding Using the HYB PinPoint Extruder. Procedia Struct. Integr. 2020, 28, 2253–2260. [Google Scholar] [CrossRef]
  24. Leoni, F.; Grong, Ø.; Ferro, P.; Berto, F. A Semi-Analytical Model for the Heat Generation during Hybrid Metal Extrusion and Bonding (HYB). Materials 2021, 14, 170. [Google Scholar] [CrossRef]
  25. Berto, F.; Sandnes, L.; Abbatinali, F.; Grong, Ø.; Ferro, P. Using the Hybrid Metal Extrusion & Bonding (HYB) Process for Dissimilar Joining of AA6082-T6 and S355. Procedia Struct. Integr. 2018, 13, 249–254. [Google Scholar]
  26. Bergh, T.; Sandnes, L.; Johnstone, D.N.; Grong, Ø.; Berto, F.; Holmestad, R.; Midgley, P.A.; Vullum, P.E. Microstructural and mechanical characterisation of a second generation hybrid metal extrusion & bonding aluminium-steel butt joint. Mater. Charact. 2021, 173, 110761. [Google Scholar]
  27. Sandnes, L.; Bergh, T.; Grong, Ø.; Holmestad, R.; Vullum, P.E.; Berto, F. Interface microstructure and tensile properties of a third generation aluminium-steel butt weld produced using the Hybrid Metal Extrusion & Bonding (HYB) process. Mater. Sci. Eng. A 2021, 809, 140975. [Google Scholar]
  28. Outinen, J.; Kesti, J.; Mäkeläinen, P. Fire design model for structural steel S355 based upon transient state tensile test results. J. Constr. Steel Res. 1997, 42, 161–169. [Google Scholar] [CrossRef]
  29. Grong, Ø. Recent Advances in Solid-State Joining of Aluminium. Weld. J. 2012, 91, 26–33. [Google Scholar]
  30. Grong, Ø.; Sandnes, L.; Berto, F. A status report on the hybrid metal extrusion & bonding (HYB) process and its applications. Mat. Des. Process 2019, 1, e41. [Google Scholar]
  31. Goldak, J.; Chakravarti, A.; Bibby, M. A new finite element model for welding heat sources. Metall. Trans. B 1984, 15, 299–305. [Google Scholar] [CrossRef]
  32. Goyal, R.; Johnson, E.; El-Zein, M.; Goldak, J.; Coulombe, M.; Tchernov, S. A model equation for the convection coefficient for thermal analysis of welded structures. In Trends in Welding Research, Proceedings of the 8th International Conference on Trends in Welding Research, Pine Mountain, GA, USA, 1–6 June 2008; ASM International: Materials Park, OH, USA, 2009; pp. 321–327. [Google Scholar]
  33. Grong, Ø. Metallurgical Modelling of Welding; Institute of Materials: London, UK, 1997. [Google Scholar]
Figure 1. Geometry of the welded plates.
Figure 1. Geometry of the welded plates.
Materials 14 06644 g001
Figure 2. Material properties of S355 steel [28].
Figure 2. Material properties of S355 steel [28].
Materials 14 06644 g002
Figure 3. Material properties of AA6082 aluminum alloy.
Figure 3. Material properties of AA6082 aluminum alloy.
Materials 14 06644 g003
Figure 4. Schematic illustration of the HYB process applied to dissimilar metals with its main parts being involved in the joining process.
Figure 4. Schematic illustration of the HYB process applied to dissimilar metals with its main parts being involved in the joining process.
Materials 14 06644 g004
Figure 5. Flowchart showing the steps involved in calculating the net power input using the heat generation model and the subsequent Finite Element and analytical/numerical analysis for each corresponding welding condition.
Figure 5. Flowchart showing the steps involved in calculating the net power input using the heat generation model and the subsequent Finite Element and analytical/numerical analysis for each corresponding welding condition.
Materials 14 06644 g005
Figure 6. (a) Mesh used for the FE simulations. (b) Heat flux at top layer (z = 0) and in the local reference system using the double ellipsoid power density distribution function, as proposed by Goldak and using the following parameters: p 0 = 1000 W, k 1 = 4 mm, k 2 = 4 mm, k 3 , f = k 3 , r = 5 mm.
Figure 6. (a) Mesh used for the FE simulations. (b) Heat flux at top layer (z = 0) and in the local reference system using the double ellipsoid power density distribution function, as proposed by Goldak and using the following parameters: p 0 = 1000 W, k 1 = 4 mm, k 2 = 4 mm, k 3 , f = k 3 , r = 5 mm.
Materials 14 06644 g006
Figure 7. Schematic representation of the welded plates with lateral and out-of-plane constraints (uniformly applied to the top and bottom surfaces of the plates) assumed in the FE model.
Figure 7. Schematic representation of the welded plates with lateral and out-of-plane constraints (uniformly applied to the top and bottom surfaces of the plates) assumed in the FE model.
Materials 14 06644 g007
Figure 8. Schematic representation of the welded plates with lateral constraints and weld line highlighted in red.
Figure 8. Schematic representation of the welded plates with lateral constraints and weld line highlighted in red.
Materials 14 06644 g008
Figure 9. Schematic representation of the deformed and undeformed configuration for the welded plates.
Figure 9. Schematic representation of the deformed and undeformed configuration for the welded plates.
Materials 14 06644 g009
Figure 10. Temperature at different time steps for the welding speed of (a,c) 12 mm/s and (b,d) 16 mm/s and power input of (a,c) 1000 W and (b,d) 1100 W, respectively, calculated via (a,b) FE model and (c,d) via rapid analytical method. Temperatures values were taken considering the mid-thickness of a cross section taken orthogonal to the welding direction in the middle length of the plates.
Figure 10. Temperature at different time steps for the welding speed of (a,c) 12 mm/s and (b,d) 16 mm/s and power input of (a,c) 1000 W and (b,d) 1100 W, respectively, calculated via (a,b) FE model and (c,d) via rapid analytical method. Temperatures values were taken considering the mid-thickness of a cross section taken orthogonal to the welding direction in the middle length of the plates.
Materials 14 06644 g010
Figure 11. Temperature–time plots exported from two points at the same distance from the weld line (i.e., 6 mm, see schematic representation at the bottom graphs) considering the two welding conditions (i.e., welding speeds of 12 and 16 mm/s and power inputs of 1000 and 1100 W, respectively) for both (a) aluminum and (b) steel sides.
Figure 11. Temperature–time plots exported from two points at the same distance from the weld line (i.e., 6 mm, see schematic representation at the bottom graphs) considering the two welding conditions (i.e., welding speeds of 12 and 16 mm/s and power inputs of 1000 and 1100 W, respectively) for both (a) aluminum and (b) steel sides.
Materials 14 06644 g011
Figure 12. Longitudinal (y direction) stress at different time steps for the welding speed of (a,c) 12 mm/s and (b,d) 16 mm/s and power input of (a,c) 1000 W and (b,d) 1100 W, respectively, calculated via (a,b) FE model and (c,d) via rapid analytical method. Values of stress refer to the mid-thickness of a cross section taken orthogonal to the welding direction in the middle length of the plates.
Figure 12. Longitudinal (y direction) stress at different time steps for the welding speed of (a,c) 12 mm/s and (b,d) 16 mm/s and power input of (a,c) 1000 W and (b,d) 1100 W, respectively, calculated via (a,b) FE model and (c,d) via rapid analytical method. Values of stress refer to the mid-thickness of a cross section taken orthogonal to the welding direction in the middle length of the plates.
Materials 14 06644 g012
Figure 13. (a) Comparison between the residual longitudinal (y direction) stress calculated with the two methods (i.e., FEM and analytical model), considering a cross section taken at the mid-length of the plates and for a welding speed of 12 mm/s and a power input of 1000 W. (b) Comparison between the residual longitudinal stress calculated with the two methods (i.e., FEM and analytical model), considering a cross section taken at the mid-length of the plates and for a welding speed of 16 mm/s and a power input of 1100 W.
Figure 13. (a) Comparison between the residual longitudinal (y direction) stress calculated with the two methods (i.e., FEM and analytical model), considering a cross section taken at the mid-length of the plates and for a welding speed of 12 mm/s and a power input of 1000 W. (b) Comparison between the residual longitudinal stress calculated with the two methods (i.e., FEM and analytical model), considering a cross section taken at the mid-length of the plates and for a welding speed of 16 mm/s and a power input of 1100 W.
Materials 14 06644 g013
Table 1. Welding conditions modelled.
Table 1. Welding conditions modelled.
MaterialsWelding TechniqueWelding SpeedPin Rotational Speed
AA6082–SS355HYB12 and 16 mm/s350 RPM
Table 2. Computational time required by the methods used.
Table 2. Computational time required by the methods used.
FEMAnalytical
64,111 s (~18 h)~10 s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Leoni, F.; Fjær, H.G.; Ferro, P.; Berto, F. Rapid Calculation of Residual Stresses in Dissimilar S355–AA6082 Butt Welds. Materials 2021, 14, 6644. https://doi.org/10.3390/ma14216644

AMA Style

Leoni F, Fjær HG, Ferro P, Berto F. Rapid Calculation of Residual Stresses in Dissimilar S355–AA6082 Butt Welds. Materials. 2021; 14(21):6644. https://doi.org/10.3390/ma14216644

Chicago/Turabian Style

Leoni, Francesco, Hallvard Gustav Fjær, Paolo Ferro, and Filippo Berto. 2021. "Rapid Calculation of Residual Stresses in Dissimilar S355–AA6082 Butt Welds" Materials 14, no. 21: 6644. https://doi.org/10.3390/ma14216644

APA Style

Leoni, F., Fjær, H. G., Ferro, P., & Berto, F. (2021). Rapid Calculation of Residual Stresses in Dissimilar S355–AA6082 Butt Welds. Materials, 14(21), 6644. https://doi.org/10.3390/ma14216644

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