Research Article  Open Access
Zhaoyan Cui, Liuhui Tu, Ming Xu, Zhongfan Chen, Qingfeng Xu, "Experimental Investigation on the LoadCarrying Capacity of SteeltoLaminated Bamboo Dowel Connection I: Single Fastener with SlottedIn Steel Plate under Tension", Advances in Civil Engineering, vol. 2021, Article ID 6683589, 10 pages, 2021. https://doi.org/10.1155/2021/6683589
Experimental Investigation on the LoadCarrying Capacity of SteeltoLaminated Bamboo Dowel Connection I: Single Fastener with SlottedIn Steel Plate under Tension
Abstract
The doweltype connection is widely applied in timber and bamboo structures. It is ambiguous regarding the calculation method of engineered bamboo connections completely referred to the timber design codes. The steeltolaminated bamboo dowel connections with slottedin steel plate tests were conducted to investigate the mechanical performance under tension based on the ASTMD565215. The effects of the thickness, dowel diameter, and end distance on the yield load, ultimate load, initial stiffness, and ductility of the connections were studied. The difference in the yield load for different end distance is negligible. With the same thickness of the connections, the lower the thickness to dowel diameter, the larger the loadcarrying capacity. The three typical yield modes and corresponding loaddisplacement curves of the connections are observed. By considering the rigidplastic model, the theoretical equation for the connections is proposed and proven to fit well with the experimental results. It presents a better prediction for the loadcarrying capacity of steeltolaminated bamboo dowel connections with slottedin steel plate.
1. Introduction
Due to the advantages of a simple manufacturing process, low adhesive content, and beautiful surface, laminated bamboo has a broad prospect as an engineered material in the construction industry [1, 2]. Compared with other biological structural materials, the laminated bamboo has higher tensile strength and stiffness with less variability [1], meaning that, for construction purposes, it would be a better choice than the others.
For the engineered structures, the reliability of the connections is the key to the structural design. The connection performance of laminated bamboo could directly affect the strength, life expectations, and robustness of the structure, making the relevant research extremely important. The dowel connection of the timber and bamboo structure can be divided into singleshear connection and doubleshear connection from the number of shears faces of each dowel. The multishear connections can be divided into a woodtowood connection and a steeltowood connection from the different materials of the connected main and side materials, and the steeltowood connection includes a steel splint connection and a steel filler connection. The simple form and reliable transmission of force make the steeltowood connection the most important form of connections in the modern timber structures, especially those doweltype connections with slottedin steel plates for having no steel plate coverage and good aesthetics; and the occurrence of cracks can be easily observed to avoid potential safety hazards [3].
Researchers have explored the calculation theories and methods of wood bolt connection for years. The European yield theory was originated from Johansen [4] theory of timber connection, which set the foundation for Eurocode 5 [5]. In the United States, the design method has been changed from empirical design to method based on the yield theory since 1991. Several factors, such as the dowel diameter, the thickness to diameter ratio, and the end distance have been investigated to study the possible effect on the loadcarrying capacity of the connection. Daudeville et al. [6] carried out the dowelbearing test on the wood, indicating that the dowel diameter has an influence on the bearing strength. In a certain range, it has a positive correlation with the dowelbearing strength, and the result is consistent with the European wood structure design specification. Solitis et al. [7] proposed the concept of thicknesstodiameter ratio in 1987 and believed that, with the same dowel diameter, the greater the thickness of the bolted joint of wood, the better the ductility. Cesar [8] pointed out in the study that when the dowel end distance to diameter ratio is less than 4, the brittle fracture is easy to occur. Doyle [9] found that an increase in the spacing of the bolts within a certain range helps to enhance the bearing capacity of the joints. The requirements on these bolted constructions directly affect the form of damage to the joints. The loading direction has also been found to have a significant effect on the performance of the connection. The failure modes of the timber loaded parallel and perpendicular to the grain are significantly different. Patel and Hindman [10] studied the bearing properties of the bolted joints in the transverse direction and found that the failure mode is related to the fracture toughness of the wood. In addition, moisture content [11] and temperature can affect the loadcarrying capacity.
At present, the existing equations of the steeltowood connection in the national design codes and standards focus on different failure modes: Eurocode 5: design of timber structures evaluates the bearing capacity of the connection with the thickness of the side members, the dowelbearing capacity of the material, and the diameter of the dowel. The minimum value of predictions under each failure mode is taken as the calculation result to predict the bearing capacity of the connection. Similar methods have been applied in the National Design Specification for Wood Construction in the US [12] and the Canadian Standard [13].
Although the bamboo material shares great similarities with wood at the aspect of physical and mechanical properties, the difference between the two materials is not negligible. There are relatively few studies on the engineered bamboo material. Cui et al. [14] tested the bearing capacity of bolted steelbamboo scrimbersteel connections and found that Eurocode 5 had better prediction accuracy. Hover [15] conducted an experimental study on the bolt connection of the laminated bamboo loaded parallel to the grain and analyzed its failure mechanism and bearing capacity. In general, the research on the bearing performance of laminated bamboo connections is still scarce at home and abroad, and further exploration is urgently needed.
The purpose of the study is to analyze the loadcarrying capacity of the steeltolaminated bamboo connection loaded parallel to the grain. The dowel diameter, the thickness, and the end distance of the connections have been taken into consideration. The experimental results have been compared with the existing equations, and a more accurate equation for laminated bamboo has been developed based on the analysis in the paper.
2. Materials and Test Methods
2.1. Materials
The fabrication process of the raw laminated bamboo material has been described specifically in previous research [16]. Based on the Chinese standard [17], ASTM D14314 [18], and ASTM D576497a [19], the physical and mechanical properties of the laminated bamboo have been tested. The mean values of its airdry and ovendry density are 609 kg/m^{3} and 592 kg/m^{3} with a moisture content of 5.9%. The compressive strength is 59.63 MPa parallel to the grain with the MOE of 12087 MPa. The corresponding tensile strength is 104.16 MPa with the MOE of 10820, and the shear strength is 17.26 MPa. The tensile strength of laminated bamboo perpendicular to grain is 3.35 MPa. The yield strength of the dowel is 480 MPa and the tensile strength is 600 MPa. The length of the unthreaded area in the middle of the screw was not less than 120 mm. The steel plate is Q345 grade steel with a thickness of 10 mm; and the thickness of the groove is 12 mm.
2.2. The Test Specimens and Method
As mentioned before, many factors may influence the bearing capacity of the connections. In this study, the dowel diameter, the thickness, and end distance are selected as reference factors. The test setup and front and lateral views of the connections are shown in Figure 1. The length of laminated bamboo specimen is 900 mm, and the width is 120 mm. It is noted that the end distances of the connections are 5D, 6D, 7D, and 8D except that the control group is 7D. The detailed parameters of different groups are shown in Table 1. Each test was repeated 3 times. The test method was based on ASTMD565215 [20] and was appropriately adjusted in combination with the test conditions. The test was carried out with a 100ton MTS Fatigue testing machine. As shown in Figures 1(b) and 1(c), 1 denotes the MTS actuator, 2 the fixed device for the steel plate, 3 the steel plate, 4 the displacement transducer support, 5 the dowel, 6 the displacement transducer roof, 7 the fixed device for the specimens, and 8 the MTS ground anchor. After preloading to eliminate the effect of the initial gap and deformation, the specimens were subjected to monotonic uniform loading in the tension direction with displacement control. The loading rate was 1 mm/min. The loading process stopped when significant damage appeared or the bearing capacity dropped below 80% of the maximum load. The load of the connections was measured by the test setup. The displacement was measured by two displacement transducers on both sides.
(a)
(b)
(c)
 
D denotes the dowel diameter, t denotes the thickness of the specimens, t_{s} denotes the thickness of the side members, l_{e} denotes the end distance, and γ denotes the thickness to diameter ratio. 
The average value measured by the two displacement meters was chosen as the relative displacement between the dowel and the laminated bamboo, and the load data are combined to generate the loaddisplacement curves. A straight line was fit to the initial linear portion of the loaddeformation curve. The line was offset by a deformation equal to 5% of the fastener diameter. The load at which the offset line intersects the loaddeformation curve was selected as the yield load, F_{y}. The elastic stiffness, K_{1}, is calculated by the slope of the initial linear portion of the loaddeformation curve. The postyield stiffness, K_{2}, is calculated by the tangent of the loaddisplacement curve through the intersection of the initial linear portion of the loaddeformation curve and the horizontal line through the yield load.
3. Results and Analysis
3.1. The Failure Modes and LoadDisplacement Curves of the Connections
According to the Eurocode yield theory, there are three yield modes for doubleshear connections in Figure 2. Mode I represents the case where fasteners do not bend, and wood fibers crushing occurres beneath the dowel without through cracks. Mode II represents the case where fasteners yield in bending at one plastic hinge point per shear plane, with bearingdominated yield of wood fibers in contact with the fasteners in side member(s), respectively. Mode III represents the case where fasteners yield in bending at two plastic hinge points per shear plane, with the limited localized crushing of wood fibers near the shear plane(s).
(a)
(b)
(c)
As shown in Figure 3, all of the three yield modes of the dowels occur in the connections. When the thickness of side members is small and the thickness to diameter ratio is also small, Mode I is prone to occur, for example, T2D1L1. It is mainly because the bending strength of the dowel is larger and the shear strength of the laminated bamboo parallel to grain is lower. But when the thickness to diameter ratio is large enough, Mode III occurs, specimen T1D2L1. It is because the bearing zones become larger and dowel bending is restrained. For the other connections, the failure modes belong to Mode II.
(a)
(b)
(c)
The failure modes of side members under the yield mode of one plastic hinge (Mode II) are mainly shear failure and splitting failure and they are accompanied by embedding failure. When the thickness of side members is small, the ductility is lower and the specimen exhibits splitting failure. As the thickness increases, the plastic displacement increases, and the dowel hole is elongated. The specimens exhibit splitting failure and embedding failure. But when the dowel diameter is larger, for example, M14 and M16, the specimens exhibit shear failure and embedding failure. When the end distance of the connections is in the range of 6D–8D, the connections exhibit better ductility and embedding failure. Therefore, the sufficient end distance of the connections is essential in the design of the laminated bamboo connections.
The typical loaddisplacement curves of the connections are shown in Figure 4. The curves can be divided into a linear branch, a nonlinear branch within and beyond the proportional limit, and falling branch rapidly. The influences of the dowel diameter, thickness, and end distance on the mechanical properties of the connections parallel to grain were analyzed in the paper. As shown in Figure 4(a), with the dowel diameter increasing from 10 mm to 16 mm, the ultimate load, yield load, and initial stiffness increase gradually. When the dowel diameter is 10 mm, the ultimate displacement is about 20 mm and it exhibits better ductility. As the thickness increases, the ultimate load and displacement increase in Figure 4(b). When the thickness is 30 mm, the ultimate load and displacement are too small, just 27.15 kN and 2.49 mm. The ultimate displacement is close in the range of 90–120 mm. The connections have good ductility and the curve change is similar to the same diameter and thickness when the end distance is above 6D as shown in Figure 4(c). Therefore, an end distance of 6 to 8 dowel diameters and a thickness equal to or greater than 7 dowel diameters will provide reasonable results for laminated bamboo connections.
(a)
(b)
(c)
3.2. The Effect of Different Factors
The test results are shown in Table 2. The mean values of a series of measurements of different groups have been collected. The initial stiffness of the connections, K_{1}, yield displacement, Δ_{y}, and yield load, F_{y}, were obtained by the 5% diameter offset method. μ is the ductility ratio, μ = Δ_{u}/Δ_{y}. and COV is the coefficient of variation of F_{u}. If the displacement at the point of intersection is larger than the ultimate displacement, the ultimate load is defined as the yield load. The ultimate load, F_{u}, was obtained as the loadcarrying capacity of the connections.

The influences of the dowel diameter, thickness, and end distance on the yield and ultimate load, initial stiffness, and ductility ratio are compared and analyzed in Figure 5. For T1D1L1 as the control group, the corresponding yield and ultimate load are 32.33 kN and 50.56 kN; and the initial stiffness and ductility ratio are 43.56 kN/mm and 10.26. For the dowel diameter, with the diameter increasing from 10 mm to 16 mm, the yield and ultimate load and the initial stiffness increase gradually in Figure 5(a). But the ductility ratio slightly reduces. It is because the dowel diameter determines the bearing area between the dowel and the bamboo material but the increase of the diameter of the connections could cause the failure mode to change from combined embedding failure and splitting failure to shear failure.
(a)
(b)
(c)
As for the thickness factor, it has the most significant effect on the initial stiffness, ductility ratio, and ultimate load of the connections in Figure 5(b). As the specimen increases from 30 mm to 120 mm, the ultimate load increases rapidly from 27.15 kN to 50.56 kN and the ductility ratio increases from 1.98 to 10.26. But when the thickness of side members is in the range of 90–120 mm, the thickness variation has no significant effect on yield load. It is mainly because the failure modes are similar, and the yield mode is determined by the compressing buckling of the bamboo fiber, which belongs to the local material property and is independent of the overall size of the specimens.
When the end distance of the connections increases from 5D to 8D, the shear area of side members enlarges as well, leading to the ultimate load increasing slightly. There is no obvious correlation between end distance with the yield load and initial stiffness because the yield mode depends on the material properties of the laminated bamboo. The initial stiffness is around 38–45 kN/mm. It is noted that when the end distance is 60 mm, the ductility ratio is very small, just 3.7. It should be avoided in the design of engineered bamboo structures. This phenomenon proved that connections with sufficient end distance could guarantee the ductility of the connections.
4. Calculation
4.1. The Existing Calculation Methods
Regarding the loadcarrying capacity and yield resistance of the connections, different national standards of timber structure have different calculation methods. In Eurocode 5 [5], the three different equations were used to calculate the corresponding failure modes based on the Eurocode yield theory; and the characteristic loadcarrying capacity of the connections can be obtained as the minimum value of the three equations. It can be expressed as follows:where F_{u} is the characteristic loadcarrying capacity per shear plane per fastener, kN, f_{h,1, k} is the embedding strength of timber members, MPa, t_{1} is the thickness of side member, mm, d is the dowel diameter, mm, and M_{y, Rk} is the yield moment of the dowel, N·mm.
In Canadian engineering design in the wood [13], the yield resistance of the connections rather than bearing capacity was calculated according to different failure modes. Like Eurocode 5, the minimum value of the connections was obtained to calculate the yield resistance. It can be expressed as follows:where F_{y} is yield resistance per shear plane per fastener, kN, f_{1} is the embedding strength of the side member, MPa, f_{2} is the embedding strength of the main member, MPa, t_{1} is the thickness of side member, mm, d is the dowel diameter, mm, f_{y} is the yield strength of the dowel, MPa, and ϕ_{y} is resistance factor for yielding failure, 0.8. For the connection with slottedin steel plate, f_{2} can be calculated aswhere K_{sp} is the parameter, and the value of 3.0 is taken for mild steel; ϕ_{steel} is the resistance factor for steel plates in connections with dowels, 0.8 for mild steel, and f_{u} is the specific tensile strength of steel, MPa.
4.2. The Equation for the Laminated Bamboo
The Eurocode yield theory is originated from Johansen theory of timber connection, which set the foundation for national standards of timber structure above. The basic assumptions of the theory are as follows: the failure mode of the connections is ductile and the embedding stressstrain relationship is the ideal rigidplastic model. Then, based on the equilibrium condition of forces and three yield modes of fasteners, the equations to express the loadcarrying capacity of the connections can be determined. The stress distributions under three different yield modes in the rigidplastic model are shown in Figure 6. The red circle denotes plastic hinge, R is the loadcarrying capacity per shear plane per fastener, kN, f_{e} is the embedding strength of timber members, MPa, t is the thickness of side member, mm, and M_{y} is the yield moment of the dowel, N·mm.
(a)
(b)
(c)
The following equations of the connections under three yield modes can be developed to calculate the loadcarrying capacity:
Regarding the yield moment of the dowel, the calculation method is different in national standards of the timber structure. In Chinese timber structure design code [21], the yield moment of the dowel is related to the yield strength and the equation can be expressed aswhere is the ratio of the plastic section modulus to the elastic section modulus, when the plastic behavior of the dowel is fully developed, 1.7, if not, 1.4. W is the moment of section resistance in mm^{3}. F_{y} is the yield strength of the dowel in MPa.
In American design specification for wood construction, the yield moment is related to bending strength of the dowel. It can be determined based on the yield load of 5% diameter deviation method or the average of the tensile strength and yield strength. Generally, it is taken as 1.3 times of yield strength. The equation can be expressed as
In Eurocode 5, the yield moment is related to tensile strength, and the equation can be expressed aswhere f_{t} is the tensile strength of the bolt in MPa.
As mentioned before, M10, M12, M14, and M16 are tested and the yield moments are 72.30 kN mm, 167.70 kN mm, 271.80 kN mm, and 394.80 kN mm, respectively. By comparison with different national standards, the results calculated in American code are the closest with experimental results; and the relative difference is less than 15%. Next, equation (8) can be substituted into equations (5) and (6) to calculate the singleshear bearing capacity of the connections. So, the equation can be expressed as
According to three yield modes of boundary conditions, equation (10) can be further solved, which includes the thickness to diameter ratio interval. It can be shown in Figure 7. So, the upper and lower limits of Mode II are determined, and the equation can be expressed as
The dowelbearing capacity of laminated bamboo is an essential parameter to evaluate the loadcarrying capacity of this kind of connection. Previous tests and researches [22, 23] suggest that the dowelbearing capacity of the laminated bamboo parallel to grain can be expressed aswhere f_{e, par} and f_{c, 0} stand for embedding strength and compressive strength of the laminated bamboo parallel to grain in MPa, respectively.
Based on the above equations, the theoretical and experimental results of loadcarrying capacity and failure modes are shown in Table 3. Compared with the ultimate load, the theoretical results are conservative but the predicted failure modes are accurate. In order to obtain an accurate solution of the loadcarrying capacity of steeltolaminated bamboo dowel connections with slottedin steel plate, the formulas would need further discussion and optimized analysis. The modified coefficient can be introduced and combined with equation (11) to predict the loadcarrying capacity of the connections. For further analysis, the modified coefficient is obtained to evaluate three yield modes. The corresponding values are 1.14 in Mode I, 1.08 in Mode II, and 1.50 in Mode III. When the connection occurs, Mode III, the modified coefficient is the largest. Through the previous analysis, when the end distance of the connections is 6 to 8 dowel diameters, they exhibit better ductility and bearing capacity. Plus, the effect of end distance on the bearing capacity of the connections is not considered in the calculation method. The sufficient end distance could guarantee the accuracy of the results and structural safety. It is noted that, considering the limited number of the connections, the modified coefficient suggested in the paper needs to be further verified. To sum up, the modified equation has better prediction and safety for the steeltolaminated bamboo dowel connection with slottedin steel plate.

5. Conclusion
The study investigated the loadcarrying performance of one doweltype steeltobamboo connection with a slottedin steel plate. The effects of dowel diameter, thickness, and end distance on the yield load, the elastic stiffness, the ultimate load, the plastic stiffness, and the ductility rate of the connections have been studied. The following conclusions can be drawn:(1)The typical loaddisplacement curves and yield modes of the connections parallel to the grain have been studied. The failure modes of side members are mainly shear failure and splitting failure, and they are accompanied by embedding failure. As for the dowel, the onehinge yield mode has been witnessed the most.(2)The influence of different factors on the loadcarrying capacity of the connections has been analyzed as well. The dowel diameter significantly influences the yield and ultimate load and the initial stiffness. When the thickness of side members is in the range of 90–120 mm, the thickness variation has no significant effect on yield load. As for the end distance, when it is in the range of 6D–8D, the connections exhibit better ductility and embedding failure.(3)Compared with experimental results, the theoretical results are conservative but the predicted failure modes are accurate. Further, the modified coefficient is introduced and combined with the theoretical equation to predict the loadcarrying capacity of the connections, which has better prediction and safety for the steeltolaminated bamboo dowel connection with slottedin steel plate.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (NSFC) (51778299) and Open Fund of Shanghai Key Laboratory of Engineering Structure Safety (2019KF05).
References
 A. Sinha, D. Way, and S. Mlasko, “Structural performance of glued laminated bamboo beams,” Journal of Structural Engineering, vol. 140, no. 1, Article ID 04013021, 2013. View at: Google Scholar
 Y. Xiao and J. Ma, “Fire simulation test and analysis of laminated bamboo frame building,” Construction and Building Materials, vol. 34, no. 34, pp. 257–266, 2012. View at: Publisher Site  Google Scholar
 J. Fueyo, M. Dominguez, J. A. Cabezas et al., “Design of connections with metal doweltype fasteners in double shear,” Materials and Structures, vol. 42, no. 3, pp. 385–397, 2009. View at: Google Scholar
 K. Johanson, “Theory of timber connections,” International Association for Bridge and Structural Engineering, vol. 9, pp. 249–262, 1949. View at: Google Scholar
 The European Standard EN 199511: 2004 (E), Eurocode 5: Design of Timber StructuresPart 11: GeneralCommon Rules and Rules for Buildings, 2004.
 L. Daudeville, L. Davenne, and M. Yasumura, “Prediction of the load carrying capacity of bolted timber joints,” Wood Science and Technology, vol. 33, no. 1, pp. 15–29, 1999. View at: Publisher Site  Google Scholar
 L. A. Soltis, F. K. Hubbard, and T. L. Wilkinson, “Bearing strength of bolted timber joints,” Journal of Structural Engineering, vol. 112, no. 9, pp. 2141–2154, 1986. View at: Publisher Site  Google Scholar
 E. César, “Capacity Predictions for Bolted Timber Joints Failing by splitting,” The 10th World Connection Timber Engineering, Miyazaki, Tokyo, Japan, 2008. View at: Google Scholar
 L. Doyle, Performance of Joints with Eight Bolts in Laminated DouglasFir (No. FSRNFPL10), Forest Products Lab Madison Wis., Madison, WI, USA, 1964.
 M. Patel and D. Hindman, “Comparison of singleand twobolted LVL perpendiculartograin connections. II: fracture models,” Journal of Materials in Civil Engineering, vol. 24, no. 4, pp. 347–355, 2011. View at: Google Scholar
 W. Mengjie, “Effect of humidity change on the performance of bolted joints of wood structural steel fillings,” Harbin Industrial University, Harbin, China, 2013, Master’s thesis. View at: Google Scholar
 American Forest & Paper Association, ANSI/AF&PA NDS2012 National Design Specification for Wood Construction, AF&PA, Washington, D.C., USA, 2012.
 Canadian Standard Association, O8601, Engineering Design in Wood, Canadian Standard Association, Toronto, Canada, 2001.
 Z. Cui, M. Xu, Z. Chen et al., “Experimental study on bearing capacity of bolted steelPSBsteel connections,” Engineering Mechanics, vol. 36, no. 1, pp. 96–103, 2019, (In Chinese). View at: Google Scholar
 J. H. P. Hover, “Load carrying capacity of multiple fastener steel to bamboo joints loaded parallel to the fibre direction,” 2017. View at: Google Scholar
 M. Xu, Z. Cui, L. Tu, Q. Xia, and Z. Chen, “The effect of elevated temperatures on the mechanical properties of laminated bamboo,” Construction and Building Materials, vol. 226, pp. 32–43, 2019. View at: Publisher Site  Google Scholar
 Jg/T 1992007, Testing Methods for Physical and Mechanical Properties of Bamboo Used in Building, Standards Press of China, Beijing, China, 2007.
 American Society for Testing and Materials, Standard Test Methods for Small Clear Specimens of Timber, American Society for Testing and Materials, West Conshohocken, PA, USA, 2014.
 American Society for Testing and Materials, Standard Test Method for Evaluating DowelBearing Strength of Wood and WoodBased Products, American Society for Testing and Materials, West Conshohocken, PA, USA, 2013.
 American Society for Testing and Materials, Standard Test Method for SingleBolt Connections in Wood and WoodBased Products, American Society for Testing and Materials, West Conshohocken, PA, USA, 2015.
 China Architecture Industry Press G. B. 50005―2003, Timber Structures Design Code, China Architecture Industry Press, Beijing, China, 2017.
 Z. Cui, L. Tu, M. Xu, Z. Chen, and C. Wang, “The evaluation of dowelbearing properties of laminated bamboo parallel to grain,” Structures, vol. 25, pp. 956–964, 2020. View at: Publisher Site  Google Scholar
 Z. Cui, M. Xu, L. Tu, Z. Chen, and B. Hui, “Determination of dowelbearing strength of laminated bamboo at elevated temperatures,” Journal of Building Engineering, vol. 30, Article ID 101258, 2020. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 Zhaoyan Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.