平衡小车动力学建模

文章详细分析了两轮自平衡小车的力学模型,建立了车轮和车体的动力学方程,包括正向运动和转向运动的状态方程。通过牛顿第二定律和刚体转动定律,推导出车轮和车体的加速度与电机扭矩、摩擦力之间的关系,最终得到了系统的状态空间表达式,为小车的控制提供了理论基础。

前言

两轮自平衡小车由车体和双轮两部分组成,可以看成一个移动的倒立摆,分别对车轮和车体进行力学分析,建立动力学模型,最后,通过对两者的分析给出系统的状态空间表达式。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-pURItxa1-1679556611143)(C:\Users\DELL\AppData\Roaming\Typora\typora-user-images\image-20230322110521545.png)]

车轮模型

以右轮为例进行受力分析,如图所示
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-JGfjgfh3-1679556611143)(C:\Users\DELL\AppData\Roaming\Typora\typora-user-images\image-20230322110550417.png)]

将车轮的运动分解为平动和转动,
由牛顿第二定律可得
mxR¨=HfR−HR(1)m\ddot{x_R}=H_{fR}-H_{R} {(1)}mxR¨=HfRHR(1)

由刚体定轴转动定律可得
IωR˙=TR−HfRr(2)I\dot{\omega_R}=T_R-H_{fR}r {(2)}IωR˙=TRHfRr(2)

其中,

物理量 描述 单位
mmm 车轮的质量 kgkgkg
rrr 车轮的半径 mmm
xRx_RxR 右轮的水平位移 mmm
HfRH_{fR}HfR 右轮受到地面的摩擦力的大小 NNN
HRH_{R}HR 右轮受到车体作用力的水平分力的大小 NNN
TRT_{R}TR 右轮电机输出的转矩的大小 N⋅mN\cdot mNm
III 车轮的转动惯量 kg⋅m2kg\cdot m^2kgm2
ωR\omega_RωR 右轮的角速度的大小 rad/srad/srad/s

联立(1)(2)(1)(2)(1)(2),消去HfRH_{fR}HfR,可得

mxR¨=TR−IωR˙r−HR(3)m\ddot{x_R}=\frac{T_R-I\dot{\omega_R}}{r}-H_{R} {(3)}mxR¨=rTRIωR˙HR(3)

在车轮不打滑的情况下,车轮的移动速度的大小和转动速度的大小成比例关系,即
{ ωR=xR˙rωR˙=xR¨r(4) \begin{cases} \omega_R=\frac{\dot{x_R}}{r} \\ \dot{\omega_R}=\frac{\ddot{x_R}}{r} \end{cases} {(4)} { ωR=rxR˙ωR˙=rxR¨(4)
将方程(4)(4)(4)代入(3)(3)(3)中,消去ωR˙\dot{\omega_R}ωR˙,可得
(m+Ir2)xR¨=TRr−HR(5) (m+\frac{I}{r^2})\ddot{x_R}=\frac{T_R}{r}-H_{R} {(5)} (m+r2I)xR¨=rTRHR(5)
由于左右轮的参数相同,对左轮有类似的结果,即
(m+Ir2)xL¨=TLr−HL(6) (m+\frac{I}{r^2})\ddot{x_L}=\frac{T_L}{r}-H_{L} {(6)} (m+r2I)xL¨=rTLHL(6)

车轮模型

车体的运动也可以分解为正向运动(前向、俯仰)和侧向运动(转向、偏航)。

正向运动

小车的正向运动可以分解为前向运动和绕车体质心PPP的相对转动(俯仰)。小车底盘中心OOO的水平位移为x=xL+xR2(7)x=\frac{x_L+x_R}{2} {(7)}x=2xL+xR(7)

将方程(5)(6)(5)(6)(5)(6)相加,等式两边同时除以222可得
(m+Ir2)xL¨+xR¨2=TL+TR2r−HL+HR2(8) (m+\frac{I}{r^2})\frac{\ddot{x_L}+\ddot{x_R}}{2}=\frac{T_L+T_R}{2r}-\frac{H_{L}+H_{R}}{2}{(8)} (m+r2I)2xL¨+xR¨=2rTL+TR2HL+HR(8)
(7)(7)(7)代入(8)(8)(8)可得
(m+Ir2)x¨=TL+TR2r−HL+HR2(9) (m+\frac{I}{r^2})\ddot{x}=\frac{T_L+T_R}{2r}-\frac{H_{L}+H_{R}}{2}{(9)} (m+r2I)x¨=2rTL+TR2HL+HR(9)

在这里插入图片描述

对车体,应用牛顿第二定律可得

在水平方向上,有
Md(x+lsinθP)dt2=HL+HR(10) M\frac{d(x+lsin\theta_P)}{dt^2}=H_L+H_R{(10)} Mdt2d(x+lsinθP)=HL+HR(10)
在竖直方向上,有
Md(lcosθP)dt2=VL+VR−Mg(11) M\frac{d(lcos\theta_P)}{dt^2}=V_L+V_R-Mg{(11)} Mdt2

评论 18
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

u小鬼

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值