Next Article in Journal
The Rise of the Xenes: From the Synthesis to the Integration Processes for Electronics and Photonics
Previous Article in Journal
Influence of Pozzolan, Slag and Recycled Aggregates on the Mechanical and Durability Properties of Low Cement Concrete
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Steel-Reinforced Reactive Powder Concrete Beam Based on Bond-Slip

1
College of Civil Engineering and Architecture, Guangxi University, Nanning 530004, China
2
Key Laboratory of Disaster Prevention and Structural Safety of Chinese Ministry of Education, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(15), 4176; https://doi.org/10.3390/ma14154176
Submission received: 21 June 2021 / Revised: 23 July 2021 / Accepted: 23 July 2021 / Published: 27 July 2021

Abstract

:
In this study, based on the concrete damaged plasticity (CDP) model in the ABAQUS software, various plastic damage factor calculation methods were introduced to obtain CDP parameters suitable for reactive powder concrete (RPC) materials. Combined with the existing tests for the bending performance of steel-reinforced RPC beams, the CDP parameters of the RPC material were input into ABAQUS to establish a finite element model considering the bond and slip between the steel and RPC for numerical simulation. The load-deflection curve obtained by the simulation was compared with the measured curve in the experiment. The results indicated that on the basis of the experimentally measured RPC material eigenvalue parameters, combined with the appropriate RPC constitutive relationship and the calculation method of the plastic damage factor, the numerical simulation results considering the bond-slip were in good agreement with the experimental results with a deviation of less than 10%. Thus, it is recommended to select a gentle compressive stress-strain curve in the descending section, an approximate strengthening model of the tensile stress-strain curve, and to use the energy loss method and Sidoroff’s energy equivalence principle to calculate the RPC plastic damage parameters.

1. Introduction

The concrete damaged plasticity (CDP) model in the ABAQUS finite element software is widely used for the simulation analysis of the mechanical properties of ordinary concrete structural members. At present, several calculation methods for CDP model parameters of concrete have been proposed [1,2,3,4,5], including the calculation of concrete damage based on energy loss, specific strain method or graphic method, and a combination of elastic modulus damage and energy equivalence principle of Sidoroff. However, these methods utilize the CDP model parameters of normal strength concrete in ABAQUS. Although the method to calculate the CDP model parameters for ordinary concrete has been established, its applicability for reactive powder concrete (RPC) has not been studied. Thus, the calculation of the CDP model parameters for RPC is of utmost importance for the numerical simulation of RPC structures using ABAQUS software. However, only a few studies have been conducted on CDP model parameters for RPC [6,7,8]. In addition, a national standard for the RPC tension compression constitutive equation in China is lacking, which is an important factor for the calculation of CDP model parameters for RPC.
In this paper, we introduce several methods for the calculation of CDP model parameters in ABAQUS software. These methods were applied to calculate the CDP model parameters of RPC. The calculated parameters were compared and analyzed to determine their applicability. This paper reviews the constitutive relationship of RPC. Furthermore, the relationships were compared and classified to determine the appropriate constitutive relationship according to the strength characteristic value of RPC measured using existing tests. The calculated CDP model parameters were applied to the ABAQUS finite element software to simulate the bending resistance of steel-reinforced RPC beams. Comparing the results with the measured load-deflection curve, the rationality of the selected calculation method and RPC constitutive relationship was verified, thereby providing a reference for numerical analysis.

2. Parameter Calculation Method of Plastic Damage Model

2.1. Specific Strain Method or Graphic Method

The definition of uniaxial compression and tension behavior of concrete in ABAQUS is depicted in Figure 1.
As shown in Figure 1, the expression of the stress-strain relationship of concrete under uniaxial compression and tension can be defined as follows:
σ k = ( 1 d k ) E 0 ( ε k ε k p l ) , ( k = t , c )
The damage factor in the above formula is calculated as follows:
d k = 1 σ k ( ε k ε k p l ) E 0 , ( k = t , c )
where σ k represents tension stress, ε k denotes tension strain, ε k p l represents equivalent plastic strain, ε k i n represents inelastic strain, E 0 denotes the initial elastic modulus of the material.
Substituting ε k = ε k i n + ε 0 k e l , ε 0 k e l = σ k / E 0 into (2), we obtain
d k = 1 σ k ( σ k + ε k p l ( 1 / b k 1 ) E 0 )
where k = t , c represents tension and compression, respectively; b k denotes the ratio of inelastic strain ε k i n to equivalent plastic strain ε k p l . According to the cyclic loading and unloading stress path of the concrete specimen, Birtel et al. [1] demarcated b k and suggested that b c = 0.7 and b t = 0.7 . Zhang Jin et al. [2] suggested that the value of b c was within the range of 0.35–0.7 under compression and the value of b t was within the range of 0.5–0.95 under tension. The expression of equivalent plastic strain can also be obtained from Figure 1
ε k p l = ε k i n d k σ k ( 1 d k ) E 0
The value of the equivalent plastic strain calculated by Equation (4) is expected to increase with an increase in the inelastic strain. Additionally, a negative value should not be obtained as the ABAQUS software will consider it as an error and terminate the calculation.

2.2. Energy Loss Method

According to the Najar damage theory [3], the uniaxial damage of brittle solid materials such as concrete can be defined as
d = W 0 W ε W 0
where W 0 = E 0 ε 2 2 denotes the strain energy density function of a non-destructive material and W ε = σ d ε = f ( ε ) d ε denotes the strain energy density function of the material in a damaged state. E 0 denotes the initial elastic modulus of the material. In the Najar damage theory, the strain energy of the damaged state is obtained by linear simplification, W ε = σ ε 2 . However, the damage value calculated using this equation is larger than that of the actual value. When the stress-strain curve equation f ( ε ) of concrete is known, the damage factor of concrete under uniaxial loading is calculated accurately using a numerical integration method.

2.3. Combination of Elastic Modulus Damage and Energy Loss

Ding et al. [4] assumed that the complete curve of concrete under monotonic loading under uniaxial stress was consistent with the skeleton curve under cyclic loading and unloading. Furthermore, deformation of concrete occurs under uniaxial stress. In this case, the strain energy of concrete is denoted as W , and the irreversible deformation after unloading is caused by energy dissipation. The energy consumption is considered as W I . The damage of concrete results in energy dissipation and stiffness reduction. The unloading stiffness of concrete can be defined as ( 1 d 1 ) E 0 and d 1 , which represent the damage variables of the elastic modulus. The energy loss variable d 2 can be defined as
d 2 = W I W = 1 W ε W
where W = σ d ε = f ( ε ) d ε ; W e denotes the recoverable elastic strain energy, and W e = σ 2 2 ( 1 d 1 ) E 0 . To simplify the calculation, it is considered that the damage variables defined are equal in value. Hence, the damage variable can be defined as follows:
d = d 1 = d 2 = 1 σ 2 W E 0

2.4. Energy Equivalence Principle of Sidoroff

According to Sidoroff’s energy equivalence principle, the elastic residual energy generated by the stress acting on the damaged material is similar to that of the elastic residual energy produced by the stress acting on the non-destructive material. Therefore, the equivalent stress or elastic modulus is replaced in the equation with the elastic modulus at the time of damage [5].
The elastic complementary energy of an undamaged material is as follows:
W 0 e = σ 2 2 E 0
The elastic residual energy of the damaged material is as follows:
W d e = σ ˜ 2 E 0 = σ 2 2 E d
Effective stress:
σ ˜ = σ 1 d
Combining the abovementioned equations, we obtain E d = ( 1 d ) 2 E 0 . Hence, the stress-strain relationship of the damaged material can be obtained by
σ = ( 1 d ) 2 E 0 ε
When the stress-strain relationship of concrete under uniaxial loading is known, the expression of the damage factor can be deduced by combining Equation (11).

2.5. Weibull Distribution Considering Damage Threshold

In [9], based on the continuous damage mechanics theory and statistical strength theory, the damage threshold parameter was introduced into the traditional two-parameter Weibull distribution function, and the damage evolution equation was derived as
D = 0 , ( ε < γ ) 1 exp ( ε γ α ) β , ( ε γ )
The damage constitutive relationship of concrete can be obtained by considering the effective stress after damage instead of the nominal stress:
σ = ( 1 D ) E 0 ε
By substituting Equation (13) into Equation (12), the expression of the damage constitutive relationship of concrete can be obtained:
σ = E 0 ε , ( ε < γ ) E 0 ε exp ( ε γ α ) β , ( ε γ )
According to the boundary conditions of the concrete stress-strain curve, the stress-strain value at the beginning of the curve is zero, and the slope at the top of the curve is zero:
β = ε p k γ ε p k ln E 0 ε p k σ p k , α = ε p k γ ε p k γ β ε p k 1 β
In Equation (15), ε p k and σ p k denote the peak strain and stress on the stress-strain curve of concrete, respectively. According to the mechanical tests conducted on ordinary concrete, fiber-reinforced concrete, and ultra-high performance concrete, the elastic limit strain is γ ( 0.7 0.9 ) ε p k . According to [8], the damage variable D in Equation (12) is different from that of the definition of the damage factor of the CDP model parameters in ABAQUS. Therefore, it cannot be directly input into the software. However, the value of D was directly input into the software in [6].

2.6. Parameter Identification Method Based on Test Data

In [7], according to the damage test, the damage factor of ordinary concrete and the normalised inelastic strain follow the first-order exponential decay function:
d k = A 0 e ( ε n o r m i n t 0 ) + B 0
where t , c represent tension and compression respectively, and ε n o r m i n represents the normalised inelastic strain. Since Equation (16) passes through two points (0,0) and (1,1), A 0 , B 0 , and t 0 can be represented by an unknown quantity. If t 0 is used to represent A 0 and B 0 , then A 0 = 1 e 1 t 0 1 , B 0 = 1 e 1 t 0 1 . The physical meaning of t 0 implies the evolution rate of the damage factor with inelastic strain. According to the test data in [7], t 0 = 0.65 is considered appropriate. To establish a general law for damage factor evolution, the relationship table of damage factor, strain ratio, and stress ratio was arranged by assuming that there exists a one-to-one correspondence between the damage factor and stress ratio or strain ratio.

3. Constitutive Relationship of RPC under Uniaxial Tension and Compression

3.1. Stress-Strain Relationship under Uniaxial Compression

Various constitutive equations have been established based on the uniaxial compression tests of ultra-high performance concrete. The ascending part of the stress-strain curve under uniaxial compression can be classified into a straight line [10,11,12], a polynomial curve, a rational fraction curve, and an exponential curve. It is an empirical approach, based on a good description of experimental data. But excepting the first part, without a certain model to sustain them, only the Ax term obeys to HOOK’s law.

3.1.1. Ascending Part of the Curve as a Polynomial

In [13], based on the clear physical quantity of the elastic modulus, the ascending part of the stress-strain curve of RPC under uniaxial compression was defined as a polynomial consisting of only one parameter, while the descending part of the curve was fitted with a rational fraction according to the test data:
y = A x + B x a + C x b , ( 0 x 1 ) x α ( x 1 ) 2 + x , ( x 1 )
where x = ε c / ε c 0 , ε c denotes the strain of concrete under compression, and ε c 0 denotes the peak strain; y = σ c / f c , σ c denotes the compressive stress of concrete and f c denotes the axial compressive strength of the prism. The parameter A = E 0 ε c 0 / f c denotes the ratio of the initial elastic modulus E 0 to the secant modulus E p at the peak point. In the following, x , y , and A have the same meaning, and α is a parameter in the descending section. B = 5 4 A , a = 4 , C = 3 A 4 , b = 5 in [13]; B = 6 5 A , a = 5 , C = 4 A 5 , b = 6 in [14].

3.1.2. Ascending Part of the Curve as a Rational Fraction

Based on the expression of the stress-strain curve of the steel-fiber-reinforced high-strength concrete, [15] simplified the stress-strain curve according to the characteristics of the entire stress-strain curve. After simplification, the ascending part of the stress-strain curve applicable to RPC is given by Equation (18), and the descending part is the same as that in [13].
y = A x x 2 1 + ( A 2 ) x , ( 0 x 1 )
The constitutive equation of ultra-high performance concrete under uniaxial compression proposed in [16] is the same as that in [15], and the descending part is shown in Formula (19), where B = 0.25 ( 0.125 A ( 0.125 A + 0.8 ) 2 0.8 0.2 ) 1.5 .
y = 1 1 + B ( x 1 ) 1.5 , ( x 1 )

3.1.3. Ascending Part of the Curve as an Exponential

Considering the experiment conducted in [17], only the ascending part of the curve was measured, and different functions were used for fitting based on the least squares principle. The fitting effect of the exponential function was found to be the best. Therefore, only the fitting function expression of the ascending part of the curve is given:
y = A x ( 1.000176 0.000176 e x 0.174 )
The above compression constitutive curve was drawn on a graph, as illustrated in Figure 2. In the ascending part of the curve, the larger the parameter value, the steeper the stress-strain curve and the shorter the elastic section; During parameter selection, the proportion, maintenance conditions, and physical characteristic values of the test block, such as peak value, limit value, and elastic modulus should be considered. The descending part of the curve fluctuates significantly because the data are complex to measure during the test and are considerably affected by the test conditions. The actual data of the test should be fitted while selecting the parameters.

3.2. Stress-Strain Relationship under Uniaxial Tension

Based on the uniaxial tensile test data of the RPC test block, in combination with the theoretical analysis, Li et al. [13] uses two equations, which are continuous at the peak point, to describe the two equations respectively:
y = A x + ( 3 2 A ) x 2 + ( A 2 ) x 3 , ( 0 x 1 ) x α ( x 1 ) β + x , ( x 1 )
where x = ε t / ε t 0 , ε t denotes the tensile strain of concrete, and ε t 0 denotes the peak strain; y = σ t / σ t 0 , σ t denotes the tensile stress of concrete and σ t 0 denotes the axial tensile strength of the prism. Parameter A = E 0 ε t 0 / f t denotes the ratio of the initial modulus of elasticity E 0 to the secant modulus E p at the peak point, and α is a parameter in the descending part of the curve, among which α = 5.5 , β = 2.2 .
Based on the experimental results, the stress-strain relationship proposed in [18] is as follows:
σ ( ε ) = f c t ε c a ε , ( 0 ε ε c a ) f c t , ( ε c a ε ε c p )
where ε c a denotes the linear deviation from the initial crack strain, ε c p denotes the ultimate strain, and f c t denotes the average stress in the strain hardening stage.
The authors of [19] provide a micro-scale calculation method for residual tensile strength of ultra-high performance concrete, which can consider fiber distribution, buried depth, and orientation. The tensile stress-strain relationship is as follows:
σ t = E 0 ε , ( ε ε p ) 0.304 f c u ρ f l f d f , ( ε p < ε < ε t lim )
where ε p = f t / E 0 denotes the ultimate elastic strain and ε t lim denotes the ultimate tensile strain, considering 0.01; d f denotes the fiber diameter, ρ f denotes the fiber volume content, and l f denotes the fiber length.
In [20], the expression of the uniaxial tensile stress-strain relationship of ultra-high performance concrete is as follows:
y = A x B x C + 1 , ( 0 x 1 ) x α ( x 1 ) β + x , ( x 1 )
where B = A 1 , C = A A 1 , parameters of the descending part α = 0.025 , β = 1.7 .
The expression of uniaxial tension relation mentioned in [14] is as follows:
y = x 0.92 x 1.09 + 1 , ( 0 x 1 ) x 0.1 ( x 1 ) 2.4 + x , ( x 1 )
Based on the tensile constitutive relationship curve depicted in Figure 3, it can be observed that the ascending part of the curve resembles a straight line, which is in good agreement with the measured data in the literature. Similar to the uniaxial tensile stress-strain constitutive relationship, the larger the value of A, the steeper the ascending part of the curve and the shorter the linear elastic section. The descending part of the curve varies significantly in different literature, reflecting the strain strengthening and softening. The descending part of the curve can be smoothed by adjusting the descent parameters.

4. Experimental Simulation of Steel-Reinforced Reactive Powder Concrete Beams under Bending

4.1. Test Data Reference

The three-point bending test of the SRC-1 steel-reinforced RPC beam in [21] was selected for the numerical simulation. The length of the test beam was 2000 mm, and the section size was B × h = 150 × 200 mm. The effective span was 1800 mm. A Q235-B section steel was used for the built-in I-beam, and the I-beam specification was I14. The section steel length was 1950 mm, which was placed in the middle of the beam, and the thickness of the protective layer was 30 mm. HPB300 grade reinforcement was used for the stirrups, which were spot-welded on the surface of the section steel. The spacing of the stirrups in the pure bending section was 200 mm, and the spacing of the other parts was 100 mm. The cross-section of the test beam is illustrated in Figure 4, and the material performance characteristic values are presented in Table 1.

4.2. Constitutive Relationship of Materials

In general, the brittleness of concrete increases with its strength. However, the addition of fiber into RPC leads to ultra-high strength and high toughness. Therefore, based on the abovementioned uniaxial tension compression constitutive relationship curve of RPC, the gentle descending section of the constitutive relationship curve was selected. In particular, the descending part of the uniaxial tension constitutive relationship curve was selected. If the tensile deformation capacity is poor, it does not match with the ultra-high compressive strength. Therefore, the RPC in the tension zone cracks, and the deformation is significantly large. However, if the RPC in the compression zone does not reach the peak stress, the section steel buckles locally, making the calculation divergent. In this study, the constitutive relationship of uniaxial compression in [13] was selected. The uniaxial tension relationship curve represents an approximate strengthening model, and the constitutive relationship curve mentioned in [14] or [20] can be used. Therefore, the convergence of the calculation is good. In this study, the relationship curve in [14] was selected. According to the measured material parameters in Table 1 and the conversion formula of the RPC basic mechanical property index mentioned in [22], the corresponding peak strain and ultimate strain can be obtained. The tensile constitutive equation is given by Equation (25), and the concrete expression of the compressive constitutive equation is shown in Equation (26):
y = 1.28 x 0.12 x 4 0.16 x 5 , ( 0 x 1 ) x 8 ( x 1 ) 2 + x , ( x 1 )
Comparing the six methods used to calculate the CDP model parameters in ABAQUS with the parameters calculated by Equations (25) and (26), as shown in Figure 5, this study uses the complex Simpson quadrature formula to obtain the damage factor.
The bilinear strengthening model was adopted for the section steel, and the elastic-plastic material properties input into the ABAQUS software are listed in Table 2. The bond-slip constitutive relationship between the section steel and RPC is based on the interface bond constitutive relationship adjusted in [23]; the discrete element method (DEM) is used, as is the bond-slip method in modeling with CDP in ABAQUS. According to the characteristic value of the interface bond strength, the relationship curve of the interface bond stress and slip was obtained, as depicted in Figure 6.

4.3. Establishment of Finite Element Model

The finite element model of the test beam is depicted in Figure 7. A rigid cushion block was set at the loading point to prevent the convergence due to stress concentration. The constraint between the stirrup and RPC was embedded, and the bond-slip was neglected. Spring 2 elements in three directions were used to simulate the interaction between the stirrup and RPC. The spring stiffness in the longitudinal and tangential direction were calculated according to the relationship curve shown in Figure 6, and the force displacement curve of the spring was obtained according to the surface area of each spring. The normal longitudinal direction considered that the two were completely bonded. The stiffness setting was larger, which was equal to the elastic modulus of the section steel or RPC, and the transverse tangential spring stiffness was set to half that of the longitudinal tangential spring stiffness. The elastic modulus of the rigid cushion block was set to 10 times that of the section steel. The rigid cushion block was bound with RPC, and two reference points were established above the cushion block to apply the displacement load. The C3D8R linear reduced integral element was used for the RPC, section steel, and rigid block, whereas the T3D2 element was used for the stirrup simulation.
The mesh size in the cross-section is fixed, while the longitudinal mesh sizes were set to 10mm, 15 mm, 20 mm, 25 mm, 30 mm, 35 mm, 40 mm, 45 mm, 50 mm, respectively. The results calculated by the FE proposed are shown in Figure 8. It can be seen from Figure 8 that the influence of the mesh size on the results is very slight, and the peak load decreases slightly as the mesh size increases. When the mesh size increases from 10 mm to 15 mm, 20 mm, 25 mm, 30 mm, 35 mm, 40 mm, 45 mm and 50 mm, the maximum load decreases by 0.29%, 0.51%, 0.65%, 1.82%, 2.32%, 3.06%, 4.15% and 5.60%, respectively. It should be recognized that mesh sizes that are too small will result in high computational costs, while mesh sizes that are too large will result in lower sensitivity. The mesh size is related to many factors. Therefore, the mesh number and size should be determined according to the specific problem.

4.4. Comparative Analysis of Results

The load-deflection curve, which could appropriately represent the flexural performance of the test beam, was selected, and the finite element simulation results of the tensile constitutive model in [14] with and without bond-slip were obtained, and the simulation results of the tensile constitutive model in [20] considering bond-slip were compared with the test results, as shown in Figure 9 and Table 3. It can be observed from the figure that the results of the tensile constitutive model in [14,20] considering the bond-slip between the steel and RPC were in good agreement with the experimental results. The difference between the simulation results and test values was large, cracking load was large, error was 24.3%, ultimate bearing capacity error was 10.7%, stiffness of the entire simulation process was large, and ultimate displacement error was 10%. Moreover, the convergence of the model was good. However, there was an error between the results with or without the bond-slip and the experimental values due to the difference between the RPC constitutive relationship curve and the actual value. The dilatancy angle, eccentricity, the ratio of initial equibiaxial compressive yield stress to initial uniaxial compressive yield stress, the ratio of the second stress invariant on the tensile meridian, and viscosity coefficient were set as 30°, 0.1, 1.16, 0.6667 and 0.0005, respectively [24,25,26,27,28,29,30,31,32,33,34,35,36]. Thus, the bond-slip constitutive relationship between the section steel and RPC was not accurate.
To verify the consistency between the ultimate load of test and finite element simulation, the ultimate load ratio was calculated, as depicted in Table 3. It could be revealed that the average value was 0.94, the standard deviation was 0.031, and the coefficient of variation was 0.033. The numerical simulation results were in good agreement with the experimental results.
The simulation results revealed that the mode of failure of the beam was bending failure. When the beam was damaged, the value of the compression strain reached 4500, which was on the edge of the compressive zone of RPC in the pure bending section. The stress-strain curve of the maximum longitudinal compressive strain point in the compressive zone of the midspan was extracted, as depicted in Figure 10. When the tensile equivalent plastic strain and maximum principal plastic strain of a point were positive, an initial crack occurred at that point, and the normal vector of the crack surface was parallel to the direction of the maximum plastic strain. As illustrated in Figure 11, the cracks were concentrated in the pure bending section in the middle of the span, forming three main cracks, with one inclined crack in each bending shear section on both sides. The tension flange and the webs of the section steel yielded under tension, and the compression flange yielded under local compression, as depicted in Figure 12. Figure 13 and Figure 14 reflect the tensile and compressive damage degrees, respectively, when the finite element simulation beam approaches the ultimate load.

5. Conclusions

(1)
The energy loss method was suitable for all types of concrete materials and exhibited a high accuracy. Furthermore, the Sidoroff energy equivalence principle demonstrated good applicability. In this study, the energy loss method and Sidoroff energy equivalent principle were used to calculate the CDP model parameters of RPC materials in ABAQUS.
(2)
The ascending part of the constitutive relationship curve of RPC and other ultra-high performance concrete under uniaxial compression adopted polynomial, rational fraction, and exponential forms. There is a minor difference among these forms, and the ascending part parameters have obvious physical meanings. However, the parameters of the descending part of the curve were significantly affected by the mix proportion, test equipment, and maintenance conditions. Therefore, there were obvious differences, which were obtained by fitting the test data. The uniaxial tension constitutive model could be classified into strain-softening, strain-hardening, and approximate strain-hardening types. Owing to the fiber crack resistance and fiber bearing tension, RPC exhibited good toughness and energy dissipation capacity. Thus, it is suggested that when the compressive constitutive model is selected, the value of parameters in the descending part of the curve should result in a gradual decline in the curve. Additionally, when the tensile constitutive model is selected, the approximate strain hardening model should be adopted.
(3)
The numerical simulation results considering the bond-slip were in good agreement with the experimental results with a deviation of less than 10%, which verified the validity of the selected RPC constitutive relationship, and the CDP parameter calculation method was suitable for RPC materials. The numerical simulation performed in this study can provide a reference for the future use of ABAQUS software to simulate the mechanical properties of RPC and ultra-high performance concrete.

Author Contributions

H.L. and J.D. conceived the methodology; H.L. and X.G. wrote the initial draft of the manuscript; H.L. analyzed the data and wrote the final manuscript. In addition, X.G. participated in the revision of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research report was financially supported by the National Natural Science Foundation of China (Grant nos. 51738004, 51878186), by the Innovation Driven Major Projects of Guangxi (Grant no. AA18118055), and by High-level Innovation Team and Outstanding Scholar Plan of Guangxi High Colleges.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The testing and analysis data used to support the findings of this study are included within the article.

Acknowledgments

The authors are very grateful for the support of the above project funds.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Birtel, V.; Mark, P. Parameterized finite element modelling of RC beam shear failure. In Proceedings of the ABAQUS Users’ Conference, Cambridge, UK, 23–25 May 2006; pp. 95–107. [Google Scholar]
  2. Zhang, J.; Wang, Q.; Hu, S. Parameter verification of ABAQUS concrete damage plasticity model. Build. Struct. 2008, 8, 127–130. [Google Scholar] [CrossRef]
  3. Wang, Z.; Yu, Z. Concrete damage model based on energy loss. J. Build. Mater. 2004, 4, 365–369. (In Chinese) [Google Scholar] [CrossRef]
  4. Ding, F.; Yu, Z.; Ou, J. Constitutive model of concrete uniaxial stress damage. J. Chang. Univ. (Nat. Sci. Ed.) 2008, 4, 70–73. (In Chinese) [Google Scholar]
  5. Xiao, Y.; Chen, Z.; Zhou, J.; Leng, Y.; Xia, R. Concrete plastic-damage factor for finite element analysis: Concept, simulation, and experiment. Adv. Mech. Eng. 2017, 9. [Google Scholar] [CrossRef]
  6. Zhu, H. Plastic Damage Analysis of RPC Materials under Uniaxial and Biaxial Compressive Stress States; Beijing Jiaotong University: Beijing, China, 2008. (In Chinese) [Google Scholar]
  7. Li, H. Parameter Identification and Finite Element Verification of Plastic Damage Constitutive Model of RPC Materials; Beijing Jiaotong University: Beijing, China, 2009. (In Chinese) [Google Scholar]
  8. Yang, K.; Sun, L.; Li, W.; Chen, L.; Xie, Z. Sectional design and mechanical properties of RPC hollow compression members. Eng. Mech. 2016, 33, 166–175. (In Chinese) [Google Scholar] [CrossRef]
  9. Xue, Y.; Li, S.; Lin, F.; Xu, H.B. Research on the damage constitutive model of steel fiber concrete considering the influence of damage threshold. Rock Soil Mech. 2009, 30, 1987–1992. [Google Scholar]
  10. Yoo, D.-Y.; Yoon, Y.-S. Structural performance of ultra-high-performance concrete beams with different steel fibers. Eng. Struct. 2015, 102, 409–423. [Google Scholar] [CrossRef]
  11. Chen, S.; Zhang, R.; Jia, L.-J.; Wang, J.-Y. Flexural behaviour of rebar-reinforced ultra-high-performance concrete beams. Mag. Concr. Res. 2018, 70, 997–1015. [Google Scholar] [CrossRef]
  12. Yan, J.-B.; Hu, H.; Wang, T. Flexural behaviours of steel-UHPC-steel sandwich beams with J-hook connectors. J. Constr. Steel Res. 2020, 169, 106014. [Google Scholar] [CrossRef]
  13. Li, L. Research on the Mechanical Properties and Design Methods of Reactive Powder Concrete Beams; Harbin Institute of Technology: Harbin, China, 2010. (In Chinese) [Google Scholar]
  14. Su, J.-Z.; Ma, X.-L.; Chen, B.-C.; Sennah, K. Full-scale bending test and parametric study on a 30-m span prestressed ultra-high performance concrete box girder. Adv. Struct. Eng. 2019, 23, 1276–1289. [Google Scholar] [CrossRef]
  15. Wu, Z.; Shi, C.; He, W.; Wu, L. Effects of steel fiber content and shape on mechanical properties of ultra high performance concrete. Constr. Build. Mater. 2016, 103, 8–14. [Google Scholar] [CrossRef]
  16. Singh, M.; Sheikh, A.H.; Mohamed Ali, M.S.; Visintin, P.; Griffith, M.C. Experimental and numerical study of the flexural behaviour of ultra-high performance fibre reinforced concrete beams. Constr. Build. Mater. 2017, 138, 12–25. [Google Scholar] [CrossRef]
  17. Xu, H.; Deng, Z. Experimental study on flexural performance of prestressed ultra-high performance steel fiber concrete beams. J. Build. Struct. 2014, 35, 58–64. (In Chinese) [Google Scholar] [CrossRef]
  18. Zhang, Z.; Shao, X.; Li, W.; Zhu, P.; Chen, H. Axial tensile performance test of ultra-high performance concrete. China J. Highw. Transp. 2015, 28, 50–58. [Google Scholar] [CrossRef]
  19. Qi, J.; Wang, J.; Ma, Z.J. Flexural response of high-strength steel-ultra-high-performance fiber reinforced concrete beams based on a mesoscale constitutive model: Experiment and theory. Struct. Concr. 2017, 19, 719–734. [Google Scholar] [CrossRef]
  20. Hu, A.; Liang, X.; Yu, J.; Shi, Q.; Li, L. Experimental study on the mechanical properties of ultra-high performance concrete under axial tension. J. Hunan Univ. (Nat. Sci. Ed.) 2018, 45, 30–37. [Google Scholar]
  21. Liu, D. Outsourcing Type Reactive Powder Concrete Beam Bending Performance Test Research; Hunan University: Changsha, China, 2017. (In Chinese) [Google Scholar]
  22. Lv, X.; Wang, Y.; Fu, C.; Zheng, W. Reactive powder concrete basic mechanics performance index value. J. Harbin Inst. Technol. 2014, 46, 1–9. (In Chinese) [Google Scholar]
  23. Bu, L.; Luo, K. Steel outsourcing the stick performance of reactive powder concrete (RPC). Sci. Technol. Eng. 2018, 18, 307–312. [Google Scholar]
  24. Huang, Z. Experimental Study and Optimal Design of UHPC Beam Model; Hunan University: Changsha, China, 2018. (In Chinese) [Google Scholar]
  25. Jian, T. Research on Mechanical Behavior of UHPC Sandwich Floor; Hunan University: Changsha, China, 2019. (In Chinese) [Google Scholar]
  26. Liu, B.; Bai, G.L.; Xu, Z.H.; Ma, J.F.; Han, Y.Y. Experimental study and finite element modelling of bond behavior between recycled aggregate concrete and the shaped steel. Eng. Struct. 2019, 201, 109840. [Google Scholar] [CrossRef]
  27. Bai, L.; Zhang, M.; Yang, L.; Hu, S.; Bao, Z. Experimental study and finite element analysis on bond properties of high ductility steel reinforced cementitious materials. Eng. Mech. 2021, 38, 98–111. [Google Scholar]
  28. ABAQUS Theory Manual; Version 6.13; ABAQUS Inc.: Palo Alto, CA, USA, 2014; Volumes 1–3.
  29. Chou, C.-C.; Wu, S.-C. Cyclic lateral load test and finite element analysis of high-strength concrete-filled steel box columns under high axial compression. Eng. Struct. 2019, 189, 89–99. [Google Scholar] [CrossRef]
  30. Nasiri, E.; Liu, Y. The out-of-plane behaviour of concrete masonry infills bounded by reinforced concrete frames. Eng. Struct. 2019, 184, 406–420. [Google Scholar] [CrossRef]
  31. Farzad, M.; Shafieifar, M.; Azizinamini, A. Experimental and numerical study on bond strength between conventional concrete and Ultra High-Performance Concrete (UHPC). Eng. Struct. 2019, 186, 297–305. [Google Scholar] [CrossRef]
  32. Rabi, M.; Cashell, K.A.; Shamass, R. Flexural analysis and design of stainless steel re-inforced concrete beams. Eng Struct 2019, 198, 109432. [Google Scholar] [CrossRef]
  33. Bypour, M.; Gholhaki, M.; Kioumarsi, M.; Kioumarsi, B. Nonlinear analysis to investigate effect of connection type on behavior of steel plate shear wall in RC frame. Eng. Struct. 2019, 179, 611–624. [Google Scholar] [CrossRef]
  34. Xiao, J. Recycled Aggregate Concrete; China Architecture & Building Press: Beijing, China, 2008. [Google Scholar]
  35. Liu, H.; Wang, Y.; He, M.; Shi, Y.; Waisman, H. Strength and ductility performance of concrete-filled steel tubular columns after long-term service loading. Eng. Struct. 2015, 100, 308–325. [Google Scholar] [CrossRef]
  36. Mirzazadeh, M.M.; Green, M.F. Non-linear finite element analysis of reinforced concrete beams with temperature differentials. Eng. Struct. 2017, 152, 920–933. [Google Scholar] [CrossRef]
Figure 1. Uniaxial compression and tensile behaviour of concrete. (a) Stress-strain relationship of plastic damage under compression; (b) Stress-strain relationship of plastic damage under tension.
Figure 1. Uniaxial compression and tensile behaviour of concrete. (a) Stress-strain relationship of plastic damage under compression; (b) Stress-strain relationship of plastic damage under tension.
Materials 14 04176 g001
Figure 2. Uniaxial compressive stress-strain curve.
Figure 2. Uniaxial compressive stress-strain curve.
Materials 14 04176 g002
Figure 3. Uniaxial tensile stress-strain curve.
Figure 3. Uniaxial tensile stress-strain curve.
Materials 14 04176 g003
Figure 4. Size of beam cross-section (unit: mm).
Figure 4. Size of beam cross-section (unit: mm).
Materials 14 04176 g004
Figure 5. Damage factor and inelastic strain. (a) Compression damage factor inelastic strain. (b) Tension damage factor inelastic strain.
Figure 5. Damage factor and inelastic strain. (a) Compression damage factor inelastic strain. (b) Tension damage factor inelastic strain.
Materials 14 04176 g005
Figure 6. Interface bond strength-slip relationship (Bu Liang-tao, 2018).
Figure 6. Interface bond strength-slip relationship (Bu Liang-tao, 2018).
Materials 14 04176 g006
Figure 7. Finite element model of SRC beam.
Figure 7. Finite element model of SRC beam.
Materials 14 04176 g007
Figure 8. Load-deflection curve of different mesh size.
Figure 8. Load-deflection curve of different mesh size.
Materials 14 04176 g008
Figure 9. Load-deflection curve of numerical simulation and test results.
Figure 9. Load-deflection curve of numerical simulation and test results.
Materials 14 04176 g009
Figure 10. Longitudinal stress-strain at the edge of the midspan compression zone.
Figure 10. Longitudinal stress-strain at the edge of the midspan compression zone.
Materials 14 04176 g010
Figure 11. Maximum plastic strain cloud diagram and crack direction.
Figure 11. Maximum plastic strain cloud diagram and crack direction.
Materials 14 04176 g011
Figure 12. Von-Mises stress diagram of section steel.
Figure 12. Von-Mises stress diagram of section steel.
Materials 14 04176 g012
Figure 13. RPC tensile damage factor cloud diagram.
Figure 13. RPC tensile damage factor cloud diagram.
Materials 14 04176 g013
Figure 14. RPC compression damage factor cloud diagram.
Figure 14. RPC compression damage factor cloud diagram.
Materials 14 04176 g014
Table 1. Material performance parameters (unit: MPa).
Table 1. Material performance parameters (unit: MPa).
Section steelfy246
fu392
Es2.06   ×   105
RPCfcu127.3
fc113.3
fts8.31
Ec4.7   ×   104
Table 2. Parameters of section steel.
Table 2. Parameters of section steel.
Elastic Model/MPaPoisson’s RatioYield Stress/MPaInelastic Strain
206369.10.3246.30
399.50.0169
Table 3. Comparison of finite element simulation results and test results.
Table 3. Comparison of finite element simulation results and test results.
ModelCracking Load/kNYield Load /kNUltimate Load/kNUltimate Deflection/mmUltimate Load Ratio
Test62.2108.013118.13/
Sliding is considered in reference [14]65.3117.1134.819.30.97
Sliding is not considered in reference [14]77.3124.9145.316.50.90
Sliding is considered in reference [20]65.4111.3136.618.30.96
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Guo, X.; Duan, J. Numerical Simulation of Steel-Reinforced Reactive Powder Concrete Beam Based on Bond-Slip. Materials 2021, 14, 4176. https://doi.org/10.3390/ma14154176

AMA Style

Li H, Guo X, Duan J. Numerical Simulation of Steel-Reinforced Reactive Powder Concrete Beam Based on Bond-Slip. Materials. 2021; 14(15):4176. https://doi.org/10.3390/ma14154176

Chicago/Turabian Style

Li, Haoxu, Xiao Guo, and Jiqiang Duan. 2021. "Numerical Simulation of Steel-Reinforced Reactive Powder Concrete Beam Based on Bond-Slip" Materials 14, no. 15: 4176. https://doi.org/10.3390/ma14154176

APA Style

Li, H., Guo, X., & Duan, J. (2021). Numerical Simulation of Steel-Reinforced Reactive Powder Concrete Beam Based on Bond-Slip. Materials, 14(15), 4176. https://doi.org/10.3390/ma14154176

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