Research Article  Open Access
Zheheng Chen, Shanwen Zhang, "EMBased Monitoring and Probabilistic Analysis of Prestress Loss of Bonded Tendons in PSC Beams", Advances in Civil Engineering, vol. 2018, Article ID 4064362, 9 pages, 2018. https://doi.org/10.1155/2018/4064362
EMBased Monitoring and Probabilistic Analysis of Prestress Loss of Bonded Tendons in PSC Beams
Abstract
The prestress level is a key factor of prestressed concrete (PSC) beams, affecting their longterm serviceability and safety. Existing monitoring methods, however, are not effective in obtaining the force or stress of embedded tendons. This paper investigates the feasibility of elastomagnetic (EM) sensors, which have been used for external tendons, in monitoring the longterm prestress loss of bonded tendons. The influence of ambient temperature, water, eccentricity ratio, plastic duct, and cement grouts on the test results of EM sensors is experimentally examined. Based on the calibrated EM sensors, prestress loss of a group of PSC beams was monitored for one year. In order to further consider the high randomness in material, environment, and construction, probabilistic analysis of prestress loss is conducted. Finally, the variation range of prestress loss with a certain confidence level is obtained and is compared with the monitored data, which provides a basis for the determination of prestress level in the design of PSC beams.
1. Introduction
Prestressed concrete (PSC) beams have been widely used in highway and railway bridges, in which the effective prestress is one key factor influencing the bridge serviceability and safety. However, it is well known that the timedependent prestress loss is a very complex process, influenced by inherent material properties, external loads or environmental effects, and construction details. [1]. Moreover, there are abundant uncertainties in such a process [2], resulting in difficulties in accurate prediction on prestress loss. As a result, excessive deflections and cracks are often observed in practice [3–6], which highlights the importance of improved design of prestress level of PSC beams.
In order to learn the authentic status of prestress loss in existing PSC beams, various longterm monitoring measures have been investigated. The vibrating wire strain gauges (VWSGs) can obtain the prestress loss due to shrinkage and creep of concrete [7–10], while they would not monitor the loss of stress relaxation. Similar problem exists when the optical fiber sensing is used [11, 12]. Due to the friction or bond of concrete, the stress along the tendon usually changes, and load cell at the beam ends can only obtain the total loss of prestress [13]. As to traditional resistance strain gauges, they are usually adhered to the surface of PT steels in the plastic corrugated pipes, while the narrow working space and the moisture may seriously affect the survival of gauges and their working performance. The elastomagnetic (EM) sensing is a nondestructive and noncontact monitoring technique that is based on the magneticelastic effects of ferromagnetic material [14]. At present, it has been usually used in stress monitoring of external tendons (i.e., unbonded tendons), such as the stayed cables, suspenders, and external tendons of PSC bridges, as shown in Figure 1. However, the application of EM sensing in bonded tendons requires further investigation due to the possible influence of plastic duct and cement grouts on measurement results.
(a)
(b)
(c)
In this paper, the feasibility of elastomagnetic (EM) sensors in monitoring the longterm prestress loss in bonded tendons is experimentally investigated, and the influence of ambient temperature, water, eccentricity ratio, plastic duct, and cement grouts on the test results is discussed. Based on the calibrated EM sensors, prestress loss of two PSC beams was monitored for one year. In order to further consider the high randomness in material, environment, and construction, probabilistic analysis of prestress loss is conducted based on the prediction model that can take into account the influence of nonprestressed steels and interaction among concrete creep, shrinkage, and steel relaxation. Finally, the variation range of prestress loss with a certain confidence level is obtained and is compared with the monitored data, which provides a basis for rational determination of prestress level in the design of PSC beams.
2. ElastoMagnetic (EM) Sensing
2.1. Principle of EM Sensing
The EM sensor is a nondestructive and noncontact test method, based on the inherent magneticelastic effect in the ferromagnetic material. The change in tension force can be obtained by measuring variations in magnetic permeability of prestressed tendons. As shown in Figure 2, the EM sensor consists of two layers of coils (i.e., the excitation and measurement coils) and the temperature sensor, and Figure 3 shows the measuring principle of EM sensor. The steel tendon will be magnetized when the pulse current flows into excitation coil, and pulsed magnetic field can be generated along the steel tendon. Due to electromagnetic induction, voltage is induced in the measurement coil. Magnitude of the induced voltage depends on the magnetic permeability of steel tendon, which is a function of the stress of the tendon, as shown in Equation (1).
The increment of magnetic permeability of the test specimen with respect to its zero stress state is expressed as follows:where represents the integral value of voltage when there is no test specimen in the EM sensor, denotes the integral value of voltage when there is a test specimen, is the area of section surrounded by measurement coil, and represents the net area of test specimen. The relationship between tension force and the increment of magnetic permeability can be expressed by the following equation:where ,, , and are the coefficients of the cubic polynomials, which can be obtained from the calibration tests of EM sensors. Note that the effect of temperature on the magnetic permeability of ferromagnetic materials should be eliminated during the test, and the variation of magnetic permeability due to temperature effect can be modified by the following formula:where and are the test and the reference temperatures, respectively; and are the increments of magnetic permeability corresponding to the two temperatures; and is the slope of the magnetic permeability vs. temperature relationship.
2.2. Tests on Influencing Factors of EM Sensor
A group of tests were designed to investigate the factors that may influence the performance of EM sensors. First, the influence of ambient temperature on the tests is experimentally investigated. It is known that the permeability of steel specimens changes with the temperature, while the slope of the relationship between the permeability and stress of steel specimens does not change within −30°C to 60°C [15]; therefore, the temperature calibration is carried out under the zero stress condition. In this study, a steel tendon with zero stress and the diameter of 15.24 mm was inserted in one EM sensor (type CCT43B) with the inner diameter of 43 mm, and they were put in a constant temperature box. The temperature started from 20°C and gradually rose to 50°C with an increment of 5°C. Under each temperature, the specimen was kept inside a temperature box for 1 hour before reading the relative magnetic permeability. The function curve of temperature and permeability is obtained, as shown in Figure 4, where linear relationship between temperature and the relative magnetic permeability is observed with the slope of 0.013.
Secondly, the influence of ambient media on the test results was analyzed, where there were four tests. In the test No. 1, the EM sensor was put in air while in the test No. 2, the EM sensor was immerged in water. No specimens were put inside the sensor in the test Nos. 1 and 2. In the test Nos. 3 and 4, the EM sensors were put in air and immerged in water, respectively, with a steel tendon inside the sensor. In each test, the integral voltages were measured for six times, and the ambient temperature was 16°C. According to the test results in Table 1, the integral voltages in water are almost the same as those in air, no matter there are steel specimens inside or not, showing that the water has no influence on the measurement results.

In addition, the inner diameter of EM sensor is usually significantly larger than the specimen; therefore, in some cases the specimen is difficult to be completely located in the centroid of EM sensor. To investigate the influence of eccentricity, another three tests were conducted, as shown in Figure 5. In test No. 5, the steel tendon with the diameter of 15.24 mm was located in the centroid of the sensor; in test No. 6 the tendon was located at the lowest position in the sensor (having the largest vertical eccentricity) and in test No. 7 the tendon was inclined. The test results presented in Table 2 indicate that the value in test No. 6 is lower than those in tests No. 5 and 7, showing that the eccentricity has certain influence on the measurement results; therefore, the sensor should avoid excessive eccentricity when used.

In practice, the steel tendons are usually embedded in corrugated ducts with or without grouts. In order to study the influence of tendon ducts and grouts, the test No. 8 and test No. 9 were conducted, where in the test No. 8, there is a corrugated plastic duct outside the tendon and in the test No. 9, the duct was grouted. As compared with the results in the test No. 3, the two tests in Table 3 are almost the same. Therefore, the corrugated ducts and cement grouts, being nonmetal materials, do not have any influence on the measurement results of EM sensor.

Finally, the relative magnetic permeability of the EM sensor was measured under various levels of loads, as illustrated in Figure 6, where the relative magnetic permeability vs. tension force relationship is basically linear, and the cubic polynomial is used to describe the relation.
3. LongTerm Monitoring of Prestress Loss of Bonded Tendons in PSC Beams
Based on the feasibility of EM sensors in monitoring the longterm prestress loss in bonded tendons, as depicted in Section 2, a series of tests were conducted, and two of them are taken for analysis in this paper. The simply supported concrete beams were made using the C55 concrete, which has the 28day cubic compressive strength of 55 MPa. The beams have a clear span of 4.2 m and a rectangular cross section of 200 × 300 mm, as shown in Figure 7. Highstrength and lowrelaxation 1860 grade steel strand was adopted as prestressed reinforcement, and the diameter of the reserved duct was 35.4 mm. The HRB335 grade twisted steel was used as longitudinal main reinforcement, the HPB235 grade round steel with diameter of 8 mm was adopted as transverse stirrups for the rest of the beams, and the stirrup spacing was 150 mm in all beams. The prestressing forces were applied directly to the target values at the age of 30 days. The initial tension forces of PC1 and PC4 are 100.1 kN and 135.9 kN, respectively, which were applied at one end of the beams with load cell installed to control the jacking force. Plain concrete blocks were applied as the secondary dead load at the two 1/3 beam spans immediately after the jacking, each with the weight of 3.125 kN. The detailed design parameters of the beams are shown in [16].
(a)
(b)
Figure 8 illustrates the timevarying prestress loss of beams PC1 and PC4, monitored using EM sensors, where the development of prestress loss with time is similar in pattern to concrete creep and shrinkage. The longterm loss developed rapidly in the early monitoring period, and the growth rate of timevarying loss becomes gradually gentle after about 250 days. Due to the higher initial tension of the tendon in PC4, the longterm loss at 1/3 section of PC4 is 62.9 MPa after one year of monitoring, which is about 1.42 times of that of PC1. In addition, the monitored longterm loss exhibits some fluctuations during the whole monitoring period, which is due to the variation of environmental factors such as humidity.
4. Probabilistic Analysis of LongTerm Prestress Loss of Bonded Tendons
4.1. Probabilistic Model of LongTerm Prestress Loss
Based on the method of ageadjusted effective modulus [17], the random prediction model for longterm prestress loss which takes into account nonprestressed steels and interaction among concrete creep, shrinkage, and steel relaxation is expressed in Equation (4), of which the detailed derivation can also be found in [16]:where , , , , and ; and represent the instantaneous elastic stresses of concrete at the centroids of PT tendons and nonprestressed steels at first loading time, respectively; , , and are the elastic moduli of nonprestressed steels, PT tendons, and concrete, respectively; and are the creep coefficient and shrinkage strain and aging coefficient at time , respectively.
In order to consider the uncertainties in creep model and shrinkage model, two coefficients, and , are adapted in this study. For the CEBFIP90 model, the statistical characteristic values of and can be taken as , and , [18], respectively; is the nominal creep coefficient; is the coefficient to describe the development of creep with time; is the nominal shrinkage strain; is the coefficient to describe the development of shrinkage with time; is the reduced relaxation loss and equals to ; denotes the reduction factor, and a value of 0.75 [19] is adopted in this study, and the calculation method of intrinsic relaxation, by Youakim et al. [20], is adopted in this paper.
The coefficients , , , and are as shown in the following equation:where and represent the ratio of nonprestressed steels and PT tendons, respectively; , where or represents the nonprestressed steels () or PT tendons (), and and are the eccentricity of nonprestressed steels and the PT tendons with respect to the centroid of concrete, respectively, and is the radius of gyration of the concrete section; is the aging coefficient, which varies between 0.6 and 1.0 [21], and the value of 0.8 is adopted in this study for simplicity.
4.2. Statistical Properties of Random Parameters
The mechanism and influence factors of concrete shrinkage and creep are very complicated, and there are usually large deviations between the predicted and measured values of effects of concrete creep and shrinkage. There are many reasons leading to this phenomenon, and one of the main reasons is that the selected prediction model of shrinkage and creep itself is not perfect. At present, most of the shrinkage and creep models are the empirical formulas obtained through statistical regression based on experimental data, and the models themselves have the fitting errors; therefore, model uncertainty factors are used to describe the deviation. Some scholars have conducted extensive research on uncertainty of shrinkage and creep models [22–24]. For the CEBFIP90 model, the statistical characteristic values of the model uncertainty factor are shown in Table 4.

Meanwhile, it is necessary to consider the randomness of parameters in the prediction model during the analysis of longterm loss of prestress, where the concrete compressive strength, ambient relative humidity, and loading age are treated as random variables. According to [25], the concrete strength used in this study follows the normal distribution with a mean value of 63.9 MPa and a coefficient of variation (COV) of 0.089. The ambient relative humidity is another significant factor which influences the timevarying effects of concrete structure, and the relative humidity in the laboratory exhibits seasonal fluctuations. According to the meteorological data of Nanjing, the probability density function (PDF) of ambient relative humidity is obtained, which follows a normal distribution with the mean value of 75% and the COV of 0.161. Other random factors include the density of concrete, initial loading age, elastic modulus of prestressing steel, areas of nonprestressed, and prestressed steel, which are described in Table 4. The geometric sizes of specimens are treated as design constants due to its minor variation. In order to facilitate stochastic analysis, random variables are assumed to be independent of each other.
4.3. Probabilistic Analysis and Discussion
According to the probabilistic model in Section 4.1 and the random factors in Table 4, probabilistic analysis of longterm prestress loss in beams PC1 and PC4 is carried out using the Latin hypercube sampling technique (5000 samples), and the timevarying PDFs are obtained, as shown in Figure 9. It can be seen that the PDFs of longterm prestress loss approximately follow the normal distribution at different ages, and the probability density curves gradually become flat as time goes by, indicating that at the same confidence level, the width of confidence interval of longterm prestress loss becomes larger. As presented in Figure 10, the standard deviation () of longterm loss increases continuously in each calculation condition from 1 month to 30 years after jacking, which shows that the discreteness of prestress loss due to the variation of influencing parameters gradually becomes larger with time, and the increase is most significant in the first ten years.
(a)
(b)
As shown in Figure 11, the measured longterm loss exhibits some fluctuations during the whole monitoring period, due to variations of ambient temperature and humidity. Although the test results of PC1 are basically consistent with the mean values of stochastic analysis, for PC2, the differences between test results and mean analyzed results become larger, which is due to the uncertainties in analysis of timevarying effect. Note that the mean values of prestress loss are usually adopted in structural design, which might underestimate their actual values and lead to unsafe design. Considering the prestress loss approximately obeying normal distribution at each time, the 95% confidence interval can be supposed as [], where represents the mean value of probability analysis of timevarying loss at time based on Equation (4), and represents the of timevarying loss at time . As shown in Figure 10, during the monitored period, the measured longterm prestress loss for beams PC1 and PC4 are located within the 95% confidence interval of probabilistic analyses. Taking PC4 as an example, the measured value of longterm loss is 70.6 MPa after about 1 year of monitoring, and the upper bound of 95% confidence interval is 97.1 MPa, which is about 37.5% larger than the measured value. Meanwhile, the width of the confidence interval of longterm loss gradually increased with time, varying from [4.3, 34.1] at one month of monitoring to [13.2, 97.1] at one year of monitoring, which indicates that the discreteness of timevarying loss becomes more obvious due to the larger variability of concrete shrinkage and creep as time goes by. Therefore, it is suggested that the upper bound of 95% confidence interval, instead of the mean value, should be used as the maximum longterm loss which may occur in the actual PSC structures.
(a)
(b)
5. Conclusions
This paper provides a preliminary investigation on the feasibility of EM sensors in monitoring the longterm prestress loss of bonded tendons. According to the presented study, conclusions are drawn as follows:(1)The performance tests of EM sensor show that the external media such as water, cement grouting, and plastic duct have no influence on the measurement results, showing that the EM sensor can be well used in steel members bonded in concrete. The temperature has certain influence on the magnetic permeability of the tested specimens, and the correction coefficient can be obtained through temperature calibration at the zero stress state. The eccentricity of the steel specimens relative to the EM sensor has some impact on the test results, and therefore, the steel specimens should be located in the centroid of the EM sensor when it is installed.(2)The longterm prestress loss due to shrinkage, creep of concrete, and stress relaxation is complicated, and there are significant uncertainties. Therefore, the uncertainty analysis of longterm loss should be carried out in order to make the design of bridges more safe and realistic. Based on the probabilistic analysis model adapted in this study, the probability analyses were conducted, which indicate that the longterm loss of prestress does not reject the normal distribution at each time, and its discreteness gradually become larger with the extension of time. The bandwidth surrounded by the upper and lower bound of 95% confidence interval also gradually increase with time. In addition, the predicted prestress loss in bridge design from deterministic analysis may be significantly different from the authentic monitored values. Therefore, to conduct a safe design, it is suggested that some suitable confidence limits such as 95% should be adopted instead of the mean value.(3)For the internal prestressed tendons of widespread PSC box girder bridges, the EM sensors installed in the construction stage can be used for the realtime monitoring of effective prestress force in the service process of bridges, which would provide the basis for the longterm performance assessment of the structures. These will be done in the next study.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The support from the Open Foundation of National Engineering Laboratory for High Speed Railway Construction under the Grant No. HSR2013029, the Education Department of Jiangsu under Grant No. JHB20121, and the Priority Academic Program Development of Jiangsu Higher Education Institutions under the Grant No. CE02138 is gratefully acknowledged.
References
 L. A. Caro, J. R. MartiVargas, and P. Serna, “Prestressed losses evaluation in prestressed concrete prismatic specimens,” Engineering Structures, vol. 48, pp. 704–715, 2013. View at: Publisher Site  Google Scholar
 E. P. Steinberg, “Probabilistic assessment of prestress loss in pretensioned prestressed concrete,” PCI Journal, vol. 40, no. 6, pp. 76–85, 1995. View at: Publisher Site  Google Scholar
 R. Malm and H. Sundquist, “Timedependent analyses of segmentally constructed balanced cantilever bridges,” Engineering Structures, vol. 32, no. 4, pp. 1038–1045, 2010. View at: Publisher Site  Google Scholar
 Z. P. Bazant, Q. Yu, and G. Li, “Excessive longtime deflections of prestressed box girders, I: recordspan bridge in Palau and other paradigms,” Journal of Structural Engineering, vol. 138, no. 6, pp. 676–686, 2012. View at: Publisher Site  Google Scholar
 T. Guo, R. Sause, D. M. Frangopol, and A. Q. Li, “Timedependent reliability of PSC boxgirder bridge considering creep, shrinkage and corrosion,” Journal of Bridge Engineering, vol. 16, no. 1, pp. 29–43, 2011. View at: Publisher Site  Google Scholar
 T. Guo, D. M. Frangopol, D. Z. Han, and Y. W. Chen, “Probabilistic assessment of deteriorating prestressed concrete boxgirder bridges under increased vehicle loads and aggressive environment,” Journal of Performance of Constructed Facilities, vol. 25, no. 6, pp. 564–576, 2011. View at: Publisher Site  Google Scholar
 J. J. Roller, H. G. Russell, R. N. Bruce, and B. T. Martin, “Longterm performance of prestressed, pretensioned high strength concrete bridge girders,” PCI Journal, vol. 40, no. 6, pp. 48–59, 1995. View at: Publisher Site  Google Scholar
 T. M. Ahlborn, C. E. French, and C. K. Shield, “High strength concrete prestressed bridge girders: long term and flexural behavior,” Minnesota Department of Transportation, Saint Paul, MN, USA, 2000, Report MN/RC–2000–32. View at: Google Scholar
 M. Lewis and V. M. Karbhari, “Experimental verification of the influence of timedependent material properties on longterm bridge characteristics,” SSRP, vol. 6, p. 20, 2006. View at: Google Scholar
 P. J. Barr, B. M. Kukay, and M. W. Halling, “Comparison of prestress losses for a prestress concrete bridge made with high performance concrete,” Journal of Bridge Engineering, vol. 13, no. 5, pp. 468–475, 2008. View at: Publisher Site  Google Scholar
 Z. Y. Liang, “Monitoring of a prestressed high performance concrete bridge using an embedded optical fiber sensor system,” New Mexico State University, Las Cruces, NM, USA, 2008, Ph.D. Dissertation. View at: Google Scholar
 C. G. Lan, Z. Zhou, and J. P. Ou, “Monitoring of structural prestress loss in RC beams by inner distributed Brillouin and fiber Bragg grating sensors on a single optical fiber,” Structural Control and Health Monitoring, vol. 21, no. 3, pp. 317–330, 2014. View at: Publisher Site  Google Scholar
 X. M. Jin, Y. L. Du, B. C. Sun, and X. Y. Li, “Study on optical fiber sensing technology for longterm real time monitoring of prestressed reinforcement,” China Safety Science Journal, vol. 15, no. 11, pp. 108–112, 2005. View at: Google Scholar
 G. D. Wang and M. L. Wang, “The utilities of Ushape EM sensors in stress monitoring,” Structural Engineering and Mechanics, vol. 17, no. 3–4, pp. 291–302, 2004. View at: Publisher Site  Google Scholar
 M. L. Wang, G. Wang, and Y. Zhao, “Application of EM stress sensors in large steel cables,” Sensing Issues in Civil Structural Health Monitoring, pp. 145–154, 2005. View at: Publisher Site  Google Scholar
 T. Guo, Z. H. Chen, S. Lu et al., “Monitoring and analysis of longterm prestress losses in posttensioned concrete beams,” Measurement, vol. 122, pp. 573–581, 2017. View at: Google Scholar
 Z. P. Bazant, “Prediction of concrete creep effects using ageadjusted effective modulus method,” ACI Journal Proceedings, vol. 69, no. 4, pp. 212–217, 1972. View at: Publisher Site  Google Scholar
 Z. P. Bazant and S. Baweja, “Justification and refinement of model B3 for concrete creep and shrinkage1, statistics and sensitivity,” Materials and Structures, vol. 28, no. 7, pp. 415–430, 1995. View at: Publisher Site  Google Scholar
 T. J. Trevino and A. Ghali, “Relaxation of steel in prestressed concrete,” PCI Journal, vol. 30, no. 5, pp. 82–94, 1985. View at: Publisher Site  Google Scholar
 S. A. Youakim, A. Ghali, S. E. Hida, and V. M. Karbhari, “Prediction of longterm prestress losses,” PCI journal, vol. 52, no. 2, pp. 116–130, 2007. View at: Publisher Site  Google Scholar
 M. F. Granata, P. Margiotta, and M. Arici, “Simplified procedure for evaluating the effects of creep and shrinkage on prestressed concrete girder bridges and the application of European and North American prediction models,” Journal of Bridge Engineering, vol. 18, no. 12, pp. 1281–1297, 2013. View at: Publisher Site  Google Scholar
 B. S. Choi, A. Scanlon, and P. A. Johnson, “Monte Carlo simulation of immediate and timedependent deflections of reinforced concrete beams and slabs,” ACI Structural Journal, vol. 101, no. 5, pp. 633–641, 2004. View at: Publisher Site  Google Scholar
 D. Diamantidis, H. O. Madsen, and R. Rackwitz, “On the variability of the creep coefficient of structural concrete,” Materials and Structures, vol. 17, no. 4, pp. 321–328, 1984. View at: Publisher Site  Google Scholar
 E. H. Khor, D. V. Rosowsky, and M. G. Stewart, “Probabilistic analysis of time dependent deflections of RC flexural members,” Computers & Structures, vol. 79, no. 16, pp. 1461–1472, 2001. View at: Publisher Site  Google Scholar
 T. Guo and Z. H. Chen, “Deflection control of longspan PSC boxgirder bridge based on field monitoring and probabilistic FEA,” Journal of Performance of Constructed Facilities, vol. 30, no. 6, Article ID 04016053, 2016. View at: Publisher Site  Google Scholar
 GBT502831999, Unified Standard for Reliability Design of Highway Engineering Structures, China Planning Press, Beijing, China, 1999, in Chinese.
 C. G. Gilbertson and T. M. Ahlborn, “A probabilistic comparison of prestress loss methods in prestressed concrete beams,” PCI Journal, vol. 49, no. 5, pp. 52–69, 2004. View at: Publisher Site  Google Scholar
 T. Guo, T. Liu, and A. Q. Li, “Deflection reliability analysis of PSC boxgirder bridge under highspeed railway loads,” Advances in Structural Engineering, vol. 15, no. 11, pp. 2001–2011, 2012. View at: Publisher Site  Google Scholar
 Y. H. Li, W. G. Bao, X. W. Guo et al., Highway Bridge Structure Reliability and Probability Limit Sate Design, China Communications Press, Beijing, China, 1997, in Chinese.
 A. S. ALHarthy and D. M. Frangopol, “Reliability assessment of prestressed concrete beams,” Journal of Structural Engineering, vol. 120, no. 1, pp. 180–199, 1994. View at: Google Scholar
Copyright
Copyright © 2018 Zheheng Chen and Shanwen Zhang. 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.