Short-Term Flexural Stiffness Prediction of CFRP Bars Reinforced Coral Concrete Beams

FRP (Fiber Reinforced Polymer) Bar reinforced coral concrete beam is a new type of structural member that has been used more and more widely in marine engineering in recent years. In order to study and predict the flexural performance of CFRP reinforced coral concrete beams, the flexural rigidity, crack morphology and failure mode of concrete were studied in detail. The results show that under the condition of similar reinforcement ratio, the flexural rigidity of CFRP reinforced coral concrete beam is significantly lower than that of ordinary reinforced concrete beam. Increasing the cross-section reinforcement ratio within a certain range can increase the bending stiffness of the test beam or reduce the deflection, but the strength utilization rate of CFRP reinforcement is greatly reduced. The short-term bending stiffness of the CFRP reinforced coral concrete beam calculated by the existing standard formula is obviously higher. This paper proposes a modified formula for introducing the strain inhomogeneity coefficient (ψ) of CFRP bars and considers the relative slip between CFRP bars and coral concrete to predict the short-term flexural stiffness of coral concrete beams reinforced by CFRP bars. The formula was verified with the test results, and it was proved that the formula has a good consistency with the test results.


Introduction
Recently, coral concrete has been gradually used in the coastal area and islands for building construction due to its environment friendly feature. coral concrete with the use of coral reefs and other island reef deposits as aggregates and seawater instead of freshwater can greatly reduce the production and transportation costs. It is major composition, coral debris, can be easily found from the ocean. This new type of product therefore has a high economic value without damaging the natural environment. However, corrosion of reinforcing bars caused by the humidity and hot marine environment in the coastal area and the large amount of salt in coral debris cause the main durability problem of coral concrete [1,2]. Fiber Reinforced Polymer Bar (FRP Bar) has been widely used due to its excellent mechanical properties and excellent corrosion resistance, which provides a new solution for the durability of coral concrete structures [3][4][5]. However, due to the unique composition of coral concrete material, the research on FRP reinforced coral concrete structure is rare. The structural behavior and performance of FRP reinforced coral concrete structure are not fully understood.
FRP reinforced concrete beams are characterized by large deflection, wide crack width and no obvious yield stage [6][7][8][9]. Gdoutos et al. [10] pointed out that the lower elastic modulus and unique bonding characteristics of FRP bars are the main reasons for the larger Materials 2021, 14, 467 2 of 18 deformation of FRP bars reinforced concrete members than that of reinforced concrete members with the same concrete strength, load, component size and reinforcement ratio. Sandor et al. [11,12] showed that the direction of the radial force component of FRP bars in concrete is quite different from that of steel bars, and the FRP bar specimens are more prone to split failure. Zhu et al. [13] discovered that short-term flexural rigidity of FRP bars reinforced concrete beams is easy to be overestimated because the potential influence of bond-slip of FRP bars is not taken into account. Wang et al. [14] compared the calculation results of 745 crack data with different national codes and found that the existing codes overestimated the crack width of FRP reinforced coral concrete beams, which affected the evaluation of the flexural performance of the members. Issa et al. [15] suggest that the deformability of FRP reinforced concrete beams need to be checked. Gravina et al. [16] predict the flexural capacity of FRP reinforced concrete beams with the consideration of the spacing and width of flexural cracks, and the bonding performance between FRP bars. Raed et al. [17] believe that the assumption of plane section remain plane in the calculation of deflection of FRP reinforced concrete beams and slabs under high-level load will lead to a significant reduction in the rigidity of FRP reinforced concrete members, resulting in significant deflection. Wang et al. [18] found that the flexural performance of concrete beams of different types of reinforcement is significantly different, and the application of FRP reinforcement to different types of concrete is also different. The existing calculation model is difficult to accurately predict their flexural performance. Sun et al. [19] pointed out that the ultimate displacement of FRP bars test beam is 1.6-1.9 times greater than that of steel bars test beam due to the difference in bond behavior. Xu et al. [13,20] analyzed different calculation models based on a large number of experimental data and deduced the calculation formulas. The results show that the existing calculation models are still ignored the unique feature of high tensile strength, low elastic modulus, brittle failure mode and low bond performance of this unique type of structure. Therefore, there are inaccurate to predict the short-term stiffness for this particular type of structure.
In order to understand the changing characteristics of the flexural performance of the test beam under different loads, accurately calculate and predict the flexural rigidity of the FRP bar-reinforced coral concrete beam, and establish a calculation model for the flexural stiffness, the flexural behavior, short-term stiffness and failure modes of CFRP reinforced coral concrete beams are studied in detail through full scale flexural tests.

Test Program
Six coral-concrete beams reinforced with fiber bars and two coral-concrete beams reinforced with steel were tested by four-point bending test. All test beams were designed according to existing fiber-reinforced concrete codes and were load-graded using force control. Coral concrete is a kind of lightweight aggregate concrete.

Test Specimens
The CFRP bars used in the tests were manufactured by Zhejiang Haining Anjie Composite Material Co., Ltd. (Haining, Zhejiang, China). Their basic mechanical properties of shown in Table 1. The cement is P.042.5 ordinary Portland cement produced in Xing'an, Guangxi. The coarse aggregate is natural continuous graded coral debris (Figure 1) from North Harbor of China South Sea. The fine aggregate is made of ordinary natural river sand. The water is artificial seawater and polycarboxylic acid superplasticizer is added. At the same time, six cubic and six prism concrete specimens were casted and tested after 28 days of curing under the same conditions as the test beams. The axial compressive strength, tensile spitting strength and elastic modulus of coral concrete were measured. The mix ratio and basic mechanical parameters of coral concrete were shown in Table 2. Eight beams were fabricated with the length of 2400 mm. Their net span is 2100 mm. The sectional dimension is 120 × 250 (mm 2 ), the cover of the concrete is 25 mm, and the diameters of the longitudinal bars are 8 mm, 10 mm and 12 mm (Figure 2), respectively. The basic parameters are shown in Table 3.      Eight beams were fabricated with the length of 2400 mm. Their net span is 2100 mm. The sectional dimension is 120 × 250 (mm 2 ), the cover of the concrete is 25 mm, and the diameters of the longitudinal bars are 8 mm, 10 mm and 12 mm (Figure 2), respectively. The basic parameters are shown in Table 3. Eight beams were fabricated with the length of 2400 mm. Their net span is 2100 mm. The sectional dimension is 120 × 250 (mm 2 ), the cover of the concrete is 25 mm, and the diameters of the longitudinal bars are 8 mm, 10 mm and 12 mm (Figure 2), respectively. The basic parameters are shown in Table 3.

Instrumentation and Loading Process
Similar to test set up in [21][22][23][24][25], as it shown in Figure 3a,b, strain gauges are symmetrically arranged at 200 mm intervals on the CFRP bars and five strain gauges (BA120-05AA) are pasted at equal intervals on the front and back surface of concrete in the middle of the span. All strain gauges are connected to a strain tester (DONGHUA DH3816N, Jingjiang, Jiangsu, China) for measurement. It can be also seen that dial gauges are placed at the midspan, loading point and support of the beam respectively to measure the deformations at these positions. In addition, dial gauges are placed at the end of the longitudinal reinforcement at beam end to measure the end slip of the beams.
Four-point bending is adopted, which is loaded by a reaction frame and a hydraulic jack (PATELI YEYA QF50T-20k, Hongkou, Shanghai, China) with a load transducer (HUADONG ELEC-TECH BHR-4 30T, Jingan, Shanghai, China). During the tests, the load is increased with 3 kN increment at each step. Each loading step is held for 5 min to ensure the instrumentation data is stabilized and the test data are recorded before the load was continuously to increase until the specimen is destroyed. In addition, the electronic crack monitor device (ZBL-F103, Xicheng, Beijing, China) is used to measure the crack width of the concrete at the height of the longitudinal rebars, and the crack development height and crack development speed are measured and recorded.

Instrumentation and Loading Process
Similar to test set up in [21][22][23][24][25], as it shown in Figure 3a,b, strain gauges are symmetrically arranged at 200 mm intervals on the CFRP bars and five strain gauges (BA120-05AA) are pasted at equal intervals on the front and back surface of concrete in the middle of the span. All strain gauges are connected to a strain tester (DONGHUA DH3816N, Jingjiang, Jiangsu, China) for measurement. It can be also seen that dial gauges are placed at the midspan, loading point and support of the beam respectively to measure the deformations at these positions. In addition, dial gauges are placed at the end of the longitudinal reinforcement at beam end to measure the end slip of the beams.
Four-point bending is adopted, which is loaded by a reaction frame and a hydraulic jack (PATELI YEYA QF50T-20k, Hongkou, Shanghai, China) with a load transducer (HUADONG ELEC-TECH BHR-4 30T, Jingan, Shanghai, China). During the tests, the load is increased with 3 kN increment at each step. Each loading step is held for 5 min to ensure the instrumentation data is stabilized and the test data are recorded before the load was continuously to increase until the specimen is destroyed. In addition, the electronic crack monitor device (ZBL-F103, Xicheng, Beijing, China) is used to measure the crack width of the concrete at the height of the longitudinal rebars, and the crack development height and crack development speed are measured and recorded.

Test Result
According to the test results, the bending performance of the test beams with different reinforcements is analyzed and compared. The development rate of deflection and crack width of the fiber reinforced test beam at different stages is different from that of the reinforced test beam. The fiber reinforcement test beam has problems such as large crack width, relative slippage of the reinforcement, and low utilization rate of the reinforcement. The stiffness changes under different load levels are worthy of further discussion.

Load Midspan Deflection Relationship
From the load-midspan deflection curve of the beams in Figure 4, it can be seen that the mid-span deflection of the beams increased with the increase of load, but the deflection curves of the beams with different type of reinforcement and different reinforcement ratios are significantly different. The load-deflection curves of CFRP reinforced coral concrete beams are not significantly different from those of conventional reinforced concrete beams. The conventional reinforced concrete beams can be roughly divided into three working stages: pre-cracking, post-cracking elasticity and reinforcement yielding. Since CFRP bar is a linear elastic material, the load-deflection curve of CFRP reinforced coral concrete beams can be divided into two working stages: no obvious yield stage and relatively sudden failure.
Before concrete crack formed in tension zone, the increase of deflection is slower with the increase of load. After concrete cracking, the stiffness of CFRP reinforced coral concrete beams decreased rapidly, and the increase rate of deflection is obviously accelerated. Under the same load conditions, the deflection of CFRP reinforced coral concrete beams is obviously larger than that of reinforced coral concrete beams. The ultimate bearing capacity of C-8 test beams is not much different from that of steel reinforced coral concrete beams, but the mid-span deflection of CFRP reinforced coral concrete beams under the same load is about 3-4 times of that of latter. Within a certain range, increasing the reinforcement ratio of CFRP bars can significantly improve the flexural rigidity and reduce the deflection. The maximum deflection of C-10 and C-12 beams is 10.8% and 27% lower than that of C-8 beams, respectively. However, with the increase of CFRP reinforcement ratio, the low strength utilization of FRP bars becomes more serious. According to GB 50608-2010 [26], the maximum deflection of l 0 /200 (l 0 is the net span of the test beam) is allowed. When the midspan deflection of CFRP reinforced coral concrete beams reach l 0 /200 (10.5 mm), the load withstood by the CFRP reinforced coral concrete beam is only 23.5-38.5% of their ultimate bearing capacity. Under the same deflection, the steel reinforced coral concrete beams bear about 90% of their ultimate load. The test results are shown in Table 4.

Test Result
According to the test results, the bending performance of the test beams with different reinforcements is analyzed and compared. The development rate of deflection and crack width of the fiber reinforced test beam at different stages is different from that of the reinforced test beam. The fiber reinforcement test beam has problems such as large crack width, relative slippage of the reinforcement, and low utilization rate of the reinforcement. The stiffness changes under different load levels are worthy of further discussion.

Load Midspan Deflection Relationship
From the load-midspan deflection curve of the beams in Figure 4, it can be seen that the mid-span deflection of the beams increased with the increase of load, but the deflection curves of the beams with different type of reinforcement and different reinforcement ratios are significantly different. The load-deflection curves of CFRP reinforced coral concrete beams are not significantly different from those of conventional reinforced concrete beams. The conventional reinforced concrete beams can be roughly divided into three working stages: pre-cracking, post-cracking elasticity and reinforcement yielding. Since CFRP bar is a linear elastic material, the load-deflection curve of CFRP reinforced coral concrete beams can be divided into two working stages: no obvious yield stage and relatively sudden failure.
Before concrete crack formed in tension zone, the increase of deflection is slower with the increase of load. After concrete cracking, the stiffness of CFRP reinforced coral concrete beams decreased rapidly, and the increase rate of deflection is obviously accelerated. Under the same load conditions, the deflection of CFRP reinforced coral concrete beams is obviously larger than that of reinforced coral concrete beams. The ultimate bearing capacity of C-8 test beams is not much different from that of steel reinforced coral concrete beams, but the mid-span deflection of CFRP reinforced coral concrete beams under the same load is about 3-4 times of that of latter. Within a certain range, increasing the reinforcement ratio of CFRP bars can significantly improve the flexural rigidity and reduce the deflection. The maximum deflection of C-10 and C-12 beams is 10.8% and 27% lower than that of C-8 beams, respectively. However, with the increase of CFRP reinforcement ratio, the low strength utilization of FRP bars becomes more serious. According to GB 50608-2010 [26], the maximum deflection of l0/200 (l0 is the net span of the test beam) is allowed. When the midspan deflection of CFRP reinforced coral concrete beams reach l0/200 (10.5 mm), the load withstood by the CFRP reinforced coral concrete beam is only 23.5-38.5% of their ultimate bearing capacity. Under the same deflection, the steel reinforced coral concrete beams bear about 90% of their ultimate load. The test results are shown in Table 4.  Load (kN)  (1) ρ f = Reinforcement ratio. (2) w = Average crack width at time of cracking. (3) w = Maximum crack width when test beam is destroyed. (4) w = Average crack width when the test beam is broken. (5) P = When the deflection of the test beam reaches l 0 /200, the corresponding external load value.

Crack Width
From Figure 5, it can be seen that under the same load conditions, the crack width of the CFRP bar reinforced beam is obviously larger than that of the steel reinforced beam, and the crack propagation height is also higher than that of the reinforced concrete beam. The average crack height is 40 mm for steel reinforced beam, and the average crack height accounts for 70% of the beam height when it reaches ultimate capacity, while the average crack height of CFRP bar reinforced beam is 140 mm, and the maximum crack height accounts for more than 90% of the beam height. The main reason is that after cracking of the concrete in the tension zone at the bottom of the test beam, the tensile stress at the cracks will be borne by the steel bar or CFRP bar. Because the elastic modulus of CFRP bar is obviously smaller than that of the steel bar, the bonding performance between CFRP bar and concrete is relatively weak [27][28][29][30][31], incurs larger deflection. With the increase of reinforcement ratio of CFRP bars, the crack development speed is slowed down, the average crack width and maximum crack width are obviously reduced, but the influence on crack spacing is not obvious. The number of cracks in each beam is between 11 and 15. It should be noted that the maximum crack width and average crack width of all test beams exceed the maximum width limit of 0.5 mm stipulated by current code [26].

Crack Width
From Figure 5, it can be seen that under the same load conditions, the crack width of the CFRP bar reinforced beam is obviously larger than that of the steel reinforced beam, and the crack propagation height is also higher than that of the reinforced concrete beam. The average crack height is 40 mm for steel reinforced beam, and the average crack height accounts for 70% of the beam height when it reaches ultimate capacity, while the average crack height of CFRP bar reinforced beam is 140 mm, and the maximum crack height accounts for more than 90% of the beam height. The main reason is that after cracking of the concrete in the tension zone at the bottom of the test beam, the tensile stress at the cracks will be borne by the steel bar or CFRP bar. Because the elastic modulus of CFRP bar is obviously smaller than that of the steel bar, the bonding performance between CFRP bar and concrete is relatively weak [27][28][29][30][31], incurs larger deflection. With the increase of reinforcement ratio of CFRP bars, the crack development speed is slowed down, the average crack width and maximum crack width are obviously reduced, but the influence on crack spacing is not obvious. The number of cracks in each beam is between 11 and 15. It should be noted that the maximum crack width and average crack width of all test beams exceed the maximum width limit of 0.5mm stipulated by current code [26].

Load and Strain Relation for Longitudinal Bar
The development trend of load-strain curve of CFRP reinforced coral concrete beams is similar to that of ordinary reinforced concrete beams [11][12][13][14]. For ordinary reinforced concrete beams, it can be divided into three stages: slow growth of strain before cracking, uniform growth of longitudinal reinforcement after cracking, yield of reinforcement to failure, as shown in Figure 6. The load-strain curve of CFRP reinforced coral concrete beams can be divided into two stages: pre-cracking stage and post-cracking stage [7][8][9].
Before the concrete cracks in the bottom tension zone, the strain values of steel reinforced coral concrete beams and CFRP reinforced coral concrete beams are similar because concrete withstand more tension stress. After cracking, the strain of CFRP bars increases sharply with the increase of load, which is obviously larger than that of steel bars under the same load value. The main reason is that the bond strength between CFRP bars and coral concrete is weaker than that of reinforcing bars. The relative slip between CFRP bars and coral concrete results in a greater redistribution of the tensile stress of CFRP bars. With the increase of reinforcement ratio of CFRP bars, the slope of load-strain curve increases, and the strain of CFRP bars decreases under the same load. When the midspan deflection of CFRP reinforced coral concrete beams reaches l 0 /200 (10.5 mm), the strain of CFRP rebars under tension is about 30-41.9% of their ultimate strain, which indicates that the strength utilization ratio of CFRP bars is low. In addition, bond slip makes the strain of CFRP bar increase nonlinearly.

Load and Strain Relation for Longitudinal Bar
The development trend of load-strain curve of CFRP reinforced coral concrete beams is similar to that of ordinary reinforced concrete beams [11][12][13][14]. For ordinary reinforced concrete beams, it can be divided into three stages: slow growth of strain before cracking, uniform growth of longitudinal reinforcement after cracking, yield of reinforcement to failure, as shown in Figure 6. The load-strain curve of CFRP reinforced coral concrete beams can be divided into two stages: pre-cracking stage and post-cracking stage [7][8][9].
Before the concrete cracks in the bottom tension zone, the strain values of steel reinforced coral concrete beams and CFRP reinforced coral concrete beams are similar because concrete withstand more tension stress. After cracking, the strain of CFRP bars increases sharply with the increase of load, which is obviously larger than that of steel bars under the same load value. The main reason is that the bond strength between CFRP bars and coral concrete is weaker than that of reinforcing bars. The relative slip between CFRP bars and coral concrete results in a greater redistribution of the tensile stress of CFRP bars. With the increase of reinforcement ratio of CFRP bars, the slope of load-strain curve increases, and the strain of CFRP bars decreases under the same load. When the midspan deflection of CFRP reinforced coral concrete beams reaches l0/200 (10.5 mm), the strain of CFRP rebars under tension is about 30-41.9% of their ultimate strain, which indicates that the strength utilization ratio of CFRP bars is low. In addition, bond slip makes the strain of CFRP bar increase nonlinearly.

Bond Slip of CFRP Reinforced Coral Concrete
It should be noted that almost all CFRP bars reinforced coral concrete beams have slip between CFRP bars and coral concrete in varying degrees. The larger the load, the more obvious the slip. The relative slip curve of CFRP bars at the end of the test beams is shown in Figure 7 (There is no obvious slip in test beam C-10-2). It can be seen that the smaller the diameter of CFRP bars is, the larger the slip is. The slip of CFRP bars has a great influence on the flexural behavior of coral concrete beams. For the same reinforced beams, the flexural capacity of C-8-1 with larger slip is 3.6% lower than that of C-8-2, C-10-1 is 8.8% lower than that of C-10-2, and C-12-1 is 7.7% lower than that of C-12-2. It shows that the larger the diameter of CFRP bars, the smaller flexural capacity due to the slip. However, the chloride ions in artificial seawater and coral have no obvious effect on the surface of CFRP bars.

Bond Slip of CFRP Reinforced Coral Concrete
It should be noted that almost all CFRP bars reinforced coral concrete beams have slip between CFRP bars and coral concrete in varying degrees. The larger the load, the more obvious the slip. The relative slip curve of CFRP bars at the end of the test beams is shown in Figure 7 (There is no obvious slip in test beam C-10-2). It can be seen that the smaller the diameter of CFRP bars is, the larger the slip is. The slip of CFRP bars has a great influence on the flexural behavior of coral concrete beams. For the same reinforced beams, the flexural capacity of C-8-1 with larger slip is 3.6% lower than that of C-8-2, C-10-1 is 8.8% lower than that of C-10-2, and C-12-1 is 7.7% lower than that of C-12-2. It shows that the larger the diameter of CFRP bars, the smaller flexural capacity due to the slip. However

Failure Mode of the Beams
The failure modes of CFRP reinforced coral concrete beams are significantly different from those of steel reinforced coral concrete beams due to the different properties of reinforcement materials. The crack width and deflection of CFRP reinforced coral concrete beams are obviously larger than that of steel reinforced coral concrete beams. In addition, CFRP reinforced coral concrete beams have a large number of deep cracks along the longitudinal reinforcement direction. The slip between CFRP bars and coral concrete is also very serious. As shown in Figure 8, when the beam C-10-1 is damaged, the CFRP bars was exposed at the ends of the beams for about 16 mm. There is no slip phenomenon in steel reinforced coral concrete beams and obvious cracks along longitudinal reinforcement, which indicates that good bonding performance between reinforcement and coral concrete is maintained. However, when the deflection of the test beam increases, the concrete in the compression zone is damaged more obviously due to stress concentration. And when the CFRP bar and coral concrete slip locally, the radial component of the FRP bar causes the concrete protective layer to crack or peel off along the longitudinal direction of the bar. Moreover, the brittleness of coral concrete is greater than that of ordinary concrete, the concrete in the upper compressive zone of the middle-span interface of the test beam has a larger compressive damage range and obvious longitudinal cracks along the test beam, as shown in Figure 8a,b.
Due to the slip of CFRP bars has a significant effect on the deflection and stiffness calculation of test beams, the relatively weak bond between FRP bars and concrete makes the strain of CFRP bars more uniform than that of steel bars [28]. When calculating flexural stiffness, the effect of relative slip of CFRP bars and coral concrete on stiffness reduction should be fully considered. In addition, the crack distribution of the test beam is shown in Figure 9.

Failure Mode of the Beams
The failure modes of CFRP reinforced coral concrete beams are significantly different from those of steel reinforced coral concrete beams due to the different properties of reinforcement materials. The crack width and deflection of CFRP reinforced coral concrete beams are obviously larger than that of steel reinforced coral concrete beams. In addition, CFRP reinforced coral concrete beams have a large number of deep cracks along the longitudinal reinforcement direction. The slip between CFRP bars and coral concrete is also very serious. As shown in Figure 8, when the beam C-10-1 is damaged, the CFRP bars was exposed at the ends of the beams for about 16 mm. There is no slip phenomenon in steel reinforced coral concrete beams and obvious cracks along longitudinal reinforcement, which indicates that good bonding performance between reinforcement and coral concrete is maintained. However, when the deflection of the test beam increases, the concrete in the compression zone is damaged more obviously due to stress concentration. And when the CFRP bar and coral concrete slip locally, the radial component of the FRP bar causes the concrete protective layer to crack or peel off along the longitudinal direction of the bar. Moreover, the brittleness of coral concrete is greater than that of ordinary concrete, the concrete in the upper compressive zone of the middle-span interface of the test beam has a larger compressive damage range and obvious longitudinal cracks along the test beam, as shown in Figure 8a,b.
Due to the slip of CFRP bars has a significant effect on the deflection and stiffness calculation of test beams, the relatively weak bond between FRP bars and concrete makes the strain of CFRP bars more uniform than that of steel bars [28]. When calculating flexural stiffness, the effect of relative slip of CFRP bars and coral concrete on stiffness reduction should be fully considered. In addition, the crack distribution of the test beam is shown in Figure 9.

Failure Mode of the Beams
The failure modes of CFRP reinforced coral concrete beams are significantly different from those of steel reinforced coral concrete beams due to the different properties of reinforcement materials. The crack width and deflection of CFRP reinforced coral concrete beams are obviously larger than that of steel reinforced coral concrete beams. In addition, CFRP reinforced coral concrete beams have a large number of deep cracks along the longitudinal reinforcement direction. The slip between CFRP bars and coral concrete is also very serious. As shown in Figure 8, when the beam C-10-1 is damaged, the CFRP bars was exposed at the ends of the beams for about 16 mm. There is no slip phenomenon in steel reinforced coral concrete beams and obvious cracks along longitudinal reinforcement, which indicates that good bonding performance between reinforcement and coral concrete is maintained. However, when the deflection of the test beam increases, the concrete in the compression zone is damaged more obviously due to stress concentration. And when the CFRP bar and coral concrete slip locally, the radial component of the FRP bar causes the concrete protective layer to crack or peel off along the longitudinal direction of the bar. Moreover, the brittleness of coral concrete is greater than that of ordinary concrete, the concrete in the upper compressive zone of the middle-span interface of the test beam has a larger compressive damage range and obvious longitudinal cracks along the test beam, as shown in Figure 8a,b.
Due to the slip of CFRP bars has a significant effect on the deflection and stiffness calculation of test beams, the relatively weak bond between FRP bars and concrete makes the strain of CFRP bars more uniform than that of steel bars [28]. When calculating flexural stiffness, the effect of relative slip of CFRP bars and coral concrete on stiffness reduction should be fully considered. In addition, the crack distribution of the test beam is shown in Figure 9.

Existing Formula for Short Term Flexural Rigidity
This section mainly compares some standard calculation formulas or calculation models with the test results and discusses the differences between different calculation formulas and cause analysis.

Short Term Flexural Rigidity Formula
Because there are few studies on the flexural behavior of FRP reinforced coral concrete, therefore, no corresponding formula for calculating flexural stiffness. There is some formula for FRP bar reinforced normal concrete can be referred to. They are from ACI 440.1R-15 [32], Zhu Hong [13] and Chinese FRP design code [26].

Formula from ACI 440.1R-15 [32]
The stiffness calculation formula of ACI 440.1R-15 is based on the code ACI 440.1R-06 [33] and ACI 318-05 [34]. Considering the load and boundary conditions and accounts for the length of the uncracked regions of the member and for the change in stiffness in the cracked regions. the stiffness factor γ is introduced in calculating the effective moment of inertia. It is as follows:

Existing Formula for Short Term Flexural Rigidity
This section mainly compares some standard calculation formulas or calculation models with the test results and discusses the differences between different calculation formulas and cause analysis.

Short Term Flexural Rigidity Formula
Because there are few studies on the flexural behavior of FRP reinforced coral concrete, therefore, no corresponding formula for calculating flexural stiffness. There is some formula for FRP bar reinforced normal concrete can be referred to. They are from ACI 440.1R-15 [32], Zhu Hong [13] and Chinese FRP design code [26] The stiffness calculation formula of ACI 440.1R-15 is based on the code ACI 440.1R-06 [33] and ACI 318-05 [34]. Considering the load and boundary conditions and accounts for the length of the uncracked regions of the member and for the change in stiffness in the cracked regions. the stiffness factor γ is introduced in calculating the effective moment of inertia. It is as follows: where I e is the effective moment of inertia of FRP bar reinforced beam after cracking, M cr is the moment when concrete cracks, M a is applied moment, I g is the moment of inertia of uncracked section, I cr is the moment of inertia of uncracked section, γ is the parameter to account for the variation in stiffness along the length of the member. β d is stiffness reduction factor.

GB 50608-2010 (Chinese Standard 2011)
where, B s is the Test beam bending stiffness, ψ is strain inhomogeneity coefficient of CFRP bars, E f is the elastic modulus of fiber, A f is effective cross-sectional area, h of is the distance from the center of the longitudinal bar to the edge of the compression surface, α f E is ratio of elastic modulus of fiber to concrete, ρ f is ratio of FRP bars in longitudinal tension: Zhu [13] further modified Formula (3) where f tk is the tensile strength of concrete, σ f k is the tensile stress of longitudinal bars, ρ te is the ratio of longitudinal tensile FRP bars calculated according to the effective tensile concrete section area, M k is the bending moment calculated by standard combination of load effects. The other symbols have the same meaning as Formula (3).

Comparison to Test Results
The formula for deflection is as follows: where f is the midspan deflection, P is the load, l 0 is the clear span, E c is the Young's module of coral concrete, I e is the moment of inertia of the section, E c I e is the flexural rigidity of the beam. From Figure 10, it can be seen that the short-term stiffness of CFRP reinforced coral concrete is obviously overestimated by the calculation formulas given by Chinese and American codes, and the difference of deflection calculation increases obviously with the increase of load. Among them, the effective moment of inertia calculation formula in ACI 440.1R-15 is modified on the basis of ACI 440.1R-06 and the stiffness reduction coefficient β d is not considered, which makes the pre-load deflection close to the calculated value, but with the increase of load, the deflection discrepancy becomes more serious. Although the calculation curve obtained by the modified formula proposed by Zhu Hong et al. [13] effectively reduces the error, the calculation value is still large when applied to the deflection calculation of coral concrete flexural members. There are two main reasons for the large discrepancy in the calculation. One is that the current formulas for calculating the stiffness of flexural members with FRP bars are mostly based on the formulas for design the normal reinforced concrete. There is not enough correlation between the elastic modulus of FRP bars and bonding performance to and that of steel reinforced concrete, this result in that the calculation are quite different to the actual result; secondly, the elastic modulus of coral concrete is lower than that of ordinary concrete of the same grade, and the bonding of FRP bars are relatively weak. The relative slip between FRP bars and coral concrete is larger. In addition, the compressive strength coral concrete in failure is higher than that of normal concrete. The brittleness of concrete is greater than that of ordinary concrete.
It should be noted that a large number of experimental studies show that the loaddeflection curve of FRP reinforced concrete beams under bending is not an ideal linear relation, but presents certain non-linear characteristics, especially in the stage of large load [35][36][37]. This is because the bond between FRP bars and concrete gradually loses with the increase of load, and the inhomogeneity coefficient of FRP bars increases gradually. In this paper, when the test beam is close to failure, the phenomenon of "false yield" mentioned above may even occur due to the large slip of FRP bars. It can be seen that the formula for calculating short-term flexural stiffness of FRP reinforced concrete beams needs further improvement. above may even occur due to the large slip of FRP bars. It can be seen that the formula for calculating short-term flexural stiffness of FRP reinforced concrete beams needs further improvement.

The Modified Formula for Short Term Flexural Rigidity
At present, there are big differences in the specifications or calculation models in different regions, and the factors considered in the formula are also different, which affect the accuracy of the calculation results of the fiber reinforced test beam. In addition, there are few studies on the stiffness of fiber reinforced coral concrete beams, and it is impossible to conduct a more comprehensive evaluation of its mechanical properties. Therefore, this article compares and analyzes the calculation results of different codes or models and proposes a stiffness calculation model for coral concrete beams on the basis of it.

Formula to Calculate ψ from Chinese Code (50010-2010)
According to the existing calculation theory of flexural rigidity of concrete beams, the elastic modulus of materials and the crack characteristics of members have a significant impact on the flexural behavior and stiffness calculation of concrete beams. For reinforced concrete beams with cracks, the concrete segment between cracks actually bears part of the tension stress, and the tensile strain of the longitudinal bar in the pure bending sector is not uniform. Therefore, the determination of strain inhomogeneity coefficient of longi-

The Modified Formula for Short Term Flexural Rigidity
At present, there are big differences in the specifications or calculation models in different regions, and the factors considered in the formula are also different, which affect the accuracy of the calculation results of the fiber reinforced test beam. In addition, there are few studies on the stiffness of fiber reinforced coral concrete beams, and it is impossible to conduct a more comprehensive evaluation of its mechanical properties. Therefore, this article compares and analyzes the calculation results of different codes or models and proposes a stiffness calculation model for coral concrete beams on the basis of it.

Formula to Calculate ψ from Chinese Code (50010-2010)
According to the existing calculation theory of flexural rigidity of concrete beams, the elastic modulus of materials and the crack characteristics of members have a significant impact on the flexural behavior and stiffness calculation of concrete beams. For reinforced concrete beams with cracks, the concrete segment between cracks actually bears part of the tension stress, and the tensile strain of the longitudinal bar in the pure bending sector is not uniform. Therefore, the determination of strain inhomogeneity coefficient of longitudinal tension bars between cracks is very important for the stiffness calculation. Previous studies have shown that the tensile strength of concrete is almost proportional to the bond strength.
Most of the formulas are fully expressed about the important influence of concrete tensile strength and actual reinforcement ratio on the inhomogeneity coefficient of reinforced concrete beams, which is in accordance with the actual stress state. The inhomogeneity coefficient ψ of steel strain in GB50010-2010 [38] is as follows: Obviously, the surface hardness and shear strength of steel bars are much better than that of FRP bars, and the shear failure of concrete is the main manifestation of interfacial slip. When the reinforcement ratio and the stress of steel bars are fixed, the inhomogeneity coefficient in Equation (8) decreases gradually with the increase of the tensile strength of concrete, which is in line with the actual situation. However, the bond-slip failure mechanism between FRP bars and concrete is more complex, and the bond performance between FRP bars and concrete is significantly different [23][24][25][26]. When the bond interface between FRP bars and concrete slips, not only the concrete is damaged by longitudinal shear, but also the surface of FRP bars [31]. When the strength of concrete increases to a certain extent, the damage focuses on the surface of FRP bars. Therefore, when the strength of concrete exceeds a certain range, the bond strength of concrete is only slightly improved with the increase of strength. There is no doubt that the inhomogeneity strain coefficient of FRP bars in flexural members will not keep a linear relationship with the tensile strength of concrete. The derived formula based on the linear relationship cannot correctly reflect the actual stress situation, especially when the concrete strength is high.
Since all the FRP bars of the flexural members have non-uniform strain and all change in a close range, in this paper, forty eight test results of simply supported concrete beams with CFRP, BFRP and GFRP bars of different strengths are investigated based on the experimental results and the available reference [13,[39][40][41][42][43][44][45][46][47][48], and the relationships between the concrete tensile strength and the inhomogeneity coefficient is developed based on these results and is shown in Figure 11 [14]. The strain inhomogeneity coefficient of longitudinal reinforcement between cracks are derived into following formula: 9) tudinal tension bars between cracks is very important for the stiffness calculation. Previous studies have shown that the tensile strength of concrete is almost proportional to the bond strength. Most of the formulas are fully expressed about the important influence of concrete tensile strength and actual reinforcement ratio on the inhomogeneity coefficient of reinforced concrete beams, which is in accordance with the actual stress state. The inhomogeneity coefficientψ of steel strain in GB50010-2010 [38] is as follows: Obviously, the surface hardness and shear strength of steel bars are much better than that of FRP bars, and the shear failure of concrete is the main manifestation of interfacial slip. When the reinforcement ratio and the stress of steel bars are fixed, the inhomogeneity coefficient in Equation (8) decreases gradually with the increase of the tensile strength of concrete, which is in line with the actual situation. However, the bond-slip failure mechanism between FRP bars and concrete is more complex, and the bond performance between FRP bars and concrete is significantly different [23][24][25][26]. When the bond interface between FRP bars and concrete slips, not only the concrete is damaged by longitudinal shear, but also the surface of FRP bars [31]. When the strength of concrete increases to a certain extent, the damage focuses on the surface of FRP bars. Therefore, when the strength of concrete exceeds a certain range, the bond strength of concrete is only slightly improved with the increase of strength. There is no doubt that the inhomogeneity strain coefficient of FRP bars in flexural members will not keep a linear relationship with the tensile strength of concrete. The derived formula based on the linear relationship cannot correctly reflect the actual stress situation, especially when the concrete strength is high.
Since all the FRP bars of the flexural members have non-uniform strain and all change in a close range, in this paper, forty eight test results of simply supported concrete beams with CFRP, BFRP and GFRP bars of different strengths are investigated based on the experimental results and the available reference [13,[39][40][41][42][43][44][45][46][47][48], and the relationships between the concrete tensile strength and the inhomogeneity coefficient is developed based on these results and is shown in Figure 11 [14]. The strain inhomogeneity coefficient of longitudinal reinforcement between cracks are derived into following formula:  Figure 11. Relationship between tensile strength and strain inhomogeneity coefficient of concrete [14].   Figure 11. Relationship between tensile strength and strain inhomogeneity coefficient of concrete [14].

The Modified Formula for Flexural Rigidity
When developing the modified short-term stiffness calculation formula of FRP bars reinforced coral concrete beams, two key factors need to be considered. Firstly, the influence of material properties and surface conditions of FRP bars on flexural stiffness should be considered. The influence factors of bonding performance of different FRP bars to coral concrete should be introduced. Based on the experimental data and the bonding characteristic coefficients β s of reinforced concrete in the code, the bonding characteristic coefficients of reinforcing bars should be taken into account. The value β s is taken as 0.9, and the influencing factor K1 reflecting the surface of FRP is introduced. (The value of K1 is influenced by the diameter, height, spacing and rib type of FRP bars, which is quite complex, and needs further study, however it is conservative to takes as 0.9 here [27,28,30].) Secondly, considering the effect of low modulus of elasticity of CFRP bars and coral concrete, and the crack height of concrete beams is higher resulting lever arm is larger, and the lever arm η is taken as 0.9 from Chinese code GB50608-2010, therefore, for the short-term stiffness Formula (10), the coefficient before ψ should be taken as 1/0.9 = 1.11.
Based on theoretical analysis and experimental data, the formulas for calculating the short-term flexural stiffness of CFRP reinforced coral concrete beams are presented as follows: where β s is the effect coefficient of bond property between FRP bars and coral concrete, K 1 is the effect factor of FRP apparent condition, and the other symbols have the same meaning as Formula (3). Figure 12 is a comparison between the load-midspan deflection curves of CFRP reinforced coral concrete beams calculated by the modified formula and the test results. It can be seen that the calculated values are in good agreement with the experimental values within the range of service loads and tend to be safe. It is noteworthy that when approaching the ultimate load, the deflection of the test beam increases rapidly and destroys due to the larger slip of CFRP bars in coral concrete. Especially when it is close to destruction, which makes the calculation curve of the revised formula deviate to a certain extent, but the difference is still small. Further research is needed if the bending stiffness of the test beam is to be calculated more accurately when it approaches failure.

Discussion
Analyzed from the test results, the existing calculation model tends to overestimate the flexural rigidity of the fiber reinforced coral concrete beam and have certain accuracy only before and after the cracking of the test beams. However, as the load increases, the error becomes larger, and the design requirements for different load levels cannot be met. In the case of test beam failure, the errors of Zhu et al. [13], GB50608-2010 [26], and ACI [48] are 17.1-20.8%, 27.6-28.7%, and 35-36% respectively. The model proposed in this paper has an error of no more than 3% within 80% of the working range, even if the error is only between 7% and 10% when it is broken.
In fact, the calculation model of fiber reinforcement basically refers to the reinforced concrete specification, and does not consider the apparent condition of the reinforcement, the relative slippage of the reinforcement and the influence of the force characteristics at different stages on the components, so it is impossible to accurately calculate the components. From another point of view, if one does not understand the force changes of the components at different loads or at different stages, when the components are damaged during use, they cannot be accurately evaluated and repaired. Therefore, in-depth study of stiffness changes has far-reaching significance for the application of fiber reinforced concrete members.

Discussion
Analyzed from the test results, the existing calculation model tends to overestimate the flexural rigidity of the fiber reinforced coral concrete beam and have certain accuracy only before and after the cracking of the test beams. However, as the load increases, the error becomes larger, and the design requirements for different load levels cannot be met. In the case of test beam failure, the errors of Zhu et al. [13], GB50608-2010 [26], and ACI [48] are 17.1-20.8%, 27.6-28.7%, and 35-36% respectively. The model proposed in this paper has an error of no more than 3% within 80% of the working range, even if the error is only between 7% and 10% when it is broken.
In fact, the calculation model of fiber reinforcement basically refers to the reinforced concrete specification, and does not consider the apparent condition of the reinforcement, the relative slippage of the reinforcement and the influence of the force characteristics at different stages on the components, so it is impossible to accurately calculate the components. From another point of view, if one does not understand the force changes of the components at different loads or at different stages, when the components are damaged during use, they cannot be accurately evaluated and repaired. Therefore, in-depth study of stiffness changes has far-reaching significance for the application of fiber reinforced concrete members.
The calculation model proposed in this paper takes into account the force characteristics and influencing factors of the test beam at different stress stages and can predict the The calculation model proposed in this paper takes into account the force characteristics and influencing factors of the test beam at different stress stages and can predict the stiffness change and deflection deformation of the test beam under different loads. However, there is a lack of long-term mechanical performance test research on test beams and long-term load degradation and damage research in this field, so long-term stiffness prediction cannot be made. In the future work, Finite Element modelling [49] can be further performed, as well consider other parameters such as effect of high strength concrete [50] 6. Conclusions

1.
The flexural behavior of CFRP reinforced coral concrete beams is similar to that of steel reinforced concrete beams with good ductility, but the flexural stiffness of CFRP reinforced coral concrete beams is significantly lower than that of ordinary reinforced concrete beams. There is no obvious yield stage. After the cracking of concrete in tension zone, the stiffness of CFRP reinforced coral concrete beams decreases rapidly and the increase rate of deflection accelerates obviously.

2.
Increasing the reinforcement ratio of CFRP bars in a certain range can significantly improve the flexural rigidity and reduce the deflection of CFRP reinforced coral concrete beams. However, by increasing the reinforcement ratio to reduce the deflection of the test beam, it will reduce the utilization of the strength of the reinforcement.

3.
The bond performance between CFRP bars and coral concrete has a significant impact on the flexural performance and failure mode of CFRP bars. Most of the tested beams show obvious slip phenomenon under tension, which results in serious loss of flexural capacity of the flexural-shear section. 4.
The revised formula for strain inhomogeneity coefficient of CFRP bars in coral concrete beams proposed in this paper has higher accuracy.

5.
A new short-term stiffness formula of CFRP reinforced coral concrete beams is developed, which considering the low elastic modulus of FRP bars, the stress-strain characteristics of coral concrete and the bond performance. The relative slip coefficient is taken as 0.9, and the surface condition coefficient is 0.9. The calculated results are close to the experimental values.