Abstract

The study aims to research the influence of salt fog corrosion cycles on seismic performance of reinforced concrete (RC) frame beam-column joints in coastal atmosphere. Based on low cyclic loading tests of six RC frame beam-column joint specimens, this study analyses the failure patterns, hysteresis loops, load carrying capacity, displacement, backbone curves, and energy dissipation capacity of corrosion-damaged RC frame beam-column joints. The effect of salt fog corrosion cycles and axial compression ratios are tested repeatedly. The results show that with the same level of axial compression of the frame joint specimens, as the increase of salt fog cycles, the strength, ductility, energy dissipation, bearing capacity, and deformation capacity of joints degenerated to different degrees. When the corrosion level is the same, the stiffness degradation appeared to be more apparent as the increase of axial compression ratio. Then, the behavior degeneration rule of the RC frame beam-column joints is analyzed and formed according to the results of the test; the degeneration restoring force models of corroded RC frame beam-column joints is formed and verified based on Clough’s three-line degenerate restoring force model and the introduction to cyclic degeneration index. The results show that the restoring force model can better describe the hysteresis characteristics of the beam-column joints of corroded RC frames. The research is a theoretical reference for the seismic analysis of the RC frame structure affected by coastal atmospheric environment.

1. Introduction

Reinforcement corrosion is the main cause of deterioration of RC structures [1, 2]. Longitudinal reinforcement corrosion will result in the decrease of effective cross-sectional area, the degradation of mechanical properties, and weakening the bond performance between the steel bar and the concrete [3]. Stirrups corrosion will reduce the restraint on the core concrete [4], which will deteriorate the performance of RC structures [5]. Because of the higher level of air humidity and chloride ion, the steel cumulative mass loss rate in coastal atmosphere is faster than that in the general atmosphere. As a result, chloride ion erosion is the main cause of steel bars corrosion in coastal atmosphere [6]. Moreover, China is a country with a long coastline and frequent earthquakes, and RC structures in coastal atmosphere are located in the earthquake regions as well. Therefore, it is necessary to study the effect of chloride ion ingress on the seismic performance of RC structures in coastal atmosphere.

Currently, domestic and foreign scholars have conducted extensive researches on the deterioration of the seismic performance of corroded RC structures. Ou and Nguyen [7, 8] have experimented the seismic performance of corroded reinforced concrete beams. It showed that corrosion in longitudinal tension reinforcement had significant detrimental effects on the yield drift, yield load, and peak load [7]. The corrosion of hoops adversely affected the deformation capacity of the beams significantly yet did not significantly influence the load-carrying capacity of the beams [8]. Zheng et al. [9] experimented on the seismic behaviors of corroded RC frame beams, which failed in flexure-shear mode under the artificial climate. It showed that the ductility, bearing capacity, stiffness, and hysteretic energy of the specimens gradually decreased with the increase of the corrosion degree of steel bars. The failure mode would transit from flexure-shear failure to shear failure. Ou et al. [10] studied the performance of a large proportion of corroded RC beams subjected to cyclic loading. The experimental results were similar to those of Zheng et al. [9]. A group of studies have researched the influence of the corrosion degrees and the axial compression ratios on the seismic performance of RC columns [1114]. Zheng et al. [11] indicated that the bearing capacity, deformation capacity, and energy dissipation capacity of the specimens decreased as the steel corrosion increases.

Gan et al. [15] proposed a joint system that connected the thin-walled square reinforced concrete–filled steel tube (RCFST) column to RC beams at both sides of the joint. Zhou et al. [16] conducted the cyclic test of square-tubed reinforced concrete column to the RC beam joints and studied the seismic performance of the joints. Kanchanadevi and Ramanjaneyulu [17] investigated the seismic response of corrosion-damaged beam-column subassemblages. The results showed that the corroded specimens displayed poor hysteretic performance in the form of huge in-cyclic strength degradation. Yang et al. [18] conducted the cyclic tests of ten corroded beam-column joints of SCOPE system, and the studies indicated that corrosion of steel bars significantly affected seismic behaviors of the structure. From the above studies, the research on seismic performance on corroded concrete structures mainly focused on flexural and compressive members, and the RC frame joint specimens were still less. It should be noted that most research works used electrochemical methods to accelerate the steel corrosion [7, 8, 1214, 19, 20], while Yuan et al. [21] and Zhang et al. [22] indicated that the corrosion mechanism and the surface characteristics of steel bars after corrosion under the condition of electricity were obviously different from those under natural conditions, but the artificial weather environment to accelerate corrosion was basically the same as natural environment corrosion. Artificial climate environment simulation technology is an important way to develop durability testing methods for RC structural members [21].

This study adopts the artificial climate environment accelerating corrosion technology (salt spray corrosion test) to simulate the coastal atmospheric environment to conduct the accelerated corrosion test of six RC frame beam-column joint specimens and the pseudo-static tests in order to research the influence of the steel corrosion degree and the axial compression ratio n (, where N is the axial force, is the concrete compressive strength, and and are the concrete column section width and height, respectively) on the seismic behavior of the RC frame beam-column joints.

In recent years, a lot of research works have been carried out on the restoring force models of structures at home and abroad. The typical restoring force models include Clough and Johnston degenerate bilinear model [23], Ramberg–Osgood restoring force model [24], Takeda degenerate three-line model [25], Bouc and Wen restoring force model [26, 27], etc. However, these studies are based on the intact specimens and do not consider the impact of environmental factors, for instance, the study on the influence of the salt-fog corrosion in coastal atmosphere on the seismic performance of RC-frame structures is less, and the study on the corresponding restoring force model is almost blank.

The purpose of this study is to research the seismic behavior of reinforced concrete frame joints in coastal atmosphere, and the backbone curve model of corrosion joints and the cyclic degradation index based on salt fog damage are obtained by the seismic experimental results. Then the restoring force model of joint considering the corrosion damage of steel bar is acquired and verified. The study can provide a basis for evaluating the seismic performance of the existing concrete structures in coastal atmosphere.

2. Experimental Program

2.1. Design and Fabrication of Specimens

Based on the design principle of “strong component, weak joint,” this experiment designed and fabricated six RC frame beam-column joint specimens, which were taken from the beam-column assembly between the conventional frame anti-bending points. Considering the limitation of the dimensions of the environmental chamber used in this study, the geometric reduction ratio of the joint specimen was 1 : 3. The main research parameters were the axial compression ratio and the degree of steel corrosion, and the design parameters of joint specimens are shown in Table 1. The joint specimens were 2.2 m long, 1.45 m high, and the cross-section of the beam was 150 × 250 mm2, the cross-section of the column was 200 × 200 mm2, the concrete strength grade was C30, the longitudinal bars of the beam-column members were all HRB335, and the stirrups were all HPB235. Cross-section and reinforcement details of the six joint specimens are shown in Figure 1.

The mixing ratio for C30 concrete was 320 kg/m³ of cement, 870 kg/m³ of sand, 870 kg/m³ of fine gravel, and135 kg/m³ of water. Ordinary Portland cement P.O 32.5R and medium sand were used. Three standard cubic specimens of concrete were prepared under the same conditions as the RC beam-column joint specimens. The real strengths of the concrete and the reinforcement bars were determined via material property tests. The average cubic compressive strength of the concrete was 27.03 MPa, and the properties of the reinforcement bars are shown in Table 2.

2.2. Specimen Corrosion Program

Neutral salt spray (NSS) tests [28] were carried out in the large-scale climate simulation laboratory at the Xi’an University of Architecture and Technology (Figure 2). To accelerate corrosion, 5% of NaCl was added to the specimens when pouring concrete except 0.3-No. After 28 days of natural conservation, the specimens were moved into the artificial climate laboratory for rapid salt spray corrosion. In order to make the indoor salt fog concentration constant, intermittent spray mist spray was used. In a one-hour cycle, the salt spray was stopped for 40 minutes after spraying for 20 minutes. This process was carried out for a total of 3 hours, and then the temperature in the climate chamber was raised to 60°C ± 2°C for drying for 2 hours. Afterwards, the temperature was reduced to 45°C at a rate of 0.5°C per minute, and the room humidity was controlled at about 90%. The salt fog drying conversion time was set to 30 minutes. Therefore, a salt spray corrosion cycle takes 6 hours and 4 cycles per day. The schematic diagram of the wet-dry cycles is shown in Figure 3.

2.3. Instrumentation and Structural Test Procedures

The cyclic loading test setup and measurement configuration are shown in Figure 4. The target vertical load was exerted on the specimens via hydraulic jacks and a pressure sensor was placed between the top of the column and the jack. At the same time, two load sensors were arranged at the beam ends to measure the tensile and compressive loads generated at the beam ends during the loading process. Linear variable displacement transducers (LVDTs) were installed to monitor lateral displacements. The lateral displacements at the loading point were measured by a horizontal LVDT, which had a height of 1650 mm to the upper surface of the strong floor. In order to prevent local compression failure of the concrete during loading, steel pads were provided at the loading point and the support position.

According to “specification of testing methods for earthquake resistant building” [30], a force-displacement hybrid control system was adopted. A load-controlled step before longitudinal reinforcement yield and a displacement-controlled after longitudinal reinforcement yield were adopted in the lateral loading procedure, as shown in Figure 5. The MTS actuator stopped when residual lateral load decreased to 85% of the peak load.

3. Experimental Results and Discussion

3.1. Corrosion Phenomenon

(1)The artificial climate corrosion process is complicated and the control method for the corrosion degree of the components is still not perfected. Jin et al. [31] pointed out that the corrosion crack width can be used to characterize the corrosion degree of the component. Therefore, a crack observer with an accuracy of 0.01 mm and a measuring range of 0 to 10 mm was used to periodically observe the rust expansion crack on the surface of the test pieces to control the crack width. To reduce the dispersion, the cracks formed in the original curing process were removed, and multiple main cracks were recorded and averaged. When the specimens’ rust expansion cracks reached the expected design width, the test pieces were taken out of the artificial climate chamber in batches.(2)The concrete of specimens with different degrees of rust were partially damaged, and then, several sections of stirrups and longitudinal reinforcements were taken out at the ends and middle portions of the test specimens, respectively. The corrosion state was observed, the cumulative mass loss rate was determined, and the material properties (take the average values) were tested. The results showed that the corrosion degrees of different specimens were not the same, and the differences in the corrosion degrees of stirrups were obvious. The corrosion degrees of longitudinal bars in beams and columns were relatively light, but the corrosion of the longitudinal reinforcements in columns were more serious than beams. The corrosion degrees of the stirrups in the core areas were significantly greater than that of the longitudinal in beams and columns.(3)Took stirrups and beam and column longitudinal reinforcements on both sides 3 for each with a length of 30 cm, after numbering, scraped the corrosion products, washed the rusted products in 12% hydrochloric acid solution, then rinsed in clean water, and neutralized in lime water, finally rinsed in water and dry, and recorded the weight of the uncorroded part. In order to reduce the discrete type of test, the cumulative mass loss rates of steel bars were averaged to the nearest 0.01% and calculated according to equation (1). The calculated values are shown in Table 3.Here, is the cumulative mass loss rate of steel bars; and are the weight of the un-corroded and corroded steel bars, respectively.(4)Mechanical properties of rusted steel bars were tested, and yield strength, ultimate strength, and elastic modulus were obtained, as shown in Table 4.

From Tables 3 and 4, it can be seen that the mechanical properties of steel bars with different corrosion degrees have been reduced to varying degrees. The experiment showed that under the salt spray corrosion, pitting corrosion of steel bars occurred, and the effective section areas reduced. When the reinforcements were stressed, the stress was highly concentrated in the pitted regions, which caused the early destruction of the steel bars.

3.2. Cracking and Failure Patterns of RC Frame Joint Specimens

Figure 6 shows the final failure modes of the specimens. The failure process of all joint specimens was the beam ends concrete cracked first and then yielded, after that, with the further increase of horizontal load, the core concrete was rapidly cracked, and the stirrups in the core area yielded until failed. However, with the increase of the design corrosion cracks widths (such as 0.3-No, 0.3-Sl, 0.3-Mo, 0.3-Se), the earlier the cracking and yielding of concrete at the beam ends occurred, the relatively rapid crack development of the core concrete, and the ductility of the specimens was reduced. When the specimens have the same crack widths and different axial compression ratios (such as 0.1-Mo, 0.45-Mo), the larger the axial compression ratio, the higher the bearing capacity.

3.3. Hysteretic Loops

Figure 7 shows the load-lateral displacement hysteresis of joint specimens under low cyclic loads. Observations made based on Figure 7 are summarized in the following:(1)In the virgin loading, the specimens were in the elastic stage. There was a linear relationship between load and displacement, and the paths of loading and unloading were basically in coincidence. As the load amplitude increased, the slope of hysteresis loops decreased gradually; stiffness degradation during the reloading paths and the residual deformation at the end of unloading paths began to appear, and specimens were in the elastic-plastic stage. After reaching the peak value, the strength and stiffness degradation became more and more obvious with the increase of lateral displacement, and the hysteresis loops exhibited obvious “pinch-in” phenomenon.(2)There were big differences in the hysteresis loops’ shape between heavily rusted specimens and un-corroded. The concrete manifestation was that as the corrosion degree increase, the peak value and slope of the hysteresis loops gradually decreased, the plumpness of the hysteresis curves gradually decreased, and the pinching phenomenon became more and more obvious, which indicated that the energy consumption of the specimens was gradually reduced.(3)For specimens 0.3-Mo, 0.1-Mo, and 0.45-Mo, it can be seen from the hysteresis curves that if the axial pressure of the joint is relatively large (e.g., 0.45-Mo), its ultimate displacement is decreased and the ductility is reduced.

3.4. Pseudo-Static Test Results

The load and displacement values of the specimens are shown in Table 5. Based on Table 5, the following observations and interpretations can be made: (1) when the axial compression ratio was the same, with the increase of the stirrup cumulative mass loss rate, the cracking, yielding, and peak loads were decreased. (2) As the cumulative mass loss rate of stirrups increased, the yield displacements of specimens 0.3-Sl, 0.3-Mo, and 0.3-Se changed little compared with that of undamaged specimen 0.3-No, while the peak and ultimate displacements basically showed the decreasing trend. (3) When the corrosion degrees were similar, the cracking, yielding, and peak loads of specimens 0.1-Mo, 0.3-Mo, and 0.45-Mo increased gradually with the increase of axial compression ratio, but the increase was small. (4) When the corrosion degrees were similar, the variation of yielding, peak, and ultimate displacements of specimens 0.1-Mo, 0.3-Mo, and 0.45-Mo were not obvious with the increase of axial compression ratio.

The backbone curves of specimens are shown in Figure 8. As Figure 8(a) shows, the yielding, peak, and ultimate load of the corrosion specimens were all lower than the intact specimen. Before yielding, the backbone curves basically coincided and the rigidity changed little; after yielding, the straight section of the backbone curves began to shorten; and after the peak load, the descending stage gradually became steeper. With the increase of the corrosion degree, the bearing capacity and ductility of the specimens were gradually degraded. As Figure 8(b) shows, the higher the axial compression ratio is, the higher the peak value of the backbone curve is; however, the smaller the corresponding displacement, the steeper the descending branch of the backbone curve become. It indicated that for the salt spray damaged RC frame joints, the load carrying capacity increased, whereas the deformability capacity decreased and the strength decay become more abrupt with the increase of axial compression.

3.5. Hysteretic Dissipation Energy

The energy dissipation capacity of RC frame specimens was analyzed in this study by using the equivalent viscosity damping coefficient and work index [32, 33].

3.5.1. Equivalent Viscous Damping Coefficient

The equivalent viscous damping coefficient reflects the energy dissipation capacity of the structure or component under an loading and unloading cycle, which can be calculated from.where is the area of enclosed hysteresis loop ABCD and and are areas of triangles OBE and ODF, respectively, as shown in Figure 9.

The equivalent viscous damping coefficients and , which corresponds to the maximum load and the ultimate load, respectively, are listed in Table 6. It could be seen that compared with un-corroded specimen, the equivalent viscous coefficients of the corrosion specimens in the limit state decrease with the increase of the cumulative mass loss rate, which demonstrates that the energy dissipation capacity of the corroded components are reduced in the limit state.

3.5.2. Work Index

The work index defined by Gosain [33] can be calculated as follows:where is the peak load in the ith cycle, is the displacement corresponding to , and and are the yielding load and displacement, respectively. The calculation results are shown in Table 6. It could be seen that as the cumulative mass loss rate increased, both the equivalent viscosity coefficient and the work index gradually decreased, which indicated that the energy dissipation capacity of the specimens were reduced. So corrosion could degrade the energy consumption of components.

4. Restoring Force Model Establishment

To study the influence of salt spray corrosion on the deformation performance of the RC frame joints and the entire assembly, this study first calculated the shear deformation in the joint core region based on the experimental data of the diagonal length, and then the corresponding shear force was calculated to establish the shear-shear deformation relationship curve. Meanwhile, it can be seen that the restoring force model of the corroded specimen and the uncorroded specimen were similar in shape, but as the reinforcement corrosion severely weakened the bond performance between the steel bars and the concretes, the stress between them was changed, as well as some feature points of the backbone curves. So after analyzing the experimental data and taking into account the restoring force model of the three-fold line shear deformation, the backbone curve model of the corrode RC frame joint was established, as shown in Figure 10. The backbone curve needs to determine five characteristic points, namely, the cracking shear , the peak shear , the elastic shear stiffness , the hardening shear stiffness , and the softening shear stiffness , respectively.

4.1. The Feature Points of Uncorroded Components
4.1.1. The Cracking Shear

The cracking shear can be calculated as follows [34]:where is the comprehensive coefficient of influence that takes the value 0.67 and is the constraint coefficient of the beam to the joint. For the joints without the orthogonal beams,  = 1.0; and are the effective width and height of the joint core section, respectively. is the stirrup stress development coefficient, which takes the value 0.1; is the yield strength of stirrups in the joint; is the sectional area of all stirrups in the same section; is the spacing of stirrups; is the effective height of the beam section; is the distance from the resultant center of the tensile reinforcement of the beam to the edge of the nearest section; and when , .

4.1.2. The Peak Shear

It was often considered that the joint reached the split state as its failure criterion, the shear strength of the joint was known as the peak value at this time. After the concrete in the joint core areas were cracked, the shear force was assumed by the horizontal stirrups in the core areas and the longitudinal reinforcements on the front and back sides of the column, thereby which form a “truss mechanism.” Reference on the current specification [34], the peak shear of the joint could be determined as follows:where is the constraint coefficient of the orthogonal beam to the joint, is the design value of the axial force at the bottom of the column at the joint considering the design value of seismic shear force. The above value can refer to the specification [34]. At the same time, the peak shear force must satisfy the condition:

4.1.3. Calculation of Shear Stiffness of Joints under Different Stress Stages

(1)Elastic Shear Stiffness at Segment 01where G is the elastic shear modulus and is calculated by: , is the elastic modulus of the joint concrete. A is the shear area of the joint core area, is Poisson ratio (for concrete materials, is 0.2).(2)Hardened Shear Rigidity of Segment 12At this stage, the joints are in an elastic-plastic state and their stiffness are degraded. The expression of is as follows:where is the stiffness degradation coefficient of the joint, the axial compression ratio, and the shear pressure ratio of the joint have an effect on the hardening stiffness. Based on references [3537], the expression of is as follows:where is the shear ratio of the joint , n is the axial compression ratio of the column, .(3)Softening Shear Rigidity of Segment 23According to Bazant [38], the ultimate shear capacity of a joint was about 1.2 times for its cracking capacity. Considering the importance of the joint, it is conservatively assumed that the joint would enter the ideal elastic-plastic state after it reached the cleft. In the actual projects, the designs are based on the principle of “strong joint,” and it is impossible for the joint to reach the limit state before the structure collapsed. Hence,

From the stiffness calculation formulas above, the shear deformation and corresponding to the crack state and the peak state could be calculated.

4.2. Calculation of Characteristic Points of Corroded Specimens

According to the deformations of the joint core areas measured in the tests, the horizontal shear force and shear deformation in the core areas of specimens were calculated. The fitting analysis is performed to obtain the degradation of characteristic parameters of the restoring force model of corroded specimens with the variation of the stirrup corrosion ratios and the axial compression ratios. The formulas calculated the characteristic parameters considering the effects of rust damage are shown in Table 7. At the same time, fitting plots of equations (11)–(14) are shown in Figure 11.

Note: In Table 7, and are the cracking shear and peak shear of the rusted specimens, respectively; and are the cracking shear deformation and the peak shear deformation of the rusted components, respectively; ,, and are the elastic stiffness, hardening stiffness, and ultimate stiffness of the rusted component, respectively; is the average cumulative mass loss rate of stirrups in the joint area; is axial compression ratio.

From the correlation coefficient in Table 7 and the fitting plots in Figure 11, it could be seen that the fitted values are similar to the actual values. Using the fitting formulas and the feature point values of the good components, the characteristic point values of the specimens with different degrees of corrosion can be obtained. The results are shown in Table 8.

4.3. Establishment of Hysteresis Rules

From the hysteretic curves of specimens, we can conclude that the strength and stiffness of specimens were deteriorated during each cyclic loading with the increase of corrosion degrees. Ibarra et al. [39] researched hysteretic models that incorporated strength and stiffness deterioration, and in this research, three of the basic hysteretic models used in seismic demand evaluation were modified to include deterioration properties: bilinear, peak-oriented, and pinching. In this, the peak-oriented model incorporate an energy-based deterioration parameter [40] that controls four cyclic deterioration modes: basic strength, post-capping strength, unloading stiffness, and reloading stiffness deterioration. Four cyclic deterioration modes had been defined by the equations and interpreted with the figures in details. Therefore, referring to Ibarra et al. [39], this study chooses the peak-oriented model as the hysteretic rules of corrosion specimens, the strength and stiffness deterioration modes are shown in Figure 12, and the equations of individual cyclic deterioration modes are listed as shown in Table 9.

Note: In Table 9, is cyclic deterioration parameter [40]; and are the cracking shear force after and before excursion i, respectively; and are the hardening stiffness after and before excursion i, respectively; and are the peak shear after and before excursion i, respectively; and are the unloading stiffness after and before excursion i, respectively; “” “represents the loading direction; “+” represents the positive loading; and “” represents the negative loading.

5. Restoration Force Model Verification

To verify the practicability of the restoration force model proposed in this study, taking the specimens 0.3-Se and 0.1-Mo as examples, the calculated hysteresis curve is compared with the test. The comparison results are shown in Figure 13.

As seen from Figure 13, the model curve and the test curve are close to each other in the positive direction, but in the negative direction, there is a large deviation. The probable cause is that corrosion of reinforced concrete specimens is a very complicated process, and the steel reinforcements are likely to be unevenly rusted, but when establishing the theoretical model, it is assumed that the specimens are uniformly rusted, so there is a difference between the experimental model and the theoretical model. 0.3-Se and 0.1-Mo represent the severe and moderate corrosion of the specimens, respectively. The calculated results are close to the test results. The pinching effect is caused by the crack development, the bond-slip of the longitudinal beams in the loading process reaches to a certain degree. This shows that the restoring force model established has a good applicability, which can reflect the degradation of the performance of corroded components during loading truly.

6. Conclusions

Six RC frame beam-column joint specimens were subjected to rapid artificial climate corrosion test in the coastal environment, and then they were tested under lateral cyclic loads. The effects of cumulative mass loss rates and axial compression ratios on strength, ductility and energy dissipation capacity of beam-column joints were analyzed, and the restoring force model of corroded RC frame beam-column joints was studied. The following conclusions can be drawn from the test results:(1)For the RC frame beam-column joints, the corrosion of the stirrups was more severe than the longitudinal reinforcements. Corroded stirrup weakened the restraint effect of the concrete in the core area of the joints, and the destruction of the joint area was aggravated.(2)With the increase of cumulative mass loss rate of steel bars, the strength, ductility, and energy dissipation capacity of RC frame joint components were all degraded to varying degrees, manifested as the reduction of bearing capacity and deformation capacity of the components. When the corrosion degree was the same, the stiffness of the joint specimens decreased more significantly with the increase of the axial compression ratio. The stiffness of specimens with a relatively small axial pressure was significantly degraded before cracking and tended to be stable after yielding.(3)The restoration force model of RC frame beam-column joints based on the tri-line restoring force model and considering the salt spray corrosion in coastal atmosphere accurately described the cyclic degradation and pinch effect caused by cyclic loading and corrosion damage.(4)The calculated results of the restoring force model of corroded RC frame beam-column joints are in good agreement with the experimental results. It shows that the restoring force model would better reflect the hysteresis characteristics of the RC frame beam-column joints in coastal environment than un-corroded the RC frame beam-column joints, which can provide a theoretical reference for the elastic-plastic seismic response analysis of RC frame structures in the coastal environment.

Data Availability

All data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors are grateful for the financial support received from the National Natural Science Foundation of China (no. 51678475), National Key R&D Project (no. SQ2019YFC150111-02), Industrialization Project of Shaanxi Education Department (no. 18JC020), and Xi’an Science and Technology Plan Project (no. 2019113813CXSF016SF026). The authors would like to thank the reviewers’ efforts in reading the manuscript.