Design of an explicit crack bridging constitutive model for engineered cementitious composites using polyvinyl alcohol (PVA) or polyethylene (PE) fiber

: This paper aims to establish an explicit crack bridging model that can link engineered cementitious composites (ECC) behavior from single fiber to single crack scale, which is great of designing ECC featuring pseudo tensile strain hardening and multiple cracks expanding by tailoring microstructure and materials selection. In this study, fiber bridging stress was divided into three parts including fiber bridging stress with no rupture, fiber debonding fracture stress, and fiber pullout fracture stress. Subsequently, the fundamental crack bridging model was emerged when fiber bridging stress with no rupture subtracted from the fiber rupture stress in debonding and pullout stage. Moreover, two-way pullout and Cook-Gordon effect were also considered to establish the complete model. It was found that the two-way pullout situation of polyvinyl alcohol (PVA) fiber has a significant influence on the crack opening width due to its slip hardening property, while the Cook-Gordon effect presents a faint crack width increment for PVA-ECC. However, the Cook-Gordon effect makes a significant contribution to the crack opening width of the composite produced by polyethylene (PE). This building intelligible model presents a better prediction for ECC through the comparison of experimental data or real average crack width in previous models, thus confirming the validity of this model.


INTRODUCTION
Fibers with characters of high elastic modulus, high elongation rate, and high ultimate tensile strength are usually used to overcome the crack sensibility and intrinsic brittleness of concrete (1,2).ECC composited by PE or PVA, cement, and fine aggregates exhibit an applausive tensile pseudo strain hardening property, meanwhile, multiple tight cracking with a mean crack width under 100 μm was usually accompanied in the tensile plane (3)(4)(5)(6).In addition, the fiber-reinforced concrete containing polyolefin and steel fibers also could exhibit great toughness and crack resistance if the fibers achieve a suitable distribution (7)(8)(9), while it is difficult to enhance the ductility obtained by the strain hardening and multiple cracks expanding action compared to PE-ECC and PVA-ECC.Hence, the superior tensile properties of ECC contributed to its wide application prospects in infrastructural construction.
On basis of fracture mechanics and micromechanical principles, ECC was designed by linking from single fiber to single crack scale to achieve the tensile pseudo strain hardening properties (10,11).The crack bridging model of composites expresses the relation of fiber bridging stress versus crack opening distance, which was of primary importance to implement the multiple crack expanding and tensile pseudo strain hardening characters by energy-based criterion and strength-based criterion (12).Hence, tailoring raw microstructure including fiber, matrix, and fiber/matrix by fiber bridging law can capture a superior tensile property of ECC (13,14).Moreover, modeling the micromechanical fiber bridging stress vs crack opening width relation also could evolve the macroscopic tensile properties of ECC (15,16).
To date, continuous efforts have been made to establish the crack bridging relation based on the fiber failure models and interface properties.Li et al. (17) first established the fiber pullout model (FPM) which considered the interfacial frictional stress as a constant during the fiber pullout process.Thereafter, Lin et al. (18) introduce the interfacial slip-hardening coefficient to advance the precision of FPM where interfacial frictional strength is enhanced by the increasing fiber pullout distance.However, fiber fracture was neglected in above two models when fiber bridging stress is greater than the fiber apparent tensile strength.Based on FPM, Maalej et al. (19) first considered fiber rupture in debonding process and proposed the fiber pullout rupture model (FPRM).Afterward, Kanda et al. (20) included a chemical bond stress to advance FPRM where the chemical bond stress and interfacial frictional stress control the fiber-matrix interface properties.Thereafter, Lin et al.
Design of an explicit crack bridging constitutive model for engineered cementitious composites using polyvinyl alcohol... • 3 (21) further introduced the interfacial slip-hardening coefficient to the above advanced FPRM.Noticeably, fiber bridging stress during pullout stage could exceed the fiber tensile strength due to the slip-hardening character of single fiber (22).However, the situation of fiber rupture in slipping stage was no considered in FPRM and two advanced FPRMs.
In recent years, fiber rupture in slipping stage was gradually introduced to advance the fiber bridging model.Based on the advanced FPRM (21),Yang et al. (23) considered fiber rupture in the debonding and slipping process, matrix micro-spalling, and Cook-Gordon effect to improve the precision accuracy of the fiber bridging model.However, no explicit expression was presented in the model.Yu et al. (24) created the implicit micromechanical modeling of crack-bridging relations of hybrid-fiber strain-hardening cementitious composites considering the interaction between different fibers based on Yang's model.In addition, based on the advanced FPRM, Huang et al. ( 14) developed an explicit crack bridging model that included fiber rupture in debonding and slipping stage (21), while the fiber two way pullout was explicated in this model.However, the crack opening width at the fiber complete debonding point was neglected in the fiber limiting embedment and bridging stress, while Cook-Gordon effect was also excluded.Cook-Gordon effect presents a premature fiber/matrix interface debonding normal to the fiber axis caused by a tensile stress located ahead of a blunt matrix crack implementing towards a fiber under remote tensile load, which makes fiber debonding take place ahead of the matrix crack, resulting in stretching of a free fiber segment and additional crack opening (23).Thus, excluding Cook-Gordon effect would reduce the crack width and mitigate the prediction accuracy of fiber bridging model.Based on the reviewed fiber bridging model, it is more reasonable to establish an advanced crack bridging model including the overall properties of the fibers, slip-hardening interfacial coefficient, fiber rupture in debonding and slipping stage, Cook-Gordon effect, and fiber two-way pullout.
The main objective of this paper was to establish an explicit fiber bridging model of ECC that could be easily adopted for all ECC researchers.On basis of Lin et al.'s FPRM model (21), the fiber rupture zone in this study in debonding and slipping stage was separately calculated.Then, the fiber bridging model was derived when complete fiber bridging stress with no fracture subtracted the bridging stress in the fiber debonding and slipping zone.In addition, Cook-Gordon effect and fiber two-way pullout were included to modify the crack bridging model, thus evolving the ultimate relation of ECC.Finally, the current relation was verified against experimental fiber bridging behavior in a series of ECCs.This analytical model was expected to favor ECC material design in terms of pseudo strain hardening and multiple crack opening behavior.

FIBER BRIDGING THEORY
To obtain pseudo strain hardening and multiple crack expanding behavior of ECC under tensile load, the fiber bridging theory was fundamentally framed based on the energy-based criterion and strengthbased criterion (20,25).The typical relation is shown in Figure 1.
(1) Strength-based criterion [1] Where and are the matrix tensile cracking stress and peak fiber bridging stress in curve, respectively.Strength-based criterion is a prerequisite for multiple cracking development, indicating that fiber bridging stress could be effectively transferred from initial crack to the next crack only when the maximum peak fiber bridging stress exceeds the matrix tensile cracking stress depended on the fracture toughness and flaw size of the matrix (26), which control the start-up of cracks.
(2) Energy-based criterion [2] Where and are the crack tip toughness and complementary energy, respectively.The specific formula of and are shown in Figure 1.Energy-based criterion was derived by the reference (27) using Jintegral analysis, which guides saturated multiple cracking development only if the crack tip toughness is lower than the complementary energy.Moreover, crack tip toughness is approximately equal to the formula , where is tensile elastic modulus.Moreover, matrix cracking can absorb energy in the cracking process, which can affect the stress field intensity factor (KI).With KI increases, the stress at the crack tip is large enough to reach the fracture strength of the material, and the crack will expand and lead to fracture.So, the Km used in this work is the matrix fracture toughness, which is the critical value of KI.The matrix fracture toughness could be obtained by single crack three-point bend test per ASTM E 1820 -05a (28).
Simultaneous requirement of Equation [1] and [2] is indispensable to achieve the characteristics of pseudo strain hardening and multiple cracking of ECC.Both criterions are well employed to guide the materials tailoring to achieve high ductility performances by strain-hardening strength indices and strain-hardening energy indices (29).For PVA-ECC, the and should exceed 1.5 and 3, respectively, while that of PE-ECC should greater than 1.2 and 3, respectively (13).

SINGLE FIBER PULLOUT BEHAVIOR
Based on the crack plane, the fiber pullout model was designed during debonding and slipping stage.Figure 2 shows the fiber inclined angle and embedment length scheme, where is the inclined angle are the embedment length of the long end and short end, respectively.
are the crack opening distance corresponding to the complete debonding point of the long end and short end, respectively.The fiber starts to slip when exceeds .The composite bridging stress can be computed by summing all single fiber stressing fibers crossing the crack plane (18).The relation between fiber bridging stress and crack opening distance was established, which mainly included fiber, matrix, and fiber-matrix properties, as shown in Equation [3].Meanwhile, the crack displacement at full-debonding point was defined in Equation [4] (18,21).[4] Where and presents the fiber bridging stress during debonding and slipping periods, respectively.
is the matrix/fiber frictional stress, is the chemical bond stress, is the modulus of elasticity of fiber, is the slip hardening coefficient, presents fiber diameter.is the fiber embedment length, which is less than Moreover, is the snubbing coefficient (30).
, where is the fiber volume fraction and is the matrix tensile elastic modulus.
Additionally, the fiber apparent tensile strength exhibits a degradation effect, which decreases with the increasing inclined angle (31).For PVA fiber, the strength reduction effect could be accounted by Equation [5].[5] Where is fiber tensile fracture strength, is the fiber tensile strength reduction coefficient.

Averaging fiber bridging stress-crack opening width
Fiber location and orientation are the main factors to influence the relation.According to (7,32), the model could be accounted by summing single stress from individual bridging fibers considering probability density functions, as shown in Equation [6].[6] Where is the probability density function of centroid distance and orientation angle , is the fiber length.For a 2D or 3D uniform random distribution, corresponding formulas are defined as follows: [7] [8] [9] Where is the fiber volume fraction.It is noteworthy that is used with limits since excessive fi-  [3] bers would lead to the agglomeration phenomenon that mitigates the tensile properties and workability.is the function of embedment length and orientation angle .After changing with and replacing the Equation [7][8][9] into (6), Equation [6] could be changed as follows: [10]

No Fiber fractures
Firstly, assuming that all fiber in crack plane would not fracture under axial tensile force, the embedment length at full debonding point can be derived when crack opening width reaches (21), as shown in Equation [11].[11] Figure 3 illustrates the fiber embedment at full debonding point with varying crack opening width.Fibers between to are in slipping stage, in which fibers between to are in debonding behavior.Moreover, the crack opening width can be obtained when fibers having a peak embedment length of are in complete debonding point, which is given in Equation [12].[12] Hence, Fiber states without fiber rupture could be divided into two kinds, as plotted in Figure 4.The crack opening width could be distinguished into two types according to the value of , , and .Then, the crack bridging stress with no fibers rupture for ECC can be calculated by integrating all contributory fibers bridging stress, which is given in Equation [13]. [13]

Fiber rupture in debonding stage
Fibers will be fractured when the maximum fiber bridging stress is higher than its apparent tensile strength during debonding process.Specially, the fiber critical embedment length  is derived by making the fiber bridging stress in complete debonding point equal the fiber tensile fracture strength (21), as given in Equation [14]. [14] Figure 5 illustrates the fiber critical embedment and debonding rupture zone with varying crack opening width.The fiber embedment increases as the crack opening width increases, which could reach when crack opening width is up to .Moreover, the debonding rupture zone increases with the crack opening width increases up to .As the crack opening width increases, three scenarios could be divided to calculate the fiber bridging stress at debonding rupture zone, as shown in Figure 6.From Figure 6a, when lies between 0 and , three main zones could be classified for all fibers, which include slipping zone (SZ), debonding zone (DZ), and debonding rupture zone (DRZ).Fibers in debonding zone (DZ1) will be fractured under tensile load, while fibers in debonding zone (DZ2) suffer full debonding and slipping process without rupture.Moreover, fibers in debonding rupture zone have fractured, which have reached the embedment length at full debonding point.
produced by the intersection point between and can be calculated by Equation [15].Moreover, could be derived by making be equal to (21), as presented in Equation [16] and [17]. [15] [16] [17] When lies between and , the debonding rupture zone (DRZ) and slipping zone (SZ) were presented in Figure 6b.Fibers in DRZ have fractured and reached the embedment length at full debonding point.Additionally, Figure 6c shows that the debonding rupture zone of fibers reduces as increases from to .The inclined angle could be calculated by making equal (21), as defined by Equation [18].Hence, the fiber rupture bridging stress in bebonding stage versus the increasing crack opening width relation can be evolved by combining the contributory fractured fibers, as given in Equation [19].

Fiber fracture in slipping stage
Fiber bridging stress can be strengthened since the slip hardening feature of PVA fiber emerges in slipping stage.Hence, fiber slipping fracture also can be presented when the maximum fiber bridging stress grow higher than the fiber apparent  tensile strength .The maximum fiber bridging stress could be deduced by Equation [20], corresponding to (33).Furthermore, another fiber critical embedment length could be derived by making equal , as given in Equation [21].[18] [19] [20] [21] As crack opening width expands, a new critical embedment length resulting in fracturing for the slipping fibers can also be derived by letting equal as shown in Equation [22], indicating that fibers with an embedment length higher than has fractured in the current crack opening width (21).[22] Fiber slipping rupture zone is determined by the combined curve of , , and .However, is an inapprehensible curve since is the quadratic function of .Huang et al. (14) strictly omitted in Equation [22].Actually, an essential should be derived to achieve the better prediction accuracy of σ(δ) relation.Hence, a new assuming in this study, embedment length all fibers conservatively lie between and and presents a random uniform distribution in the slipping rupture zone, which was proposed to identify, and account for the slipping fracture section.Afterward, the probability density function of is proposed to get the average , as defined in Equation [23] and [24].The revised cur-rent critical embedment length can be shown as in Equation [25].[23] [24] [25] Figure 7 illustrates the limiting embedment length at slipping state and slipping rupture zone with varying crack opening distance.The slipping rupture zone gradually increases as crack opening expands.Figure 8 shows the fiber slipping rupture zone.It is noteworthy to see that the slipping rupture zone starts to occur when crack plane expands a minimum width that can be approximately calculated by making equal using Equation [26].Thus, the slipping rupture zone is zero when crack opening width is lower than , shown in Figure 8a.[26] Moreover, the slipping rupture zone increases as the crack plane increases from to (see Figure 8 b).The inclined angle could be calculated by letting equal (21), as given in Equation [27] and [28].The curve of and almost coincide when the crack opening width expand to (see Figure 8 c).Generally, the crack opening width corresponding to the peak fiber bridging stress is lower than (around 750 μm).Hence, the fiber rupture bridging stress in slipping stage versus increasing crack opening width (lower than ) for composites can be deduced by integrating the contributory fractured fibers, as presented in Equation [29].

Effective fiber bridging stress
On basis of above derivation, the fiber bridging stress without fiber fracture deducts the fiber debonding and slipping rupture stress, thus derivating the effective fiber bridging stress of composites, as presented in Equation [30]. [30]

FIBER TWO-WAY PULLOUT CONSIDER-ATION
Fiber two way pullout was first reported with respect to PVA fiber due to slip-hardening interface behavior (34).The exclude the fiber two-way pullout situation.In this study, the so-called nominal fiber was adopted to supply a reasonably approximate crack opening width contributed by fiber two-way pullout (14,23).
The single fiber bridging stress is the function of embedment length and inclined angle , thus Equation ( 10) can be rewritten as Equation [31] (23) Where is the efficiency of fiber bridging, which represents the fiber contributing stress as increasingly fracture across a crack plane.Thus, will gradually reduce to zero.It is noteworthy that the efficiency fiber contributing stress for fibers containing the rupture fiber in debonding and slipping stage is defined as which is given in Appendix I. Hence, according to Equation [31], is presented in Equation [33].For fibers having a long embedment length , under the effective bridging stress , the displacement of fiber in debonding and pullout stage can be deduced (14), as given in Equation [34] and [35], respectively.
[31]  Where is the nominal fiber bridging stress, can be derived in Equation [32].[32] [33] [33] [34] [35] Where . The displacement of is the function of maximum fiber bridging stress and embedment length of the long side in the loading process, as shown in Equation [36].[36] Assuming a uniform random distribution in the interval from to thus the probability density function of is defined as by Equation [37].Afterwards, the remedial crack opening width is approximately averaging by the integral of , presented in Equation [38].
[37] Moreover, it is known from Equation [34], the maximum fiber bridging stress in debonding stage has a critical value (see Equation [39]), which deduces the current fiber critical embedment length , as given in Equation [40]. [39] [40] Where .A long embedment length of fibers with higher than is still at the debonding stage.
According to the above definition, two scenarios could be generated according to namely , implying that the long embedded side of fibers is in debonding process, is given by Huang et al. (14).
[41] , namely , which indicates that the long embedded side of fibers is in the debonding and pullout process, is given as follows [42]: [42]

COOK-GORDON EFFECT
Cook-Gordon effect will contribute the additional width to the , describes a precocious fiber debonding normal process.Meanwhile, the fiber-matrix is separated under the fiber bridging stress related to the elastic crack tip field of the adjoining matrix crack in the horizontal plane, thus leading to an attached crack opening generated from the fiber elastic stretching of the fiber segment (23).Therefore, the fiber debonding process initiates ahead of the matrix cracking, leading to stretch a free fiber segment , as given in Equation [43] (23). [43] Where is the Cook-Gordon parameter.In general, for PVA fiber was suggested, while α is set as 15 for polyethylene (PE), and polypropylene (PP) fibers (35).
Overall, the total crack opening width including the contribution of fiber two-way pullout and Cook-Gordon effect is given as [44]: [44]

MODEL VALIDATION
Figure 9 shows the computation flow scheme of final .The experimental data including fiber, matrix, and fiber/matrix parameters from Yang et al. (23) were employed to verify the final .The micro-parameters for PE-ECC and PVA-ECC are listed in Table 1.It is well known that PVA fiber has a strong chemical bond strength and slip hardening behavior in hardened pastes.Table 1 illustrates the micro-parameters for PVA-ECC and PE-ECC.
Figure 10 shows theoretical relation of PVA-ECC analyzed by the current model.When only one way was considered, the predicted peak bridging stress and corresponding crack opening width were 64 μm and 1.50 MPa, respectively.If two ways are accounted, the same peak stress is obtained at the crack opening width of 165.8 μm.Moreover, the Cool-Gordon effect can add 2.1 μm to the crack opening width and reach 167.9 μm, indicating that the contribution of Cool-Gordon effect to crack opening width is faint due to the small assuming Cook-Gordon parameter for PVA fiber.Figure 11 illustrates the predicted composite relation of PVA-ECC.The predicted and experimental composite characteristics of PVA-ECCs and PE-ECCs are listed in Table 2. From Figure 11, it was difficult to research the variation of pseudo elastic properties of three models in the elastic part.The approximate first crack strength from three models was presented, which was lower than the experimental data.Afterward, cracks gradually develop during strain in three models, which show a different maximum fiber bridging and corresponding crack opening width.Moreover, the accuracy coefficients of the two key parameters by three models are listed in Table 2.For PVA-ECC with a volume fraction of 0.1%, the experimental curves scattered inside the grey zone referenced by (23), where its peak stress and corresponding crack opening width was 0.295 MPa and 165.0 μm, respectively, while that of Huang al. (14) presented a peak stress of 0.31 MPa at the crack width of 131.0 μm, and that of Yang et al. (23) has a peak stress of 0.35 Mpa at the crack width of 107.0 μm.Hence, employing the current model could further improve the accuracy of crack opening width that reaches 175.1μm with a proximal peak bridging stress of 3.0 MPa.In addition, for PVA-ECC with a volume fraction of 0.5 %, the maximum fiber bridging stress obtained by Huang's model, Yang's model, and the current model were 1.52 MPa, 1.72 MPa, and 1.50 MPa, respectively, while the corresponding crack opening width were 141.5μm, 125.6 μm, and 171.9 μm, respectively.Hence, compared to the real experimental maximum fiber bridging stress of 1.65MPa and crack opening width of 160.1 μm, the current model has a higher prediction accuracy than the other two models (see Table 2).It was clearly found that a proximate crack opening presented higher precision by using the current model than the other two predicted models.In fact, the current models use the opposite logic calculation compared to the Huang's model, in which fiber bridging stress with no rupture subtracts the fiber fracture stress in debonding and pullout stage.In Huang's model, the operative partial fiber bridging stress in the debonding and pullout stage are severalty picked, thereafter summing those partial stresses to form the fiber bridging stress of composites.Therefore, the current model provides an effective accuracy compared to the experimental data, which is also easy to understand for researchers focusing on the fiber reinforced cement composites field.Based on the micro-parameters in the previous study (36), Figure 12 demonstrates the comparison of the predicted σ − δ relation for PE-ECC by three models.It can be seen from Figure 12 and Table 2 that, the maximum bridging stress and corresponding crack opening width of the M-C group by using Yang's model are around 12.99 MPa and 99.0 μm, respectively, while that of Huang's model presented a maximum stress of 13.1 Mpa at the crack width of 56.8 μm.Noticeably, the Cool-Gordon effect was significant in verifying the real width of PE-ECC due to the big Cool-Gordon effect of α = 15 which was not included in Huang's model from Figure 12 (b), (c) and (d).Employing the current model could further improve this accuracy of crack opening width that reaches 97.9 μm with a slightly bigger peak bridg-ing stress of 14.0 MPa.Hence, Yang's model and the current model have a higher average crack width than Huang's model for PE-ECC since the Cool-Gordon effect of Huang's model could not effectively be implemented.On the other hand, Yang's model and current model show a proximal value for M-25, M-50, and M-100 compared to the real average crack width calculated from the crack patterns of the tensile test.A near-peak bridging stress was illustrated for three models.However, the current model is easy to calculate compared to the implicit Yang's model.However, the corresponding crack opening width at maximum bridging stress by three models is lower than the experimental data for PE-ECC since there is amount of inelastic strain produced by micro-cracks near the main crack (15).
Fiber length (m) Fiber bridging stress during debonding periods (Pa) V f Fiber volume fraction Fiber bridging stress during slipping periods (Pa) E f Fiber modulus of elasticity (Pa) Current debonding ruptured length (m) E m Elastic modulus of the matrix (Pa) Crack opening width when fibers are fully debonding (m) β Slip hardening coefficient L d (θ) Critical debonding embedment length (m) Fiber tensile strength reduction coefficient L p (θ) Critical slipping embedment length (m) f Snubbing coefficient Current slipping ruptured length (m) Fiber tensile fracture strength (Pa) Inclined angle when intersects with L d (θ) Matrix tensile cracking stress (Pa) Inclined angle when L f /2 intersects with L d (θ) Le Fiber embedment length (m) Nominal fiber bridging stress (Pa) L L Fiber embedment length of long end (m) Efficiency of fibers bridging Ls Fiber embedment length of short end (m) Cook-Gordon parameter τ 0 Matrix/fiber frictional stress (Pa) Fiber stretch segment according to Cook-Gordon effect G d Chemical bond stress (Pa) Matrix tensile cracking stress in

Figure 2 .
Figure 2. Fiber inclined angle and embedment length schematic diagram.

Figure 3 .
Figure 3.The fiber embedment at full debonding point with a varying crack opening width.

Figure 5 .
Figure 5.The fiber critical embedment and debonding rupture zone with a varying crack opening width.

Figure 7 .
Figure 7.The limiting embedment length in slipping stage and slipping rupture zone with varying crack opening distance.

Figure 12 .
Figure 12.Comparison of predicted σ − δ relation for PE-ECC by three models.

table 1 .
The micro-parameters for PVA-ECC and PE-ECC.

table 2 .
Predicted and experimental composite characteristics of PVA-ECCs and PE-ECCs.
a : 0.1% and 0.5% represent the volume fraction of PVA fibers.b : σ peak_mean is the mean peak stress of σ(δ) curve from experimental data.