A novel semi-active suspension design based on decoupling skyhook control

Yao, Jialing; Taheri, Saied; Tian, Songmei; Zhang, Zhongnan; Shen, Liang

Journal of Vibroengineering

Browse Journal

Submit article

Published: 15 May 2014

Check for updates

A novel semi-active suspension design based on decoupling skyhook control

Jialing Yao¹

Saied Taheri²

Songmei Tian³

Zhongnan Zhang⁴

Liang Shen⁵

^{1, 4, 5}College of Automobile and Traffic Engineering, Nanjing Forestry University, Nanjing, 210037, China

^{2, 1}Center for Tire Research, Mechanical Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA 24060, USA

³Research Institute of Automotive Engineering, Nanjing Automobile Group Limited Company, Nanjing, 210028, China

Corresponding Author:

Jialing Yao

Cite the article Download PDF

Downloads 1311

WoS Core Citations 15

Abstract

A semi-active suspension design based on the traditional method of skyhook control is not capable of effectively controlling the attitude of the vehicle. However, an innovative approach called decoupling skyhook control allows the attitude of the vehicle body and its vibration characteristics to be effectively controlled. In this paper, a new decoupling skyhook controller for semi-active suspension is presented. Vehicle body motions in the three directions of vertical, pitch, and roll have been adopted to develop three skyhook controllers and directly control the vehicle body attitude. Furthermore, three orientation skyhook control forces are converted into actual damping forces of four adjustable dampers through the input decoupling transformation. The simulation results show that the developed controller is more effective than the traditional skyhook control in improving ride comfort.

1. Introduction

Since the skyhook control scheme was introduced in the early 1970s by D. Karnopp, this method has been used in many vibration isolation applications [1]. In active and semi-active suspension control, the skyhook control scheme is an extremely effective control strategy and has been used in some commercial applications. Moreover, the skyhook scheme is often used as the reference model for other control strategies because of its simplicity and its prominent effect in improving ride comfort [2-4].

To develop the semi-active control algorithm for a full-vehicle system, two different methodologies can be used [5]. One methodology uses the quarter car model to develop the control algorithm, which is then implemented in a full vehicle model for controlling and analysis [6-8]. However, this methodology does not directly consider the pitch and roll motions of the vehicle. Furthermore, four accelerometers on the four corner points of the sprung mass are usually used to design controllers. As the aim of the suspension is to control the acceleration felt by the passengers, this acceleration should be considered at the vehicle’s center of gravity and not over the wheels. Thus, the control of the body posture may be inadequate.

Another methodology uses a full-vehicle model to directly control both the vertical motion (heave) and the angular motion (pitch and roll) of the vehicle body [9, 10]. The advantage of this methodology is that a control strategy can be designed that effectively controls the attitude of the sprung mass and suppresses the vibration of the suspension. In [9], three fuzzy controllers are designed to produce three forces $f_{z}$ , $f_{θ}$ and $f_{φ}$ to suppress the heave, pitch, and roll movements. The input variables for the controllers are velocities and accelerations in three orientations. Then, the three forces are decoupled into four forces generated by active force generators. A double-loop control scheme is presented to control the active suspension in [10]. The inner control loop adopts the active filtering of the spring and damping coefficient to design the ride controller, while the outer control loop maintains load-leveling and load distribution during vehicle maneuvers using skyhook scheme. Finally, three forces, which control the heave, pitch, and roll movements, are also decoupled into the four active generators’ forces. Reference [11] also proposes a decoupled quarter car active suspension model that introduces a new abstract variable denoted by $γ$ . For semi-active suspension, a coupling design using skyhook damping is mentioned in [12]. However, no further explanation is given. According to the theory of skyhook damping forces that damp heave and angular motion proposed in [1], this study develops semi-active suspension controllers through a direct use of the skyhook damping. First, three ideal skyhook controllers are designed to command forces that damp the vehicle body’s heave, pitch, and roll motions. Next, the three orientation forces are decoupled into damping forces of four adjustable dampers, which are constrained by passivity. The developed strategy is simple and effective.

2. Full vehicle model

The vehicle is represented as a linearized seven-degree-of-freedom system, as shown in Fig. 1. This model is similar to the model used in [13]. It consists of a sprung mass (car body) and four unsprung masses. For simplicity, all pitch and roll angles are assumed to be small. The suspension is also modeled as linear spring and damper, while the tires are modeled as simple linear springs without damping. The damper is divided into two parts: the damper with a fixed damping coefficient and the damper with a variable damping coefficient.

Newton’s second law applied to the system results in the following differential equations of motion:

1

m_{s} \ddot{z} = - m_{s} g - K_{f l} (z_{f l} - z_{u f l}) - K_{f r} (z_{f r} - z_{u f r}) - K_{r l} (z_{r l} - z_{u r l})

- K_{r r} (z_{r l} - z_{u r r}) + \dots - C_{f l} ({\dot{z}}_{f l} - {\dot{z}}_{u f l}) - C_{f r} ({\dot{z}}_{f r} - {\dot{z}}_{u f r}) - C_{r f} ({\dot{z}}_{r f} - {\dot{z}}_{u r f})

- C_{r r} ({\dot{z}}_{r r} - {\dot{z}}_{u r r}) + f_{f l} + f_{f r} + f_{r l} + f_{r r},

I_{y} \ddot{θ} = K_{f l} a (z_{f l} - z_{u f l}) + K_{f r} a (z_{f r} - z_{u f r}) - K_{r l} b (z_{r l} - z_{u r l}) - K_{r r} b (z_{r l} - z_{u r r}) + . . . + C_{f l} a ({\dot{z}}_{f l} - {\dot{z}}_{u f l}) + C_{f r} a ({\dot{z}}_{f r} - {\dot{z}}_{u f r}) - C_{r l} b ({\dot{z}}_{r l} - {\dot{z}}_{u r l}) - C_{r r} b ({\dot{z}}_{r r} - {\dot{z}}_{u r r}) - f_{f l} a - f_{f r} a + f_{r l} b + f_{r r} b,

I_{x} \ddot{φ} = - 0.5 K_{f l} w (z_{f l} - z_{u f l}) + 0.5 K_{f r} w (z_{f r} - z_{u f r}) - 0.5 K_{r l} w (z_{r l} - z_{u r l})

+ 0.5 K_{r r} w (z_{r r} - z_{u r r}) + \dots - 0.5 C_{f l} w ({\dot{z}}_{f l} - {\dot{z}}_{u f l}) + 0.5 C_{f r} w ({\dot{z}}_{f r} - {\dot{z}}_{u f r}) - 0.5 C_{r l} w ({\dot{z}}_{r l} - {\dot{z}}_{u r l}) + 0.5 C_{r r} w ({\dot{z}}_{r r} - {\dot{z}}_{u r r}) + . . . 0.5 w f_{f l} - 0.5 w f_{f r} + 0.5 w f_{r l} - 0.5 w f_{r r},

m_{u f l} {\ddot{z}}_{u f l} = - m_{u f l} g + K_{f l} (z_{f l} - z_{u f l}) + C_{f l} ({\dot{z}}_{f l} - {\dot{z}}_{u f l}) - K_{u f r} (z_{u f l} - z_{r f l}) - f_{f l},

m_{u f r} {\ddot{z}}_{u f r} = - m_{u f r} g + K_{f r} (z_{f r} - z_{u f r}) + C_{f r} ({\dot{z}}_{f r} - {\dot{z}}_{u f r}) - K_{u f r} (z_{u f r} - z_{r f r}) - f_{f r},

m_{u r l} {\ddot{z}}_{u r l} = - m_{u r l} g + K_{r l} (z_{r l} - z_{u r l}) + C_{r l} ({\dot{z}}_{r l} - {\dot{z}}_{u r l}) - K_{u r l} (z_{u r l} - z_{r r l}) - f_{r l},

m_{u r r} {\ddot{z}}_{u r r} = - m_{u r r} g + K_{r r} (z_{r r} - z_{u r r}) + C_{r r} ({\dot{z}}_{r r} - {\dot{z}}_{u r r}) - K_{u r r} (z_{u r r} - z_{r r r}) - f_{r r} .

The kinematic equations representing the displacement at each corner are given as follows:

2

z_{f l} = z - a θ + 0.5 w φ,

z_{f r} = z - a θ - 0.5 w φ,

z_{r l} = z + b θ + 0.5 w φ,

z_{r l} = z + b θ - 0.5 w φ,

where: $z$ – heave displacement of the vehicle body’s mass center, $θ / φ$ – pitch/roll angle, $z_{f l}$ , $z_{f r}$ , $z_{r l}$ , $z_{r r}$ – heave displacement of the vehicle body at the suspension locations, $z_{u f l}$ , $z_{u f r}$ , $z_{u r l}$ , $z_{u r r}$ – heave displacement of the unsprung mass, $z_{r f l}$ , $z_{r f r}$ , $z_{r r l}$ , $z_{r r r}$ – excitation from the road surface, $f_{f l}$ , $f_{f r}$ , $f_{r l}$ , $f_{r r}$ – variable damping force, and $g$ – the acceleration due to gravity.

Fig. 1Vehicle model with seven degrees of freedom

The full-vehicle mode parameters were adopted from [8] and are summarized in Table 1.

Table 1Model parameters

Symbol	Description	Units
$m_{s}$	Sprung mass	1465 kg
$m_{u f l}, m_{u f r}, m_{u r l}, m_{u r r}$	Unsprung mass	40 kg
$K_{f l}, K_{f r}$	Front suspension stiffness	19960 N/m
$K_{r l}, K_{r r}$	Rear suspension stiffness	17500 N/m
$K_{u f l}, K_{u f r}, K_{u r l}, K_{u r r}$	Tire stiffness	175500 N/m
$C_{f l}, C_{f r}$	Fixed damping of the front dampers	258 N/m/s
$C_{r l}, C_{r r}$	Fixed damping of the rear dampers	324 N/m/s
$I_{y}$	Pitch moment of inertia	2460 kg/m²
$I_{x}$	Roll moment of inertia	460 kg/m²
$a$	Longitudinal distance from vehicle c.g. to front axle	1.4 m
$b$	Longitudinal distance from vehicle c.g. to rear axle	1.7 m
$w$	Track width	3 m

3. Input decoupling transformation

The actuators used to implement the semi-active suspension are continuously adjustable dampers (such as MR dampers), which are installed in all four corners of the vehicle. Each damper is independently controlled to improve the vehicle ride. The sprung mass acceleration at the body’s center of gravity should be considered rather than that over the wheels, because the suspension control target is to reduce the acceleration in the location of the occupants. Thus, the damping force of four dampers should be associated with the control forces for three sprung mass accelerations. In this project, the vehicle body posture is directly controlled, considering the heave, pitch, and roll motions, and the three sprung mass accelerations ( $\ddot{z}$ , $\ddot{θ}$ and $\ddot{φ}$ ) at the center of gravity are also controlled.

From Eq. (1), the equivalent forces for heave, pitch, and roll can be defined by:

3

f_{z} (t) = f_{f l} (t) + f_{f r} (t) + f_{r l} (t) + f_{r r} (t),

f_{θ} (t) = - a f_{f l} (t) - a f_{f r} (t) + b f_{r l} (t) + b f_{r r} (t),

f_{φ} (t) = 0.5 w (f_{f l} (t) - f_{f r} (t) + f_{r l} (t) - f_{r r} (t)),

and expressed in a matrix form:

4

[\begin{matrix} f_{z} \\ f_{θ} \\ f_{φ} \end{matrix}] = W [\begin{matrix} f_{f l} \\ f_{f r} \\ f_{r l} \\ f_{r r} \end{matrix}] = [\begin{matrix} 1 & 1 & 1 & 1 \\ - a & - a & b & b \\ 0.5 w & - 0.5 w & 0.5 w & - 0.5 w \end{matrix}] [\begin{matrix} f_{f l} \\ f_{f r} \\ f_{r l} \\ f_{r r} \end{matrix}] .

Eq. (4) takes the form ‘ $b = A x$ ’, where $b$ is the $m \times 1$ vector, $x$ is the $n \times 1$ vector, and $A$ is the $n \times m$ matrix. The pseudo right inverse can be obtained using $W = A^{T} {(A A^{+})}^{- 1}$ , $W A = I$ . Thus, Eq. (4) can be expressed as:

5

[\begin{matrix} f_{f l} \\ f_{f r} \\ f_{r l} \\ f_{r r} \end{matrix}] = [\begin{matrix} \frac{b}{2 (a + b)} & - \frac{1}{2 (a + b)} & \frac{1}{2 w} \\ \frac{b}{2 (a + b)} & - \frac{1}{2 (a + b)} & - \frac{1}{2 w} \\ \frac{a}{2 (a + b)} & \frac{1}{2 (a + b)} & \frac{1}{2 w} \\ \frac{a}{2 (a + b)} & \frac{1}{2 (a + b)} & - \frac{1}{2 w} \end{matrix}] [\begin{matrix} f_{z} \\ f_{θ} \\ f_{φ} \end{matrix}] .

By entering the decoupling transformation expressed by (5), the control for heave, pitch, and roll of the car body is converted into the control for the four dampers.

4. Skyhook control design

In a full-vehicle suspension system, the rigid body of the vehicle (i.e., the sprung mass) has three DOFs of heave, pitch, and roll, as HPR in the rest of this paper.

According to the theory of skyhook damping presented in [1], the ideal skyhook control forces in three orientations (i.e., heave, pitch, and roll) for the suspension system of a full car model can be expressed as:

6

F_{j} = C_{j} {\dot{Z}}_{j}, j = z, θ, φ,

where $C_{j}$ is the skyhook control gain, ${\dot{Z}}_{j}$ consists of the HPR velocities of the vehicle body, and $z$ , $θ$ and $φ$ represent the heave, pitch, and roll orientations, respectively.

We obtain an over-actuated system with three degrees of freedom controlled by four actuators. Thus, the three ideal skyhook control forces must be decoupled. Meanwhile, the actual adjustable dampers cannot produce active forces, so the semi-active suspension controller can only be conducted in an on-off control. Thus, the actual skyhook controller needs the velocities and relative suspension velocities of the sprung mass. Given the problem of measuring relative velocity in a practical application, velocities of sprung mass and relative velocities at each suspension corner are adopted to build a skyhook controller to achieve a direct input decoupling. In this decoupling skyhook control, seven sensors are needed: three sprung mass accelerometers (measuring $\ddot{z},$ $\ddot{θ}$ and $\ddot{φ}$ ) and four displacement sensors (measuring suspension relative displacement).

The semi-active skyhook control forces that are constrained by passivity can be expressed as:

7

f_{j} = \{\begin{array}{l} C_{j} \cdot {\dot{Z}}_{j}, & {\dot{Z}}_{j} {\dot{Z}}_{i} \geq 0, \\ 0, & otherwise, \end{array}

where $i = f l$ , $f r$ , $r l$ , $r r$ and $j = z$ , $θ$ and $φ$ .

${\dot{Z}}_{j}$ consists of the HPR velocities of the vehicle body, and ${\dot{Z}}_{i}$ of the relative velocities at each suspension corner.

From Eq. (2), the HPR velocities of the four corner points of the vehicle body can be obtained:

8

{\dot{z}}_{f l} = \dot{z} - a \dot{θ} + 0.5 w \dot{φ},

{\dot{z}}_{f r} = \dot{z} - a \dot{θ} - 0.5 w \dot{φ},

{\dot{z}}_{r l} = \dot{z} + b \dot{θ} + 0.5 w \dot{φ},

{\dot{z}}_{r r} = \dot{z} + b \dot{θ} - 0.5 w \dot{φ} .

Finally, the actual control force demanded from each damper can be calculated by Eqs. (5) and (7).

5. Simulation

The simulation focuses on the vehicle ride comfort and is representative of the vehicle body accelerations. To validate the decoupling skyhook control (DSKY) scheme proposed in this study, a passive suspension is also simulated with fixed damping forces acting on the four dampers. The passive suspension is treated as the baseline for DSKY control. Furthermore, a traditional skyhook control (TSKY) for a full-vehicle system is simulated too. It is controlled from its four corners according to the skyhook control scheme for the quarter-car model. This controller uses eight sensors (i.e., four acceleration and four deflection sensors). For simplicity, the MRDs are considered as ideal linear adjustable dampers, although a polynomial model has been proposed in our previous research [14].

The suspension system in this study has been simulated for two cases, in which a delay is introduced between the front and rear axles. In the first case, the periodic road inputs have three different values of frequency, that is, $ω =$ 8 rad/s, $ω =$ 69 rad/s, and $ω =$ 140 rad/s. The input terrain disturbances have been defined as sinusoids given by $z_{r f l} (t) = A_{f l} s i n (ω t),$ $z_{r f r} (t) = A_{f r} s i n (ω t)$ , $z_{r r l} (t) = A_{r l} s i n (ω (t + τ))$ , $z_{r r r} (t) = A_{r r} s i n (ω (t + τ))$ , where $ω$ is the terrain disturbance frequency and $A_{f l}$ , $A_{f r}$ , $A_{r l}$ and $A_{r r}$ are the terrain disturbance amplitudes at the front-left, front-right, rear-left, and rear-right wheels, respectively. The time delay is given by $τ = (a + b) / v$ and the speed at which the vehicle travels is given by $v$ . The vehicle is driving along a straight longitudinal trajectory at a constant speed $v = 24 m/s$ . Hence, longitudinal and lateral dynamics do not interfere with the vertical dynamics. Other simulation parameter values are selected as $A_{f l} = A_{r l} = 0.05 m$ , $A_{f r} = A_{r r} = 0.15 m$ . In the second case, the road disturbance is the “speed bump” used in [15], the profile and dimension of which are shown in Fig. 2. Only the left side of the vehicle meets a bump, so that both pitch and roll motion are produced at the same time. Moreover, the skyhook gains for heave, pitch, and roll are $C_{z} =$ 2000 N/m/s, $C_{θ} =$ 3000 N/m/s, and $C_{φ} =$ 3000 N/m/s, respectively.

Fig. 2‘Speed bump’ profile

Figs. 3-5 show the heave, pitch, and roll accelerations at low and high frequencies when the vehicle is given periodic road inputs. As can be seen, responses from low to high frequencies are considerably better for the new configuration compared with those from the traditional skyhook and passive cases.

Figure 6 shows the vehicle acceleration response for a speed bump. The vertical acceleration is shown in Fig. 6(a). The DSKY configuration clearly yields the best results. The reduction of peak acceleration is approximately 20 % more than that in the passive case and 10 % more than that in the traditional skyhook case. Both the DSKY and TSKY configurations attenuate the vertical motions faster than the passive configuration does. The pitch response is shown in Fig. 6(b). As can be seen, the DSKY configuration yields the smallest peak and the fastest time in eliminating the pitch motions. Fig. 6(c) shows the roll response. As can be seen, the DSKY configuration also yields the smallest peak. However, the passive configuration converges faster than do the DSKY and TSKY cases.

Fig. 3Vehicle acceleration responses for input frequency ω= 8 rad/s. ––––––– DSKY, - - - - - TSKY, -∙-∙-∙-∙- Passive

a) Heave acceleration

b) Pitch acceleration

c) Roll acceleration

Fig. 4Vehicle acceleration responses for input frequency ω= 69 rad/s. ––––––– DSKY, - - - - - TSKY, -∙-∙-∙-∙- Passive

a) Heave acceleration

b) Pitch acceleration

c) Roll acceleration

Fig. 5Vehicle acceleration responses for input frequency ω= 140 rad/s. ––––––– DSKY, - - - - - TSKY, -∙-∙-∙-∙- Passive

a) Heave acceleration

b) Pitch acceleration

c) Roll acceleration

Fig. 6Vehicle acceleration responses for ‘speed bump’ road exciting. ––––––– DSKY, - - - - - TSKY, -∙-∙-∙-∙- Passive

a) Heave acceleration

b) Pitch acceleration

c) Roll acceleration

6. Conclusions

A semi-active suspension using a novel methodology called DSKY was developed to control the attitude of the vehicle body and to suppress the vibration of the suspension. In accordance with the vehicle body motion (vertical, pitch, and roll), three skyhook controllers were designed. An input decoupling transformation converted the three orientation skyhook control forces into actual damping forces of four adjustable dampers. The simulations were conducted under two kinds of road inputs: periodic road inputs with three different frequencies and a “speed bump” input. The results showed that the proposed semi-active system provides better ride control in the heave, pitch, and roll motions of the vehicle body compare with the traditional semi-active and passive systems.

The time delay of system is an important issue that needs careful treatment to avoid poor performance or even possible instability of the control closed-loop system. For semi-active full-vehicle model systems, unavoidable time delays may appear in the controlled channel, in sensors, and in MRD where the delays are taken by the dampers to build up the required control force. Although this paper doesn’t consider time delay, it doesn’t affect the effectiveness of the proposed methodology and the simulation results. Usually, the stability and performance of the system due to a retarded control force have been studied only after the controller was designed [16]. The time delay of system will be considered before an experiment study in future research.

References

Karnopp D., Crosby M. J., Harwood R. A. Vibration control using semi-active force generators. J. of Engineering for Industry, 1974, p. 619-626.

Search CrossRef
Narges Maleki, Ali Khaki Sedigh, Batool Labibi Robust model reference adaptive control of active suspension system. Proceedings of the Control and Automation-MED, 14th Mediterranean Conference, 2006, p. 1-6.

Search CrossRef
Ardeshir Karami Mohammadi Vehicle active suspension control, using a variable structure model reference adaptive controller. Proceedings of the ACSE 05 Conference, CICC, Cairo, Egypt, 2005, p. 616.

Search CrossRef
Parthasarathy S., Srinivasa Y. G. Design of an active suspension system for a quarter-car road vehicle model using model reference control. J. Systems and Control Engineering, Vol. 220, p. 91-108.

Search CrossRef
Zifan Fang, Wenhui Shu, Daojia Du, et al. Semi-active suspension of a full-vehicle model based on double-loop control. Procedia Engineering, Vol. 16, 2011, p. 428-473.

Search CrossRef
Sedeh R. S., Ahmadian M. T., Abdollahpour R., et al. Application of car semi-active suspension systems to achieve desired performance on decreasing effect of road excitation on human health. Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Vol. 6, 2005, p. 375-382.

Search CrossRef
Yagiz N., Sakman L. E. Fuzzy logic control of a full vehicle without suspension gap degeneration. International Journal of Vehicle Design, Vol. 42, Issue 1-2, 2006.

Search CrossRef
Emmanuel D. Blanchard On the control aspects of semiactive suspensions for automobile applications. Master Dissertation, Blacksburg, Virginia, USA, Virginia Polytechnic Institute and State University, 2003.

Search CrossRef
Laoufi M., Eskandarian A. Fuzzy logic control for active suspension of a non-linear full-vehicle model. IEEE Intelligent Vehicles Symposium Proceedings, 2009, p. 677-684.

Search CrossRef
Ikenaga S., Lewis F. L., Campos J., et al. Active suspension control of ground vehicle based on a full-vehicle model. Proceedings of the American Control Conference, Vol. 6, 2000, p. 4019-4024.

Search CrossRef
Lakehal-Ayat M., Diop S., Fenaux E. Development of a full active suspension system. Proceedings of the IFAC World Congress, Barcelona, Spain, 2002, p. 2658.

Search CrossRef
Xubin Song Cost-effective skyhook control for semiactive vehicle suspension applications. The Open Mechanical Engineering Journal, Vol. 3, 2009, p. 17-25.

Search CrossRef
Chalasani R. M. Ride performance potential of active suspension system – part 2: comprehensive analysis based on a full-car model. Symposium on Simulation and Control of Ground Vehicles and Transportation Systems, Anaheim CA, ASME AMD-Vol. 80, DSC-Vol. 2, 1986, p. 205-234.

Search CrossRef
Yao Jialing, Zheng Jiaqiang, Gao Weijie, Zhang Zhenjun Sliding mode control of vehicle semi-active suspension with magnetorheological dampers having polynomial model. Journal of System Simulation, Vol. 21, Issue 8, 2009, p. 2400-2404, (in Chinese).

Search CrossRef
Dong Xiaomin Human simulated intelligent control of automobile magnetorheological semi-active suspension. Ph. D. Dissertation of Chongqing University, Chongqing, China, 2006.

Search CrossRef
Haiping Du, Nong Zhang H^∞ control of active vehicle suspensions with actuator time delay. Journal of Sound and Vibration, Vol. 301, 2007, p. 236-252.

Search CrossRef

About this article

Received

18 November 2013

Accepted

25 December 2013

Published

15 May 2014

Keywords

skyhook control

input decoupling

semi-active suspension

Acknowledgements

This work is funded by The Overseas Training Program Foundation for University Excellent Young and Middle-aged Teachers and Presidents by Jiangsu Province, and by China Postdoctoral Science Foundation funded project (2011M500935).

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.