Numerical Simulation of Steel Fiber Pull-Out Process Based on Cohesive Zone Model and Uniﬁed Phase-Field Theory

: In steel ﬁber reinforced concrete, the interface is a very complex and weak structure. It is because of the weak interface layer between the steel ﬁber and the matrix that the reinforcing and toughening properties of the steel ﬁber cannot be fully exerted. The interface bond performance is the core of the meso-mechanical properties of steel ﬁber reinforced concrete. To study its inﬂuence on the mechanical properties of steel ﬁber reinforced concrete, three-phase ﬁnite element models of steel ﬁber pull-out are established based on the cohesive zone model and uniﬁed phase-ﬁeld theory by means of FEM in this paper. The interface bond is simulated by a zero-thickness cohesive element, and the pull-out process of steel ﬁber in the concrete matrix is analyzed to provide a basis for the fracture research of steel ﬁber reinforced concrete. In this paper, the inﬂuence of factors such as the embedment depth, length–diameter ratio, embedment angle, and interface properties of steel ﬁbers on the pull-out mechanical properties of steel ﬁbers are considered, and the relevant ﬁnite element models are established to conduct numerical simulations of the pull-out process of steel ﬁbers. The numerical simulation results are in good agreement with the experimental results, and this veriﬁes the reliability of the model. The results show that the steel ﬁber pull-out ﬁnite element model established by the cohesive zone model and phase-ﬁeld regularized cohesive zone model (PF-CZM) has a certain reliability; the peak pull-out load of the steel ﬁber increases with an increase in the embedment depth of the steel ﬁber, and decreases with an increase in the length–diameter ratio and embedment angle of the steel ﬁber; by controlling the strength of the interface layer and the concrete matrix, the reinforcement effect of the steel ﬁber on the concrete matrix can be improved; and the PF-CZM has a good characterization of the damage and failure evolution process of the concrete matrix.


Introduction
As a widely used construction material, concrete has high compressive strength and good durability.However, with new developments over time, higher requirements are put forward for engineering materials, and these requirements make the application of ordinary concrete somewhat limited because of its low tensile strength and high brittleness [1].The incorporation of steel fibers in the cement matrix can significantly alleviate these deficiencies.As a composite material with excellent performance, not only does steel fiber reinforced concrete inherit the high compressive strength of ordinary concrete [2], but steel fibers also inhibit the development of cracks, so that the ductility and toughness of steel fiber reinforced concrete structures are enhanced when the structures are damaged [3,4].Therefore, steel fiber reinforced concrete is increasingly used in most buildings, roads, bridges, tunnels, dams, and other important civil engineering structures [5].
However, due to the presence of the interface layer between steel fibers and the concrete matrix, the reinforcing and toughening properties of steel fiber on the concrete matrix cannot be fully exploited.As a key factor for steel fibers to play a role in the reinforcing and toughening properties of the concrete matrix, the interface bond performance is generally measured by steel fiber pull-out tests.Naaman et al. [6] conducted pull-out tests on single fibers around different shapes of steel fibers, different strengths of the cement matrix, and different diameters and embedment depths of steel fibers in order to investigate the relationship between interfacial bond shear stress and bond slip.The results show that with an increase in the strength of the matrix, the bond between the fiber and the matrix increases correspondingly, and the fiber diameter and fiber embedment length do not affect the characteristics of the pull-out load versus end slip relationship.Tuyan M et al. [7] conducted a steel fiber pull-out test on slurry infiltrated fiber concrete (SIFCON) and investigated the influence of the length-diameter ratio of steel fiber, the embedment depth, and the fiber type on the pull-out behavior of single fibers in the SIFCON matrix.The results show that the peak pull-out load of fibers increases with an increase in fiber embedment depth and fiber diameter, and it was found that end-hooked fibers have better interface bond performance than straight fibers.Bhutta A et al. [8] compared three different shapes of steel fibers pulled out of the concrete matrix through their experiments, and concluded that end-hooked steel fibers have the best bond performance compared with straight and wavy steel fibers.Shannag et al. [9] designed a single fiber pull-out device to investigate the interface bond properties and the reinforcing and toughening mechanism of steel fibers when they are pulled out from the concrete matrix.The test results show that the dense small particle matrix can improve the interface bond performance better than the traditional mortar matrix and that increasing the embedment depth of steel fibers and the volume fraction of steel fibers can improve the pull-out load-bearing capacity of steel fibers.Abu-Lebdeh et al. [10,11] discussed the influence of loading rates on pull-out behavior by conducting a test of steel fiber-matrix bond characteristics; their results show that the peak pull-out load increases with an increase in the pull-out rate, and they also investigated the influence of matrix strength on the pull-out behavior of steel fiber reinforced concrete and found that for all unbroken deformed fibers, the maximum pull-out load and total pull-out energy increased with an increase in matrix strength.Yoo et al. [12] investigated the influence of fiber type and fiber spacing on the pull-out behavior of steel fibers in ultra-high-performance concrete (UHPC).The results show that for straight steel fibers, the matrix damage is small and does not affect the pull-out performance of adjacent fibers; however, for end-hooked steel fibers, the matrix displays severe damage, and as the distance between the fibers decreases, a larger spalling area occurs between them.Ye J D et al. [13] investigated the bond properties of the helical steel fiber-UHPC matrix interface by conducting single pull-out tests rectangular helical steel fibers with different embedded depths.Lee et al. [14] investigated the influence of the steel fiber embedment angle on the bond properties between steel fibers and the matrix through a pull-out test of steel fibers with five different embedment angles, namely 0 • , 15 • , 30 • , 45 • , and 60 • .The authors of [15,16] also found that the peak pull-out load of steel fibers decreases with an increase in the inclined angle, and the peak displacement increases accordingly.
Single experiments to study the pull-out process of steel fibers often have certain limitations.In recent years, the numerical simulation method has gradually become an important direction in civil engineering research [17,18], and using numerical simulation methods to simulate the pull-out process of steel fibers has become an important means to study the mechanical properties of the interface [19].Numerical simulation can effectively save costs and time for making specimens, and the process of fiber pull-out and concrete damage can be better observed at the meso-level.Hsueh et al. [20] and Kim et al. [21] regarded the contact surface between the fiber and the concrete matrix as an interface layer of extremely thin thickness, and analyzed the interface's shear stress, but their analysis results do not show the relationship between the fiber pulling displacement and the interface debonding length.Huang et al. [22] then established three-phase and four-phase finite element analysis models on this basis, and their research results show that the threephase model is more practical.Liu F [23] established three-phase finite element models based on the cohesive zone model and the concrete damage plasticity (CDP) model and simulated several tests, which verified the correctness and feasibility of using the cohesive zone model to simulate the interface bond performance.Zong Y W et al. [24] established a three-phase finite element model for steel fiber pull-out and analyzed the influence of different factors such as the steel fiber embedment depth, interface layer thickness, and steel fiber length-diameter ratio on the pull-out characteristics of steel fiber.Yang C C et al. [25] found that the fiber falling off due to the failure of the matrix itself is the main failure form of steel fiber pull-out according to a real bending-pulling failure fracture electron microscope scan image, and established two pull-out models based on the cohesive zone model to analyze the influence of different strengths of the interface layer and the matrix itself on the pull-out properties of the fiber.
According to the large number of numerical simulation works on fiber pull-out, it is known that the use of the cohesive zone model can simplify the analysis of the complex interfacial separation process, and the use of the cohesive zone model to simulate the interfacial interaction has been verified to have certain accuracy and simplicity.Based on this, in the finite element model of steel fiber pull-out established in this paper, the interface layer between steel fiber and the concrete matrix is simulated with zero-thickness cohesive elements to simulate the bond performance between the interfaces.In addition, the selection of the constitutive model of the concrete matrix has an important influence on the accurate simulation of the steel fiber pull-out process, but most scholars assume that the concrete matrix is undamaged and is represented by linear elasticity, which appears to be inconsistent with the real situation.Wu [26] proposed a unified phase-field theory based on the theoretical framework of thermodynamics and subsequently proposed the phase-field regularized cohesive zone model (PF-CZM) applicable to quasi-brittle materials.This model provides great convenience for studying the generation and evolution of concrete cracks and is very suitable for simulating the damage of concrete.Zhuang L J [27] verified the applicability of the PF-CZM in a simulation analysis of the whole concrete damage process by simulating several experiments with the PF-CZM.On this basis, Qiu J F et al. [28] discussed the similarities of, differences between, and advantages of the PF-CZM and the extended finite element method (XFEM) for simulating the whole process of concrete nonlinear damage.Li B B et al. [29] investigated the meso-mechanical model of plain concrete based on the PF-CZM and confirmed the feasibility of the PF-CZM in heterogeneous materials.Considering the performance of the PF-CZM in the numerical simulation of concrete damage, the concrete matrix part of the finite element model of steel fiber pull-out established in this paper is represented by this damage model.
In summary, a large number of engineering practices and studies have shown that adding an appropriate number of steel fibers into concrete will have a significant strengthening and toughening effect, but the contact area between steel fibers and concrete is still a weak area where damage and fracture occur in steel fiber reinforced concrete; this has a decisive influence on the load-bearing capacity of the structure.Pull-out experiments on concrete specimens containing single steel fibers and their numerical simulation are important means to investigate the interface bonding properties between steel fibers and concrete.This paper combines the latest unified phase-field theory [30] to establish a finite element model for a steel fiber pull-out test with a cohesive zone model for the interface part and the PM-CZM for the concrete matrix part, in order to analyze the influences of different factors and different failure forms on the steel fiber pull-out process and to investigate the interface bonding properties of steel fiber reinforced concrete, so as to provide a relevant theoretical basis for practical engineering and application.
The remainder of this paper is outlined as follows.The theoretical basis and main equations related to the cohesive zone model and the unified phase-field theory used in this work are recalled in Section 2. In Section 3, the reliability of the combination of the cohesive model and the PF-CZM is verified using a simulation of the established numerical calculation model of steel fiber pull-out.In Section 4, the influence of different factors on the pull-out process of steel fiber is systematically studied.The most relevant conclusions are drawn in Section 5, closing this paper.

Cohesive Zone Model
In this work, the zero-thickness layer of the cohesive element is embedded between the steel fiber and the concrete matrix to simulate the interface bond.As a functional element, the cohesive element is located between two adjacent elements, and the use of a reasonable constitutive model is key to describing the interface damage accurately and effectively.The commonly used cohesive zone models include bilinear, trapezoidal, exponential, and polynomial forms [31], and the bilinear cohesive zone model is used in this paper.The bilinear cohesive zone model is the most effective and commonly used model for characterizing the interface between fibers and the concrete matrix.Taking the normal direction as an example, the model's constitutive curve is shown in Figure 1.

Cohesive Zone Model
In this work, the zero-thickness layer of the cohesive element is embedded between the steel fiber and the concrete matrix to simulate the interface bond.As a functional element, the cohesive element is located between two adjacent elements, and the use of a reasonable constitutive model is key to describing the interface damage accurately and effectively.The commonly used cohesive zone models include bilinear, trapezoidal, exponential, and polynomial forms [31], and the bilinear cohesive zone model is used in this paper.The bilinear cohesive zone model is the most effective and commonly used model for characterizing the interface between fibers and the concrete matrix.Taking the normal direction as an example, the model's constitutive curve is shown in Figure 1.   is the critical bond strength, representing the maximum bond force that the interface can withstand before separation; max  is the critical displacement of the element at the time of failure, representing the maximum separation distance that can transmit the interface bond force, when max    , which means that the interface is completely debonding and the elements on both sides are separated; the area intersected by the   curve and the horizontal axis is the interface fracture energy C G ; and K is the tangent stiffness of the element before its degradation.The interface fracture energy C G , the bond strength  , and the displacement  are three important parameters of the cohesive zone model, and the relationship between them can be expressed as The bilinear cohesive element has its judgment basis from damage initiation to damage evolution.The damage initiation criterion is used to determine when the damage begins to occur and when the stiffness begins to decay.In this paper, the maximum nominal stress criterion is selected as the initial initiation criterion, and the initial damage is satisfied as follows: where n  and t  represent the actual stress values of the cohesive element perpendicular to and parallel to the crack direction, respectively; and 0 n  and 0 t  represent the maximum stress that can be borne perpendicular to the crack and parallel to the crack, respectively.
After the initial damage occurs, the damage evolution process follows.The damage evolution criterion in the cohesive zone model is used to express the relationship between the separation displacement and the traction force, and the degradation of its mechanical properties is characterized by stiffness degradation.Damage evolution can be expressed as τ max is the critical bond strength, representing the maximum bond force that the interface can withstand before separation; δ max is the critical displacement of the element at the time of failure, representing the maximum separation distance that can transmit the interface bond force, when δ > δ max , which means that the interface is completely debonding and the elements on both sides are separated; the area intersected by the τ ∼ δ curve and the horizontal axis is the interface fracture energy G C ; and K is the tangent stiffness of the element before its degradation.The interface fracture energy G C , the bond strength τ, and the displacement δ are three important parameters of the cohesive zone model, and the relationship between them can be expressed as The bilinear cohesive element has its judgment basis from damage initiation to damage evolution.The damage initiation criterion is used to determine when the damage begins to occur and when the stiffness begins to decay.In this paper, the maximum nominal stress criterion is selected as the initial initiation criterion, and the initial damage is satisfied as follows: where τ n and τ t represent the actual stress values of the cohesive element perpendicular to and parallel to the crack direction, respectively; and τ 0 n and τ 0 t represent the maximum stress that can be borne perpendicular to the crack and parallel to the crack, respectively.
After the initial damage occurs, the damage evolution process follows.The damage evolution criterion in the cohesive zone model is used to express the relationship between the separation displacement and the traction force, and the degradation of its mechanical properties is characterized by stiffness degradation.Damage evolution can be expressed as The stiffness degradation is represented by D, and its value range is 0~1.If D = 0, it means that no damage has occurred; if D = 1, it means complete damage has occurred.
After damage occurs, the stiffness of the element during the degradation is represented by K 0 :

Unified Phase-Field Theory
According to the existing research [25], in the process of steel fiber pull-out, not only does interface damage between the steel fiber and the concrete matrix lead to the fiber falling off, but the fiber falling off due to damage to the concrete matrix is the main damage form.To more intuitively and accurately characterize concrete matrix damage, this paper adopts the phase-field regularized cohesive zone model (PF-CZM), which was proposed based on the unified phase-field theory.The basic idea and main equations of the unified phase-field theory are briefly introduced below [26,30,32,33].
The reference configuration of a cracking solid, Ω ⊂ R n dim (n dim = 1, 2, 3), with a set of sharp cracks, S, is shown in Figure 2a.In addition, the material particles of the solid are labeled by their spatial coordinates, x; thus, u(x) represents the displacement field of the solid and ε(x) := ∇ sym u(x) represents the strain field of the solid.The solid boundary is divided into two parts, ∂Ω u and ∂Ω t , which satisfies ∂Ω u ∩ ∂Ω t = ∅ and ∂Ω u ∪ ∂Ω t = ∂Ω.There is a displacement boundary condition, u * (x), at x ∈ ∂Ω u and a surface tractions boundary condition, t * (x), at x ∈ ∂Ω t .Thus, the displacement space U u can be expressed as The stiffness degradation is represented by D , and its value range is 0~1.If 0 D  , it means that no damage has occurred; if 1 D  , it means complete damage has occurred.After damage occurs, the stiffness of the element during the degradation is represented by 0

Unified Phase-Field Theory
According to the existing research [25], in the process of steel fiber pull-out, not only does interface damage between the steel fiber and the concrete matrix lead to the fiber falling off, but the fiber falling off due to damage to the concrete matrix is the main damage form.To more intuitively and accurately characterize concrete matrix damage, this paper adopts the phase-field regularized cohesive zone model (PF-CZM), which was proposed based on the unified phase-field theory.The basic idea and main equations of the unified phase-field theory are briefly introduced below [26,30,32,33].
The reference configuration of a cracking solid, After the cracking solid is stressed, once the requirements of the cracking criterion are met, the sharp crack  will occur in the solid.As shown in Figure 2b, in the phase- After the cracking solid is stressed, once the requirements of the cracking criterion are met, the sharp crack S will occur in the solid.As shown in Figure 2b, in the phase-field method, the sharp crack S will be regularized over a localization band, B ⊆ Ω, there is a continuous damage field, d(x), in the localization band, and there is no damage to the area outside the localization band.The damage space U d can be expressed as where the derivative .d(x) ≥ 0 of damage to time represents the irreversibility of the damage.According to the first and second laws of thermodynamics, the energy dissipation inequality of the solid deformation process is as follows:

Constitutive Relation and Damage Evolution Law
.
where ψ represents the local free energy density of the solid, and considering the simplest isotropic case, the local free energy density ψ can be expressed as the following general form: where ω(d) is the energy degradation function; ψ 0 (ε) is the initial strain energy density; σ = E 0 : ε is the effective stress tensor; and E 0 and C 0 are the elastic stiffness and compliance, respectively.At this time, the driving force Y of the damage is given by where Y := ∂ψ/∂ω is the effective damage driving force.In summary, the following constitutive relation can be obtained: (2) Damage evolution law Following Griffith's fracture mechanics theory, the sharp cracks can be expressed as regularized cracks within the localization band B: According to the law of conservation of energy, there is the following relationship between the damage evolution of the localization band and the original sharp crack: where G f is the fracture energy, usually regarded as a material property; γ(d) is the crack surface density function, expressed as follows: where the length scale b characterizes the width of the localization band B, which is an important parameter of the unified phase-field damage theory, and the geometric crack function α(d) ∈ [0, 1] determines the damage distribution in the localization band B.
Based on the irreversibility of damage and Gauss's divergence theorem, the damage evolution law can be obtained and its Neumann-type boundary conditions are as follows: where n B is the outward normal of the boundary ∂B; and δ d γ is the derivative of the crack surface density function shown in Equation ( 14).

Phase-Field Regularized Cohesive Zone Model
There are two important characteristic functions in the unified phase-field damage theory, which are the geometric crack function α(d), used to normalize the topology of physical cracks, and the energy degradation function ω(d), used to reflect the softening process of the material.Wu [26,30,32,33] combined the unified phase-field theory with the cohesive zone model, clarified the suggested values of the geometric crack function α(d) and the energy degradation function ω(d) in practical applications, and then proposed the PF-CZM.The specific suggested forms are as follows: 16) where the coefficients a 1 , a 2 , and a 3 are derived from the following material: where k 0 is the initial slope, b is the length scale, and ûu is the ultimate crack opening.When p and a i (i = 1, 2, 3) take specific values, the specific softening criterion can be fitted.The Cornelissen softening criterion [26,30,32,33] commonly used in concrete is adopted in this paper.It is calibrated according to the relevant test data of concrete and is most suitable for concrete materials.At this time, p = 2, a 2 = 1.3868, and a 3 = 0.6567.

Numerical Implementation (1) Governing equations
Based on the constitutive relations and damage evolution laws, the governing equations in strong form can be expressed as and the following is the Neumann-type boundary condition: Using the weighted residual method, the governing equations and boundary conditions can be transformed into the following weak forms: where δP is the virtual work by the external force, and δγ is the variation of the crack surface density function.The specific expression is as follows: In conclusion, the unified phase-field damage theory starts from thermodynamics and transforms the complex structural damage problem and fracture mechanics problem into two governing equations, one for solving the displacement field problem and the other for solving the damage phase-field problem.
(2) Finite element discretization In the two-dimensional (2D) case, the calculation domain Ω can be divided into meshes, M h , for discretization, and the superscript h represents the mesh size of the finite element domain M h .It is worth noting that the mesh size h should be much smaller than the localization bandwidth b, which is generally taken as h/b ≤ 1/5 to ensure the accuracy of the calculation results.Each element node has three nodal degrees of freedom (dofs): two degrees of freedom of displacement and one degree of freedom of phase-field damage.
Based on the above settings, the displacement field u h and the strain field ε h of the structure can be approximated by interpolation of the nodal displacements χ := {χ 1 , . . . ,χ i , . ..}T : In the 2D case, the interpolation function matrix N i (x) and the displacement-strain matrix B i (x) at node i are the following, respectively: where N i (x) is the interpolation function at node i.Similar to the discretization of the displacement field, the damage field d h and the damage gradient ∇d h can be approximated by the interpolation of the node damage χ := { χ 1 , . . . ,χ i , . ..}T : Accordingly, the interpolation function matrix N i (x) and the displacement-strain matrix B i (x) at node i are the following, respectively: With the above finite element discretization, the weak form (Equation ( 24)) can be transformed into a residual form of the equilibrium equation: (3) Solving algorithm To solve the nonlinear algebraic equations system shown in Equation (29), this paper's research is implemented using a user-defined subroutine.The user-defined subroutine is based on the corresponding interface of the finite element software, and the user writes the code according to the Fortran language.Based on the implementation of the user-defined subroutine [32], the PF-CZM achieves the simulation of the whole process of structural damage and destruction using the alternate minimization (AM) algorithm [34,35] of the above two control equations.The specific solution process is as follows (Figure 3): To solve the nonlinear algebraic equations system shown in Equation (29), this p per's research is implemented using a user-defined subroutine.The user-defined subro tine is based on the corresponding interface of the finite element software, and the us writes the code according to the Fortran language.Based on the implementation of t user-defined subroutine [32], the PF-CZM achieves the simulation of the whole process structural damage and destruction using the alternate minimization (AM) algorithm [34, of the above two control equations.The specific solution process is as follows (Figure 3):

Model Building
In this section, numerical simulations are carried out on the experiments conduct by Shannag et al. [9] to investigate the pull-out process of steel fiber; the schematic d gram and dimensions of the steel fiber pull-out specimen are shown in Figure 4a.The s of the concrete matrix is 23 mm × 25 mm, and the fiber diameter is 0.19 mm.In this paper, ABAQUS finite element software was used to establish the numerical calculation model.For this finite element model, firstly, two parts of steel fiber and a concrete matrix are established, which are respectively assigned to different sets and assembled into one part.The concrete matrix part adopts the PF-CZM, and the CPS4T element is selected; the CPS4T is a four-node plane stress full-integration element with both dis- In this paper, ABAQUS finite element software was used to establish the numerical calculation model.For this finite element model, firstly, two parts of steel fiber and a concrete matrix are established, which are respectively assigned to different sets and assembled into one part.The concrete matrix part adopts the PF-CZM, and the CPS4T element is selected; the CPS4T is a four-node plane stress full-integration element with both displacement and temperature dofs built in; the steel fiber part adopts the bilinear cohesive zone model, and the CPS4R element is selected; the CPS4R is a four-node plane stress reduced-integration element, and when the steel fiber is bent and deformed, the linear reduced-integration element is selected to avoid the shear self-locking phenomenon.Then, the zero-thickness cohesive element is embedded between the two sets to represent the interface bond between the steel fiber and the concrete matrix.The cohesive element adopts the bilinear cohesive zone model, and the COH2D4 element is selected; the COH2D4 is the 2D cohesive element.
The steel fiber is subjected to a vertical downward concentrated force and restrained to vertical displacement on the upper side of the concrete matrix and horizontal displacement on both sides of the concrete matrix.The loading method adopts displacement loading, and the loading and restraint positions are shown in the 2D steel fiber pull-out finite element model shown in Figure 4b.The incorporated length scale parameter is b = 1 mm and the edge length of the quadrangle is h = 0.5 mm, and the 2D finite element mesh consists of 14,528 nodes and 14,211 quadrilateral linear elements.
The material parameters of the steel fiber and the concrete matrix were obtained from ref. [9], and the material parameters are shown in Table 1.The two parameters of the bond strength and fracture energy of cohesive elements in the interface layer are unknown in the finite element simulation, the parameters of the cohesive elements are determined using the trial calculation method, and the optimal parameters are determined by adjusting the two parameters of bond strength and fracture energy and comparing the pull-out load-displacement curves obtained by the numerical simulation with the experimental results.To investigate the influence of different factors on the mechanical properties of steel fiber pull-out during the steel fiber pull-out process and to consider the calculation efficiency of the model, the finite element model shown in Figure 5 is established in this paper.The size of the concrete matrix is 12.7 mm × 25.4 mm, and the fiber diameter is 0.5 mm.The model establishment, boundary conditions, and loading methods are the same as above.The finite element model of steel fiber pull-out vertical to the concrete matrix shown in Figure 5a is used to simulate and analyze the influence of different embedment depths, length-diameter ratios, and interface properties of steel fibers on the pull-out mechanical properties of steel fibers.The influence of different embedment angles of steel fibers on the pull-out mechanical properties of steel fibers is simulated and analyzed using the finite element model of steel fiber pull-out inclined to the concrete matrix shown in Figure 5b.
as above.The finite element model of steel fiber pull-out vertical to the concrete matrix shown in Figure 5a is used to simulate and analyze the influence of different embedment depths, length-diameter ratios, and interface properties of steel fibers on the pull-out mechanical properties of steel fibers.The influence of different embedment angles of steel fibers on the pull-out mechanical properties of steel fibers is simulated and analyzed using the finite element model of steel fiber pull-out inclined to the concrete matrix shown in Figure 5b.

Model Validation and Analysis
Figure 6 shows the comparison between the pull-out load-displacement curve obtained from numerical simulation and the results of the steel fiber pull-out test in ref. [9], where the vertical axis is the vertical downward pull-out load applied to the steel fiber, and the horizontal axis is the relative displacement value of the steel fiber and the concrete matrix.When the bond strength  of the cohesive element is taken as 0.5 MPa and the fracture energy C G is taken as 0.25 N/mm, the peak pull-out load obtained by the simulation is 10.148 N, the displacement corresponding to the peak pull-out load is 0.0187 mm, and the experiment results are 10.183N and 0.0181 mm, with errors of 0.4% and 3.3%, respectively.Thus, the numerical simulation results are in good agreement with the experiment results.
It can be seen from Figure 6 that the numerical simulation results in this paper are in good agreement with the experiment results.The whole process of steel fiber pull-out can be roughly divided into four stages: in the first stage, the steel fiber is pulled at the beginning, and the pull-out load and the displacement grow linearly, while the steel fiber and the concrete matrix are in an elastic bond state; in the second stage, the increase in the pull-out load, the shear stress between the steel fiber and the concrete matrix reaches the interface bond strength, and the structure begins debonding, at which time the pull-out load and the displacement grow nonlinearly and increase to the maximum value; in the third stage, as the pull-out load continues to increase, the bonding area disappears, the steel fiber is completely debonding, and the occlusal force between the steel fiber and the concrete matrix gradually decreases; in the fourth stage, the pull-out load is completely borne by friction at this time and gradually decreases as the steel fiber is pulled out.

Model Validation and Analysis
Figure 6 shows the comparison between the pull-out load-displacement curve obtained from numerical simulation and the results of the steel fiber pull-out test in ref. [9], where the vertical axis is the vertical downward pull-out load applied to the steel fiber, and the horizontal axis is the relative displacement value of the steel fiber and the concrete matrix.When the bond strength τ of the cohesive element is taken as 0.5 MPa and the fracture energy G C is taken as 0.25 N/mm, the peak pull-out load obtained by the simulation is 10.148 N, the displacement corresponding to the peak pull-out load is 0.0187 mm, and the experiment results are 10.183N and 0.0181 mm, with errors of 0.4% and 3.3%, respectively.Thus, the numerical simulation results are in good agreement with the experiment results.From this analysis, it can be seen that the finite element model established in this paper can better reflect the whole process of pulling out the steel fiber from the concrete matrix, which also shows that the model established in this paper has a certain accuracy and reliability.

Effect of Steel Fiber Embedment Depth
To investigate the influence of steel fiber embedment depth on the pull-out mechanical properties of steel fiber, finite element models with different steel fiber embedment depths of 6, 8, and 10 mm are established to simulate the pull-out process of steel fibers, as shown in Figure 5a.Other parameters are the same as above, and Figures 7 and 8 show the pull-out load-displacement curves, peak pull-out load, and the displacement corresponding to the peak pull-out load for steel fibers with different embedment depths, respectively.It can be seen from Figure 6 that the numerical simulation results in this paper are in good agreement with the experiment results.The whole process of steel fiber pull-out can be roughly divided into four stages: in the first stage, the steel fiber is pulled at the beginning, and the pull-out load and the displacement grow linearly, while the steel fiber and the concrete matrix are in an elastic bond state; in the second stage, the increase in the pull-out load, the shear stress between the steel fiber and the concrete matrix reaches the interface bond strength, and the structure begins debonding, at which time the pull-out load and the displacement grow nonlinearly and increase to the maximum value; in the third stage, as the pull-out load continues to increase, the bonding area disappears, the steel fiber is completely debonding, and the occlusal force between the steel fiber and the concrete matrix gradually decreases; in the fourth stage, the pull-out load is completely borne by friction at this time and gradually decreases as the steel fiber is pulled out.
From this analysis, it can be seen that the finite element model established in this paper can better reflect the whole process of pulling out the steel fiber from the concrete matrix, which also shows that the model established in this paper has a certain accuracy and reliability.

Effects of Different Factors on the Mechanical Properties of Steel Fiber Pull-Out 4.1. Effect of Steel Fiber Embedment Depth
To investigate the influence of steel fiber embedment depth on the pull-out mechanical properties of steel fiber, finite element models with different steel fiber embedment depths of 6, 8, and 10 mm are established to simulate the pull-out process of steel fibers, as shown in Figure 5a.Other parameters are the same as above, and Figures 7 and 8 show the pull-out load-displacement curves, peak pull-out load, and the displacement corresponding to the peak pull-out load for steel fibers with different embedment depths, respectively.
From this analysis, it can be seen that the finite element model established in this paper can better reflect the whole process of pulling out the steel fiber from the concrete matrix, which also shows that the model established in this paper has a certain accuracy and reliability.

Effect of Steel Fiber Embedment Depth
To investigate the influence of steel fiber embedment depth on the pull-out mechanical properties of steel fiber, finite element models with different steel fiber embedment depths of 6, 8, and 10 mm are established to simulate the pull-out process of steel fibers, as shown in Figure 5a.Other parameters are the same as above, and Figures 7 and 8 show the pull-out load-displacement curves, peak pull-out load, and the displacement corresponding to the peak pull-out load for steel fibers with different embedment depths, respectively.The results show that the distribution of the three groups of pull-out load-displacement curves is basically the same, and the ascending sections partially overlap, indicating that the embedment depth does not work before the steel fibers are debonding.With the increase in embedment depth, the contact area between the steel fiber and the concrete matrix increases, so the peak pull-out load also increases accordingly.The pull-out loadbearing capacities of steel fibers with embedment depths of 8 mm and 10 mm are 30.7% and 61.37%, respectively, these are higher than that of steel fibers with an embedment depth of 6 mm, which indicates that with the increase in the embedment depth of steel fiber, the pull-out load-bearing capacity of steel fiber increases greatly, and the displacement corresponding to the peak pull-out load also has a small increase.Therefore, the embedment depth of steel fiber significantly affects the mechanical properties of steel fiber pull-out.

Effect of Length-Diameter Ratio of Steel Fiber
The length and diameter of steel fibers are important dimensional parameters.To investigate the influence of the length-diameter ratio of the steel fiber on the mechanical The results show that the distribution of the three groups of pull-out load-displacement curves is basically the same, and the ascending sections partially overlap, indicating that the embedment depth does not work before the steel fibers are debonding.With the increase in embedment depth, the contact area between the steel fiber and the concrete matrix increases, so the peak pull-out load also increases accordingly.The pull-out load-bearing capacities of steel fibers with embedment depths of 8 mm and 10 mm are 30.7% and 61.37%, respectively, these are higher than that of steel fibers with an embedment depth of 6 mm, which indicates that with the increase in the embedment depth of steel fiber, the pull-out load-bearing capacity of steel fiber increases greatly, and the displacement corresponding to the peak pull-out load also has a small increase.Therefore, the embedment depth of steel fiber significantly affects the mechanical properties of steel fiber pull-out.

Effect of Length-Diameter Ratio of Steel Fiber
The length and diameter of steel fibers are important dimensional parameters.To investigate the influence of the length-diameter ratio of the steel fiber on the mechanical properties of steel fiber pull-out, finite element models with steel fiber diameters of 0.5 mm, 0.2 mm, and 0.1 mm and an embedment depth of 6 mm are established to simulate the steel fiber pull-out process, as shown in Figure 5a.With all other parameters unchanged, the peak pull-out load and the displacement corresponding to the peak pull-out load for steel fibers with different length-diameter ratios are shown in Figure 9.
and 61.37%, respectively, these are higher than that of steel fibers with an embedment depth of 6 mm, which indicates that with the increase in the embedment depth of steel fiber, the pull-out load-bearing capacity of steel fiber increases greatly, and the displacement corresponding to the peak pull-out load also has a small increase.Therefore, the embedment depth of steel fiber significantly affects the mechanical properties of steel fiber pull-out.

Effect of Length-Diameter Ratio of Steel Fiber
The length and diameter of steel fibers are important dimensional parameters.To investigate the influence of the length-diameter ratio of the steel fiber on the mechanical properties of steel fiber pull-out, finite element models with steel fiber diameters of 0.5 mm, 0.2 mm, and 0.1 mm and an embedment depth of 6 mm are established to simulate the steel fiber pull-out process, as shown in Figure 5a.With all other parameters unchanged, the peak pull-out load and the displacement corresponding to the peak pull-out load for steel fibers with different length-diameter ratios are shown in Figure 9.It can be seen that with the same steel fiber embedment depth, the larger the diameter, the greater the pull-out force required to pull out the steel fibers.This also indicates that thicker steel fibers provide better reinforcement to steel fiber reinforced concrete.Therefore, steel fiber reinforced concrete with better mechanical properties can be obtained by using thick steel fiber in engineering practices.However, it is not the case that thicker steel fibers are always better; steel fibers that are too thick are not conducive to structural load-bearing and would affect the ductility and toughness of the structure [36].

Effect of Steel Fiber Embedment Angle
Steel fibers are randomly distributed in concrete materials, so they are mixed in concrete at different angles in practical engineering.To study the influence of the embedment angle of steel fibers on the mechanical properties of steel fiber pull-out, finite element models with different embedment angles of steel fibers, such as 0 • , 15 • , 30 • , 45 • , and 60 • , are established and simulated as shown in Figure 5b.The displacement loads at all embedment angles are perpendicular to the concrete matrix surface and the embedment depth of the steel fibers is 10 mm.
The results shown in Figure 10 were obtained with only the embedment angle of the steel fibers changed.The simulation results show the following: (1) the magnitude of the peak pull-out load decreases with the increase in the embedment angle of the steel fiber; (2) the displacement corresponding to the peak pull-out load of steel fiber increases with the increase in the embedment angle of steel fiber.The numerical simulation results in this section are in agreement with the conclusions reached by Long X [19].
Figure 11 shows the stress diagram of steel fiber pulled out at 30 • .Unlike the case in which the steel fibers are pulled out vertically from the concrete matrix, there is a certain contact pressure between the steel fiber and the left side of the matrix during the pullout process when the steel fibers have an embedment angle.With the increase in the embedment angle of steel fibers, the pressure of steel fibers on the left side of the matrix increases and the matrix is deformed continuously, which leads to a crack between the steel fiber and the right side of the matrix; debonding occurs first, and the pull-out channel of the steel fiber changes.Thus, the pull-out load-bearing capacity of steel fiber decreases.models with different embedment angles of steel fibers, such as 0°, 15°, 30°, 45°, and 60°, are established and simulated as shown in Figure 5b.The displacement loads at all embedment angles are perpendicular to the concrete matrix surface and the embedment depth of the steel fibers is 10 mm.
The results shown in Figure 10 were obtained with only the embedment angle of the steel fibers changed.The simulation results show the following: (1) the magnitude of the peak pull-out load decreases with the increase in the embedment angle of the steel fiber; (2) the displacement corresponding to the peak pull-out load of steel fiber increases with the increase in the embedment angle of steel fiber.The numerical simulation results in this section are in agreement with the conclusions reached by Long X [19].Figure 11 shows the stress diagram of steel fiber pulled out at 30°.Unlike the case in which the steel fibers are pulled out vertically from the concrete matrix, there is a certain contact pressure between the steel fiber and the left side of the matrix during the pull-out process when the steel fibers have an embedment angle.With the increase in the embedment angle of steel fibers, the pressure of steel fibers on the left side of the matrix increases and the matrix is deformed continuously, which leads to a crack between the steel fiber and the right side of the matrix; debonding occurs first, and the pull-out channel of the steel fiber changes.Thus, the pull-out load-bearing capacity of steel fiber decreases.

Effect of Interface Properties
According to previous research, the debonding damage of the fibers pulled from the matrix mainly depends on the strength of the matrix itself and the strength of the interface layer between the fibers and the matrix.There are two main forms of damage: (1) when the bond strength between the fiber and the matrix is less than the load of the fiber pulled out, the interface between fiber and matrix will be damaged, resulting in fiber falling out, the interface between the fiber and the matrix is damaged, resulting in the fiber falling off; (2) when the bond strength between the fiber and the matrix is greater than the load that the fiber is pulled out of, the concrete matrix itself is damaged, which causes the fiberbonded part of the concrete to fall off.In this section, the influence of interface properties on the mechanical properties of steel fiber pull-out is investigated using the finite element model with a steel fiber embedment depth of 6 mm, as shown in Figure 5a.

Failure of the Interface Layer between Steel Fiber and Concrete Matrix
When steel fiber reinforced concrete material is subjected to an external load, part of the stress of the concrete matrix is transferred to the steel fiber through the interface layer.Since the strength of steel fibers is much greater than the strength of the concrete matrix, whether the steel fibers play a role mainly depends on whether the interface layer fails.This section mainly studies the influence of the bond strength  and fracture energy C G of the cohesive element of the interface layer on the mechanical properties of steel fiber pull-out.

Effect of Interface Properties
According to previous research, the debonding damage of the fibers pulled from the matrix mainly depends on the strength of the matrix itself and the strength of the interface layer between the fibers and the matrix.There are two main forms of damage: (1) when the bond strength between the fiber and the matrix is less than the load of the fiber pulled out, the interface between fiber and matrix will be damaged, resulting in fiber falling out, the interface between the fiber and the matrix is damaged, resulting in the fiber falling off; (2) when the bond strength between the fiber and the matrix is greater than the load that the fiber is pulled out of, the concrete matrix itself is damaged, which causes the fiber-bonded part of the concrete to fall off.In this section, the influence of interface properties on the mechanical properties of steel fiber pull-out is investigated using the finite element model with a steel fiber embedment depth of 6 mm, as shown in Figure 5a.

Failure of the Interface Layer between Steel Fiber and Concrete Matrix
When steel fiber reinforced concrete material is subjected to an external load, part of the stress of the concrete matrix is transferred to the steel fiber through the interface layer.Since the strength of steel fibers is much greater than the strength of the concrete matrix, whether the steel fibers play a role mainly depends on whether the interface layer fails.This section mainly studies the influence of the bond strength τ and fracture energy G C of the cohesive element of the interface layer on the mechanical properties of steel fiber pull-out.
(1) Effect of bond strength τ When analyzing the influence of interface bond strength on the mechanical properties of steel fiber pull-out, G C = 0.25 N/mm is taken as a certain value and τ max is taken as 0.2 MPa, 0.5 MPa, 0.8 MPa, 1.1 MPa, and 1.4 MPa for five different groups of values, respectively, while all other parameters remain unchanged.The peak pull-out load and displacement corresponding to the peak pull-out load for steel fibers are obtained by numerical simulation, as shown in Figure 12.
whether the steel fibers play a role mainly depends on whether the interface layer fails.This section mainly studies the influence of the bond strength  and fracture energy C G of the cohesive element of the interface layer on the mechanical properties of steel fiber pull-out.
(1) Effect of bond strength  When analyzing the influence of interface bond strength on the mechanical properties of steel fiber pull-out, C G = 0.25 N/mm is taken as a certain value and max  is taken as 0.2 MPa, 0.5 MPa, 0.8 MPa, 1.1 MPa, and 1.4 MPa for five different groups of values, respectively, while all other parameters remain unchanged.The peak pull-out load and displacement corresponding to the peak pull-out load for steel fibers are obtained by numerical simulation, as shown in Figure 12.The results show that when τ max is 0.5~1.4MPa, the pull-out capacity of steel fiber increases by 66.6~244.5% compared with 0.2 MPa.It can be seen that the peak tensile force increases greatly with the increase in interface bonding strength, and the reinforcement effect of steel fiber is better.This is because τ max represents the stress value required to be reached by the interface element for the debonding between the steel fiber and the concrete matrix.Therefore, the higher the τ max , the greater the load required for the debonding of the steel fiber.It can be seen that the bonding strength of the interface layer has a great influence on the mechanical properties of steel fiber pulled out from the concrete matrix.The bearing capacity of the member or structure can be effectively improved by enhancing the bonding strength of the interface layer.
(2) Effect of fracture energy G C of the interface layer To analyze the influence of the fracture energy of the cohesive element on the mechanical properties of steel fiber pull-out, τ max = 0.5 MPa is taken as a certain value and G C is taken as 0.15, 0.25, 0.35, 0.45, and 0.55 N/mm for five different groups, respectively, while the other parameters remain unchanged.The peak pull-out load and displacement corresponding to the peak pull-out load for steel fibers are obtained by numerical simulation, as shown in Figure 13.The results show that when max  is 0.5~1.4MPa, the pull-out capacity of steel fiber increases by 66.6~244.5% compared with 0.2 MPa.It can be seen that the peak tensile force increases greatly with the increase in interface bonding strength, and the reinforcement effect of steel fiber is better.This is because max  represents the stress value required to be reached by the interface element for the debonding between the steel fiber and the concrete matrix.Therefore, the higher the max  , the greater the load required for the debonding of the steel fiber.It can be seen that the bonding strength of the interface layer has a great influence on the mechanical properties of steel fiber pulled out from the concrete matrix.The bearing capacity of the member or structure can be effectively improved by enhancing the bonding strength of the interface layer.The results show that when G C takes 0.25~0.55N/mm, the pull-out load-bearing capacity of steel fiber is increased by 0.24~0.46%compared with when G C takes 0.15 N/mm.Because the fracture energy G C represents the energy consumed by the debonding be-tween the steel fiber and the concrete matrix caused by the stiffness degradation of interface elements, the external load also increases with the increase in the load required" to "increases with the increase of G C , which increases the load required for the debonding between the steel fiber and the concrete matrix; that is, the bearing capacity of steel fiber reinforced concrete specimens is improved.However, according to the simulation results, the influence is much smaller than the influence of bonding strength on the pull-out performance of steel fibers, and the peak tensile strength increases gently.

Failure of the Concrete Matrix Itself
According to the scanning electron microscope image of the fiber damage form in ref. [25], it can be seen that the matrix's damage is the main damage form for fiber pull-out.Therefore, in this section, it is assumed that the interface bond strength between steel fibers and the concrete matrix is high (τ max = 3 MPa) and no debonding damage occurs, and the influence of the mechanical properties of the concrete matrix itself on the mechanical properties of steel fiber pull-out is analyzed. (

1) Effect of tensile strength f t
To analyze the influence of the tensile strength of concrete matrix material on the mechanical properties of steel fiber pull-out, G C = 0.25 N/mm is taken as a certain value and f t is taken as 1.9, 2.4, and 2.9 MPa for three different groups of values, respectively, while the other parameters remain unchanged.
When f t = 2.9 MPa and the steel fiber pull-out displacement is 0.1 mm, the damage distribution diagram of the model is shown in Figure 14.Damage is applied to CPS4T (concrete matrix) elements with values between 0 and 1, and when the value of damage is close to 1, it indicates that the element is damaged completely.It can be seen from Figure 14 that the damage in the localization zone of the crack evolves from 0 to 1, and the damaged part can basically be regarded as a real crack.From the damage diagram, it can be seen that the steel fibers carry part of the concrete to fall off in a conical shape at this time.Therefore, in terms of both appearance and specific values, the phase-field regularized cohesive zone model can accurately reflect the real physical cracks and is suitable for the study of steel fiber pull-out.To analyze the influence of the tensile strength of concrete matrix material on t mechanical properties of steel fiber pull-out, C G = 0.25 N/mm is taken as a certain val and t f is taken as 1.9, 2.4, and 2.9 MPa for three different groups of values, respective while the other parameters remain unchanged.When t f = 2.9 MPa and the steel fiber pull-out displacement is 0.1 mm, the dama distribution diagram of the model is shown in Figure 14.Damage is applied to CPS (concrete matrix) elements with values between 0 and 1, and when the value of dama is close to 1, it indicates that the element is damaged completely.It can be seen from Figu 14 that the damage in the localization zone of the crack evolves from 0 to 1, and the dam aged part can basically be regarded as a real crack.From the damage diagram, it can seen that the steel fibers carry part of the concrete to fall off in a conical shape at this tim Therefore, in terms of both appearance and specific values, the phase-field regulariz cohesive zone model can accurately reflect the real physical cracks and is suitable for t study of steel fiber pull-out.The peak pull-out load and displacement corresponding to the peak pull-out load f steel fibers are obtained by numerical simulation, as shown in Figure 15.The results sho that when the tensile strength t f is taken as 2.4 and 2.9 MPa, the pull-out load-bearing c pacity of steel fiber is increased by 20.84 and 40.43%, respectively, compared with wh t f takes 1.9 MPa.It can be seen that the peak pull-out load increases with the increase the tensile strength of the matrix, and the corresponding displacement of the peak pu out load also increases slightly.It is also suggested that the bearing capacity of steel fib reinforced concrete can be improved by improving the strength of concrete in engineerin practices.The peak pull-out load and displacement corresponding to the peak pull-out load for steel fibers are obtained by numerical simulation, as shown in Figure 15.The results show that when the tensile strength f t is taken as 2.4 and 2.9 MPa, the pull-out load-bearing capacity of steel fiber is increased by 20.84 and 40.43%, respectively, compared with when f t takes 1.9 MPa.It can be seen that the peak pull-out load increases with the increase in the tensile strength of the matrix, and the corresponding displacement of the peak pull-out load also increases slightly.It is also suggested that the bearing capacity of steel fiber reinforced concrete can be improved by improving the strength of concrete in engineering practices.t pacity of steel fiber is increased by 20.84 and 40.43%, respectively, compared with when t f takes 1.9 MPa.It can be seen that the peak pull-out load increases with the increase in the tensile strength of the matrix, and the corresponding displacement of the peak pullout load also increases slightly.It is also suggested that the bearing capacity of steel fiber reinforced concrete can be improved by improving the strength of concrete in engineering practices.(2) Effect of fracture energy G C of the concrete matrix To analyze the influence of the fracture energy of concrete matrix material on the mechanical properties of steel fiber pull-out, f t = 2.9 MPa is taken as a certain value and G C is taken as 0.15, 0.25, and 0.35 N/mm for three different groups of values, respectively, while the other parameters remain unchanged; the numerical simulation results are shown in Figure 16.The results show that when the fracture energy C G is taken as 0.25 and 0.35 N/mm, the pull-out load-bearing capacity of steel fiber is increased by 5.3 and 8.7%, respectively, compared with when C G takes 0.15 N/mm.It can be seen that the peak pull-out load increases with the increase in the fracture energy of the matrix material and the corresponding displacement of the peak pull-out load can also increase slightly.Therefore, the bearing capacity of steel fiber reinforced concrete can be improved by increasing the fracture energy of the concrete matrix material.

Conclusions
At present, steel fiber reinforced concrete is widely used in civil engineering because of its excellent mechanical properties, but there are still many problems to be studied in terms of the damage, fracture failure mechanism, and evolution process of steel fiber reinforced concrete.The pull-out test and numerical simulation of concrete specimens containing single steel fibers are important means for investigating the bond strength of the interface between the steel fibers and the concrete.Compared with related published numerical simulation work, this paper adopts the newly proposed unified phase-field theory, which has the following advantages in the characterization of concrete materials: it The results show that when the fracture energy G C is taken as 0.25 and 0.35 N/mm, the pull-out load-bearing capacity of steel fiber is increased by 5.3 and 8.7%, respectively, compared with when G C takes 0.15 N/mm.It can be seen that the peak pull-out load increases with the increase in the fracture energy of the matrix material and the corresponding displacement of the peak pull-out load can also increase slightly.Therefore, the bearing capacity of steel fiber reinforced concrete can be improved by increasing the fracture energy of the concrete matrix material.

Conclusions
At present, steel fiber reinforced concrete is widely used in civil engineering because of its excellent mechanical properties, but there are still many problems to be studied in terms of the damage, fracture failure mechanism, and evolution process of steel fiber reinforced concrete.The pull-out test and numerical simulation of concrete specimens containing single steel fibers are important means for investigating the bond strength of the interface between the steel fibers and the concrete.Compared with related published numerical simulation work, this paper adopts the newly proposed unified phase-field theory, which has the following advantages in the characterization of concrete materials: it can accurately describe the formation, propagation, and bifurcation of cracks; there is no stress oscillation, which makes it easy to implement numerically; and the grid has low sensitivity, which can ensure high calculation accuracy.This paper calculates and analyzes the influences of different factors and different failure forms on the steel fiber pull-out process, and can thus provide a more comprehensive and in-depth study on the bonding properties of the interface between steel fibers and concrete.The following main conclusions are drawn: (1) Based on the cohesive zone model and the unified phase-field theory, the pull-out loaddisplacement curve obtained by the simulation of the steel fiber pull-out test is in good agreement with the test results, which indicates that the finite element model of steel fiber pull-out established by the cohesive zone model and the phase-field regularized cohesive zone model has a certain accuracy and reliability.In addition, according to the results, the phase-field regularized cohesive zone model can represent the damage and failure evolution process of a concrete matrix intuitively and accurately.(2) By studying the influence of different factors on the pull-out mechanical properties of steel fiber, it is found that the pull-out load-bearing capacity of steel fiber when the embedment depth of steel fiber is 8mm and 10mm is increased by 30.7 and 61.37%, respectively, compared with that when the embedment depth is 6 mm.However, the pull-out load-bearing capacity of steel fiber is reduced with the increase in the length-diameter ratio and the embedment angle of the steel fiber.Therefore, these factors can be adjusted in actual projects to achieve the purpose of improving the mechanical properties of steel fiber reinforced concrete.(3) By studying the influence of interface properties on the pull-out mechanical properties of steel fiber, it is found that the bond strength of the interface layer and the strength of the concrete matrix have great influences on the pull-out mechanical properties of steel fiber.Among these influences, when the interface bond strength is 0.5~1.4MPa, the pull-out load-bearing capacity of steel fiber is increased by 66.6~244.5% compared with 0.2 MPa; and when the matrix strength is 2.4 and 2.9 MPa, the pull-out loadbearing capacity of steel fiber is increased by 20.84 and 40.43%, respectively, compared with that when the matrix strength is 1.9 MPa.Therefore, the reinforcing effect of steel fiber on the concrete matrix can be improved by controlling the bond strength of the interface layer and the strength of the concrete matrix.
Considering factors such as computational resources and costs, numerical simulations were conducted in this paper for only a 2D finite element model of steel fiber reinforced concrete.The 3D finite element model should be able to reflect the steel fiber pull-out test more realistically and comprehensively.Therefore, future efforts could establish the 3D and meso-level finite element models based on unified phase-field theory under all kinds of conditions.

Figure 2 .
Figure 2. Fracture diagram of the phase-field model: (a) sharp crack and (b) localized representation.

Figure 2 .
Figure 2. Fracture diagram of the phase-field model: (a) sharp crack and (b) localized representation.

3 .
Establishment and Verification of Numerical Calculation Model for Steel Fiber Pull-Out 3.1.Model Building In this section, numerical simulations are carried out on the experiments conducted by Shannag et al. [9] to investigate the pull-out process of steel fiber; the schematic diagram and dimensions of the steel fiber pull-out specimen are shown in Figure 4a.The size of the concrete matrix is 23 mm × 25 mm, and the fiber diameter is 0.19 mm.Sustainability 2023, 15, x FOR PEER REVIEW 10 of 21 (a) (b)

Figure 4 .
Figure 4.The pull-out specimen of Shannag: (a) schematic diagram of the pull-out specimen and (b) finite element model of the pull-out specimen.

Figure 4 .
Figure 4.The pull-out specimen of Shannag: (a) schematic diagram of the pull-out specimen and (b) finite element model of the pull-out specimen.

Figure 5 .
Figure 5.The pull-out model of steel fiber: (a) geometric conditions and boundary conditions of steel fiber pull-out model and (b) pull-out model of steel fiber with inclination.

Figure 5 .
Figure 5.The pull-out model of steel fiber: (a) geometric conditions and boundary conditions of steel fiber pull-out model and (b) pull-out model of steel fiber with inclination.

Figure 7 .
Figure 7. Pull-out load-displacement curve of steel fiber with different embedment depths.Figure 7. Pull-out load-displacement curve of steel fiber with different embedment depths.

Figure 7 . 21 Figure 8 .
Figure 7. Pull-out load-displacement curve of steel fiber with different embedment depths.Figure 7. Pull-out load-displacement curve of steel fiber with different embedment depths.Sustainability 2023, 15, x FOR PEER REVIEW 13 of 21

Figure 8 .
Figure 8. Peak pull-out load and displacement of steel fiber with different embedment depths.

Figure 9 .
Figure 9. Peak pull-out load and displacement of steel fibers with different length-diameter ratios.Figure 9. Peak pull-out load and displacement of steel fibers with different length-diameter ratios.

Figure 9 .
Figure 9. Peak pull-out load and displacement of steel fibers with different length-diameter ratios.Figure 9. Peak pull-out load and displacement of steel fibers with different length-diameter ratios.

Figure 10 .
Figure 10.Peak pull-out load and displacement of steel fiber with different embedment angles.

Figure 10 .Figure 11 .
Figure 10.Peak pull-out load and displacement of steel fiber with different embedment angles.Sustainability 2023, 15, x FOR PEER REVIEW 15 of 21

( 1 )
Effect of bond strength  When analyzing the influence of interface bond strength on the mechanical properties of steel fiber pull-out, C G = 0.25 N/mm is taken as a certain value and max  is taken as 0.2 MPa, 0.5 MPa, 0.8 MPa, 1.1 MPa, and 1.4 MPa for five different groups of values, respectively, while all other parameters remain unchanged.The peak pull-out load and

Figure 11 .
Figure 11.Pull-out stress diagram of steel fiber.

Figure 12 .
Figure 12.Variation of peak pull-out load and displacement with interface bond strength.Figure 12. Variation of peak pull-out load and displacement with interface bond strength.

Figure 12 .
Figure 12.Variation of peak pull-out load and displacement with interface bond strength.Figure 12. Variation of peak pull-out load and displacement with interface bond strength.

( 2 )
Effect of fracture energy C G of the interface layerTo analyze the influence of the fracture energy of the cohesive element on the mechanical properties of steel fiber pull-out, max  = 0.5 MPa is taken as a certain value and C G is taken as 0.15, 0.25, 0.35, 0.45, and 0.55 N/mm for five different groups, respectively, while the other parameters remain unchanged.The peak pull-out load and displacement corresponding to the peak pull-out load for steel fibers are obtained by numerical simulation, as shown in Figure13.

Figure 13 .
Figure 13.Variation of peak pull-out load and displacement with interface fracture energy.

Figure 13 .
Figure 13.Variation of peak pull-out load and displacement with interface fracture energy.

Figure 15 .
Figure 15.Variation of peak pull-out load and displacement with matrix strength.

( 2 )
Effect of fracture energy C G of the concrete matrix

Figure 15 .
Figure 15.Variation of peak pull-out load and displacement with matrix strength.

Sustainability 2023 ,
15, x FOR PEER REVIEW 18 of 21To analyze the influence of the fracture energy of concrete matrix material on the mechanical properties of steel fiber pull-out, t f = 2.9 MPa is taken as a certain value and C G is taken as 0.15, 0.25, and 0.35 N/mm for three different groups of values, respectively, while the other parameters remain unchanged; the numerical simulation results are shown in Figure16.

Figure 16 .
Figure 16.Variation of peak pull-out load and displacement with matrix fracture energy.

Figure 16 .
Figure 16.Variation of peak pull-out load and displacement with matrix fracture energy.