Effect of Fractal-Based Fiber Clustering on Tensile Properties of BFRP

To explore the clustering phenomenon of discontinuous fibers in composite materials, this paper deduces the fiber uniform distribution coefficient and analytical expressions of fiber clustering content based on fractal theory and establishes a tensile strength prediction model of fiber/epoxy resin composite materials containing cluster fibers. With basalt fiber/epoxy resin composites (BFRP) as an example, this paper analyzes the tensile strength law of BFRP under fiber clustering effect. The results show that when the fiber volume fraction is constant, the tensile strength of the composite in the presence of agglomerated fibers is only related to the fractal dimension of the circumference and cross-sectional area of the inner fiber agglomerate. The calculated value of the composite tensile strength based on fractal theory is lower than the experimental value, but closer to the experimental value than the approximate method. The research conclusions can provide theoretical support for strength prediction of fiber/epoxy resin composites.


Introduction
Basalt fiber reinforced epoxy resin composites (BFRP) can be formed by mixing basalt fiber into epoxy resin materials at a certain ratio. BFRP materials are widely used in civil engineering, water conservancy, and hydropower and energy engineering fields because of their excellent performance, mechanical properties, and durability. At present, the research on BFRP materials focuses on the mechanical properties of the materials and their components, which are mainly in the following aspects.
First, the performance and mechanism of BFRP materials are analyzed, aimed at improving their mechanical properties and durability and exploring their internal mechanism. For example, [1][2][3][4] used different methods to treat and modify the fiber surface and studied the influence of fiber surface modification on the tensile properties of BFRP. e results showed that the mechanical parameters, such as tensile strength and elastic modulus, were significantly improved after the modification. Shokrieh et al. [5,6] studied the sulfuric acid corrosion resistance of BFRP, analyzed the damage mechanics characteristics under bending load and corrosive medium, and obtained the life curve of BFRP materials. Bulut et al. [7] studied the influence of graphene nanoparticles on the damping and vibration characteristics of BFRP. e results indicated that when the mass fraction of graphene nanoparticles was 0.1% and 0.2%, the damping and vibration characteristics of samples were significantly improved. Chelliah et al. [8] studied the influence of titanium carbide on the wear resistance of BFRP. Experiment results showed that the mass loss rate of BFRP decreased significantly when titanium carbide was added to BFRP materials and that it reached the lowest when the content of titanium carbide was 2%. Sampath [9] analyzed the influence of nano-fly ash filler on the mechanical and chemical corrosion resistance of BFRP. Experiment results demonstrated that nano-fly ash filler could improve the properties of BFRP to varying degrees. Kim et al. [10] modified the fracture properties of BFRP with different concentrations of natural graphite sheets. e results suggested that natural graphite particles act as crack nails, thus improving the mechanical and thermal properties of BFRP.
It can be seen that the performance of BFRP materials is achieved mainly by treating the fiber surface and adding admixtures to the substrate. e working mechanism is as follows. (1) By treating the fiber surface, the bond strength between the fiber and the substrate is increased, and the sliding damage is reduced, so that most fibers show "pulloff" damage. When the fiber could reach the yield strength, the composite will be damaged. us, the fiber utilization rate is improved and the mechanical properties of the composite are improved. (2) By adding some tiny particles into the epoxy resin substrate, the microstructure of the substrate is affected in some way. e defects such as small voids in the substrate are reduced, and the stress concentration is reduced. erefore, the mechanical properties of the substrate are improved, so as to enhance the mechanical properties of the composite.
Second, the properties of BFRP components are studied. It mainly refers to research on the mechanical properties of BFRP material as reinforcement, plate, and net. For example, Zha et al. [11] applied BFRP as reinforcement to concrete beam members and studied the impact bending performance of BFRP reinforced concrete beams. Experiment results indicated that the impact resistance of BFRP reinforced concrete beams was equivalent to that of ordinary reinforced concrete beams, and BFRP reinforced concrete beams could be applied to beam members instead of steel bars. Zhang et al. [12] used BFRP as reinforcement to strengthen concrete column members. ey found that BFRP reinforcement significantly improved the ductility of column members, and the combination of inner and outer rectangular BFRP reinforcement was the best method to improve members' ductility. Rui [13] studied the fatigue bond performance of BFRP reinforcement and RAC by drawing test and revealed the influence of design parameters on bond stiffness, bond strength, and residual bond strength. Kim et al. [14] applied BFRP as grid cloth to the reinforcement and repair of existing concrete members and conducted relevant experiments. e results showed that BFRP grid cloth could effectively delay the failure time of damaged concrete members. Liu [15] applied BFRP as a laminated plate to strengthen prestressed reinforced concrete beams and compared it with the traditional CFRP laminated plate. It was found that BFRP had greater failure strain and strength and showed significant creep fracture behavior. It can be seen that the BFRP component has excellent mechanical properties and enjoys great potentials in the field of existing structural reinforcement engineering.
Based on the above analysis, scholars at home and abroad have extensively studied the properties of BFRP materials and components and obtained fruitful achievements, which played a positive role in promoting the development of BFRP materials. e existing studies are largely based on the modification of fiber surface and substrate microstructure to improve the properties of the composite. However, there is little research on the impact of fiber distribution on the mechanical properties of BFRP. In practice, due to the different molding processes of BFRP, it is difficult to achieve uniform distribution of fibers in the matrix material, and fiber agglomeration may occur [16][17][18]. erefore, based on the "fiber clustering model" proposed in [19], this paper analyzes the geometric shape of the cluster fiber section on the fracture surface of BFRP materials by using the fractal theory and obtains the actual fiber uniform distribution coefficient and the clustering fiber content of BFRP materials. Besides, this paper compares the above values with those obtained from the approximate calculation method and proposes a composite tensile strength model considering fiber clustering effect based on fractal theory. It provides theoretical support for the strength prediction of fiber/epoxy resin composites.

Tensile Fractal Model of Composite
Materials considering Fiber Clustering Effect

Fractal Geometric Characteristics of Cluster Fibers.
is section analyzes the geometric characteristics of one fiber agglomerate. In Figure 1, the black part is the true cross-sectional area of the fiber agglomerate. Figure 2 shows its internal void area, and Figure 3 presents the cross-sectional perimeter. To analyze the fiber beating effect based on the "fiber beating model" proposed in [19], it is necessary to know the indexes such as fiber uniform distribution coefficient and fiber beating content. In this case, this paper needs to obtain the important parameters for calculating the fiber nominal distribution coefficient and fiber beating content, including the true cross-sectional area, nominal cross-sectional area (the sum of the true cross-sectional area and the internal void area of the fiber agglomerate), and the cross-sectional perimeter of the fiber agglomerate. In practice, the cross section of fiber agglomerates is composed of irregular circular bodies arranged randomly, so it is difficult to calculate the real area of fiber agglomerates by Euclidean geometry theory. Reference [19] assumes that fibers in the agglomerate are arranged in squares and triangles, and the fiber uniform distribution coefficient and fiber cluster content are calculated by the approximate method, which is different from the reality. erefore, how to calculate the true cross-sectional area and the true crosssectional perimeter of fiber agglomerates is the key to analyzing the fiber clustering effect.
In this case, this paper analyzes the cross-sectional features of real fiber agglomerates. ey can be considered to be formed by arranging fibers according to certain physical rules (as shown in Figure 1), and their cross-sectional perimeters are formed by arcs of different lengths (as shown in Figure 3). It can be seen that the cross-sectional area of each fiber agglomerate is similar in shape, and the cross-sectional perimeter is self-similar in shape. According to the fractal theory [20][21][22], the set of shapes with self-similarity has fractal features. erefore, the cross-sectional area and crosssectional perimeter of fiber agglomerates have fractal characteristics.

Analysis of Fiber Clustering Effect Based on Fractal eory.
In this section, the relationship between fiber uniform distribution coefficient and fractal dimension of fiber agglomerates is analyzed by using relevant principles of fractal theory. e box counting dimension method [23] is used to calculate the cross-sectional perimeter C and cross-sectional area S of the fiber agglomerate, as shown in the following equations: where δ L and δ S are length and area measurement yardsticks, respectively. According to the box counting method, the circumference of the cross section of the fiber agglomerate contains is the number of δ L ), and the cross section contains M 2 (δ S )δ S (M 2 (δ S ) is the number of δ S ); D L and D S represent the fractal dimensions of the circumference and area of the section of the fiber agglomerate. According to the fiber uniform distribution coefficient defined in [19], the number (N) of fibers contained in the fiber agglomerate can be calculated according to equation (3), and the fiber uniform distribution coefficient can be calculated according to equation (4): where τ is the bonding strength between fibers and the substrate, μ indicates the friction coefficient between fibers and the substrate, d f is the diameter of fibers, l f is the length of fibers, and σ e represents the shrinkage stress of the substrate. When calculating the circumference and area of a fiber agglomerate by box counting method, it is assumed that the side length of its square box is δ 0 . en, Substitute equations (5) and (6) into (4), and the following formula can be obtained: When fibers and the substrate are determined, τ, μ, d f , l f , and σ e are all constants. For a certain fiber agglomerate, M 1 (δ L ) and M 2 (δ S ) are also constants. erefore, let en, equation (7) can be simplified as Equation (9) is the fractal expression of the uniform distribution coefficient of the i th fiber agglomerate. If the composite contains j fiber agglomerates, the fiber uniform   Advances in Civil Engineering distribution coefficient of the composite is the average value of j fiber agglomerates. As shown in equation (10), the fiber clustering content w can be calculated according to equation (11).
where V e is the volume fraction of agglomerated fiber and V f is the fiber volume fraction.

Fractal Prediction Model of Tensile Strength of Composite
Materials considering Fiber Clustering Effect. According to the composite mechanics theory [24,25], FRP is regarded as a multi-item composite elastic system with fiber reinforced resin. e strength and elastic modulus of the composite are equal to the elastic superposition of the properties of each component. e mechanical model is shown in Figure 4. erefore, the stress relationship in BFRP elastic stage can be expressed as where n � E f /E m , V m + V f � 1, E fc indicates elastic modulus of FRP, and σ fc is the stress of FRP in elastic stage; σ m represents the stress of epoxy resin matrix and σ f is the stress of fiber; V m is the volume ratio of epoxy resin. It can be seen that the tensile initial cracking strength of FRP depends on the tensile strength of epoxy resin. e larger the elastic modulus ratio of fiber to epoxy resin, the greater the initial cracking strength of FRP. Assuming that epoxy resin cracks after reaching its tensile strength under tensile load, and all load is transferred to the fiber until it is broken, then where σ u fc is FRC tensile strength, σ u f is the fiber tensile strength, and η f is the fiber effective coefficient. en, the FRC tensile strength considering fiber clustering effect is expressed as Substitute equations (11) and (12) into (14), and the FRP fractal analytical equation of ultimate tensile strength can be obtained: According to equation (15), if the fiber orientation, length, and volume fraction are fixed, the FRP tensile strength under fiber clustering effect is only related to the fractal dimension of the fiber agglomerate's perimeter and area.

Fractal Geometric Characteristics of Fiber Agglomerates in BFRP Tensile Specimens.
Taking BFRP materials prepared in [19] as an example, this article selects the BFRP samples with fiber volume fraction of 0.6%, 0.9%, 1.2%, and 1.5% to represent the cross section. After cutting and grinding, they are observed at 500 times magnification, and the arrangement of fiber agglomerate number and internal fibers under different fiber volume fractions is obtained.
For example, Figure 5 shows the gray-scale image of a fiber agglomerate when the volume fraction is 1.2%. e grayscale image is binarized (an image containing only black and white colors). at is, in the gray-scale image, there are significant differences in gray-scale values among fibers, base materials, and fiber agglomerates. e grayscale values of cluster fiber cross section are taken as the threshold. en, the regions with gray-scale values larger than the threshold are all white by default; otherwise, they are black, thus obtaining the cross-sectional arrangement and distribution of clustered fibers. As shown in Figure 6, the white part is the fiber cross section, and the black part is the gap between base materials and fiber clusters. Taking Figure 6 as an m × n matrix digital image with a height of m and a width of n, different pixel images can be obtained by changing the values of m and n. erefore, there is oneto-one correspondence between pixel positions in the image and element positions in the matrix. Since the binarized image only contains two colors (black and white), the elements of m × n matrix can be expressed by two numbers. Here, the element value "1" represents black color and "0" represents white color. e fractal dimension of the area of white part in Figure 6 is calculated by box counting method. e principle is that square boxes with different area sizes (δ Si ) are used to cover the images of the white area. At the same time, the number of boxes contained in the cross-sectional area of each fiber agglomerate M 2i (δ S ) was recorded. en, M 2 (δ S ) can be expressed as where a is a constant. en, take the logarithm of equation (16): According to equation (17), the fractal dimension of the cross-sectional area of the clustering fiber, D S , can be obtained from the linear regression analysis on the set (ln M 2 (δ S ), ln δ Si ) generated by the box counting method. It is the absolute value of the slope of the curve.
Divide the gray image of the fiber agglomerate in Figure 6 to obtain the contour map of the fiber agglomerate section, as shown in Figure 7. Likewise, the fractal dimension of contour line in Figure 7 is calculated by the box counting method. en, the fractal dimension of fiber agglomerate cross-sectional perimeter is obtained. Record the number of boxes included in the contour line image of the fiber agglomerate cross section as M 1i (δ L ); then where b is a constant. en, the set (lgM1(δ L ), lgδLi) is generated by the box counting method, and the method of calculating D S is utilized to calculate D L . e data points with self-similar characteristics in a limited scale range are fitted, and the slope of the obtained curve is the fractal dimension D L of the cross-sectional profile of the cluster fiber (cross-sectional perimeter of the cluster fiber).
Based on the above methods, the fractal dimension of fiber agglomerate cross-sectional area and perimeter on the fracture surface of BFRP samples with different fiber volume doping rates can be obtained. As shown in Figure 8, the fractal dimension of fiber agglomerate cross-sectional area is about 1.8, and the fractal dimension of cluster fiber crosssectional perimeter is about 1.1.
Take the square box area with a side length of 0.01d f as the benchmark, that is, δ S � (0.01d f ) 2 , and record the number of boxes, M 2 (δ S ). e area of the white part, S f , can be obtained through equation (2). Besides, start measuring with δ L � 0.01d f as the side length of the box and record the number of boxes M 1 (δ L ). en, the perimeter C f of different fiber agglomerate cross-sectional contours can be obtained through equation (1).     circumference of fiber agglomerate section is about 0.9-2.3 mm. Compared with the results of the approximate method, the circumference of fiber agglomerate section based on fractal theory has increased. It is due to the fact that different benchmarks can be used for actual measurement when calculating by the box counting method, and the smaller the scale, the longer the circumference of the fiber agglomerate. erefore, the perimeter calculated by fractal theory is accurate in strict sense. e cross-sectional area of the fiber agglomerate calculated by this method is about 0.03-0.17 mm 2 , which is smaller than that calculated by the approximate method. is is because the calculated area of fractal theory is based on the binarization result of the enlarged image of the fiber cluster, which reflects the true arrangement of fibers in the agglomerate. erefore, the calculated results are more objective.

Analysis of BFRP Fiber Clustering Effect Based on Fractal
eory. After substituting the data in Table 1 into equations (8) and (10), the value of fiber uniform distribution coefficient of BFRP samples based on fractal theory can be obtained. Combined with the results of [19] (calculated by approximate method), it is found that the results based on fractal theory have similar rules with the calculated results of [19]: (1) the fiber uniform distribution coefficient decreases with the increase of fiber volume fraction; (2) if the fiber volume fraction is constant, the fiber uniform distribution Table 1: Calculated values of c j and S j with fractal method.  coefficients calculated by fractal theory are all lower than those calculated by the approximate method, as shown in Figure 9. Substitute the data in Table 1 into equation (11), and the value of fiber clustering content of BFRP samples based on fractal theory can be obtained. Combining the above value with the results of [19] (calculated by approximate method), it can be found that the fiber clustering content increases with the increase of fiber volume fraction. Comparing the above value with the calculated results by the approximate method, it can be concluded that when the fiber volume fraction is constant, the calculated value of fiber clustering content by fractal theory is higher than that by the approximate method. As shown in Figure 10, when the fiber volume fraction is greater than 1.2%, the increase of the calculated values under fractal theory has slowed down.

Tensile Strength Analysis of BFRP Based on Fractal eory.
Based on the analysis in Section 3.2, this paper substitutes the calculated values of M 1 (δ L ), M 2 (δ S ), c j , and S j into equation (15) and obtains the composite tensile strength based on fractal theory. en, the results are compared with those of [19], as shown in Figure 11.
It can be seen that, with the increase of fiber volume fraction, the calculated value of composite tensile strength displays a similar trend with the experimental value. When the fiber volume fraction is constant, the calculated value of composite tensile strength under fractal theory is slightly lower than the experimental value. rough the regression analysis between the results obtained by the two methods and the experimental values, it is found that the calculated values of BFRP tensile strength based on fractal theory are in good agreement with the experimental values, and the correlation coefficient is as high as 0.941, which is closer to the experimental results than the calculated results of [19] (approximate method). It can provide theoretical basis for the design and strength prediction of BFRP composites.

Conclusion
In this paper, the fractal characteristics of agglomerated fibers were analyzed and established a prediction model of tensile strength of fiber/epoxy resin composites containing agglomerated fibers. With BFRP materials as an example, the following conclusions were reached: (1) When the fiber orientation, length, and volume fraction are fixed, the tensile strength of the composite containing agglomerated fibers is only related to the fractal dimension of the circumference and area of the inner fiber agglomerate. The fiber volume fraction V f (%) Figure 9: Fiber equidistribution coefficient as a function of fiber volume fraction.
Calculation results of approximate method [20] Calculation results of fractal method Calculation results of approximate method [20] Calculation results of fractal method Test results [20]  e calculated values of fiber uniform distribution coefficient and fiber clustering content based on fractal theory are similar to those calculated in [19]. It means the fiber uniform distribution coefficient decreases with the increase of fiber volume fraction, while the fiber clustering content increases with the increase of fiber volume fraction. e difference is that the fiber uniform distribution coefficient calculated by fractal theory is lower than that calculated by the approximate method, and the fiber clustering content calculated by fractal theory is higher than that calculated by the approximate method.
(3) e calculated value of BFRP tensile strength based on fractal theory is lower than the experimental value, but closer to the experimental value than the calculated value of the approximate method. e research conclusion can provide theoretical support for BFRP strength prediction.

Data Availability
Some or all data, models, or code generated or used during the study are available from the corresponding author upon request.

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