Next Article in Journal
Nitration of 2,6,8,12-Tetraacetyl-2,4,6,8,10,12-Hexaazaisowurtzitane Derivatives
Previous Article in Journal
Effect of High-Current Pulsed Electron Beam on Properties of Graphene-Modified Aluminum Titanium Carbide Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Microstructural Approach to the Strength of an Explosion Weld

by
Alexander G. Kolpakov
1,* and
Sergei I. Rakin
2
1
SysAn (System Analysis in Engineering), 630075 Novosibirsk, Russia
2
Mathematics Department, Siberian Transport University, 630049 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Materials 2022, 15(22), 7878; https://doi.org/10.3390/ma15227878
Submission received: 6 October 2022 / Revised: 30 October 2022 / Accepted: 4 November 2022 / Published: 8 November 2022
(This article belongs to the Section Mechanics of Materials)

Abstract

:
In this paper, the local stress–strain state in an explosion weld was investigated and the local strength of the welded materials near the weld analyzed. It follows from the experimental data that the explosion weld at the microlevel looks like a wavy line. In the first approximation, this wavy line may be assumed to be periodic. We used the two-scale method to analyze the corresponding interface elasticity problem. We carried out numerical computations for three of the most referenced types of weld geometry: the symmetric wave, the asymmetric wave, and the wave with crest. We found that the wave geometry of the weld leads to increase in local stress in the weld zone. The stress concentration varied from 20% to 200% in dependence on the weld geometry and the macroscopic loading. Explosion welding is accompanied by strain hardening of the materials in the welding zone. In some cases, the strain hardening may compensate for the increasing local stress. As a result, the weld may be both stronger and weaker than the welded materials.

1. Introduction

An explosion weld occurs when metal plates collide due to plastic deformation in the impact zone. The first works on explosion welding were carried out in the USA and the USSR and date back to the period 1958–1963, see [1,2,3,4,5,6,7,8,9,10,11]. Interest in explosion welding continues unabated to this day.
The strength of a structure/material produced by the welding is determined by the strength of both the welded materials and the weld. Many publications indicate that the bonding strength of explosion welded metals exceeds the strength of the weakest material of the pair (see, e.g., [12]) and destruction occurs at a certain distance from the plane of the connection along the least durable metal (see [11,12,13,14]). In [13,14], it is explained by the strengthening of the material in/near the weld as a result of its plastic deformation. At the same time, some publications report destruction along the weld, see, e.g., [15]. In this paper, we investigated the local stress–strain state in an explosion weld and analyzed the local strength of the welded materials in the vicinity of the weld.
When welding with an explosion of bodies, at the macro level, the weld looks like a straight line. However, at the micro level, the weld is wavy, see Figure 1. We call the wave-like weld and its surroundings the welding zone P 0 , see Figure 1. The parts P and P + in Figure 1d of the bodies outside the welding zone are called the main parts of the materials.
At a microlevel, the weld D ε has wavy geometry, which can be taken as periodic. The length of the period and the wave amplitude are of the order of several hundred micrometers, see Figure 2. Wave shapes vary from sinusoidal to asymmetric and to a wave with a crest (see Figure 1a–c). The wavelength and amplitude of the weld wave are of the order of hundred micrometers. These dimensions are compared with the dimension of the microstructure in metals [13]. The grain size of metals varies from 2.7 μm (0.0027mm.) to 1 mm. It is considered that large grains have a diameter from 62 μm (0.062 mm) to 1 mm, while small grains have a diameter from 2.7 μm to 44 μm, see, e.g., [14]. For example, aluminum grain diameters range from 10 μm to 100 μm. In the welds obtained by thermal welding of metals, a granular structure is different from the structure of welded metals, see [15] for details.
According to Figure 2a–c, the microstructure of the materials in an explosion weld is similar to the microstructure of the materials in the main parts of the welded materials, see Figure 2. The explosion weld differs from a fusion weld (for the fusion weld see, e.g., [13,15]).
Using the same scale in all figures in Figure 2, in Figure 2a–c (figures are taken from [16,17]), the microstructure size of the materials is compared to the wavelength of the weld. In Figure 2d–f, the microstructure size is small compared to the amplitude of the weld. In Figure 2a–c, the non-uniform microstructure is clearly seen, while in Figure 2d–f, the microstructure is practically invisible. Figure 2a displays the microstructure of the pair steel–steel, Figure 2b displays the microstructure of the pair copper–aluminum. For the materials displayed in Figure 2a–c, the first-order approximation is two homogeneous bodies connected along a wavy line (the model displayed in Figure 1c). The second-order approximation accounts for the inhomogeneous structure of the connected bodies. For the materials displayed in Figure 2d,e, the model displayed in Figure 1c may be accepted as an accurate model.
The strength of the weld is determined by the following factors:
  • stress concentration caused by the wavy geometry of the joint
  • change of the strength of the welded materials as a result of the intensive plastic deformation in the process of the explosion welding.

2. Formulation and Two-Scale Analysis of the Problem

In this paper, we analyzed in detail the local stress–strain state caused by the wavy geometry of the explosion weld. We distinguished the weld region and the main parts of the materials, which are placed far from the weld.
Although explosion welding shows a clear boundary between the welded bodies, it is impossible to separate the local stress–strain state in the weld zone D ε and the adjacent “left” P and “right” P + main parts of the bodies. To describe the stress–strain state in the entire region P D ε P + , we used the “local perturbations” method, developed in [18,19]. The “local disturbances” are a set of functions defined in the domain P D ε P + that are fast decaying with distance from the weld zone D ε . The “local disturbances” describe both the classical boundary layers near the weld and the deformations of general form in the weld zone D ε .
The real weld D ε is a wavy surface (or line, if a planar elasticity problem is considered). The characteristic dimension of the wave period and amplitude is denoted by ε . The parameter ε is small: ε 1 , and is formalized as ε 0 . We assume that the weld is periodic with period ε T with respect to variable x 2 . Due to periodicity of the weld, the entire problem is also periodic with respect to the variable x 2 with period ε T . In the “fast” variables, the problem is periodic with respect to the variable y 2 with period T (Figure 1c). We assume, in the first approximation, that the elastic constants of the materials in the welded zone are the same as in the main parts of the materials.
At the macrolevel, the weld looks like a surface (or line). We assume that this surface (line) is flat, see Figure 1a,b.
Denote x = (s, t, n) = (x1, x2, x3). the “slow” coordinates in the orthonormal coordinate system. The notation ( s , t , n ) is used to separate the coordinates into the “tangential” coordinates ( s , t ) in the plane of the weld D and is the “normal” coordinate n in the direction perpendicular to the weld, see Figure 1a. Following [20], we introduce “fast” variables y = ( y 1 , y 2 , y 3 ) = x / ε . We assume that the connecting area has a periodic structure in the variables ( y 1 , y 2 ) with the periodicity cell T , see Figure 2b. In the macroscopic coordinate system ( x 1 , x 2 , x 3 ) , the periodicity cell T turns into ε T .
The explosion weld is a microscopic object [21]. This means that, although the weld is small in size, it may be investigated by the methods of classical continuum mechanics [22]. The approaches based on the molecular and atomic models (see such models in [23,24,25,26]) are not intended for this paper. In this paper, we take the elasticity theory problem as the starting point of our investigation. The corresponding problem has the form:
{ ( a i j k l ( x ) u k , l y ) , i y = q i ( x )   in   Q , u ( x ) = 0   on   Q .  
The wavy structure of the weld is taken into account by setting the local elastic constants a i j k l ( x ) in (1) in the form:
a i j k l ( x ) = { a i j k l   on   the   left   from   D ε , a i j k l +   on   the   right   from   D ε .
The local elastic constants a i j k l ( x ) are periodic with respect to the variable x 2 with the period ε T . It means that a i j k l ( x ) are the high-oscillating functions.
If welded bodies are subjected to a macroscopic force, at the macro-level, boundary layers [27,28,29] arise on the left and on the right of the weld D . At the micro-level, the weld D turns into a weld zone D ε and, in addition to the “left” and the “right” boundary layers, the local deformation occurs in the weld zone D ε placed between the “left” and the “right” boundary layers. The local stress–strain state inside the weld zone D ε is of general form (see Figure 3, Figure 4 and Figure 5 below) and demonstrates no similarity with the boundary layers. The boundary layers have classical form at a distance from the weld zone D ε , in the domains between the weld zone D ε and the main parts P and P + of the connected bodies. The local stress–strain state in the weld zone D ε is a micro-structural solution, while the boundary layers are macro-structural solutions.
Following [18,19], we seek the solution to the problem in the form:
u ( x ) = u 0 ( x ) + ε u 1 ( x / ε ) .
In (3), u 0 ( x ) is the solution to the macroscopic problem (see Figure 2a,b):
{ ( a i j k l ( x ) u 0 k , l y ) , i y = q i ( x )   in   P 0 , u 0 ( y 1 0 , y 2 ) = u 0 ( y 1 + 0 , y 2 )   on   D ,   σ n ( y 1 0 , y 2 ) = σ n ( y 1 + 0 , y 2 )   on   D ,   u 0 ( x ) = 0   on   Q  
and corrector ε u 1 ( x / ε ) is the “local perturbations”, i.e., ε u 1 ( x / ε ) satisfied the following conditions:
u 1 ( y ) 0   as   | y 3 |   and   u 1 ( y ) 0   is   periodic   in   ( y 1 , y 2 )   with   PC   T
The local perturbation
  • decays in y 3 with distance from the weld D ε ;
  • oscillates in the variables ( y 1 , y 2 ) with PC T .
Substituting equality (3) into the elasticity theory equations, we arrive at the following periodicity cell problem (see for details [18,19]):
{ ( a i j k l ( y ) u 1 k , l y + a i j α β ( y ) u 0 α , β x ( x ) ) , i y = 0   in   P 0 , u 1 ( y )   is   periodic   in   y 1 , y 2   with   PC   T ,   u 1 ( y ) 0   as   | y 3 |
In (6), the local elastic constant is given by (2).
Problem (6) looks similar to the periodicity cell problem of the homogenization theory [20]. The difference is the decay condition: u 1 ( y ) 0 as | y 3 | in (6), while in the homogenization theory, the solution to the periodicity cell problem is periodic in all the variables [20]. This difference leads to significant differences between the limit problems. The first one is that the individual properties of the weld do not influence the macroscopic conditions on the weld surface, see for details [18,19]. The second difference is the discontinuity of several components of the macroscopic stress–strain state at the weld surface (the term u 0 α , β x ( x ) in (6) corresponds to the macroscopic stress–strain state). The “jump” of macroscopic deformations is seen in Figure 5a, Figure 6c and Figure 7c. As a result of this “jump”, separation of the “fast” and the “slow” variables becomes impossible. The essence of the separation of the “fast” and the “slow” variables is treating u 0 α , β x ( x ) as a parameter and “removing” it from the problem (6). This separation of the “fast” and the “slow” variables is the basis for the classical homogenization theory technique and leads to the periodicity cell problem arising in the homogenization theory [20]. In the case under consideration, we have to solve problem (6) in the original form. Note that a compatibility condition arises for the problem (6), see [18,19] for details. The existence of the compatibility condition also distinguishes (6) from the periodicity cell problem of the homogenization theory (it is known, see, e.g., [20], that there is no compatibility condition for the periodicity cell problem of the homogenization theory).

3. Numerical Computation of the Local Stresses in the Vicinity of the Waved Weld

As above, s , t mark the tangential and n marks the normal direction in the coordinates related to the welding seam. In Figure 1, ( s , t , n ) = ( x 1 , x 2 , x 3 ) .
It is important for the following that at the macroscopic interface surface D , see Figure 2c, the strains e s s ,   e t t ,   e s t and stresses σ n n ,   σ n s ,   σ n t are continuous [22]. The remaining strains and stresses may be discontinuous [22].
The total number of the continuous components of the macroscopic stress–strain state is six, as the number of stress (or strain) tensor components. Thus, one can express the microscopic stresses in the weld zone through the six continuous components e s s ,   e t t ,   e s t and σ n n ,   σ n s ,   σ n t of the macroscopic stress–strain tensors.
By using Hooke’s law and Formula (3), we write the local (microscopic) stresses σ i j l o c ( x ) in the following form:
σ i j l o c ( x ) = a i j k l ( x / ε ) ( u 0 ( x ) + ε u 1 ( x , x / ε ) ) k , l y a i j k l ( x / ε ) u 1 k , l y + a i j α β ( x / ε ) u 0 α , β x ( x ) .
As follows from [18,19], one can write the Formula (7) in the terms of the six continuous components of the macroscopic stress and strain tensors ( { σ n t } = { σ 33 ,   σ 31 ,   σ 32 } and { e s t } = { e 11 ,   e 22 ,   e 12 } ) in the form:
σ i j l o c ( x ) = K i j n n ( x / ε ) σ n n ( x ) + K i j n t ( x / ε ) σ n t ( x ) + K i j s t ( x / ε ) e s t ( x ) ,
where K i j n n ( y ) and K i j n t ( y ) are the functions similar to the local stress concentration tensors [30,31], K i j s t ( y ) are stresses in problem (6) with unit macroscopic strains. All these functions are expressed through the elastic stresses in (6).
Consider the planar problem. An argument in favor of such consideration is the geometry of the explosion weld. The geometry of the explosion weld changes significantly in two directions and changes slowly in the third direction. If we neglect slow changes in the third direction, we arrive at a problem that depends on two spatial variables. Such a three-dimensional problem of elasticity theory can be reduced to a planar problem. The planar periodicity cell problem corresponds to the planar strain state (i.e., the local strains ε x z ( y ) = ε y z ( y ) = 0 and ε z z ( y ) = 0 ). The local stresses σ x z ( y ) = σ y z ( y ) = 0 and σ z z ( y ) = ν ( y ) ( σ x x ( y ) + σ y y ( y ) ) .
Following, we present the result of our numerical calculations. The problem (6) was solved using the FEM program ANSYS [32]. In the computations, the welded material was as follows: on the left, there was copper: Young’s modulus E = 1 · 10 11 Pa, Poisson’s ratio ν = 0.35, tensile strength σ C u 0 = 180 × 10 6 Pa; on the right there was aluminum: Young’s modulus E = 0.5 · 10 11 , Poisson’s ratio ν = 0.34, tensile strength σ A l 0 = 100 × 10 6 Pa. We carried out the numerical computations for the three typical shapes of the weld wave: symmetric sinusoidal-like wave (Figure 1a), asymmetric wave (Figure 1b) and wave with a crest (Figure 1c). The weld geometries were constructed based on the measurement of the photograph shown in Figure 1a–c. The computations were done for three basic planar strains: the tensions along O x - and O y -axes and the shift in O x y -plane. Note that not all of these strains were continuous.
Tension along O x -axis. We considered three widely observed experimental types of weld waves: symmetric wave, Figure 2a, asymmetric wave, Figure 2b, and wave with crest, Figure 3c. The results of the numerical computations for the three types of weld waves subjected to the tension along O x -axis are presented in Figure 3. Figure 3 displays the actual von Mises stress in the periodicity cells.
Tension along O y -axis. Figure 4 displays the actual von Mises stress in the periodicity cell for three types of the weld waves.
Shift in O x y -plane. Figure 5 displays the actual von Mises stress in the periodicity cell for three types of the weld waves.

4. The Problem of the Weld Strength

We assumed that the strength criteria for the connected materials had the form f i ( σ i j l o c ) < σ i * (the index i being the index of the material). After the problem (6) was solved, one could verify the strength criterion at every point of the connected materials. This approach correlated with the idea of the homogenization theory and with the approach used in the sub-structural analysis [33,34].
Two approaches (densely related one to another) are possible. The first is the computation of the local stresses in the vicinity of the weld, and the following fracture analysis. This is the problem related to the geometry of the weld. This approach is densely connected with the idea of the local stress concentration in inhomogeneous materials and structures [35]. The explosion not only forms a weld of complex geometry, but also changes the properties of the welded materials. It is possible to separate, at least theoretically, the problems related to the geometry of the weld and the problems of material properties in the weld. The second approach assumes the construction of the so-called homogenized strength criterion (if it exists for the problem under consideration).
We present solutions to the mentioned two problems for the planar problem, where the numerical computations are not as massive as in three-dimensional problems. In the planar problems, the interface strength criterion depends on following three macroscopic continuous stress–strain components: the macroscopic stresses σ n n ,   σ n s and the macroscopic strain ε s s . In both approaches, it is necessary to calculate the functions, K i j n n ( y ) , K i j n 3 ( y ) and K i j 33 ( y ) (the direction normal to the weld is x 3 n ).
The functions K i j n n ( y ) , K i j n 3 ( y ) and K i j 33 ( y ) ( i , j = 1 , 2 ) were obtained by dividing the stresses σ i j n n ( y ) , σ i j n 3 ( y ) , σ i j 33 ( y ) by the macroscopic stresses (for σ i j n n ( y ) and σ i j n 3 ( y ) ) or macroscopic strain (for σ i j 33 ( y ) ) at the long distance from the weld. At the distance “far from the weld”, the microscopic stresses/strains coincide with macroscopic stresses/strains. In accordance with our computations, the “far from the weld” distance was 3–5 widths of PC. Thus, all the data necessary for computation K i j n n ( y ) , K i j n 3 ( y ) and K i j 33 ( y ) are determined from the solution to the problem (6). Note that, although the microscopic and macroscopic stresses/strains coincide, the values may be different on the left and on the right from the weld. For this reason, K i j n 3 ( y ) are not exactly the classical stress concentration tensors.

4.1. Stress Concentrations in the Weld Zone

When the periodicity cell is subjected to tension along O x -axis or shift, the stresses to the left and right of the weld are continuous and the macroscopic deformations are discontinuous. When the periodicity cell is subjected to tension along the O y -axis, the stresses to the left and right of the weld are discontinuous and the macroscopic deformation is continuous. This situation shows a similarity with the Voigt- and Reiss-type deformations arising in a layered material. In the case under consideration, the deformations in the materials were localized near the weld. In the case of macroscopic tension along O x -axis or shift, the von Mises stress concentration in the weld zone is calculated according to the formula:
k ( x ) = σ M ( x ) / σ M ( ) ,
In the case of the macroscopic tension along the O y -axis, the von Mises stress concentration is determined as follows:
  • to the left of the weld, k ( x ) = σ M ( x ) / σ M | Z 1 = ;
  • to the right of the weld, k ( x ) = σ M ( x ) / σ M | Z 1 = + .
Here σ M means the stress intensity, and σ M | Z 1 = ± are stresses to the left and right of the weld, at the distances “far from the weld”. According to the results of numerical calculations from [17], these distances are ( 3 5 ) T ( T is the characteristic dimension of the weld waves in the “fast” variables, see Figure 1).
The maximum Mises stress concentration is determined as k = max x   k ( x ) .
Following, we present the results of our numerical computations. Figure 6 displays the von Mises stress concentrations in the periodicity cell corresponding to the symmetric weld.
Figure 7 displays the von Mises stress concentrations in the periodicity cell corresponding to the asymmetric weld.
Figure 8 displays the von Mises stress concentrations in the periodicity cell corresponding to the weld with crest.
We collected the maximum computed von Mises stress concentrations in Table 1.

4.2. Constructing the Macroscopic Strength Criterion of the Weld

In this section, we discuss the following question: “Is the weld an independent object with its own strength characteristics, or is it a sub-structural fragment of its encompassing structure without its own physical properties?” Note that in the classical elasticity theory, the interface does not form an independent object. Although the interface conditions in the classical elasticity theory are written in terms of the characteristics of the connected bodies [22], the interface is not a sub-structural fragment of the connected bodies.
In this section, we present an answer to the question under the assumption that the local strength limits are constant over the connected material. This assumption does not take into account the phenomenon of strain hardening. We discuss accounting for the strain hardening below. Note that the analysis presented in the previous section provided us with the stress concentration coefficients in the weld and did not involve the strength of the materials.
Let the strength criteria of the welded materials have the form
f ( σ i j ) < σ * ( y ) ,   f + ( σ i j ) < σ * + ( y )
(“−” marks the material to the left of the weld, “+” marks the material to the right of the weld).
The strength criterion in the arbitrary point y may be written in the form:
f ( y , σ i j l o c ( y ) ) < σ * ( y )
where the local stress limit is:
σ * ( y ) = { σ *   in   P , σ * +   in   P +
and the function is:
f ( y , σ i j ) = { f ( σ i j )   in   P , f + ( σ i j )   in   P +
Substituting σ i j l o c ( y ) following (8) into (9), we obtain:
f ( y , K i j n t ( x / ε ) σ n t ( x ) + K i j s t ( x / ε ) e s t ( x ) ) < σ * ( y )
As explained in [35], there is no failure at any points of the PC if the following inequality is satisfied:
F ( σ n t ( x ) , e s t ( x ) ) max y P f ( x / ε , K i j n t ( x / ε ) σ n t ( x ) + K i j s t ( x / ε ) e s t ( x ) ) σ * ( y ) < 1 .
The failure starts when the condition F ( σ n t ( x ) , e s t ( x ) ) = 1 is satisfied. The failure occurs at the points y 0 at which:
f ( y 0 , K i j n t ( y 0 ) σ n t ( x ) + K i j s t ( y 0 ) e s t ( x ) ) = σ * ( y 0 ) .
Both the “no failure” criteria (12) and the “first failure” condition F ( σ n t ( x ) , e s t ( x ) ) = 1 are written in terms of macroscopic strains e s s ,   e t t ,   e s t and macroscopic stresses σ n n ,   σ n s ,   σ n t . This means that the “no failure” and the “first failure” criteria at the micro-level can be written in terms of the macroscopic stress–strain state.
When considering a planar problem, the macroscopic strength criterion of the interface can be written in the terms of two macroscopic stresses σ n n ,   σ n t and one strain ε t t . It has the form F ( σ n t ( x ) , e s t ( x ) ) = 1 .
In our computations, we used the von Mises strength criterion at the micro-level. This strength criterion is often used for plastic materials, in particular, metals (see, e.g., [36]). Then f ( y , σ i j l o c ( y ) ) = σ M i s e s l o c ( y ) .
The macroscopic strength criterion of the weld can be constructed in the following steps:
1. Compute K i j n n ( y ) , K i j n 3 ( y ) and K i j 33 ( y ) , and save the computed values.
2. Form the linear combinations:
σ i j ( y ) = K i j 33 ( y ) ε x x + K i j n n ( y ) σ y y + K i j n 3 ( y ) σ x y .
and select from the combinations in (14) ones for which the condition in (15) is satisfied:
max y P σ M i s e s ( y ) σ * ( y ) = 1   ( or   1 ) .
In (15), σ M i s e s ( y ) means von Mises stress corresponding to the local stresses (14).
Step 1 may be done with ANSYS or any similar FEM software. Step 2 requires additional programming. The program must process the ANSYS data files and make the computations indicated in (14) and (15).
The coefficients of the linear combinations (14) satisfying equality (15) form the macroscopic failure surface of the weld.
One can make calculations (14) and (15) in each of the welded bodies: in the left domain p and in the right domain p + (see Figure 1c) separately and compute σ = max y P σ M i s e s l o c ( y ) σ * ( y ) and σ + = max y P + σ M i s e s l o c ( y ) σ * ( y ) . The strength criterion (15) may be written in the form max ( σ M i s e s , σ M i s e s + ) 1 . The failure starts at max ( σ M i s e s , σ M i s e s + ) = 1 . This equality means σ = 1 or σ + = 1 . Thus, it is possible to indicate in which materials the failure starts (in other words, it is possible to indicate the “weakest” material).
The authors performed computations described in item 1 by using ANSYS FEM software and developed a computer program in order to perform the computations mentioned in item 2.
Biaxial tension. Consider the case of biaxial tension. In this case the macroscopic shift stress σ 12 = 0 , and the strength criterion of the weld is written in terms of the continuous macroscopic axial tension ε 22 and axial stress σ 11 . The failure surface is a curve in plane O σ x x ε y y . It is suitable to express the strength criteria of the welded materials, aluminum and copper, also in the terms of σ 11 and ε 22 . The domains S A l = { ( σ 11 , ε 22 ) : σ M i s e s ( σ 11 , e 22 ) σ A l * } and S C u = σ M i s e s ( σ 11 , e 22 ) σ C u * display the failure surfaces of pure aluminum and copper, correspondingly. The domains S = { ( σ 11 , ε 22 ) : max σ M i s e s l o c ( y ) / y p σ A l * 1 } and S + = { ( σ 11 , ε 22 ) : max σ M i s e s l o c ( y ) / y p + σ C u * 1 } display the copper and the aluminum failure surfaces in the vicinity of the weld (copper is to the left of the weld, aluminum is to the right of the weld).
The failure surface of the weld is the intersection S = S S + . We developed computer programs for the computation of the domains S C u , S A l , S , S + and S . Figure 9 displays all the mentioned domains in the plane O σ 11 ε 22 .
It is seen in Figure 9a,b that aluminum and copper as the elements of the weld were less strong than the pure aluminum and copper, correspondingly. This was the result of the stress concentration in the weld zone, see Figure 6, Figure 7 and Figure 8. At the same time, we observed no direct relationship between the strength of materials in the weld zone and the strength of the pure materials. In our computations, copper was approximately twice as strong as the aluminum: σ C u * = 180 × 10 6 Pa and σ A l * = 100 × 10 6 Pa. In the weld, on the top and the bottom of S , copper was the strongest element of the weld. On the left and the right sides of S , aluminum was the strongest element of the weld, see Figure 9c.
Planar deformation of general form. In the general case, σ x x ,   ε y y ,   σ x y take non-zero values, thus, the failure surfaces are three-dimensional objects in the coordinate system O σ 11 ε 22 σ 12 . Denote S A l = { ( σ 11 , ε 22 , σ 12 ) : σ M i s e s ( σ 11 , ε 22 , σ 12 ) σ A l * }, S C u = { ( σ 11 , ε 22 , σ 12 ) : σ M i s e s ( σ 11 , ε 22 , σ 12 ) σ C u * }, S = { ( σ 11 , ε 22 , σ 12 ) : max y p + σ M i s e s l o c ( y ) / σ C u * 1 } and S + = { ( σ 11 , ε 22 , σ 12 ) : max y p σ M i s e s l o c ( y ) / σ A l * 1 }. The domains S A l and S C u are displayed in Figure 10. The surfaces are colored to improve the visibility of these objects.
The domains S and S + are displayed in Figure 11a,b. The weld failure surface S is the intersection of S and S + : S = S S + .
In the case under consideration, a large part of S lay inside S + , but not all S . The domain S cut off the left and right parts of S + . The remaining part of S was the S = S S + . It is displayed in Figure 11c. The failure started on the left (copper) or the right (aluminum) side of the weld. It depended on the macroscopic stress–strain state of the weld.
The failure surfaces S A l and S C u are compared in Figure 10 and the weld failure surface S = S S + in Figure 11b. This comparison made it possible to estimate the decrease in the strength of the weld compared with the strength of the welded materials numerically. In the case displayed in Figure 10 and Figure 11, the strength of the weld was between two to four times less than the strength of the welded materials.

5. Prospective

The effect of the microstructure on the mechanical properties of explosive welding joints is an actual problem nowadays, see, e.g., [37,38,39], and references in these papers. The problem requires both experimental and theoretical research. The number of experimental researches prevail over the number of theoretical researches. In the opinion of the authors, the two-scale method of local perturbations, developed in [18,19] and partially inspired by the homogenization method in composite material [20,27,28,29,30,31], is an adequate tool for both the theoretical analysis and numerical computation of the effect of the microstructure on the mechanical properties of explosive welding joints.
Above, we presented the theory of the two-scale method, as applied to weld joints, and demonstrated the capabilities of this method for some specific problems. Now, we discuss several prospective problems. The prospects for research in the field of strength of welds based on the analysis of the microscopic stress–strain state in them, by the methods of the mechanics of solids, are very diverse. Below, we discuss several prospective problems directly related to our research.

5.1. Accounting for Strain Hardening in a Weld

The photos of explosion weld in Figure 1 lead to the conclusion that during welding, significant deformations occur in the region of wave formation of the weld, and fade at a distance equal to one weld wave amplitude. Similar information may be found in [11,12,13,14].
The plastic deformations were non-zero in the weld zone and zero at a distance from the weld zone. This meant that the plastic deformations were functions of the spatial variable y . Correspondingly, the strength limit of the welded materials was the function of the spatial variable y and had the following form:
σ * ( y ) = { σ *   in   P , σ * 0 ( y )   in   P 0 , σ * +   in   P +
Replacing (10) with (16) does not bring a significant change to the mathematical algorithm described above. As for the stress concentrations, this replacement does not affect the calculations related to Section 4.2 at all.
The only problem is obtaining the function σ * 0 ( y ) , which describes the strength of materials after strain hardening. If the function σ * 0 ( y ) is known, we arrive at the problem, which is reminiscent of the problems arising in connection with the graded materials discussed in [40,41,42,43]. One can construct an approximation of the function, σ * 0 ( y ) , if one knows the characteristic magnitude of the plastic deformation during the welding process and the form of the functions in (16). Little is known about the exact values of these parameters. One can estimate the deformations in the weld zone by using the photographs of the weld zone. Regarding the form of the functions in (16), the papers on the graded material often assume that the distribution may be described by power functions, see, e.g., [40,41,42,43,44].

5.2. Accounting for Interlayer in the Weld Zone

Between the welded metals there can be an interlayer which is both artificial (for example, an additional thin technological layer) and natural (for example, the “mixing zone” that occurs during the welding, Figure 12a).
This interlayer can be taken into account within the framework of the local perturbation approach if we introduce the interlayer region, P i and define the elastic constants and ultimate strength as follows:
a i j k l ( y ) = { a i j k l   in   P , a i j k l   in   P i , a i j k l +   in   P + . σ * ( y ) = { σ * ( y )   in   P , σ * i   in   P i , σ * + ( y )   in   P + ,
The layer shown in Figure 12b looks like a typical composite material of random structure. So, then you can use the homogenization method to calculate (estimate) the elasticity constants a i j k l   in   P i and the strength limit σ * i   in   P i . Undoubtedly, this is a non-trivial problem, since the material in the domain P i has a random structure and its components subjected to large plastic deformations. However, such problems are being studied, and some results can be found in the literature on composite materials, see [45,46] for homogenization of random composites and [47] for homogenization in a plastic problem, see also reference in [45,46].

6. Conclusions

We presented the theory of the two-scale method, as applied to weld joints, and used this method for construction of the strength criterion of the joint. We applied the developed method to the computation of the local stress concentrations in the explosion weld. In accordance with our computation, the stress concentration in the explosion welding zone changed from 1.2 to 1.7, depending on the weld wave geometry and the macroscopic deformation. If strain hardening exceeded stress concentration, the strength of the explosion weld metals exceeded the strength of the weakest material of the pair and destruction occurred at a certain distance from the plane. In accordance with our considerations, the destruction did not necessarily occur along the least durable metal. In accordance with our considerations, the weakest point is determined not by the strength limit (the strength limit after the strain hardening), but by the ratio:
strength   limit   after   the   strain   hardening stress   concentration   ratio
At the same time, if the strain hardening in the weld was small. compared with the stress concentration in the weld, the destruction occurred along the weld.
In any case, at the macrolevel, the weld may be treated as an independent solid mechanical object, possessing its own macroscopic strength criterion.

Author Contributions

Conceptualization, methodology, validation, formal analysis, writing—original draft preparation, writing—review and editing, visualization, supervision, A.G.K.; data curation, writing—original draft preparation, visualization, project administration, funding, S.I.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available by request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carl, L.R. Brass welds made by detonation impulse. Metal Progr. 1944, 46, 102–103. [Google Scholar]
  2. Walsh, J.M.; Shreffler, R.G.; Willig, E.J. Limiting conditions for jet formation in high velocity collisions. J. Appl. Phys. 1953, 24, 349–359. [Google Scholar] [CrossRef] [Green Version]
  3. Allen, W.; Mapes, J.; Wilson, G. An effect produced by oblique impact of a cylinder on a thin target. J. Appl. Phys. 1954, 25, 675–676. [Google Scholar] [CrossRef]
  4. Pearson, J. Recent Advances in Explosive Pressing and Welding. In Proceedings of the 2nd Metals Energy Conference of Explosive, Chicago, IL, USA, 3 November 1959; 1959; pp. 32–36. [Google Scholar]
  5. Pearson, J. Metal working with explosives. JOM 1960, 12, 673–681. [Google Scholar] [CrossRef]
  6. Cowan, G.; Douglas, J.; Holtzman, A.H. Explosive Bonding. U.S. Patent 3,137,937, 26 October 1960. [Google Scholar]
  7. Philipchuk, V.; Le Roy, B.F. Explosive Welding. U.S. Patent 3,024,526, 31 August 1960. [Google Scholar]
  8. Philipchuk, V. Explosive welding status. In Creative Manufacturing Seminar; Techn. Paper; American Society of Tool & Manufacturing Engineers: Dearborn, MI, USA, 1965; pp. 65–109. [Google Scholar]
  9. Sedikh, V.S.; Deribas, A.A.; Bichenkov, E.I.; Trishin, Y.A. Welding by explosion. Weld. Prod. 1962, 5, 3–6. [Google Scholar]
  10. Bayce, A. Explosive welding. Stanf. Res. Inst. J. 1967, 13, 2–8. [Google Scholar]
  11. Deribas, A.A.; Kudinov, V.M.; Matveenkov, F.I.; Simonov, V.A. Explosive welding. Combust. Explos. Shock Waves 1967, 3, 69–72. [Google Scholar] [CrossRef]
  12. Available online: http:www.hydro.nsc.ru/technology/techn1.php (accessed on 6 October 2022).
  13. Degarmo, E.P.; Black, J.T.; Kohser, R.A. Materials and Processes in Manufacturing, 9th ed.; Wiley: Chichester, UK, 2003. [Google Scholar]
  14. GOST 5639-82; Steels and Alloys. Methods Detecting and Identifying Grain Size. GOST: Moscow, Russia, 1983.
  15. Smith, W.F.; Javad, H. Foundations of Materials Science and Engineering; McGraw-Hill: New York, NY, USA, 2006. [Google Scholar]
  16. Maliutina, I.N.; Stepanova, N.V.; Cherkov, A.G.; Chuchkova, L.V. Welding of dissimilar materials with interlayer employment containing copper and tantalum. Obrab. Met. 2015, 4, 61–71. [Google Scholar]
  17. Kolpakov, A.G.; Rakin, S.I. Estimation of stress concentration in a welded joint formed by explosive welding. J. Appl. Mech. Technol. Phys. 2018, 59, 569–575. [Google Scholar] [CrossRef]
  18. Kolpakov, A.G.; Andrianov, I.V.; Rakin, S.I.; Rogerson, G.A. An asymptotic strategy to couple homogenized elastic structures. Int. J. Eng. Sci. 2018, 131, 26–39. [Google Scholar] [CrossRef]
  19. Kolpakov, A.G.; Andrianov, I.V.; Prikazchikov, D.A. Asymptotic strategy for matching homogenized structures. Conductivity problem. Q. J. Mech. Appl. Math. 2018, 71, 519–535. [Google Scholar] [CrossRef]
  20. Sanchez-Palencia, E. Non-Homogeneous Media and Vibration Theory; Springer: Berlin/Heidelberg, Germany, 1980. [Google Scholar]
  21. Born, M.; Huang, K. Dynamical Theory of Crystal Lattices; Clarendon Press: Oxford, UK, 1954. [Google Scholar]
  22. Timoshenko, S.; Goodier, J.N. Theory of Elasticity; McGraw-Hill: New York, NY, USA, 1951. [Google Scholar]
  23. Wagner, G.J.; Liu, W.K. Coupling of atomistic and continuum simulations using a bridging scale decomposition. J. Comput. Phys. 2003, 190, 249–274. [Google Scholar] [CrossRef]
  24. Xiao, S.P.; Belytschko, T. A bridging domain method for coupling continua with molecular dynamics. Comput. Methods Appl. Mech. Eng. 2004, 193, 1645–1669. [Google Scholar] [CrossRef]
  25. Parks, M.L.; Bache, P.B.; Lehoucq, R.B. Connecting atomistic-to-continuum coupling and domain decomposition. Multiscale Model. Simul. 2008, 7, 362–380. [Google Scholar] [CrossRef]
  26. Pfaller, S.; Possart, G.; Steinmann, P.; Rahimi, M.; Müller-Plathe, F.; Böhm, M.C. A comparison of staggered solution schemes for coupled particle-continuum systems modeled with the Arlequin method. Comp. Mech. 2012, 49, 565–579. [Google Scholar] [CrossRef]
  27. Oleynik, O.A.; Shamaev, A.S.; Yosifian, G.A. Mathematical Problems in Elasticity and Homogenization; North-Holland: Amsterdam, The Netherlands, 1992. [Google Scholar]
  28. Bakhvalov, N.S.; Panasenko, G.P. Averaging Processes in Periodic Media. Mathematical Problems in Mechanics of Composite Materials; Kluwer: Dordrecht, The Netherlands, 1989. [Google Scholar]
  29. Sanchez-Palencia, E. Boundary Layers and Edge Effects in Composites: Homogenization Techniques for Composite Materials; Sanchez-Palencia, E., Zaoui, A., Eds.; Springer: Berlin/Heidelberg, Germany, 1987; pp. 122–193. [Google Scholar]
  30. Dvorak, G. Micromechanics of Composite Materials; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  31. Buryachenko, V.A. Micromechanics of Heterogeneous Materials; Springer: New York, NY, USA, 2007. [Google Scholar]
  32. ANSYS Inc. ANSYS 5.7 Structural Analysis Guide; ANSYS Inc.: Tokyo, Japan, 2000. [Google Scholar]
  33. Vankan, W.J.; van den Brinkm, W.M.; Maas, R. Multi-Level Structural Analysis for Sub-Component Validation in Aircraft Composite Fuselage Structures; NLRNLR-TP-2015-275; Nederlands Lucht- en Ruimtevaartcentrum: Amsterdam, The Netherlands, 2015. [Google Scholar]
  34. Zhou, K.; Cashany, M.; Wu, Z.Y. Substructure Analysis Framework for Accelerated Finite Element Modeling; Technical Report for Bentley Systems; Bentley Systems: Watertown, CT, USA, 2015. [Google Scholar]
  35. Kalamkarov, A.L.; Kolpakov, A.G. Analysis, Design and Optimization of Composite Structures; Wiley: Chichester, UK, 1997. [Google Scholar]
  36. Yu, M.-H. Unified Strength Theory and Its Applications; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  37. Wang, K.; Kuroda, M.; Chen, X.; Hokamoto, K.; Li, X.; Zeng, X.; Nie, S.; Wang, Y. Mechanical Properties of Explosion-Welded Titanium/Duplex Stainless Steel under Different Energetic Conditions. Metals 2022, 12, 1354. [Google Scholar] [CrossRef]
  38. Mahmood, Y.; Dai, K.; Chen, P.; Zhou, Q.; Bhatti, A.A.; Arab, A. Experimental and numerical study on microstructure and mechanical properties of Ti-6Al-4V/Al-1060 Explosive Welding. Metals 2019, 9, 1189. [Google Scholar] [CrossRef] [Green Version]
  39. Zhou, Q.; Liu, R.; Zhou, Q.; Ran, C.; Fan, K.; Xie, J.; Chen, P. Effect of microstructure on mechanical properties of titanium-steel explosive welding interface. Mater. Sci. Eng. A 2022, 830, 142260. [Google Scholar] [CrossRef]
  40. Erdogan, F. Fracture mechanics of functionally graded materials. Compos. Eng. 1995, 5, 753–770. [Google Scholar] [CrossRef]
  41. Tan, C.; Chew, Y.; Bi, G.; Wang, D.; Ma, W.; Yang, Y.; Zhou, K. Additive manufacturing of steel–copper functionally graded material with ultrahigh bonding strength. J. Mater. Sci. Technol. 2021, 72, 217–222. [Google Scholar]
  42. Kolednik, O. The yield stress gradient effect in inhomogeneous materials. Int. J. Solids Struct. 2000, 37, 781–808. [Google Scholar] [CrossRef]
  43. Angshuman, C.; Gopinath, M.; Sagar, S.; Vikranth, R.; Ashish Kumar, N. Mitigation of cracks in laser welding of titanium and stainless steel by in-situ nickel interlayer deposition. J. Mater. Process. Technol. 2022, 300, 117403. [Google Scholar]
  44. Baytak, T.; Bulut, O. Thermal stress in functionally graded plates with a gradation of the coefficient of thermal expansion only. Exp. Mech. 2022, 62, 655–666. [Google Scholar] [CrossRef]
  45. Mityushev, V.; Rylko, N. Effective properties of two-dimensional dispersed composites. Part I. Schwarz’s alternating method. Comput. Math. Appl. 2022, 111, 50–60. [Google Scholar] [CrossRef]
  46. Mityushev, V. Effective properties of two-dimensional dispersed composites. Part II. Revision of self-consistent methods. Comput. Math. Appl. 2022, 121, 74–84. [Google Scholar] [CrossRef]
  47. Xia, D.; Zhang, X.; Oskay, C. Large-deformation reduced order homogenization of polycrystalline materials. Comput. Methods Appl. Mech. Eng. 2021, 387, 114–119. [Google Scholar] [CrossRef]
Figure 1. The explosion weld: (b) macroscopic view, (c) microscopic view, (a,d) models (macroscopic and microscopic).
Figure 1. The explosion weld: (b) macroscopic view, (c) microscopic view, (a,d) models (macroscopic and microscopic).
Materials 15 07878 g001
Figure 2. View of the weld at a microlevel: (a) symmetric wave; (b,d) asymmetric wave; (c,e,f) wave with a crest.
Figure 2. View of the weld at a microlevel: (a) symmetric wave; (b,d) asymmetric wave; (c,e,f) wave with a crest.
Materials 15 07878 g002
Figure 3. The local von Mises stresses corresponding to the tension along O x -axis for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Figure 3. The local von Mises stresses corresponding to the tension along O x -axis for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Materials 15 07878 g003
Figure 4. The local von Mises stresses corresponding to the tension along Oy-axis for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Figure 4. The local von Mises stresses corresponding to the tension along Oy-axis for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Materials 15 07878 g004
Figure 5. The local von Mises stresses corresponding to the shift for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Figure 5. The local von Mises stresses corresponding to the shift for (a) symmetric wave, (b) asymmetric wave, (c) wave with crest, Figure 3c.
Materials 15 07878 g005
Figure 6. The von Mises stress concentration in the vicinity of symmetric weld: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Figure 6. The von Mises stress concentration in the vicinity of symmetric weld: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Materials 15 07878 g006
Figure 7. The von Mises stress concentration in the vicinity of asymmetric weld: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Figure 7. The von Mises stress concentration in the vicinity of asymmetric weld: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Materials 15 07878 g007
Figure 8. The von Mises stress concentration in the vicinity of the weld with crest: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Figure 8. The von Mises stress concentration in the vicinity of the weld with crest: (a) tension along O x -axis, (b) tension along O y -axis, (c) shift.
Materials 15 07878 g008
Figure 9. The failure curves for biaxial tension ((a) cooper, (b) aluminum, (c) weld).
Figure 9. The failure curves for biaxial tension ((a) cooper, (b) aluminum, (c) weld).
Materials 15 07878 g009
Figure 10. Failure surfaces: (a) pure copper, (b) pure aluminum.
Figure 10. Failure surfaces: (a) pure copper, (b) pure aluminum.
Materials 15 07878 g010
Figure 11. Failure surfaces: (a) for the left (copper) and (c) the right (aluminum) sides of the weld separately; (b) for the weld as a whole, (a,b) S and S + together.
Figure 11. Failure surfaces: (a) for the left (copper) and (c) the right (aluminum) sides of the weld separately; (b) for the weld as a whole, (a,b) S and S + together.
Materials 15 07878 g011
Figure 12. The interlayer in the weld zone melted zone ((a) view at the distance, (b) zoomed).
Figure 12. The interlayer in the weld zone melted zone ((a) view at the distance, (b) zoomed).
Materials 15 07878 g012
Table 1. Maximum von Mises stress concentration k for different welds and different macroscopic strains.
Table 1. Maximum von Mises stress concentration k for different welds and different macroscopic strains.
Weld Typek
Tension along the Ox-AxisTension along the Oy-AxisShift in the
Oxy-Plane
Symmetric wave1.211.491.25
Asymmetric wave1.311.541.28
Wave with a crest1.141.701.20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kolpakov, A.G.; Rakin, S.I. A New Microstructural Approach to the Strength of an Explosion Weld. Materials 2022, 15, 7878. https://doi.org/10.3390/ma15227878

AMA Style

Kolpakov AG, Rakin SI. A New Microstructural Approach to the Strength of an Explosion Weld. Materials. 2022; 15(22):7878. https://doi.org/10.3390/ma15227878

Chicago/Turabian Style

Kolpakov, Alexander G., and Sergei I. Rakin. 2022. "A New Microstructural Approach to the Strength of an Explosion Weld" Materials 15, no. 22: 7878. https://doi.org/10.3390/ma15227878

APA Style

Kolpakov, A. G., & Rakin, S. I. (2022). A New Microstructural Approach to the Strength of an Explosion Weld. Materials, 15(22), 7878. https://doi.org/10.3390/ma15227878

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