#### Abstract

The fatigue performance of the bridge deck significantly affects the safety and durability of the overall steel-concrete composite beam bridge. Based on the vehicle flow information of the highway within 10 years, the fatigue performance of a two-way four-lane steel-concrete composite continuous beam bridge deck is studied in this research. The results indicate that the effect of the wheel track position is negligible for two-way four-lane bridge when the wheel track sways laterally, and the fatigue stress of bridge deck concrete is the most unfavorable while the loading position is 7.0 m away from the bridge center line. The fatigue damage decreases by 30%–40% when the centerline of the lane deviates from the most unfavorable stress position by 1 m. The punching fatigue of the concrete is more sensitive to the changes in slab thickness, and the thickness of the deck concrete slab is recommended to be ≥35 cm.

#### 1. Introduction

Compared with concrete bridges, steel-concrete composite structure bridges have the advantages of a lower self-weight and a larger span; compared with steel bridges, they have the advantages of less steel consumption, better structural stability, higher bending rigidity, and higher ductility [1–3]. However, steel-concrete composite beams are prone to fatigue problems under long-term repeated loading [2, 4–6]. The fatigue life of the welded joint—an important connecting part—has been studied extensively, and the fatigue problem is no longer a challenge [7, 8]. Owing to the discontinuous stiffness at the crack, internal force redistribution and an uneven reinforcement force occur in the cracked composite beam. Therefore, fatigue problems tend to occur when the bridge deck operates with cracks [9, 10].

According to the existing research results, there is no relatively mature crack calculation theory for composite beams, and the test data about cracks of composite beams are also limited [11–20]. For example, Leonhardt et al. [11] performed fatigue tests on two composite beams with different strength connectors to study the slippage of welding nails in the negative-bending moment zone and the occurrence of cracks in concrete slabs. Zanuy et al. [14] obtained the calculation formula of the average strain of concrete cracks under fatigue load based on the tie-bar model and recommended the adoption of CEB-FIP and European standards. El-Zohairy et al. [20] showed that the degree of shear connection between steel beam and concrete deck played a major role in controlling longitudinal fatigue cracks of concrete deck.

In current research, structural fatigue assessment methods are mainly divided into two categories. The first category is assessment methods that do not calculate fatigue damage [21, 22]. This approach based on the maximum fatigue-load effect value obtained by loading the most unfavorable position is compared with the resistance value to assess the safety of the structure. The results obtained through this design method were conservative.

The second type of the assessment method is the one that requires calculating the degree of fatigue damage [22–26]. It can be divided into two subcategories. The first design method is stress history process without calculating the fatigue details, such as the simple fatigue assessment method in the European Code [25], the AASHTO method in the United States [26], and the bike method in the British BS5400 [22]. This kind of method is similar to assessment without the fatigue damage in calculating the fatigue-load effect. The most unfavorable load-effect amplitude is determined by applying the fatigue load to bridges, the number of stress cycles is determined according to the traffic volume, and then the fatigue damage is calculated according to the fatigue cumulative damage principle to obtain the fatigue life corresponding to fatigue details. The second design method is a design method that needs to calculate the stress history process of fatigue details, e.g., the linear cumulative damage method [25] in the European Code, the fatigue assessment guidelines for road bridges in Japan, and the vehicle load spectrum method [22] in the British code BS5400. In this method, the stress spectrum of fatigue details is obtained by loading the influence line according to the fatigue load, and the stress frequency spectrum is obtained through counting methods (such as the rain-flow method or water-discharge method). Then, the fatigue damage degree is calculated according to Miner’s linear cumulative damage principle, and a fatigue check is performed [27]. This loading mode reveals the process of fatigue damage and is more sophisticated than other methods.

The bridge deck directly bears vehicle loads; thus, its design affects the safety and durability of the overall bridge structure [28–30]. The position of the wheel load significantly influences the total strain/stress response, as well as the overall effect of the bridge deck. The local effect of the wheel load contributes significantly to the total strain/stress response, and the global effect may be negligible or significant, depending on the location [31, 32].

At present, research on the effect of the wheel position on the fatigue performance of bridge decks mainly focuses on steel decks [33–37]. For example, Zhou et al. [34] showed that the lateral position of a vehicle was very sensitive to the transverse stress of the steel bridge deck and the fatigue damage of key parts. Zeng et al. [35] showed that the cumulative effect of plastic deformation should be considered when analyzing the fatigue damage of the steel deck under a vehicle wheel load. Zhu et al. [33] studied the stress behavior and fatigue life estimation of the composite system of orthotropic steel deck (OSD) and ultrahigh performance concrete (UHPC) under concentrated wheel load by means of field monitoring and finite-element analysis. The results indicated that the stress value corresponding to the fatigue detail was sensitive to the transverse position of a single wheel load. Although in-depth studies have been conducted on the mechanical performance of concrete bridge decks [32, 38–40], the fatigue behavior of concrete slab of steel-concrete composite beam bridge under wheel load still needs further study.

Previous studies indicated that the reinforcement ratio of the bridge deck significantly affects the fatigue of the bridge deck [41–45]. Amir et al. [41] studied the effect of the reinforcement ratio on the transverse early cracking of a GFRP-RC bridge deck to evaluate the effect of noncorrosive glass fiber-reinforced polymer (GFRP) on the shrinkage performance of reinforced concrete (RC). Ghatefar et al. [42] conducted tests on a GFRP-RC bridge deck with a reinforcement ratio of 0.7% under freeze-thaw and dry-wet conditions. Su et al. [45] studied the effects of steel fiber-reinforced concrete panels with different reinforcement ratios on the negative-bending performance of orthotropic steel-concrete composite beams. However, no in-depth studies have been performed on the fatigue performance of steel-concrete composite beams.

Thus, in the present study, the fatigue performance of the deck of the Dazhenggang steel-concrete composite continuous beam bridge on the Shanghai-Hangzhou Expressway was investigated. First, the fatigue assessment process of the bridge deck was determined. Second, a finite-element analysis model was established, and the fatigue evaluation parameters of the bridge deck were determined according to traffic flow data of the expressway for 10 years. Finally, the effects of the wheel transverse position, steel and concrete damage, lane location, deck reinforcement ratio, and deck thickness on the fatigue performance of steelconcrete composite beams were analyzed.

#### 2. Fatigue Assessment Process of Bridge Deck

Various fatigue assessment methods are also applicable to punching fatigue of bridge decks [23, 24]. To accurately reflect the fatigue damage process, the method of calculating the fatigue detail stress history [46] was used to evaluate the fatigue of the bridge deck in this study, as shown in Figure 1. The steps for evaluating the fatigue damage of the bridge deck included the determination of the fatigue load, the calculation of the fatigue-load effect, and the calculation of the fatigue damage. The fatigue calculation of the steel bar was based on the fatigue vehicle, and the punching fatigue calculation of the concrete was based on the axle.

#### 3. Finite-Element Analysis

##### 3.1. Project Overview

As shown in Figure 2, the Shanghai-Hangzhou Expressway is a part of Shanghai-Kunming Expressway G60. It starts from Shanghai Xinzhuang; passes through Minhang, Songjiang, Jinshan in Shanghai, Jiashan, Jiaxing, Tongxiang, Haining, and Yuhang in Zhejiang; and finally connects Pengbu town in the eastern suburb of Hangzhou to Hangzhou-Ningbo Expressway.

The Dazhenggang steel-concrete composite beam bridge (2 × 75 m) of the Shanghai-Hangzhou Expressway is located between the Songjiang and Dayun toll stations. A total of 163,369,169 passenger cars and 65,998,076 trucks passed the Dayun toll station of the Shanghai-Hangzhou Expressway in the 10-year study period. To facilitate the calculation and analysis, the vehicles are classified, as shown in Table 1. Model 1 is a two-axle vehicle weighing <3 tons, model 2 is a two-axle vehicle weighing >3 tons, model 3 is a three-axle vehicle, model 4 is a four-axle vehicle, model 5 is a five-axle vehicle, and model 6 is a six-axle vehicle. As shown in Figure 3, two-axle vehicles weighing <3 tons were most common (78.1%), followed by two-axle vehicles weighing >3 tons (9.8%). Additionally, the number of five-axle large trucks was large (6.4%).

##### 3.2. Finite-Element Model

The bridge deck slab not only bears the local action of the vehicle but also participates in the force of the whole part as part of the composite beam. It is difficult to consider the effect of transverse beam on the supporting stiffness and shear lag of bridge deck when the simplified model is adopted, so it is impossible to accurately simulate the coupling effect of the first and second system models. The spatial finite-element model can reflect the actual stress state of the bridge [47]. The calculation model simulates the steel beam, concrete bridge deck, floor concrete, steel beam longitudinal and horizontal stiffener systems, and steel beam diaphragm systems.

The spatial force characteristics of the structure cannot be reflected accurately through the simplified method of effective distribution width [48, 49], and the finite-element method is used to simulate the spatial force characteristics of bridge in this study (Figure 4). The total of 132,503 units was divided in the finite-element model, and 22,634 spring units in each direction of *X*, *Y*, and *Z* to simulate bond slip were included, and 64,601 concrete units and 45268 shell unites was used to simulate the main beam and profiled steel plate. The steel bars were dispersed in the concrete, and both steel bars and concrete adopt linear elastic constitutive relations. The modeling method was consistent with references [50], and the correctness of the model is ensured.

#### 4. Fatigue Assessment Parameters of Bridge Deck

##### 4.1. Fatigue Load

According to a traffic investigation and statistical analysis of Jiangyin Bridge over the Yangtze River, the Third Bridge of the Nanjing Yangtse River, Humen Bridge, and the Second Nanjing Yangtze Bridge, the fatigue design loads were determined, as shown in Figure 5. There were 202,793,690 vehicles in the traffic volume survey data, 9,770,757 in the vehicle weight survey data, 30,431,682 axles in the axle load survey data, and 11057903 axles spacing in the wheelbase survey data.

**(a)**

**(b)**

To better reflect the characteristics of the vehicle load on the Shanghai-Hangzhou Expressway, the traffic load of the Dazhenggang continuous composite beam bridge was converted into the fatigue design load, as shown in Figure 5. According to Miner’s cumulative damage criterion [51],where is the equivalent number of vehicles under the fatigue design load, represents the vehicle weight of the fatigue design load model, which is taken as 330 kN, and represents the number of vehicles with weight .

Let *N* represent the total traffic flow, *P*_{ti} represent the proportions of different vehicle types, and *P*_{di} represent the probability of the axle load distribution. and can be determined using the following formulas:

According to the vehicle weight statistics for 474,660 vehicles that traveled on Jiangyin Bridge over the Yangtze River, the distribution and statistical characteristics of the vehicle weight for different vehicle types are presented in Figure 6 and Table 2. The equivalent number of vehicles corresponding to the fatigue design load for different vehicle types is presented in Table 3.

In this study, the distribution of vehicle types passing through the Dayun toll station is used to approximately replace the distribution of vehicle types passing through the Dazhenggang composite girder bridge in the Shanghai-Hangzhou Expressway. The vehicle distribution is as follows: Model 1, 78.1%; Model 2, 9.8%; Model 3, 2.5%; Model 4, 1.1%; Model 5, 6.4%; Model 6, 2.1%; and vehicles weighing >3 tons, 21.9%. Neglecting the effects of vehicles weighing <3 tons on the bridge fatigue, the equivalent vehicle conversion coefficient of the fatigue design load for vehicles weighing >3 tons can be obtained as follows:

##### 4.2. Axle Load Spectrum of Concrete Punching Fatigue Vehicle

For developing the axle load spectrum model, Jiangyin Bridge over Yangtze River, the Third Bridge of Nanjing Yangtse River, Humen Bridge, and the Second Nanjing Yangtze Bridge were taken as the research objects. The vehicle flow, vehicle axle load, wheelbase, and lane flow distribution were investigated and statistically analyzed. The axle load spectra for various vehicles (single axle and double wheel, single axle four-wheel, double axle, and triple-axle) were obtained, as shown in Figures 7–10.

The total axle load spectrum for all types of axles obtained from the statistical data is shown in Figure 11.

In this study, single-axle double-wheel, single axle four-wheel, double-axle, and three-axle units were converted to equivalent single-axle units according to the axle load. Models 1 and 2 are two-axle vehicles, while models 3–6 can be regarded as three-axle vehicles. According to the vehicle survey of the Dazhenggang Bridge, the number of axles (as shown in Table 4) can be determined. After the total vehicle axle number is determined, each axle load can be calculated according to the spectral value of the axle load.

#### 5. Analysis of Results

##### 5.1. Selection of Most Unfavorable Position of Lateral Force and Effect of Wheel Transverse Position

The fatigue failure of the RC bridge deck is mainly caused by the direct action of the wheels, and the stresses of the bridge deck at different transverse positions under the wheel load are not identical. To determine the most unfavorable position of the bridge deck, unit loads were arranged as close to each other as possible along the transverse direction of the bridge. The loading area was selected according to the provisions of the “General Design Specification for Highway Bridges” (JTGD60-015). A concentrated load of 100 kN was applied, and the spacing was 0.5 m along the transverse direction of the bridge. The stress of the concrete slab under the wheel load was calculated under loading, and the accuracy was sufficient to satisfy the calculation requirements. The transverse load distribution is shown in Figure 12.

The bridge deck at the top of the bearing in the negative-bending moment area was selected as the research object, as shown in Figure 13.

Vehicles do not drive exactly along the track line. Additionally, the distribution of the wheel track line along the transverse direction varies significantly among different countries; thus, it cannot be simply applied [52]. The Dazhenggang Bridge is a one-way four-lane bridge. According to statistics, the lateral swing range of centerline of the wheel track is approximately 1.5 m in the same lane. Considering the different stresses of the bridge deck caused by the wheel load in different transverse positions, the bridge deck 1.0 m from the centerline of the bridge was taken as the research object, and the horizontal spacing of 300 mm was taken as a working condition. This was divided into five working conditions: path 0, path –1, path 1, path –2, and path 2. The transverse position of the bridge is shown in Figure 14. The frequency diagram of the corresponding vehicle transverse distribution is shown in Figure 15, and the effect of the wheel transverse position on the stress amplitude of the bridge deck is shown in Figure 16.

The influence lines of different positions of the wheels, i.e., the influence surface of the bending moment of the concrete slab, were used for calculation. To compare the effects of various wheel lateral distribution models on the roof stress amplitude, we refer to the calculation method for the equivalent stress amplitude [53] and consider the equivalent stress amplitude of the vehicle lateral distribution according to Miner’s linear cumulative damage principle [54]. The equivalent calculation formula is as follows:where represents the equivalent force amplitude; represents the probability corresponding to the *i*^{th} wheel path (Figure 15); represents the stress amplitude corresponding to the *i*^{th} wheel path (Figure 16); and represents the slope of the S-N curve, which is taken as 3.

To better understand the differences between various wheel lateral distribution models and determine whether the effect of the bridge deck pavement on the lateral distribution effect of the wheels should be considered, the ratio of the equivalent stress amplitude calculated using the lateral distribution of the wheels to that calculated using the most unfavorable loading position was examined, as shown in Table 5. Because the performance of concrete bridge panels under a shared load is satisfactory and the swing of the wheel track has little effect on the stress of the bridge deck, the effect of the wheel track position was not considered in this study.

##### 5.2. Damage to Reinforcement and Concrete of Bridge Deck

###### 5.2.1. Checking Calculation of Fatigue of Reinforcement

The traffic volume used for the fatigue checking of the reinforcing bars was the predicted traffic volume in Table 4. There were four lanes in one direction. In the absence of lane distribution statistics, each lane was evenly allocated. During fatigue checking, it was assumed that vehicles weighing <3 tons did not cause fatigue damage to the steel bars. Therefore, these vehicles were not considered in the fatigue checking. A total of 21.9% of the vehicles in the traffic flow weighed >3 tons. According to the calculations and analysis presented in Section 3.1, each vehicle weighing >3 tons is equivalent to 0.948 fatigue vehicles.

The fatigue life of the bridge-deck reinforcement at different locations is presented in Table 6. According to calculations based on ECCS100 [25] and BS5400 [22], the fatigue life of the reinforcement at different locations is >100 years when the loading position is 1.0 m from the bridge centerline. According to ECCS100 [25], when the loading position is 7.0 m from the centerline of the bridge, the fatigue life of the cross-beam reinforcement is >100 years, while that of the midspan reinforcement is only 15 years. According to BS5400 [22], when the loading position is 7.0 m from the centerline of the bridge, the fatigue life of the cross-member reinforcement is 44 years and that of the midspan reinforcement is only 7 years. A comparative analysis revealed that the fatigue stress of the steel bar in the middle of the span is the most disadvantageous when the loading position is 7.0 m from the centerline of the bridge.

###### 5.2.2. Checking Calculation of Concrete Fatigue

For checking the punching fatigue of the concrete, the axle load spectrum used was the number of axle loads converted from the traffic volume predicted in Table 4. There were four lanes in one direction. At present, there are few provisions regarding the S-N curve of the punch-shear fatigue of concrete in national codes. It is generally considered that the punch-shear fatigue resistance is related to the ratio of the section shear force to the maximum shear capacity. The results of Matsui [55] were used in the present study.

Here, represents the punching shear and represents the punching fatigue resistance of the sections. According to formula (6),where represents the slab width, represents the shear strength of the concrete, represents the tensile strength of the concrete, represents the height of the compression zone, and represents the thickness of the protective layer of the steel bar. The tensile and compressive strengths of concrete are specified in the “Code for Design of Highway Reinforced Concrete and Prestressed Concrete Bridges and Culverts” (JTG 3362–2018). The shear strength was calculated using the formula developed by Guo and Xudong [56]:

According to the S-N curve of formula (6), the fatigue damage degree of the concrete at the center pivot position was calculated. The fatigue damage of the concrete was considered when the damage degree reached 1.0.

The fatigue life of the bridge-deck concrete at different locations is presented in Table 7. The fatigue life of the concrete at different locations was >100 years when the loading position was 1.0 m from the centerline of the bridge. When the loading position was 7.0 m from the centerline of the bridge, the fatigue life of the cross-beam concrete was >100 years, while that of the midspan concrete was only 37 years. Therefore, the fatigue stress of the concrete at the midspan was the most disadvantageous when the loading position was 7.0 m from the centerline of the bridge.

##### 5.3. Effect of Lane Position on Fatigue of Reinforcement Bars in Bridge Deck

According to the analysis presented in Section 4.2, the distance between the centerline of the lane and the most unfavorable fatigue position should be considered when deciding the layout of the lanes because of the different stress properties of the bridge deck in the transverse direction. Three different lane positions were selected to analyze the effect of the lane layout on the fatigue of bridge-deck reinforcement. Dazhenggang is a double-span bridge with four lanes in one direction. Half of the bridge structure was selected for the analysis. If the width of each lane is 3.75 m and the total width of lanes that can be arranged on the half-deck is 9 m, there are two lanes in each half-deck. Three different transverse positions were selected for the lane layout, as shown in Figures 17(a)–17(c).

**(a)**

**(b)**

**(c)**

The wheel positions and wheel distribution ranges for various working conditions are presented in Table 8.

According to the analysis presented in Section 4.2, the locations of the left wheel of Lane 1 and the right wheel of Lane 2 are close to the most unfavorable position of the lateral force. The bridge deck 1.0 m from the centerline of the bridge does not experience fatigue damage. Therefore, the bridge deck under the left wheel of Lane 1 was not analyzed; only the fatigue damage of the right wheel of Lane 2 was analyzed. The highest position of the wheel distribution frequency under various working conditions is presented in Table 9.

The fatigue damage at the most unfavorable transverse position of the bridge deck, that is, 7.0 m from the centerline of the bridge, and that at the highest frequency of wheels were compared. The results are presented in Table 10.

In working condition 1, the fatigue damage was high owing to the overlap of the lane centerline and the most disadvantageous position of the deck. However, in working conditions 2 and 3, the probability of driving in the most unfavorable position was reduced because the distance between the centerline of the lane and the most unfavorable position of the deck under stress was approximately 1.0 m, which significantly reduced the fatigue damage of the reinforcement. At the position 7.0 m from the centerline of the bridge, the damage degrees of conditions 2 and 3 were 56.2% and 61.8% of that of condition 1, respectively. The amount of damage at the most disadvantageous position under the static stress of the bridge deck was smaller than the amount of fatigue damage at the highest wheel frequency. Therefore, if the influence of lane position on fatigue is considered, the service life of bridge deck can be significantly improved.

##### 5.4. Effects of Reinforcement Ratio and Thickness of Bridge Deck on Bridge-Deck Fatigue

The tensile fatigue of the steel bar in the deck and the punching fatigue of the concrete are significantly affected by the thickness of the deck and reinforcement ratio of the steel bar [57]. According to the analysis results presented in Section 4.2, both the reinforcement and the concrete of the concrete deck 7.0 m from the centerline of the bridge will have suffered fatigue damage within 100 years. In this section, the bridge deck is taken as the research object. The S-N curve of BS5400 [22] is used for calculations, and the parameters shown in Table 11 are used. Three reinforcement ratios and fatigue damage degrees of the concrete and reinforcement under three slab thicknesses are selected for comparison, which provides a useful reference for the design of the bridge deck.

As shown in Figures 18–20 and Table 12, the fatigue life of the deck reinforcement significantly increased with an increase in the reinforcement ratio. At a slab thickness of <30 cm, even when the reinforcement ratio reached 3%, the fatigue life of the deck reinforcement was <50 years, which was significantly shorter than the 100-year design life of the bridge (Table 12). At a slab thickness of 35 cm, when the reinforcement ratio reached 3%, the fatigue life of the steel bar in the deck reached 101 years, which was longer than the design life of the bridge (Table 12).

The relationships between the punching fatigue life of the bridge deck and the reinforcement ratio and slab thickness are presented in Table 13. Increasing the reinforcement ratio was effective for restraining the punching fatigue damage to the bridge-deck concrete. When the slab thickness was 25 cm and the reinforcement ratio was 1.76%, the punching fatigue life of the bridge-deck concrete was only 5 years. Even when the reinforcement ratio reached 3%, the punching fatigue life of the deck concrete was only 34 years, which was significantly lower than the 100-year design life of the bridge (Table 13). When the slab thickness was ≥30 cm and the reinforcement ratio was ≥2.5%, the punching fatigue life of the bridge-deck concrete was >100 years and was longer than the bridge design life.

In conclusion, when the slab thickness of the bridge deck is <30 cm, the effects of increasing the reinforcement ratio on the fatigue life of the bridge-deck reinforcement and the punching shear fatigue life of the concrete are not obvious, and the increase in fatigue life is negligible in comparison with the 100-year design life of the bridge. For a 35 cm thick bridge deck, even a reinforcement ratio of 1.76% will not cause fatigue damage. Therefore, the concrete slab thickness of the bridge deck is recommended to be ≥35 cm, which complies with the strict requirement for the thickness of the concrete bridge deck in the Japanese Road specification [58].

#### 6. Conclusion

According to a load investigation, a fatigue assessment method for a composite beam bridge deck was proposed, and the fatigue damage and life for steel bar tension fatigue and concrete punching shear fatigue of a composite beam RC bridge deck were calculated. The following conclusions are drawn:(1)When the centerline of the wheel track swings laterally on the same lane for approximately 1.5 m, it has little influence on the force of the deck. The effect of the wheel track position can be ignored in the fatigue checking of the deck.(2)For the deck of the steel-concrete composite beam, the fatigue properties of the steel bar and concrete in the deck depend significantly on the transverse loading position. The fatigue life of the deck at different positions is >100 years when the loading position is 1.0 m from the centerline of the bridge. The fatigue life of the concrete across the middle span is only 37 years when the loading position is 7.0 m from the centerline of the bridge. Therefore, the fatigue stress of the concrete on the bridge deck is the most disadvantageous when the loading position is 7.0 m from the centerline of the bridge.(3)The lane arrangement significantly affects the fatigue damage of the bridge deck; thus, the overlap of the lane centerline and the most unfavorable stress position should be avoided to the greatest extent possible. The fatigue damage is reduced by 30%–40% when the centerline of the lane deviates from the most unfavorable stress position by 1 m.(4)The reinforcement ratio and thickness of the deck have similar effects on the fatigue life of deck reinforcement, and the punching fatigue of the concrete is more sensitive to changes in the thickness of the deck. When the deck thickness is <30 cm, the effects of increasing the reinforcement ratio on the fatigue life of the deck reinforcement and the punching fatigue life of the concrete are not obvious, and the increase in fatigue life is negligible in comparison with the 100-year design life of the bridge. For 35 cm thick decks, even a reinforcement ratio of 1.76% will not cause fatigue damage. The thickness of the concrete slab of the bridge deck is recommended to be ≥35 cm.

#### Data Availability

The data used to support the findings of this study are not available.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 51878623), Innovative Research Team (in Science and Technology) in University of Henan Province (no. 20IRTSTHN009), Program for Young Backbone Teachers in Colleges and Universities in Henan (no. 2018GGJS005), and Foundation for Postdoctoral Students in Henan Province (no. 1901024).