Mechanical Performance and Numerical Simulation of Basalt Fiber Reinforced Concrete (BFRC) Using Double-K Fracture Model and Virtual Crack Closure Technique (VCCT)

0is paper mainly investigates the fracture parameters of Basalt Fiber Reinforced Concrete (BFRC) with various fiber lengths and dosages using Double-K fracture model. 0e model was developed by fracture criterion using ABAQUS Virtual Crack Closure Technique (VCCT), and the results of themodel and experiments were compared.0e basalt fiber with length of 6mm and 12mmwas added into concrete in the dosage of 0.0%, 0.1%, 0.2%, 0.3%, 0.4%, and 0.5% by volume of concrete, respectively. Concrete specimens were cast into three dimensions, i.e., 60mm× 180mm× 480mm, 80mm× 240mm× 640mm, and 100mm× 300mm× 800mm. 0en, three-point bending test was conducted on precast-notched beams. 0e load versus cracking mouth opening displacement (PCMOD curve) was developed in order to evaluate cracking and breaking load. 0e initial fracture toughness and unstable fracture toughness were derived from the Double-K fracture model aimed to optimize the fiber length and dosage.0e results showed that the initial fracture toughness and unstable fracture toughness increased first and then decreased with the increase in fiber dosage, and basalt fiber with length of 6mm and dosage of 0.2% performed the best toughening effect on concrete.0e comparison results showed that numerical simulation can better simulate the initiation and propagation of BFRC fractures and achieve the dynamic propagation process of fractures.


Introduction
e concrete panel rock-fill dam is extensively developed due to its low cost, short construction period, and strong adaptability, but concrete panel is vulnerable to cracks, which restricted its application. As a quasibrittle material, concrete has low tensile strength (tensile strength is 1/10 to 1/20 of compressive strength). In the process of casting, the concrete panel is in high risk to yield cracks due to uneven settlement, hydration reaction, and dry shrinkage deformation [1][2][3]. ese cracks can cause seepage and remove fine particles in the panel, which seriously impacts the structural stability, reduces the antiseepage performance and durability of the concrete panel, and becomes a great threat to the safety of the dam [4][5][6][7].
According to previous research, adding fibers, such as steel fiber, carbon fiber, and polypropylene fiber into concrete, can improve the basic mechanical properties, ductility, toughness, crack resistance, and durability of concrete structures [8][9][10]. Steel fiber was a commonly used fiber type to enhance the basic mechanical properties of concrete. Pan et al. found that appropriate addition of steel fiber reduced the fragility of the concrete, improved the shape of the major crack, and increased the toughness and ductility of the concrete matrix [11][12][13]. However, addition of steel fiber raises the weight of concrete itself, which can result in difficulty of construction and differential settlement in terms of mass concrete panels. Moreover, as a water retaining structure, concrete panel perennially contacts with water which makes steel fiber rust easily.
us, it induces microcracks in the panel, causes the panel seepage, and further reduces its durability. Consequently, steel fiber is not suitable in concrete panel [14][15][16].
Carbon fiber has the advantages of high tensile strength, low density, large elastic modulus and strong toughness, so it was believed as a good fiber for improving concrete performance. Tong et al. put forward that proper amount of short-cut carbon fiber can enhance the tensile splitting strength, compressive strength, and flexural strength of concrete [17,18]. Deng and Wang et al. studied the fracture parameters of carbon fiber concrete and concluded that carbon fiber could greatly improve the fracture energy, fracture toughness, and critical crack mouth opening displacement of concrete [19,20]. However, carbon fiber is too expensive for a large amount of concrete panel pouring. erefore, carbon fiber concrete has not been widely adopted in water conservancy construction.
Polypropylene fiber has been widely used in recent years. Liang et al. found that polypropylene fiber enhanced the bending performance of concrete and greatly improved the initiation toughness [21]. Nevertheless, the addition of polypropylene fiber declined the strength of concrete, and polypropylene fiber concrete was less effective in reducing dry shrinkage cracks [22]. It is not feasible for Qinghai area along with strong radiation and large temperature difference. erefore, looking for a more suitable fiber to improve the crack resistance and durability of concrete in Qinghai province is necessary.
In recent years, both domestic and foreign scholars have found that basalt fiber can better improve the performance of concrete. Basalt fiber is formed by melting and drawing of basalt ore at high temperature. It has advantages such as high-temperature resistance, chemical resistance, good mechanical properties, and natural compatibility with concrete [23][24][25][26][27]. Research studies on BFRC have been carried out. Wang et al. found that basalt fiber could improve the splitting tensile and flexural strength of concrete, but it had a delayed effect on the early compressive strength of concrete [28]. Basalt fiber significantly enhanced the freezethaw resistance and durability of concrete under corrosion conditions [29]. Yu  eir results showed that, same as plain concrete, the fracture process of BFRC could also be divided into three stages, i.e., crack nonpropagation, crack stable propagation, and crack unstable propagation. Moreover, basalt fiber greatly improved the fracture parameters [34][35][36]. Wang et al. found that adding basalt fiber had effects on compressive strength, tensile splitting strength, flexural tensile strength and size effect, among which compressive strength had the greatest effect [37]. Based on previous studies, basalt fiber had great improvement on the durability and crack resistance of plain concrete.
Up to date, most of relevant research programs focused on the basic mechanical properties of BFRC. However, it is rarely employed to analyze fracture performance with Double-K fracture model. In this paper, the BFRC fracture parameters were analyzed by the Double-K fracture model. Double-K fracture model was put forward by Xu and Reinhardt, which described the fracture process of concrete by two stress intensity factors: initiation toughness and unstable fracture toughness. It introduced the unstable fracture toughness parameters and believed that the aggregate and mortar in concrete caused the increase of unstable fracture toughness. And the fracture of concrete could be divided into three stages: crack nonpropagation, crack stable propagation, and crack unstable propagation [34,[38][39][40][41][42][43][44]. Ruiz et al. investigated the effects of crack size and bond on Double-K fracture parameters of concrete [45]. e feasibility of using peak loads to determine Double-K parameters was demonstrated by Kumar et al. [46]. Wang et al. studied the Double-K fracture parameters of steel fiber reinforced concrete [47].
e Double-K fracture parameters of carbon fiber reinforced concrete were evaluated by Zhang et al. [48]. However, few researchers have applied Double-K fracture parameters to investigate the fracture performance of BFRC. us, it is significant to study the Double-K fracture parameters of BFRC.
For the constitutive model study of fracture, Bernardi et al. used nonlinear finite element on double-parameter simulation for polypropylene fiber concrete. eir results showed that the simulated results were highly coherent with the experimental results [49]. Chi et al. evaluated the value of damage constitutive model on steel and polypropylene fiber concrete by ABAQUS [50]. Hu et al. found that the deviations of the whole-process simulation on crack initiation, stable propagation, and unstable failure were less than 10% for three-point bending beam using extended finite element method (XFEM) [51]. Roth  e results indicated that the variations of the strain energy release rate from test and simulation were between 0.4% and 11.4% for type I and type II fracture, respectively [53]. Valvo accessed energy release rate of type I, II, and III compound cracks using VCCT [54]. e constitutive model has been widely studied in concrete materials, but VCCT has hardly been used to study the fracture behavior of BFRC. Comparing with XFEM, VCCT can better simulate the crack initiation and crack propagation of concrete and reduce time on calculations. erefore, VCCT was selected to simulate and analyze fracture performance of BFRC. 2 Advances in Civil Engineering

Double-K Fracture Parameters and Calculation.
As a quasibrittle material, BFRC can be composed of base phase and dispersed phase. Base phase is a mixture of cement, water, fly ash, water-reducing agent, air entraining agent, and basalt fiber, while dispersed phase is the aggregate. e fracture failure of BFRC can be divided into three levels, the first level: the bonding surface of mortar and aggregate is a weak surface, and cracks are generated on the bonding surface of mortar and aggregate firstly; the second level: with the development of cracks between mortar and aggregate binding surfaces, cracks spread into mortar and the interface of sand and hardened cement slurry begins to crack; third level: as cracks between sand and hardened cement slurry develop, hardened cement slurry starts to crack, accompanied by the pulling-out failure of fibers. Xu and Reinhardt put forward the Double-K fracture model to evaluate the fracture properties of concrete more accurately. In Double-K fracture model, the fracture of concrete consists of three processes, namely, the nonpropagation crack, the stable-propagation crack, and the unstablepropagation crack. In order to better describe the fracture process, the initiation toughness (K Q IC ) and unstable fracture toughness (K S IC ) were introduced. When K < K Q IC , crack will not propagate; when K Q IC � K, crack starts stable propagation; when K Q IC < K < K S IC , crack is in stable propagation; and when K S IC < K, crack is in unstable propagation. In this paper, Double-K fracture parameters were determined according to P-CMOD curve.

Initiation Toughness.
As a quasibrittle material, the crack propagation starts with linear elastic propagation, where crack has not yet appeared. Next stage is elastoplastic propagation, where the crack is stably propagating until the fracture load. Specimen is damaged in this stage. In P-CMOD curve, the load at turning point from the linear part to the nonlinear part is called the crack initiation load (F Q ). According to DL/T 5332-2005 [55], the initiation toughness can be expressed as where K Q IC is initiation toughness, MPa·m 1/2 ; F Q is crack initiation load, kN; a 0 is length of precast crack, m; t is thickness of specimen, m; h is height of specimen, m; S is midspan distance, m; m is mass of midspan specimen, kg; and g is acceleration of gravity, 9.81 m/s 2 .

Critical Effective Crack Length.
e fracture process of BFRC specimen is not all linear due to the bond strength between aggregates and binder and the bridging effect of basalt fiber after matrix cracking. erefore, the nonlinear characteristic is equivalent to elastic fracture when fracture occurs for calculation of unstable fracture toughness. Hereby, the critical effective crack length can be calculated following equation (3) as per DL/T 5332-2005 [55]: where h 0 is the thickness of the edge sheet of the clip extensometer, m; V C is the critical value of crack opening displacement, μm; and E is the elasticity (GPa) as expressed by where a 0 is the length of precast crack, m, and c i is the ratio of displacement to load in the linear stage of load displacement, μm/kN.

Unstable Fracture Toughness.
When the load reaches the maximum, the specimen is in unstable state. e corresponding fracture toughness at this moment is the unstable fracture toughness. Based on DL/T 5332-2005 [55], unstable fracture toughness can be determined by equations (2) and (5): where K S IC is unstable fracture toughness, MPa·m 1/2 ; F max is maximum load, kN; ac is critical effective crack length, m; t is thickness of specimen, m; h is height of specimen, m; S is midspan distance, m; m is mass of midspan specimen, kg; and g is acceleration of gravity, 9.81 m/s 2 .
In this paper, three-point bending tests were carried out for precast notched BFRC beams with different lengths, admixtures, and specimen sizes. e load versus cracking mouth opening displacement curves (P-CMOD curves) were developed in order to evaluate cracking and fracture load. e initiation toughness and unstable toughness parameters of BFRC were calculated by using the Double-K fracture model. e effects of basalt fiber length and content on fracture properties of BFRC were analyzed as well. Due to the limited test conditions, the entire process of BFRC fracture is difficult to determine. erefore, VCCT technique was used to simulate the fracture performance of BFRC and analyze the entire process of BFRC fracture.
II grade fly ash was used in this paper, and its properties are listed in Table 2 as per DL/T 5055-2007 [60].
Advances in Civil Engineering e properties of used water-reducing admixture and air entraining admixture are shown in Tables 3 and 4 based on DL/T 5100-2014 [61].
Coarse aggregate and fine aggregate are natural aggregates (containing artificial crushed pebbles) from the right bank of Yangqu Hydropower Station. eir properties are summarized in Tables 5 and 6 according to DL/T 5144-2015 [62] and DL/T 5151-2014 [63]. e properties of basalt fibers used in the test are shown in Table 7, and the samples are shown in Figure 1.

Specimen Preparation.
In this study, a forced mixer was used to mix BFRC. Coarse aggregate and fine aggregate were mixed for 30 s first. In order to ensure a uniform distribution of the basalt fibers in concrete, a dry fiber mixing method was used to improve the uniformity of fiber distribution based on [66]. Cement, fly ash, and superplasticizers were poured to mix for another 30 s. en, basalt fibers were spread uniformly into mixer followed by 60 s mixing. After water and AEA adding, all materials were mixed for 60 s at last. Next, the mixed BFRC was cast and then consolidated using a table vibrator, followed by surface finishing. After 24 h, specimens were demolded and cured in standard curing room for 28 days. Before three-point bending test, specimens were cut to preform crack as required using cutting machine. In order to avoid error caused during specimen preparation, all specimens were prepared once.
ree different dimensions of specimens used for threepoint bending test were 60 mm × 180 mm × 480.6 mm, 80 mm × 240 mm × 640.8 mm, and 100 mm × 300 mm × 801 mm, respectively. Basalt fibers of 6 mm length were added to the mixture in accordance with the dosage of 0%, 0.1%, 0.2%, 0.3%, 0.4%, and 0.5%, respectively, the same for the basalt fiber of 12 mm length. e detailed information of specimens is listed in Table 9.

Test Equipment.
ree-point bending test was carried out at continuous loading rate of 0.4 MPa/min using Kaiwen universal hydraulic machine. A linear variable differential transformer (LVDT) sensor was installed vertically at the bottom of beam two centimeters from the midspan to measure vertical displacement. An extensometer was mounted at precast crack incision to determine CMOD. Resistance strain gauges (model: S2150-150AA, resistance: 149.5 ± 0.2%, sensitivity coefficient: 2.032 ± 0.32%, grid width: 5 mm, grid length: 150 mm) were pasted every 2 cm at the top of the precast crack.
e specimen under loading is illustrated in Figure 2.

P-CMOD Curves.
e effects of various basalt fiber dosage on P-CMOD curve were obtained from the tests as shown in Figures 3 and 4 for 6 mm fiber and 12 mm fiber, respectively.
According to the P-CMOD curve analysis, in general, the expansion of concrete cracks can be divided into three stages (stage I, stage II, and stage III). Stage I: at the beginning of load application, the load and displacement were in linear relationship. At this moment, BFRC had not started to crack. In this stage, BFRC elastically expanded. Stage II: with load continuing, the microcracks occurred under steady propagation until the ultimate load was reached. BFRC displayed a nonlinear expansion, with the growth of microcracks. It stepped into matrix fracture Stage III: After ultimate load, unstable fracture took over in crack propagation until the cracks produced over the whole specimen. In this stage, the specimen was completely fractured along with the fibers pulling out. In this test, the fracture process of most tested specimens was exhibited these three stages, while some type A and B specimens only had stage I and stage III. e nonlinear rising stage is not apparent enough. is might be due two reasons: (1) the yield process is relatively short and lose a lot of yield values during force load; (2) the sizes of type A and B specimens are relatively small, and the existed size effect leads to small yield strength for type A and B specimens and the relatively short yield process. In the future tests, displacement loading should be selected to better collect BFRC fracture data and analyze BFRC fracture   Table 10. According to Table 10, the initiation toughness and unstable toughness showed irregular growth. However, compared with PC, the growth rates of the initiation toughness and unstable toughness are 1.023-1.736 and 1.015-1.541, respectively. Because concrete will shrink during the forming process, the shrinkage will cause relative displacement of the concrete matrix. e scattered distribution of basalt fibers could prevent the movement of the matrix aggregate and withstand the tensile stress caused by shrinkage, reducing the stress concentration, preventing the expansion of the crack, and enhancing the fracture toughness of concrete. In the process of fracture failure of basalt fiber concrete, the fiber will be pulled out from the matrix, prolonging the stable expansion stage of the basalt fiber concrete specimen and enhancing the fracture toughness of the concrete.

Initiation Toughness
Discussion. BFRC fracture propagation is a complex process. e defects such as microcracks, voids, and weak contact surface between mortar and aggregates exist in BFRC matrix among which most of them are interfacial cracks and very few are mortar cracks. At the beginning of loading, it exhibits the elastic stage because these defects are stable. In the process of increasing the load until it reaches the critical load, these microcracks, voids, and weak surfaces between mortar and aggregate begin to extend or expand, which shows the linear stage for BFRC fracture. e development process of these defects depends on the magnitude of the crack initiation load and initiation toughness.
us, crack initiation load and initiation toughness play an important role in evaluating crack resistance of concrete materials. Figure 4 shows that with the increase in dosage for 6 mm long basalt fiber from 0.0% to 0.5%, the initiation toughness     Advances in Civil Engineering increased first and decreased after 0.2% dosage. e initiation toughness of BFRC was greater than that of plain concrete. Regardless of specimen dimension, initiation toughness was maximum at 0.2% dosage. However, initiation toughness of specimens A increased most which is up to 61.55%. When the basalt fiber content was 0.4%, the initiation toughness was relatively small. e reason might be that the fiber agglomeration caused defects such as microcracks and voids in the concrete. As the raise of dosage for 12 mm long basalt fiber, initiation toughness increased first and then decreased as well. Compared with the plain concrete, the addition of basalt fiber increased initiation toughness. For different dimensions of specimen, the maximum initiation toughness appeared at different fiber dosage. For instance, the peak values of initiation toughness for A and B specimens are at 0.2% fiber dosage, but it was 0.3% for C specimens. e initiation toughness of A specimens was improved most up to 73.6% among the three dimensions of specimen.
Generally, considering the time of crack development, the optimal dosage of basalt fiber was 0.2%. Regarding the increment of initiation toughness, the BFRC with 6 mm basalt fiber displayed relatively stable increase in initiation toughness. Among the BFRC specimens with 12 mm basalt fiber, only one specimen performed a large raise in initiation toughness probably due to randomness. Overall, the optimal combination of basalt fiber was 6 mm in length and 0.2% in dosage. Basalt fiber could increase the initiation toughness because it had higher strength and ductility than the plain concrete. Basalt fibers scattered and connected with each other in concrete matrix, which could inhibit the development of microcracks, voids, and the weak interface of mortar and aggregate, prolong the linear propagation and increase the initiation toughness. However, it should be noted that if the basalt fibers are distributed inhomogeneously and clustered together, they cannot well restrain the development of microcracks, voids, and weak interface in concrete and even lead to the decline of initiation toughness. is is associated with the voids among the fibers and the very low friction among the fibers.

Unstable Fracture Toughness Discussion.
After linear propagation, BFRC started unstable failure until it reached the ultimate load and failed. Figure 5 presents the trends of unstable fracture toughness with various basalt fiber dosages.
As the load increases, the crack in BFRC ends linear stage and steps into the nonlinear stage. In this stage, the   Note. All specimens are not listed here due to tremendous quantity. e specimen is identified by X-Y-Z, where X represents size of specimen, Y is length of basalt fiber, and Z stands for fiber content. A is the specimen in dimension of 60 mm × 180 mm × 480 mm, B is the specimen in dimension of 80 mm × 240 mm × 640 mm, and C is the specimen in dimension of 100 mm × 300 mm × 800 mm. 6 Advances in Civil Engineering Advances in Civil Engineering defects such as the microcracks, voids, and weak contact surface between substrate and aggregate continue to develop. Cracks and other defects occurred in the mortar, meanwhile sand and interface of hardened cement slurry were damaged, resulting in crack propagation into the slurry. When the slurry damaged, the specimen failed. e fracture failure limit of BFRC depends upon the unstable load and unstable fracture toughness, which plays an important role in the evaluation of fracture performance. Figure 5, similar trends of unstable fracture toughness for both 6 mm and 12 mm fibers can be obtained comparing with Figure 5. It can be observed that the maximum unstable fracture toughness was at 0.3% fiber dosage for A specimens and at 0.2% for B, C specimens when introducing 6 mm fibers. While for the specimens with 12 mm fibers, the peak value of unstable fracture toughness was achieved at 0.2% fiber dosage for A, B specimens and at 0.3% for C specimens. Specimens of A dimension had the most increment up to 62.5% at 0.2% dosage. According to Figure 5, the optimal content of basalt fiber is 0.2% for 12 mm considering the failure of the specimen. However, the long fiber will wrap a large amount of water during mixing process, which will adversely affect slump of concrete and construction. erefore, 6 mm basalt fiber is recommended in practical project. Basalt fiber can increase unstable fracture toughness of concrete. Its mechanism is that basalt fiber can be scattered in the concrete, connecting with each other. When crack propagates to cement slurry and cement damages, accompanied by tensile failure of basalt fiber (>2000 MPa), basalt fiber delays the crack propagation and the slurry destruction owing to large consumption of energy during crack propagation which explains that the fracture process of BFRC presents a certain plastic zone.

Numerical Simulation and Analysis
BFRC is a quasibrittle material with complex fracture mechanism. e three-point bending fracture test in the laboratory is difficult to keep tracking the fracture behaviors during the loading process beyond the ultimate load. With the development of computer technology, the macromechanical properties can be characterized using mesh division, finite element calculation, and fracture behavior simulation for concrete. It is very helpful to access fracture performance of BFRC under the correct numerical simulation by ABAQUS, a commercial software, which was used for the finite element modeling. Virtual crack closure technique (VCCT) was first proposed by Rybicki and Kanninen when they studied linear crack problems [67]. Shivakumar et al. applied VCCT to investigate surface crack of concrete [68]. e accuracy and applicability of VCCT in nonlinear materials have been verified by Cao et al. [69]. Additionally, in finite element analysis, VCCT is almost unaffected by mesh size in the calculation process, and shorten the finite element analysis steps, which greatly improves the calculation efficiency [70]. erefore, VCCT was used to numerically simulate the fracture performance of basalt fiber reinforce concrete.

Virtual Crack Closure Technique (VCCT).
VCCT is based on Irwin energy theory, and its core idea is to assume that the energy released from crack propagation equates to the energy for crack closure. When the crack develops from Figure 6(a) to 6(b), assuming the front shape of the crack is constant, the opening size of the crack after propagation is the same as before. During crack propagation, assume that the energy release rate is G I , and the critical energy release rate required for crack propagation is G IC . When G I > G IC , crack propagates. Equation (6) can be used to assess if crack starts to initiates: If f > 1, crack starts to propagate, where b and d represent the width and length of the unit crack tip. As shown in Figure 7, 3,4 is the normal force between nodes 3 and 4; v 1,2 represents the normal displacement between nodes 1 and 2.
However, this rule of propagation is only suitable for type I fracture. e equivalent strain energy release rate G e has to be introduced. en, the criterion of crack propagation is f � (G e /G eC ) ≥ 1. ABAQUS provides BK criterion, POWER criterion, and Reeder criterion to calculate G e .
BK criterion is the best criterion in this case that critical strain energy release rates on two shear directions are the same. erefore, BK criterion was adopted in the simulation as shown in the following equation: where G eC is equivalent strain energy release rate; G I , G II , and G III are crack strain energy release rates of type I, II, and III     10 Advances in Civil Engineering crack, respectively; and G IC , G IIC , and G IIIC are critical crack strain energy release rates of type I, II, and III crack, respectively. η is usually between 0.5 and 2.5 for damage factor.

Comparison between Experimental and Numerical
Results. e finite element method was used to simulate the experiment in ABAQUS software. Units are mm, N, tonne, MPa, tonne/mm 3 , and MPa·mm 1/2 . Hexahedral mesh was selected. e dimension of model and experiment was exactly the same. In order to better simulate the BFRC fracture process, enhanced VCCT was selected because it can simulate the crack initiation, which is more consistent with the real process. e parameters used in this simulation include    density, elastic modulus, Poisson's ratio, crack initiation energy release rate, fracture energy release rate, and damage factor η. e density was measured by the experiments, and the elastic modulus, the crack initiation energy release rate, and the fracture energy release rate were calculated according to the DL/T5332-2005 [55]. Poisson's ratio was obtained based on the GB 50010-2010 [71]. e damage factor η is 1.75 according to empirical research. e specific parameters are shown in Table 11. e mesh division and stress pattern are shown in Figure 9. It can be seen from Figure 8 that the cracks initiated using VCCT to simulate the three-point bending and the stress were completely symmetrical, which is coherent with the real condition, and stress is concentrated in the three supports. e anticipated result was relatively acceptable as well. e three-point bending test on the specimens with 0.2% dosage of basalt fiber in 6 and 12 mm was selected to be simulated. e comparison results from simulation and experiments are presented in Figure 9.
Based on Figure 9, it indicates that under loading, the curves first rise linearly and then go up nonlinearly until reaching the peak value. e softening phenomenon occurred at some point in the descending stage. Regarding the coherence of the results from experiments and numerical simulation, it is relatively high before the peak value but it becomes low for B-6-20 after peak value. Nevertheless, the length of crack can be obtained up to 0.7 mm for numerical simulation but only about 0.4 mm for experiments. e agreement of critical displacement is not good due to displacement load. In other words, the loading rate is relatively high for displacement load. It is difficult to collect data after peak value.
e optimum dosage of basalt fiber, initiation toughness, and unstable fracture toughness were evaluated. e conclusions are as follows: (1) Basalt fiber can enhance the initiation toughness and unstable fracture toughness of concrete. Moreover, it can be uniformly distributed in the concrete matrix, preventing the development of microcracks, voids, weak interface due to basalt fiber bridging effect and enhancing the crack resistance of concrete. (2) e P-CMOD curves were established by experiments, and the initial crack load was obtained based on the state of crack propagation. e initiation toughness was calculated by initial crack load. Basalt fiber increased initiation toughness. e optimal dosage of basalt fiber was 0.2%, and the toughening effect of 6 mm basalt fiber was better than that of 12 mm basalt fiber. Meanwhile, the slump of BFRC was also good to meet the requirement of engineering construction. (3) Based on P-CMOD curves, when the load reached the maximum, the specimens failure. e buckling load and unstable fracture toughness were determined. Basalt fiber enhanced unstable fracture toughness. e optimal dosage of basalt fiber was 0.2%, while the toughening effect of 12 mm basalt fiber was better compared to 6 mm basalt fiber. (4) e results of three-point bending test for the BFRC were analyzed using VCCT. e overall trend from numerical simulation results and the experiment results showed a high degree in coherence. Numerical simulation can better present the whole process of BFRC fracture. However, the numerical simulation and experimental results have some errors in the critical displacement fitting.

12
Advances in Civil Engineering

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.