汽车侧垂向耦合动力学模型(汽车动力学建模进阶知识)

一.汽车侧垂向耦合动力学模型

1.1 侧向轮胎模型

      这里的汽车侧向轮胎模型,采取的是魔术轮胎半经验轮胎模型,该模型现在在汽车动力学建模领域的应用最为广泛,输入为轮胎法向反力 F z F_z Fz和车轮侧偏角 α \alpha α,输出为车轮侧向力 F y F_y Fy,该经验公式如下:[1]       F y ( α ) = D y sin ⁡ ( C y a r c t a n ( B y α − E y ( B y α − arctan ⁡ ( B y α ) ) ) ) F_y\left(\alpha\right)=D_y\sin{(}C_yarctan\left(B_y\alpha-E_y\left(B_y\alpha-\arctan{\left(B_y\alpha\right)}\right)\right)) Fy(α)=Dysin(Cyarctan(ByαEy(Byαarctan(Byα)))) C y = 1.30 C_y=1.30 Cy=1.30

D y = b 1 F z 2 + b 2 F z D_y=b_1F_z^2+b_2F_z Dy=b1Fz2+b2Fz

B y = b 3 s i n ( b 4 a r c t a n ( b 5 F z ) ) / C y D y B_y=b_3sin{\left(b_4arctan{\left(b_5F_z\right))}\right.}/{C_yD_y} By=b3sin(b4arctan(b5Fz))/CyDy

E x = b 6 F z 2 + b 7 F z + b 8 E_x=b_6F_z^2+b_7F_z+b_8 Ex=b6Fz2+b7Fz+b8
      

拟合参数[1]
C y C_y Cy 1.30
b 1 b_1 b1 -23.4
b 2 b_2 b2 981
b 3 b_3 b3 1078
b 4 b_4 b4 1.67
b 5 b_5 b5 0.202
b 6 b_6 b6 0.000
b 7 b_7 b7 -0.387
b 8 b_8 b8 0.707

      确定好各种参数后,搭建模型即可,这里建议用Function函数,在Matlab内编写程序控制输入输出即可。
汽车侧垂向耦合动力学模型(汽车动力学建模进阶知识)_第1张图片

1.2 侧垂向耦合动力学模型

汽车侧垂向耦合动力学模型(汽车动力学建模进阶知识)_第2张图片
汽车侧垂向耦合动力学模型(汽车动力学建模进阶知识)_第3张图片
      如上图一,为标准的汽车侧向动力学模型,其中左右两个方块代表车轮模型, x x x代表汽车纵向运动方向, y y y代表汽车侧向的运动方向, ψ ˙ \dot{\psi} ψ˙为汽车的横摆角速度; x ˙ \dot{x} x˙为汽车速度纵向分量, y ˙ \dot{y} y˙为汽车的速度侧向分量; β \beta β为质心侧偏角; F y f F_{yf} Fyf为汽车前轮侧向力, F y r F_{yr} Fyr为汽车后轮侧向力; a a a b b b分别为质心到前后轴的距离; α f \alpha_f αf α r \alpha_r αr分别为前后轮的侧偏角; Y a w Y_{aw} Yaw为汽车横摆角速度; δ \delta δ为前轮转角输入。
      如上图二,即为简化的半车模型。该图形中各个参数的含义为 c s r c_{sr} csr c s f c_{sf} csf为前后悬架的阻尼系数; F M R r F_{MRr} FMRr F M R f F_{MRf} FMRf是可控阻尼力; k s f k_{sf} ksf k s r k_{sr} ksr为前后悬架的弹簧刚度; q f q_f qf q r q_r qr为前后轮的路面激励, F z f F_{zf} Fzf F z r F_{zr} Fzr为地面给轮胎的法向反力;了解这两张原理图以后,根据牛顿第二定律,可列出以下微分方程[1]。
      整车侧向刚体运动学微分方程:
M t a y = F x f cos ⁡ δ + F y r M_ta_y=F_{xf}\cos{\delta}+F_{yr} Mtay=Fxfcosδ+Fyr      上式中 M t M_t Mt为整车质量, a y a_y ay为总的纵向加速度 ( a y = y ¨ + x ˙ ψ ˙ ) (a_y=\ddot{y}+\dot{x}\dot{\psi}) (ay=y¨+x˙ψ˙) F x f F_{xf} Fxf为前轮纵向力。在该微分方程中 F y r F_{yr} Fyr为后轮侧向力, δ \delta δ为前轮转角输入控制量。
      整车横摆刚体运动学微分方程: I z ψ ¨ = a c o s δ F y f − b F y γ I_z\ddot{\psi}=acos{\delta F_{yf}}-bF_{y\gamma} Izψ¨=acosδFyfbFyγ      簧载质量垂向动力学运动学微分方程: M b ( z ¨ − x θ ˙ ) = − k s f ( z − a θ − z u f ) − k s r ( z + b θ − z u r ) − C s f ( z ˙ − a θ ˙ − z ˙ u f ) − c s r ( z + b θ ˙ − z ˙ u r ) + F M R f + F M R r M_b\left(\ddot{z}-x\dot{\theta}\right)=-k_sf\left(z-a\theta-z_{uf}\right)-k_{sr}\left(z+b\theta-z_{ur}\right)-C_{sf}\left(\dot{z}-a\dot{\theta}-{\dot{z}}_{uf}\right)-c_{sr}\left(z+b\dot{\theta}-{\dot{z}}_{ur}\right)+F_{MRf}+F_{MRr} Mb(z¨xθ˙)=ksf(zaθzuf)ksr(z+bθzur)Csf(z˙aθ˙z˙uf)csr(z+bθ˙z˙ur)+FMRf+FMRr      簧载质量俯仰动力学运动微分方程: I y θ ¨ = a k s f ( z − a θ − z u f ) − b k s r ( z + b θ − z u r ) + a c s f ( z ˙ − a θ ˙ − z ˙ u f ) − b c s r ( z ˙ + b θ ˙ − z ˙ u r ) − a F M R f + b F M R r + h p F y f s i n δ + M b g h d s i n θ I_y\ddot{\theta}=ak_{sf}\left(z-a\theta-z_{uf}\right)-bk_{sr}\left(z+b\theta-z_{ur}\right)+ac_{sf}\left(\dot{z}-a\dot{\theta}-{\dot{z}}_{uf}\right)-bc_{sr}\left(\dot{z}+b\dot{\theta}-{\dot{z}}_{ur}\right)-aF_{MRf}+bF_{MRr}+h_pF_{yf}sin{\delta}+M_bgh_dsin{\theta} Iyθ¨=aksf(zaθzuf)bksr(z+bθzur)+acsf(z˙aθ˙z˙uf)bcsr(z˙+bθ˙z˙ur)aFMRf+bFMRr+hpFyfsinδ+Mbghdsinθ      车轮侧偏角的用下面两个公式求。
α f l = δ − a r c t a n y ˙ + a ψ ˙ x ˙ , α r l = − arctan ⁡ y ˙ − b ψ ˙ x ˙ \alpha_{fl}=\delta-arctan{\frac{\dot{y}+a\dot{\psi}}{\dot{x}}},\alpha_{rl}=-\arctan{\frac{\dot{y}-b\dot{\psi}}{\dot{x}}} αfl=δarctanx˙y˙+aψ˙αrl=arctanx˙y˙bψ˙      上述为汽车侧垂向耦合的动力学模型。本文的内容为学习卢少波博士论文经验所得。本文篇幅有限,如果想深入贯通了解需要在看一下侧向和纵向还有耦合的模型,可以去知网搜索一下该论文。

二.参考文献

[1] 卢少波. 汽车底盘关键子系统及其综合控制策略研究[D].重庆大学,2009.

你可能感兴趣的:(汽车系统动力学及控制,数学建模,matlab,经验分享,算法,矩阵)