Similar articles

Abstract: The fracture processes of concrete were described by a cohesive crack model based on initial toughness criterion. The corresponding analytical method to predict the instability state was proposed. In this model, the initial toughness was adopted as the crack propagation criterion and the weight function method was used to calculate the stress intensity factor and the crack opening displacement caused by the cohesive stress. The unstable toughness can be easily obtained using the proposed method without measuring parameters at the critical state that was necessary in traditional methods. The proposed method was verified by existing experimental data of wedge splitting specimens with different grades of concrete and the sensitivity of the results on the tensile softening curve was discussed. The results demonstrate that the proposed method can well predict the peak load, the critical effective crack length, and the unstable toughness of concrete specimens. Moreover, the calculated unstable toughness is not sensitive to the tensile softening curve.

基于起裂韧度准则的混凝土失稳韧度预测研究

研究目的：研究混凝土失稳韧度的理论预测方法。
创新要点：1.提出混凝土失稳韧度的理论预测方法；2.利用楔入劈拉试件计算不同级配混凝土的失稳韧度；3.研究失稳韧度受拉伸软化曲线的影响。
研究方法：1.基于起裂韧度扩展准则，采用理论分析手段研究混凝土的失稳韧度计算方法；2.利用楔入劈拉试件（见图3）计算不同级配混凝土的失稳韧度。
重要结论：1.基于起裂韧度准则可合理地预测峰值荷载状态及失稳韧度；2.混凝土失稳韧度受断裂能的影响较小；3.拉伸软化曲线对混凝土失稳韧度的影响较小。

关键词：混凝土；起裂韧度；失稳韧度；楔入劈拉实验

1.  Introduction

 A large number of experimental studies on concrete fracture have shown that an obvious fracture process zone (FPZ) exists at the crack tip (Bazant and Planas, 1998). The traditional linear elastic fracture mechanics (LEFM), is no longer applicable to the analysis of the concrete fracture process (Jenq and Shah, 1985a) when the size of FPZ is comparable to the structure. To date, numerous nonlinear models for concrete fracture have been developed from different perspectives and based on various assumptions. The primary models and their characteristics are briefly summarized herein.

 The fictitious crack model (FCM) (Hillerborg et al., 1976) and crack band model (CBM) (Bazant and Oh, 1983) both considering the softening relation in the FPZ, were primarily solved numerically. The size effect law (SEL) developed by Bazant (1984) analyzes the size effect of the fracture by comparing nominal stress in the failure of structures with different sizes. The effective crack model (ECM) (Karihaloo and Nallathambi, 1990), known as the equivalent LEFM, has an explicit analytical solution. The two-parameter fracture model (TPFM) proposed by Jenq and Shah (1985a) was based on the assumption that crack propagation starts when the stress intensity factor and crack tip opening displacement (CTOD) reach their critical values. The critical stress intensity factor in both ECM and TPFM was calculated from the peak load P _max and critical effective crack length a _c using a formula developed in LEFM.

 Three different stages in the concrete fracture processes have been identified, namely, crack initiation, stable crack propagation, and unstable crack propagation (Xu and Reinhardt, 1999a). To depict the three stages of crack propagation without compromising simplicity, Xu and Reinhardt (1999a) developed the double-K model, which employed the initial toughness and the unstable toughness as control parameters of the two important instantaneous states of crack initiation and instability, respectively. Macroscopic damage in concrete is believed to be initiated once the stress intensity factor at the crack tip reaches (Zhang and Xu, 2011). Crack instability is the catastrophic point between the stable and unstable crack propagation stages and is significant in investigating the mechanisms of concrete fracture. The double-K model has been playing an important role in engineering practice. For instance, and have been used in safety warning systems (DL/T 5332-2005, 2006) and in evaluating the performance of concrete structures with cracks.

 The initial toughness can be easily obtained using formulae of the initial cracking load and initial crack length based on LEFM (Xu and Reinhardt, 1999a). Several measurement methods for the initial cracking load, such as tests using photo-elastic coating, laser speckle (Xu and Reinhardt, 1999a) and strain gauge (Zhang and Xu, 2011), and analytical methods have been formulated. Xu and Reinhardt (1999b; 1999c; 2000) proposed the double-K method and developed a simplified method later. Kumar and Barai (2009) developed a weight function method. The comparison of the above analytical methods can be found in (Zhang and Xu, 2011). Recently, Qing and Li (2013) proposed a theoretical method to obtain based on experimental peak load. Alternatively, as an easy method, Reinhardt and Xu (1999) and Zhang et al. (2010) determined the initial cracking load by the transition point that separates the linear and nonlinear segments of the load P-CMOD (crack mouth opening displacement) curve.

 Nevertheless, predicting the unstable toughness is difficult without making any assumption. In the experimental aspect, can be calculated by the peak load P _max and the critical effective crack length a _c, which were mainly measured at the peak load state through the three point bending (TPB) specimens or wedge splitting (WS) specimens (DL/T 5332-2005, 2006). However, the critical effective crack length a _c was usually difficult to measure accurately. In traditional methods, a _c was calculated by critical crack mouth opening displacement CMOD_c through empirical formulae (Xu and Reinhardt, 1999b; 1999c; Kumar and Barai, 2009; Zhang and Xu, 2011). Furthermore, in the numerical aspect, a similar step-by-step procedure must be followed until the maximum load was reached (Wu et al. 2007; Dong et al., 2013). Therefore, it is of significance to establish a simple theoretical method for predicting the unstable toughness .

 Based on the analysis of concrete fracture mechanisms, the current study attempts to develop an easy-to-use theoretical method for predicting the unstable toughness . In this method, the initial toughness was adopted as the control parameter of crack propagation and a weight function method was used to calculate the stress intensity factor and the crack opening displacement caused by the cohesive stress. Then the proposed method was verified using experimental data on wedge splitting specimens. The sensitivity of the results to the tensile softening curve was discussed.

2.  Model development

2.1.  Fracture processes of concrete

 As mentioned in the previous section, the fracture processes in concrete structures include three different stages: crack initiation, stable crack propagation, and unstable fracture. To identify the three different states, the cohesive crack model with crack tip singularity (Elices and Planas, 1991) was adopted in the present work. Fig. 1 shows the different characteristics of the three stages, including the initial crack length a ₀, crack propagation length Δa, external load P, stress strength factor K _I, and fracture toughness K _IC. Take type I load for example. When concrete is subjected to monotonic loading, cracks will not propagate until the external load P reaches the crack initiation load P _ini (P=P _ini), and (Fig. 1a). After crack initiation and as the external load P continuously increases, the crack starts to propagate, and a cohesive zone forms ahead of the initial crack, as shown in Fig. 1b. This period is commonly known as the stable crack propagation stage. Then the external load P reaches the peak (critical) value P _max (Fig. 1c). After that, unstable crack propagation may occur and the external load P decreases.




 A typical P-a/D (a is the effective crack length, and D is the specimen height) curve is shown in Fig. 2 (Reinhardt and Xu, 1999; Kumar and Barai, 2008; Qing and Li, 2013). When P reaches the initial cracking load P _ini, crack begins to grow and it gradually and nonlinearly increases with a. When P reaches the peak load P _max, a=a _c. Then P gradually decreases with a. The derivation of P to a at P=P _max can be assumed to be continuous (Qing and Li, 2013). The maximum theory can also be adopted to predict the crack instability.

2.2.  Criterion for concrete fracture

 As discussed above, a crack is initiated in concrete once the external load reaches the crack load P _ini or the stress intensity factor increases to the initial toughness . Thus, can be regarded as the toughness of the structural material to crack growth attributable to external forces. In the present study, the following criterion for crack initiation and propagation was employed (Xu and Reinhardt, 1999a; Wu et al., 2007; Dong et al., 2013): (1) , where K _I is the stress intensity factor at the tip of the effective crack tip in a mode I fracture.

 Before crack initiation, LEFM can be applied, and K _I can be calculated by (2a) , where K _I ^P is the stress intensity factor attributable to the external load in a mode I fracture.

 After the crack initiation, K _I can be expressed by the superposition scheme: (2b) , where K _I ^C is the stress intensity factor attributable to cohesive stress.

2.3.  Cohesive stress distribution

 It is assumed that the cohesive stress distribution on the FPZ can be expressed by Eq. (3) (Li et al., 2012). The cohesive stress at the tip of the effective crack equals the tensile strength. The cohesive stress at the tip of the initial crack and the crack tip opening displacement CTOD satisfy the tensile softening curve. (3) , where σ _s is the stress on the initial crack tip, f _t is the tensile strength, and m is the cohesion distribution index.

 According to Eq. (3), the distribution of cohesive stress changes with the index m. When m=1, the cohesive stress distribution is linear. Eq. (3) is similar to the cohesive force distribution adopted by Reinhardt (1985). The difference is that the cohesive stress at the initial crack tip is not zero in the present model. Based on the theoretical model, the corresponding analytical approach is developed in the following section.

3.  Analytical method for fracture

 Take the WS specimen for example, the configuration of the WS specimen is shown in Fig. 3 (DL/T 5332-2005, 2006).




  K _I ^P can be expressed by (Xu and Reinhardt, 1999c): (4) , where α=a/D, and k(α) is a geometric factor, which can be calculated by (DL/T 5332-2005, 2006): (5) . Eq. (5) is valid for 0.2≤α≤0.8 with 2% accuracy (Xu and Reinhardt, 1999c).

 To obtain the value of K _I ^C, a distribution of cohesive stress is assumed on the fictitious crack in an infinite strip (Jenq and Shah, 1985b). Previous studies have shown that the descending segment of the P-CMOD curve is affected by the cohesive stress distribution, while the peak load and critical crack mouth opening displacement are almost not affected (Li et al., 2012). Therefore, to simplify the analysis of the fracture processes, a linear cohesive stress assumption (m=1) was adopted for the propagation state of concrete cracks in the current study. The cohesive stress distribution is shown in Fig. 4, where σ _s(CTOD) is the function of CTOD.




 The weight function method proposed by Kumar and Barai (2009) was adopted in this study to calculate K_I ^C: (6) , where g(a) can be expressed using the following fourth-order weight function: (7) , where A ₁=σ _s(CTOD), A ₂=[f _t−σ(CTOD)]/(a−a ₀), s=1−a ₀/a, M ₁, M ₂, and M ₃ can be expressed by α (Kumar and Barai, 2009).

 The CTOD and σ_s(CTOD) values are required to calculate K _I ^C using Eqs. (6) and (7), where σ_s(CTOD) can be expressed by (Reinhardt et al., 1986): (8) , where c ₁, c ₂, and w ₀ are parameters.

 CTOD in Eq. (8) needs to be calculated by adding up two displacements caused by the external load P and cohesive stress. Using Paris’ displacement formula (Tada et al., 2000; Mai, 2002), CTOD can be expressed by (9) , where E is the elasticity modules, ξ is the integration variable, m(x, a) can be expressed using the following fourth-order weight function (Kumar and Barai, 2009): (10) .

 Thus, using Eqs. (7)–(9), g(a) can be derived, and K_I ^C can be obtained using Eq. (6). Then, substituting Eq. (4) for K _I ^P and Eq. (6) for K _I ^C into Eq. (2b) leads to: (11) , where

 Eq. (11) obviously shows that the external load P can be explicitly expressed as a function of the effective crack length a. The derivation of this equation with respect to a at the moment of the crack instability can be expressed by the following partial differential equation: (12) , and is calculated by (13) , where

  g′(a) can be expressed as (14) , where

  k′(α) can be expressed as (15) .

 If the initial toughness is given, the critical effective crack length a _c can be calculated using Eq. (12) by iterating method. Then, the critical crack tip opening displacement CTOD_c and the peak load P _max can be obtained by substituting a=a _c into Eq. (9) and Eq. (11), respectively. The unstable toughness can be calculated by substituting a=a _c and P=P _max into Eq. (4). Hence, an analytical approach for calculating was developed.

4.  Results and discussion

4.1.  Verification of the analytical method

 In this section, the proposed method for concrete fracture was verified using the data from Xu et al. (2006)’s experiments on WS specimens with different maximum aggregate sizes. The parameters of the specimens are shown in Table 1. In the experiments, the crack initiation load P _ini was measured by resistance strain gauges and the initial toughness was obtained by the LEFM formula.




 According to the experimental conditions, the parameters in Eq. (8) for the proposed method were taken as follows: c ₁=3, c ₂=6.93, and w ₀=160 μm (Reinhardt et al., 1986).

 The calculated and measured values of the peak load P _max are compared in Tables 2–4. It can be seen from the comparison that the values of the peak load P _max calculated by the proposed method are generally in good agreement with the experimental values for different maximum aggregate sizes of concrete specimens. It can be concluded that the method is effective for predicting the instability of concrete.










 Moreover, the values of unstable toughness are also shown in Tables 2–4. The calculated values of by the proposed method generally agree well with those by experimental method. Some of the calculated results are slightly different from those obtained by the experimental method, for example, the ratio of the predicted result to experimental result of WS15 is 1.21. According to Eq. (4), the discrepancy in the unstable toughness results from the critical relative effective crack length a _c/D. The primary reason is that, as mentioned above, an empirical formula was adopted in the experimental method to calculate the critical effective length a _c. The proposed method, differing from the experimental method, is capable of considering the effects of both external load and cohesive force when calculating a _c.

4.2.  Effect of tensile softening curve

 To study the sensitivity of the method to the shape of the tensile softening curve, two groups of parameters characterizing different tensile softening curves were used to predict fracture using the proposed theoretical method.

 The parameters of first group were c ₁=3, c ₂=6.93, and w ₀=160 μm, which were adopted in the last section. For comparison, the second set of the softening curve parameters was taken as as c ₁=1.5, c ₂=6.3, and w ₀=90 μm. The tensile softening curves represented by these two sets of parameters are shown in Fig. 5.




 With these different parameters, the values of P _max, a _c/D, and were calculated using the proposed method, and the results are shown in Figs. 6–8. As can be seen from Figs. 6–8, the calculated P _max, a _c/D, and using the first group of parameters are slightly different from those of the other one. Generally, this disparity in the calculated results using these two sets of parameters is considered to be insignificant in view of the great difference between the two tensile softening curves. For example, the fracture energy G _F (area under the tensile softening curve) corresponding to the first softening curve is about twice that for the second one.

5.  Conclusions

 In this study, the instability state of fracture in concrete was predicted by a theoretical model which adopted the initial toughness as the crack propagation criterion, and used the weight function method to calculate the stress intensity factor and the crack opening displacement caused by the cohesive stress. The applicability of the proposed method was verified by experimental data obtained on WS specimens, and the parameters at the peak load state, such as the peak load P _max, critical effective crack length a _c and unstable toughness , were calculated using the proposed method. The good agreement between the calculated results and the experimental results demonstrates that the proposed method can accurately predict the unstable toughness . In addition, the sensitivity of the results to the tensile softening curve was discussed. The results showed that the proposed method for fracture is not sensitive to the tensile softening curve, which verifies the reasonability of the proposed method. Future studies can be conducted from many aspects, such as the investigation of measurement methods for the effective crack length a _c.

---

* Project supported by the National Natural Science Foundation of China (Nos. 51309073 and 51309203), the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20131317120012), and the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering (No. sklhse-2014-C-02), China