Abstract

This study is to systematically analyze the influences of time-varying meshing stiffness (TVMS) and meshing impact on the dynamic characteristics of high-speed gear transmission in the two-stage pure electric vehicle (PEV) gearbox, as well as the effect of tooth surface modification on the vibration control. First, the dynamic model was established, and the TVMS and meshing impact were calculated. Then, the vibration characteristics of single-stage and two-stage helical gear transmission were analyzed under three different excitation conditions, excitation of TVMS, excitation of meshing impact, and excitation of both. The results show that the effect of rotating speed on the system vibration is not significant outside the resonant region under the excitation of TVMS, while the effect of meshing impact becomes the main exciting component with the increasing rotating speed. The vibrations of the two gear pairs interact with each other; the vibration frequency of one gear pair contains both its meshing frequency and the coupling frequency of the other gear pair. Tooth surface modification in the input-stage gear pair can reduce the vibration of both the input- and the output-stage obviously; that is, more attention should be paid to the input-stage gear pair in the modification design of PEV gearbox.

1. Introduction

With the aggravation of the energy and environment issue, the PEVs have shown a broad application prospect and achieved rapid development because of the advantages in energy conservation and environmental protection, as well as the performance and the potential for intelligent driving [13]. In the PEV power transmission system, the combination of electric motor and two-stage gearbox is widely used [4, 5], such as in Tesla Model S and NIO ES8. For increasing the power density ratio, the electric motor speed is relatively high. Thus, the gear pair of the input-stage often works in the high-speed off-resonance region. High rotating speed will necessarily bring a stern challenge to the dynamic performance of the gear transmission in PEV gearbox, which has a significant effect on vehicle reliability, stability, and vibration noise.

The dynamic characteristics of gear transmission system have always been a hot topic of scholars [68]. There are many excitation factors affecting the dynamic characteristics of gear transmissions. Among them, the meshing stiffness [9], errors [10], and the meshing impact [11] are the main internal excitations, and the frictional excitations along the line of action direction are the secondary excitations, when the aforementioned two kinds of excitations are minimized, the shuttling excitation along the axial direction may be dominant [12].

In early researches, TVMS was replaced by the mean meshing stiffness, while the accurate methods for calculating the TVMS have been proposed in recent years [13], and the influence of TVMS on the gear transmission has attracted wide attention. Saxena et al. calculated the gear meshing stiffness of misaligned shaft based on the potential energy method and studied the effect of misalignment and friction on the cracked gear pair [14]. Later, a calculation method for meshing stiffness of helical gear considering the actual pitch deviation was proposed, and the influence mechanism of pitch deviation on the meshing stiffness and vibration noise of helical gears were analyzed by Yin et al. [15], and another accurate calculation method for the helical gear meshing stiffness was proposed by Tang et al. [16]. Han and Qi also proposed a TVMS of helical gears by slicing and discrete integral method, and the effects of gear tooth defects, such as cracks, pitting, spalling, and breakage, on TVMS were analyzed as well [17]. In a recent research, Sánchez et al. proposed the meshing stiffness equations including Hertzian effects for spur gears, as well as the approximate equations for the meshing stiffness and the load sharing ratio [18]. For the meshing impact, Chen et al. calculated the impact force of the sun gear in a planetary gear system under different loads, and the meshing stiffness and the dynamic transmission error of integrated planetary gear were investigated when the planetary gear meshed with the sun or ring gear [19]. In Wang’s work, a dynamic model of herringbone gear including TVMS, meshing impact, and tooth surface friction was established and the vibration characteristics of the system were analyzed [20]. Zhou and Chen also proposed a meshing impact calculation method considering tooth surface friction [21].

Most of the above studies, however, are about the medium and low rotating speed conditions for single-stage gear transmission. On the one hand, as mentioned above, the input rotating speed of most of the gear transmissions in PEV field has reached more than 10000 rpm. The meshing impact has become a non-negligible internal excitation factor. On the other hand, the effects of excitation transfer and superposition between multistage gear pairs on the dynamic performance of the overall system are not negligible.

This paper focuses on finding out the effects of TVMS and meshing impact on the dynamic characteristics of high-speed helical gear transmissions in PEV gearbox, as well as the effect of tooth surface modification on the vibration control. A two-stage high-speed helical gear transmission system in a PEV reducer is taken as the research object in this study. First, the bending-torsion-axis dynamic model is established. Next, the TVMS, the impact time, and the amplitude of the impact force are calculated. Afterward, the dynamic characteristics of the single-stage and two-stage helical gear systems are analyzed under three kinds of excitation conditions, and the dynamic characteristics of the single-stage are compared to the two-stage gear systems as well. Finally, the effects of two modification schemes on the dynamic characteristics are studied, and some concluding remarks regarding this research are provided.

2. Dynamics Model

The 3D model and basic parameters of the two-stage high-speed helical gear transmission in the PEV reducer are shown in Figure 1 and Table 1, respectively.

In two-stage helical gear transmission, the angle of the two planes of action αm can be calculated by equation (1) according to Figure 2:where indicate the input and load torque, the tangents A1A2 and B1B2 represent the planes of action for each gear pair, the base radii of each gear are , the working pressure angles are α1 and α2, and the angle of two axis planes is γa. For the two-stage gear transmission given in Table 1, the calculation result of αm is 0°. That is, the two planes of action are parallel. Therefore, the following dynamic model is established by the lumped-mass method, as shown in Figure 3. and are the mass and the inertia moment of the gears; the values of the equivalent support stiffness and damping , as well as the torsional stiffness and damping of the intermediate shaft, are given in Table 2.

Then, the differential equations of this helical gear transmission system are derived:where , indicate the vibration displacements of the input shaft, the vibration displacements of intermediate and output shaft are , and , , the torsional displacements of the gears are , , , , the relative angle of gear 1 and pinion 2 is , and the meshing impacts of two gear pairs are represented by . The displacement matrix of this system can be expressed as

To compare the dynamic characteristics of the single-stage gear system with the two-stage gear system, the differential equations of the input gear pair are as follows:

The tangential and axial dynamic meshing forces of the system are expressed inwhere indicates the TVMS; the helix angles and pressure angle are given in Table 1.

The formula for calculating the meshing damping is shown inwhere is the damping ratio.

To eliminate the influence of the angular displacement of the gear rigid body on the subsequent analysis, the relative angular displacement is introduced to convert the angular displacement into linear displacement, as shown in

The subscripts 1 and 2 represent the input and output gear pair, respectively. Substituting equation (7) into equations (4) or (3),where and are the equivalent torsional qualities of the gear pairs:

3. Calculation of the TVMS

The TVMS is calculated based on the improved loaded tooth contact analysis (LTCA) method [22]. This LTCA method combines the gear geometry analysis with the mechanical analysis organically, and many factors such as actual coincidence degree and error are considered comprehensively [23], and the tooth surface load distribution as well as the loaded transmission error (LTE) under the current contact position load can be obtained. The LTE is then converted into the meshing normal deformation according towhere , , and represent the normal modulus, the number of driven teeth, and the end pressure angle, respectively.

The TVMS is then obtained:

The LTE of the gear pair is mainly caused by the geometric transmission error (TE), the gear tooth bending deformation, and tooth surface contact deformation. The TE is determined by tooth surface design and machining. When the material, the geometric parameters, and the meshing position are confirmed, the tooth bending deformation and the contact deformation are determined by the load, which can be obtained according to the following equation [24]:where represent the underdetermined coefficients, while represent the geometric transmission error, the tooth bending deformation, and the contact deformation.

The meshing period is divided into parts, and represents the deformation of the k-th () meshing point under the nominal load , which can be obtained by the gear LTCA program.

Thus, the load deformation of the k-th meshing position under the nominal load can be expressed as

The LTE of meshing positions under a given load can be obtained by LTCA, and the meshing stiffness of each point is calculated by equations (10) and (11). By the curve fitting method, the TVMS function is then derived under a certain load.

The LTEs of a meshing point under three different loads can be obtained by the same method. And then the coefficients can be determined; the relationship between load and deformation of the gear pair under different loads can be obtained. The meshing stiffness of a certain meshing point under any load can then be obtained. Finally, the TVMS can be obtained under any load conditions by computations on each meshing point in the whole meshing period.

According to the previous method, the TVMS curve of the input-stage and output-stage gear pairs of the PEV reducer under torque is calculated and shown in Figure 4:

4. Calculation and Analysis of the Meshing Impact

The meshing points of the ideal gear pair (without error deformation) are on the theoretical meshing line. Due to the meshing synthesis pitch error caused by the gear machining error, the installation error, and the deformation of the gear, the base pitch of two gears is no longer equal, so the meshing point of the gear pair deviates from the theoretical meshing line at the meshing-in and meshing-out position. Thus, the gear transmission ratio and meshing force change abruptly, which leads to vibration and noise.

The formation mechanisms of the above two meshing impacts are the same, but the effect of the meshing-in impact is more significant than the meshing-out impact. Thus, only the meshing-in impact is considered in this study. The schematic of the meshing-in impact is shown in Figure 5.

In Figure 5, the points , indicate the centre of the pinion and the wheel, the pitch point is , the point is the starting position of the out-of-line meshing, while point is the starting position of normal meshing, the inversion point of point is , and the combined deformation amount of the meshing teeth can be obtained by LTCA method. According to the geometric relationship of and , the following equations can be obtained:where , , , , , are the tooth tip circle radius, base circle radius, and pitch circle radius, and the centre distance is .

Substituting equations (14)–(19) into equation (20), can then be obtained, as well as the impact time :where is the angular velocity of the pinion.

According to the results of the above theoretical derivation, the meshing impact time is only related to the meshing teeth comprehensive deformation when the gear size parameters are determined. The load torque is an important factor affecting the . Figure 6 shows the change of the impact time under different loads.

It is seen from Figure 6 that the percentage of the meshing impact in the meshing period is within 3%–8.4%, and as the load torque increases, the percentage gradually levels off, because the increase of load torque will enlarge the contact area to increase the load-carrying capacity of the tooth surface.

According to the theory of impact dynamics, the amplitude of the impact force can be obtained by the following equation:

Here, indicates the impact velocity of the point of meshing-in, and is the single-tooth compliance of the point meshing-in, and the variables , , and indicate the rotational inertia of driving and driven wheels and the gear transmission ratio. The driving wheel base radius is , while the instantaneous base circle radius of the driven wheel is .

Since the impact time is relatively short, the meshing impact force is usually reduced to a sawtooth wave function. Based on the percentage of the meshing impact and the amplitude of the impact force obtained above, the meshing impact force curve is shown in Figure 7.

5. Dynamic Analysis at High Speed

The vibration accelerations in all directions of the gear are synthesized as the relative vibration acceleration in the line of action on the end face of the gear, which is regarded as the main index to measure the vibration of the system, as is shown in

Then, the root mean square (RMS) value can be calculated, and the vibration acceleration diagram can be obtained. The dynamic characteristics of the helical gear transmission system are then analyzed under three kinds of excitation conditions, excitation of TVMS, excitation of meshing impact, and excitation of both of them.

5.1. Excitation of TVMS

The TVMS curve obtained above is substituted into the dynamic equations, and the vibration responses of the single-stage gear system are obtained under the input torques of , , and . The vibration acceleration diagrams are shown in Figure 8.

As is shown in Figure 8, the 2nd- and 3rd-order superharmonic resonance occur when the rotating speed approaches 1/2 and 1/3 resonance speed N0. The vibration increases with the increase of load torque. The RMS value of the vibration acceleration is not significantly increased with the increase of the rotating speed excluding the resonance peak; on the contrary, there is a slight downward trend after the resonance region. That is, the effect of speed increasing on the vibration characteristics is not significant under the excitation of TVMS when the speed is far away from the resonance speed. This is because the TVMS is a parameter excitation. The vibration of the system under parametric excitation is called parametric vibration. The system has no external force input except the driving force, so the system vibration is caused by the periodic change of the system internal parameters. As the rotating speed increases to the off-resonance zone, the frequency of parameter excitation gradually moves away from the system resonance frequency, so the influence of this parameter change on the system also decreases. Therefore, the system vibration does not increase significantly with the rotating speed increase under the excitation of TVMS in the off-resonance zone.

To analyze the dynamic characteristics of the two-stage helical gear transmission system, and compare the dynamic characteristic differences between the single-stage and two-stage helical gear transmission system, the vibration responses of the above two systems are analyzed under the excitation of TVMS at 3,000 rpm and 15,000 rpm, as is shown in Figures 914.

Figures 9 and 10 show the vibration responses of a single-stage helical gear transmission system at 3,000 rpm and 15,000 rpm. As can be seen, the phase diagram of the single-stage helical gear transmission system is a closed curve under the excitation of TVMS. And the Poincaré section is a single point. It is an obvious periodic appearance for the vibration response. And the system response frequency only has its meshing and multiplication frequency; the frequency component is relatively simple. Figures 1114 show the vibration responses of the input-stage gear pair of the two-stage helical gear transmission system at 3,000 rpm and 15,000 rpm. represents the meshing frequency of the input gear pair, and represents the meshing frequency of the output gear pair. As can be seen from (b) and (d) in Figures 1114, the phase diagram of the two-stage helical gear transmission system under the excitation of TVMS is a series of curves intersecting with each other, and the Poincaré section is a series of discrete points that can form a closed curve. Therefore, the vibration responses show distinct quasiperiodic characteristics. The nonperiodic vibration response of the two-stage helical gear transmission system is significantly enhanced, compared with the single-stage helical gear transmission system.

As can be seen from (c) in Figures 1114, the frequency domain response component of the two-stage gear system becomes complicated. In addition to its own meshing frequency component, the vibration frequency of the input-stage gear pair also includes the meshing frequency of the output-stage gear pair and the frequencies of the overlapping and combination of the two meshing frequencies. This phenomenon also exists in the vibration frequency domain diagram of the output-stage gear pair. Since the input-stage gear pair is the power source of the output-stage gear pair, this phenomenon is more obvious in the output-stage gear pair, as shown in Figure 12(c); in the frequency response of the output-stage gear pair, the higher amplitudes are mostly related to the meshing frequency of the input-stage gear pair; however, the response frequency amplitude related to the meshing frequency of the output-stage gear pair is not that large.

This phenomenon indicates that the vibrations of the two gear pairs in the two-stage helical gear transmission system affect each other, which makes the vibration of the system complicated. This is because the two-stage helical gear transmission system can be regarded as a multifrequency excitation system. In addition to the excitation frequency generated by its own, the vibration frequencies generated by the input-stage gear pair also have an excitation effect on the output-stage gear pair. At this time, the response frequency of system vibration is overlapping and a combination of the two gear pairs meshing frequency and its multiple frequencies. Due to the fact that the input-stage gear pair is the power source of the output-stage gear pair, the influence of the input-stage gear pair on the output-stage gear pair is more obvious. This phenomenon is a prominent feature of the vibration response of two-stage helical gear transmission systems. Therefore, the response of system vibration becomes more complicated and the system’s aperiodicity is enhanced.

5.2. Excitation of Meshing Impact

When substituting the meshing impact obtained above into the dynamic equations, the vibration response of the single-stage gear system is obtained under the input torques of , , and . The vibration acceleration diagram is shown in Figure 15.

Figure 15 is a vibration acceleration diagram of three different loads under the excitation of meshing impact. It shows that the vibration of the gear increases with the increase of the load. The 2nd- and 3rd-order superharmonic resonance occur when the rotating speed approaches 1/2 and 1/3 its resonance speed N0. The RMS value of the vibration acceleration is significantly increased with the increase of the rotating speed excluding the resonance peak. The influence of rotating speed increasing on the system vibration is significantly far away from the resonance speed. This is because the meshing impact force increases with the increase of rotating speed, and the growth trend is approximately a linear relation. This is consistent with the conclusion in a previous literature that the dynamic load coefficient increases with the increase of line speed [25].

To analyze the time domain characteristics of the helical gear system under the excitation of meshing impact, the time domain diagrams of the single-stage helical gear system at two different rotation speeds are shown in Figures 16 and 17.

In Figure 16, the main manifestation of the meshing impact zone is the abrupt change in the negative direction of vibration acceleration. It is because the rotating speed of the driven wheel will jump sharply when the meshing impact occurs, and the vibration acceleration of line of action will decrease sharply, combining with equation (23). In Figure 16, the vibration acceleration of the system is gradually attenuated by Td, and 1/Td is the natural frequency of the system. The meshing period Tz decreases gradually with the increase of rotational speed. The system will resonate when the meshing period Tz equals the vibration acceleration attenuation period Td.

Figure 17 shows a time domain diagram of single-stage helical gear transmission system at 8,900 rpm under the excitation of meshing impact, where 8,900 rpm is the resonance speed of the single-stage helical gear transmission system.

Figure 17 clearly shows that the next impact is entered about one period of Td after the previous impact, which is consistent with the conclusion above.

For further comparison of the dynamic characteristics of two-stage and single-stage helical gear transmission systems, the time domain diagrams of the input-stage gear pair and the output-stage gear pair are analyzed under the excitation of meshing impact at 3,000 rpm and 15,000 rpm, as shown in Figures 18 and 19.

As can be seen in Figure 18(a), the two-stage helical gear transmission system under the excitation of meshing impact shows an abrupt change of vibration acceleration, which is consistent with the result in Figure 16. In addition to the abrupt change of vibration acceleration caused by its impact force, the output-stage gear pair has an abrupt change at the meshing impact positions corresponding to the input-stage gear pair. This is due to the influence of the input-stage gear pair on the output-stage gear pair. Because of the opposite rotational speeds of the two gear pairs, the directions of the acceleration abrupt change are opposite. In addition, the effect of the output-stage gear pair on the input-stage gear pair is not obvious as shown in Figure 18.

Figure 19 shows that, at the meshing time corresponding to the meshing impact of the input-stage gear pair, the output-stage gear pair has an obvious abrupt change of vibration acceleration, which is also consistent with the analysis in Figure 16. The difference is that the influence of the meshing impact of the output-stage gear pair on the input-stage gear pair appears in Figure 19. As can be seen in Figure 19, the influence of the input-stage gear pair on the output-stage gear pair is significantly greater; this is mainly because the input-stage gear pair is the power source of the output-stage gear pair, and the input-stage gear pair has a relatively high rotation speed. Furthermore, the higher the input rotation speed is, the more obvious the interaction between the two gear pairs will be.

5.3. Excitation of Both

Substituting the TVMS and meshing impact obtained above into the dynamic equations simultaneously, the vibration responses of the single-stage gear system are obtained under three input conditions of , , and . The vibration acceleration diagram is shown in Figure 20.

Figure 20 is a vibration acceleration diagram of the comprehensive excitation under three kinds of loads. As can be seen, the vibration of the gear increases with the increase of the load. And 2nd- and 3rd-order superharmonic resonance occur when the rotating speed approaches 1/2 or 1/3 its resonance speed N0; and the RMS value of the vibration acceleration is significantly increased with the increase of the rotational speed excluding the resonance peak. However, the tendency of the increase is slower than under the excitation of meshing impact. This is due to the slight decrease in system vibration with the increase of the rotation speed under the excitation of TVMS in the high-speed off-resonance region.

The vibration acceleration diagrams of the single-stage helical gear transmission system under three different excitations with the input torque of are compared in Figure 21. The system vibration under the excitation of meshing impact is larger than that under the excitation of TVMS. With the increase of the rotating speed, the system vibration under the excitation of meshing impact tends to be consistent with the excitation of both in the off-resonance zone. That is, the meshing impact is the main excitation component at high speed.

To compare the system responses under the three excitation modes, the diagrams of the frequency domain under three kinds of excitations for single-stage helical gear system at 12,000 rpm are shown in Figure 22, while under the excitation of TVMS independent, the third harmonic frequency is attenuated to zero. When the impact excitation component exists, the twentieth harmonic frequency still has a certain amplitude and its attenuation is slow. It is because the excitation frequency of the meshing impact is complicated; the system generates a lot of high-frequency harmonics. That is, the presence of the excitation of meshing impact makes the system vibration more complicated. This phenomenon is more prominent in the situation of high rotation speed. Therefore, more attention should be paid to the impact component in the optimization of vibration in the high-speed gear system.

6. Effect of Tooth Surface Modification

Tooth surface modification is one of the commonly used methods to reduce vibration and noise in gear transmission. After proper modification of the gear, the meshing performance of the gear system can be effectively improved, the vibration and noise of the system can be reduced, and the service life of the gear system can be improved.

To analyze the effects of the two modification schemes, modification of the input-stage pinion and of the output-stage pinion, on the vibration reduction of the PEV high-speed two-stage helical gear transmission, the dynamic model is established by Romax software and shown in Figure 23, and the corresponding tooth surface modifications of the pinions are shown in Figures 24 and 25.

Then, the LTE curves of the two gear pairs, with and without modification, are simulated by Romax software, as shown in Figures 26 and 27.

In Figure 26, the peak-to-peak value of the LTE is reduced from 1.47 to 0.75, which is reduced by 49% after the modification of the input-stage pinion, while in Figure 27, the peak-to-peak value of the LTE is reduced from 0.95 to 0.69, which is reduced by 27.4% after the modification of the output-stage pinion.

To analyze the vibration reduction, the TVMS curves of two gear pairs are calculated based on LTEs. The TVMS curves are then substituted into the dynamic equations (3), and the system response is compared and analyzed after solving the equation. Figures 28 and 29 show the time domain responses at input speed 12,000 rpm under two kinds of modification schemes, separately.

Figure 28 shows the vibration response under the modification of the input-stage pinion. As can be seen from Figure 28(a), the fluctuation of the vibration acceleration in the input-stage gear pair is reduced, and the abrupt change of the vibration acceleration caused by the impact force is also greatly reduced. In Figure 28(b), the fluctuation of the vibration acceleration of the output-stage gear pair is reduced, and the abrupt change of the vibration acceleration caused by the impact of the input-stage gear pair is also reduced to a certain extent. However, the abrupt change of the acceleration caused by its impact force is not changed obviously.

As can be seen from Figure 29, after the modification of the output-stage pinion, the abrupt change of the vibration acceleration caused by the meshing impact of the output-stage gear pair is reduced. However, the vibration acceleration of the input-stage gear pair is not significantly changed.

The phenomena in Figures 28 and 29 coincide with the results of the previous analysis that the vibration between the input-stage gear pair and the output-stage gear pair affects each other, and the influence of the former on the latter is obviously greater than that of the latter on the former.

To compare the effect of vibration reduction under the two kinds of modifications, the dynamic equations are solved, and the RMS values of the vibration acceleration are calculated at various input rotation speeds, as shown in Figures 30 and 31, where the vibration acceleration reduction rate represents the percentage reduction of RMS value of vibration acceleration on each stage before and after modification.

Figure 30 shows that, after modification of the input-stage pinion, not only does the input-stage gear pair have a damping effect (41.2%∼48.1%), but also the vibration of the output-stage gear pair is significantly reduced (10.3%∼37.8%). Figure 31 shows that, after modification of the output-stage pinion, the vibration of the output-stage gear pair is significantly reduced (4.2%∼23.3%), but the vibration of the reduction of the input-stage gear pair is not significant (0.3%∼1.2%).

According to the analysis of the system dynamic characteristics in the previous section, the phenomena in Figures 30 and 31 are due to the fact that the vibration of the input-stage gear pair has a great influence on the output-stage gear pair. Thus, the modification of the input-stage pinion will significantly weaken the vibration on both stages. Therefore, more attention should be paid to the input-stage gear pair in the modification design of multistage gear transmission.

7. Conclusions

The RMS value of the vibration acceleration is not significantly increased with the increase of the rotating speed excluding the resonance region under the excitation of TVMS. The influence of rotating speed increasing on vibration is not significantly away from the resonance speed. The RMS value of the vibration acceleration is significantly increased with the increase of the rotating speed under the excitation of meshing impact and the excitation of both. However, the increasing tendency of the vibration acceleration under the excitation of both is slower than that of under the excitation of meshing impact.

In the high-speed off-resonance region, the influence of meshing impact on the system vibration is greater than TVMS. That is, the meshing impact becomes the main excitation component at high speed.

The vibration of the two gear pairs has a mutual effect in the two-stage gear transmission. The main manifestation is that the vibration response frequency of one gear pair contains not only its own meshing frequency, but also the coupling frequencies. This phenomenon makes the vibration response of the high-speed gear transmission system more complicated.

Modification of the input-stage pinion will lead to significant vibration reduction for both stages in the two-stage high-speed helical gear transmission, about 41.2%∼48.1% reduction for the input-stage and 10.3%∼37.8% for the output-stage, while it only has a limited effect on its own stage when the output-stage pinion is modified, about 4.2%∼23.3% reduction for output-stage but only 0.3%∼1.2% for the input-stage. Therefore, more attention should be paid to the input-stage gear pair in the modification design of multistage gear transmissions.

Data Availability

The data used to support the conclusions of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the financial aid and support from the Key Laboratory of Shaanxi Province for Development and Application of New Transportation Energy (Chang’an University) Open Foundation (grant no. 300102220507), the National Natural Science Foundation of China (grant nos. 51705419 and 61701397), the China Postdoctoral Science Foundation (grant no. 2019M653702), the Postdoctoral Research Project of Shaanxi Province (grant no. 2018BSHEDZZ10), and the Shaanxi Provincial Department of Education Service Local Special (Industrialization) Plan Project (19JC030).