Abstract
In order to study the dynamic mechanical properties of recycled concrete with ferronickel slag, dynamic experiments on cylindrical specimens (φ 73 mm × 36.5 mm) were performed using a split Hopkinson pressure bar (SHPB) system. The effects of loading rate and recycled aggregate replacement rate on the compressive strength, toughness index, peak strain, and failure pattern of the concrete were investigated, and four different loading rates and aggregate replacement rates were selected. In addition, cohesive elements were embedded between C3D8R elements in the numerical model using Python language to simulate the failure process of the recycled concrete, and the reliability of the simulation was validated by the obtained experimental results. It was found that the peak strain of the recycled concrete was significantly higher than that of the ordinary concrete; however, its peak stress was lower than that of the ordinary concrete. In comparison to the recycled concretes containing 70% and 100% recycled coarse aggregates (RCAs), the toughness index amplification of the recycled concrete containing 30% RCAs was the most obvious. The failure pattern and mechanical properties of the numerical model were consistent with those of the test specimens, thus validating the reliability of the numerical simulation.
1. Introduction
Ferronickel slag (FNS) is an industrial metallurgy slag produced in ferronickel manufacturing processes. It is estimated that about 14 tons of slag are generally produced per ton of nickel [1]. China produces about 4 million tons of slag every year [2]; however, less than 10% of slag is actually utilized. FNS accumulates in the open air, and its harmful substances cause soil quality to decline. Additionally, the dust dispersion of waste during transportation reduces air quality. Therefore, it is pointed out that the recovery and utilization of FNS have good economic and social benefits [3–5]. Moreover, recycled concretes with ferronickel slag (RCFNS) have good strength [6].
The chemical composition of FNS mainly includes SiO2, MgO, CaO, Al2O3, and Fe2O3 [7], and crystalline minerals in nickel-iron slag consist of enstatite, forsterite, and diopside [8]. The content of amorphous minerals in FNS is generally higher than 50% [9]. FNS powder has a pozzolanic activity and is a potential mineral admixture for concrete. Dourdounis et al. [10] studied the utilization of ferronickel electroreduction furnace (FeNi-ERF) slag for the production of high-alumina cement and found that the reduction smelting process of ERF slag with limestone and low-grade diasporic bauxite had considerable potential for optimization. Katsiotis et al. [11] studied the performance and hydration of a mixed cement containing FNS powder and noticed that the addition of FNS reduced the water requirement of the mixed cement and prolonged the setting time. Ding et al. [12] found that the FNS powder-containing concrete exhibited good performance at a suitable replacement rate.
Although the application of recycled aggregate concretes (RAC) can reduce construction waste and save natural resources, their strength is lower than that of natural aggregate concretes [13]. Etxeberria et al. [14] studied the production of RAC with four different recycled aggregate replacement rates and found that the RAC with 25% of recycled coarse aggregates (RCA) achieved the same mechanical properties as those of conventional concretes. Limbachiya et al. [15] studied the effects of RCA content on the ceiling strength, bulk engineering, and durability of RAC and reported that 30% coarse RCA had no effect on the concrete strength; however, the concrete strength decreased with the increase of RCA content.
Nowadays, due to the increase of extreme incidents, such as chemical explosions, terrorist attacks, high-speed car crashes, and explosions, widespread concern about the security of building structures is raised. In an explosion, tremendous energy is released in a very short period and causes great damage. Therefore, in order to apply RAC members in structural engineering, it is necessary to study their dynamic behavior under impact loading. Rao et al. [16] investigated the behavior of RAC beams with varied RAC ratios subjected to low-velocity impacts and found that 25% RCA did not influence the concrete strength. Li et al. [17] investigated the effects of nanoparticles on the dynamics behavior of RAC members under impact loading by a split Hopkinson pressure bar (SHPB) system and found that nano-SiO2 particles were more effective to improve the dynamic compressive strength of RAC. Mastali et al. [18] studied the impact resistance and mechanical properties of self-compacting concrete reinforced with recycled glass fiber-reinforced polymers (GFRP) and noticed that the impact resistance and mechanical properties of the concrete were linearly correlated. Li et al. [19] investigated the behavior of a recycled aggregate concrete-filled steel tube (RACFST) at elevated temperatures by experimental and numerical studies and reported that the impact behavior of the tube was improved with the increasing steel ratio.
In view of the above analysis, the strength of recycled concrete is generally lower than that of natural aggregate concrete, and the ferronickel slag is helpful to improve the strength of concrete. Therefore, in order to promote the application of recycled concrete in impact engineering, ferronickel slag was used to make up for the lack of bearing capacity of recycled concrete, and the impact tests of recycled concrete containing ferronickel slag (RCFS) were carried out to determine its dynamic mechanical properties. In the present work, dynamic experiments of RCFS specimens were performed using a SHPB system at different replacement rates of RCA and loading rates. The compressive strength, toughness index, peak strain, and failure pattern of the concrete specimens were analyzed. In addition, the failure process of recycled concrete with ferronickel slag was simulated in ABAQUS software. The failure mechanism of the numerical model was analyzed, and the reliability of the numerical simulation was validated by experimental results.
2. Experiment
2.1. Materials
Both the ordinary concrete and the recycled concrete were composed of P·O 42.5 cement, and their strengths were 60 MPa. Blast furnace FNS prepared by water quenching was used in this research. The density, specific surface area, average particle diameter, and particle size distribution of FNS were 2.86 g/cm3, 425.6 m2/kg, 13.953 μm, and 1–100 μm, respectively. The primary components of FNS and cement are presented in Table 1. The X-ray diffraction (XRD) curve of FNS is displayed in Figure 1. The main components of FNS were CaCO3, MgAlO4, and (Mg, Fe)2SiO4. The main physical properties of cement are presented in Table 2.

The water-cement ratio of the ordinary concrete specimen (N0-60) was 0.32, and no RCA was used in this concrete. In order to ensure the fluidity and compactness of the concrete mixing process, a water reducing agent and silica fume are added into the ordinary concrete and the recycled concrete. Moreover, fine aggregates of natural river sand with a fineness modulus of 2.7 were added to both types of concretes. RCAs of the recycled concrete were crushed particles (size = 16 mm) of C60 ordinary concrete. The mixing proportions of these two types of concretes are depicted in Table 3.
2.2. Specimens
Impact tests of the two types of concrete were carried out by a SHPB system with a diameter of 80 mm. The main parameters of the tests were the loading rate and the replacement rate of RCA. The diameter and height of the test specimens were 73 mm and 36.5 mm, respectively (Figure 2). The propagation of one-dimensional stress waves in the specimens was assumed, and the inertial and dispersion effects under an impact were eliminated [20]. After pouring the specimens, the standard curing was carried out for 28 days, and then the cross sections of the specimens were polished to ensure the flatness requirements of the International Society for Rock Mechanics (ISRM).

In order to facilitate the analysis of test data, the specimens were numbered. For example, in the R30-60-1400 specimen, R represents recycled concrete, 30 represents 30% RCAs, 60 represents the concrete strength of 60 MPa, and 1400 represents the loading rate of 1400 GPa/s.
2.3. Experimental Equipment
SHPB systems have a history of more than 100 years and are widely used to determine the dynamic mechanical properties of materials [21–25].
The pressure bar of the SHPB system used in this experiment was made of high-strength alloy steel (elastic modulus = 210 GPa and Poisson’s ratio = 0.25). The diameter of the bar was 80 mm, and the length of the incident and transmitted bars was 3000 mm and 2040 mm, respectively. A damper was installed behind the transmission bar. The air pressure of the launching device of the SHPB system was electronically controlled to ensure a stable launching speed of the impact bar. In addition, an infrared velocimeter was used to capture the speed of the impact bar. The test device is displayed in Figure 3.

In order to investigate the influence of loading rate on the mechanical properties of the specimens, four different loading rates of 1400 GPa/s, 1500 GPa/s, 1600 GPa/s, and 1700 GPa/s were selected.
A DH3840 ultradynamic strain indicator was used to capture voltage signals of the pressure bar based on the half-bridge method. Strain gauges on the incident and transmission bars were connected to the superdynamic strain indicator. Voltage fluctuations of the strain gauges were caused by stress waves, and the detected voltage signals were transformed into force versus time by equations (1) and (2) (Figure 4):where K denotes the sensitivity coefficient of the strain gauges, U0 denotes the output voltage signal, E0 denotes the bridge voltage, n denotes the gain coefficient (n = 300), E is Young’s modulus of the incident and transmission bars, A is the cross-sectional area of the bars, and , , and are the recorded strains.

(a)

(b)

(c)

(d)
It is noticeable from Figure 4 that superposition waves were formed by incident waves, and reflected waves were in good agreement with transmitted waves. Therefore, the dynamic force balance of the specimens could be realized when stress waves propagated across the specimens several times, and the quasistatic theory was used to analyze the specimens.
3. Test Results and Analysis
In this research, four types of concrete specimen (N0-60, R30-60, R70-60, and R100-60) were investigated. For each type, four loading rates were selected, and three specimens were tested under each loading rate. A total of 48 specimens were prepared in this test. Due to a large number of test specimens, some typical test results were selected for the analysis.
3.1. Failure Phenomenon and SEM Micrograph
Taking the R30-60 specimen as an example, the detected failure patterns were divided into four grades over a wide range of loading rates (Figure 5). Grade 1 indicates that the specimen was split into two parts under the loading rate of 1400 GPa/s; however, these two parts experienced little damage and no spalling. Grade 2 indicates that the specimen was split into numerous fragments under the loading rate of 1500 GPa/s. Grade 3 indicates that the specimen still had numerous fragments under the loading rate of 1600 GPa/s; however, their volume was smaller than that of Grade 2. Grade 4 indicates that the specimen was completely crushed under the loading rate of 1700 GPa/s.

(a)

(b)

(c)

(d)
Figure 6 displays the SEM micrographs of four types of specimens (N0-60-1400, R30-60-1400, and R70-60-1400) under a local magnification of 6000 times. The N0-60-1400 specimen contained a large amount of hydration products (Figure 6(a)). A large amount of Al2O3 and SO2 was present in FNS, and these substances reacted with Ca(OH)2 produced by cement. The C-S-H gel could fill porosity in the specimen and combine aggregates and cementitious materials more closely. In the R30-60-1400 specimen, hydration products generated on the surface contained a large amount of acicular ettringite (Figure 6(b)). The formation of ettringite led to an expansion of this concrete and affected the interface strength between aggregates and the mortar. In the R70-60-1400 specimen, numerous granular substances and holes were present on the surface, the boundary between the mortar and aggregates was prominent, and the adhesion between particles was weak (Figure 6(c)).

(a)

(b)

(c)
3.2. Relationship between Loading Rate and Stress-Strain
Figure 7 displays the mechanical properties of the concrete specimens (N0-60, R30-60, R70-60, and R100-60) under the loading rates of 1400 GPa/s and 1700 GPa/s. The stress-strain curves were linear in the early stage and nonlinear in the later stage [26].

(a)

(b)

(c)
It is clear from Figure 7(a) that when the loading rate was 1400 GPa/s, the peak stress of the R30-60-1400 specimen was 0.1% lower than that of N0-60-1400. When the loading rate reached 1700 GPa/s, the peak stress of the N0-60-1700 specimen was 9% higher than that of R30-60-1700. In comparison to the N0-60 specimen, the peak stress of the R30-60 concrete increased more slowly with increasing loading rate.
It is noticeable from Figure 7(b) that when the loading rate was 1400 GPa/s, the peak stress of the R70-60-1400 specimen was higher than that of N0-60-1400. When the loading rate reached 1700 GPa/s, the peak strain of the R70-60-1700 specimen was higher than that of N0-60-1700; however, its peak stress was lower than that of N0-60-1700. In comparison to the N0-60 specimen, the peak stress of the R70-60 concrete increased more slowly with the increasing loading rate.
When the RCA content in the recycled concretes varied from 30% to 100%, the peak strain of the recycled concretes was significantly higher than that of the N0-60 specimen; however, their peak stress was lower than that of N0-60. When the replacement rate of RCA reached 100%, the peak stress decreases by 21% under the loading rate of 1700 GPa/s.
3.3. Toughness Index
When a specimen is destroyed under impact, the ability of the specimen to absorb the impact energy is called toughness, which is used to evaluate the strength and ductility of materials. The area enclosed by a stress-strain curve and coordinate axes is defined as toughness [27–30]. The toughness values of four types of concretes are presented in Figure 8, and the toughness index was expressed by T.

The toughness index of the N0-60 specimen increased from 85.76 J/cm3 to 138.24 J/cm3 (increment of 52.48 J/cm3) with the increasing loading rate from 1400 GPa/s to 1700 GPa/s. The toughness index of the R30-60 specimen increased from 71.37 J/cm3 to 161.32 J/cm3 (increment of 89.95 J/cm3) as the loading rate increased from 1400 GPa/s to 1700 GPa/s. The toughness index of the R70-60 specimen increased from 55.21 J/cm3 to 126.25 J/cm3 (increment of 71.04 J/cm3) when the loading rate increased from 1400 GPa/s to 1700 GPa/s. The toughness index of the R100-60 specimen increased by only 55.32 J/cm3. Therefore, the difference of the toughness index of the recycled concrete with the increasing loading rate from 1400 GPa/s and 1700 GPa/s decreased with the increase of RCA content.
Under the same replacement rate of RCA, the toughness index of the recycled concrete increased with the increase of the loading rate. In comparison to the recycled concretes containing 70% and 100% RCAs (R70-60 and R100-60 specimens, respectively), the toughness index increase of the R30-60 specimen was the most obvious.
3.4. Effect of the Replacement Rate of RCA on Peak Strain
Peak strain is an important index to indicate the deformation and instability of a specimen after an impact. Figure 9 displays the relationship between the peak strain of each specimen and the loading rate.

It is clear that the peak strain of each concrete increased with the increase of the loading rate. When the loading rate increased from 1400 GPa/s to 1700 GPa/s, the peak strains of N0-60, R30-60, and R100-60 increased by 1.9 × 10−3, 8.67 × 10−3, and 2.166 × 10−2, respectively.
Table 4 presents a summary of the obtained test results. The average peak strain, average toughness index, and average peak stress of the specimens had an increasing trend with the increase of the loading rate. When the loading rate was 1700 GPa/s, the N0-60 specimen had the largest average peak stress, whereas the average peak strain of R100-60 was found to be the largest. When the loading rate was 1400 GPa/s, all four types of specimens experienced only cracking phenomenon. When the loading rate varied between 1500 GPa/s and 1600 GPa/s, large particle fragments were generated in all four types of specimens. When the loading rate was 1700 GPa/s, all four types of specimens were crushed.
4. Numerical Simulation
Taking the R100-60-1700 specimen as an example, the secondary development of ABAQUS was carried out by Python language to realize the batch embedding of cohesive elements and the random generation of coarse aggregates in the numerical model, and by this method, the random crack propagation and random element shedding of the model were realized. Moreover, the reliability of the numerical simulation was validated by the obtained experimental results.
4.1. Cohesive Model
The cohesive model is first introduced by Barenblatt [31] to simulate the decohesion of atomic lattices by the traction-separation law. The cohesive model of ABAQUS also adopts this damage mechanics model. The linear triangular cohesive model shown in Figure 10 is generally used to simulate the fracture of brittle materials. , , and GC in Figure 10 represent the initial damage displacement, the effective displacement of traction force, and the energy release rate, respectively.

4.2. Damage Initiation Criteria
When the contact traction or separation reaches the peak point A in Figure 10, the damage of the cohesive element initiates. In the present research, the damage initiation criteria of ABAQUS cohesive elements adopted the maximum nominal stress criterion.where is normal traction force on the fracture surface, and are two tangential traction forces on the fracture surface, and , , and are the peak values of anisotropic nominal stress. When f = 1.0, the damage starts.
4.3. Damage Evolution
The damage evolution law describes that the cohesive model reaches the corresponding initiation criterion during the impact, and the cohesive stiffness begins to degenerate. The degradation process is indicated by the straight line AB in Figure 9. The damage process of cohesive elements can be expressed aswhere , , and are the contact stress components predicted by the elastic traction-separation behavior for the current separation without damage D (damage factor). If damage evolution is modeled, D monotonically evolves from 0 to 1 upon further loading after the initiation of damage.
4.4. Numerical Model
ABAQUS, a commercial finite element software, provides COH3D8 elements to simulate crack propagation. Crack propagation along any path can be realized by the batch embedding of COH3D8 elements and is affected by the cohesive strength. The functions of ABAQUS software cannot implement the batch embedding of COH3D8 elements into a model. However, the batch embedding of COH3D8 elements can be realized by processing the INP file in ABAQUS using Python language. Therefore, the node number of C3D8R elements in a model can be obtained, and COH3D8 elements with an extremely thin thickness are formed by the node number of C3D8R elements (shown in region A of Figure 11). By this method, cracks can propagate randomly between C3D8R elements during an impact, and C3D8R elements in the model can be separated.

Figure 11 shows the 8-node hexahedral linear reduced integral element (C3D8R). Compared with the complete integration element, the linear reduced integration element only contains one integration point in the center of the element, while the number of integration points of the quadratic reduced integration element is the same as that of the linear complete integration element. The linear reduced integral element has the following advantages: (1) it is not easy to shear self-locking under bending load; (2) the result of displacement is more accurate; and (3) when the mesh is distorted, the accuracy of the analysis will not be greatly affected.
Therefore, C3D8R elements are widely used to establish the numerical model of concrete. When the cohesive elements are embedded, the C3D8R elements can fall off freely when the concrete model is damaged, which is very similar to the concrete particle falling off in the test.
Two materials were included in the concrete numerical model: mortar and coarse aggregate. The random generation of coarse aggregate was realized using Python language to process C3D8R element numbers in the INP file of the ABAQUS code, and the length, width, and height of coarse aggregate model are between 12 mm and 19 mm (Figure 11). C3D8R elements were used for both mortar and coarse aggregate, and COH3D8 elements were embedded between C3D8R elements. The total numbers of COH3D8 and C3D8R elements in the model were 164,400 and 56,000, respectively. For brittle heterogeneous materials, especially concrete materials, only linear elastic behavior is considered in C3D8R elements (Table 5). The elastic modulus, Poisson’s ratio, and density of COH3D8 elements were similar to those of surrounding C3D8R elements; however, COH3D8 elements also had tensile and shear strengths, which could be determined by splitting and shear tests (Table 5). In addition, related concrete studies have reported that the spacing between concrete particles is approximately 0.01 mm [32]; hence, the thickness of COH3D8 elements in the model was set as 0.01 mm.
According to the schematic view of the SHPB system in Figure 3, its simulation model was established. The linear elastic model was adopted for the incident and transmission bars. A hard contact without friction between the pressure bar and the specimen was considered. Stress waves in the R100-60 specimen are displayed in Figure 4.
4.5. Numerical Simulation Results
The numerical simulation results of the R100-60-1700 specimen are displayed in Figure 12, and the failure pattern of the numerical model was consistent with the test result.

During the dynamic loading process at t = 673–688 μs, stress concentration occurred at the incident loading end of the model. When the loading time was t = 792 μs, an extrusion protrusion appeared in the middle of the model. During the dynamic loading process at t = 847–870 μs, the element damage of the model was more serious, and the elements appeared to fall off. When the loading time was t = 910 μs, the volume of the model was greatly compressed, and a large number of elements at the loading end and middle part of the model fall off. When the loading time was t = 1024 μs, the mortar region in the model was completely crushed. At this moment, the failure phenomenon of the coarse aggregate region is exhibited in Figure 13. The failure phenomenon of residual coarse aggregate in the numerical model was consistent with that of the test specimen.

In order to gain further insight into the failure mechanism, the stress-strain curve obtained by the numerical simulation was compared with that obtained from the experiment (Figure 14(a)). It is noticeable that the stress-strain curve obtained by the numerical simulation was consistent with that obtained from the experiment, and the peak stress difference between the two curves was only 1.7%, confirming the reliability of the numerical simulation.

(a)

(b)
Figure 14(b) displays the peak stress versus loading rate curves obtained from the simulation and the experiment. The maximum and minimum peak stress differences between the two curves were 3.5% and 0.6%, respectively.
5. Conclusions
(1)The failure of the recycled concrete was divided into four grades under different loading rates. When the loading rates were 1400 GPa/s, 1500 GPa/s, 1600 GPa/s, and 1700 GPa/s, the failure patterns of the concrete were graded as 1 to 4, respectively.(2)The peak strain of the recycled concrete was significantly higher than that of the ordinary concrete; however, its peak stress was lower than that of the ordinary concrete. When the replacement rate of RCA reached 100%, the peak stress decreased by 21%.(3)Under the same replacement rate of RCA, the toughness index of the recycled concrete increased with the increasing loading rate. In comparison to the recycled concretes containing 70% and 100% RCAs, the toughness index amplification of the recycled concrete containing 30% RCAs was the most obvious.(4)In the numerical model, cohesive elements of the mortar region were completely discrete, whereas those of the coarse aggregate region were only partially discrete. The failure pattern and mechanical properties of the numerical model were consistent with those of the test specimens, thus validating the reliability of the numerical simulation.Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This work was financially supported by the National Natural Science Foundation of China (Grant nos. 11672194 and U19A2098). The authors thank LetPub (https://www.letpub.com) for its linguistic assistance during the preparation of this manuscript.