# Crack Extension Resistance of Normal-Strength Concrete Subjected to Elevated Temperatures

Authors: Jing Chen; Zhoudao Lu
Journal: Advances in Materials Science and Engineering
 (2014)

---

## Abstract

Determination of the residual crack extension resistance curves (K
R-curves) associated with cohesive force distribution on fictitious crack zone of complete fracture process is implemented in present research. The cohesive force distributes according to bilinear softening traction-separation law proposed by Petersson. Totally ten temperatures varying from 20°C to 600°C and the specimen size of 230
×
200
×
200 mm with initial-notch depth ratios 0.4 are considered. The load-crack mouth opening displacement curves (P-CMOD) of postfire specimens are obtained by wedge-splitting method from which the stress intensity factor curves (K-curves) are calculated. In each temperature, with the distribution of cohesive force along the fracture process zone, the residual fracture toughness K
R (Δ
a) increases with increasing crack length Δ
a, whereas the K
R-curves decrease with increasing temperatures T
m for the thermal damage induced. The stability analysis on crack propagation demonstrates that when the residual K
R-curve is higher than K-curve, the crack propagates steadily; otherwise, the crack propagates unsteadily.


---

## Body

## 1. Introduction

To describe the crack propagation in concrete structures, the crack extension resistance in terms ofK
R has been measured and the related characteristics were investigated by Hilsdorf and Brameshuber in 1984 [1], Mai in 1984 [2], and Karihaloo and Shah in 1987 [3] as well as Xu and Reinhardt in 1999 [4] which was calculated using the conventional approach proposed by Irwin et al. in the early 1960s. Xu and Reinhardt in 1998 [5] proposed an analytical method to determine the crack extension resistance (K
R-curve) according to the cohesive force on the fictitious crack zone that is described directly by the softening traction-separation law. The basic principle of the approach is that the crack extension resistance is composed of two parts. One part is the inherent toughness K
Ic
ini, which resists the initial propagation of an initial crack under loading. This means that a crack does not propagate when the stress intensity factor at the initial crack tip is less than the inherent toughness K
Ic
ini. Another part is induced by the cohesive force distributed on the fictitious crack during crack propagation. Therefore, it is a function of the cohesive force distribution f
(
σ
), the tensile strength f
t of the material, and the length of the propagating crack.The main requirement for determining crack extension resistance curve based on cohesive force distribution during crack propagation is to know the load-crack mouth opening displacement (P-CMOD) curve a priori. Some characteristics of K
R-curve were investigated numerically on standard TPBT specimen for different concrete strength and specimen sizes [6]. It was observed that K
R-curve was dependent on compressive strength of concrete and had almost the same S-shape. It increased with increasing crack length and increasing strength of concrete. The obtained K
R-curves were almost independent of the specimen sizes. However, it was found that some difference could be noticed on the gained K
R-curve by using bilinear [7] and nonlinear [8] softening functions of concrete. The influence of specimen geometry (TPB and CT specimens) on the K
R-curves is considered by Kumar and Barai [9], and the influence of the specimen geometry on the K
R-curves was not observed for the specified specimen size and initial-crack length/depth ratio.The influence of temperature on the fracture properties was considered by several researchers, mainly on the fracture energy and material brittleness [10–14], relative fewer discussions on the fracture toughness [15, 16], not mention the crack extension resistance of the complete fracture process. For the stability analyses of crack propagation, the comparison between K
R-curves and the corresponding stress intensity factor curves (K-curves) can be taken as a crack propagation criterion to judge the stability of the crack in the loaded structures. This criterion could also be used in the analysis of the stability of structures that suffered high temperature or fire in a real situation.The main objective of this research is to determine the residual crack extension resistance (K
R-curve) of postfire concrete based on the cohesive force distributed on the fictitious crack zone, and the influence of temperatures on the K
R-curve is discussed. The wedge-splitting experiments of a total of ten temperatures varying from 20°C to 600°C and the specimens size 230 × 200 × 200 mm with initial-notch depth ratios of 0.4 are implemented.


## 2. Determination of ResidualK
R-Curve Based on Cohesive Stress Distribution

### 2.1. Background

According toK
R-curve criterion [5], the crack extension resistance of a cracked solid consists of the inherent toughness K
Ic
ini and the cohesive toughness K
c
(
Δ
a
) which increases with the increasing amount of crack extension. The cohesive toughness depends upon cohesive stress distribution f
(
σ
) which is a function of crack opening displacement w and tensile strength of concrete f
t and the propagating crack length a. At the onset of unstable crack propagation, the stress intensity factor K at the tip of the propagating crack is expressed as
(1)
K
=
K
R
(
Δ
a
)
,
where K
R
(
Δ
a
) is the crack extension resistance at crack extension length Δ
a
=
a
-
a
0. Also K
R
(
Δ
a
) is expressed in the following relation:
(2)
K
R
(
Δ
a
)
=
K
Ic
ini
+
K
c
(
Δ
a
)
,
(3)
K
c
(
Δ
a
)
=
F
1
(
f
t
,
f
(
σ
)
,
Δ
a
)
.In order to develop theK
R-curve for complete fracture process considering the cohesive stress in fictitious fracture zone, the value of cohesive toughness K
c at every stage of loading is important to determine. During crack propagation, four different stages are considered with the help of three characteristic crack lengths (a
0, a
c, and a
w
c) as represented in Figure 1 in which r
u is the undamaged portion of the ligament, a
0 is initial crack length, a
c is crack length at critical condition of unstable crack propagation, and a
w
c is the length of fully developed fictitious fracture zone after which the stress-free crack propagation will begin.Four stages of crack propagation during the fracture process.
(a)
Initial notch length
(b)
Critical notch extension
(c)
Fully developed crack zone
(d)
Formation of new stress-free crack length


### 2.2. Softening Traction-Separation Law of Postfire Concrete

The softening traction-separation law is a priori to determine theK
R-curve, at room temperature, and many expressions have been proposed based on direct tensile test [8, 17–20]. Based on numerical studies, simplified bilinear expressions for the softening traction-separation law (Figure 2) were suggested by Petersson in 1981 [17], Hilsdorf and Brameshuber in 1991 [19], and Phillips and Zhang in 1993. The area under the softening curve was defined as the fracture energy G
F by Hillerborg et al. in 1976 [21]. Therefore, one can get the following equation:
(4)
G
F
=
1
2
(
f
t
w
s
+
σ
s
w
0
)
.Figure 2
The bilinear softening traction-separation law.As a consequence, a general form of the simplified bilinear expression of the softening traction-separation law is given as follows:(5)
σ
=
{
f
t
-
(
f
t
-
σ
s
)
w
w
0
0
≤
w
≤
w
s
,
σ
s
(
w
0
-
w
)
(
w
0
-
w
s
)
w
s
≤
w
≤
w
0
,
0
w
≥
w
0
.Different values of the break point (σ
s, w
s) and the crack width w
0 at stress-free point were used for the expression proposed by different researchers. In present work, the bilinear softening function of concrete proposed by Petersson is used for postfire specimens as follows:
(6)
σ
s
=
f
t
3
,
w
s
=
0.8
G
F
f
t
,
w
0
=
3.6
G
F
f
t
.


### 2.3. Analytical Method

The standard Green function [22] for the edge cracks with finite width of plate subjected to a pair of normal forces is used to evaluate the value of cohesive toughness. The general expression for the crack extension resistance for complete fracture associated with cohesive stress distribution in the fictitious fracture zone for Mode I fracture is given as follows:
(7)
K
R
(
Δ
a
)
=
K
Ic
ini
+
∫
a
0
a
2
σ
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
,
where
(8)
F
(
x
a
,
a
h
)
=
3.52
(
1
-
x
/
a
)
(
1
-
a
/
h
)
3
/
2
-
4.35
-
5.28
x
/
a
(
1
-
a
/
h
)
+
{
1.30
-
0.30
(
x
/
a
)
3
/
2
1
-
(
x
/
a
)
2
+
0.83
-
1.76
x
a
}
×
{
1
-
(
1
-
x
a
)
a
h
}
.
Equation (7) is employed according to the conditions of four stages of crack propagation as mentioned below.


#### 2.3.1. Case 1: Whena
=
a
0

There is no advancement in the initial notch length at this stage of loading and the body remains in elastic condition, subjected to small load (up toP
ini) without any slow crack growth. Hence, cohesive stress σ
(
x
) = 0; the crack growth resistance remains equal to initiation toughness of the material. From (7), it is expressed as
(9)
K
R
=
K
Ic
ini
.


#### 2.3.2. Case 2: Whena
0
≤
a
≤
a
c

The stable slow crack growth will take place until the effective crack extensiona
c corresponding to the maximum load P
u is achieved. The cohesive force will start acting across the fictitious fracture zone resulting in the increase of crack extension resistance.(a) For specimens subjected to temperatures less than 120°C, the critical CTODc corresponding to maximum load P
u is less than w
s, as shown in Figure 3(a). The distribution of cohesive stress along the fictitious fracture zone is approximated to be linear, as shown in Figure 4(a). The variation of cohesive stress along the fictitious fracture zone for this loading condition, that is, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(10)
σ
(
x
)
=
σ
(
w
)
+
(
f
t
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
-
a
0
)
,
where σ
(
w
) and w are the values of cohesive stress and crack opening displacement, being at the tip of initialnotch, respectively. The value of σ
(
w
) is determined by using bilinear softening function as follows:
(11)
σ
(
w
)
=
σ
s
(
w
s
)
+
w
s
-
w
w
s
(
f
t
-
σ
s
(
w
s
)
)
.
The crack extension resistance in this case is evaluated using (7) and (8).Two different situations for CTODc and w
s.
(a)
When CTODc
≤
w
s
(b)
When CTODc
>
w
sThe cohesive distribution of crack propagation stage II.
(a)
The linear distribution of cohesive force
(b)
The bilinear distribution of cohesive force(b) For specimens subjected to temperatures higher than 120°C, the critical CTODc corresponding to maximum load P
u is wider than w
s, as shown in Figure 3(b). The distribution of cohesive stress along the fictitious fracture zone is approximated to be bilinear, as shown in Figure 4(b). The variation of cohesive stress along the fictitious fracture zone for this loading condition, also, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(12)
σ
1
(
x
)
=
σ
(
w
)
+
(
σ
s
(
w
s
)
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
s
-
a
0
)
a
s
≤
x
≤
a
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.
The value of σ
(
w
) is determined by using bilinear softening function as follows:
(13)
σ
(
w
)
=
w
0
-
w
w
0
-
w
s
σ
s
(
w
s
)
.
The limits of integration of (7) should be taken in two steps: a
0
≤
x
≤
a
s for cohesive stress σ
1
(
x
) and a
s
≤
x
≤
a for cohesive stress σ
2
(
x
), respectively. The same Green’s function F
(
x
/
a
,
a
/
h
) for a given effective crack extension will be determined using (8). The calculated formula is listed as follows:
(14)
K
R
=
K
Ic
ini
+
∫
a
0
a
s
2
σ
2
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
+
∫
a
s
a
2
σ
3
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
.
The effective crack length at break point a
s (as shown in Figure 4(b)) is computed from the following nonlinear expression [23] by substituting COD
(
a
s
), CMOD, a, and h:
(15)
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
,
where COD
(
a
s
) is the crack opening displacement at a
s, a is the effective crack length (according to (18)), and h is the specimen height.


#### 2.3.3. Case 3: Whena
c
≤
a
≤
a
w
0

During this stage of the applied load for all temperatures, corresponding CTOD and effective crack length have increased more than maximum loadP
u, CTODc, and critical effective crack extension a
c, respectively. The cohesive stress distribution for this case will take a bilinear shape, as shown in Figure 4(b). It is the same situation as case 2 part (b), and the crack extension resistance K
R at this stage would be calculated according to (14). This case is shown in Figure 5.Figure 5
Crack propagation stage III.


#### 2.3.4. Case 4: Whena
≥
a
w
0

This situation of loading corresponds to the descending portion ofP-CMOD curve. At the effective crack extension, a
w
0, full shape of cohesive stress distribution is allowed to develop, and beyond further extension in crack during loading, a new stress-free crack in front of the initial notch tip will form. This case is shown in Figure 6 and the stress distribution is expressed with the following relation:
(16)
σ
1
(
x
)
=
0
a
0
≤
x
≤
a
w
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
(
x
-
a
w
0
)
(
a
s
-
a
w
0
)
a
s
≤
x
≤
a
w
0
,
σ
3
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.Figure 6
Crack propagation stage IV.Similar to (15), the effective crack length a
w
0 corresponding to zero stress of new fictitious fracture zone as shown in Figure 6 is computed from the following nonlinear expression:(17)
COD
(
a
w
0
)
=
CMOD
{
(
1
-
a
w
0
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
w
0
a
-
(
a
w
0
a
)
2
]
}
1
/
2
,
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
.
The evaluation of crack extension resistance in this case is also done using (7) and (8).


### 2.4. Determination of Equivalent Crack Extension

#### 2.4.1. Assumption

Linear asymptotic superposition assumption [5, 6] is considered to introduce the nonlinearity effect inP-CMOD curves during loading of concrete test specimens. This assumption enables us to introduce linear elastic fracture mechanics (LEFM) for calculating fracture parameters during every stage of loading in concrete structures. The hypotheses of the assumption are given as follows.(1)
The nonlinear characteristic of theP-CMOD curve is caused by fictitious crack extension in front of a stress-free crack.
(2)
An effective crack consists of an equivalent-elastic stress-free crack and equivalent-elastic fictitious crack extension.


#### 2.4.2. Equivalent Crack Extension for WS Specimen

The equivalent-elastic crack length for WS specimen is expressed as [24]
(18)
a
=
(
h
+
h
0
)
{
1
-
(
13.18
CMOD
E
t
/
P
+
9.16
)
1
/
2
}
-
h
0
.
The empirical expression (18) is valid within 2% accuracy for 0.2 ≤ α ≤ 0.8. The residual Young’s modulus E is calculated using theP-CMOD curve as
(19)
E
=
1
t
c
i
[
13.18
×
(
1
-
a
0
+
h
0
h
+
h
0
)
2
-
9.16
]
,
where c
i = CMOD/P, the segment compliance; t is the specimens thickness; h is the specimens height; and h
0 is the thickness of clip gauge holder. The value of equivalent-elastic crack length a and residual E is listed in Table 1.Table 1
Experimental results of fracture parameters.
Specimen
Temperature
P
ini  (kN)
CMODini  (mm)
P
u  (kN)
CMOD
c  (mm)
CTOD
c  (mm)
E  (GPa)
a
c  (h)
K
Ic
ini   (MPa·m1/2)
K
Ic
un   (MPa·m1/2)
G
F  (N/m)
WS1
20°C
6.19
0.068
8.33
0.174
0.065
15.30
0.53
0.505
1.061
234.15
WS2
6.28
0.047
9.81
0.120
0.039
20.51
0.48
0.523
1.070
483.66
WS3
7.26
0.063
10.40
0.210
0.079
20.66
0.57
0.610
1.497
438.22
WS4
7.02
0.086
7.92
0.152
0.060
18.88
0.56
0.357
1.091
219.39
WS5
5.65
0.060
9.39
0.237
0.096
15.45
0.54
0.503
1.213
321.05
Average
6.55
0.065
9.17
0.178
0.068
18.16
0.54
0.498
1.186
339.30
WS6
65°C
6.98
0.055
11.31
0.195
0.078
21.73
0.56
0.550
1.594
425.91
WS7
3.88
0.050
8.23
0.163
0.100
24.79
0.66
0.303
1.664
482.62
WS8
6.88
0.078
10.41
0.212
0.087
19.43
0.57
0.557
1.518
487.75
WS9
7.94
0.052
10.71
0.164
0.087
23.25
0.60
0.511
1.685
480.51
WS10
6.32
0.056
11.67
0.229
0.086
16.60
0.54
0.562
1.507
522.36
Average
6.42
0.058
10.47
0.193
0.088
21.16
0.59
0.480
1.594
479.83
WS11
120°C
5.03
0.064
8.37
0.191
0.056
10.65
0.47
0.518
0.900
396.52
WS13
4.69
0.093
8.25
0.224
0.084
11.87
0.53
0.417
1.058
517.82
WS12
4.71
0.070
7.53
0.357
0.152
9.48
0.60
0.419
1.202
654.73
WS14
2.79
0.030
7.53
0.198
0.083
15.42
0.58
0.249
1.107
345.46
WS15
—
—
—
—
—
—
—
—
—
—
Average
4.31
0.064
7.92
0.243
0.094
11.86
0.55
0.401
1.067
478.63
WS21
300°C
1.89
0.182
3.40
0.653
0.283
2.45
0.61
0.168
0.556
437.92
WS22
3.48
0.185
5.53
0.667
0.280
3.49
0.59
0.309
0.841
611.47
WS23
1.82
0.121
3.38
0.672
0.271
1.91
0.57
0.162
0.480
341.77
WS24
2.61
0.194
4.97
0.577
0.262
1.99
0.52
0.232
0.589
564.12
WS25
2.03
0.096
4.17
0.651
0.361
4.03
0.68
0.175
0.913
549.99
Average
2.37
0.156
4.29
0.644
0.291
2.78
0.59
0.209
0.676
501.05
WS36
450°C
1.52
0.126
3.37
1.009
0.544
1.41
0.62
0.135
0.582
611.53
WS37
—
—
—
—
—
—
—
—
—
—
WS38
1.52
0.163
3.26
1.419
0.660
1.46
0.62
0.135
0.527
482.45
WS39
1.12
0.296
3.07
1.348
0.617
1.34
0.64
0.100
0.563
663.10
WS40
0.99
0.105
2.94
1.394
0.666
1.58
0.68
0.088
0.659
678.79
Average
1.29
0.172
3.16
1.293
0.622
1.48
0.64
0.115
0.583
608.97
WS46
600°C
0.76
0.443
1.13
1.482
0.684
0.47
0.65
0.067
0.221
228.23
WS47
0.53
0.139
1.48
2.082
0.684
0.48
0.64
0.063
0.277
395.06
WS48
0.81
0.324
1.65
1.908
0.813
1.14
0.76
0.072
0.550
539.22
WS49
0.58
0.436
1.14
1.687
0.973
0.38
0.65
0.052
0.225
331.99
WS50
0.62
0.279
1.48
2.082
0.727
0.38
0.62
0.068
0.213
273.07
Average
0.62
0.324
1.38
1.848
0.799
0.57
0.67
0.064
0.297
353.51


#### 2.4.3. Calculation of Crack Opening Displacement

The crack opening displacement at positionx along the fictitious crack line COD
(
x
) is computed from the known value of CMOD using the following expression [23]:
(20)
COD
(
x
)
=
CMOD
{
(
1
-
x
a
)
2
+
(
1.018
-
1.149
a
h
)
[
x
a
-
(
x
a
)
2
]
}
1
/
2
.
Also, the value of cohesive stress along the fictitious fracture zone corresponding to crack opening displacement at all stages of loading is evaluated using bilinear stress-displacement softening law as given in (3).


### 2.5. Calculation of Double-K Fracture Parameters

The parameters like initiation toughnessK
Ic
ini and stress intensity factors (SIF) are required to be calculated to judge the stability of a propagating crack in a loaded structure using K
R
(
Δ
a
) curve analysis. The LEFM formula for the corresponding test specimen geometries is used for this purpose.The SIF for WS test specimens is written as [24]
(21)
K
(
P
,
a
)
=
P
×
1
0
-
3
t
h
1
/
2
f
(
α
)
,
(22)
f
(
α
)
=
3.675
×
[
1
-
0.12
(
α
-
0.45
)
]
(
1
-
α
)
3
/
2
,
α
=
a
h
.
The empirical expression (21) is valid within 2% accuracy for 0.2 ≤ α ≤ 0.8.


## 2.1. Background

According toK
R-curve criterion [5], the crack extension resistance of a cracked solid consists of the inherent toughness K
Ic
ini and the cohesive toughness K
c
(
Δ
a
) which increases with the increasing amount of crack extension. The cohesive toughness depends upon cohesive stress distribution f
(
σ
) which is a function of crack opening displacement w and tensile strength of concrete f
t and the propagating crack length a. At the onset of unstable crack propagation, the stress intensity factor K at the tip of the propagating crack is expressed as
(1)
K
=
K
R
(
Δ
a
)
,
where K
R
(
Δ
a
) is the crack extension resistance at crack extension length Δ
a
=
a
-
a
0. Also K
R
(
Δ
a
) is expressed in the following relation:
(2)
K
R
(
Δ
a
)
=
K
Ic
ini
+
K
c
(
Δ
a
)
,
(3)
K
c
(
Δ
a
)
=
F
1
(
f
t
,
f
(
σ
)
,
Δ
a
)
.In order to develop theK
R-curve for complete fracture process considering the cohesive stress in fictitious fracture zone, the value of cohesive toughness K
c at every stage of loading is important to determine. During crack propagation, four different stages are considered with the help of three characteristic crack lengths (a
0, a
c, and a
w
c) as represented in Figure 1 in which r
u is the undamaged portion of the ligament, a
0 is initial crack length, a
c is crack length at critical condition of unstable crack propagation, and a
w
c is the length of fully developed fictitious fracture zone after which the stress-free crack propagation will begin.Four stages of crack propagation during the fracture process.
(a)
Initial notch length
(b)
Critical notch extension
(c)
Fully developed crack zone
(d)
Formation of new stress-free crack length


## 2.2. Softening Traction-Separation Law of Postfire Concrete

The softening traction-separation law is a priori to determine theK
R-curve, at room temperature, and many expressions have been proposed based on direct tensile test [8, 17–20]. Based on numerical studies, simplified bilinear expressions for the softening traction-separation law (Figure 2) were suggested by Petersson in 1981 [17], Hilsdorf and Brameshuber in 1991 [19], and Phillips and Zhang in 1993. The area under the softening curve was defined as the fracture energy G
F by Hillerborg et al. in 1976 [21]. Therefore, one can get the following equation:
(4)
G
F
=
1
2
(
f
t
w
s
+
σ
s
w
0
)
.Figure 2
The bilinear softening traction-separation law.As a consequence, a general form of the simplified bilinear expression of the softening traction-separation law is given as follows:(5)
σ
=
{
f
t
-
(
f
t
-
σ
s
)
w
w
0
0
≤
w
≤
w
s
,
σ
s
(
w
0
-
w
)
(
w
0
-
w
s
)
w
s
≤
w
≤
w
0
,
0
w
≥
w
0
.Different values of the break point (σ
s, w
s) and the crack width w
0 at stress-free point were used for the expression proposed by different researchers. In present work, the bilinear softening function of concrete proposed by Petersson is used for postfire specimens as follows:
(6)
σ
s
=
f
t
3
,
w
s
=
0.8
G
F
f
t
,
w
0
=
3.6
G
F
f
t
.


## 2.3. Analytical Method

The standard Green function [22] for the edge cracks with finite width of plate subjected to a pair of normal forces is used to evaluate the value of cohesive toughness. The general expression for the crack extension resistance for complete fracture associated with cohesive stress distribution in the fictitious fracture zone for Mode I fracture is given as follows:
(7)
K
R
(
Δ
a
)
=
K
Ic
ini
+
∫
a
0
a
2
σ
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
,
where
(8)
F
(
x
a
,
a
h
)
=
3.52
(
1
-
x
/
a
)
(
1
-
a
/
h
)
3
/
2
-
4.35
-
5.28
x
/
a
(
1
-
a
/
h
)
+
{
1.30
-
0.30
(
x
/
a
)
3
/
2
1
-
(
x
/
a
)
2
+
0.83
-
1.76
x
a
}
×
{
1
-
(
1
-
x
a
)
a
h
}
.
Equation (7) is employed according to the conditions of four stages of crack propagation as mentioned below.


### 2.3.1. Case 1: Whena
=
a
0

There is no advancement in the initial notch length at this stage of loading and the body remains in elastic condition, subjected to small load (up toP
ini) without any slow crack growth. Hence, cohesive stress σ
(
x
) = 0; the crack growth resistance remains equal to initiation toughness of the material. From (7), it is expressed as
(9)
K
R
=
K
Ic
ini
.


### 2.3.2. Case 2: Whena
0
≤
a
≤
a
c

The stable slow crack growth will take place until the effective crack extensiona
c corresponding to the maximum load P
u is achieved. The cohesive force will start acting across the fictitious fracture zone resulting in the increase of crack extension resistance.(a) For specimens subjected to temperatures less than 120°C, the critical CTODc corresponding to maximum load P
u is less than w
s, as shown in Figure 3(a). The distribution of cohesive stress along the fictitious fracture zone is approximated to be linear, as shown in Figure 4(a). The variation of cohesive stress along the fictitious fracture zone for this loading condition, that is, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(10)
σ
(
x
)
=
σ
(
w
)
+
(
f
t
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
-
a
0
)
,
where σ
(
w
) and w are the values of cohesive stress and crack opening displacement, being at the tip of initialnotch, respectively. The value of σ
(
w
) is determined by using bilinear softening function as follows:
(11)
σ
(
w
)
=
σ
s
(
w
s
)
+
w
s
-
w
w
s
(
f
t
-
σ
s
(
w
s
)
)
.
The crack extension resistance in this case is evaluated using (7) and (8).Two different situations for CTODc and w
s.
(a)
When CTODc
≤
w
s
(b)
When CTODc
>
w
sThe cohesive distribution of crack propagation stage II.
(a)
The linear distribution of cohesive force
(b)
The bilinear distribution of cohesive force(b) For specimens subjected to temperatures higher than 120°C, the critical CTODc corresponding to maximum load P
u is wider than w
s, as shown in Figure 3(b). The distribution of cohesive stress along the fictitious fracture zone is approximated to be bilinear, as shown in Figure 4(b). The variation of cohesive stress along the fictitious fracture zone for this loading condition, also, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(12)
σ
1
(
x
)
=
σ
(
w
)
+
(
σ
s
(
w
s
)
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
s
-
a
0
)
a
s
≤
x
≤
a
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.
The value of σ
(
w
) is determined by using bilinear softening function as follows:
(13)
σ
(
w
)
=
w
0
-
w
w
0
-
w
s
σ
s
(
w
s
)
.
The limits of integration of (7) should be taken in two steps: a
0
≤
x
≤
a
s for cohesive stress σ
1
(
x
) and a
s
≤
x
≤
a for cohesive stress σ
2
(
x
), respectively. The same Green’s function F
(
x
/
a
,
a
/
h
) for a given effective crack extension will be determined using (8). The calculated formula is listed as follows:
(14)
K
R
=
K
Ic
ini
+
∫
a
0
a
s
2
σ
2
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
+
∫
a
s
a
2
σ
3
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
.
The effective crack length at break point a
s (as shown in Figure 4(b)) is computed from the following nonlinear expression [23] by substituting COD
(
a
s
), CMOD, a, and h:
(15)
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
,
where COD
(
a
s
) is the crack opening displacement at a
s, a is the effective crack length (according to (18)), and h is the specimen height.


### 2.3.3. Case 3: Whena
c
≤
a
≤
a
w
0

During this stage of the applied load for all temperatures, corresponding CTOD and effective crack length have increased more than maximum loadP
u, CTODc, and critical effective crack extension a
c, respectively. The cohesive stress distribution for this case will take a bilinear shape, as shown in Figure 4(b). It is the same situation as case 2 part (b), and the crack extension resistance K
R at this stage would be calculated according to (14). This case is shown in Figure 5.Figure 5
Crack propagation stage III.


### 2.3.4. Case 4: Whena
≥
a
w
0

This situation of loading corresponds to the descending portion ofP-CMOD curve. At the effective crack extension, a
w
0, full shape of cohesive stress distribution is allowed to develop, and beyond further extension in crack during loading, a new stress-free crack in front of the initial notch tip will form. This case is shown in Figure 6 and the stress distribution is expressed with the following relation:
(16)
σ
1
(
x
)
=
0
a
0
≤
x
≤
a
w
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
(
x
-
a
w
0
)
(
a
s
-
a
w
0
)
a
s
≤
x
≤
a
w
0
,
σ
3
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.Figure 6
Crack propagation stage IV.Similar to (15), the effective crack length a
w
0 corresponding to zero stress of new fictitious fracture zone as shown in Figure 6 is computed from the following nonlinear expression:(17)
COD
(
a
w
0
)
=
CMOD
{
(
1
-
a
w
0
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
w
0
a
-
(
a
w
0
a
)
2
]
}
1
/
2
,
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
.
The evaluation of crack extension resistance in this case is also done using (7) and (8).


## 2.3.1. Case 1: Whena
=
a
0

There is no advancement in the initial notch length at this stage of loading and the body remains in elastic condition, subjected to small load (up toP
ini) without any slow crack growth. Hence, cohesive stress σ
(
x
) = 0; the crack growth resistance remains equal to initiation toughness of the material. From (7), it is expressed as
(9)
K
R
=
K
Ic
ini
.


## 2.3.2. Case 2: Whena
0
≤
a
≤
a
c

The stable slow crack growth will take place until the effective crack extensiona
c corresponding to the maximum load P
u is achieved. The cohesive force will start acting across the fictitious fracture zone resulting in the increase of crack extension resistance.(a) For specimens subjected to temperatures less than 120°C, the critical CTODc corresponding to maximum load P
u is less than w
s, as shown in Figure 3(a). The distribution of cohesive stress along the fictitious fracture zone is approximated to be linear, as shown in Figure 4(a). The variation of cohesive stress along the fictitious fracture zone for this loading condition, that is, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(10)
σ
(
x
)
=
σ
(
w
)
+
(
f
t
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
-
a
0
)
,
where σ
(
w
) and w are the values of cohesive stress and crack opening displacement, being at the tip of initialnotch, respectively. The value of σ
(
w
) is determined by using bilinear softening function as follows:
(11)
σ
(
w
)
=
σ
s
(
w
s
)
+
w
s
-
w
w
s
(
f
t
-
σ
s
(
w
s
)
)
.
The crack extension resistance in this case is evaluated using (7) and (8).Two different situations for CTODc and w
s.
(a)
When CTODc
≤
w
s
(b)
When CTODc
>
w
sThe cohesive distribution of crack propagation stage II.
(a)
The linear distribution of cohesive force
(b)
The bilinear distribution of cohesive force(b) For specimens subjected to temperatures higher than 120°C, the critical CTODc corresponding to maximum load P
u is wider than w
s, as shown in Figure 3(b). The distribution of cohesive stress along the fictitious fracture zone is approximated to be bilinear, as shown in Figure 4(b). The variation of cohesive stress along the fictitious fracture zone for this loading condition, also, a
0
≤
a
≤
a
c or 0 ≤ CTOD ≤ CTODc, is written as
(12)
σ
1
(
x
)
=
σ
(
w
)
+
(
σ
s
(
w
s
)
-
σ
(
w
)
)
(
x
-
a
0
)
(
a
s
-
a
0
)
a
s
≤
x
≤
a
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.
The value of σ
(
w
) is determined by using bilinear softening function as follows:
(13)
σ
(
w
)
=
w
0
-
w
w
0
-
w
s
σ
s
(
w
s
)
.
The limits of integration of (7) should be taken in two steps: a
0
≤
x
≤
a
s for cohesive stress σ
1
(
x
) and a
s
≤
x
≤
a for cohesive stress σ
2
(
x
), respectively. The same Green’s function F
(
x
/
a
,
a
/
h
) for a given effective crack extension will be determined using (8). The calculated formula is listed as follows:
(14)
K
R
=
K
Ic
ini
+
∫
a
0
a
s
2
σ
2
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
+
∫
a
s
a
2
σ
3
(
x
)
F
(
x
/
a
,
a
/
h
)
π
a
d
x
.
The effective crack length at break point a
s (as shown in Figure 4(b)) is computed from the following nonlinear expression [23] by substituting COD
(
a
s
), CMOD, a, and h:
(15)
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
,
where COD
(
a
s
) is the crack opening displacement at a
s, a is the effective crack length (according to (18)), and h is the specimen height.


## 2.3.3. Case 3: Whena
c
≤
a
≤
a
w
0

During this stage of the applied load for all temperatures, corresponding CTOD and effective crack length have increased more than maximum loadP
u, CTODc, and critical effective crack extension a
c, respectively. The cohesive stress distribution for this case will take a bilinear shape, as shown in Figure 4(b). It is the same situation as case 2 part (b), and the crack extension resistance K
R at this stage would be calculated according to (14). This case is shown in Figure 5.Figure 5
Crack propagation stage III.


## 2.3.4. Case 4: Whena
≥
a
w
0

This situation of loading corresponds to the descending portion ofP-CMOD curve. At the effective crack extension, a
w
0, full shape of cohesive stress distribution is allowed to develop, and beyond further extension in crack during loading, a new stress-free crack in front of the initial notch tip will form. This case is shown in Figure 6 and the stress distribution is expressed with the following relation:
(16)
σ
1
(
x
)
=
0
a
0
≤
x
≤
a
w
0
,
σ
2
(
x
)
=
σ
s
(
w
s
)
(
x
-
a
w
0
)
(
a
s
-
a
w
0
)
a
s
≤
x
≤
a
w
0
,
σ
3
(
x
)
=
σ
s
(
w
s
)
+
(
f
t
-
σ
s
(
w
s
)
)
(
x
-
a
s
)
(
a
-
a
s
)
a
s
≤
x
≤
a
.Figure 6
Crack propagation stage IV.Similar to (15), the effective crack length a
w
0 corresponding to zero stress of new fictitious fracture zone as shown in Figure 6 is computed from the following nonlinear expression:(17)
COD
(
a
w
0
)
=
CMOD
{
(
1
-
a
w
0
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
w
0
a
-
(
a
w
0
a
)
2
]
}
1
/
2
,
COD
(
a
s
)
=
CMOD
{
(
1
-
a
s
a
)
2
+
(
1.018
-
1.149
a
h
)
×
[
a
s
a
-
(
a
s
a
)
2
]
}
1
/
2
.
The evaluation of crack extension resistance in this case is also done using (7) and (8).


## 2.4. Determination of Equivalent Crack Extension

### 2.4.1. Assumption

Linear asymptotic superposition assumption [5, 6] is considered to introduce the nonlinearity effect inP-CMOD curves during loading of concrete test specimens. This assumption enables us to introduce linear elastic fracture mechanics (LEFM) for calculating fracture parameters during every stage of loading in concrete structures. The hypotheses of the assumption are given as follows.(1)
The nonlinear characteristic of theP-CMOD curve is caused by fictitious crack extension in front of a stress-free crack.
(2)
An effective crack consists of an equivalent-elastic stress-free crack and equivalent-elastic fictitious crack extension.


### 2.4.2. Equivalent Crack Extension for WS Specimen

The equivalent-elastic crack length for WS specimen is expressed as [24]
(18)
a
=
(
h
+
h
0
)
{
1
-
(
13.18
CMOD
E
t
/
P
+
9.16
)
1
/
2
}
-
h
0
.
The empirical expression (18) is valid within 2% accuracy for 0.2 ≤ α ≤ 0.8. The residual Young’s modulus E is calculated using theP-CMOD curve as
(19)
E
=
1
t
c
i
[
13.18
×
(
1
-
a
0
+
h
0
h
+
h
0
)
2
-
9.16
]
,
where c
i = CMOD/P, the segment compliance; t is the specimens thickness; h is the specimens height; and h
0 is the thickness of clip gauge holder. The value of equivalent-elastic crack length a and residual E is listed in Table 1.Table 1
Experimental results of fracture parameters.
Specimen
Temperature
P
ini  (kN)
CMODini  (mm)
P
u  (kN)
CMOD
c  (mm)
CTOD
c  (mm)
E  (GPa)
a
c  (h)
K
Ic
ini   (MPa·m1/2)
K
Ic
un   (MPa·m1/2)
G
F  (N/m)
WS1
20°C
6.19
0.068
8.33
0.174
0.065
15.30
0.53
0.505
1.061
234.15
WS2
6.28
0.047
9.81
0.120
0.039
20.51
0.48
0.523
1.070
483.66
WS3
7.26
0.063
10.40
0.210
0.079
20.66
0.57
0.610
1.497
438.22
WS4
7.02
0.086
7.92
0.152
0.060
18.88
0.56
0.357
1.091
219.39
WS5
5.65
0.060
9.39
0.237
0.096
15.45
0.54
0.503
1.213
321.05
Average
6.55
0.065
9.17
0.178
0.068
18.16
0.54
0.498
1.186
339.30
WS6
65°C
6.98
0.055
11.31
0.195
0.078
21.73
0.56
0.550
1.594
425.91
WS7
3.88
0.050
8.23
0.163
0.100
24.79
0.66
0.303
1.664
482.62
WS8
6.88
0.078
10.41
0.212
0.087
19.43
0.57
0.557
1.518
487.75
WS9
7.94
0.052
10.71
0.164
0.087
23.25
0.60
0.511
1.685
480.51
WS10
6.32
0.056
11.67
0.229
0.086
16.60
0.54
0.562
1.507
522.36
Average
6.42
0.058
10.47
0.193
0.088
21.16
0.59
0.480
1.594
479.83
WS11
120°C
5.03
0.064
8.37
0.191
0.056
10.65
0.47
0.518
0.900
396.52
WS13
4.69
0.093
8.25
0.224
0.084
11.87
0.53
0.417
1.058
517.82
WS12
4.71
0.070
7.53
0.357
0.152
9.48
0.60
0.419
1.202
654.73
WS14
2.79
0.030
7.53
0.198
0.083
15.42
0.58
0.249
1.107
345.46
WS15
—
—
—
—
—
—
—
—
—
—
Average
4.31
0.064
7.92
0.243
0.094
11.86
0.55
0.401
1.067
478.63
WS21
300°C
1.89
0.182
3.40
0.653
0.283
2.45
0.61
0.168
0.556
437.92
WS22
3.48
0.185
5.53
0.667
0.280
3.49
0.59
0.309
0.841
611.47
WS23
1.82
0.121
3.38
0.672
0.271
1.91
0.57
0.162
0.480
341.77
WS24
2.61
0.194
4.97
0.577
0.262
1.99
0.52
0.232
0.589
564.12
WS25
2.03
0.096
4.17
0.651
0.361
4.03
0.68
0.175
0.913
549.99
Average
2.37
0.156
4.29
0.644
0.291
2.78
0.59
0.209
0.676
501.05
WS36
450°C
1.52
0.126
3.37
1.009
0.544
1.41
0.62
0.135
0.582
611.53
WS37
—
—
—
—
—
—
—
—
—
—
WS38
1.52
0.163
3.26
1.419
0.660
1.46
0.62
0.135
0.527
482.45
WS39
1.12
0.296
3.07
1.348
0.617
1.34
0.64
0.100
0.563
663.10
WS40
0.99
0.105
2.94
1.394
0.666
1.58
0.68
0.088
0.659
678.79
Average
1.29
0.172
3.16
1.293
0.622
1.48
0.64
0.115
0.583
608.97
WS46
600°C
0.76
0.443
1.13
1.482
0.684
0.47
0.65
0.067
0.221
228.23
WS47
0.53
0.139
1.48
2.082
0.684
0.48
0.64
0.063
0.277
395.06
WS48
0.81
0.324
1.65
1.908
0.813
1.14
0.76
0.072
0.550
539.22
WS49
0.58
0.436
1.14
1.687
0.973
0.38
0.65
0.052
0.225
331.99
WS50
0.62
0.279
1.48
2.082
0.727
0.38
0.62
0.068
0.213
273.07
Average
0.62
0.324
1.38
1.848
0.799
0.57
0.67
0.064
0.297
353.51


### 2.4.3. Calculation of Crack Opening Displacement

The crack opening displacement at positionx along the fictitious crack line COD
(
x
) is computed from the known value of CMOD using the following expression [23]:
(20)
COD
(
x
)
=
CMOD
{
(
1
-
x
a
)
2
+
(
1.018
-
1.149
a
h
)
[
x
a
-
(
x
a
)
2
]
}
1
/
2
.
Also, the value of cohesive stress along the fictitious fracture zone corresponding to crack opening displacement at all stages of loading is evaluated using bilinear stress-displacement softening law as given in (3).


## 2.4.1. Assumption

Linear asymptotic superposition assumption [5, 6] is considered to introduce the nonlinearity effect inP-CMOD curves during loading of concrete test specimens. This assumption enables us to introduce linear elastic fracture mechanics (LEFM) for calculating fracture parameters during every stage of loading in concrete structures. The hypotheses of the assumption are given as follows.(1)
The nonlinear characteristic of theP-CMOD curve is caused by fictitious crack extension in front of a stress-free crack.
(2)
An effective crack consists of an equivalent-elastic stress-free crack and equivalent-elastic fictitious crack extension.


## 2.4.2. Equivalent Crack Extension for WS Specimen

The equivalent-elastic crack length for WS specimen is expressed as [24]
(18)
a
=
(
h
+
h
0
)
{
1
-
(
13.18
CMOD
E
t
/
P
+
9.16
)
1
/
2
}
-
h
0
.
The empirical expression (18) is valid within 2% accuracy for 0.2 ≤ α ≤ 0.8. The residual Young’s modulus E is calculated using theP-CMOD curve as
(19)
E
=
1
t
c
i
[
13.18
×
(
1
-
a
0
+
h
0
h
+
h
0
)
2
-
9.16
]
,
where c
i = CMOD/P, the segment compliance; t is the specimens thickness; h is the specimens height; and h
0 is the thickness of clip gauge holder. The value of equivalent-elastic crack length a and residual E is listed in Table 1.Table 1
Experimental results of fracture parameters.
Specimen
Temperature
P
ini  (kN)
CMODini  (mm)
P
u  (kN)
CMOD
c  (mm)
CTOD
c  (mm)
E  (GPa)
a
c  (h)
K
Ic
ini   (MPa·m1/2)
K
Ic
un   (MPa·m1/2)
G
F  (N/m)
WS1
20°C
6.19
0.068
8.33
0.174
0.065
15.30
0.53
0.505
1.061
234.15
WS2
6.28
0.047
9.81
0.120
0.039
20.51
0.48
0.523
1.070
483.66
WS3
7.26
0.063
10.40
0.210
0.079
20.66
0.57
0.610
1.497
438.22
WS4
7.02
0.086
7.92
0.152
0.060
18.88
0.56
0.357
1.091
219.39
WS5
5.65
0.060
9.39
0.237
0.096
15.45
0.54
0.503
1.213
321.05
Average
6.55
0.065
9.17
0.178
0.068
18.16
0.54
0.498
1.186
339.30
WS6
65°C
6.98
0.055
11.31
0.195
0.078
21.73
0.56
0.550
1.594
425.91
WS7
3.88
0.050
8.23
0.163
0.100
24.79
0.66
0.303
1.664
482.62
WS8
6.88
0.078
10.41
0.212
0.087
19.43
0.57
0.557
1.518
487.75
WS9
7.94
0.052
10.71
0.164
0.087
23.25
0.60
0.511
1.685
480.51
WS10
6.32
0.056
11.67
0.229
0.086
16.60
0.54
0.562
1.507
522.36
Average
6.42
0.058
10.47
0.193
0.088
21.16
0.59
0.480
1.594
479.83
WS11
120°C
5.03
0.064
8.37
0.191
0.056
10.65
0.47
0.518
0.900
396.52
WS13
4.69
0.093
8.25
0.224
0.084
11.87
0.53
0.417
1.058
517.82
WS12
4.71
0.070
7.53
0.357
0.152
9.48
0.60
0.419
1.202
654.73
WS14
2.79
0.030
7.53
0.198
0.083
15.42
0.58
0.249
1.107
345.46
WS15
—
—
—
—
—
—
—
—
—
—
Average
4.31
0.064
7.92
0.243
0.094
11.86
0.55
0.401
1.067
478.63
WS21
300°C
1.89
0.182
3.40
0.653
0.283
2.45
0.61
0.168
0.556
437.92
WS22
3.48
0.185
5.53
0.667
0.280
3.49
0.59
0.309
0.841
611.47
WS23
1.82
0.121
3.38
0.672
0.271
1.91
0.57
0.162
0.480
341.77
WS24
2.61
0.194
4.97
0.577
0.262
1.99
0.52
0.232
0.589
564.12
WS25
2.03
0.096
4.17
0.651
0.361
4.03
0.68
0.175
0.913
549.99
Average
2.37
0.156
4.29
0.644
0.291
2.78
0.59
0.209
0.676
501.05
WS36
450°C
1.52
0.126
3.37
1.009
0.544
1.41
0.62
0.135
0.582
611.53
WS37
—
—
—
—
—
—
—
—
—
—
WS38
1.52
0.163
3.26
1.419
0.660
1.46
0.62
0.135
0.527
482.45
WS39
1.12
0.296
3.07
1.348
0.617
1.34
0.64
0.100
0.563
663.10
WS40
0.99
0.105
2.94
1.394
0.666
1.58
0.68
0.088
0.659
678.79
Average
1.29
0.172
3.16
1.293
0.622
1.48
0.64
0.115
0.583
608.97
WS46
600°C
0.76
0.443
1.13
1.482
0.684
0.47
0.65
0.067
0.221
228.23
WS47
0.53
0.139
1.48
2.082
0.684
0.48
0.64
0.063
0.277
395.06
WS48
0.81
0.324
1.65
1.908
0.813
1.14
0.76
0.072
0.550
539.22
WS49
0.58
0.436
1.14
1.687
0.973
0.38
0.65
0.052
0.225
331.99
WS50
0.62
0.279
1.48
2.082
0.727
0.38
0.62
0.068
0.213
273.07
Average
0.62
0.324
1.38
1.848
0.799
0.57
0.67
0.064
0.297
353.51


## 2.4.3. Calculation of Crack Opening Displacement

The crack opening displacement at positionx along the fictitious crack line COD
(
x
) is computed from the known value of CMOD using the following expression [23]:
(20)
COD
(
x
)
=
CMOD
{
(
1
-
x
a
)
2
+
(
1.018
-
1.149
a
h
)
[
x
a
-
(
x
a
)
2
]
}
1
/
2
.
Also, the value of cohesive stress along the fictitious fracture zone corresponding to crack opening displacement at all stages of loading is evaluated using bilinear stress-displacement softening law as given in (3).


## 2.5. Calculation of Double-K Fracture Parameters

The parameters like initiation toughnessK
Ic
ini and stress intensity factors (SIF) are required to be calculated to judge the stability of a propagating crack in a loaded structure using K
R
(
Δ
a
) curve analysis. The LEFM formula for the corresponding test specimen geometries is used for this purpose.The SIF for WS test specimens is written as [24]
(21)
K
(
P
,
a
)
=
P
×
1
0
-
3
t
h
1
/
2
f
(
α
)
,
(22)
f
(
α
)
=
3.675
×
[
1
-
0.12
(
α
-
0.45
)
]
(
1
-
α
)
3
/
2
,
α
=
a
h
.
The empirical expression (21) is valid within 2% accuracy for 0.2 ≤ α ≤ 0.8.


## 3. Experimental Research

### 3.1. Experimental Program and Experimental Phenomena

To obtain the completeP-CMOD curves, the wedge-splitting tests were implemented. A total of 50 concrete specimens with the same dimensions 230 × 200 × 200 mm were prepared, the geometry of the specimens is shown in Figure 7 (b
=
200 mm, d
=
65 mm, h
=
200 mm, f
=
30 mm, a
0 = 80 mm, and θ = 15°). The concrete mix ratios (by weight) were Cement : Sand : Coarse aggregate : Water = 1.00 : 3.44 : 4.39 : 0.80, with common Portland cement-mixed medium sand and 16 mm graded coarse aggregate. All the specimens had a precast notch of 80 mm height and 3 mm thickness, achieved by placing a piece of steel plate into the molds prior to casting. Each wedge splitting specimen was embedded with a thermal couple in the center of the specimens for temperature control.Figure 7
Geometry of specimens.Nine heating temperatures, ranging from 65°C to 600°C (T
m = 65°C, 120°C, 200°C, 300°C, 350°C, 400°C, 450°C, 500°C, and 600°C), were adopted with the ambient temperature as a reference. Because it was recognized that the fracture behavior measurements were generally associated with significant scatter, five repetitions were performed for each temperature. An electric furnace with net dimensions 300 × 300 × 900 mm was used for heating. When the designated T
m was reached, the furnace was shut down, and the specimens were naturally cooled for 7 days prior to the test.A closed-loop servocontrolled hydraulic jack with a maximum capacity of 1000 kN was employed to conduct the wedge splitting test. Two clip-on extensometers were suited at the mouth and the tip of the crack to measure the crack mouth opening displacement (CMOD) and crack tip opening displacement (CTOD). To obtain the completeP-CMOD curves (as shown in Figure 8), the test rate was fixed at 0.4 mm/min, such that it took approximately 20 minutes to complete a single test of specimens subjected to less than 300°C and 30 minutes for beyond 300°C.Figure 8
Load versus CMOD curves of specimens with temperatures.


### 3.2. Experimental Results and Discussion

The recorded maximum loadP
u, the recorded crack mouth opening displacement CMO
D
c at P
u, the calculated crack tip opening displacement CTO
D
c based on (20), the initial cracking load P
ini determined by graphical method, the calculated residual Young’s modulus E based on (18), the double-K fracture parameters, that is, K
Ic
ini and K
Ic
un, and the residual fracture energy G
F are summarized in Table 1.Figure8 shows typical complete load-displacement curves for different heating temperatures up to 600°C. The figure shows that the ultimate load P
u decreases significantly with increasing temperatures T
m, whereas the crack mouth opening displacement (CMOD) increases with T
m. The initial slope of ascending branches decreases with heating temperatures and the curves become gradually shorter and more extended.From Table1, it is found that the initial load P
ini, ultimate load P
u, the residual Young’s modulus E, and the double-K fracture parameters decrease with the increasing temperatures, whereas the CMODini, CMODc, CTODc, and a
c
/
h increase with T
m. The G
F sustains a hold-increase-decrease tendency with T
m; the detailed explanation could be found in our previous work [25].


## 3.1. Experimental Program and Experimental Phenomena

To obtain the completeP-CMOD curves, the wedge-splitting tests were implemented. A total of 50 concrete specimens with the same dimensions 230 × 200 × 200 mm were prepared, the geometry of the specimens is shown in Figure 7 (b
=
200 mm, d
=
65 mm, h
=
200 mm, f
=
30 mm, a
0 = 80 mm, and θ = 15°). The concrete mix ratios (by weight) were Cement : Sand : Coarse aggregate : Water = 1.00 : 3.44 : 4.39 : 0.80, with common Portland cement-mixed medium sand and 16 mm graded coarse aggregate. All the specimens had a precast notch of 80 mm height and 3 mm thickness, achieved by placing a piece of steel plate into the molds prior to casting. Each wedge splitting specimen was embedded with a thermal couple in the center of the specimens for temperature control.Figure 7
Geometry of specimens.Nine heating temperatures, ranging from 65°C to 600°C (T
m = 65°C, 120°C, 200°C, 300°C, 350°C, 400°C, 450°C, 500°C, and 600°C), were adopted with the ambient temperature as a reference. Because it was recognized that the fracture behavior measurements were generally associated with significant scatter, five repetitions were performed for each temperature. An electric furnace with net dimensions 300 × 300 × 900 mm was used for heating. When the designated T
m was reached, the furnace was shut down, and the specimens were naturally cooled for 7 days prior to the test.A closed-loop servocontrolled hydraulic jack with a maximum capacity of 1000 kN was employed to conduct the wedge splitting test. Two clip-on extensometers were suited at the mouth and the tip of the crack to measure the crack mouth opening displacement (CMOD) and crack tip opening displacement (CTOD). To obtain the completeP-CMOD curves (as shown in Figure 8), the test rate was fixed at 0.4 mm/min, such that it took approximately 20 minutes to complete a single test of specimens subjected to less than 300°C and 30 minutes for beyond 300°C.Figure 8
Load versus CMOD curves of specimens with temperatures.


## 3.2. Experimental Results and Discussion

The recorded maximum loadP
u, the recorded crack mouth opening displacement CMO
D
c at P
u, the calculated crack tip opening displacement CTO
D
c based on (20), the initial cracking load P
ini determined by graphical method, the calculated residual Young’s modulus E based on (18), the double-K fracture parameters, that is, K
Ic
ini and K
Ic
un, and the residual fracture energy G
F are summarized in Table 1.Figure8 shows typical complete load-displacement curves for different heating temperatures up to 600°C. The figure shows that the ultimate load P
u decreases significantly with increasing temperatures T
m, whereas the crack mouth opening displacement (CMOD) increases with T
m. The initial slope of ascending branches decreases with heating temperatures and the curves become gradually shorter and more extended.From Table1, it is found that the initial load P
ini, ultimate load P
u, the residual Young’s modulus E, and the double-K fracture parameters decrease with the increasing temperatures, whereas the CMODini, CMODc, CTODc, and a
c
/
h increase with T
m. The G
F sustains a hold-increase-decrease tendency with T
m; the detailed explanation could be found in our previous work [25].


## 4. Residual Crack Extension Resistance Curves (K
R-Curves) and Stability Criterion

### 4.1. Crack Extension Resistance at Various Temperatures

Since theK
R-curve is considered as a criterion for complete description of crack propagation in structure, it is regarded as the material properties of the complete fracture process. The procedure to calculate the crack extension resistance curves (K
R-curves) is programmed using the analytical expressions given in Section 3. TheP-CMOD curves shown in Figure 8 are used. The K
R-curves at different temperatures with crack extension from 65°C to 600°C are plotted in Figures 9(a)~9(c).K
R-curves with crack extension length of various temperatures.
(a)
K
R -curves at 65°C
(b)
K
R-curves at 200°C
(c)
K
R -curves at 450°CWhen temperatures are less than 200°C, theK
R-curves have almost the same S-shape, which have a good coincidence with existed literature [5].


### 4.2. The Influence of Temperatures on ResidualK
R-Curves

Figure10 plots the K
R-curves of all temperatures; it is concluded that the K
R-curves decrease with increasing temperatures. Generally, temperatures under 120°C appear not to induce much thermal damage to concrete; the K
R-curves of 20°C, 65°C, and 120°C are close. Between 200°C and 600°C, higher temperatures cause more damage to the concrete, and the K
R-curves drop significantly. When temperature reaches 600°C, the K
R-curve almost increases linearly with a small margin.Figure 10
K
R-curves with crack extension length of all temperatures.


### 4.3. The Stability Analysis of Crack Propagation

TheK
R-curve presents the material properties of the complete fracture process. So, the K
R-curve can be taken as a criterion for describing the crack propagation in a structure or a structural component.On the contrary, the stress intensity factor curve (K-curve) during the crack propagation in the structure or the structural component must be calculated already. At an arbitrary loading stage on a wedge-splitting specimen, the stress intensity factor at the tip of a propagating crack can be evaluated by inserting the load P and the length of the propagating crack a into formula (21). For the complete fracture process, the stress intensity factor at the tip of the propagating crack caused by the external load P can be plotted as a curve using formula (21).Herein, the length of the propagating crack is taken as a horizontal axis, and the crack extension resistanceK
R
(
Δ
a
), the stress intensity factor K, and the corresponding load P are plotted on the vertical axis in Figures 11(a)~11(d). From those figures, the stability analysis of the crack propagation during the complete fracture process could be carried out.Stability analysis of crack propagation of various temperatures.
(a)
Stability analysis of crack propagation at room temperature
(b)
Stability analysis of crack propagation at 120°C
(c)
Stability analysis of crack propagation at 300°C
(d)
Stability analysis of crack propagation at 450°CFrom the initial pointB to the critically unstable pointC during the crack propagation, the toughness of the material increases with increasing crack extension from the inherent toughness K
Ic
ini to the value of the crack extension resistance at the critically unstable point which is denoted by K
R
(
Δ
a
c
). The difference between K
R
(
Δ
a
c
) and K
Ic
ini is caused by the cohesive force on the fictitious crack zone during the crack propagation.A common character of Figures11(a)~11(d) is that when the maximum load P
max
⁡ is reached, the stress intensity factor curve (K-curve) coincides with the crack extension resistance curve (K
R-curve) at the point at which the value of the vertical coordinate is equal to the value of the unstable fracture toughness K
Ic
un. All of them correspond to the related critical crack length a
c.Furthermore, it can be seen that when the curve of the stress intensity factor is lower than the crack extension resistance curve, the crack propagates steadily which can be observed in the region between pointB and pointC. Otherwise, when the K-curve coincides with or is higher than the K
R-curve, the crack propagates unsteadily. The stability analysis can be expressed mathematically as given below:
(23)
K
Ic
ini
<
K
(
P
,
a
)
<
K
Ic
un
crack
propagates
steadily
,
K
(
P
,
a
)
>
K
Ic
un
crack
propagates
unsteadily
.


## 4.1. Crack Extension Resistance at Various Temperatures

Since theK
R-curve is considered as a criterion for complete description of crack propagation in structure, it is regarded as the material properties of the complete fracture process. The procedure to calculate the crack extension resistance curves (K
R-curves) is programmed using the analytical expressions given in Section 3. TheP-CMOD curves shown in Figure 8 are used. The K
R-curves at different temperatures with crack extension from 65°C to 600°C are plotted in Figures 9(a)~9(c).K
R-curves with crack extension length of various temperatures.
(a)
K
R -curves at 65°C
(b)
K
R-curves at 200°C
(c)
K
R -curves at 450°CWhen temperatures are less than 200°C, theK
R-curves have almost the same S-shape, which have a good coincidence with existed literature [5].


## 4.2. The Influence of Temperatures on ResidualK
R-Curves

Figure10 plots the K
R-curves of all temperatures; it is concluded that the K
R-curves decrease with increasing temperatures. Generally, temperatures under 120°C appear not to induce much thermal damage to concrete; the K
R-curves of 20°C, 65°C, and 120°C are close. Between 200°C and 600°C, higher temperatures cause more damage to the concrete, and the K
R-curves drop significantly. When temperature reaches 600°C, the K
R-curve almost increases linearly with a small margin.Figure 10
K
R-curves with crack extension length of all temperatures.


## 4.3. The Stability Analysis of Crack Propagation

TheK
R-curve presents the material properties of the complete fracture process. So, the K
R-curve can be taken as a criterion for describing the crack propagation in a structure or a structural component.On the contrary, the stress intensity factor curve (K-curve) during the crack propagation in the structure or the structural component must be calculated already. At an arbitrary loading stage on a wedge-splitting specimen, the stress intensity factor at the tip of a propagating crack can be evaluated by inserting the load P and the length of the propagating crack a into formula (21). For the complete fracture process, the stress intensity factor at the tip of the propagating crack caused by the external load P can be plotted as a curve using formula (21).Herein, the length of the propagating crack is taken as a horizontal axis, and the crack extension resistanceK
R
(
Δ
a
), the stress intensity factor K, and the corresponding load P are plotted on the vertical axis in Figures 11(a)~11(d). From those figures, the stability analysis of the crack propagation during the complete fracture process could be carried out.Stability analysis of crack propagation of various temperatures.
(a)
Stability analysis of crack propagation at room temperature
(b)
Stability analysis of crack propagation at 120°C
(c)
Stability analysis of crack propagation at 300°C
(d)
Stability analysis of crack propagation at 450°CFrom the initial pointB to the critically unstable pointC during the crack propagation, the toughness of the material increases with increasing crack extension from the inherent toughness K
Ic
ini to the value of the crack extension resistance at the critically unstable point which is denoted by K
R
(
Δ
a
c
). The difference between K
R
(
Δ
a
c
) and K
Ic
ini is caused by the cohesive force on the fictitious crack zone during the crack propagation.A common character of Figures11(a)~11(d) is that when the maximum load P
max
⁡ is reached, the stress intensity factor curve (K-curve) coincides with the crack extension resistance curve (K
R-curve) at the point at which the value of the vertical coordinate is equal to the value of the unstable fracture toughness K
Ic
un. All of them correspond to the related critical crack length a
c.Furthermore, it can be seen that when the curve of the stress intensity factor is lower than the crack extension resistance curve, the crack propagates steadily which can be observed in the region between pointB and pointC. Otherwise, when the K-curve coincides with or is higher than the K
R-curve, the crack propagates unsteadily. The stability analysis can be expressed mathematically as given below:
(23)
K
Ic
ini
<
K
(
P
,
a
)
<
K
Ic
un
crack
propagates
steadily
,
K
(
P
,
a
)
>
K
Ic
un
crack
propagates
unsteadily
.


## 5. Conclusions

The residual crack extension resistanceK
R
(
Δ
a
) associated with cohesive stress distribution in fictitious crack zone is evaluated for wedge-splitting specimens using analytical method. The analytical formulations of the residual K
R
(
Δ
a
) for the complete fracture process are presented which consist of two parts. One part is the initiation toughness of the material K
Ic
ini, and the other part is an overall value of stress intensity factor K
Ic
c caused by the cohesive stress along the fictitious crack zone that increases with the crack extension Δ
a. The distributions of the cohesive forces along the fictitious crack zone for varied loading stages are determined according to the softening traction-separation law of the concrete materials.In the calculation of residual crack extension curves (K
R-curves) of wedge-splitting specimens subjected to different temperatures, it is found that in each temperature, the fracture toughness K
R
(
Δ
a
) increases with increasing crack length Δ
a, while the K
R-curves decrease with increasing temperatures. Generally, temperatures under 120°C appear not to induce much thermal damage to concrete; the K
R-curves of 20°C, 65°C, and 120°C are close. Between 200°C and 600°C, high temperatures cause more thermal damage to the concrete, and the K
R-curves drop significantly. When temperature reaches 600°C, the K
R-curve increases almost linearly with a small margin.For the stability analyses of crack propagation, the comparison betweenK
R-curves and the corresponding stress intensity factor curves (K-curves) shows that K
R-curves can be taken as a crack propagation criterion to judge the stability of the crack in the loaded structures; that is, when the stress intensity factor K
(
P
,
a
) is smaller than the crack extension resistance K
R
(
Δ
a
), the crack propagates steadily. Contrarily, when K
(
P
,
a
) is larger than K
R
(
Δ
a
), the crack propagates unsteadily.



---

*Source: 683756-2014-06-26.xml*