自动驾驶-车辆运动学模型


1、车辆运动学模型

车辆运动学模型是描述车辆位置、速度和加速度等运动状态的模型,主要关注车辆的位置、速度、前轮转角等的关系,而不考虑任何力的影响。一般考虑运动学模型时,将车辆模型简化成单车模型 (bicycle model) 。

在这里插入图片描述

自行车模型(Bicycle Model)的相关假设条件包括:

  1. 不考虑车辆在垂直方向(Z轴方向)的运动,即假设车辆的运动是一个二维平面(X, Y)上的运动。
  2. 假设车辆左右侧轮胎在任意时刻都拥有相同的转向角度和转速;车辆左右两个轮胎的运动可以合并为一个轮胎来描述。
  3. 假设车辆行驶速度变化缓慢,忽略前后轴载荷的转移。
  4. 假设车身和悬架系统都是刚性系统。

单车模型中的符号和定义如下:

A A A:前轮中心点; B B B:后轮中心点;(左右轮都合并为单个轮子)
C C C:车辆重心点;
δ f \delta _{f} δf:前轮转向角; δ r \delta _{r} δr:后轮转向角;
l f {l}_{f} lf:前轮中心点到车辆重心点的距离; l r {l}_{r} lr:后轮中心点到车辆重心点的距离;
L L L:车辆轴距, L = l f + l r L={l}_{f}+{l}_{r} L=lf+lr
V V V:车辆重心的速度;
β \beta β:滑移角,车辆重心速度方向与车辆纵向轴的夹角;
O O O:车辆的瞬时旋转中心;
R R R:车辆的轨迹半径;
Ψ \Psi Ψ:车辆的航向角;

假设速度矢量 V V V 的方向在 A A A 点和 B B B 点的方向与转向角的方向相同,换句话说,在 A A A 点的速度矢量与车辆纵轴的夹角为 δ f \delta _{f} δf,同样 B B B 点的速度矢量与车辆纵轴的夹角为 δ r \delta _{r} δr,也就是说前后轮的滑移角都为 0 0 0。该条件假设成立前提的是车辆速度很低 ( < 5 m / s <5m/s <5m/s),此时轮胎产生的横向力很小,可以忽略。

车辆运动模型基于单车模型推导,推导过程不考虑车辆受到的横向力,故该模型只适用于车辆速度很低的情形。

1.1 运动方程推导

姿态信息 (偏航角)

三角形 O C A OCA OCA 根据正弦定理可得:
sin ⁡ ( δ f − β ) l f = sin ⁡ ( π 2 − δ f ) R (1) \begin{aligned} \frac{\sin\left ( \delta _{f}-\beta \right )}{{l}_{f}}=\frac{\sin\left ( \frac{\pi }{2}- \delta _{f} \right )}{R} \end{aligned}\tag{1} lfsin(δfβ)=Rsin(2πδf)(1)
三角形 O C B OCB OCB 根据正弦定理可得:
sin ⁡ ( β − δ r ) l r = sin ⁡ ( π 2 + δ r ) R (2) \begin{aligned} \frac{\sin\left ( \beta - \delta _{r} \right )}{{l}_{r}}=\frac{\sin\left ( \frac{\pi }{2}+ \delta _{r} \right )}{R} \end{aligned}\tag{2} lrsin(βδr)=Rsin(2π+δr)(2)
对等式 ( 1 ) 、 ( 2 ) (1)、(2) (1)(2) 分别运用和差定理得:
sin ⁡ δ f cos ⁡ β − cos ⁡ δ f sin ⁡ β l f = cos ⁡ δ f R (3) \begin{aligned} \frac{\sin \delta _{f} \cos \beta -\cos \delta _{f}\sin \beta }{{l}_{f}}=\frac{ \cos \delta _{f}}{R} \end{aligned}\tag{3} lfsinδfcosβcosδfsinβ=Rcosδf(3)
sin ⁡ β cos ⁡ δ r − cos ⁡ β sin ⁡ δ r l r = cos ⁡ δ r R (4) \begin{aligned} \frac{\sin\beta \cos\delta _{r} -\cos\beta \sin \delta _{r} }{{l}_{r}}=\frac{ \cos \delta _{r}}{R} \end{aligned}\tag{4} lrsinβcosδrcosβsinδr=Rcosδr(4)
等式 ( 3 ) 、 ( 4 ) (3)、(4) (3)(4) 两边分别同乘 l f cos ⁡ δ f \frac{{l}_{f}}{\cos \delta _{f}} cosδflf l r cos ⁡ δ r \frac{{l}_{r}}{\cos \delta _{r}} cosδrlr得:
tan ⁡ δ f cos ⁡ β − sin ⁡ β = l f R (5) \begin{aligned} \tan \delta _{f} \cos \beta -\sin \beta =\frac{{l}_{f}}{R} \end{aligned}\tag{5} tanδfcosβsinβ=Rlf(5)
sin ⁡ β − tan ⁡ δ r cos ⁡ β = l r R (6) \begin{aligned} \sin \beta - \tan \delta _{r} \cos \beta =\frac{{l}_{r}}{R} \end{aligned}\tag{6} sinβtanδrcosβ=Rlr(6)
对等式 ( 5 ) 、 ( 6 ) (5)、(6) (5)(6) 相加得:
( tan ⁡ δ f − tan ⁡ δ r ) cos ⁡ β = l f + l r R (7) \begin{aligned} \left(\tan \delta _{f} - \tan \delta _{r}\right) \cos \beta=\frac{{l}_{f}+{l}_{r}}{R} \end{aligned}\tag{7} (tanδftanδr)cosβ=Rlf+lr(7)
根据假设条件,在车辆速度很低时,车辆航向角的变化率 ( Ψ ˙ \dot{\Psi} Ψ˙) 可近似等于车辆的角速度 ( ω \omega ω),根据车辆的角速度公式 ω = V R \omega = \frac{V}{R} ω=RV,可得
Ψ ˙ = V R (8) \begin{aligned} \dot{\Psi}= \frac{V}{R} \end{aligned}\tag{8} Ψ˙=RV(8)
将等式(7)带入等式(8)中,消除 R R R 项得:
Ψ ˙ = V cos ⁡ β l f + l r ( tan ⁡ δ f − tan ⁡ δ r ) (9) \begin{aligned} \dot{\Psi}= \frac{V\cos \beta}{{l}_{f}+{l}_{r}}\left(\tan \delta _{f} - \tan \delta _{r}\right) \end{aligned}\tag{9} Ψ˙=lf+lrVcosβ(tanδftanδr)(9)
等式(9)中共有三个输入变量: δ f 、 δ r 、 V \delta _{f}、\delta _{r}、V δfδrV

位置信息

X ˙ = V cos ⁡ ( β + Ψ ) (10) \begin{aligned} \dot{X} =V\cos \left ( \beta +\Psi \right ) \end{aligned}\tag{10} X˙=Vcos(β+Ψ)(10)
Y ˙ = V sin ⁡ ( β + Ψ ) (11) \begin{aligned} \dot{Y} =V\sin \left ( \beta +\Psi \right ) \end{aligned}\tag{11} Y˙=Vsin(β+Ψ)(11)

滑移角 β \beta β 推导

等式 ( 5 ) 、 ( 6 ) (5)、(6) (5)(6) 两边分别同乘 l r {l}_{r} lr l f {l}_{f} lf得:
l r tan ⁡ δ f cos ⁡ β − l r sin ⁡ β = l f ∗ l r R (12) \begin{aligned} {l}_{r}\tan \delta _{f} \cos \beta - {l}_{r}\sin \beta =\frac{{l}_{f} \ast {l}_{r}}{R} \end{aligned}\tag{12} lrtanδfcosβlrsinβ=Rlflr(12)
l f sin ⁡ β − l f tan ⁡ δ r cos ⁡ β = l f ∗ l r R (13) \begin{aligned} {l}_{f}\sin \beta - {l}_{f}\tan \delta _{r} \cos \beta =\frac{{l}_{f} \ast {l}_{r}}{R} \end{aligned}\tag{13} lfsinβlftanδrcosβ=Rlflr(13)
等式 (12)、(13) 相减得:
cos ⁡ β ( l f tan ⁡ δ r + l r tan ⁡ δ f ) = ( l f + l r ) sin ⁡ β (14) \begin{aligned} \cos \beta \left ( {l}_{f}\tan \delta _{r}+{l}_{r}\tan \delta _{f} \right ) =\left ( {l}_{f} + {l}_{r} \right ) \sin \beta \end{aligned}\tag{14} cosβ(lftanδr+lrtanδf)=(lf+lr)sinβ(14)
等式 (14) 两端同乘 1 ( l f + l r ) cos ⁡ β \frac{1}{\left ( {l}_{f} + {l}_{r} \right ) \cos \beta} (lf+lr)cosβ1得:
tan ⁡ β = l f tan ⁡ δ r + l r tan ⁡ δ f l f + l r (15) \begin{aligned} \tan \beta =\frac{ {l}_{f}\tan \delta _{r}+{l}_{r}\tan \delta _{f} }{ {l}_{f} + {l}_{r} } \end{aligned}\tag{15} tanβ=lf+lrlftanδr+lrtanδf(15)
等式 (15) 取反三角函数得:
β = arctan ⁡ ( l f tan ⁡ δ r + l r tan ⁡ δ f l f + l r ) (16) \begin{aligned} \beta =\arctan \left (\frac{ {l}_{f}\tan \delta _{r}+{l}_{r}\tan \delta _{f} }{ {l}_{f} + {l}_{r} } \right) \end{aligned}\tag{16} β=arctan(lf+lrlftanδr+lrtanδf)(16)

微分方程形式

根据等式(9)、(10)、(11)、(16),对于前轮转向系统, δ r = 0 \delta _{r}=0 δr=0,可得:
X ˙ = V cos ⁡ ( β + Ψ ) Y ˙ = V sin ⁡ ( β + Ψ ) Ψ ˙ = V cos ⁡ β tan ⁡ δ f L β = arctan ⁡ l r tan ⁡ δ f L \begin{aligned} \dot{X} =V\cos \left ( \beta +\Psi \right ) \\ \dot{Y} =V\sin \left ( \beta +\Psi \right ) \\ \dot{\Psi}= \frac{V\cos \beta\tan \delta _{f} }{L} \\ \beta =\arctan \frac{ {l}_{r}\tan \delta _{f} }{ L } \\ \end{aligned} X˙=Vcos(β+Ψ)Y˙=Vsin(β+Ψ)Ψ˙=LVcosβtanδfβ=arctanLlrtanδf
当车辆重心点在后轴中心点时, β ≅ 0 \beta \cong 0 β0,上述微分方程变为:
X ˙ = V cos ⁡ Ψ Y ˙ = V sin ⁡ Ψ Ψ ˙ = V tan ⁡ δ f L \begin{aligned} \dot{X} =V\cos \Psi \\ \dot{Y} =V\sin \Psi \\ \dot{\Psi}= \frac{V\tan \delta _{f} }{L} \\ \end{aligned} X˙=VcosΨY˙=VsinΨΨ˙=LVtanδf

1.2 补充-阿克曼转向几何

在这里插入图片描述
在单车模型中,假设转向时的左、右前轮偏角为同一角度,但实际两个角度并不完全相等。通常情况下,内侧轮胎转角更大,如上图所示。
上图中, l w {l}_w lw为车轴宽度; δ o 、 δ i \delta _{o}、\delta _{i} δoδi分别表示车辆前轮外侧、内侧转角; L = l f + l r L = {l}_f+{l}_r L=lf+lr为车辆轴距,其值远远小于车辆的轨迹半径 R R R,滑移角 β \beta β接近于 0 0 0
在车辆后轴为固定轴的情况下, δ r = 0 \delta _{r}=0 δr=0,等式(9)可以近似为:
Ψ ˙ = V tan ⁡ δ f L (17) \begin{aligned} \dot{\Psi}= \frac{V\tan \delta _{f} }{L} \end{aligned}\tag{17} Ψ˙=LVtanδf(17)
由于 δ f \delta _{f} δf很小,故 tan ⁡ δ f ≈ δ f \tan\delta _{f}\approx\delta _{f} tanδfδf,根据等式(8)、(17)可得:
Ψ ˙ V ≈ δ f L = 1 R (18) \begin{aligned} \frac{\dot{\Psi}}{V}\approx \frac{\delta _{f} }{L}=\frac{1}{R} \end{aligned}\tag{18} VΨ˙Lδf=R1(18)
不区分车辆前后轴时,等效转向角为:
δ = L R (19) \begin{aligned} \delta= \frac{L}{R} \end{aligned}\tag{19} δ=RL(19)
由于内外轮转弯半径不同,根据等式 (19),车辆外轮转角、内轮转角分别为:
δ o = L R + l w 2 (20) \begin{aligned} \delta_o= \frac{L}{R+\frac{{l}_w}{2}} \end{aligned}\tag{20} δo=R+2lwL(20)
δ i = L R − l w 2 (21) \begin{aligned} \delta_i= \frac{L}{R-\frac{{l}_w}{2}} \end{aligned}\tag{21} δi=R2lwL(21)
车辆前轮平均转向角为:
δ = δ o + δ i 2 = L R − l w 2 4 R (22) \begin{aligned} \delta= \frac{\delta_o+\delta_i}{2}= \frac{L}{R-\frac{{l}_w^{2}}{4R}} \end{aligned}\tag{22} δ=2δo+δi=R4Rlw2L(22)
由于 l w {l}_w lw远远小于 R R R,故 l w 2 4 R ≈ 0 \frac{{l}_w^{2}}{4R}\approx0 4Rlw20,等式(22)可近似为:
δ = L R (23) \begin{aligned} \delta= \frac{L}{R} \end{aligned}\tag{23} δ=RL(23)
车辆内外轮转角之差 Δ δ \Delta\delta Δδ为:
Δ δ = δ i − δ o = L ⋅ l w R 2 − l w 2 4 ≅ L R 2 l w = δ 2 l w L (24) \begin{aligned} \Delta\delta= \delta_i-\delta_o= \frac{L\cdot {l}_w}{R^{2}-\frac{ {l}_w^{2}}{4}}\cong \frac{L}{R^{2}}{l}_w=\delta^{2}\frac{{l}_w}{L} \end{aligned}\tag{24} Δδ=δiδo=R24lw2LlwR2Llw=δ2Llw(24)
根据等式 (24) 可知,前轮内外转向角的差值 Δ δ \Delta\delta Δδ 与平均转向角 δ \delta δ 的二次方成正比,所以前轮转向角越大,内外轮的转向角误差就越大。


参考文献:

1、控制算法 - 运动学模型
2、自动驾驶中的车辆运动学模型

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值