力学与实践, 2020, 42(1): 116-118 DOI: 10.6052/1000-0879-19-027

身边力学的趣话

再谈无人自行车的自稳定性

刘延柱1)

上海交通大学工程力学系,上海 200240

REVISIT OF SELF-STABILIZATION OF UNMANNED BICYCLE

LIU Yanzhu1)

Department of Engineering Mechanics, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 1)E-mail:liuyzhc@163.com

收稿日期: 2019-01-16   网络出版日期: 2020-02-20

Received: 2019-01-16   Online: 2020-02-20

作者简介 About authors

摘要

本文讨论无人自行车稳定行驶的实现可能性,对前叉受控的自行车数学模型做稳定性分析。分析表明,要保证自行车直立行驶的渐近稳定性,控制系统必须考虑前叉转轴的阻尼因素,且控制规律必须计入车体侧向倾斜的角速度信息。

关键词: 无人自行车 ; 渐近稳定性 ; Hurwitz判据

Abstract

The realization of an unmanned self-stabilizing bicycle is discussed. A stability analysis is made for a mathematical model of a bicycle controlled by the rotation of the front set. It is shown that the self- running bicycle can be asymptotically stabilized only when the damping factor in the rotation axis of the front set is considered, and the input signals of the control system should contain the angular velocity of the side-leaning of the bicycle.

Keywords: unmanned bicycle ; asymptotic stability ; Hurwitz's criterion

PDF (643KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘延柱. 再谈无人自行车的自稳定性. 力学与实践[J], 2020, 42(1): 116-118 DOI:10.6052/1000-0879-19-027

LIU Yanzhu. REVISIT OF SELF-STABILIZATION OF UNMANNED BICYCLE. MECHANICS IN ENGINEERING[J], 2020, 42(1): 116-118 DOI:10.6052/1000-0879-19-027

2016年Google宣布已造出能自动行驶的无人自行车,引起不小的轰动。但随后声称只是个愚人节玩笑。此前有报导称,2011年荷兰Delft大学的Kooijman曾制造出无人自行车。但他的实验对象是特殊设计的小车,与常规的自行车截然不同[1]。2015年,清华大学自动化系学生在赵明国老师指导下造出了真正的自平衡无人自行车,在全国虚拟仪器大赛中获得唯一的特等奖。2017年在美国奥斯汀的NI WEEK 2017大会中展示,入围 NI Engineering Impact Awards,受到科技界的极大关注[2]

关于自行车稳定性问题的讨论已延续了两百多年。早在19世纪提出的离心力效应至今仍是自行车稳定性的经典解释。车轮的陀螺效应有助于自行车的稳定,但1971年英国人Jones[3]的实验证实陀螺效应的影响极其微弱,其作用远不如与前轮接触点有关的脚轮效应。1982 年,英国人Lowell和McKell建立的自行车数学模型考虑了离心力效应、陀螺效应和脚轮效应等所有因素,还增加了前叉转角随车身侧向倾斜角变化的简单控制规律。基于此模型导出的特征根存在正实部,数值计算也得出振幅不断增大的不稳定振荡的运动轨迹。表明即使考虑所有稳定因素且增加对前叉的自动控制,也不能使自行车保持稳定[4-5]

如上所述,自稳定无人自行车已在近期内成功实现。上述不稳定性的分析结论已被实践彻底否定。实际上凭常识也能判断,自行车能否稳定的重要因素是驾车人的控制技巧。关键在于,控制规律是否真正体现人的驾车技术。本文仍采用相同的数学模型,仅对控制规律作一些修改,对自行车自稳定性问题继续做些探讨。

将自行车视为由车架、前叉和前后轮组成的刚体系。设$O_1$和$O_2$为后轮和前轮与地面的接触点,以$O_1$为原点建立参考坐标系($O_1$-$xyz$),其中$x$轴沿$O_1 O_2$,$y$轴为垂直轴。设绕轴逆时针转过$\theta$角后的位置为($O_1$-$x_1 y_1 z_1)$,$\left( {x_1 ,y_1 }\right)$为车架平面。设$O_1$与$O_2$的距离为$a$,质心$O_{\rm c}$的直立高度为$h$,在$x$轴上的投影与$O_1$的距离为$b$。设车架随$O_1$点以速度$v$沿$x$ 轴匀速平动(图\,1),前轮偏角$\psi$的出现使自行车转为曲线运动(图\,2)。图中$O$为曲率中心,$R = a/\psi $为曲率半径。近似将车体的全部质量$m$集中在$O_{\rm c}$ 点,仅保留$\theta$和$\varphi$的一次项,列写车体在重力$mg$和离心力${mv^2}/R$作用下绕水平轴$x$ 转动的动力学方程。引入参数$\alpha = {bv}/{ha}$, $\beta = {v^2}/{ga}$, $\gamma = g/h$,写作

$\begin{equation} \label{eq1} \ddot {\theta } + \alpha \dot {\psi } + \gamma \left( {\beta \psi - \theta } \right) = 0 \end{equation}$

列写前叉与前轮的组合件绕转轴转动的动力学方程时,为简化数学推导,仍按文献[4]中的假定将倾斜的前叉转轴$y^*$拉直,使与$O_1 O_2$正交。将($O_1$-$x_1 y_1 z_1$)的原点移至$O_2$,绕$y_1 $轴顺时针转过$\psi $角后的位置为($O_2$-$x_2 y_2 z_2$)。 其中$y_2$轴平行于前叉转轴$y^ * $,$x_2$轴沿轮缘在$O_2$点处的切线,$\left( {x_2 ,y_2 }\right)$为前叉的对称平面。考虑地面在$O_1$和$O_2$处对车体和前叉施加沿$z_1 $和$z_2$轴的摩擦力与离心力平衡。且考虑前叉转轴$y^*$与地面的交点$O^*$与$O_2$点不重合的距离$\varDelta$形成的脚轮效应(图\,3)。设车架与前叉之间存在绕转轴$y^*$ 的阻尼力矩

$\begin{equation} \label{eq2} M_{\text{d}y} = C\dot {\psi } \end{equation}$

其中$C$为黏性摩擦系数。考虑控制系统的作用,将陀螺仪量测到的车体侧向倾斜角$\theta$和角速度$\dot {\theta}$的信息输入电机,产生绕前叉转轴的控制力矩$M_{\text{c}y} $,控制规律为

$\begin{align} M_{\text{c}y} = J_y ( {k_1 \theta + k_2 \dot {\theta }}) \end{align}$

其中$J_y $为前叉连同前轮绕转轴$y^*$的惯量矩。引入参数$\mu ={mgb\varDelta }/({J_y a})$,$c = C/{J_y}$,列出前叉的动力学方程

$\begin{equation} \label{eq3} \ddot {\psi } + c\dot {\psi } + \mu \left( {\beta \psi - \theta } \right) - k_1 \theta - k_2 \dot {\theta } = 0 \end{equation}$

方程(1)和(4)的具体推导过程可参阅文献[4,5]。仅对文献中的动力学方程略去影响微弱的陀螺力矩,增加了阻尼力矩$M_{\text{d}y}$和控制力矩$M_{\text{c}y}$。

图1

图1   参考坐标系


图2

图2   运动轨迹的曲率半径


图3

图3   前叉转轴与前轮触地点


方程组(1), (4)的特征方程为

$\begin{align} \lambda ^4 + a_1 \lambda ^4 + a_2 \lambda ^3 + a_3 \lambda ^2 + a_4 = 0 \end{align}$

其中的系数定义为

$\begin{align} \left.\begin{array}{ll} & a_1 = c \\ & a_2 = \alpha k_2 - \gamma + \mu \beta \\ & a_3 = \alpha ({\mu + k_1}) + \gamma({\beta k_2-c}) \\ & a_4 = \beta \gamma k_1 \end{array}\right\} \end{align}$

利用Hurwitz判据判断线性系统的稳定性,所有系数为正值是方程组(1)和(4)零解渐近稳定的必要条件。此条件要求控制系数$k_1$和$k_2$满足

$\begin{align} \left.\begin{array}{ll} & k_1 > 0 \\ & \alpha k_2 > \gamma - \beta \mu \\ & \alpha k_1 + \gamma \beta k_2 > \gamma c - \alpha \mu \end{array}\right\} \end{align}$

如$a_1 a_2 - a_3 > 0$条件也满足,即

$\begin{align} \label{eq4} & c\left( {\alpha k_2 + \mu \beta - \gamma } \right)\left[ {\alpha \left( {\mu + k_1 } \right) + \gamma \left( {\beta k_2 - c} \right)} \right] > \notag\\ & c^2\beta \gamma k_1 + \left[ {\alpha \left( {\mu + k_1 } \right) + \gamma \left( {\beta k_2 - c} \right)} \right]^2\end{align}$

则渐近稳定性充分成立。在$c \ne 0$,$k_2 \ne 0$条件下,借助控制系数$k_1$和$k_2$的调整可使上述诸条件均得到满足。从而证实无人自行车直立行驶状态存在渐进稳定性的可能性。

以上分析表明,利用对前叉的自动控制,实现无人自行车直立行驶的稳定性在理论上并无问题。但必须考虑前叉转轴的阻尼因素,且控制规律中除车体倾斜的角度信息以外,还必须计入角速度信息。Lowell和McKell的论文未考虑上述因素,才导致不稳定的计算结果[4-5]

直立稳定性只是无人自行车的关键问题之一。要真正实现无人自行车,还必须攻克转弯、过障碍、跟踪等许多技术难关。清华大学的老师和学生们不畏艰难,坚持不懈,终于造出了无人驾驶自行车,震惊了世界,也开启了自行车发展史的新阶段。

参考文献

刘延柱 .

自稳定的无人自行车

力学与实践, 2015,37(1):146-148

[本文引用: 1]

Liu Yanzhu .

Self-stable bicycle without rider

Mechanics in Engineering, 2015,37(1):146-148 (in Chinese)

[本文引用: 1]

https://www.eet-china.com/news/201705231414.html

URL     [本文引用: 1]

Jones DEH .

The stability of the bicycle

Physics Today, 1970,23(4):34-40

[本文引用: 1]

Lowell J, McKell HD .

The stability of bicycles

Amer. J. Physics, 1982,50(12):1106-1112

[本文引用: 4]

刘延柱 .

关于自行车的稳定性

力学与实践, 2012,34(2):90-93

[本文引用: 3]

Liu Yanzhu .

On the stability of bicycle

Mechanics in Engineering, 2012,34(2):90-93 (in Chinese)

[本文引用: 3]

/

Baidu
map