Research Article  Open Access
Changcheng Zhou, Leilei Zhao, Yuewei Yu, Fuxing Yang, Song Wang, "Kinetodynamic Model and Simplified Model of XType Seat System with an Integrated Spring Damper", Shock and Vibration, vol. 2018, Article ID 4962436, 10 pages, 2018. https://doi.org/10.1155/2018/4962436
Kinetodynamic Model and Simplified Model of XType Seat System with an Integrated Spring Damper
Abstract
A simplified analytical model with transformation coefficients of Xtype seat system with an integrated spring damper for control strategies development is proposed on the basis of a kinetodynamic model. Firstly, based on a commercial seat of trucks, the relationship between the suspension support force and the spring force was created by using virtual work principle. The analytical formulae of the equivalent stiffness and the stiffness transformation coefficient were deduced. Based on the principle of conservation of energy, the analytical formulae of the equivalent damping coefficient and the damping transformation coefficient were deduced. Then, the motion equation of the simplified model was created. Secondly, the nonlinear dynamic equation of a complex seat model including the kinematic characteristics was established. Thirdly, the road test was conducted using a heavy truck to collect the seat vibration signals. Finally, the simplified model was validated by the tested data and compared with the complex model. The results show that the accuracy of the simplified model is acceptable. Moreover, the influence laws of kinematic parameters on and were revealed. The proposed simplified model provides an accurate and efficient tool for designing controllable seat suspension system that minimizes a necessary tuning process.
1. Introduction
Trucks play an important role in transportation and logistics industry. As products of human civilization and technology progress, they have made great contributions to the social development [1]. Because of the complex and changeable road conditions, the vibration caused by the road excitation and the power assembly not only leads to the damage of the vehicle parts but also reduces the ride comfort [2]. More seriously, it damages the driver’s physical and mental health and endangers the driving safety [3].
With the rapid development of the truck industry, improving ride comfort has become a major concern of modern truck designers [4–6]. The seat suspension system is an important part of modern trucks [7]. The seat suspension system directly affects the NVH (noise, vibration, and harshness) performance [8]. Thus, it becomes a key link to improve ride comfort of trucks [9].
The seat suspension system can be mainly divided into three categories: passive seat suspensions [10], semiactive seat suspensions [11], and active seat suspensions [12]. The vibration isolation capability of passive seat suspensions is limited. In order to improve drivers’ ride comfort as much as possible, the controllable seat suspension is becoming a research hotspot [13].
The development and implementation of most of the control strategies depend on the dynamic model of the seat suspension system. Thus, the dynamic model of the seat suspension system is an important foundation for control strategies development and electric control. To facilitate control strategies development, the seat dynamic model should be simplified as much as possible under the premise of meeting the model accuracy. Generally speaking, it is necessary to compromise the model accuracy and the solution time.
The existing dynamic models considering the seat suspension structure are mainly used for the vibration characteristics analysis. For example, Feng and Hu optimized the vibration acceleration transmissibility for seat suspension system with a parallel mechanism [14]. Wang et al. made a hierarchical optimization for scissor seat suspension [15]. Shangguan et al. optimized vehiclespecific seat suspension systems using a seat kinetodynamic model [16]. These complex models are useful for analyzing the internal and fundamental characteristics of seat mechanisms. They also benefit parameter analysis and optimization. However, from the point of view of controller design, their equations of motion are of no practical use because they cannot satisfy the realtimeliness for control strategies realization in practical engineering.
In most cases, the classical single DOF (degree of freedom) model of seat suspension system is used for control strategies development. The model is composed of one lumped mass, a linear spring, and a damper [17, 18]. In recent years, scholars have carried out related research on controllable seat suspensions based on the classical single DOF model. Moreover, some fruitful research results have been achieved. Do et al. designed a new adaptive fuzzy slidingmode controller to control a seat suspension [19]. Ning et al. researched an active seat suspension based on the TakagiSugeno fuzzy control [20]. Rajendiran et al. simulated the responses of seat suspension with PID and fuzzy logic controllers [21]. Ning et al. designed an active seat suspension for vibration control of heavyduty vehicles [22]. Although the classical single DOF model can be used for exploring various kinds of control strategies, it does not establish the relationship between the model equivalent parameters and the suspension structure. Thus, it cannot be directly used in the design of the seat suspension controller for actual vehicles [23, 24]. Actually, the seat suspension kinematic mechanism and the layout form have great influences on the dynamic responses of the seat system [25, 26]. Zhang and Zhao researched the nonlinear stiffness and vibration isolation characteristics of scissorlike structure with full types [27]; however, they did not provide the analytical formula of the equivalent damping coefficient of the seat damper. Available studies provide useful references for the design, control, and modeling of seat suspension system. However, there are few studies on providing an accurate simple model of seat suspension system with complex scissor mechanisms for control strategies development. Creating a simple seat model with reasonable accuracy is the primary motivation of this paper.
To facilitate control strategies development, this paper proposed a simplified analytical dynamic model with transformation coefficients of seat suspension system. In Section 2, according to a commercial seat of trucks, the nonlinear relationship between the suspension support force and the spring force was created by using virtual work principle. Based on the relationship, the analytical formula of the equivalent stiffness and the stiffness transformation coefficient was deduced. Based on the principle of conservation of energy, the analytical formula of the equivalent damping coefficient was deduced, and the damping transformation coefficient was proposed. The motion equation of the simplified model was created. In Section 3, to facilitate comparison, the nonlinear dynamic equation of the seat kinetodynamic model is established by using Lagrange modeling method. In Section 4, the simplified model was validated by the tested data and compared with the kinetodynamic model. In Section 5, the influence laws of kinematic parameters on and were revealed.
2. Development of a Simplified Analytical Model of Seat Suspension System
2.1. Physical Structure and Model of the Seat Suspension System
The seat suspension with an integrated spring damper has been gradually applied in trucks in recent years because of its compact structure [28]. The integrated spring damper is composed of a hydraulic damper, a spiral spring, and two bushings. The hydraulic damper and the spiral spring are connected in parallel. On the basis of the physical structure of Xtype seat, the kinetodynamic model is created as shown in Figure 1(a). The moving parts contains seven rigids: the damper rod, the damper tube, the guide wheel at point P, the guide wheel at point R, the seat pan, the linkage PS, and the linkage OR. The coordinate system is shown in Figure 1(a) and its origin is at S point. In the model, the lengths of the linkages RQ and PS are equal, so denote as . For simplicity, denote the lengths of AS, BR, and BQ as , , and , respectively. The lengths and are called as the lower installation distance and the upper installation distance, respectively. The angle between the linkage PS and the seat floor is . Table 1 shows the topological relations among moving parts. This kinetodynamic model is rather complex, and it is inconvenient for control strategies development of the seat suspension. Thus, in Figure 1(b) a simplified analytical model with single DOF is proposed, which includes three parameters: is the effective mass of the driver body and the seat pan supported by the seat suspension, is the vertical equivalent stiffness of the seat suspension, is the vertical equivalent damping coefficient of the seat suspension. For conveniently translating the complex model to the simplified DOF model, the key problem is how to obtain the analytical formulae of and . In the following, the kinetodynamic model is denoted as the complex model. The simplified model is the equivalent model.
(a)
(b)

2.2. The Analytical Formula of the Equivalent Stiffness
If the seat pan produces an upward displacement from the static equilibrium position and , the point moves to the new point . In this condition, . Thus, is defined as the vertical deformation of the seat suspension system. According to the geometric relationship, can be expressed as
Substituting and in Equation (1), we obtain the following:
From Equation (2), we obtain the following:
At the static equilibrium position, the distance between the spring installation points can be expressed as
Substituting , , , and in Equation (4), we obtain the following:
Let , , and , and Equation (5) can be expressed as
When the spring is deformed from the static equilibrium position, the distance between the spring installation points can be expressed as
Substituting , , , and in Equation (7), we obtain the following:
Equation (8) can be further expressed as
The deformation of the spring can be expressed as
According to the principle of virtual work, the following can be obtained:where is the vertical external force acting on the seat pan, is the spring force, and is the virtual displacement.
From Equation (10), can be expressed as
Substituting Equation (12) in Equation (11), we obtain the following:
Because the elastic element of the seat suspension is a spiral spring, thus, the spring force can be expressed aswhere is the precompression deformation of the spring and is the spring stiffness.
Substituting Equation (14) in Equation (13), the force can be further expressed as
The vertical equivalent stiffness of the seat suspension can be solved bywhere is the vertical elastic restoring force of the seat suspension and , is the effective mass of the driver body supported by the seat suspension, and is the effective mass of the seat pan supported by the seat suspension.
According to Equations (15) and (16), can be expressed aswhere the stiffness transformation coefficient can be expressed as
In addition, at the static equilibrium position, let . That is
According to Equations (15) and (19) and , we obtain the following:
From Equation (20), at the static equilibrium position, the precompression deformation can be determined by
2.3. The Analytical Formula of the Equivalent Damping Coefficient
According to the conservation of energy, the energy dissipated by the damper in simplified model is equal to that in the kinetodynamic model. Thus,where is the damper damping.
According to the chain derivation rule, Equation (22) can be expressed as
From Equation (23), we obtain the following:
From Equation (10), can be expressed as
Substituting Equation (25) in Equation (24), the vertical equivalent damping coefficient of the seat suspension can be expressed aswhere the damping transformation coefficient can be expressed as
2.4. The Motion Equation of the Simplified Model
Substituting in Equations (17) and (26), the vertical equivalent stiffness and the vertical equivalent damping coefficient at the static equilibrium position can be obtained. They are denoted as and , respectively. At the static equilibrium position, the stiffness transformation coefficient and the damping transformation coefficient are denoted as and , respectively. Thus, and can be further, respectively, expressed as
Using Newton’s second law, the motion equation of the simplified model can be written as
3. Motion Equation of the Complex Model
The vibration equation of the seat system is created using the Lagrange modeling method in the section. Under the input excitation of the seat floor, the angle variation is , which is designated as the generalized coordinate for the kinetodynamic model [28].
The Lagrange equation of the seat system can be expressed as
where, , , and are the kinetic energy, the potential energy, and the dissipated energy of the seat system, respectively.
From (30), we obtain the nonlinear motion equation of the kinetodynamic model:
4. Model Validation and Comparison
In order to verify the effectiveness of the simplified model, the numerical simulations and comparison of the vertical acceleration responses were carried out in this section. The parameters of the seat system are as shown in Table 2. The driver mass is 75 kg.

4.1. Verification of Random Vibration Response
To obtain the vibration signals for simulation and validation, according to the national standard GB/T 49702009, the road test was carried out and the vibration signals were collected. Based on the complex model and the simplified model, the values of the seat pan vertical frequencyweighted RMS acceleration were calculated using the measured seat base accelerations at different speeds as inputs. The calculated results are shown in Table 3. In order to make it clear to observe the differences, a histogram is used to provide a comparison of the calculated results, as shown in Figure 2.

From Figure 2, it can be seen that the values of the seat pan vertical frequencyweighted RMS acceleration calculated from the tested data are larger than those calculated from the simulated data based on the complex model and the simplified model. For the vehicle at lower speeds 50, 60, and 70 km/h, the values of calculated from the complex model and the simplified model are the same, respectively. However, for the vehicle at higher speeds 80 and 90 km/h, the values of calculated from the complex model are slightly larger than those from the simplified model. The maximum of the relative error between the simulated values of calculated from the complex model and the tested data is 8.3%. The maximum of the relative error between the simulated values of calculated from the simplified model and the tested data is 10.4%. The results show that the two models can reflect the random acceleration response of the seat pan, and the accuracy of the simplified model is slightly lower than that of the complex model. From the point of view of controller design, the accuracy of the simplified model can satisfy the control strategies realization in practical engineering.
4.2. Verification and Comparison of the Vibration Transmissibility
4.2.1. Verification of the Vibration Transmissibility
In this section, the transmissibility of the seat system was measured through the test of sine wave sweep. Before the test, the seat was fixed on the test bench. Moreover, the driver was replaced by a weight block with 55 kg which is calculated by (8 + 75) × 9.8 × 0.73. On the basis of the frequency sweep method, the transmissibility of the seat system was measured. In addition, the transmissibility of the complex model was calculated based on the frequency sweep method, and the transmissibility of the simplified model was analytically calculated using the transfer function from to . A comparison of the transmissibility curves is shown in Figure 3. Based on Figure 3, the values of the natural frequency and the maximal transmissibility were extracted. A comparison of the two characteristic parameters of the seat system is shown in Table 4.

From Figure 3, it can be seen that the simulation curves of the transmissibility coincides with that of the measured curve. There is a certain deviation in the frequency bands 0.5∼1.5 Hz and 8.0∼10.0 Hz. The transmissibility curve calculated from the simplified model can reflect the real transmission characteristics of the seat suspension system. Table 4 shows that the tested is the same as the simulated ; the absolute deviation between the tested and the simulated is 0.01, and the relative deviation is 1.0%. The results further prove that the simplified model is effective. The main reason for the deviations is that the friction and gaps between the moving parts of the scissor mechanism are ignored in the simplified model.
4.2.2. Comparison of the Vibration Transmissibility Curves at Different Values of
Figure 4 provides a comparison of the vibration transmissibility curves at different values of with the static equilibrium position unchanged. Where = 55, 65, 75, and 85 kg. From Figure 4, it can be seen that the transmissibility curves calculated from the simplified model are in good agreement with those calculated from the complex model at different values of . Moreover, as the value of the mass increases, the natural frequency decreases, however, the maximal transmissibility increases.
4.3. Comparison of Harmonic Response
The vibration response analysis of the seat system under harmonic excitation is one of the effective methods to investigate the differences of the models and the real system. Figure 5 provides a comparison of the seat pan acceleration responses for the seat system under harmonic excitations with the amplitude = 10 mm and the excitation frequency = 1.0 and 2.0 Hz. From Figure 5, it can be seen that the simplified model can be considered suitable for reproducing the real seat acceleration response.
(a)
(b)
5. Influence Analysis of Kinematic Parameters on ρ_{k}, ρ_{c}, K_{e}, and C_{e}
In this section, the influence laws of kinematic parameters on the stiffness transformation coefficient , the damping transformation coefficient , the equivalent stiffness , and the equivalent damping were revealed on the basis of the values for the seat system in Table 2.
The influences of the kinematic parameters and on the stiffness transformation coefficient and the equivalent stiffness are shown in Figure 6 and Figure 7, respectively. The influences of the kinematic parameters and on the damping transformation coefficient and the equivalent damping are shown in Figure 8 and Figure 9, respectively.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
Figure 6(a) shows that the stiffness transformation coefficient approximately linearly increases with the increase of the suspension deformation . Moreover, with the increase of the kinematic parameters , the stiffness transformation coefficient decreases. Figure 6(b) depicts the same rules of the kinematic parameters and on the equivalent stiffness . Figure 7 further illustrates the above rules. From Figure 8, it can be seen that the damping transformation coefficient also approximately linearly increases with the increase of the suspension deformation . The larger the values of the kinematic parameters and are, the smaller the value of the damping transformation coefficient is. Similar change characteristics are also observed to occur in Figure 9.
6. Conclusions
For the seat suspension system, the development and implementation of most of the control strategies depend on its dynamic model. The primary motivation of this paper is to create a simple seat model with reasonable accuracy. Thus, a simplified analytical model with transformation coefficients of seat suspension with a complex scissor mechanism for control strategies development is proposed. Firstly, the analytical formulae of the equivalent stiffness and the stiffness transformation coefficient were deduced. The analytical formulae of the equivalent damping coefficient and the damping transformation coefficient were also deduced. Then, the motion equation of the simplified model was created. Secondly, the nonlinear dynamic equation of a complex seat model including the kinematic characteristics is established. Finally, the simplified model was validated by the tested data and compared with the complex model and the tested data. Moreover, the influence laws of kinematic parameters on and were revealed.
(1) Through the accuracy verification of random vibration response, it can be seen that the accuracy of the simplified model can satisfy the control strategies realization in practical engineering. (2) By comparison of the vibration transmissibility, the results show that the transmissibility curve calculated from the simplified model can reflect the real transmission characteristics of the seat suspension system. (3) By comparison of harmonic response, the results show that the simplified model can be considered suitable for reproducing the real seat acceleration response. (4) Both the stiffness transformation coefficient and the damping transformation coefficient approximately linearly increase with the increase of the suspension deformation . With the increase of the kinematic parameters and , both and decrease.
The proposed model provides an accurate and efficient tool for designing controllable seat suspension system that minimizes a necessary tuning process. In the following study, the nonlinear stiffness transformation and the asymmetric damping transformation would be an interesting topic about getting reasonable simple models for control strategies development.
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 no conflicts of interest.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (51575325) and BUPT Excellent Ph.D. Students Foundation (CX2016206).
References
 B. B. Du, P. L. Bigelow, R. P. Wells, P. Hall, and H. W. Davies, “The impact of different seats and wholebody vibration exposures on truck driver vigilance and discomfort,” Ergonomics, vol. 61, no. 4, pp. 528–537, 2018. View at: Publisher Site  Google Scholar
 N. Raffler, J. Rissler, R. Ellegast, C. Schikowsky, and T. Kraus, “Combined exposures of wholebody vibration and awkward posture: a cross sectional investigation among occupational drivers by means of simultaneous field measurements,” Ergonomics, vol. 60, no. 11, pp. 1–12, 2017. View at: Publisher Site  Google Scholar
 L. L. Zhao, C. C. Zhou, Y. W. Yu, and F. X. Yang, “Hybrid modelling and damping collaborative optimisation of Fivesuspensions for coupling driverseatcab system,” Vehicle System Dynamics, vol. 54, no. 5, pp. 667–688, 2016. View at: Publisher Site  Google Scholar
 F. Diba and E. Esmailzadeh, “Development of hybrid electric heavyduty truck with selfpropelled trailer,” International Journal of Heavy Vehicle Systems, vol. 25, no. 2, pp. 203–222, 2018. View at: Publisher Site  Google Scholar
 G. Wang, K. Li, J. Yang, and C. Liang, “Simulation and vibration isolation optimization on the cab mount of a commercial vehicle,” Automotive Engineering, vol. 39, pp. 1081–1086, 2017. View at: Google Scholar
 A. Gupta and V. Rastogi, “Computational modelling and simulation of commercial truck with roadtyre interactions to evaluate ride comfort,” International Journal of Heavy Vehicle Systems, vol. 25, no. 2, p. 133, 2018. View at: Publisher Site  Google Scholar
 X. Wang, F. Bi, and H. Du, “Reduction of low frequency vibration of truck driver and seating system through system parameter identification, sensitivity analysis and active control,” Mechanical Systems and Signal Processing, vol. 105, pp. 16–35, 2018. View at: Publisher Site  Google Scholar
 F. X. Yang, L. L. Zhao, Y. W. Yu, and C. C. Zhou, “Analytical description of ride comfort and optimal damping of cushionsuspension for wheeldrive electric vehicles,” International Journal of Automotive Technology, vol. 18, pp. 1121–1129, 2017. View at: Publisher Site  Google Scholar
 L. L. Zhao, Y. W. Yu, C. C. Zhou, and S. Mao, “Test verification of a seatcabin system vertical dynamic model for highgrade fork lift trucks,” International Journal of Modeling, Simulation, and Scientific Computing, vol. 9, Article ID 1850003, 2018. View at: Publisher Site  Google Scholar
 L. L. Zhao, Y. W. Yu, and C. C. Zhou, “Comparative research on optimal damping matching of seat system for an offhighway dump truck,” International Journal of Engineering, Transactions B: Applications, vol. 31, no. 2, pp. 204–211, 2018. View at: Publisher Site  Google Scholar
 D. Ning, S. Sun, H. Du, and W. Li, “Integrated active and semiactive control for seat suspension of a heavy duty vehicle,” Journal of Intelligent Material Systems and Structures, vol. 29, pp. 91–100, 2018. View at: Publisher Site  Google Scholar
 J. E. D. Ekoru, O. T. C. Nyandoro, and T. Chingozha, “An UdwadiaKalaba Equation based active seat suspension controller for mining dump trucks,” IFACPapersOnLine, vol. 49, pp. 184–189, 2016. View at: Publisher Site  Google Scholar
 A. Alfadhli, J. Darling, and A. J. Hillis, “The control of an active seat with vehicle suspension preview information,” Journal of Vibration and Control, vol. 24, no. 8, pp. 1412–1426, 2018. View at: Publisher Site  Google Scholar
 X. L. Feng and J. Hu, “Numerical optimization of vibration acceleration transmissibility for seat suspension system in vehicles,” Journal of Vibroengineering, vol. 18, pp. 1115–1132, 2016. View at: Google Scholar
 C. Wang, X. Zhang, K. Guo, J. Lv, and Y. Yang, “Hierarchical optimisation on scissor seat suspension characteristic and structure,” Vehicle System Dynamics, vol. 54, no. 11, pp. 1538–1553, 2016. View at: Publisher Site  Google Scholar
 W. B. Shangguan, Y. Shui, and S. Rakheja, “Kinetodynamic design optimisation for vehiclespecific seatsuspension systems,” Vehicle System Dynamics, vol. 55, no. 11, pp. 1–22, 2017. View at: Publisher Site  Google Scholar
 L. L. Zhao, C. C. Zhou, and Y. W. Yu, “Identification and optimisation of surrogate model parameters of the bus driver’s seat suspension,” International Journal of Vehicle Systems Modelling and Testing, vol. 11, no. 4, pp. 330–342, 2016. View at: Publisher Site  Google Scholar
 D. X. Phu, S. M. Choi, and S. B. Choi, “A new adaptive hybrid controller for vibration control of a vehicle seat suspension featuring MR damper,” Journal of Vibration and Control, vol. 23, no. 20, pp. 3392–3413, 2016. View at: Publisher Site  Google Scholar
 P. X. Do, S. B. Choi, Y. S. Lee, and M. S. Han, “Vibration control of a vehicle’s seat suspension featuring a magnetorheological damper based on a new adaptive fuzzy slidingmode controller,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 230, no. 4, pp. 437–458, 2016. View at: Publisher Site  Google Scholar
 D. Ning, S. Sun, F. Zhang, H. Du, and W. Li, “Disturbance observer based TakagiSugeno fuzzy control for an active seat suspension,” Mechanical Systems and Signal Processing, vol. 93, pp. 515–530, 2017. View at: Publisher Site  Google Scholar
 S. Rajendiran and P. Lakshmi, “Simulation of PID and fuzzy logic controller for integrated seat suspension of a quarter car with driver model for different road profiles,” Journal of Mechanical Science and Technology, vol. 30, no. 10, pp. 4565–4570, 2016. View at: Publisher Site  Google Scholar
 D. Ning, S. Sun, F. Zhang, H. Du, and W. Li, “An active seat suspension design for vibration control of heavyduty vehicles,” Journal of Low Frequency Noise Vibration and Active Control, vol. 35, no. 4, pp. 264–278, 2016. View at: Publisher Site  Google Scholar
 L. L. Zhao, C. C. Zhou, and Y. W. Yu, “Hybrid modeling of driver’s seatcushion coupled system for metropolitan bus,” Journal of Low Frequency Noise, Vibration and Active Control, vol. 36, no. 3, pp. 214–226, 2017. View at: Publisher Site  Google Scholar
 L. L. Zhao, Y. W. Yu, C. C. Zhou, S. F. Wang, and F. X. Yang, “A new hydraulic semiactive suspension based on road statistical properties and its road identification,” Applied Science, vol. 8, no. 5, pp. 1–13, 2018. View at: Publisher Site  Google Scholar
 R. D. Nandurkar, A. Mohanty, and P. Barath, “Design and simulation of automotive seat height adjuster driving mechanism based on multibody dynamics,” Applied Mechanics and Materials, vol. 592–594, pp. 2282–2286, 2014. View at: Publisher Site  Google Scholar
 K. Hill and A. Dhingra, “Modeling, analysis and optimization of a scissors linkage seat suspension,” Engineering Optimization, vol. 35, no. 4, pp. 341–357, 2003. View at: Publisher Site  Google Scholar
 W. Zhang and J. Zhao, “Analysis on nonlinear stiffness and vibration isolation performance of scissorlike structure with full types,” Nonlinear Dynamics, vol. 86, no. 1, pp. 1–20, 2016. View at: Publisher Site  Google Scholar
 L. L. Zhao, F. X. Yang, Y. W. Yu et al., “Influence of shock absorber asymmetry on vibration responses of seatoccupant system under impact condition,” in Proceedings of International Conference on Modelling, Identification and Control (ICMIC 2017), pp. 394–399, IEEE Xplore, Kunming, China, July 2017. View at: Google Scholar
Copyright
Copyright © 2018 Changcheng Zhou 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.