Research Article | Open Access
Wenjun Hu, Quanman Zhao, Yao Liu, Zhigang Li, Xianghui Kong, "Damage Evaluation of the Paving around Manholes under Vehicle Dynamic Load", Advances in Materials Science and Engineering, vol. 2020, Article ID 6510439, 11 pages, 2020. https://doi.org/10.1155/2020/6510439
Damage Evaluation of the Paving around Manholes under Vehicle Dynamic Load
It is very important for pavement engineers to know which factors are the main reasons for the damage of the paving around manholes. Based on the investigation on the damage of paving around manholes, a vibration model with multidegree of freedom for the vehicle-manhole cover was established and analyzed. After that, the Matlab software was used to obtain the variation law of impact load over time, and the 95% fourth power of the aggregate force was used as the index to evaluate the pavement damage. Finally, many influencing factors on pavement damage were analyzed by the method of grey correlation entropy. The results indicated that the impact load reached the maximum for the first time when the vehicle reached the top of the manhole cover, which was 1.29 times that of the static load, and the pavement damage coefficient was 2.12 times that of the static load. The influencing factors had different degrees of influence on the pavement damage; from large to small, they were change of road longitudinal slope > driving speed > damping of tire > stiffness of tire > height difference from the pavement damage > height difference from the manhole settlement > stiffness of the manhole cover.
Manholes play an important role in the installation, inspection, and maintenance of urban facilities. Many manholes are located on urban roads, and they are weak parts. Under repeated vehicle load, the paving around manholes easily sustain serious damage, such as settlement, cracks, and pits . Surveys were conducted on 300 paving around manholes in Guangzhou  and 212 in Tianjin . It was found that there were different types of damage, including settlement, cracks, spalling, and pits (Figure 1) .
When the paving around manhole sustains damage, it not only affects the beauty of the pavement but also worsens pavement roughness. When a vehicle passes over uneven pavement, a significant “bump” occurs. At the same time, as a rigid plate with limited thickness, the manhole cover undergoes obvious deformation and vibration under the vehicle load, which would lead to increased vibration of the vehicle, resulting in greater vehicle impact load, accelerating the destruction of the pavement, which would in turn increase the vibration of the vehicle, forming a vicious circle. In the maintenance of urban roads, the cost of repairing the paving around manholes is very high each year, but there remains the problem of further destruction very soon .
At present, research on manholes and the paving around manholes mostly focuses on manhole settlement mechanisms , new manhole structures [7–9], damage analysis of the paving around manholes , and treatment technology for pavement damage [11, 12], and the control of damage to paving around manholes has improved to some extent.
In the above studies, the dynamic load of a vehicle was a key factor. There were three main ways to obtain the dynamic load of a vehicle: software analysis, mathematical modeling, and the field test. Liu et al. , Chen et al. , and Liu et al.  used a multibody dynamics software SIMPACK to build a heavy-duty vehicle model and studied the dynamic load characteristics of vehicles under random excitation of pavement unevenness. The 1/4 vehicle model [16, 17], 1/2 vehicle model , and full car model [19, 20] were, respectively, established based on the mathematical method, and the dynamic load characteristics, driving comfort, and their influencing factors were studied. Liu et al.  used an accelerometer to monitor the acceleration of the front and rear axle and obtained the dynamic load.
When the vehicle passed over the paving around the manhole, the vibration characteristics of the vehicle were significantly different from those when passing over the general pavement. This was because the vehicle and the pavement were a weakly coupled system , and the deformation and vibration of the pavement were often neglected. But, when the vehicle passed over the manhole, there was obvious deformation and vibration of the manhole cover, which could not be ignored. However, a static load was often used in the study of the failure mechanism for the paving around manholes, while dynamic load was often ignored, meaning the failure mechanism of the paving around manholes is not yet comprehensively understood.
Therefore, based on the investigation on the damage of paving around manholes, this paper studies the dynamic characteristics of a vehicle when it passes over a manhole and lays the foundation for research on the failure mechanism of the paving around manholes. In the research, the deformation and vibration of the manhole cover is considered, a vibration model with multidegree of freedom for the vehicle-manhole cover is established, and the 95% fourth power of the aggregate force is used to evaluate the pavement damage. At the same time, the variation law of impact load over time is studied, and many influencing factors on pavement damage are analyzed with the method of grey correlation entropy, the key factors are identified, and provide reference for maintenance and design of paving around manholes.
From March 10 to 25, 2019, investigations including on damage type to the paving around manholes and settlement of manholes were carried out on four roads—Jingshidong (300 manholes), Xinluo (100 manholes), Tianchen (100 manholes), and Xueshan (100 manholes)—in Jinan, Shandong Province, China. During the investigations, the damage was divided into four categories: manhole settlement, pavement cracks, pavement damage, and pavement deformation. Pavement damage included pit slot, looseness, and repaired areas. The deformation included settlement, rutting, and waves.
During the investigations, it was found that the settlement in the “Up” point was often less than that in the “Down” point (Figure 2), resulting in a change of road longitudinal slope, and had great influence on the impact load of passing vehicles. The investigation data of pavement damage and manhole settlement are shown in Table 1. The settlement of manholes and the changes of road longitudinal slope (Δi) are shown in Table 2. In Table 1, the region of the damage is equivalent to a circular region, represented by radius r.
As shown in Table 1, when the road had been used for a short time (Tianchen and Xinluo road had been used for 3 years; Xueshan road had been used for 5 years), the radius of the damaged region was small, mainly concentrated within the range 0.4–0.8 m. As the time increased (Jingshidong road had been used for 15 years), the radius of the damaged region gradually increased (range 0.7–1.2 m), and the proportion of the large radius of the damaged region gradually increased. Since the Jingshidong road had been used for 15 years and had reached the maximum design service life, it could be considered that the radius of the damaged region had reached the maximum of 1.2 m. Therefore, the paving around the manhole was defined as a circular area, with the circle center being the manhole center with a radius of 1.2 m.
It can be seen from Table 1 that 98% of the paving around manholes on Xueshan road had crack damage, and 56% had pavement damage; 100% of that in Jingshidong road had pavement damage (all had been repaired); 91% of that on Tianchen road had manhole settlement, and 51% had pavement damage; and 92% of that on Xinluo road had manhole settlement. Combined with the time that the road had been used, it was found that the probability of damage occurrence for manholes and the paving around manholes was 100%. The occurrence order of damage was manhole settlement ⟶ pavement cracks ⟶ pavement deformation ⟶ pavement damage (in need of repair).
As shown in Table 2, the settlement of manholes was generally small and concentrated in the range of 0–10 mm, mainly because the larger settlement had been treated, but there was still a small amount of settlement (>20 mm) that had not been treated. The value of was mainly distributed between −1% and 1%: the maximum value was 1.8% and the minimum value was −2.1%.
3. Vehicle-Manhole Cover Vibration Model with Multidegree of Freedom
3.1. Model Building
As a complex system of multimass vibration, a vehicle is generally considered to have 18 degrees of freedom (DOF). It is very complicated to build a model according to the actual situation; so, it needed to be simplified. The current simplified models are the 1/2 vehicle vibration model with 5 DOF  and the 1/4 vehicle vibration model with 3 or 2 DOF . In this paper, when the vehicle passed over the paving around manhole, a 1/4 vehicle vibration model with 3 DOF was built; when the vehicle passed over the manhole cover, the vibration and deformation of the manhole cover were considered to build the 1/4 vehicle vibration model with 4 DOF (Figure 3).
As shown in Figure 3, the process of a vehicle passing over manhole covers and the paving around manholes was divided into four stages: (1) on manhole paving ⟶ (2) on the manhole cover ⟶ (3) on manhole paving ⟶ (4) on normal pavement. In stage 2, considering the deformation and vibration of the manhole cover, a 4 DOF vehicle vibration model was established, and a 3 DOF vehicle vibration model was established in the other stages. In Figure 3, m1, m2, m3, and m4 are the weights of the vehicle seat, frame, wheel, and manhole cover; y1, y2, y3, and y4 are the displacements; k1, k2, k3, and k4 are the stiffness coefficients; c1, c2, c3, and c4 are the damp coefficients; is the pavement unevenness; and is the vehicle speed.
It was assumed that the pavement roughness was good except for the manhole settlement and any damage to the paving around the manhole. Considering the damage as an unevenness incentive in the calculation process, the 3-DOF vibration model was established:
The 4 DOF vibration model was established:
When the vehicle was driving on the manhole cover, the impact load factor (f) was calculated by equation (3). When the vehicle was driving on other areas, the calculation of f still used equation (3) but :where is the impact load and calculated by equation (3), and F is the static load:
3.2. Calculation Example
We assumed that the maximum height difference caused by the damage to paving around the manhole was H1 = 1 cm, and the position was 0.6 m from the edge of the manhole cover; the diameter of the manhole cover was 0.7 m, and its settlement was H2 = 1 cm; if and = 10 m/s, the driving direction was uphill. In stage (1), the initial conditions of equation (1) were , y1 = 0.01, y2 = 0.01, and y3 = 0.01. In stage (2), the initial conditions of equation (2) were , , , , y1 = y1 (t1) + 0.01, y2 = y2 (t1) + 0.01, y3 = y3 (t1) + 0.01, y4 = 0 (t1 was the duration of the vehicle driving on the paving around the manhole), and so on for the remaining stages. The transfer matrix method  was used to solve equations (1) and (2). Combined with equation (3), equation (4), and the initial conditions and the programmed Matlab software, the change of f with time was obtained (Figure 4).
It can be seen from Figure 4 that under the excitation of pavement roughness, the vehicle load was a simple harmonic vibration of amplitude attenuation. When there was a new unevenness excitation, the vehicle dynamic load changed obviously. The f was 1.12 at 0.06 s when the vehicle was on the “Up” point, 1.29 at 0.13 s when the vehicle was on the “Down” point, and reached a maximum value of 1.33 at 0.64 s. This was the reason why the settlement of the “Down” point was bigger than the “Up” point. Considering the subgrade around the manhole is difficult to compact, and the paving around the manhole was most likely to be damaged.
4. Evaluation of Pavement Damage and Influencing Factors
4.1. Evaluation of Pavement Damage
At present, there are three main indexes for the evaluation of pavement damage: dynamic load coefficient (D), road stress factor (θ), and 95% fourth power of the aggregate force (J) . In D, the relation between pavement damage loads is linear, but the repeatability of the vehicle dynamic load in the spatial distribution is not considered, which was quite different from the actual situation. Later, through many experiments, it was proved that the fatigue damage of a semirigid base asphalt pavement was directly proportional to the fourth power of vehicle load. Therefore, Potter et al. proposed θ  but still did not consider the spatial distribution of the vehicle dynamic load. Based on this, Cole and Cebon  put forward J (equation (5)). In J, the repeatability of the vehicle dynamic load in spatial distribution is considered and greatly improves the accuracy of pavement damage evaluation. Therefore, J was used as the pavement damage evaluation index in this paper. Under a static load, the index was equal to 1, and the larger the value, the greater the damage to the pavement:where A4 is the fourth power of the aggregate force, is the standard deviation of A4, and is the average of A4.
Since the vibration caused by the vehicle passing over the manhole was transient vibration, referring to the literature , we calculated J within 1 s and started timing when the vehicle entered the paving around the manhole area. With the basic parameters described in Table 3 and Section 3.2, J = 2.12 > 1 calculated by equations (3) and (5). When the vehicle passed over the paving around the manhole and the manhole cover, the damage of the vehicle load to the pavement was significantly greater than that of a static load.
4.2. Influencing Factors
There are many factors affecting the dynamic load and pavement damage . Seven factors, the vehicle speed (), the height difference caused by pavement damage (H1), the manhole settlement (H2), the changes of road longitudinal slope (), the stiffness coefficient of the manhole cover (k4), the tire stiffness coefficient (k3), and the tire damping coefficient (c3), were analyzed. For each factor, was determined by the design speed of the urban road, H1 was determined by the pavement damage, H2 was determined by the manhole settlement, was determined by the uneven settlement of the manhole, k4 was determined by the material type of the manhole cover, and k3 and c3 were determined by the tire model and material. By equations (1)–(4), the time history curves for f under different influencing factors were calculated (Figure 5). By equation (5), J under different influencing factors was calculated (Table 4):(1)As can be seen from Figure 5(a) and Table 4, with the increase of vehicle speed, the maximum f and J increased gradually. When the vehicle speed increased from 10 km/h to 80 km/h, J increased by 52.6%, which was obvious.(2)As can be seen from Figures 5(b) and 5(c) and Table 4, with the increase of H1 and H2, the maximum f and J were significantly increased. When H1 increased from 0.5 cm to 8 cm, J increased by 48.5%; When H2 increased from 0.5 cm to 8 cm, J increased by 36.5%. Although the two factors of H1 and H2 were the same as the unevenness of the pavement, their influences on J were different. Especially, when H was larger than 1 cm, the influence of H1 was greater than the influence of H2. At this time, the presence of the manhole cover with elastic characteristics buffered the vehicle vibration and reduced J under the vehicle load.(3)As can be seen from Figure 5(d) and Table 4, with k4 increased, the maximum f decreased gradually and J decreased slightly; when k4 increased from 102 N/m to 109 N/m, J only decreased by 15.3%. Therefore, k4 had little influence on J.(4)As can be seen from Figure 5(e) and Table 4, with the increase of , the maximum f and J increased significantly. When increased from 1% to 8%, J increased by 33.1%. In fact, also belonged to the “unevenness” of pavement. Under the same excitation of pavement roughness, the larger , the bigger J, but this factor was often ignored.(5)As can be seen from Figures 5(f) and 5(g) and Table 4, with the increase of c3, the maximum f and J decreased gradually. When c3 increased from 103 to 104, J decreased by 13.2%. As k3 increased, the maximum f increased significantly, the vehicle vibration period shortened, and J increased first and then decreased. When k3 was about 8 × 105, J reached the maximum value of 2.31. It can be seen that k3 and c3 had a certain degree of influence on the pavement damage, but in fact, k3 and c3 were mainly affected by the tire pressure, and as the tire pressure was basically stable, k3 and c3 varied little; therefore, J changed little.
Through the above analysis, it is clear that all the above factors had a certain degree of influence on J. However, among these factors, which factors were more important, played a major role, and which factors played a secondary role? As the units of the above factors were different, comparing the change of J with the change of each factor was not significant. It was necessary to select a suitable evaluation method to compare the influence degree of each factor.
5. Grey Correlation Entropy Method to Compare the Degree of Influence of Each Factor on J
5.1. Grey Correlation Entropy Method
Grey correlation refers to the uncertain relationship between different things. It is a systematic analysis method to measure the degree of correlation between factors and systems and to compare the influence degree between each factor. However, this method can easily generate the problem that the local point correlation value controls the overall point correlation value, causing loss. To solve this problem, the grey correlation entropy method was proposed .
For grey correlation entropy analysis, let X be the grey correlation factor set. Take J as the reference set (), and take the 7 influencing factors as the comparison set . In the calculation process, the magnitude and unit of each factor were different. Therefore, in order to eliminate the influence on the calculation result, the data needed to be processed to obtain the dimensionless set Xi and the reference set X0, which was .
The process of grey correlation entropy analysis is summarized into the following 5 steps: Step 1. The grey correlation coefficient of Xi–X0 was obtained from the following equation: where is the grey correlation coefficient, and is the resolution coefficient ( ∈ [0, 1] and is 0.5 generally). Step 2. The distribution density of grey correlation entropy was obtained according to equation (7). Let , and the defined distribution density of grey correlation entropy is as follows: where , and , . Step 3. The grey correlation entropy was obtained from equations (8) and (9). Let X = (x1, x2, …, xn), , xi ≥ 0, and , and grey correlation entropy is as follows: where xi is attribute information. The grey correlation entropy of Xi is Step 4. The grey entropy correlation degree of Xi was obtained from the following equation: where is the maximum value of the grey entropy and represents the maximum value of the difference information column composed of n elements. Step 5. According to the magnitude of the grey entropy correlation degree, which factor was more important was determined. The larger the value, the more important the factor.
As can be seen from Table 5, the order of degree of influence of each factor on J was > > c3 > k3 > H1 > H2 > k4. Among them, c3 and k3 were mainly determined by the tire pressure, but tire pressure was basically stable; so, the range of c3 and k3 was small and had little influence on J . Therefore, in the analysis of J, , , H1, and H2 were the 4 major factors that needed consideration.
Considering the deformation and vibration of the manhole cover, the coupled vibration model of the vehicle-manhole cover was established, and the road damage was evaluated by J. The research results provide the theoretical basis for the maintenance of the paving around manholes, and the following conclusions can be drawn:(1)When there was a new unevenness excitation, the vehicle dynamic load changed obviously. With the basic parameters, f reached the maximum value of 1.33 at 0.64 s, and J was 2.12. The unevenness of the pavement made the paving around the manhole easier to damage.(2)When there was damage in the paving around the manhole and the manhole cover, which made the pavement roughness worse, the impact load of the vehicle was greatly increased, and the damage to the pavement was accelerated, and the repair of this pavement should be timely. Although the two factors of H1 and H2 were the same as the unevenness of the pavement, their influences on J were different. Especially, when H was larger than 1 cm, the influence of H1 was greater than the influence of H2. At this time, the presence of the manhole cover with elastic characteristics buffered the vehicle vibration and reduced J under the vehicle load.(3)According to the analysis of grey correlation entropy, the order of degree of influence of each factor on J was > > c3 > k3 > H1 > H2 > k4. Among them, c3 and k3 were mainly determined by the tire pressure, but tire pressure was basically stable, so the range of c3 and k3 was small and had little influence on J. Therefore, , , H1, and H2 were the 4 major factors that had significant influence on the J.
However, during the analysis of pavement damage, the influence of load time was not considered, which affected the accuracy of the analysis results. It is suggested that in subsequent research, the influence of load magnitude and time on pavement damage should be considered, and the pavement damage problem should be studied in depth.
The test data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This research was supported by the project ZR2018BEE039 supported by Shandong Provincial Natural Science Foundation and project 2019GSF109067 supported by the Key Research and Development Program of Shandong Province. The authors gratefully acknowledge the financial support.
- J. Zhou, “Settlement law and improvement measures of the manhole on the roadway,” Chongqing University, Chongqing, China, 2013, Master’s thesis.
- Z. H. Jin, “Analysis of numerical calculation of urban road inspection manhole installation and pavement smoothness,” South China University of Technology, Guangzhou, China, 2018, Master’s thesis.
- H. Li, “Prevention and treatment research on pavement distresses around inspection wells,” Hebei University of Technology, Tianjin, China, 2014, Master’s thesis.
- Q. Li and J. Q. Liu, “Asphalt pavement evenness deterioration analysis based on the vehicle-pavement interaction,” Jounal of Vibration and Shock, vol. 37, no. 6, pp. 76–81+116, 2018.
- M. Q. He and L. Y. Wei, “Mechanical analysis of urban road inspection shaft,” Journal of Hebei University of Technology, vol. 3, pp. 109–112, 2014.
- T. Tobita, G. C. Kang, and S. Iai, “Estimation of liquefaction-induced manhole uplift displacements and trench-backfill settlements,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 138, no. 4, pp. 491–499, 2012.
- A. Abbas, F. Ruddock, R. Alkhaddar, G. Rothwell, and R. Andoh, “Improving the geometry of manholes designed for separate sewer systems,” Canadian Journal of Civil Engineering, vol. 46, no. 1, pp. 13–25, 2019.
- Z. Luo, K. Dong, M. Guo et al., “Composite manhole covers prepared with recycled printed circuit boards as a reinforcing filler,” Polymer Composites, vol. 39, no. 12, pp. 4532–4541, 2018.
- J. Stoikes, “Manhole risers cut costs and preserve pavements,” Pavement Maintenance & Reconstruction, vol. 29, no. 3, pp. 41–50, 2014.
- L. Y. Chen, M. Y. Wang, C. Chao, and W. Lo, “Assessment of asphalt concrete pavement quality by using infrared thermal imaging technology,” Journal of Marine Science and Technology, vol. 23, no. 3, pp. 331–338, 2015.
- J. R. Chang, S. R. Yang, and K. H. Lin, “Performance of new resin cement in rapid manhole vertical adjustment,” Materials Research Innovations, vol. 18, no. 2, pp. 857–862, 2014.
- Y. Wang, L. Kong, Q. Chen, B. Lau, and Y. Wang, “Research and application of a black rapid repair concrete for municipal pavement rehabilitation around manholes,” Construction and Building Materials, vol. 150, pp. 204–213, 2017.
- D. W. Liu, Z. H. Dai, Y. Chen, and H. M. Chen, “Dynamic stress response of flexible pavement under multiple wheel dynamic loads of vehicle,” China Journal of Highway and Transport, vol. 30, no. 11, pp. 36–44, 2017.
- Y. Chen, Z. H. Dai, H. M. Chen, and D. W. Liu, “Strain analysis of a flexible asphalt pavement under multi-wheel random dynamic loads of vehicles,” Jounal of Vibation and Shock, vol. 35, no. 19, pp. 15–19+44, 2016.
- X. Y. Liu, Y. Chen, P. Gao, and D. W. Liu, “Simulation analysis for multiaxial heavy vehicles’ dynamic load characteristics under bilateral tracks’ road excitation,” Jounal of Vibation and Shock, vol. 34, no. 13, pp. 48–52, 2015.
- X. Shao, F. Naghdy, H. Du, and Y. Qin, “Coupling effect between road excitation and an in-wheel switched reluctance motor on vehicle ride comfort and active suspension control,” Journal of Sound and Vibration, vol. 443, no. 43, pp. 683–702, 2019.
- H. Zeng, H. Park, and B. L. Smith, “Impact of vehicle dynamic systems on a connected vehicle-enabled pavement roughness estimation,” Journal of Infrastructure Systems, vol. 25, no. 1, pp. 426–434, 2018.
- S. Li, S. Yang, and L. Chen, “A nonlinear vehicle-road coupled model for dynamics research,” Journal of Computational and Nonlinear Dynamics, vol. 8, no. 2, pp. 101–114, 2013.
- S. H. Li and S. P. Yang, “Investigation on dynamics of a three-directional coupled vehicle-road system,” Journal of Vibroengineering, vol. 17, no. 7, pp. 125–136, 2015.
- Y. Lu, S. Yang, and S. Li, “Research on dynamics of a class of heavy vehicle-tire-road coupling system,” Science China Technological Sciences, vol. 54, no. 8, pp. 2054–2063, 2011.
- D. W. Liu, J. Chen, W. Huo, and B. C. Chen, “Experimental study on random dynamic load of vehicles on pavement,” Transactions of the Chinese Society for Agricultural Machinery, vol. 7, pp. 12–14+19, 2005.
- Z. L. Zheng and P. M. Lv, “Dynamic load of large trucks acting on pavement,” Journal of Chang’an University (Natural Science Edition), vol. 29, no. 2, pp. 101–106, 2009.
- H. L. Zhang and W. Q. Yang, “Evaluation method of pavement roughness based on 5-DOF human-vehicle-road vibration model,” Jounal of Traffic and Tranportation Engoneering, vol. 10, no. 4, pp. 16–22, 2010.
- Q. M. Zhao, H. L. Zhang, and W. Jia, “Study on the contact force between vehicle and road and the influence factors in punch-out area,” Jounal of Hefei Univesity of Technology (Nature Science), vol. 38, no. 3, pp. 300–304, 2015.
- Q. M. Zhao, “Study on mechanism and prediction of the distress of punchout in continuously reinforced concrete pavement based on two scale model,” Chang’an University, Xi’an, China, 2015, Ph.D. thesis.
- D. W. Liu, H. M. Chen, W. Liu, and R. C. Jiang, “Evaluation of road damage based on a full heavy vehicle multi-body dynamic model,” Jounal of Vibation and Shock, vol. 29, no. 4, pp. 92–96+233, 2010.
- T. E. C. Potter, D. Cebon, and D. J. Cole, “Review paper: assessing “road-friendliness”: a review,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 211, no. 6, pp. 455–475, 1997.
- D. J. Cole and D. Cebon, “Truck suspension design to minimize road damage,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 210, no. 2, pp. 95–107, 1996.
- P. H. Zhao, “Study on the effects of formation and stability for emulsified asphalt,” China University of Petroleum (East China), Qingdao, China, 2013, Ph.D. thesis.
- X. T. Nie, L. X. Shi, H. Gu, and S. H. Zhu, “Research on the radial stiffness and damping of tractor coefficient tires through test,” Journal of Nanjing Agricultural University, vol. 34, no. 5, pp. 139–143, 2011.
Copyright © 2020 Wenjun Hu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.