A journal of IEEE and CAA , publishes high-quality papers in English on original theoretical/experimental research and development in all areas of automation
Volume 7 Issue 6
Oct.  2020

IEEE/CAA Journal of Automatica Sinica

  • JCR Impact Factor: 15.3, Top 1 (SCI Q1)
    CiteScore: 23.5, Top 2% (Q1)
    Google Scholar h5-index: 77, TOP 5
Tinghua Li, Qinghua Yang, Xiaowei Tu and Bin Ren, "An Improved Torque Sensorless Speed Control Method for Electric Assisted Bicycle With Consideration of Coordinate Conversion," IEEE/CAA J. Autom. Sinica, vol. 7, no. 6, pp. 1575-1584, Nov. 2020. doi: 10.1109/JAS.2020.1003360
Citation: Tinghua Li, Qinghua Yang, Xiaowei Tu and Bin Ren, "An Improved Torque Sensorless Speed Control Method for Electric Assisted Bicycle With Consideration of Coordinate Conversion," IEEE/CAA J. Autom. Sinica, vol. 7, no. 6, pp. 1575-1584, Nov. 2020. doi: 10.1109/JAS.2020.1003360

An Improved Torque Sensorless Speed Control Method for Electric Assisted Bicycle With Consideration of Coordinate Conversion

doi: 10.1109/JAS.2020.1003360
Funds:  This work was supported by the National Natural Science Foundation of China (51775325) and Hong Kong Scholars Program of China (XJ2013015)
More Information
  • In this paper, we propose an improved torque sensorless speed control method for electric assisted bicycle, this method considers the coordinate conversion. A low-pass filter is designed in disturbance observer to estimate and compensate the variable disturbance during cycling. A DC motor provides assisted power driving, the assistance method is based on the real-time wheel angular velocity and coordinate system transformation. The effect of observer is proved, and the proposed method guarantees stability under disturbances. It is also compared to the existing methods and their performances are illustrated through simulations. The proposed method improves the performance both in rapidity and stability.

     

  • WITH the improvement of modern life, people pay more and more attention to health and personal safety. Almost every country has been calling for green travel and the shared bicycles are very popular in all big cities around the world, like ofo bicycle and mobike [1]. So the development of intelligent bicycles must be an important trend in the future [2].

    Traditional bicycles are widely used in our daily life as an individual transportation tool, which has a merit to relieve the traffic pressure. Although bicycle travel is very environmental friendly, there are still some limitations. Since traditional bicycle is completely driven by cyclist, if the road is in a poor condition with various geometric complexity, the torque output by the cyclist may decrease as well as the speed of bicycle. Besides, it will be very laborious for people to ride on an uphill road.

    Once lacking of the power of human input, the bicycle will not maintain running. As far as this major drawback is concerned, man-powered bicycles are not as good as electric bicycles.

    Electric bicycle equipped with motor can solve this problem. The energy expenditure of riding an electric bicycle was 31% lower than the conventional bicycle on the uphill [3]. It has greatly extended the scope of application for traditional bicycle by providing assisted power with a motor [4]−[6]. Nowadays, electric bicycle has become one of the main vehicles in the world. Compared with the other vehicles powered by gasoline, electric bicycle can reduce the energy consumption, the air pollution and noise [7]. As we can see from Li et al. [8], the emissions of pollutants by electric bicycles per person per kilometer are several times smaller than those by motorcycles and cars, more or less equivalent to those by buses, and higher than those by bicycles.

    For electric bicycles, cyclist can control the output of motor by adjusting the power controller to adapt different road conditions or individual preferences. However, electric bicycle also has its limitations. Since the cyclist can control the motor freely, the speed of electric bicycles may be too fast to be safe. It not only affects the service life of the motor, but also increases pollution.

    Moreover, electric bicycles are heavy, when the motor runs out of electricity, it will be very hard for people to ride with his own human power.

    In order to protect the environment as well as to achieve the assisted riding, the electric bicycle equipped with torque sensor has been produced. The sensor detects the pedal torque and converts it into an electrical signal, which is the input for the motor controller [9]. However, the torque sensor may be easily damaged since it contacts with the cyclist’s feet directly and the damage costs are not negligible [10].

    Besides, the emerging development of connected and automated vehicles imposes a significant challenge on current vehicle control systems [11].

    To solve the problems mentioned above, the best way is to design a torque estimation method to measure the pedaling torque in real time and provide assistance [12], [13].

    Niki et al. [14] proposed a torque separation method based on high pass filter, they used high pass filter (HPF) to separate the human torque input from running friction. In 2015, Sankaranarayanan and Ravichandran [15] designed a torque sensorless controller for bicycles, but they did not take the disturbance force during cycling into account. The Fourier filtering was also applied to decouple the load torque and pedaling torque [16].

    In 2017, Fukushimal and Fujimoto [17] presented a method for pedaling torque estimation using the recursive least square algorithm with multiple forgetting factors and they considered the condition for travelling on upward as well. In 2018, Rallo et al. [18] measured the wheel speed with a magnetic wheel encoder and showed that the pedaling cadence is correlated with the speed oscillation.

    In addition, Bertucci et al. [19] analyzed the crank torque in road cycling on level and uphill, and the results show that when the pedaling cadences are the same, the crank torque profile differences are minimal between level ground and uphill road. In this way, we can see that the main factor affecting pedaling torque is cadences. The inertial measurement unit is also used to detect the angular velocity [20], [21].

    Nevertheless, none of them focused on the speed control. The purpose of electric assistance is to keep the speed within a certain range in all circumstances. And it is worthwhile to note that although the inertial measurement unit was used, none of the above papers considered the coordinate conversion of angular velocity.

    Similarly to the inertial navigation, in order to obtain the accurate positions, angular velocity based on the inertial coordinate system needs to be converted to another coordinate, which is consistent with the movement of the bicycle. For instance, Nilsson et al. [22] presented an open-source wireless foot-mounted inertial navigation module with an intuitive and significantly simplified dead reckoning interface, which can display the data after coordinate conversion.

    Disturbance observer is often used to eliminate the effects of interference on system control, considering that there are many kinds of disturbances during cycling [23].

    Therefore, in order to achieve a safe and convenient riding, it is very important to design an improved highly accurate speed control method for electric bicycles.

    This paper proposes an improved speed control method for the electric assisted bicycle. A second-order low-pass filter is designed in disturbance observer to estimate and compensate disturbances during cycling. The quaternion method is used to convert the coordinate and users can set a motor auxiliary ratio which is variable with the real-time speed of bicycle. The result finally comes out that the speed can be self-adaptively stabilized in a set safety range despite of the disturbances during cycling.

    The rest of the paper is organized as follows. Section II describes the dynamic model of the electric assisted bicycle system. In Section III, the quaternion method is used to re-coordinate the angular velocity which is measured by inertial measurement units (IMU). Besides, we simplify the calculation of pedaling torque. Section IV presents the design of low-pass filter (LPF) in disturbance observer to estimate and compensate the disturbances during riding. Section V is the introduction of motor assisted control. The stability of the proposed method is analyzed in Section VI. Simulation and results comparison are presented in Section VII. Finally, some conclusions are drawn and summarized in Section VIII.

    Table I lists the acronyms and symbols explanations used in this paper.

    Table  I.  List of Parameter
    Parameter Description
    ${E_t}$ Translational kinetic energy
    ${E_r}$ Rotational kinetic energy
    $M$ Mass of an electric assisted bicycle and a cyclist
    $V$ Velocity of bicycle
    $J$ Moment of inertia of wheel
    $\omega $ Angular velocity of wheels
    $T$ Total kinetic energy of bicycle
    $r$ Radius of rear wheel
    ${P_h}$ Gravitational potential energy of a bicycle
    $h$ Vertical height of cycling
    $l$ Length of a wheel movement
    $\varphi$ Slope angle of uphill
    $\theta$ Rotation angle of wheel
    ${F_{\rm{man}} }$ Pedaling force of cyclist
    ${F_{\rm{motor}} }$ Assisted power of motor
    ${T_{\rm{fric}} }$ Disturbances during cycling
    $B$ Friction coefficient
    ${T_{\rm{man}} }$ Pedaling torque
    ${J_m}$ Moment of inertia of crank shaft
    ${\omega _0}$ Pedaling angular velocity
    ${\omega _1}$ Angular velocity of the large gear
    ${\omega _2}$ Angular velocity of pinion
    ${\omega _3}$ Angular velocity of rear wheel
    $D$ Quaternion from inertial coordinate to geographic coordinate
    $N$ Quaternion from geographic coordinate to carrier coordinate
    $C$ Quaternion from inertial coordinate to carrier coordinate
    $\Omega $ Absolute angular velocity vectors of carrier coordinate
    ${\Omega _g}$ Absolute angular velocity vectors of geographic coordinate
    ${\omega _X}$, ${\omega _Y}$, ${\omega _Z}$ Absolute angular velocity vectors on pedal carrier coordinate
    $i$, $j$, $k$ Three axial unit vectors
    ${\omega _{Xg}}$, ${\omega _{Yg}}$, ${\omega _{Zg}}$ Absolute angular velocity vectors on geographic coordinate
    $ \lambda $ Longitude
    $\alpha $ Latitude
    ${\omega _t}$ Absolute angular velocity in carrier coordinate system
    ${\omega _{gt}}$ Angular velocity of carrier coordinate relative to the geographic coordinate
    $d$ Length of crank shaft
    $m$ Mass of crank shaft
    ${G_1}$ Model of bicycle
    $\tau $ Disturbance during cycling
    $\hat \tau $ Estimated disturbance
    ${G_{\tau {\omega _{\rm{3}}}}}$ Transfer function from $\tau $ to ${\omega _3}$
    ${G_{T{\omega _{\rm{3}}}}}$ Transfer function from $T$ to ${\omega _3}$
    ${Q_1}$ Second-order low-pass filter (LPF)
    ${i_1}$, ${i_2}$, ${i_3}$, ${i_4}$ Branch current
    ${V_1}$, ${V_2}$, ${V_3}$ Branch voltage
    ${A_u}$ Pass-band gain
    $Q$ Quality factor
    ${\omega _c}$ Cut-off angular frequency of LPF
    ${f_c}$ Cut-off frequency of LPF
    $c$ Auxiliary coefficient
    $k$ Transmission ratio between gears
    ${K_t}$ Motor constant
    $I$ Current of motor winding
    $U$ Average voltage of motor
    ${K_e}$ Motor constant
    $R$ Winding resistance
     | Show Table
    DownLoad: CSV

    Assuming that an electric assisted bicycle is running on an uphill road, the longitudinal dynamic model of the electric bicycle is shown in Fig. 1 and the forces during riding are indicated. We consider a rear-wheel-assisted electric bicycle, the motor is installed on the rear wheel, so as to collect the angular velocity and provide assisted power directly.

    Figure  1.  Electric assisted bicycle model

    The kinetic energy of a moving rigid body consists of translational kinetic energy and rotational kinetic energy.

    The expression of translational kinetic energy is

    $$ {E_t} = \frac{1}{2}M{V^2} $$ (1)

    where $ M $ is the mass of the electric assisted bicycle and a cyclist, $ V $ is the velocity of the bicycle.

    Considering there are two wheels on the bicycle and the expression of rotational kinetic energy is

    $$ {E_r} = 2 \left(\frac{1}{2}J{\omega ^2}\right) $$ (2)

    where $ J $ is the moment of inertia of the wheel and $ \omega $ is the angular velocity of the wheel.

    Therefore, the total kinetic energy of the electric assisted bicycle can be expressed as

    $$ T = {E_t} + {E_r} = \frac{1}{2}M{V^2} + 2 \left(\frac{1}{2}J{\omega ^2}\right) $$ (3)

    where $ T $ is the total kinetic energy of the bicycle.

    Considering $V = \omega \times r$, where $ r $ is the radius of rear wheel, here we suppose that the front wheel has the same radius as the rear wheel, T can be written as

    $$ T = \frac{1}{2}M{(\omega \times r)^2} + J{\omega ^2}. $$ (4)

    Due to the uphill road, the gravitational potential energy of the bicycle $ {P_h} $ can be expressed as

    $$ {P_h} = Mgh $$ (5)

    where $ h $ is the vertical height of cycling.

    Fig. 2 shows the uphill road riding, $ l $ is the length of the wheel movement, $ \varphi $ is the slope angle of the uphill, and $ \theta $ is the rotation angle of the wheel.

    Figure  2.  Uphill road riding

    With the equation of arc length calculation, we have

    $$ l = r \times \theta . $$ (6)

    The vertical height of cycling can be rewritten as

    $$ h = (r \times \theta ) \times \sin \varphi . $$ (7)

    The potential energy of the bicycle can be rewritten as

    $$ {P_h}{\rm{ = }}Mg(r \times \theta ) \times \sin \varphi . $$ (8)

    Based on the kinetic equation and the Lagrange function $ L = T - {P_h} $, we get the Lagrange function

    $$ L = \frac{1}{2}M{\omega ^2}{r^2} + J{\omega ^2} - Mgr\theta \sin \varphi . $$ (9)

    The Lagrange function satisfies $ \dfrac{d}{{dt}}\left(\dfrac{{\partial L}}{{\partial \theta }}\right) - \dfrac{{\partial L}}{{\partial \theta }} = Q $. Considering there are two major powers during cycling, which are the pedaling force of cyclist ${F_{\rm{man}}}$ and the assisted power of motor ${F_{\rm{motor}}}$, we have

    $$ Q = {T_{\rm{man}}} + {T_{\rm{motor}}} . $$ (10)

    At the same time, considering the disturbances during cycling ${T_{\rm{fric}}}$ and the impedance force between tires and ground, $ B $ is friction coefficient, the final Lagrange equation can be written as follows:

    $$ (M{r^2} + 2J)\ddot \theta + B\dot \theta + Mgr\sin \varphi + {T_{\rm{fric}}} = {T_{\rm{man}}} + {T_{\rm{motor}}} . $$ (11)

    Based on $ \ddot \theta = \dot \omega $ and $ \dot \theta = \omega $, the equation (11) can be rewritten as

    $$ (M{r^2} + 2J)\dot \omega + B\omega + Mgr\sin \varphi + {T_{\rm{fric}}} = {T_{\rm{man}}} + {T_{\rm{motor}}} . $$ (12)

    Then the transfer function can be expressed as

    $$ (M{r^2} + 2J)\dot \omega s + B\omega = {T_{\rm{man}}} + {T_{\rm{motor}}} - {T_{\rm{fric}}} - Mgr\sin \varphi . $$ (13)

    Assuming $ {J_{\rm{1}}} =M{r^2}{\rm{ + 2}}J $, the block diagram of a transfer function of the bicycle is shown in Fig. 3.

    Figure  3.  Block diagram of a transfer function

    To get the pedaling torque ${T_{\rm{man}}}$, we can simplify the pedaling process into a circular motion and use moment formula $ M = {J_m}\dfrac{{d\omega }}{{dt}} $ to replace it, where $ {J_m} $ is the moment of inertia and $ \omega $ is the angular velocity of the pedal rotation.

    Inertial measurement unit (IMU) is used to collect the pedal angular velocity, it is mainly composed of a three axes of gyroscope and an three directional accelerometer, which can measure the data of the three axes of the $ OXYZ $ coordinate system.

    In order to measure the angular velocity of the pedal and the rear wheel, one IMU is installed directly on the right pedal of the bicycle, and the other IMU is firmly connected with the rear wheel to detect the angular velocity of the electric assisted bicycle.

    Firstly, the coordinate must be mentioned. In this paper, there are three kinds of coordinate systems. The first one is the earth centered inertial coordinate system ($ i $ system). The origin of coordinate is located at the earth’s center, $ {O_i}{Z_i} $ is along the polar axis of the earth, $ {O_i}{X_i} $ and $ {O_i}{Y_i} $ are on the equatorial plane of the earth.

    The second one is the geographic coordinate system ($ t $ system), where the coordinate origin is also the barycenter of the carrier. $ {O_t}{Z_t} $ axis points to the zenith along the local reference ellipsoid’s normal line, $ {O_t}{X_t} $ axis is along the local latitude line to the east and the $ {O_t}{Y_t} $ axis points to the north along the local meridian.

    The third is the carrier coordinate system ($ b $ system), which is firmly connected with a carrier. In this paper, the inertial measurement unit is loaded on the pedal, so the pedal is the carrier. The origin of this coordinate is located at the barycenter of pedal. The $ {O_b}{X_b} $ axis points to the direction of bicycle heading on the plane of the carrier, the $ {O_b}{Y_b} $ axis is perpendicular to the $ {O_b}{X_b} $ axis and points to the right side. The $ {O_b}{Z_b} $ axis is perpendicular to the plane that consists of the $ {O_b}{X_b} $ axis and $ {O_b}{Y_b} $ axis.

    $ {\omega _0} $ is the pedaling angular velocity, $ {\omega _1} $ is the angular velocity of the large gear and $ {\omega _2} $ is the angular velocity of the pinion. The angular velocity $ \omega $ that we mentioned in Section II is the angular velocity $ {\omega _3} $ of the rear wheel, which can be detected by the IMU. However, IMU measures the projection of the absolute angular velocity vector on the axis of the carrier coordinate system. To be exact, they are not equal to the previous $ {\omega _0} $ and $ {\omega _3} $.

    In order to obtain the accurate angular velocity and acceleration to control the speed of the electric bicycle, it is necessary to get the real-time data of the carrier coordinate system (pedal) relative to the geographic coordinate system. In this paper, the quaternion method is used to convert the data measured in the carrier coordinate system to the geographic coordinate system [24].

    Matveev et al. [25] proposed a coordinate transformation method in 2009. Assume that the quaternion of the transformation from the inertial coordinate system to the geographic coordinate system is $ D $, the quaternion of the transformation from the geographic coordinate system to the carrier coordinate system is $ N $, and $ C $ is the quaternion that only needs to be transformed once from the inertial coordinate system to the carrier coordinate system. It can be composed of the two intermediate coordinate transformations. Based on the principle of quaternion, we get

    $$ C = D \circ N $$ (14)

    where “$ \circ $” is the product symbol of the quaternions.

    The purpose of the paper is to get the data of the coordinate system relative to the geographic coordinate system, so we need to calculate the quaternion $ N $.

    Left multiplying both sides of (14) by the conjugate quaternion $ \bar D $, we can get

    $$ \bar D \circ C = \bar D \circ D \circ N .$$ (15)

    Considering $ \bar D \circ D = 1 $ and reversing the equation, we get

    $$ N = \bar D \circ C . $$ (16)

    Take the derivative of both sides of (16)

    $$ \dot N = \dot {\bar D} \circ C + \bar D \circ \dot C . $$ (17)

    Using the dynamic equation of quaternion $ 2\dot N = N \circ \Omega $ and $ 2\dot {\bar N} = - \Omega \circ \bar N $, $ \dot C $ and $ \dot {\bar D} $ are written as follows

    $$ 2\dot C = C \circ \Omega \quad\;\; $$ (18)
    $$ 2\dot {\bar D} = - {\dot \Omega _g} \circ \bar D . $$ (19)

    Substituting (18) and (19) into (17), we get

    $$ \dot N = - \frac{1}{2}{\Omega _g} \circ \bar D \circ C + \bar D \circ \frac{1}{2}C \circ \Omega . $$ (20)

    Substituting (16) into (20), the final expression can be written as

    $$ 2\dot N = N \circ \Omega - {\Omega _g} \circ N $$ (21)

    where $ \Omega $ and $ {\Omega _g} $ are the absolute angular velocity vectors of the carrier coordinate system and geographic coordinate system. The three components of $ \Omega $ can be measured directly by the three axes gyroscopes in IMU.

    Assuming $ {\omega _X} $, $ {\omega _Y} $ and $ {\omega _Z} $ are the absolute angular velocity vectors on the three axes of the pedal carrier coordinate system, $ \Omega $ can be expressed as

    $$ \Omega = i{\omega _X} + j{\omega _Y}{\rm{ + }}k{\omega _Z} $$ (22)

    where $ i $, $ j $ and $ k $ are the three axial unit vectors, and the quaternion matrix for the corresponding $ \Omega $ can be described as

    $$ M({\omega _t}) = \left[ {\begin{array}{*{20}{c}} 0&{ - {\omega _X}}&{ - {\omega _Y}}&{ - {\omega _Z}}\\ {{\omega _X}}&0&{ - {\omega _Z}}&{{\omega _Y}}\\ {{\omega _Y}}&{{\omega _Z}}&0&{ - {\omega _X}}\\ {{\omega _Z}}&{ - {\omega _Y}}&{{\omega _X}}&0 \end{array}} \right]. $$ (23)

    Assuming that $ {\omega _{Xg}} $, $ {\omega _{Yg}} $ and $ {\omega _{Zg}} $ are the absolute angular velocity vectors on the three axes of the geographic coordinate system, $ {\Omega _g} $ can be expressed as

    $$ \qquad\qquad {\Omega _g} = i{\omega _{Xg}}{\rm{ + }}j{\omega _{Yg}} + k{\omega _{Zg}} $$ (24)
    $$ \qquad\qquad {\omega _{Xg}} = ({\omega _{ie}} + \dot \lambda )\cos \alpha $$ (25)
    $$ \qquad\qquad {\omega _{Yg}} = ({\omega _{ie}} + \dot \lambda )\sin \alpha $$ (26)
    $$ \qquad\qquad {\omega _{Zg}} = - \dot \alpha $$ (27)

    where $ {\omega _{ie}} $ is the rotational angular velocity of the earth and the value is $7.292 \times {10^{ - 5}}\;{\rm{rad/s}}$, $ \lambda $ is the longitude and $ \alpha $ is the latitude.

    Assuming $ N = {\lambda _0} + i{\lambda _1} + j{\lambda _2} + k{\lambda _3} $ and substituting $ \Omega $ and $ {\Omega _g} $ into (21), we can get

    $$ \begin{split} {\rm{2}}\left[ {\begin{array}{*{20}{c}} {{{\dot \lambda }_{\rm{0}}}}\\ {{{\dot \lambda }_{\rm{1}}}}\\ {{{\dot \lambda }_{\rm{2}}}}\\ {{{\dot \lambda }_{\rm{3}}}} \end{array}} \right]=\; & \left[ {\begin{array}{*{20}{c}} {{\lambda _{\rm{0}}}}&{-{\lambda _{\rm{1}}}}&{-{\lambda _{\rm{2}}}}&{-{\lambda _{\rm{3}}}}\\ {{\lambda _{\rm{1}}}}&{{\lambda _{\rm{0}}}}&{-{\lambda _{\rm{3}}}}&{{\lambda _{\rm{2}}}}\\ {{\lambda _{\rm{2}}}}&{{\lambda _{\rm{3}}}}&{{\lambda _{\rm{0}}}}&{-{\lambda _{\rm{1}}}}\\ {{\lambda _{\rm{3}}}}&{-{\lambda _{\rm{2}}}}&{{\lambda _{\rm{1}}}}&{{\lambda _{\rm{0}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {\rm{0}}\\ {{\omega _X}} \end{array}}\\ {{\omega _Y}} \end{array}}\\ {{\omega _Z}} \end{array}} \right] \\ & - \left[ {\begin{array}{*{20}{c}} 0&{ - {\omega _{Xg}}}&{ - {\omega _{Yg}}}&{ - {\omega _{Zg}}}\\ {{\omega _{Xg}}}&0&{ - {\omega _{Zg}}}&{{\omega _{Yg}}}\\ {{\omega _{Yg}}}&{{\omega _{Zg}}}&0&{ - {\omega _{Xg}}}\\ {{\omega _{Zg}}}&{ - {\omega _{Yg}}}&{{\omega _{Xg}}}&0 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\lambda _0}}\\ {{\lambda _1}}\\ {{\lambda _2}}\\ {{\lambda _3}} \end{array}} \right]. \end{split} $$ (28)

    Then, (28) can be simplified to

    $$ \left\{ \begin{aligned} & 2{{\dot \lambda }_0} = {\lambda _1}({\omega _{Xg}} - {\omega _X}) + {\lambda _2}({\omega _{Yg}} - {\omega _Y}) + {\lambda _3}({\omega _{Zg}} - {\omega _Z})\\ & 2{{\dot \lambda }_1} = {\lambda _0}({\omega _X} - {\omega _{Xg}}) + {\lambda _2}({\omega _Z} + {\omega _{Zg}}) - {\lambda _3}({\omega _Y} + {\omega _{Yg}})\\ & 2{{\dot \lambda }_2} = {\lambda _0}({\omega _Y} - {\omega _{Yg}}) - {\lambda _1}({\omega _Z} + {\omega _{Zg}}) + {\lambda _3}({\omega _X} + {\omega _{Xg}})\\ & 2{{\dot \lambda }_3} = {\lambda _0}({\omega _Z} - {\omega _{Zg}}) + {\lambda _1}({\omega _Y} + {\omega _{Yg}}) - {\lambda _2}({\omega _X} + {\omega _{Xg}}) .\end{aligned} \right. $$ (29)

    By solving the first order equation group (29), we can get the values of $ {\lambda _0} $,$ {\lambda _1} $,$ {\lambda _2} $,$ {\lambda _3} $, and the expression $ N = {\lambda _0} + i{\lambda _1} + j{\lambda _2} + k{\lambda _3} $.

    $ N $ is the quaternion that transforms from the geographic coordinate system to the carrier coordinate system, while our purpose is to convert the coordinate reversely. Therefore, after measuring the absolute angular velocity $ {\omega _t} $ in the carrier coordinate system, the angular velocity $ {\omega _{gt}} $ of the carrier coordinate system relative to the geographic coordinate system can be obtained by

    $$ {\omega _{gt}} = N \circ {\omega _t} \circ \bar N . $$ (30)

    For the convenience of calculation, we simplify the human pedal model in this paper and the pedaling torque formula can be described as

    $$ {T_{\rm{man}}} = {J_m}\frac{{d{\omega _0}}}{{dt}}\quad\quad\quad\quad\quad\quad\quad \quad $$ (31)
    $$ {\omega _0} = \left| {{\omega _{gt}}} \right| = \sqrt[2]{{{\omega _{gtx}}^2 + {\omega _{gty}}^2 + {\omega _{gt\textit{z}}}^2}} $$ (32)

    where $ {\omega _0} $ is the angular velocity of pedal, $ {J_m} $ is the moment of inertia of the crank shaft. $ d $ and $ m $ are the length and mass of crank shaft respectively.

    The chain transmission model of electric bicycle is shown in Fig. 4. With the principle of chain transmission, the rotational angular velocity of the rear wheel $ {\omega _3} $ can be calculated by the pedaling angular velocity $ {\omega _0} $. However, $ {\omega _3} $ is not proportional to $ {\omega _0} $ because of the disturbance factors such as inertia and friction during cycling. When the pedal stops rotating, $ {\omega _3} $ will still exist. That is why in this paper, an inertial measurement unit is loaded on the rear wheel to detect the angular velocity $ {\omega _3} $ in real time.

    Figure  4.  Chain transmission model of the electric bicycle

    As we mentioned in Section II, the electric bicycle will be disturbed by many kinds of frictions during riding. This paper mainly considers three kinds of frictions, namely the friction between tires and road, the wind resistance and the impedance force due to the gravitational potential energy.

    On the basis of kinematic model, Huang et al. [26] used a sensor fusion algorithm as an observer. In this paper, in order to compensate the disturbances of frictions during riding, a disturbance observer is used to estimate the interference, as shown in Fig. 5 [27]. Considering the variability of frictions, we set a special cut-off frequency in LPF so that the majority of disturbances can be observed and compensated.

    Figure  5.  Block diagram of the proposed method

    $ {G_1} $ is the model of the bicycle, we assume that the bicycle system is a first order model defined as $ {G_1}(s) = \dfrac{{\rm{1}}}{{{J_1}s + B}} $, $ {G_{\rm{2}}}(s) = {G^{{\rm{ -1}}}(s)} $, where $ {J_1} = M{r^2} + 2J $.

    $ \tau $ is the disturbance during cycling and $\tau = - {T_{\rm{fric}}}$, it can be observed and compensated by disturbance observer described in this paper. In this way, the transfer function from $ \tau $ to the angular velocity of the wheel $ {\omega _3} $ can be written as

    $$ {G_{\tau {\omega _{\rm{3}}}}} = \frac{{{G_1}(1 - {Q_1})}}{{1 + {G_1}{G_2}{Q_1} - {Q_1}}} . $$ (33)

    $ T $ is the input of the observer. With Mason formula, the transfer function from $ T $ to $ {\omega _3} $ can be expressed as

    $$ {G_{T{\omega _{\rm{3}}}}} = \frac{{{G_1}}}{{1 + {G_1}{G_2}{Q_1} - {Q_1}}} . $$ (34)

    $ \hat \tau $ is the estimated disturbance, it can be observed to compensate $ \tau $, the equation is as follow:

    $$ \hat \tau = [(t + \tau ){G_1}{{G}_2}{ - t]} \times {{{Q}}_1}. $$ (35)

    $ {Q_1} $ is a second-order low-pass filter (LPF), it is used in observer to filter out high-frequency [28], [29]. As shown in (36), when $ {Q_{\rm{1}}}(s) = 1 $ in the low-frequency stage, we get $\hat \tau = \tau$. The disturbance can be estimated in this way. Fig. 6 shows the circuit diagram of LPF.

    Figure  6.  The second-order low-pass filter (LPF)

    Firstly, with Kirchhoff’s law, the following can be obtained

    $$ \qquad\qquad \left\{ \begin{aligned} &{i_1} = \frac{{{V_1} - {V_2}}}{{{R_1}}}\\ & {i_2} = ({V_2} - {V_3}){C_1}s = ({V_2} - {V_p}){C_1}s\\ & {i_3} = \frac{{{V_2} - {V_p}}}{{{R_2}}}\\ & {i_4} = {V_p}{C_2}s\\ & {i_1} = {i_2} + {i_3}\\ & {i_3} = {i_4} \text{.} \end{aligned} \right. $$ (36)
    $$ \qquad\qquad \left\{ \begin{aligned} &\frac{{{V_1} - {V_2}}}{{{R_1}}} = ({V_2} - {V_3}){C_1}s + \frac{{{V_2} - {V_3}}}{{{R_2}}}\\ & \frac{{{V_{\rm{2}}} - {V_3}}}{{{R_2}}} = {V_3}{C_2}s .\end{aligned} \right. $$ (37)

    Therefore, the transfer function of LPF can be expressed as

    $$ \qquad\qquad \begin{array}{l} {Q_1}(s) = \dfrac{{{A_u} \cdot {\omega _c}^2}}{{{s^2} + \frac{{{\omega _c}}}{Q}s + {\omega _c}^2}} \end{array} $$ (38)
    $$ \qquad\qquad \left\{ \begin{aligned} &{\omega _c} = \sqrt {\frac{1}{{{R_1}{R_2}{C_1}{C_2}}}} = 2\pi {f_c}\\ & {A_u} = 1\\ & Q = \frac{{\sqrt {{R_1}{R_2}{C_1}{C_2}} }}{{({R_1} + {R_2}){C_2}}} \end{aligned} \right. $$ (39)

    where $ {A_u} $ is the pass-band gain, $ Q $ is quality factor, $ {\omega _c} $ is the cut-off angular frequency of LPF and $ {f_c} $ is the cut-off frequency.

    The filtering effect on a certain frequency section can be achieved by selecting the different cut-off angular frequencies. In the low-frequency stage, $ {Q_{\rm{1}}}(s) = 1 $, $ {G_{T{w_{\rm{3}}}}} = {G_{\rm{2}}}^{{\rm{ - 1}}} $ and $ {G_{\tau {\omega _{\rm{3}}}}} = 0 $, the effect of disturbance can be compensated by $ {G_{\rm{2}}} $.

    Similarly, in the high-frequency stage, $ {Q_{\rm{1}}}(s) = 0 $, the disturbance which is higher than $ {f_c} $ is effectively suppressed.

    In general, people often pedal 65−80 times in a minute when they ride bicycles in a normal frequency. The corresponding frequency is from 1.083 Hz to 1.33 Hz. But the frequency can not be cut off immediately in the low-pass filter, because of the transition-band. Therefore, the cut-off frequency of LPF is set to be 4 Hz in this paper.

    Taking the characteristics of the circuit into account, the resistance value cannot be too large, and the $ Q $ value is preferably less than 3. We set $ Q = 0.7 $, then the value of $ {R_1} $, $ {R_2} $, $ {C_1} $, $ {C_2} $ can be calculated.

    The expression of LPF can be described as follows:

    $$ {Q_{\rm{1}}}{\rm{(s) = }}\frac{{632}}{{{s^2} + 36s + 632}} . $$ (40)

    Fig. 7 shows the amplitude-frequency characteristic of designed LPF. We can see clearly that the curve shows a declining trend after 4 Hz. As indicated in the label, when it is at the cut-off frequency, the gain value is −2.93 dB, which is in accordance with the design of the low-pass filter.

    Figure  7.  The amplitude-frequency characteristic of LPF

    In order to test and verify the accuracy of disturbance observer $ \tau $, we consider the disturbance as sine wave, the expression is $ \tau = {\rm{2}}0\sin t $.

    The simulation results are shown in Fig. 8. The dotted line represents the simulated disturbance $ \tau $ and the solid line is the estimated disturbance $ \hat \tau $ through observation. As we can see, there are small phase and amplitude differences between the two waveforms.

    Figure  8.  Disturbance and observation

    Fig. 9 captures and amplifies the waveforms from 0.6 s to 2 s. The two dashed lines indicate the amplitude values of the waveform vertexes respectively. It is clear that the difference in amplitude is 0.4, which proves that the observer can estimate the majority of disturbance.

    Figure  9.  Difference between disturbance and observation

    To achieve the purpose of the assisted riding, a permanent magnet brushless DC motor is used in this paper and installed in the rear wheel.

    When the riding speed decreases and riding becomes more laborious, the motor can provide assistance based on the real-time angular velocity of wheel. As shown in Fig. 5, $ {\omega _{\rm{3}}} $ is fed back to the motor-assist module, and the auxiliary equation can be expressed as

    $$ {T_{\rm{motor}}} = c{T_{\rm{man}}} $$ (41)

    where $ c $ is the auxiliary coefficient determined in real-time by the wheel angular velocity.

    Japanese industry proposed the PAP method, which has been adopted by regulation. The PAP method can be formulated by the linear piecewise equation [30]. Considering the motor assistance, we rewrite the auxiliary coefficient, and the specific relationship can be expressed as

    $$ c = \left\{ {\begin{array}{*{20}{l}} 1, & \qquad\qquad {\omega < {\omega _{\min }}}\\{(0,1),} & \qquad\qquad{{\omega _{\min }} < \omega < {\omega _{\max }}} \\{ - 1,}& \qquad\qquad{\omega > {\omega _{\max }}} \end{array}} \right. $$ (42)

    where $ {\omega _{\min }} $ and $ {\omega _{\max }} $ are the pre-defined angular velocity limits. Once the wheel angular velocity $ {\omega _{\rm{3}}} $ is lower than the minimum limit $ {\omega _{\min }} $, the motor will provide assistance and the auxiliary coefficient $ c = 1 $. The motor and cyclist will provide power together with a proportion of 1:1.

    When $ {\omega _{\rm{3}}} $ exceeds the maximum limit $ {\omega _{\max }} $, the motor will be reversed and $ c = {\rm{ - }}1 $. In this way, the electric bicycle will not be overspeed and become too dangerous.

    As the permanent magnet brushless DC motor is used in this paper, it is equipped with position sensors. According to the detected difference of the rotor position, each phase coil will be electrified in turns, making the direction change of magnetic field, which is produced by stators. Therefore, we can use logic circuits to control the conduction sequence, so as to achieve motor inversions.

    In the same way, if $ {\omega _{\rm{3}}} $ is in a range of $ {\omega _{\min }} \leq \omega \leq {\omega _{\max }} $, the auxiliary coefficient will vary between 0 and 1, so as to provide assistance intelligently.

    As shown in Fig. 10, the predefined wheel angular velocity limits are $ {\omega _{\min }}{\rm{ = 12}} $ rad/s and $ {\omega _{\max }} = 16 $ rad/s respectively. The auxiliary coefficient values at different wheel angular velocities are shown in the figure.

    Figure  10.  Auxiliary coefficient

    Because of the gear transmission between motor and the rear wheel, the formula of speed transmission can be expressed as

    $$ {\omega _{\rm{wheel}}} = k{\omega _{\rm{motor}}} $$ (43)

    where ${\omega _{\rm{wheel}}}$ is the angular velocity of the rear wheel, $ k $ is the transmission ratio between gears and ${\omega _{\rm{motor}}}$ is the angular velocity of motor.

    In this paper, $ {\omega _{\rm{3}}} $ is ${\omega _{\rm{wheel}}}$, and the value of $ k $ is related to the radii of the gears.

    For the permanent magnet brushless DC motor, the electromagnetic torque of the motor can be expressed as

    $$ {T_{\rm{motor}}} = {K_t} \times I $$ (44)

    where $ {K_t} $ is the motor constant, and $ I $ is the current of the motor winding.

    The motor winding current $ I $ can be calculated as

    $$ I=\frac{{U - {K_e} \times {\omega _{\rm{motor}}}}}{R} $$ (45)

    where $ U $ is the average voltage of motor, $ {K_e} $ is the motor constant, ${\omega _{\rm{motor}}}$ is the angular velocity of motor and $ R $ is the winding resistance.

    So, ${T_{\rm{motor}}}$ can be rewritten as

    $$ {T_{\rm{motor}}} = {K_t} \times \frac{{U - {K_e}{\omega _{\rm{motor}}}}}{R} . $$ (46)

    Substituting ${T_{\rm{motor}}} = c{T_{\rm{man}}}$ into (46), we have

    $$ {\omega _{\rm{motor}}} = \frac{{U{K_t} - cR{T_{\rm{man}}}}}{{{K_t}{K_e}}} . $$ (47)

    Considering (43), the angular velocity of the bicycle wheel can be rewritten as

    $$ {\omega _{\rm{3}}} ={\omega _{\rm{wheel}}} = \frac{{k \times (U{K_t} - cR{T_{\rm{man}}})}}{{{K_t}{K_e}}}. $$ (48)

    In this way, the speed of the electric bicycle can be controlled and dynamically stabilized in a certain safe range.

    Based on the controller (41) and the coefficient function $ c(\omega_3) $, we have

    $$ \qquad T_{\rm{man}} = J_m \ddot{\omega}_0 $$ (49)
    $$ \qquad T_{\rm{motor}} = c(\omega_3)T_{\rm{man}} = c(\omega_3)J_m\dot{\omega}_0 . $$ (50)

    Combining (11), the differential equation of the system can be written as

    $$ (M{r^2} + 2J){{\dot \omega }_3} + B\omega_3+Mgr\sin\varphi +\tau -\hat{\tau } = J_m\dot{\omega }_0+\\c(\omega_3)J_m\dot{\omega}_0.$$ (51)

    The system dynamics can be written as

    $$ x = \omega_3-\omega^{*} $$ (52)
    $$ \dot{x} = \dot{\omega }_3 = \frac{J_m \dot{\omega}_0+c(\omega _3)J_m\dot{\omega }_0+\hat{\tau} -\tau -Mgr\sin\varphi -B\omega _3}{Mr^2+2J} $$ (53)

    where $ \omega^{*} $ is the pre-set expected angular velocity, and we have

    $$ J_m \dot{\omega}_0+c(\omega ^*)J_m\dot{\omega}_0+\hat{\tau} -\tau -Mgr\sin\varphi -B\omega ^* = 0. $$ (54)

    Consider the Lyapunov function $ V(x) $ as

    $$ V(x) = \frac{1}{2}x^2 . $$ (55)

    Its time derivative is given by

    $$ \begin{split} \dot{V}(x) =\; & x\dot{x}\\ =\; & x\frac{J_m \dot{\omega}_0+c(x+\omega ^*)J_m\dot{\omega }_0+\hat{\tau} -\tau -Mgr\sin\varphi }{Mr^2+2J}-x\frac{B(x+\omega ^*)}{Mr^2+2J}\\ =\; & \frac{J_m \dot{\omega }_0\left[ c(x+\omega ^*)-c(\omega ^{*}) \right ]x-Bx^2}{Mr^2+2J}\\ =\; & \frac{J_m \dot{\omega }_0Ax-Bx^2}{Mr^2+2J} \\[-17pt]\end{split} $$ (56)

    where $ A = \left[ c(x+\omega ^*)-c(\omega ^{*}) \right ] $. Since $ c $ is a function of $ \omega_3 $, $ \omega ^* $ is set between $\omega_{\rm{min}}$ and $\omega_{\rm{max}}$. Then $ x\neq 0 $, $ Ax <0 $. We also have $ Ax = 0 $ when $ x = 0 $.

    Due to the fact that $J_m > 0,\; \dot{\omega}_0>0,\;B>0$, then $ \dot{V}(x) $ is negative-definte and $ \dot{V}(x) = 0 $ only when $ x = \omega_3-\omega^* = 0 $. Therefore, the system is asymptotically stable at $ \omega^* $.

    To show the validity of the proposed method, several simulations are implemented. For the convenience of this research, the signals for the input angular velocities of the three directions $ XYZ $ are all assumed to be sinusoid.

    The expressions are depicted as $ {\omega _x} = 10\sin (2\pi t) + 10 $, $ {\omega _y} = 12\sin (2\pi t) + 12 $ and ${\omega _\textit{z}} = 14\sin (2\pi t) + 14$, respectively. The waveform is shown in Fig. 11.

    Figure  11.  Angular velocity before transformation

    Besides, the dynamic parameters used in simulation are shown in Table II.

    Table  II.  Simulation Parameter Values
    Parameter Value
    $M(man + cyclist)$ $100\;{\rm{kg}}$
    $J$ $0.15\;{\rm{N \cdot m \cdot {s^2}/rad}}$
    $r$ $0.3\;{\rm{m}}$
    $B$ $0.2$
    ${J_m}$ $0.12\;{\rm{kg \cdot {m^2}}}$
     | Show Table
    DownLoad: CSV

    Fig. 12 shows the angular velocity and three components curve diagram after coordinate transformation. We can see clearly that the value of angular velocity has changed from the original carrier coordinate system to the geographic coordinate system, so as to provide the subsequent calculation of the pedaling torque.

    Figure  12.  Angular velocity after transformation

    Take a set of data of the initial angular velocity as an example, $ {\omega _x} = 10 $, $ {\omega _y} = 12 $ and ${\omega _\textit{z}} = 14$. They make up an angular velocity matrix [0 10 12 14]T and the transformed matrix we get is [0 5.836 9.814 17.68]T. Thus, the quaternion from the geographic coordinate system to the carrier coordinate system can be calculated and expressed as N = [–0.9265 –0.216 –0.08526 –0.3018]T.

    As we have got the value of transformed angular velocity, the pedaling torque can be obtained by the torque formula. The result is shown in Fig. 13.

    Figure  13.  Pedaling torque

    With ${T_{\rm{man}}}$ and the angular velocity of the wheel, the motor module can provide assistance in real time.

    Article 58 of China’s road traffic safety law stipulates that the maximum velocity of an electric bicycle riding on a non-motorized lane must not exceed 15 km/h.

    In 2013, Zhao et al. [31] studied that the average velocity of a bicycle is 15 km/h on an isolated lanes between motor vehicles and non-motor vehicles with better riding conditions in China.

    The electric bicycle safety specification stipulates that the maximum velocity of an electric bicycle shall not exceed 25 km/h, which is issued in China in May 2018. Before this specification, the Chinese national standard stipulated that the maximum velocity of electric bicycles was 20 km/h.

    Considering that the electric assisted bicycle discussed in this paper is not a fully electric powered bicycle, we believe that velocity should be controlled according to the standard of the bicycle.

    Therefore, the safe velocity of ordinary people riding a bicycle is between 12 km/h and 16 km/h. Assuming the radius of wheel $ r $ is 0.3 m and the formula $v = \omega \times r$, we can get the range of angular velocity is 11.1−14.8 rad/s.

    Accordingly, the minimum limited angular velocity for the wheel can be set to 10 rad/s and the maximum limit is 13 rad/s. If the wheel angular velocity is lower than 10 rad/s, the motor and manpower will provide power together by the ratio of 1:1. If the wheel angular velocity is higher than 13 rad/s, motor will be reversed, so as to avoid danger.

    The angular velocity variation curves of the rear wheel are shown in Figs. 14 and 15.

    Figure  14.  Angular velocity of the wheel in previous methods
    Figure  15.  Angular velocity of wheel in previous and new proposed methods

    Three kinds of methods are simulated in our system, two of which are the previous method. A low pass filter is used in method 1. In method 2, low pass filter and band pass filter are used in combination to estimate disturbance and human pedaling torque.

    In Fig. 14, the dotted line indicates the result of previous method 1, and the solid line represents the previous method 2. After setting the same safety angular velocity range for both of methods, the rapidity of method 2 is faster but fluctuation is more turbulent.

    Fig. 15 shows the contrast curve between the results of this paper and the previous method 1. It can be seen that they can both well control the angular velocity but the rapidity and stability are not the same. Obviously, the method proposed in this paper is much faster and more stable, it takes reaction in the first 5 seconds.

    As the cyclist begins to step on, the angular velocity value is the largest at the beginning in Fig. 15. When the value exceeds 13 rad/s, motor will turn to the reverse direction to decrease the value and the curve will show a downward trend. Later, when the value is between 10 rad/s and 13 rad/s, the curve rises up due to the motor’s assistance.

    With the control and assistance, we can see clearly that the wheel angular velocity can still be effectively stabilized in a safe range, even though there are disturbances such as frictions and slopes during riding, so the same as the velocity of an electric bicycle.

    The contribution of this paper lies in the presenting of an improved speed control method for the electric assisted bicycle without torque sensor installed.

    On the basis of dynamic model system, we combine the disturbance observer, coordinate transformation and motor assistance to realize the speed control.

    Based on the simulation results above, the proposed method can well estimate and compensate the majority of disturbances during cycling, and improve the measurement and control accuracy by means of coordinate system conversion. Finally, the bicycle speed can be intelligently stabilized within a safe range with the motor assistance.

    Compared with the previous methods, this new proposed method has a higher rapidity and stability than that of the previous methods. Meanwhile, the stability of the method is also analyzed. It further validates that the speed control method can not only assist the cycling under different road conditions, but also improve the accuracy.

  • [1]
    Z. B. Li, W. Wang, C. Yang, and H. Y. Ding, “Bicycle mode share in China: A city-level analysis of long term trends,” Transportation, vol. 44, no. 4, pp. 773–788, Jul. 2017. doi: 10.1007/s11116-016-9676-8
    [2]
    E. Salmeron-Manzano and F. Manzano-Agugliaro, “The electric bicycle: Worldwide research trends,” Energies, vol. 11, pp. 1894, Jul. 2018.
    [3]
    B. C. Langford, C. R. Cherry, D. R. Bassett Jr., E. C. Fitzhugh, and N. Dhakal, “Comparing physical activity of pedal-assist electric bikes with walking and conventional bicycles,” J. Transp. Health, vol. 6, pp. 463–473, Sep. 2017. doi: 10.1016/j.jth.2017.06.002
    [4]
    X. Fan and M. Tomizuka, “Robust disturbance observer design for a power-assist electric bicycle,” in Proc. American Control Conf., Baltimore, USA, 2010, pp. 1166−1171.
    [5]
    Y. Li, J. F. Hu, F. B. Chen, S. P. Liu, Z. T. Yan, and Z. Y. He, “A new-variable-coil-structure-based IPT system with load-independent constant output current or voltage for charging electric bicycles,” IEEE Trans. Power Electron., vol. 33, no. 10, pp. 8226–8230, Oct. 2018. doi: 10.1109/TPEL.2018.2812716
    [6]
    M. Corno, D. Berretta, and S. M. Savaresi, “An IMU-driven rider-on-saddle detection system for electric-power-assisted bicycles,” IEEE Trans. Intell. Trans. Syst., vol. 17, no. 11, pp. 3184–3193, Nov. 2016. doi: 10.1109/TITS.2016.2544959
    [7]
    A. Wolf and S. Seebauer, “Technology adoption of electric bicycles: A survey among early adopters,” Trans. Res. Part A:Policy Pract., vol. 69, pp. 196–211, Nov. 2014.
    [8]
    T. Z. Li, F. Qian, and C. Su, “Energy consumption and emission of pollutants from electric bicycles,” Appl. Mech. Mater., vol. 505−506, pp. 327–333, Jan. 2014.
    [9]
    P. Spagnol, M. Corno, and S. M. Savaresi, “Pedaling torque reconstruction for half pedaling sensor,” in Proc. European Control Conf., Zurich, Switzerland, 2013, pp. 275−280.
    [10]
    D. S. Cheon and K. H. Nam, “Pedaling torque sensor-less power assist control of an electric bike via model-based impedance control,” Int. J. Autom. Technol., vol. 18, no. 2, pp. 327–333, Jan. 2017. doi: 10.1007/s12239-017-0033-5
    [11]
    F.-Y. Wang, N. N. Zheng, D. P. Cao, C. M. Martinez, L. Li, and T. Liu, “Parallel driving in CPSS: A unified approach for transport automation and vehicle intelligence,” IEEE/CAA J. Autom. Sinica, vol. 4, no. 4, pp. 577–587, Sep. 2017. doi: 10.1109/JAS.2017.7510598
    [12]
    H. Y. Guo, D. P. Cao, H. Chen, C. Lv, H. J. Wang, and S. Q. Yang, “Vehicle dynamic state estimation: State of the art schemes and perspectives,” IEEE/CAA J. Autom. Sinica, vol. 5, no. 2, pp. 418–431, Mar. 2018. doi: 10.1109/JAS.2017.7510811
    [13]
    L. Chen, X. M. Hu, W. Tian, H. Wang, D. P. Cao, and F. Y. Wang, “Parallel planning: A new motion planning framework for autonomous driving,” IEEE/CAA J. Autom. Sinica, vol. 6, no. 1, pp. 236–246, Jan. 2019. doi: 10.1109/JAS.2018.7511186
    [14]
    H. Niki, J. Miyata, and T. Murakami, “Power assist control of electric bicycle taking environment and rider’s condition into account,” in Proc. SPIE 6052, Optomechatronic Systems Control, Sapporo, Japan, 2005, pp. 60520H.
    [15]
    V. Sankaranarayanan and S. Ravichandran, “Torque sensorless control of a human-electric hybrid bicycle,” in Proc. Int. Conf. Industrial Instrumentation and Control, Pune, India, 2015, pp. 806−810.
    [16]
    T. Kurosawa, Y. Fujimoto, and T. Tokumaru, “Estimation of pedaling torque for electric power assisted bicycles,” in Proc. 40th Annu. Conf. of the IEEE Industrial Electronics Society, Dallas, USA, 2014, pp. 2756−2761.
    [17]
    N. Fukushima and Y. Fujimoto, “Estimation of pedaling torque for electric power-assisted bicycle on slope environment,” in Proc. IEEE Int. Conf. Advanced Intelligent Mechatronics, Munich, Germany, 2017, pp. 1682−1687.
    [18]
    G. Rallo, S. Formentin, M. Corno, and S. M. Savaresi, “Real-time cycling cadence estimation via wheel speed measurement,” Int. J. Adapt. Control Signal Process., vol. 32, no. 7, pp. 1052–1066, Jul. 2018. doi: 10.1002/acs.2885
    [19]
    W. Bertucci, F. Grappe, A. Girard, A. Betik, and J. D. Rouillon, “Effects on the crank torque profile when changing pedalling cadence in level ground and uphill road cycling,” J. Biomech., vol. 38, no. 5, pp. 1003–1010, May 2005. doi: 10.1016/j.jbiomech.2004.05.037
    [20]
    C. F. Huang, B. H. Dai, and T. J. Yeh, “Observer-based sensor fusion for power-assist electric bicycles,” in Proc. American Control Conf., Boston, USA, 2016, pp. 5486−5491.
    [21]
    D. M. Bevly, J. Ryu, and J. C. Gerdes, “Integrating INS sensors with GPS measurements for continuous estimation of vehicle sideslip, roll, and tire cornering stiffness,” IEEE Trans. Intell. Trans. Syst., vol. 7, no. 4, pp. 483–493, Dec. 2006. doi: 10.1109/TITS.2006.883110
    [22]
    J. O. Nilsson, A. K. Gupta, and P. Händel, “Foot-mounted inertial navigation made easy,” in Proc. Int. Conf. Indoor Positioning and Indoor Navigation, Busan, South Korea, 2014, pp. 24−29.
    [23]
    S. Y. Yu, J. Wang, Y. Wang, and H. Chen, “Disturbance observer based control for four wheel steering vehicles with model reference,” IEEE/CAA J. Autom. Sinica, vol. 5, no. 6, pp. 1121–1127, Nov. 2018. doi: 10.1109/JAS.2016.7510220
    [24]
    H. Y. Zhang, Z. H. Wang, L. S. Zhou, Q. N. Xue, L. Ma, and Y. F. Niu, “Explicit symplectic geometric algorithms for quaternion kinematical differential equation,” IEEE/CAA J. Autom. Sinica, vol. 5, no. 2, pp. 479–488, Mar. 2018. doi: 10.1109/JAS.2017.7510829
    [25]
    V. Matveev and V. Raspopov, Fundamentals of Designing Strapdown Inertial Navigation Systems. St. Petersburg, 2009, pp. 140.
    [26]
    C. F. Huang, B. H. Dai, and T. J. Yeh, “Determination of motor torque for power-assist electric bicycles using observer-based sensor fusion,” J. Dyn. Syst.,Meas. Control, vol. 14, no. 7, pp. 071019, Jul. 2018.
    [27]
    Z. L. Zhao, B. Z. Guo, and Z. P. Jiang, “A new extended state observer for output tracking of nonlinear MIMO systems,” in Proc. 56th IEEE Annu. Conf. Decision and Control, Melbourne, Australia, 2017, pp. 6732−6737.
    [28]
    F. F. Zhao, C. Q. Chen, W. He, and S. S. Ge, “A filtering approach based on MMAE for a SINS/CNS integrated navigation system,” IEEE/CAA J. Autom. Sinica, vol. 5, no. 6, pp. 1113–1120, Nov. 2018. doi: 10.1109/JAS.2017.7510445
    [29]
    M. Chen, “Robust tracking control for self-balancing mobile robots using disturbance observer,” IEEE/CAA J. Autom. Sinica, vol. 4, no. 3, pp. 458–465, Jul. 2017. doi: 10.1109/JAS.2017.7510544
    [30]
    C. L. Chen and M. C. Tsai, “Torque control of power-assisted bike via planetary gear mechanism,” in Proc. 20th Int. Conf. Electrical Machines and Systems, Sydney, Australia, 2017, pp. 1−5.
    [31]
    Y. Zhao and X. D. Pan, “Experimental research on the bicycle speed and turning radius on urban road intersection,” J. Transport Science and Engineering, vol. 29, no. 1, pp. 81–85, Mar. 2013.
  • Related Articles

    [1]Changjun He, Hua Geng, Kaushik Rajashekara, Ambrish Chandra. Analysis and Control of Frequency Stability in Low-Inertia Power Systems: A Review[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(12): 2363-2383. doi: 10.1109/JAS.2024.125013
    [2]Zheng Chen, Shizhao Zhou, Chong Shen, Litong Lyu, Junhui Zhang, Bin Yao. Observer-Based Adaptive Robust Precision Motion Control of a Multi-Joint Hydraulic Manipulator[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(5): 1213-1226. doi: 10.1109/JAS.2024.124209
    [3]Dan Zhang, Jiabin Hu, Jun Cheng, Zheng-Guang Wu, Huaicheng Yan. A Novel Disturbance Observer Based Fixed-Time Sliding Mode Control for Robotic Manipulators With Global Fast Convergence[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(3): 661-672. doi: 10.1109/JAS.2023.123948
    [4]Zhongcai Zhang, Guangren Duan. Stabilization Controller of An Extended Chained Nonholonomic System With Disturbance:  An FAS Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(5): 1262-1273. doi: 10.1109/JAS.2023.124098
    [5]Haoxiang Ma, Mou Chen, Qingxian Wu. Disturbance Observer-Based Safe Tracking Control for Unmanned Helicopters With Partial State Constraints and Disturbances[J]. IEEE/CAA Journal of Automatica Sinica, 2023, 10(11): 2056-2069. doi: 10.1109/JAS.2022.105938
    [6]Lu Zhang, Jun Yang, Shihua Li. A Model-Based Unmatched Disturbance Rejection Control Approach for Speed Regulation of a Converter-Driven DC Motor Using Output-Feedback[J]. IEEE/CAA Journal of Automatica Sinica, 2022, 9(2): 365-376. doi: 10.1109/JAS.2021.1004213
    [7]Guowei Dong, Liang Cao, Deyin Yao, Hongyi Li, Renquan Lu. Adaptive Attitude Control for Multi-MUAV Systems With Output Dead-Zone and Actuator Fault[J]. IEEE/CAA Journal of Automatica Sinica, 2021, 8(9): 1567-1575. doi: 10.1109/JAS.2020.1003605
    [8]Zhijia Zhao, Zhijie Liu. Finite-Time Convergence Disturbance Rejection Control for a Flexible Timoshenko Manipulator[J]. IEEE/CAA Journal of Automatica Sinica, 2021, 8(1): 157-168. doi: 10.1109/JAS.2020.1003378
    [9]Tongfu He, Zhong Wu. Iterative Learning Disturbance Observer Based Attitude Stabilization of Flexible Spacecraft Subject to Complex Disturbances and Measurement Noises[J]. IEEE/CAA Journal of Automatica Sinica, 2021, 8(9): 1576-1587. doi: 10.1109/JAS.2021.1003958
    [10]Dailin Zhang, Zining Wang, Masayoshi Tomizuka. A Variable-Parameter-Model-Based Feedforward Compensation Method for Tracking Control[J]. IEEE/CAA Journal of Automatica Sinica, 2020, 7(3): 693-701. doi: 10.1109/JAS.2020.1003135
    [11]Zhenhua Zhao, Jun Yang, Shihua Li, Wen-Hua Chen. Composite Nonlinear Bilateral Control for Teleoperation Systems With External Disturbances[J]. IEEE/CAA Journal of Automatica Sinica, 2019, 6(5): 1220-1229. doi: 10.1109/JAS.2018.7511273
    [12]Aishwarya Apte, Ujjwala Thakar, Vrunda Joshi. Disturbance Observer Based Speed Control of PMSM Using Fractional Order PI Controller[J]. IEEE/CAA Journal of Automatica Sinica, 2019, 6(1): 316-326. doi: 10.1109/JAS.2019.1911354
    [13]Jiajun Wang. Speed-assigned Position Tracking Control of SRM With Adaptive Backstepping Control[J]. IEEE/CAA Journal of Automatica Sinica, 2018, 5(6): 1128-1135. doi: 10.1109/JAS.2016.7510019
    [14]Mou Chen. Robust Tracking Control for Self-balancing Mobile Robots Using Disturbance Observer[J]. IEEE/CAA Journal of Automatica Sinica, 2017, 4(3): 458-465. doi: 10.1109/JAS.2017.7510544
    [15]Songsong Cheng, Shengguo Wang, Yiheng Wei, Qing Liang, Yong Wang. Study on Four Disturbance Observers for FO-LTI Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(4): 442-450.
    [16]Lu Wang, Jianbo Su. Trajectory Tracking of Vertical Take-off and Landing Unmanned Aerial Vehicles Based on Disturbance Rejection Control[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 65-73.
    [17]Pei Jia, Fei Hao, Hao Yu. Function Observer Based Event-triggered Control for Linear Systems with Guaranteed L-Gain[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(4): 394-402.
    [18]Chaoxu Mu, Qun Zong, Bailing Tian, Wei Xu. Continuous Sliding Mode Controller with Disturbance Observer for Hypersonic Vehicles[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 45-55.
    [19]Yanlong Zhou, Mou Chen, Changsheng Jiang. Robust Tracking Control of Uncertain MIMO Nonlinear Systems with Application to UAVs[J]. IEEE/CAA Journal of Automatica Sinica, 2015, 2(1): 25-32.
    [20]Wen Qin, Zhongxin Liu, Zengqiang Chen. Formation Control for Nonlinear Multi-agent Systems with Linear Extended State Observer[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(2): 171-179.

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(15)  / Tables(2)

    Article Metrics

    Article views (1024) PDF downloads(36) Cited by()

    Highlights

    • Coordinate Conversion In order to obtain the accurate angular velocity and acceleration to control the speed of the electric bicycle, we calculate and use the real-time data of the carrier coordinate system (pedal) relative to the geographic coordinate system by quaternion-based coordinate conversion.
    • Disturbance observer Considering the variability of frictions, we build a second-order low pass filter with a special cut-off frequency , so that the majority of disturbances can be observed and compensated.
    • Velocity control A permanent magnet brushless DC motor is used for adaptively assisted riding. Based on the real-time angular velocity of the rear wheel, the motor can provide assistance power for controlling the riding speed within a certain range.

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return