Quantitative Evaluation of Carbon Fiber Dispersion in Amorphous Calcium Silicate Hydrate-Based Contact-Hardening Composites

Carbon fiber dispersion has a substantial influence on the properties of amorphous calcium silicate hydrate-based contact-hardening composites. In this study, a mixture of carbon fiber and calcium silicate hydrate powder was compressed into solid composites at 40 MPa for one minute. The mechanical properties and electrical resistivity of the solid materials were measured, and the dispersion of carbon fibers was quantitatively evaluated by digital image processing technology. The Taipalu model was used to build the correlation between the electrical resistivity of the composites and the carbon fiber dispersion. The results of the electrical resistivity showed that the down threshold of carbon fiber content in the contact-hardening composites was 1.0 wt.% and the electrical resistivity was 30,000 Ω·cm. As the fiber content increased to 2.0 wt.%, the electrical resistivity dropped to 2550 Ω·cm, which was attributed to the increase in fiber dispersion uniformity in the solid composites, and the value of the fiber distribution coefficient reached a maximum value of 0.743. A subsequent decrease in the uniformity of the fiber dispersion was observed at a high fiber content. In addition, the carbon fiber content showed a slight influence on the fiber orientation in the contact-hardening composites.


Introduction
Several decades ago, researchers discovered that calcium silicate hydrate powder ((C-S-H), where C=CaO, S=SiO 2 and H=H 2 O) in an amorphous nature can be directly compressed to produce water-resistant artificial blocks with sufficient strength [1][2][3][4]. This process is similar to powder metallurgy where the strength can be derived from the bonds resulting from the deformation of the particles when they are in contact under pressure. In this case, Glukhovsky and Pu et al. [5,6] postulated "contact-hardening cementitious properties" for unstable or metastable minerals, and amorphous calcium silicate hydrate was one of the available materials for producing solid blocks with water-resistant properties. This behavior is due to the decrease in high surface energy under compression, leading to the bridging between particles in contact or even recrystallization between the surfaces of the particles [2,5].
To date, several investigations have been conducted on the contact-hardening properties of amorphous calcium silicate hydrate powder, and it has been found that solid compacts with substantial strength can be obtained in a short period at room temperature; for instance, pressing for 1 to 10 min [7,8]. No further treatment, such as standard curing, heat curing or sintering, is required. The hardening process is somehow similar to the "cold sintering" process of ceramics where interparticulate bonding is formed [9]. This provides an innovative approach to prepare composites with other types of materials, (SiO 2 ). Calcium oxide (with a free CaO content of 88.8%) was supplied by the Shandong lime plant in Chongqing. Its Blaine fineness was 410 m 2 /kg, and the average particle size was 0.8 µm. Quartz powder (with SiO 2 content of 96.5%) was milled for 90 min from quartz sand (supplied by Hunan Magnificent Vision New Materials Technology Co. Ltd.). The materials were dried at 105 • C to remove absorbed water before hydrothermal synthesis. Calcium oxide, quartz powder and water were mixed in an autoclave (20 L, Weihai Chemical Machinery Co. Ltd., Weihai, China). The Ca/Si molar ratio of the starting materials was 1:1, and the water/solid weight ratio was 5:1. It took 5 h for the temperature to increase from room temperature to 180 • C, and it was then kept constant at 180 • C for 3 h. The mixture was stirred at a rate of approximately 400 r/min. When the reaction was finished, the mixture was quickly cooled to 60 • C before being filtered. The filtered product was oven-dried at 60 • C for approximately 24 h according to the previous study [29], and the moisture content of the powder used in this study was 11.75%, compared to the powder supercritical dried at 105 • C to the constant weight. The oven-dried powder with moisture content of 11.75% was then dispersed into powders for a duration of 10 s by a mixer. Laser granulometry supplied by Haixinrui Technology Co., Ltd. (Beijing, China) was used to evaluate the particle size of the calcium silicate hydrate powders described above. The particle size distribution curve is shown in Figure 1. The average size of the particles was approximately 1 µm.

Raw Materials
The calcium silicate hydrate powder used in this study was synthesized by a dynamical hydrothermal method from a mixture of calcium oxide (CaO) and siliceous materials (SiO2). Calcium oxide (with a free CaO content of 88.8%) was supplied by the Shandong lime plant in Chongqing. Its Blaine fineness was 410 m 2 /kg, and the average particle size was 0.8 μm. Quartz powder (with SiO2 content of 96.5%) was milled for 90 min from quartz sand (supplied by Hunan Magnificent Vision New Materials Technology Co. Ltd.). The materials were dried at 105 °C to remove absorbed water before hydrothermal synthesis. Calcium oxide, quartz powder and water were mixed in an autoclave (20 L, Weihai Chemical Machinery Co. Ltd., Weihai, China). The Ca/Si molar ratio of the starting materials was 1:1, and the water/solid weight ratio was 5:1. It took 5 h for the temperature to increase from room temperature to 180 °C, and it was then kept constant at 180 °C for 3 h. The mixture was stirred at a rate of approximately 400 r/min. When the reaction was finished, the mixture was quickly cooled to 60 °C before being filtered. The filtered product was oven-dried at 60 °C for approximately 24 h according to the previous study [29], and the moisture content of the powder used in this study was 11.75%, compared to the powder supercritical dried at 105 °C to the constant weight. The oven-dried powder with moisture content of 11.75% was then dispersed into powders for a duration of 10 s by a mixer. Laser granulometry supplied by Haixinrui Technology Co., Ltd. (Beijing, China) was used to evaluate the particle size of the calcium silicate hydrate powders described above. The particle size distribution curve is shown in Figure 1. The average size of the particles was approximately 1 μm. The mineral phases of the synthesized calcium silicate hydrate powder were determined by using X-ray diffraction (XRD, RigakuD/max-1200, Cu Kα, Osaka, Japan) with an output power of 3 kW. A step size of 0.02°, scanning rate of 2°/min and scanning range from 15° to 65° were applied. The XRD pattern of the powders is shown in Figure 2. The diffused peaks can be observed at d-spacings of approximately 1.4 nm (2θ of 6.5°), 0.303 nm (2θ of 29.4°) and 0.182 nm (2θ of 51.2°), which are typical for the amorphous calcium silicate hydrate [30]. In addition, crystalline phases, including α-SiO2, α-C2SH and calcite, were also detected in the powder. α-SiO2 with d-spacings of 0.425, 0.334, 0.1819 and 0.1541 nm was derived from the quartz powder. α-C2SH is a crystalline calcium silicate hydrate with d-spacings of 0.421, 0.324, 0.277, 0.250, 0.218, 0.179 nm [31]. The presence of calcite was due to the carbonation of calcium silicate hydrate during drying and storage process. The mineral phases of the synthesized calcium silicate hydrate powder were determined by using X-ray diffraction (XRD, RigakuD/max-1200, Cu Kα, Osaka, Japan) with an output power of 3 kW. A step size of 0.02 • , scanning rate of 2 • /min and scanning range from 15 • to 65 • were applied. The XRD pattern of the powders is shown in Figure 2. The diffused peaks can be observed at d-spacings of approximately 1.4 nm (2θ of 6.5 • ), 0.303 nm (2θ of 29.4 • ) and 0.182 nm (2θ of 51.2 • ), which are typical for the amorphous calcium silicate hydrate [30]. In addition, crystalline phases, including α-SiO 2 , α-C 2 SH and calcite, were also detected in the powder. α-SiO 2 with d-spacings of 0.425, 0.334, 0.1819 and 0.1541 nm was derived from the quartz powder. α-C 2 SH is a crystalline calcium silicate hydrate with d-spacings of 0.421, 0.324, 0.277, 0.250, 0.218, 0.179 nm [31]. The presence of calcite was due to the carbonation of calcium silicate hydrate during drying and storage process.
The nanostructure of the synthetic powder was characterized by transmission electron microscopy (TEM, Libra 200, Zeiss, Germany) with a maximum resolution of 1.4 nm. The particles consisted of globular particles with sizes of approximately 20 to 50 nm (Figure 3a), and the selected area electron diffraction (SAED) pattern in Figure 3b further demonstrated the presence of amorphous calcium silicate hydrate due to the presence of a faint halo around it.  The carbon fibers in the current study were supplied by Zhongfu Shenying Carbon Fiber Co., Ltd. The properties are shown in Table 1.

Preparation of Contact-Hardening Composites
In this study, a mixer (O-500Y, Oubo high speed pulverizer, Jinhua, China) was used to mix calcium silicate hydrate powders and carbon fibers at different weight ratios for approximately 10 s to produce a homogeneous mixture before being poured into a 40 × 100 × 160 mm 3 steel mold. A certain amount of the mixture was compressed by increasing the pressure on both sides with a compression pressure of 40 MPa at a rate of 0.4 MPa/s.   The carbon fibers in the current study were supplied by Zhongfu Shenying Carbon Fiber Co., Ltd. The properties are shown in Table 1.

Preparation of Contact-Hardening Composites
In this study, a mixer (O-500Y, Oubo high speed pulverizer, Jinhua, China) was used to mix calcium silicate hydrate powders and carbon fibers at different weight ratios for approximately 10 s to produce a homogeneous mixture before being poured into a 40 × 100 × 160 mm 3 steel mold. A certain amount of the mixture was compressed by increasing the pressure on both sides with a compression pressure of 40 MPa at a rate of 0.4 MPa/s. The carbon fibers in the current study were supplied by Zhongfu Shenying Carbon Fiber Co., Ltd. The properties are shown in Table 1.

Preparation of Contact-Hardening Composites
In this study, a mixer (O-500Y, Oubo high speed pulverizer, Jinhua, China) was used to mix calcium silicate hydrate powders and carbon fibers at different weight ratios for approximately 10 s to produce a homogeneous mixture before being poured into a 40 × 100 × 160 mm 3 steel mold. A certain amount of the mixture was compressed by increasing the pressure on both sides with a compression pressure of 40 MPa at a rate of 0.4 MPa/s. The pressure was maintained at 40 MPa for 1 min. The contacthardening specimen with size of 40 × 40 × 160 mm 3 prepared from the mixture of calcium silicate hydrate and carbon fiber was removed from the mold immediately after being compacted. The procedure of the contact-hardening composite preparation is shown in the flow chart below (Figure 4). More details of the preparation of the hardening specimens were described in [32]. The influence of fiber content (ranging from 0.5 wt.% to 4.0 wt.%) on the properties of the contact-hardening composite was investigated.
The pressure was maintained at 40 MPa for 1 min. The contact-hardening specimen with size of 40 × 40 × 160 mm 3 prepared from the mixture of calcium silicate hydrate and carbon fiber was removed from the mold immediately after being compacted. The procedure of the contact-hardening composite preparation is shown in the flow chart below (Figure 4). More details of the preparation of the hardening specimens were described in [32]. The influence of fiber content (ranging from 0.5 wt.% to 4.0 wt.%) on the properties of the contact-hardening composite was investigated.

Characterization of the Contact-Hardening Composites
The flexural strength, compressive strength, and bulk density of the specimens were measured according to the Test Methods for Concrete Blocks and Bricks (GB/T 4111-2013).
The electrical conductivity was tested according to the two-point electrical conductance test principle by using the inductance capacitance resistance (LCR) method [33]. The specimen was placed between two metal plates. Wet cotton cloths were placed between the specimen and metal plates to ensure that they were in full contact. When the test was carried out, an LCR bridge was used to alternate the current with a frequency setting of 1 Hz to reduce the influence of the direct current and low frequency on the polarization of the specimen. The test was conducted after 2 min of contact, and steady-state of the current was achieved. The electrical resistivity (Ω·cm) was calculated by using Equation (1).

=
(1) where R (Ω) is the electrical resistance of the specimen, S (cm 2 ) is the cross-section area of the specimen (cm 2 ) and l (cm) is the length of the specimen.

Image Method for Fiber Distribution and Orientation
To improve the accuracy and validity of the calculation results of carbon fiber distribution and orientation, three specimens were prepared in each mixing proportion, and each specimen was cut into 5 pieces. Initially, the surface of the cross-section was polished roughly on an electrical grinding machine by using # 300 coarse sandpaper. The specimens were then polished by using finer sandpapers, i.e., # 600 paper, # 1000 paper and then # 1500 fine sandpaper, to produce an adequately smooth surface with minimum damage. A highresolution charge-coupled device camera (CCD, EOS 80D) supplied by Canon was used to take photographs of the polished surface.
Digital image processing technology (DIPT) was used to quantify the distribution and orientation of the carbon fibers in this study. The whole process is shown in Figure 5. The original image was converted into a grey image. It was then converted into a binary image by referring to the maximum interclass variance method proposed by Otsu [34],

Characterization of the Contact-Hardening Composites
The flexural strength, compressive strength, and bulk density of the specimens were measured according to the Test Methods for Concrete Blocks and Bricks (GB/T 4111-2013).
The electrical conductivity was tested according to the two-point electrical conductance test principle by using the inductance capacitance resistance (LCR) method [33]. The specimen was placed between two metal plates. Wet cotton cloths were placed between the specimen and metal plates to ensure that they were in full contact. When the test was carried out, an LCR bridge was used to alternate the current with a frequency setting of 1 Hz to reduce the influence of the direct current and low frequency on the polarization of the specimen. The test was conducted after 2 min of contact, and steady-state of the current was achieved. The electrical resistivity ρ (Ω·cm) was calculated by using Equation (1).
where R (Ω) is the electrical resistance of the specimen, S (cm 2 ) is the cross-section area of the specimen (cm 2 ) and l (cm) is the length of the specimen.

Image Method for Fiber Distribution and Orientation
To improve the accuracy and validity of the calculation results of carbon fiber distribution and orientation, three specimens were prepared in each mixing proportion, and each specimen was cut into 5 pieces. Initially, the surface of the cross-section was polished roughly on an electrical grinding machine by using # 300 coarse sandpaper. The specimens were then polished by using finer sandpapers, i.e., # 600 paper, # 1000 paper and then # 1500 fine sandpaper, to produce an adequately smooth surface with minimum damage. A high-resolution charge-coupled device camera (CCD, EOS 80D) supplied by Canon was used to take photographs of the polished surface.
Digital image processing technology (DIPT) was used to quantify the distribution and orientation of the carbon fibers in this study. The whole process is shown in Figure 5. The original image was converted into a grey image. It was then converted into a binary image by referring to the maximum interclass variance method proposed by Otsu [34], and thus, the geometric morphological information about carbon fibers was obtained, as shown in Figure 5 In terms of the distribution of the carbon fibers, the centroids of each carbon fiber were extracted by using the "Centroid Function" in the MATLAB 2018a software. The Voronoi diagram of the distribution of the centroids (see Figure 5-5 # ) was generated based on the Delaunay triangulation method (see Figure 5-4 # ). The specific description of the Delaunay triangulation method and Voronoi diagram can be found in Bosman's study [35].
With regard to the orientation of the carbon fibers, the Sobel operator [36] was used to verify the edge profile of the carbon fiber in a binary image. The edge contour of the fibers was further quantified by using the "External Area Rectangle" in MATLAB. The orientation of the carbon fibers at the cross-section was presented in two ways. When cutting along the cross-section, fibers showed elliptical features and when cut along the longitudinal section, the fibers showed rectangular features, as shown in Figure 5-8 # . The orientation of the carbon fibers was calculated based on the study conducted by Song [37], and the colormap in MATLAB was used to characterize the fiber orientation nephogram. and thus, the geometric morphological information about carbon fibers was obtained, as shown in Figure 5 In terms of the distribution of the carbon fibers, the centroids of each carbon fiber were extracted by using the "Centroid Function" in the MATLAB 2018a software. The Voronoi diagram of the distribution of the centroids (see Figure 5-5 # ) was generated based on the Delaunay triangulation method (see Figure 5-4 # ). The specific description of the Delaunay triangulation method and Voronoi diagram can be found in Bosman's study [35].
With regard to the orientation of the carbon fibers, the Sobel operator [36] was used to verify the edge profile of the carbon fiber in a binary image. The edge contour of the fibers was further quantified by using the "External Area Rectangle" in MATLAB. The orientation of the carbon fibers at the cross-section was presented in two ways. When cutting along the cross-section, fibers showed elliptical features and when cut along the longitudinal section, the fibers showed rectangular features, as shown in Figure 5-8 # . The orientation of the carbon fibers was calculated based on the study conducted by Song [37], and the colormap in MATLAB was used to characterize the fiber orientation nephogram.

Mechanical Properties of Contact-Hardening Composites
The compressive strength, flexural strength, and bulk density of the contact-hardening composites are shown in Figure 6.

Mechanical Properties of Contact-Hardening Composites
The compressive strength, flexural strength, and bulk density of the contact-hardening composites are shown in Figure 6.
From Figure 6a, it can be seen that as the carbon fiber content increased from 0 to 4.0 wt.%, the apparent density gradually increased since the density of the carbon fibers was higher than that of the amorphous calcium silicate hydrate-based contact-hardening specimen prepared at 40 MPa. The structure became more compact with an increase in fiber content. An opposite tendency was observed in the compressive strength of the contact-hardening composites, in which it decreased slightly from 19.4 MPa to 18.0 MPa as the fiber content increased to 4.0 wt.% (Figure 6b). This is probably because calcium silicate hydrate powders were not able to cover the surface of the carbon fibers, and it is easy to form defects at the interface between fiber and C-S-H. With increasing carbon fiber content, the number of defects increased further. During compression, the microcracks expanded more rapidly and therefore caused failure to occur earlier. However, the flexural strength was greatly improved by increasing the carbon fiber content. When the carbon fiber contents were 0.5, 1.0, 2.0, 3.0, and 4.0 wt.%, the flexural strength increased by 141, 216, 367, 408, and 558%, respectively, compared with the specimen that contained no carbon fibers. This can be explained by the following reason. The flexural strength is the ability of a material to resist bending without breaking. The effective bearing area of the carbon fibers increased with increasing fiber numbers, improving the brittleness of the matrix, as described by Jia et al. [38]. From Figure 6a, it can be seen that as the carbon fiber content increased from 0 to 4.0 wt.%, the apparent density gradually increased since the density of the carbon fibers was higher than that of the amorphous calcium silicate hydrate-based contact-hardening specimen prepared at 40 MPa. The structure became more compact with an increase in fiber content. An opposite tendency was observed in the compressive strength of the contacthardening composites, in which it decreased slightly from 19.4 MPa to 18.0 MPa as the fiber content increased to 4.0 wt.% (Figure 6b). This is probably because calcium silicate hydrate powders were not able to cover the surface of the carbon fibers, and it is easy to form defects at the interface between fiber and C-S-H. With increasing carbon fiber content, the number of defects increased further. During compression, the microcracks expanded more rapidly and therefore caused failure to occur earlier. However, the flexural strength was greatly improved by increasing the carbon fiber content. When the carbon fiber contents were 0.5, 1.0, 2.0, 3.0, and 4.0 wt.%, the flexural strength increased by 141, 216, 367, 408, and 558%, respectively, compared with the specimen that contained no carbon fibers. This can be explained by the following reason. The flexural strength is the ability of a material to resist bending without breaking. The effective bearing area of the carbon fibers increased with increasing fiber numbers, improving the brittleness of the matrix, as described by Jia et al. [38].

Electrical Resistivity of Contact-Hardening Composites
Resistivity is an important physical parameter for measuring the electrical conductivity of engineering materials, especially fiber-reinforced cement-based materials. The influence of the carbon fiber content on the electrical resistivity of contact-hardening composites is shown in Figure 7.

Electrical Resistivity of Contact-Hardening Composites
Resistivity is an important physical parameter for measuring the electrical conductivity of engineering materials, especially fiber-reinforced cement-based materials. The influence of the carbon fiber content on the electrical resistivity of contact-hardening composites is shown in Figure 7. From Figure 7, it can be seen that the electrical resistivity of the contact-hardening composites decreased from 42,000 Ω·cm to 30,000 Ω·cm as the carbon fiber content increased from 0 to 1.0 wt.%, but a more dramatic decrease occurred at fiber contents between 1.0 wt.% and 2.0 wt.%. This was due to the percolation phenomenon that when the conductive phase content in the nonconductive matrix exceeded a certain threshold, the conductive phases would connect with each other and lead to a sudden change in electrical resistivity [39]. The mass fraction of conductive materials where the electrical resistivity dramatically decreased with increasing carbon fiber content is called the down threshold, which was 1 wt.% for the material studied with an electrical resistivity of 30,000 Ω·cm. This was followed by a dramatic decrease in electrical resistivity of 2550 Ω·cm at a fiber content of 2.0 wt.%, which was considered the upper threshold in this study.
In addition, when the carbon fiber content was lower than 1.0 wt.%, the value of electrical resistivity was relatively high due to the comparatively large distance between the carbon fibers. In this case, the carbon fibers failed to connect with each other to form a network. Therefore, a good electrical path was unlikely formed to provide considerable conductance. As the carbon fiber content increased, the distance between carbon fibers decreased, the number of connections increased and a network structure formed. Electrical electricity was conducted by the tunneling effect of the carbon fibers, as demonstrated by Zhang et al. [40]. In addition, the electrical resistivity decreased dramatically at fiber  Figure 7, it can be seen that the electrical resistivity of the contact-hardening composites decreased from 42,000 Ω·cm to 30,000 Ω·cm as the carbon fiber content increased from 0 to 1.0 wt.%, but a more dramatic decrease occurred at fiber contents between 1.0 wt.% and 2.0 wt.%. This was due to the percolation phenomenon that when the conductive phase content in the nonconductive matrix exceeded a certain threshold, the conductive phases would connect with each other and lead to a sudden change in electrical resistivity [39]. The mass fraction of conductive materials where the electrical resistivity dramatically decreased with increasing carbon fiber content is called the down threshold, which was 1 wt.% for the material studied with an electrical resistivity of 30,000 Ω·cm. This was followed by a dramatic decrease in electrical resistivity of 2550 Ω·cm at a fiber content of 2.0 wt.%, which was considered the upper threshold in this study.
In addition, when the carbon fiber content was lower than 1.0 wt.%, the value of electrical resistivity was relatively high due to the comparatively large distance between the carbon fibers. In this case, the carbon fibers failed to connect with each other to form a network. Therefore, a good electrical path was unlikely formed to provide considerable conductance. As the carbon fiber content increased, the distance between carbon fibers decreased, the number of connections increased and a network structure formed. Electrical electricity was conducted by the tunneling effect of the carbon fibers, as demonstrated by Zhang et al. [40]. In addition, the electrical resistivity decreased dramatically at fiber contents of between 1.0 wt.% and 2.0 wt.%. When the carbon fiber content was higher than 2.0 wt.%, the electrical resistivity started to slow down. This was because the conductive network of carbon fibers was already formed and an excess of carbon fiber content would probably increase the number of fiber clusters, leading to decreasing uniformity of the fiber distribution in the matrix. As a result, only a minor change was observed in the electrical resistivity. More information about the fiber distribution is discussed in the following section.
According to the above results, an optimum dispersion of carbon fibers in contacthardening composites appeared at a carbon fiber content of 2.0 wt.%.

Analysis of the Carbon Fiber Distribution
Carbon fibers in the contact-hardening composite can connect with each other and become electrical conductors. In addition, the distribution of the fibers was also assumed to play an important role in inhibiting the propagation of cracks and improving the flexural strength of the specimen. Figure 8 shows the original high-resolution CCD images of the cross-section of the contact-hardening composites with different carbon fiber contents, and Figure 9 presents the binary images from Figure 8 after processing by using Otsu's method [34].     Voronoi diagrams with different carbon fiber contents are shown in Figure 10. Each picture was divided into nine subsections, which are shown in Figure 11. The number of fibers in each subsection was counted. The formula of the fiber distribution coefficient (α f ) proposed by Han et al. [41] was used to quantitatively evaluate the uniformity of the fiber distribution, as described in Equation (2).
where n = 9, which reflects the number of segmented portions, X i is the number of carbon fibers positioned in the ith area, and X average is the average number of carbon fibers in all areas. α f is the fiber distribution coefficient. α f = 1 indicates that fibers are completely homogeneously distributed, and α f = 0 indicates that the agglomeration of fibers. j is the number of the cross-section. For each mixture, a section that was close to the average carbon fiber distribution coefficient was selected from 15 sections of each mixing proportion, and the section comprehensively reflected the dispersion characteristics of fibers with different mass fractions in the composites. The statistic of the number of carbon fibers for the different additions are shown in Table 2. The calculation results are shown in Figure 12. It can be seen that with the increase in carbon fiber content, the values of α f first increased and then decreased. The values were 0.613, 0.680, 0.743, 0.557 and 0.479, corresponding to the fiber content of 0.5, 1.0, 2.0, 3.0 and 4.0 wt.%, respectively. At a fiber content of 0.5 wt.%, the number of fibers was not enough to be dispersed in the whole area of the matrix. The initial increase in the value of α f with an increase in fiber content was ascribed to the decrease in the discreteness degree of the fibers in each subarea as the fibers gradually came into contact. As a result, the distribution of the fibers became more homogeneous in the matrix as the fiber content increased up to 2.0 wt.%, and α f increased to 0.743. The subsequent decrease in the α f value when the fiber content was beyond 2 wt.%, demonstrating that a cluster of fibers occurred because of the excessive amount of fibers. Such a result was consistent with Kang's research [28]. This demonstrates that a carbon fiber content of 2 wt.% was substantial in the amorphous calcium silicate hydrate-based contact-hardening composite, which also provided an optimum value of electrical resistivity and rapid growth in flexural strength.
proposed by Han et al. [41] was used to quantitatively evaluate the uniformity of the fiber distribution, as described in Equation (2).
where n = 9, which reflects the number of segmented portions, Xi is the number of carbon fibers positioned in the i th area, and Xaverage is the average number of carbon fibers in all areas. f is the fiber distribution coefficient. f = 1 indicates that fibers are completely homogeneously distributed, and f = 0 indicates that the agglomeration of fibers. is the number of the cross-section.  For each mixture, a section that was close to the average carbon fiber distribution coefficient was selected from 15 sections of each mixing proportion, and the section com prehensively reflected the dispersion characteristics of fibers with different mass fraction in the composites. The statistic of the number of carbon fibers for the different addition

Analysis of the Fiber Orientation
The orientation of the carbon fibers provided a clear and easy way to characterize the carrying capacity. When conducting an analysis of the orientation, the image processing procedures described in Section 2.4 were used (see Figure 5-7 # , Figure 5-8 # , Figure 5-9 # ). The "Sobel operator" of the MATLAB 2018-a software, which is a typical edge detection operator based on the first-order derivative of the image brightness function, was used in the current study. The algorithm is more applicable than the Prewitt operator and Roberts operator because it introduces local average operator [42]. The extraction figure is shown in Figure 13. Then, the external area rectangle (EAR) was used to measure the fiber orientation, which is shown in Figure 14. When using the EAR, each fiber was wrapped by drawing the external area rectangle.
The inclined angle coefficient θ of these two types of sections of carbon fibers can be calculated by Equation (3) [37]. Finally, the fiber orientation angle nephogram was used to illustrate the orientation of the carbon fibers, which is shown in Figure 5-9#.

= arctan
where lr is the length of the external area rectangle and d is the width of the external area rectangle.θ tending to 90 o indicates that the carbon fibers are aligned vertically, showing a weak carrying capacity, and θ = 0° means that the carbon fibers are aligned laterally, showing a strong carrying capacity.

Analysis of the Fiber Orientation
The orientation of the carbon fibers provided a clear and easy way to characterize the carrying capacity. When conducting an analysis of the orientation, the image processing procedures described in Section 2.4 were used (see Figure 5-7 # -9 # ). The "Sobel operator" of the MATLAB 2018-a software, which is a typical edge detection operator based on the first-order derivative of the image brightness function, was used in the current study. The algorithm is more applicable than the Prewitt operator and Roberts operator because it introduces local average operator [42]. The extraction figure is shown in Figure 13. Then, the external area rectangle (EAR) was used to measure the fiber orientation, which is shown in Figure 14. When using the EAR, each fiber was wrapped by drawing the external area rectangle.
The inclined angle coefficient θ of these two types of sections of carbon fibers can be calculated by Equation (3) [37]. Finally, the fiber orientation angle nephogram was used to illustrate the orientation of the carbon fibers, which is shown in Figure 5-9#.
where l r is the length of the external area rectangle and d is the width of the external area rectangle. θ tending to 90 • indicates that the carbon fibers are aligned vertically, showing a weak carrying capacity, and θ = 0 • means that the carbon fibers are aligned laterally, showing a strong carrying capacity. The carbon fiber orientation angle was further divided into every five degrees to calculate the ratio of the number of carbon fibers at that angle range to the total number of fibers in each section. The histograms of the carbon fiber orientation angles are shown in Figure 15.     The carbon fiber orientation angle was further divided into every five degrees to calculate the ratio of the number of carbon fibers at that angle range to the total number of fibers in each section. The histograms of the carbon fiber orientation angles are shown in Figure 15.  and 87.8% of the fibers were in the angle range of between 25 • and 65 • , which showed no obvious difference. This result indicates that fibers mainly oriented in the angle range of 40 • to 50 • . The fiber content has little influence on the fiber orientation of contact-hardening composites. In traditional cement paste, the fiber orientation increases gradually with the increase in viscosity and yield stress [43]. However, in the process of powder compression, the orientation of carbon fiber is only affected by the compression pressure, so the increase in fiber content has little influence on the orientation. The figure suggests that with increasing fiber content, the frequency of connection between fibers in each subarea increased, although no substantial changes in the fiber orientation angles were observed in the fiber distribution ( Figure 16). However, an excessive fiber content (e.g., 4.0 wt.% in Figure 16e) increases the possibility of fibers connecting, winding, and balling during the contact-hardening process, and the orientation angle becomes increasingly random.
From Figure 15, it can be seen that the carbon fibers in amorphous calcium silicate hydrate-based contact-hardening composites were randomly distributed. The fiber orientation angles in each group mainly ranged from 45° to 50°. Approximately 30.7, 30.1, 31.6, 29.7 and 31.1% of the fibers were in this angle range for the composites containing carbon fiber contents of 0.5, 1.0, 2.0, 3.0 and 4.0 wt.%, respectively. Meanwhile, 89.6, 88.6, 89.4, 87.7 and 87.8% of the fibers were in the angle range of between 25° and 65°, which showed no obvious difference. This result indicates that fibers mainly oriented in the angle range of 40° to 50°. The fiber content has little influence on the fiber orientation of contact-hardening composites. In traditional cement paste, the fiber orientation increases gradually with the increase in viscosity and yield stress [43]. However, in the process of powder compression, the orientation of carbon fiber is only affected by the compression pressure, so the increase in fiber content has little influence on the orientation. The figure suggests that with increasing fiber content, the frequency of connection between fibers in each subarea increased, although no substantial changes in the fiber orientation angles were observed in the fiber distribution ( Figure 16). However, an excessive fiber content (e.g., 4.0 wt.% in Figure 16e) increases the possibility of fibers connecting, winding, and balling during the contact-hardening process, and the orientation angle becomes increasingly random.

The Taipalus Model
The electrical resistivity reflects the electrical conductance properties of contact-hardening composites in different fiber contents, but it fails to intuitively reflect the real characteristics of carbon fiber dispersion. It is necessary to establish the relationship model between the electrical resistivity and dispersion. Based on the thermal and electrical conductive model of particles' filler composites [44], Weber and Kamal [45] proposed the fiber contact model to calculate the volume resistivity of FRC (fiber-reinforced composites).

The Taipalus Model
The electrical resistivity reflects the electrical conductance properties of contacthardening composites in different fiber contents, but it fails to intuitively reflect the real characteristics of carbon fiber dispersion. It is necessary to establish the relationship model between the electrical resistivity and dispersion. Based on the thermal and electrical conductive model of particles' filler composites [44], Weber and Kamal [45] proposed the fiber contact model to calculate the volume resistivity of FRC (fiber-reinforced composites). Taipalus et al. [46] improved Weber and Kamal's model by introducing the electrical resistivity of the matrix. The modified classical model has been widely used to study the electrical effect of carbon fiber on polypropylene [47], rod-filled composites [48], and hexagonal-boron nitride/methyl-vinyl-silicone rubber (h-BN/MVQ) [49].These theoretical values of resistivity are somehow dependent on the experimental results. However, there is a lack of studies on the classical model that has been conducted on contact-hardening composites. Therefore, it is necessary to investigate the carbon fiber dispersion of the novel constituent system in the framework by applying the Taipalus model. The electrical resistivity of the composite can be defined by Equation (4).
where ρ c in Ω·cm is the electrical resistivity of the composites, ρ m in Ω·cm is the electrical resistivity of the matrix, and d c is the diameter of the contact circle, the value of which is 40 µm according to Ref. [50]. d 0 is the diameter of the carbon fiber with a value of 40 µm, and l 0 is the average carbon fiber length, which equals 10 mm. σ f is the electrical resistivity of the carbon fibers, and the value is 290 Ω·cm. θ, V p and X are explained below. θ denotes the orientation angle of the carbon fibers, which can be calculated by Equation (5).
where ρ(θ) is the probability function of the carbon fiber orientation, n is the step length of the differential equation, and n = 9 in this study. θ i is the fiber orientation angle of the i th subsection, which is defined as the average of the interval between the minimum and maximum angles. ρ(θ i ) is the probability value of θ i . The parameter V p denotes the volume fraction of carbon fibers involved in forming conductive paths, which can be expressed by Equation (6) [46].
where V f is the mass fraction of carbon fibers in the composites. V per is the threshold mass fraction corresponding to the lower threshold, which means that there are not enough carbon fibers to bridge each other to form pathways, and V per was 1.0 wt. % in this study. V t is the mass fraction corresponding to the upper threshold, which means that a sufficient amount of carbon fibers can bridge into pathways so that the electrical resistivity of the composite does not increase obviously, and V t was 2.0 wt.% in this research.
The parameter X denotes the contact factor of the carbon fibers per unit volume. Weber and Kamal [45] provided an empirical formula to establish X and the contact numbers of the fibers. However, due to the difficulty of determining the maximum numbers of contacts, this formula cannot be universally applied to other complicated systems. Therefore, in this study, a modified formula is provided in Equation (7).
where α f is the fiber distribution coefficient. When the carbon fibers were evenly distributed, α f was 1. In this case, each fiber can be regarded as being completely covered by calcium silicate hydrate powder without direct contact, and as a result, X = 0. α f was zero when the carbon fibers were completely agglomerated. They may be considered to be in contact with each other. By substituting Equations (5)-(7) into the Equation (4), the empirical formula of the carbon fiber distribution can be obtained, as shown in Equation (8). The residual can be calculated based on Equation (9) Residual = ρ model − ρ experiment (9) where ρ model is the predictive electrical resistivity, and ρ experiment is the experimental electrical resistivity.

Model Validation and Analysis
To validate the Taipalu model, all the parameters involved in Equations (4)-(7) as well as R-square (R 2 ), were solved by the software MATLAB 2018-a and listed in Table 3. From Table 3, the R-square value is 0.976, which shows that there is obvious agreement between the predicted value from the model and the experimental results [51]. The parameter V p indicates that the electrical resistivity of the contact-hardening composites satisfies the percolation-threshold theory, as demonstrated in [39]. Moreover, θ has no substantial influence on increasing carbon fiber content due to the fibers bridging each other and the phenomenon of preferred orientation, which is completely different from the study of carbon-filled polymer composites reported in [45]. Moreover, when the carbon fiber content was lower than 2.0 wt.%, the distribution of carbon fibers had no direct effect on the conductivity of the composites due to the lack of adequate electrical paths in the composites. When the carbon fiber content was higher than 2.0 wt.%, the cos 2 θ was between 0.486 and 0.508, and it could be assumed that all the carbon fiber groups had a cos 2 θ value equal to 0.50. Substituting the parameters in Table 3 into Equation (8), the distribution coefficient α f for the carbon fiber content of 2.0, 3.0 and 4.0 wt.% shown in Figure 14 can be described by Equation (10). Figure 17 shows the relationship between α f and the electrical resistivity. A strong relationship between α f and ρ c can be observed in the figure, and the electrical resistivity of the contract-hardening composites decreased with an increased carbon fiber content at a content higher than 2.0 wt.%, although a poorer fiber distribution appeared. This may be because the high carbon fiber content forms more conductive paths than low fiber dosages.

Correlation between Mechanical Properties and Dispersion Parameters
In Figure 6b, it shows that the compressive strength decreased with carbon fiber content. For the compressive strength, it is greatly affected by the interface between calcium silicate hydrate powder and carbon fibers. At a higher carbon fiber content, considering the complexity of carbon fiber agglomeration and the wall effect of powder, it is not easy to establish a relationship between fiber dispersion (including fiber distribution and orientation) and compressive strength.

Correlation between Mechanical Properties and Dispersion Parameters
In Figure 6b, it shows that the compressive strength decreased with carbon fiber content. For the compressive strength, it is greatly affected by the interface between calcium silicate hydrate powder and carbon fibers. At a higher carbon fiber content, considering the complexity of carbon fiber agglomeration and the wall effect of powder, it is not easy to establish a relationship between fiber dispersion (including fiber distribution and orientation) and compressive strength.
With respect to the flexural strength, the value increased with carbon fiber content ( Figure 6b). However, it is shown in Figure 12 that the carbon fiber distribution coefficient of the contact-hardening composites increased first with increasing fiber content to 2 wt.%, and then decreased with an excessive fiber content. Meanwhile, the orientation angle ranges were not affected by the carbon fiber content as shown in Figure 15, due to the constant compression pressure applied on the specimen for one minute. Thus, there should be no direct relationship between flexural strength and carbon fiber distribution and orientation. For further study on the fiber distribution, however, Bosman et al. [35] proposed an indirect relationship between the flexural strength and fiber distribution. Gao et al. [17] investigated the influence of the spatial distribution coefficient on the flexural strength and toughness of carbon fiber-reinforced cement-based composites by the grey entropy correlation analysis. In this study, the product of the total number of the fibers in the composites and the distribution coefficient was termed as the effective carbon fiber numbers ( ), as shown in Equation (11): where Ne is the number of effective carbon fibers, αf is the distribution coefficient of carbon fibers, Xtotal is the total number of carbon fibers. The relationship between the flexural strength and the effective carbon fiber numbers is shown in Figure 18. In general, a linear relationship with R-square of 0.8872 was observed, and the relationship between flexural strength and the carbon fiber distribution coefficient can be further estimated by the following equation: = 573 + 135,155 10,000 where is the flexural strength of the carbon fiber-reinforced contact-hardening composites.
Despite the fact that the equation indicates that the distribution coefficient of carbon fiber affects the flexural properties of the matrix, it should also be noted that at a higher With respect to the flexural strength, the value increased with carbon fiber content ( Figure 6b). However, it is shown in Figure 12 that the carbon fiber distribution coefficient of the contact-hardening composites increased first with increasing fiber content to 2 wt.%, and then decreased with an excessive fiber content. Meanwhile, the orientation angle ranges were not affected by the carbon fiber content as shown in Figure 15, due to the constant compression pressure applied on the specimen for one minute. Thus, there should be no direct relationship between flexural strength and carbon fiber distribution and orientation. For further study on the fiber distribution, however, Bosman et al. [35] proposed an indirect relationship between the flexural strength and fiber distribution. Gao et al. [17] investigated the influence of the spatial distribution coefficient on the flexural strength and toughness of carbon fiber-reinforced cement-based composites by the grey entropy correlation analysis. In this study, the product of the total number of the fibers in the composites and the distribution coefficient was termed as the effective carbon fiber numbers (N e ), as shown in Equation (11): N e = α f X total (11) where N e is the number of effective carbon fibers, α f is the distribution coefficient of carbon fibers, X total is the total number of carbon fibers. The relationship between the flexural strength and the effective carbon fiber numbers is shown in Figure 18. In general, a linear relationship with R-square of 0.8872 was observed, and the relationship between flexural strength and the carbon fiber distribution coefficient can be further estimated by the following equation: where σ c f t is the flexural strength of the carbon fiber-reinforced contact-hardening composites. Despite the fact that the equation indicates that the distribution coefficient of carbon fiber affects the flexural properties of the matrix, it should also be noted that at a higher carbon fiber content, the effective carbon fiber numbers started with a slight change, but a remarked increase in flexural strength was still observed. This is due to the stack of the fiber observed by optical microscopy at a higher content. In future work, X-ray computed tomography (X-CT) will be used to characterize the spatial fiber in the composites to provide a three-dimensional fiber distribution prediction. carbon fiber content, the effective carbon fiber numbers started with a slight change, but a remarked increase in flexural strength was still observed. This is due to the stack of the fiber observed by optical microscopy at a higher content. In future work, X-ray computed tomography (X-CT) will be used to characterize the spatial fiber in the composites to provide a three-dimensional fiber distribution prediction.

Conclusions
(1) The carbon fiber content has a slight influence on the compressive strength and bulk density of amorphous calcium silicate hydrate-based contact-hardening composites, but it can greatly improve the flexural strength of the composites. When the carbon fiber content is 0.5, 1.0, 2.0, 3.0, and 4.0 wt.%, the flexural strength increases by 141, 216, 367, 408, and 558%, respectively, compared to the contact-hardening composite without fibers.
(2) The electrical resistivity of the specimens decreased with increasing carbon fiber content. A lower threshold for contact-hardening composites with carbon fibers appeared at a fiber content of 1.0 wt.%, and the electrical resistivity was 30,000 Ω·cm. A higher threshold was present at a fiber content of 2.0 wt.% with an electrical resistivity of 2550 Ω·cm. It was found that a carbon fiber content of 2.0 wt.% in the contact-hardening composite provided optimum uniformity of the carbon fiber dispersion in the contact-hardening composites.
(3) Direct analysis of the fiber distribution and orientation by using the imaging method suggests that by increasing the fiber content up to 2.0 wt.%, the value of the fiber distribution coefficient ( ) increased dramatically to 0.743, representing a more homogeneous distribution of fibers in the composite. When the fiber content was beyond 2 wt.%, the value of decreased with increasing fiber content, which is ascribed to the fiber cluster in the matrix. Approximately 30% of the fibers showed orientation angles between 45° and 50°, while most of the fibers (approximately 87%~90%) ranged from 25° to 65° according to the analysis of the orientation angles.
(4) When the fiber content was lower than 2.0 wt.%, the distribution and orientation of the carbon fibers showed no apparent effect on the electrical resistivity. As the fiber content exceeded 2.0 wt.%, a conductive network of carbon fibers was formed, and the electrical resistivity was further reduced. However, excessive carbon fiber content would cause fiber clusters, leading to a decreasing uniformity of the fiber dispersion. A quantitative relationship between the distribution coefficient (αf) and electrical resistivity (ρ ) was obtained based on the Taipalus model, and the expression was = 0.9236 − 448.0652 ⁄ , with an R-square value of 0.976. This formula can be used to predict the dispersion of carbon fibers with contents higher than 2 wt.% by measuring the electrical resistivity.

Conclusions
(1) The carbon fiber content has a slight influence on the compressive strength and bulk density of amorphous calcium silicate hydrate-based contact-hardening composites, but it can greatly improve the flexural strength of the composites. When the carbon fiber content is 0.5, 1.0, 2.0, 3.0, and 4.0 wt.%, the flexural strength increases by 141, 216, 367, 408, and 558%, respectively, compared to the contact-hardening composite without fibers.
(2) The electrical resistivity of the specimens decreased with increasing carbon fiber content. A lower threshold for contact-hardening composites with carbon fibers appeared at a fiber content of 1.0 wt.%, and the electrical resistivity was 30,000 Ω·cm. A higher threshold was present at a fiber content of 2.0 wt.% with an electrical resistivity of 2550 Ω·cm. It was found that a carbon fiber content of 2.0 wt.% in the contact-hardening composite provided optimum uniformity of the carbon fiber dispersion in the contact-hardening composites.
(3) Direct analysis of the fiber distribution and orientation by using the imaging method suggests that by increasing the fiber content up to 2.0 wt.%, the value of the fiber distribution coefficient (α f ) increased dramatically to 0.743, representing a more homogeneous distribution of fibers in the composite. When the fiber content was beyond 2 wt.%, the value of α f decreased with increasing fiber content, which is ascribed to the fiber cluster in the matrix. Approximately 30% of the fibers showed orientation angles between 45 • and 50 • , while most of the fibers (approximately 87%~90%) ranged from 25 • to 65 • according to the analysis of the orientation angles.
(4) When the fiber content was lower than 2.0 wt.%, the distribution and orientation of the carbon fibers showed no apparent effect on the electrical resistivity. As the fiber content exceeded 2.0 wt.%, a conductive network of carbon fibers was formed, and the electrical resistivity was further reduced. However, excessive carbon fiber content would cause fiber clusters, leading to a decreasing uniformity of the fiber dispersion. A quantitative relationship between the distribution coefficient (α f ) and electrical resistivity (ρ c ) was obtained based on the Taipalus model, and the expression was α f = 0.9236 − 448.0652/ρ c , with an R-square value of 0.976. This formula can be used to predict the dispersion of carbon fibers with contents higher than 2 wt.% by measuring the electrical resistivity.
(5) The flexural strength (σ c f t ) of the contact-hardening composites was dependent on the effective carbon fiber numbers, the product of the total number of the fibers (X total ) in the composites and the distribution coefficient (α f ). The expression was σ c f t = 573α f X total + 135,155 /10,000.
Funding: This research was funded by the Natural Science Foundation of China, grant number 51678093, and the Research Funds for the Central Universities, grant number 2020CDJ-LHZZ-088.

Data Availability Statement:
The data presented in this study are available in the article.