# The Effect of Compression on the Void Coalescence under Strong Dynamic Loading

Authors: Ruichen Tu; Ning Wei; Yongmao Pei; Yan Liu; Fengguo Zhang; Dongsheng Zhang
Journal: Advances in Materials Science and Engineering
 (2022)

---

## Abstract

The void coalescence under strong dynamic loading is a significant spallation process for ductile metals. Since the spallation is basically dominated by tension waves, most void coalescence studies have focused on the tension effect. However, it is known that in spallation, the material initially undergoes a strong compression wave, and then an irreversible deformation is produced by the compression wave inside the material. Therefore, in this paper, the effect of compression on the void coalescence is investigated using the molecular dynamics (MD) simulation. It was found that as the compressive strain increases, the yield strength decreases first and then increases. The results showed that due to the Bauschinger effect (BE), the yield strength decreases by 19.43% from 5.66 GPa without compressive loading to 4.56 GPa when the compressive strain is −7.5%, after which the yield strength increases. The voids do not coalesce when the compressive strain is −8%. In addition, it was found that during the compressive phase, the void surfaces would generate dislocations, which could obstruct the void coalescence in the tensile phase. Furthermore, under compressive loading, the temperature effect on the void coalescence was studied, and it was found that lower temperatures could suppress the void coalescence.


---

## Body

## 1. Introduction

Dynamic tensile fracture (spallation) of ductile metals under strong dynamic loading is a typical failure mode of material dynamic damage, which is also a fundamental problem in engineering. Under strong dynamic loading, the shock wave interaction in opposite directions generates a stress pulse in ductile metals, which induces microvoid nucleation, growth, and coalescence corresponding to the accumulation, growth, and rapid development of damage. Once the damage has accumulated to a certain extent, it causes material failures. As the rapid growth of damage is a critical process which causes material failure, the mechanism of the void coalescence needs to be further studied.The relations between the void coalescence and void sizes, the mutual position of voids, loading conditions, and material properties have been studied in many references. Tvergaard et al. [1, 2] introduced two adjustment parameters to consider the interaction between voids. A model, which is known as the Gurson–Tvergaard–Needleman (GTN) model [3, 4], has taken into account the process of rapid loss of load carrying capacity caused by the void coalescence. Gologanu et al. [5, 6] developed a model named Gologanu–Leblond–Devaux (GLD), which considered the effect of void shape evolutions on the void coalescence. Earlier studies often focused on the void coalescence under quasistatic loading, which led researchers to study the case under strong dynamic loading. In order to study dynamic tensile fracture, the split Hopkinson pressure bar (SHPB) loading technology [7] has been widely used for problems with a strain rate of less than 104 s−1 and the laser loading technology [8] for higher strain rate conditions between 105–107 s−1. A model was established by Thomason et al. [9], and they found that the effects of inertia and strain rate influenced the growth of dynamic damage. The results of this model were in good agreement with the relevant experimental results. Moreover, temperature has been found as an important parameter. For instance, the effect of temperature on void evolutions has been confirmed [9–11]; it significantly influenced the critical stress and damage growth rate, and the fracture strength had a negative correlation with the initial temperature. In addition, it is worth mentioning that the distance between two voids and the void radius [12] can also greatly affect the coalescence efficiency. A study [13] found that there was a relationship between the intervoid ligament distance (ILD) of two same voids. When the ILD was equal to the void radius, the voids were most likely to coalesce.Due to experimental and observational limits, the internal evolution of void coalescence in extreme conditions has not been universally investigated. Therefore, numerical simulation methods are used to study the void coalescence under strong dynamic loading. Jacques et al. [14] revealed the performance difference of samples under quasistatic and strong dynamic loading. The inertia effect of ductile metals subjected to dynamic loads was also studied [15], indicating that this effect on the growth and coalescence of voids is significant. Lubarda et al. [16] pointed out that at higher strain rates, the damage could only grow through surface emission dislocations, and the critical stress of the dislocation gradually decreased with increasing damage size. Horstemeyer et al. [17] proposed a critical ILD to define the distance between two voids, and they employed the micromechanical finite element method (FEM) to quantify the coalescence effect based on temperature and voids that had different spatial arrangements. It is known that the molecular dynamics (MD) method can properly simulate the processes of void nucleation, growth, and coalescence, so that it has gradually become popular for investigating void coalescence. Potirniche et al. [18] proved that the void size and strain rate of loading influence the void growth and coalescence of single crystal nickel. Tang et al. [19] by comparing the metals of face-centered cubic (FCC) crystal with hexagonal close-packed (HCP) crystal, found that the void coalescence mode of the HCP crystal largely depended on crystal orientations. Cui et al. [20] considered the void shape in coalescence and confirmed that the elliptical voids with an oblate-oblate combination would be easier to coalesce. Deng et al. [21] studied the effect of different loading directions on the void coalescence based on the [100] crystal orientation. The results illustrated that the voids were most likely to coalesce when the loading direction was at an angle of 60° to the centerline of the two voids. Zhang et al. [22] suggested that even if the dislocation of the twin crystal was not as slow as that of a single crystal, the coalescence velocity of the voids of a twin crystal is still faster than that of a single crystal. The twin boundary can reduce the critical stress and suppress dislocations. Peng et al. [23] performed tens of millions of atoms of MD simulations, in which the voids collapsed under an impact loading, while after a tensile loading was applied, the voids were re-initiated during the contraction. Diwu et al. [24] studied the compression of porous voids under high strain rates and revealed that the porous structure under compressive loading would cause the elastic modulus and rheological peak stress of FCC metals to decrease. Liu et al. [25] investigated the effects of void radius and ILD on the void coalescence, and it was found that the rapid growth stage for the void growth was dependent on the initial void radius, and smaller radii can lead to higher growth rates. The first peak value of stress triaxiality increases first and then decreases with the increase of ILD.However, in spallation, the compression wave is transformed into a tensile wave on the free surface, where it encounters the incident unloading wave and a local tensile stress region is generated. Consequently, the generation of shock waves is often accompanied by compressive and tensile phases. Note that the void coalescence studied using the compressive-tensile loading method is rarely found in the reference. Therefore, in this study, the MD method has been utilized to simulate single-crystal aluminum under uniaxial compressive-tensile plastic deformation in the presence of nanovoids along the [100] crystal orientation. In addition, the effects of temperature and strain rate are also considered in studying the influence of compression on the void coalescence. In this regard, the yield strength of specimens under different compressive strains is also compared.


## 2. Simulation Method

In this study, the large-scale atomic/molecular massively parallel simulator (LAMMPS) was used in the MD simulation [26]. Since the accuracy of the interaction potential function between molecules or atoms has a tremendous influence on the MD method, it is extremely important to choose a suitable potential function. In this simulation, the interatomic interactions of aluminum were described by the embedded atom method (EAM) potential [27, 28].First, the atomic models for single-crystal aluminum were constructed. In order to eliminate the influence of box size as much as possible, some tensile simulations of the model without voids were set at first. The stress-strain curves of these simulations are shown in Figure1. When the model size is diminutive, the stress-strain curve oscillates severely, while this curve tends to be stabilized when the size is gradually increased to 30a, in which “a” is a unit cell. The value of “a” equals the lattice constant of aluminum metal at 300 k; thus, a = 0.405 nm. Considering the stability and capacity of calculations, the side length of the simulated box was set to “50a”. The dimensions of the specimens were L = h = 50a and thickness t = 25a. Two spherical voids were prefabricated in the constructed model with a void radius of r = 5a, as shown in Figure 2, and the distance between the two voids’ centers was three times the radius [13]. All MD simulations were first relaxed at temperature 300 k and 0 Pa pressure in three directions for 20 picosecond (ps) using the Nosé–Hoover thermostat and barostat [29]. The periodic boundary conditions were set in all simulated directions. The time step was 1 femtosecond (fs) in all simulations. The loading strain rate was controlled to achieve uniaxial compressive-tensile loading (Figure 2), which was utilized throughout the simulation, and the initial compressive strain was set to eleven groups: 0, −1, −2, −3, −4, −5, −6, −6.5, −7, −7.5, and −8%. The model was subjected to strain rate controlled uniaxial compressive loading and then unloaded, and after the equilibrium of 20 ps, the strain rate controlled tensile loading was performed. In order to better simulate the extreme loading conditions [8], the loading strain rate was set to 1 × 1010 s−1. The specimens were visualized by the open visualization tool (OVITO) [30] (Figure 3), and all atoms were colored by the centrosymmetry parameter (CSP). The green lines in the visual image are the dislocations marked with CSP.Figure 1
The stress-strain curves of boxes with different sizes such as 5a, 10a,… 70a, where a = 0.405 nm.Figure 2
Loading conditions and dimensions designed for simulations.Figure 3
(a) Three moments as the analysis objectives; (b) the snapshot of section view; (c–f) the atomic snapshot of different compressive strains at the “three important moments” shown in Figure3(a). (a) Three important moments: 1. the compressive section has just ended; 2. the material has reached yield; 3. the void has coalescence. (c) εic=0%. (d) εic=−5%. (e) εic=−6%. (f) εic=−7%. (g) εic=−8%.
(a)(b)(c)(d)(e)(f)(g)


## 3. Results and Discussion

### 3.1. Comparison of the Void Coalescence under Uniaxial Compressive-Tensile Loading with Different Compressive Strains

Figure4 illustrates the stress-strain curves under different compressive strains, in which the curve could be divided into three stages. In the first stage, except for the compressive strains −7.5 and −8%, the stress varied linearly with respect to the strain at a similar rate. Meanwhile, it can be seen that during the compressive phase, the curves of the compressive strains −7.5 and −8% were nonlinear, indicating that the material has entered the local plastic region. In the second stage, the stress-strain curve reached its peak, implying that the material yields. The enlarged section of Figure 4 shows that the yield strength fluctuates within a certain range before the initial compressive strain of −7%. From the compressive strain of −7%, the yield strength decreases significantly. However, when the compressive strain is −8%, the yield strength is higher than when the compressive strain is −7.5%. In addition, as the initial compressive strain decreases, the material will yield in advance. The curves enter the platform section in the third stage, which also represents that the material enters the plastic flow section.Figure 4
Stress-strain curves of different compressive strains (the inset is a partial enlarged view of the yield strength in the stress-strain curve).To analyze the phenomena appearing in the stress-strain curves more clearly, the model was visualized by OVITO, taking the center of the two voids as theY-axis and slicing in the X-direction to obtain the profile at different moments. This is equivalent to selecting a length on the Z-axis to study the plane strain issue. Hence, there is no need to include the influence of the boundary conditions in the Z-direction on the simulation. Three important moments have been selected as the objectives of the analysis in −7% of the compressive strain sample ( Figure 3(a)). Three moments correspond to the snapshots in Figures 3(b)–3(g), which are “the compressive phase has just ended,” “the material has reached its yield (ε≈0.1),” and “the void has got coalescence ε≈0.2.” Figure 3(b) is the section view of the initial configuration, where the axis and size are marked in Figure 3(b). Figure 3(c) shows that the voids are continuously developing. After growing to a certain size under uniaxial tensile loading, the ligament that is at the center of two voids begins necking and then the voids coalesce. In the second image, dislocations are generated between two voids. After coalescence, the dislocations originating from the voids continue to propagate. Note that Figures 3(d)–3(g) show that the voids still coalesce. In addition, when the initial compressive strain is −5%, there are no obvious dislocations after the compressive phase. In the second section, symmetric dislocations are also generated between two voids and in the process of void coalescence. Similarly, dislocations form an X-shaped distribution [31]. A different phenomenon can be observed where dislocations are generated between the voids at the end of the compression when the compressive strain is greater than −7%, and then dislocations evolve into an asymmetric distribution. The result shows that although the voids ultimately coalescence, dislocations have already occurred around the voids when the compressive phase has ended, as shown in Figure 3(e). Figure 3(g) shows that the voids also generate dislocations at the end of the compressive phase, but no coalescence occurs. In addition, dislocations are generated between the voids and are emitted in an approximate direction of 43° from the tensile direction under the compressive phase.In MD simulation, the void coalescence is obtained based on statistical methods, which is generally determined by merging two voids shown in the atomic snapshot. However, so far, there are no suitable parameters to depict the void coalescence, so we apply the coalescence strain to describe the void coalescence; that is, when two voids merge in the atomic snapshot, the current strain is recorded and defined as the coalescence strain. When the strain is equal to the coalescence strain, the voids coalesce. Figure5 shows the relationship between the initial compressive strain and coalescence strain. As the initial compressive strain increases, the coalescence strain also increases, and both increase linearly before the initial compressive strain of −7%. When the initial strain is −7.5%, because dislocations have already occurred during the compressive phase, the void coalescence begins to become more difficult. Figure 6(a) illustrates the relationship between the yield strength and compressive strain. It can be clearly observed from the figure that before the compressive strain of −7%, the yield strength fluctuates around 5.6 GPa (blue dash line). After the compressive strain of −7%, there is a significant decrease in the yield stress, and when the compressive strain is −8%, the yield strength is higher than when the compressive strain is −7.5%. The decline in the yield strength is due to the existence of BE. For further clarification, Figure 6(b) illustrates the relationship between the stress difference and compressive strain. Except for the sample with a compressive strain of -6%, all critical stress differences and the minimum stress difference fluctuate in the vicinity of the blue dash line, which is similar to what is described in the BE [32]. In addition, Figure 6(a) shows that the yield strength with a compressive strain of −8% is higher than that with a compressive strain of −7.5%, while for the compressive strain of −7.5%, the voids do not coalescence. This phenomenon why the voids cannot coalescence will be explained in next section.Figure 5
Relationship between the initial compressive strain and coalescence strain.Figure 6
(a) Initial compressive strain vs. yield stress; (b) the relationship between the stress difference and compressive strain.
(a)(b)


### 3.2. Dislocation Analysis of the Void Coalescence under Compressive-Tensile Loading

In order to profoundly understand the effect of the compressive phase on the void coalescence, Figures7 and 8 show the dislocation processing obtained by different methods at the end of the compressive phase. In Figure 7, CSP is employed to depict the dislocation, in which Figures 7(a)–7(c) illustrate that the specimen with an initial compressive strain from −5 to −6.5% has few dislocations during the compressive phase. Obviously, dislocations are generated in the specimen when the compression strains are −7 and −8% (Figures 6(d)–6(f)), so that an increase of the initial compressive strain would cause the dislocation to increase. A dislocation begins to form, mainly on the void surface, which connects into a piece in the middle of the voids and finally expands to the material boundary. Moreover, it is not isotropic when it expands, and the voids are arranged in more dislocations. Figure 8 shows the dislocation lines that are extracted by the dislocation extraction algorithm (DXA) [30, 33]. The dislocation density ρd can be calculated by the length of the dislocation line:(1)ρd=lV,where V is the total volume of the box at the initial moment, and l is the length of the dislocation line. As shown in Figure 8, in the compressive phase, the dislocations are produced continuously, which results in an increase in the dislocation density that affects the void coalescence. From the surface of two spherical voids, dislocations are absorbed in the form of dislocation lines which are formed around the voids. As the strain increases, both dislocation lines merge into a large dislocation group, which is accompanied by lots of disorderly dislocation bands. In the compressive phase, the evolution of voids is also accompanied by the adsorption-mediated growth of dislocations, similar to those found by Sills [34]. Figure 8(c) shows that after applying −8% compressive strain, the dislocations grow in a large amount, and the growth has an effect on the later tensile phase. Therefore, the dislocation density with −8% compressive strain is compared with −5% compressive strain. Figure 9 shows the dislocation density-stress-time curve of the voids, in which it can be seen that both samples in the first 5 ps are under compressive stress, so that the 0–5 ps segment is coincident. In the 5–8 ps stage, dislocations begin to grow quickly, and the dislocation density increases rapidly in the sample, for which the compressive strain is −8%. The dislocation density begins to decrease during 8–20 ps because the tensile stress begins to dominate the dislocation movement. It should be noted that the stress is nonlinear, which indicates that the voids have formed a local plastic zone at this time in the stress-time diagram. At the 20–30 ps stage, dislocations begin to form in the specimen with a compressive strain of −5%, and the time 20 ps corresponds to the moment when stress reaches the yield. Compared to the compressive strain of −5%, the dislocation of the sample with a compressive strain of -8% grows in the compressive phase, and it causes the yield point hysteresis on the stress-time curve. After 30 ps, both the stress-time and density-time curves enter the platform stage.Figure 7
Dislocation emission by the void surface by CSP. (a)εic=−5%. (b) εic=−6%. (c) εic=−6.5%. (d) εic=−7%. (e) εic=−7.5%. (f) εic=−8%.
(a)(b)(c)(d)(e)(f)Figure 8
The dislocation lines for the compressive strain from the void surface (-7%–8%). (a)εic=−7%, (b) εic=−7.5%, and (c) εic=−8%.
(a)(b)(c)Figure 9
The curve of dislocation-stress-time with different compressive strains.The beginning of plastic deformation is the slip of dislocations, and the large-scale evolution of dislocations makes other internal dislocations difficult to move. The direction of dislocation movement depends on the BE. Due to the presetting compressive strain, dislocations are generated in the void surface which causes back stress in the material. The pilling up of dislocations leads to the dislocations where closed ligaments in the tensile phase cannot continue to slip. In addition, the back stress prevents the positive direction deformation, which results in a decrease in the yield strength while the voids do not coalesce.


### 3.3. Temperature Effect under Compressive Loading on the Void Coalescence

The effects of temperature and strain rate on the void coalescence should also be considered in studying extreme conditions. Figure10 is a snapshot of the specimen under different temperatures subjected to compressive-tensile loading, in which the compressive strain is set to −5%. As can be seen, the effect of temperature on the void coalescence is obvious. For instance, when a strain is 20%, the voids in the 300 and 500 k get coalescence, while the void coalescence in the 100 k sample occurs at a strain of 50%, as shown in Figure 10(c). This behavior implies that a decrease in temperature will suppress the void coalescence. In addition, as the temperature increases, the yield strength gradually decreases. Therefore, the low temperature would inhibit the void coalescence. From the energy perspective, the transition and rearrangement of atoms can only obtain energy through external deformation, which makes it difficult for voids to coalesce. Similarly, the strain rate also has an influence on the void coalescence (Figure 11). After compressive and tensile loading, a decrease of the strain rate would cause more obvious noncoalescence of the voids. When the strain is 20%, only the specimen whose strain rate is 1010 s−1 has coalescence. This shows that at a lower strain rate, the compression effect has a stronger influence on the voids, and it is necessary to consider the influence of compression on the void coalescence.Figure 10
(a–c) The snapshots of the specimen with a strain of 20%; (a) 100 k atε=20%. (b) 300 k at ε=20%. (c) 500 k at ε=20%. (d) The stress-strain curve at different temperatures.
(a)(b)(c)(d)Figure 11
The stress-strain curve at different strain rates when the compressive strain is set to −5%.


## 3.1. Comparison of the Void Coalescence under Uniaxial Compressive-Tensile Loading with Different Compressive Strains

Figure4 illustrates the stress-strain curves under different compressive strains, in which the curve could be divided into three stages. In the first stage, except for the compressive strains −7.5 and −8%, the stress varied linearly with respect to the strain at a similar rate. Meanwhile, it can be seen that during the compressive phase, the curves of the compressive strains −7.5 and −8% were nonlinear, indicating that the material has entered the local plastic region. In the second stage, the stress-strain curve reached its peak, implying that the material yields. The enlarged section of Figure 4 shows that the yield strength fluctuates within a certain range before the initial compressive strain of −7%. From the compressive strain of −7%, the yield strength decreases significantly. However, when the compressive strain is −8%, the yield strength is higher than when the compressive strain is −7.5%. In addition, as the initial compressive strain decreases, the material will yield in advance. The curves enter the platform section in the third stage, which also represents that the material enters the plastic flow section.Figure 4
Stress-strain curves of different compressive strains (the inset is a partial enlarged view of the yield strength in the stress-strain curve).To analyze the phenomena appearing in the stress-strain curves more clearly, the model was visualized by OVITO, taking the center of the two voids as theY-axis and slicing in the X-direction to obtain the profile at different moments. This is equivalent to selecting a length on the Z-axis to study the plane strain issue. Hence, there is no need to include the influence of the boundary conditions in the Z-direction on the simulation. Three important moments have been selected as the objectives of the analysis in −7% of the compressive strain sample ( Figure 3(a)). Three moments correspond to the snapshots in Figures 3(b)–3(g), which are “the compressive phase has just ended,” “the material has reached its yield (ε≈0.1),” and “the void has got coalescence ε≈0.2.” Figure 3(b) is the section view of the initial configuration, where the axis and size are marked in Figure 3(b). Figure 3(c) shows that the voids are continuously developing. After growing to a certain size under uniaxial tensile loading, the ligament that is at the center of two voids begins necking and then the voids coalesce. In the second image, dislocations are generated between two voids. After coalescence, the dislocations originating from the voids continue to propagate. Note that Figures 3(d)–3(g) show that the voids still coalesce. In addition, when the initial compressive strain is −5%, there are no obvious dislocations after the compressive phase. In the second section, symmetric dislocations are also generated between two voids and in the process of void coalescence. Similarly, dislocations form an X-shaped distribution [31]. A different phenomenon can be observed where dislocations are generated between the voids at the end of the compression when the compressive strain is greater than −7%, and then dislocations evolve into an asymmetric distribution. The result shows that although the voids ultimately coalescence, dislocations have already occurred around the voids when the compressive phase has ended, as shown in Figure 3(e). Figure 3(g) shows that the voids also generate dislocations at the end of the compressive phase, but no coalescence occurs. In addition, dislocations are generated between the voids and are emitted in an approximate direction of 43° from the tensile direction under the compressive phase.In MD simulation, the void coalescence is obtained based on statistical methods, which is generally determined by merging two voids shown in the atomic snapshot. However, so far, there are no suitable parameters to depict the void coalescence, so we apply the coalescence strain to describe the void coalescence; that is, when two voids merge in the atomic snapshot, the current strain is recorded and defined as the coalescence strain. When the strain is equal to the coalescence strain, the voids coalesce. Figure5 shows the relationship between the initial compressive strain and coalescence strain. As the initial compressive strain increases, the coalescence strain also increases, and both increase linearly before the initial compressive strain of −7%. When the initial strain is −7.5%, because dislocations have already occurred during the compressive phase, the void coalescence begins to become more difficult. Figure 6(a) illustrates the relationship between the yield strength and compressive strain. It can be clearly observed from the figure that before the compressive strain of −7%, the yield strength fluctuates around 5.6 GPa (blue dash line). After the compressive strain of −7%, there is a significant decrease in the yield stress, and when the compressive strain is −8%, the yield strength is higher than when the compressive strain is −7.5%. The decline in the yield strength is due to the existence of BE. For further clarification, Figure 6(b) illustrates the relationship between the stress difference and compressive strain. Except for the sample with a compressive strain of -6%, all critical stress differences and the minimum stress difference fluctuate in the vicinity of the blue dash line, which is similar to what is described in the BE [32]. In addition, Figure 6(a) shows that the yield strength with a compressive strain of −8% is higher than that with a compressive strain of −7.5%, while for the compressive strain of −7.5%, the voids do not coalescence. This phenomenon why the voids cannot coalescence will be explained in next section.Figure 5
Relationship between the initial compressive strain and coalescence strain.Figure 6
(a) Initial compressive strain vs. yield stress; (b) the relationship between the stress difference and compressive strain.
(a)(b)


## 3.2. Dislocation Analysis of the Void Coalescence under Compressive-Tensile Loading

In order to profoundly understand the effect of the compressive phase on the void coalescence, Figures7 and 8 show the dislocation processing obtained by different methods at the end of the compressive phase. In Figure 7, CSP is employed to depict the dislocation, in which Figures 7(a)–7(c) illustrate that the specimen with an initial compressive strain from −5 to −6.5% has few dislocations during the compressive phase. Obviously, dislocations are generated in the specimen when the compression strains are −7 and −8% (Figures 6(d)–6(f)), so that an increase of the initial compressive strain would cause the dislocation to increase. A dislocation begins to form, mainly on the void surface, which connects into a piece in the middle of the voids and finally expands to the material boundary. Moreover, it is not isotropic when it expands, and the voids are arranged in more dislocations. Figure 8 shows the dislocation lines that are extracted by the dislocation extraction algorithm (DXA) [30, 33]. The dislocation density ρd can be calculated by the length of the dislocation line:(1)ρd=lV,where V is the total volume of the box at the initial moment, and l is the length of the dislocation line. As shown in Figure 8, in the compressive phase, the dislocations are produced continuously, which results in an increase in the dislocation density that affects the void coalescence. From the surface of two spherical voids, dislocations are absorbed in the form of dislocation lines which are formed around the voids. As the strain increases, both dislocation lines merge into a large dislocation group, which is accompanied by lots of disorderly dislocation bands. In the compressive phase, the evolution of voids is also accompanied by the adsorption-mediated growth of dislocations, similar to those found by Sills [34]. Figure 8(c) shows that after applying −8% compressive strain, the dislocations grow in a large amount, and the growth has an effect on the later tensile phase. Therefore, the dislocation density with −8% compressive strain is compared with −5% compressive strain. Figure 9 shows the dislocation density-stress-time curve of the voids, in which it can be seen that both samples in the first 5 ps are under compressive stress, so that the 0–5 ps segment is coincident. In the 5–8 ps stage, dislocations begin to grow quickly, and the dislocation density increases rapidly in the sample, for which the compressive strain is −8%. The dislocation density begins to decrease during 8–20 ps because the tensile stress begins to dominate the dislocation movement. It should be noted that the stress is nonlinear, which indicates that the voids have formed a local plastic zone at this time in the stress-time diagram. At the 20–30 ps stage, dislocations begin to form in the specimen with a compressive strain of −5%, and the time 20 ps corresponds to the moment when stress reaches the yield. Compared to the compressive strain of −5%, the dislocation of the sample with a compressive strain of -8% grows in the compressive phase, and it causes the yield point hysteresis on the stress-time curve. After 30 ps, both the stress-time and density-time curves enter the platform stage.Figure 7
Dislocation emission by the void surface by CSP. (a)εic=−5%. (b) εic=−6%. (c) εic=−6.5%. (d) εic=−7%. (e) εic=−7.5%. (f) εic=−8%.
(a)(b)(c)(d)(e)(f)Figure 8
The dislocation lines for the compressive strain from the void surface (-7%–8%). (a)εic=−7%, (b) εic=−7.5%, and (c) εic=−8%.
(a)(b)(c)Figure 9
The curve of dislocation-stress-time with different compressive strains.The beginning of plastic deformation is the slip of dislocations, and the large-scale evolution of dislocations makes other internal dislocations difficult to move. The direction of dislocation movement depends on the BE. Due to the presetting compressive strain, dislocations are generated in the void surface which causes back stress in the material. The pilling up of dislocations leads to the dislocations where closed ligaments in the tensile phase cannot continue to slip. In addition, the back stress prevents the positive direction deformation, which results in a decrease in the yield strength while the voids do not coalesce.


## 3.3. Temperature Effect under Compressive Loading on the Void Coalescence

The effects of temperature and strain rate on the void coalescence should also be considered in studying extreme conditions. Figure10 is a snapshot of the specimen under different temperatures subjected to compressive-tensile loading, in which the compressive strain is set to −5%. As can be seen, the effect of temperature on the void coalescence is obvious. For instance, when a strain is 20%, the voids in the 300 and 500 k get coalescence, while the void coalescence in the 100 k sample occurs at a strain of 50%, as shown in Figure 10(c). This behavior implies that a decrease in temperature will suppress the void coalescence. In addition, as the temperature increases, the yield strength gradually decreases. Therefore, the low temperature would inhibit the void coalescence. From the energy perspective, the transition and rearrangement of atoms can only obtain energy through external deformation, which makes it difficult for voids to coalesce. Similarly, the strain rate also has an influence on the void coalescence (Figure 11). After compressive and tensile loading, a decrease of the strain rate would cause more obvious noncoalescence of the voids. When the strain is 20%, only the specimen whose strain rate is 1010 s−1 has coalescence. This shows that at a lower strain rate, the compression effect has a stronger influence on the voids, and it is necessary to consider the influence of compression on the void coalescence.Figure 10
(a–c) The snapshots of the specimen with a strain of 20%; (a) 100 k atε=20%. (b) 300 k at ε=20%. (c) 500 k at ε=20%. (d) The stress-strain curve at different temperatures.
(a)(b)(c)(d)Figure 11
The stress-strain curve at different strain rates when the compressive strain is set to −5%.


## 4. Conclusions

To understand why the void coalescence stopped after applying some compressive loading, the MD method was used in this paper to simulate single crystal aluminum under uniaxial compressive-tensile plastic deformation with nanovoids. Moreover, the effects of temperature and strain rate on the void coalescence were investigated, and the yield strength of specimens under different compressive strains was analyzed. It was found that the compression not only reduces the yield strength of the material, but also obstructs the void coalescence, which also modifies the existing damage constitutive model. The main conclusions are as follows:(i)
Before the compressive strain of -7%, the yield strength fluctuates within a certain range, and the yield strength decreases with the increasing compressive strain after -7%. Compared with the uniaxial tensile loading, this behavior can be explained by BE. Note that if the compression effect is not considered in engineering applications and material manufacturing, it will cause damage to the material before its service life, posing a great risk.(ii)
The compressive phase causes dislocation absorption at the void ligament, and the dislocation growth obstructs the dislocation movement during the tensile phase. Therefore, the compressive phase obstructs the void coalescence. This is why the yield strength ofεic=−8% is higher than εic=−7.5%. The temperature and strain rate also have an obvious effect on the void coalescence under compressive and tensile loading. A low temperature would suppress the void coalescence, and the strain rate would inhibit the void coalescence as well.



---

*Source: 9990161-2022-02-28.xml*