Brittle or QuasiBrittle Fracture of Engineering Materials: Recent Developments and New Challenges
View this Special IssueResearch Article  Open Access
Kequan Yu, Jiangtao Yu, Zhoudao Lu, "Average Fracture Energy for Crack Propagation in Postfire Concrete", Advances in Materials Science and Engineering, vol. 2013, Article ID 143208, 13 pages, 2013. https://doi.org/10.1155/2013/143208
Average Fracture Energy for Crack Propagation in Postfire Concrete
Abstract
Wedgesplitting tests of postfire concrete specimens were carried out in the present research to obtain the loaddisplacement curves. Ten temperatures varying from room temperature to 600°C were employed. In order to calculate the accurate fracture energy, the tails of loaddisplacement curves were best fitted using exponential and power functions. Three fracture energy quantities (fracture energy , stable fracture energy , and unstable fracture energy ) with their variation tendency and their mutual relationship were determined to predict energy consumption for the complete fracture propagation. Additionally, the stable fracture work was also calculated. All these fracture parameters sustain an increasedecrease tendency which means that the fracture property of postfire concrete shares the same tendency.
1. Introduction
Since the application of fracture mechanics to concrete, the energy consumption for crack propagation in concrete has been a popular topic. For concrete, the specific fracture energy has been proven to be a useful parameter in the structure design and fracture behavior modeling. The specific fracture energy of concrete was defined based on a tensile test as the energy absorbed per unit crack area in widening the crack from zero to or beyond the critical value above which no stress can transmit [1]. Based on the workoffracture principle, threepoint bending test [2], compact tension [3], and wedgesplitting method [4] were proposed as alternative methods to determine the specific fracture energy . It is computed as the area under the entire imposed load and loadline displacement curve divided by the projected area of uncracked ligament, so the fracture energy represents the average or nominal energy consumption of concrete for an entire crack propagation process.
The existence of fracture process zone FPZ ahead of a crack is now well accepted. Since the 1970s, it has been known that the evolution of the FPZ undergoes two distinct periods—precritical stable crack growth and unstable fracture process [5]. There is no doubt that crack propagation is accompanied by energy dissipation, and the motive for crack propagation comes from either work provided by the imposed load or released strain energy. Fracture energy is one appropriate consideration to describe the amount of energy consumed during crack propagation process.
It is worth noting that the fracture energy can only represent the amount of average energy dissipation for entire crack propagation from crack initiation to complete failure without characterizing crack stable propagation and unstable fracture periods. So even with , it is still not clear how much energy is dissipated during those two crack extension periods. Xu et al. [6] proposed two new concepts the stable fracture energy and unstable fracture energy to describe fracture responses for different crack propagation periods. It is found that kept constant for different ligament lengths, whereas and showed the apparent size effect. But the accurate calculation of fracture surface remains unsolved. It is known that the true path of crack extension is tortuous, not straight as expected. The projected area underestimates the true fracture area. Hence, these parameters are actually nominal values.
The fracture energy of postfire concrete has been studied by several researchers [7–12]. It is found that the residual fracture energy sustained an increasedecrease tendency with the turning point at approximately 450°C. The increase tendency is due to the energy dissipation of microcracks distributing in the concrete, whereas the thermal damage induced by high temperatures reduces the residual fracture energy. However, in these researches, the influence of loadingdisplacement tail was unknown or not considered.
In present paper, wedgesplitting experiments of under ten temperatures levels varying from 20°C to 600°C and the specimens size of 230 mm × 200 mm × 200 mm with initialnotch depth ratios 0.4 are implemented [12]. Based on the workoffracture idea, the residual fracture energy is calculated considering the influence of loaddisplacement tail. Furthermore, the fracture energy consumption for crack stable extension and unstable extension, that is, , , is investigated. However, the true fracture surface remains undetermined and extremely difficult for specimens subjected to high temperatures, so these three fracture energy parameters are still nominal values. Hence, corresponding to and , the stable and unstable fracture work which neglect the fracture surface and their variation about temperatures are thus determined. From these parameters the fracture properties of postfire concrete could be described.
2. Fracture Energies for an Entire Crack Propagation Period
Based on the global energy balance principle, the work performed by a generalized force on its displacement will be transformed into energy: one part of energy is stored in the body in the form of strain energy and the other part is used for crack propagation. It can be mathematically written as [6]: where = average energy needed for unit crack propagation during period from crack initiation to any instant of time ; = newly formed fracture area for this period; the product of these two quantities represents the energy absorbed for crack propagation during this period; = work performed by the external force during the same period; = increase of elastic strain energy of the body until time ; = crack opening displacement corresponding to . If time is approaching the failure terminal point = 1, where the load drops to zero, then in (1) will reduce to where = RILEM fracture energy; = total energy provided by the external force for crack propagation; = projected fracture area perpendicular to the tensile stress direction; = specimen thickness; = specimen height; = initial crack length.
During the stable crack propagation period, the load will increase up to its peak value with its corresponding displacement reaching ; see Figure 1. So, based on (1), the average energy absorbed during this period can be evaluated by where by naming the stable fracture energy = average energy needed for the crack to grow unit area during the crack stable propagation; = change in fracture area; energy provided for stable crack propagation as shown in Figure 1 by the shaded area; work performed by the external force for crack increases its area ; = = increased elastic strain energy of the body; and = critical crack opening displacement.
From the definition of the stable fracture energy , it is implied that the work from the imposed load is expended in two forms during the crack stable propagation: one part is stored in the body in the form of the elastic strain energy () and the other part is the energy absorbed by the fracture zone mainly to counteract the resistance caused by the cohesive forces along the FPZ. In the same way, the unstable fracture energy, denoted by , is defined based on energy equilibrium as where = energy needed for the crack unstable propagation period and = change in fracture area. During this period, the stored energy is completely released until the deformation reaches its maximum value when the load closes to zero. Similar to the RILEM fracture energy , the establishment of the stable fracture energy and unstable fracture energy rests upon the implicit assumption that no other energy consumption occurs outside the fracture zone.
3. Experimental Program and Experimental Phenomena
In this test, concrete specimens were prepared using an ordinary silicate cement PO. 42.5 produced conforming to the Chinese standard. Coarse aggregate was calcareous crushed stone with a maximum size of 16 mm, and river sand was used as the fine aggregates and its maximum diameter was 5 mm. Details of the mix proportioning (by weight) used for concrete and some mechanical properties are Cement : Sand : Limestone Coarse aggregate : Water : fly ash = 1.00 : 3.44 : 4.39 : 0.80 : 0.26.
Fracture properties of concrete were determined by means of the wedge splitting test [4]. The test setup and geometry of the specimen are schematically represented in Figure 2. Compared to threepoint bending notched beams, the wedgesplitting test has following advantages. For the threepoint bending beams, inaccurate measurement of loadpoint displacement and the selfweight of the specimen could influence the real value of the fracture energy. During the test, beams should be carefully handled due to their heavy weight. However, using the WS specimens, the recorded in a horizontal plane is not affected by the crushing of the specimen at the supports or some other factors. Besides, the WS specimens are simple and easily prepared in laboratories or on site.
A total of 50 concrete specimens with the same dimensions 230 × 200 × 200 mm were prepared; the geometry of the specimens and the test setup are shown in Figure 2 ( mm, mm, mm, mm, mm, and ). 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 specimen for temperature control.
Nine heating temperatures, ranging from 65°C to 600°C ( = 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 was reached, the furnace was shut down, and the specimens were naturally cooled for 7 days prior to the test. It averagely took 50, 95, 135, 182, 218, 254, 294, 342, and 453 minutes for the specimens to reach the final temperatures, respectively (from 65°C to 600°C). Figure 3 shows the typical temperature history for several cases with different maximum temperatures. After heating, microcracks disperse on the specimen surface, especially for temperatures higher than 200°C (see in Figure 4).
The fracture surfaces at different temperature intervals (20°C, 200°C, 350°C, 450°C, and 600°C) are shown in Figure 5, which became lighter but more tortuous with increasing temperatures.
(a) Fracture surface of 20°C
(b) Fracture surface of 600°C
(c) Fracture surfaces at typical temperatures (20°C, 200°C, 350°C, 450°C, and 600°C)
A universal machine with a maximum capacity of 1000 kN was used to conduct the wedge splitting tests. During the test, the vertically applied load and the crack opening displacement along the horizontal load line were simultaneously recorded through a data acquisition system. Referring to Figure 2, the splitting force is actually the horizontal component of the force acting on the bearing. Taking the wedge angle into consideration, its relation with the recorded applied load can be developed based on force equilibrium ignoring the small influence from the friction for different roller bearings:
To obtain the complete  curves (shown in Figure 6), 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.
(a) versus curves of specimens at 120°C
(b) versus curves of specimens at 450°C
(c) versus curves of specimens at 600°C
(d) Typical versus curves of specimens at all temperatures
The fracture of specimen is essentially due to the bending moment caused by the horizontal splitting force , vertical component , and selfweight of the specimen. Two symmetrical supports are placed below the center of gravity of each half of the specimen. In doing so, the influence of the dead weight of the specimen and part of the vertical component force on the calculation of the fracture energies could be counteracted. Each roll axis is fixed at the same horizon as the lower plane of the groove and is very close to the center of gravity of each half of the specimen in the vertical position (shown in Figure 2). Due to the carefulness in the choice of the specimen geometry, the roll axis location, and the placement of the supports, the horizontal force contributes most to fracturing the specimen. Therefore,  curves were directly used in the calculation of the RILEM fracture energy , stable fracture energy , and unstable fracture energy .
For our test results, the  curves could be easily obtained from the monitored  curves and (5). Figure 6 contains the plots of  curves for several temperatures and typical  curves for all temperatures.
From Figure 6(a) to 6(c), it is found that with the increasing of temperature (20°C–600°C), the divergence between the curves for the same temperature is more significant. In particular 600°C, the ultimate load of specimen WS50 is one time higher than the one of specimen WS47. Additionally, the whole loading process is not stable for specimens WS49 (a sudden snapback) due to the thermal damage induced by high temperature. Figure 6(d) shows the typical  curves of all temperatures. The ultimate load decreases significantly with increasing temperatures , whereas the crack mouth opening displacement () increases with . The initial slope of ascending branches decreases with heating temperatures, and the curves become gradually shorter and more extended.
It is found that the ultimate load decreases with the increasing temperatures, whereas the increases with (Figure 8). The average value of decreases from 9.17 kN at ambient temperature to 7.92 kN at 120°C, 4.29 kN at 300°C, 3.16 kN at 450°C, and finally 1.38 kN at 600°C, with a final drop of 85%. The value of increases from 0.178 mm at ambient temperature to 0.352 mm at 200°C, 0.901 mm at 400°C, and 1.848 mm at 600°C, nearly 10 times as the ambient value.
4. Experimental Results and Analysis
4.1. Determination of Residual Fracture Energy
In the calculation of RILEM fracture energy , two things should be carefully considered. One thing is that the work done by the selfweight of loading device should be taken into account. In present experiments, the loading device is attached to the testing machine and the  curve includes the selfweight; hence the work should not be calculated again. The other concern pertains to the tail part of the  curve: the recorded point (, ) is just one point when the experiment stops, not the actual point of when the load drops to zero. Therefore, cutting the loaddeflection tail may lead to noticeable inaccuracy in the true value of fracture energy. Thus, to account for these two factors, the actual RILEM fracture energy based on the workoffracture method becomes (see Figure 9) where = measured work enclosed by the  curve until = ; = fresh fracture area; is the work caused by tail effect part. For , curve fitting technology based on Excel software was used to get its approximate value. In previous research from the test results at ambient temperature [13], it was found out that the descending branch of the  curve after onethird of the peak load could be well described by power function (the coefficient of determination for each curve is close to 1). This study indicates that for specimens subject to no more than 120°C, the power function is more accurate; however, for specimens subject to higher temperatures, exponential function is more suitable (see Figure 9).
For power function , , where , = constants for fitting curves. And then the unrecorded work owning to cutting loaddeflection tail can be written as
For exponential function , , where , = constants for fitting curves. And then the unrecorded work owning to cutting loaddeflection tail can be written as
The parameters of , for 20°C~120°C and , for 200°C~600°C, and the value of , is listed in Table 1.

Some necessary test results are tabulated in Table 1, including the maximum value of the horizontal load and its corresponding crack opening displacement , the endpoint where approaches zero and the crack opening displacement arrives at , the initial slope of , that is, the initial stiffness , and the modulus of elasticity .
Hence, the residual fracture energy described by (1) could be calculated. The values are listed in Table 2 and are shown in Figure 10.

Though the residual fracture energy at each temperature has significant scatter, Figure 10 shows that the average values sustain an increasedecrease tendency with . From 20°C to 450°C, average increases from 339.3 Nm^{−1} to 609 Nm^{−1}, while the temperature reaches 600°C and the fracture energy falls back to 307.8 Nm^{−1}. The detailed explanation would be seen elsewhere [12].
4.2. Determination of Stable Fracture Energy and Unstable Fracture Energy
As an extension of fracture energy , two other energybased fracture characteristics, the stable fracture energy and unstable fracture energy , are proposed to describe fracture responses for different crack propagation periods. The analysis shows that fracture energy is actually the weighed average of and , and and could be regarded as two components of . This is very helpful in understanding the whole crack propagation process from the aspect of energy consumption.
To obtain the values of and , the critical effective crack length should first be determined for the calculation of the fracture areas and . Herein, the value was computed based on the double fracture model [14]: where is the compliance of specimens, is modulus of elasticity, is specimen thickness, is specimen height, and is the thickness of the clip gauge holder. For calculation of critical value of equivalent elastic crack length , the value of crack mouth opening displacement () and are taken as and , respectively. Equation (9) is valid for 0.2–0.8 within 2% accuracy. The value of is reported in Table 1.
In a typical loaddisplacement curve shown in Figure 1, is the initial stiffness of the loaddisplacement curve before the start of the fracture process. At critical state when the load arrives at its peak value and its corresponding displacement reaches , the stiffness may not be the same as the initial . The degradation of stiffness is the result of crack propagation. To make it simple, here, the value of at any value of displacement is assumed to be linear with increased deformation [15]. where = final deformation when the load approaches zero and = displacement before which the stiffness is still kept at . At the origin of  curve, the initial tangent stiffness is . The stiffness becomes degraded beyond the origin point; hence is presumed to be zero. It is clear from (10) that when = 0 and when = , and is assumed to be a linear function with respect to the displacement. With the crack opening displacement = , the stiffness at this point could approximately be obtained through (10). Note that  is a general term for loaddisplacement; in our case, the load is the horizontal force and its corresponding displacement is the . Therefore, with , , and , the approximate degraded stiffness can be gained from (10) and the result is shown in Table 1. Figure 11 shows that the values of stiffness decrease monotonously with temperature due to the thermal damage induced by high temperatures.
The stable fracture energy and unstable fracture energy can be derived by following (11)~(13). The specific values of and are compiled in Table 2:
Figure 12 shows the variation of stable fracture energy with temperatures. Similar to the residual fracture energy , it also keeps an increasedecrease tendency with temperatures. Temperatures less than 120°C appear not to induce much thermal damage to concrete, so the cracking resistance almost keeps constant. The values of stable fracture energy at these temperatures are 213 Nm^{−1}, 186 Nm^{−1}, and 180 Nm^{−1}, respectively. The fracture surfaces tend to be more tortuous between 200°C and 450°C than those observed at lower temperatures (see Figure 4) and there exist several cracks competing to form the final fracture, so more energy was dissipated in these specimens. Additionally, the opening of microcracks on the surface and inside the specimens also dissipates energy (see Figure 3). Finally, higher heating temperatures would cause more micro cracks, dehydration, and decomposition and would degrade the resistance. After 450°C, cracking resistance continuously decreases with .
Table 2 shows that the unstable fracture also sustains an increasedecrease tendency with , and its value is much larger than stable fracture energy. Two reasons for larger values of are provided [5]. First, energy consumption besides the main fracture zone takes places for the whole fracture process. During the unstable fracture process, the energy consumption for plasticity or other nonlinear deformation beyond the main FPZ would be much higher, especially for specimens subjected to high temperatures. The other possible reason lies in the calculation of true fracture area. It is known that the interface (transition zone) between cement and aggregate is the weak link in microstructure for normal concrete, and crack propagation would proceed in the path where the energy needs are least. So the true path of crack extension is tortuous and the higher the temperature is, the more tortuous the fracture surface is (see Figure 4), not straight as expected. The projected area is used in the calculation of fracturing surface or , which underestimates the true fracture area. Since the crack experiences much longer distance for the unstable extension, which leads to a greater underestimation of the calculation of the newly fractured area , the calculated is overestimated.
Moreover, in view of (11) and (12), another expression for the fracture energy with respect to and can be followed: Since , the fracture energy is the weighed average of and . However, for engineering application, the stable crack propagation is considered to be more important, because when the load exceeds the maximum value the whole structure would be in an unstable situation. Hence from Figures 7 and 9, it is concluded that the fracture property of postfire concrete sustains an increasedecrease tendency to 600°C, with a turning point at 450°C.
(a) Curve fitting for specimen of ambient temperature
(b) Curve fitting for specimen of 200°C
Furthermore, as mentioned above, the fracture surfaces of , , are project areas, so , , are nominal fracture energies. To avoid the violence of fracture surface, the variation tendency of stable fracture work is determined (see Figure 10). Similarly, also keeps an increasedecrease tendency at the same turning temperature of 450°C.
5. Conclusions
Energy consumption during an entire crack propagation period has been investigated, including the fracture energy , stable fracture energy , unstable fracture energy , and stable fracture work . The conclusions are as follows.(1)Wedgesplitting tests of ten temperatures levels varying from room temperature to 600°C and the specimen size of 230 mm × 200 mm × 200 mm with initialnotch depth ratios 0.4 have been presented. Complete  curves and the curve tails are obtained using exponential and power functions. For specimens subject to no more than 120°C, the power function is more accurate; for higher temperatures, exponential function is more suitable.(2)Three fracture energy quantities corresponding to different aspects of fracture are proposed. The fracture energy in a general case only represents the average energy dissipation for an entire crack propagation process. , the stable fracture energy, denotes the average energy absorption during crack stable propagation, and is used to characterize the average energy consumption for crack unstable propagation and a higher value of than is observed. is actually the weighed average of and .(3)However, for engineering applications, the stable crack propagation is considered to be more important. From the wedgesplitting tests of different temperatures, it is concluded that , sustain an increasedecrease tendency to 600°C, with a turning point at 450°C. Furthermore, the variation of stable fracture work is determined and shares the same tendency with and (Figure 13). All these three parameters mean that the fracture property of postfire concrete sustains an increasedecrease tendency.
Acknowledgments
The State Key Laboratory of Disaster Reduction in Civil Engineering (SLDRCE09D02) and Young Scientist Project of Natural Science Foundation of China (NSFC) have supported this research.
References
 A. Hillerborg, “The theoretical basis of a method to determine the fracture energy ${G}_{\text{F}}$ of concrete,” Materials and Structures, vol. 18, no. 4, pp. 291–296, 1985. View at: Publisher Site  Google Scholar
 RILEMDraftRecommendation (50FCM) (RILEM), “Determination of the fracture energy of mortar and concrete by means of threepoint bend tests on notched beams,” RILEM Materials and Structures, vol. 18, no. 4, pp. 285–290, 1985. View at: Google Scholar
 ASTM International Standard E39906, “Standard test method for linearelastic method planestrain fracture toughness ${K}_{\text{IC}}$ of Metallic Materials,” 2006. View at: Google Scholar
 E. Brühwiler and F. H. Wittmann, “The wedge splitting test, a new method of performing stable fracture mechanics tests,” Engineering Fracture Mechanics, vol. 35, no. 1–3, pp. 117–125, 1990. View at: Google Scholar
 S. Xu and H. W. Reinhardt, “Determination of doubleK criterion for crack propagation in quasibrittle fracture, Part I: experimental investigation of crack propagation,” International Journal of Fracture, vol. 98, no. 2, pp. 111–149, 1999. View at: Google Scholar
 S. Xu, Y. Zhao, and Z. Wu, “Study on the average fracture energy for crack propagation in concrete,” Journal of Materials in Civil Engineering, vol. 18, no. 6, pp. 817–824, 2006. View at: Publisher Site  Google Scholar
 Z. P. Bazand and P. C. Prat, “Effect of temperatures and humidity on fracture energy of concrete,” ACI Materials Journal, vol. 85, no. 4, pp. 262–271, 1988. View at: Google Scholar
 G. Baker, “The effect of exposure to elevated temperatures on the fracture energy of plain concrete,” Materials and Structures, vol. 29, no. 190, pp. 383–388, 1996. View at: Google Scholar
 B. Zhang, N. Bicanic, C. J. Pearce, and G. Balabanic, “Residual fracture properties of normal and highstrength concrete subject to elevated temperatures,” Magazine of Concrete Research, vol. 52, no. 2, pp. 123–136, 2000. View at: Google Scholar
 C. V. Nielsen and N. Bićanić, “Residual fracture energy of highperformance and normal concrete subject to high temperatures,” Materials and Structures, vol. 36, no. 262, pp. 515–521, 2003. View at: Publisher Site  Google Scholar
 B. Zhang and N. Bicanic, “Fracture energy of highperformance concrete at high temperatures up to 450°C: the effects of heating temperatures and testing conditions (hot and cold),” Magazine of Concrete Research, vol. 58, no. 5, pp. 277–288, 2006. View at: Publisher Site  Google Scholar
 J. Yu, K. Yu, and Z. Lu, “Residual fracture properties of concrete subjected to elevated temperatures,” Materials and Structures, vol. 45, no. 8, pp. 1155–1165, 2012. View at: Publisher Site  Google Scholar
 Y. H. Zhao, The analytical study on the energy in the fracture process of concrete [Ph.D. thesis], Dalian University of Technology, Dalian, China, 2002.
 S. Xu and H. W. Reinhardt, “Crack extension resistance and fracture properties of quasibrittle softening materials like concrete based on the complete process of fracture,” International Journal of Fracture, vol. 92, no. 1, pp. 71–99, 1998. View at: Publisher Site  Google Scholar
 R. K. Navalurkar, C. T. Hsu, S. K. Kim, and M. Wecharatana, “True fracture energy of concrete,” ACI Materials Journal, vol. 96, no. 2, pp. 213–225, 1999. View at: Google Scholar
Copyright
Copyright © 2013 Kequan Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.