很多刚接触车辆动力学仿真的朋友一上来就问我“Simulink里怎么搭一个整车模型”这个问题其实挺难回答的因为整车模型本身是个很庞大的工程涉及悬架、轮胎、转向、制动等一大堆子系统。但如果把问题聚焦在“理解车辆运动的基本规律”上那答案就简单多了——从二自由度和三自由度模型入手。这两个模型虽然结构简单但把车辆横向、纵向、横摆这几个最核心的运动状态都涵盖进去了是理解ESP、ABS、四轮转向这些上层控制算法的基础。这篇文章就是写给初学者的我会从模型本身的物理意义讲起推导过程尽量写得通俗然后给出Simulink里的具体搭建方法包括直接积分法、状态空间法以及从二自由度扩展到三自由度时要注意的问题。文章里所有内容我都用实际跑过的模型验证过参数也给了参考值你照着搭就能跑起来。1. 为什么先学二自由度模型操纵稳定性分析的最小完整系统很多初学者容易犯一个错误就是上来就追求模型“越复杂越真实”结果模型建到一半参数标定出了问题仿真结果一团糟还不知道是哪里错了。我个人的建议是先踏踏实实把二自由度模型吃透它虽然舍弃了很多细节但保留了车辆操纵稳定性分析中最核心的矛盾——轮胎侧偏力与横摆运动的耦合关系。1.1 二自由度模型到底忽略了什么保留了什举二自由度模型通常指包含横摆角速度和质心侧偏角两个自由度的“自行车模型”。它做了一系列人为简化忽略悬架运动认为车身只做平面运动忽略空气阻力、滚动阻力等纵向力把左右车轮合并为一个车轮所以叫“自行车模型”认为轮胎工作在线性区即侧偏力与侧偏角成线性关系。这些简化听起来“很不真实”但恰恰是这些简化让我们能用一组常系数线性微分方程来描述车辆运动。而线性系统的好处是可以用解析方法求解可以用状态空间理论分析稳定性这也是后续所有车辆稳定性控制算法的理论基石。举一个具体的例子ESP系统在判断车辆是否失稳时核心依据就是“实际横摆角速度与理想横摆角速度的偏差”。而这个“理想横摆角速度”从哪里来就是从二自由度模型的稳态响应里推出来的。你说这个简单模型重不重要1.2 两个自由度各自的物理含义第一个自由度是横摆角速度记为 \(\omega_r\)描述的是车辆绕垂直轴旋转的快慢单位是rad/s。你开车打方向盘车头开始转动这个转动的快慢就是横摆角速度。它是衡量车辆“灵活性”最直接的指标。第二个自由度是质心侧偏角记为 \(\beta\)描述的是车辆质心速度方向与车身纵轴方向的夹角。这个量稍微抽象一点你可以这样理解车头指向的方向和车实际上在移动的方向之间存在一个夹角。正常直线行驶时这个角接近零高速过弯时这个角会明显增大。质心侧偏角过大说明车辆正在失去轨迹保持能力这是非常危险的工况。这两个自由度不是独立的它们通过轮胎侧偏力互相耦合——横摆角速度影响前轮侧偏角前轮侧偏角产生的前轮侧偏力又会影响质心侧偏角的变化率而质心侧偏角又反过来影响后轮侧偏角。这一组耦合关系就是二自由度模型的核心内容。2. 从牛顿第二定律到状态方程二自由度模型的完整推导我见过不少同学Simulink模型搭得很熟练但你问他“这个增益模块里的数值是怎么算出来的”他就答不上来了。这是很危险的——模型的正确性无法验证参数一改就出错。所以这一节我们把推导过程完整走一遍我保证每一步都写得足够清楚。2.1 坐标系约定与关键变量定义在推导之前先约定坐标系。采用车辆坐标系x轴沿车辆纵轴向前y轴垂直于纵轴指向左侧或者右侧不同教材约定不同我这里采用向左为正但这不影响最终结果的形式z轴垂直向上。坐标原点位于质心。符号含义单位\(m\)整车质量kg\(V\)质心速度沿轨迹方向m/s\(\beta\)质心侧偏角rad\(\omega_r\)横摆角速度rad/s\(\delta\)前轮转角rad\(k_1\)前轮侧偏刚度负值N/rad\(k_2\)后轮侧偏刚度负值N/rad\(a\)质心到前轴距离m\(b\)质心到后轴距离m\(I_z\)绕z轴的转动惯量kg·m²\(L\)轴距\(L a b\)m这里需要提醒一下侧偏刚度的符号约定在不同的教材里不一样。有些教材定义为正值有些定义为负值因为侧偏力与侧偏角方向相反。我这里采用负值约定即 \(k_1 0\)这在后续推导里会体现出来。如果读者看到其他教材里的符号和你预期的不一样请先确认它的符号约定这是初学者最容易搞混的地方。2.2 侧偏角的几何关系前轮和后轮各有各的侧偏角轮胎侧偏角的定义是轮胎接地印迹中心的速度方向与轮胎旋转平面即车轮平面之间的夹角。我们要做的第一步是把前轮和后轮的侧偏角用状态变量 \(\beta\)、\(\omega_r\) 和前轮转角 \(\delta\) 表示出来。先看前轮。前轮中心点的速度由两部分组成质心速度 \(V\) 加上由于横摆运动产生的相对速度。前轮中心的横向速度为\[ v_{y1} V\sin\beta a\omega_r \]纵向速度为\[ v_{x1} V\cos\beta \]前轮中心处的速度方向角为 \(\arctan\left(\frac{V\sin\beta a\omega_r}{V\cos\beta}\right)\)。在侧偏角较小的情况下\(\beta\) 和 \(\omega_r\) 都不大可以近似为\[ \arctan\left(\frac{V\sin\beta a\omega_r}{V\cos\beta}\right) \approx \beta \frac{a\omega_r}{V} \]而前轮的速度方向与轮胎平面之间的夹角还要减去前轮转角 \(\delta\)因为轮胎平面已经相对车身转了 \(\delta\) 角。所以前轮侧偏角为\[ \alpha_1 \delta - \left(\beta \frac{a\omega_r}{V}\right) \]注意这里符号的选择当 \(\delta \beta \frac{a\omega_r}{V}\) 时\(\alpha_1 0\)此时前轮侧偏力为负指向右侧在左侧为正的约定下产生顺时针横摆力矩。这个细节很重要后面列方程时要保持符号一致。再看后轮。后轮没有转角所以\[ \alpha_2 -\left(\beta - \frac{b\omega_r}{V}\right) \]这里的负号来自速度方向的近似计算。很多初学者容易漏掉这个负号导致后轮侧偏力的方向搞反。2.3 力与力矩方程两个方程两个未知数有了侧偏角就可以计算侧偏力了\[ F_{y1} k_1\alpha_1 k_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) \]\[ F_{y2} k_2\alpha_2 k_2\left(-\beta \frac{b\omega_r}{V}\right) \]接下来对车辆列写牛顿第二定律。沿y轴方向的力平衡方程\[ m\left(\dot{V}\sin\beta V\cos\beta\cdot\dot{\beta}\right) F_{y1}\cos\delta F_{y2} \]在小角度假设下\(\cos\beta \approx 1\)\(\sin\beta \approx \beta\)且认为 \(\dot{V} 0\)匀速行驶上式简化为\[ mV\left(\dot{\beta} \omega_r\right) F_{y1} F_{y2} \]这里用到了一个重要关系质心速度方向角的变化率等于横摆角速度。具体推导是\(\dot{\beta} \omega_r\) 实际上是车辆加速度在y轴方向的分量除以速度后得到的“横向加速度变化率”形式。初学者可以先记住这个结论理解时抓住“横向动力学”的本质即可。绕z轴的力矩平衡方程\[ I_z\dot{\omega}r aF{y1}\cos\delta - bF_{y2} \]小角度下 \(\cos\delta \approx 1\)所以\[ I_z\dot{\omega}r aF{y1} - bF_{y2} \]注意后轮侧偏力那一项的符号后轮侧偏力方向和前轮相反所以对质心产生的力矩方向也相反用减号。2.4 整理成状态空间形式把 \(F_{y1}\) 和 \(F_{y2}\) 的代达式代入整理成关于 \(\dot{\beta}\) 和 \(\dot{\omega}_r\) 的方程\[ mV(\dot{\beta} \omega_r) k_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) k_2\left(-\beta \frac{b\omega_r}{V}\right) \]展开并移项\[ mV\dot{\beta} mV\omega_r -(k_1 k_2)\beta - \frac{ak_1 - bk_2}{V}\omega_r k_1\delta \]\[ \dot{\beta} \frac{k_1 k_2}{mV}\beta \left(\frac{ak_1 - bk_2}{mV^2} - 1\right)\omega_r \frac{k_1}{mV}\delta \]第二个方程\[ I_z\dot{\omega}_r ak_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) - bk_2\left(-\beta \frac{b\omega_r}{V}\right) \]展开整理\[ \dot{\omega}_r \frac{ak_1 - bk_2}{I_z}\beta \frac{a^2k_1 b^2k_2}{I_z V}\omega_r \frac{ak_1}{I_z}\delta \]写成矩阵形式就是标准的二自由度车辆模型状态方程\[ \begin{bmatrix} \dot{\beta} \\ \dot{\omega}_r \end{bmatrix}\begin{bmatrix} \frac{k_1 k_2}{mV} \frac{ak_1 - bk_2}{mV^2} - 1 \\ \frac{ak_1 - bk_2}{I_z} \frac{a^2k_1 b^2k_2}{I_z V} \end{bmatrix} \begin{bmatrix} \beta \\ \omega_r \end{bmatrix} \begin{bmatrix} \frac{k_1}{mV} \\ \frac{ak_1}{I_z} \end{bmatrix} \delta \]这套方程就是整个二自由度模型的核心。后面在Simulink里搭建时本质上就是在用积分器、增益模块和加法模块来“复现”这两个方程。3. Simulink里的第一种搭法按微分方程直接搭积分模块现在进入实操环节。我先介绍最直观、最适合初学者理解的一种方式——直接根据状态方程搭建积分模块。这种方式的好处是每个模块对应方程里的一项你能清楚地看到“哪个增益是干什么的”对理解模型本身帮助很大。3.1 从方程到方块图的映射思路先看状态方程的结构。第一个方程告诉我们 \(\dot{\beta}\) 由三部分线性叠加而成一项是 \(\beta\) 的倍数一项是 \(\omega_r\) 的倍数还有一项是 \(\delta\) 的倍数。第二个方程同理。所以整个Simulink模型的基本骨架应该是两个积分器分别输出 \(\beta\) 和 \(\omega_r\)积分器的输入分别是 \(\dot{\beta}\) 和 \(\dot{\omega}_r\)\(\dot{\beta}\) 由三个支路求和得到\(\dot{\omega}_r\) 也是由三个支路求和得到支路来自状态变量经过增益模块的输出以及外部输入 \(\delta\) 经过增益模块的输出。用一句话概括就是“状态变量的导数由状态变量和输入的线性组合决定”这本质上是线性系统理论里“状态方程决定状态导数”的直接体现。想通了这一点整个模型结构的草图就能在脑子里画出来。3.2 参数设置用一组真实感强的车辆参数在开始搭建之前先给出一组参考参数。这些参数来自我平时用的一辆类似家用轿车的简化数据量级是合理的你可以直接拿来用。参数值单位\(m\)1500kg\(I_z\)2500kg·m²\(a\)1.2m\(b\)1.4m\(k_1\)-60000N/rad\(k_2\)-80000N/rad\(V\)20m/s约72km/h\(L\)2.6m在Simulink里增益模块的参数就是直接填入根据这些值计算出的数字。你在建模时可以把这些增益的计算用MATLAB脚本算好也可以直接在增益模块里写入表达式比如(k1k2)/(m*V)。我建议在模型里用变量名然后在模型回调或MATLAB脚本里赋值这样改参数时不用改模型更灵活。3.3 搭建步骤详解打开Simulink新建一个空白模型。我按步骤带你走一遍。第一步添加积分器和输出接口。从Simulink库浏览器里拖入两个Integrator模块。为清晰起见把第一个积分器重命名为“beta”第二个重命名为“omega_r”。这不仅仅是命名规范的问题更重要的是当模型越来越复杂时清晰的命名能让你快速定位。另外给积分器设置初始条件\(\beta(0) 0\)\(\omega_r(0) 0\)对应车辆从直线行驶状态开始运动。第二步构建 \(\dot{\beta}\) 的求和单元。拖入一个Sum模块设置符号为“”三个输入都为正输入数量改为3。三条输入支路分别是\(\beta \times \frac{k_1 k_2}{mV}\)\(\omega_r \times \left(\frac{ak_1 - bk_2}{mV^2} - 1\right)\)\(\delta \times \frac{k_1}{mV}\)这三个增益模块分别连接对应的信号源beta积分器输出、omega_r积分器输出、外部输入delta输出端都接到Sum模块的输入端。第三步构建 \(\dot{\omega}_r\) 的求和单元。同样拖入一个Sum模块符号“”输入数3。增益分别为\(\frac{ak_1 - bk_2}{I_z}\)\(\frac{a^2k_1 b^2k_2}{I_z V}\)\(\frac{ak_1}{I_z}\)三个支路同样取自对应的状态量和输入。注意第二项 \(\omega_r\) 前的系数中 \(a^2k_1 b^2k_2\) 依然是负值因为 \(k_1,k_2\) 都是负值所以积分器omega_r的输出经过这个增益模块后对 \(\dot{\omega}_r\) 的贡献是负反馈。第四步连接积分器形成闭环。这是最关键的一步。Sum模块的输出是 \(\dot{\beta}\)把它接到beta积分器的输入beta积分器的输出再次分支一路作为状态变量引出一路返回给增益模块用于计算 \(\dot{\beta}\) 和 \(\dot{\omega}_r\)。omega_r同理。千万注意积分器输出到增益模块再到求和模块再到积分器输入这个闭环绝对不能断否则状态量无法更新仿真结果就是错的。第五步添加输入和显示输出。用一个Step模块作为前轮转角输入 \(\delta\)步长设为0.02rad约1.15度这是个中等强度的转向输入。输出方面可以接Scope查看 \(\beta\) 和 \(\omega_r\) 的时间历程。如果你想看车辆的行驶轨迹可以额外计算 \(x\) 和 \(y\) 方向的速度再积分。3.4 调试中的常见错误我见过很多初学者在搭建这个模型时遇到以下几个问题提前帮你排掉增益符号写反侧偏刚度用了正值导致整个系统的极点跑到右半平面仿真结果直接发散。记住在我们这套约定下\(k_1, k_2\) 必须为负。闭环没连上积分器输入悬空状态量永远不变输出始终是初始值图形是一条平线。单位搞混角度用弧度而不是度。如果你用角度制输入Step模块结果会大得离谱因为38度的转向输入在实际车辆上是不可能出现的。求解器步长太大如果默认的求解器自动步长下结果出现锯齿形的抖动手动将最大步长调到0.001s左右。4. 更工程化的第二种搭法状态空间模块与S函数的对比直接按方程搭建的方法优点是直观缺点是当你想修改参数、批量仿真时不太方便。工程上更常用的是利用MATLAB强大的矩阵运算能力把整个状态方程封装成一个状态空间模块或者用S函数实现更复杂的变参数仿真。4.1 用State-Space模块一行代码搞定MATLAB/Simulink自带State-Space模块它直接接受A、B、C、D四个矩阵作为参数内部帮你完成连续系统的仿真。对二自由度模型来说这个模块的参数设置比手动搭积木简单得多。% 参数定义 m 1500; % 整车质量 kg I_z 2500; % 转动惯量 kg.m^2 a 1.2; % 质心到前轴距离 m b 1.4; % 质心到后轴距离 m k1 -60000; % 前轮侧偏刚度 N/rad k2 -80000; % 后轮侧偏刚度 N/rad V 20; % 纵向速度 m/s % 状态矩阵 A [(k1k2)/(m*V), (a*k1-b*k2)/(m*V^2)-1; ... (a*k1-b*k2)/I_z, (a^2*k1b^2*k2)/(I_z*V)]; B [k1/(m*V); a*k1/I_z]; C eye(2); D zeros(2,1);把这些矩阵填入State-Space模块的参数对话框中输入信号仍然是 \(\delta\)输出就是 \(\beta\) 和 \(\omega_r\)。这种方式的好处是代码可读性高参数集中管理后续做参数扫描优化时非常方便——只需要在m脚本里改参数然后调用sim()函数批量仿真。4.2 用S函数实现变参数模型如果你的仿真里 \(V\)车速是变化的那二自由度模型就不能再用定常矩阵A、B来描述因为矩阵元素里含有 \(V\)。这时候两种方案一是用“变参数状态空间”的积分方式手动搭建把 \(V\) 作为一个外部输入信号接入增益模块二是用S函数在代码里动态计算 \(A\)、\(B\) 矩阵。用Level-2 MATLAB S函数实现的核心伪代码如下function update(block) % 获取输入信号delta, V delta block.InputPort(1).Data; V block.InputPort(2).Data; % 从对话框参数中获取固定参数 m block.DialogPrm(1).Data; Iz block.DialogPrm(2).Data; a block.DialogPrm(3).Data; b block.DialogPrm(4).Data; k1 block.DialogPrm(5).Data; k2 block.DialogPrm(6).Data; % 连续状态方程 f1 (k1k2)/(m*V) * block.ContStates(1) ... ((a*k1-b*k2)/(m*V^2) - 1) * block.ContStates(2) ... k1/(m*V) * delta; f2 (a*k1-b*k2)/Iz * block.ContStates(1) ... (a^2*k1b^2*k2)/(Iz*V) * block.ContStates(2) ... a*k1/Iz * delta; block.Derivatives.Data [f1; f2]; end这里我把车速 \(V\) 也当作输入信号了这样你可以模拟“匀减速过弯”“加速出弯”等更接近真实工况的场景。需要注意的是当 \(V 0\) 时方程奇异矩阵元素分母为零所以仿真起始时刻要设置一个很小的初值比如 \(V 0.1\) m/s或者用条件判断防止除零错误。4.3 两种实现方式的对比与选型建议维度直接搭建积分模块使用State-Space模块使用S函数直观性高每个增益对应一项中矩阵封装了细节低需要阅读代码适合人群初学者理解原理快速搭建、参数扫描变参数、复杂仿真扩展性低状态增加时模型膨胀中修改矩阵即可高可灵活控制调试难度低容易定位问题中矩阵数值错误难查高需要代码调试经验我的建议是学习阶段用第一种把每个增益算一遍、查一遍模型自然就懂了。做课题或项目时用第二种把参数集中到脚本里后续做批量仿真、参数优化、蒙特卡洛分析时省事很多。如果涉及变车速、变载荷等非线性情况直接用S函数一步到位。5. 三自由度模型加入纵向运动后的动力学扩展二自由度模型把车速 \(V\) 当作常数这在分析稳态转向特性时够用了。但如果你想研究车辆在驱动/制动过程中的横摆稳定性或者模拟ESP介入时的整车响应就必须要考虑纵向车速的变化。这时就需要把模型从二自由度扩展到三自由度。5.1 第三个自由度选什么纵向车速 \(V_x\)一种最常见的扩展方式是保留横摆角速度 \(\omega_r\) 和质心侧偏角 \(\beta\)再加上纵向车速 \(V_x\)作为第三个自由度。这样状态变量就是 \([V_x,\ \beta,\ \omega_r]^T\)。相比于二自由度模型它不再假设 \(V\) 为常数而是通过纵向力平衡方程显式地计算 \(V_x\) 的变化。需要注意的是严格来说这里的状态变量采用 \(V_x\)纵向速度而不是标量车速 \(V\) 更合适因为车辆存在侧偏时\(V \sqrt{V_x^2 V_y^2}\)而 \(\beta \arctan(V_y/V_x)\)。在Simulink中要么同时计算 \(V_x\) 和 \(V_y\)要么用 \(V\) 和 \(\beta\) 作为一组广义坐标。下面我以 \(V_x, V_y, \omega_r\) 三者的方式推导这样物理含义更直接。5.2 三自由度模型的运动微分方程在车辆坐标系下列写x方向、y方向的力平衡和绕z轴的力矩平衡。注意因为坐标系本身在旋转加速度项里会出现额外的“哥氏项”。x方向\[ m(\dot{V}x - V_y\omega_r) F{x1} F_{x2} \]y方向\[ m(\dot{V}y V_x\omega_r) F{y1} F_{y2} \]z方向横摆\[ I_z\dot{\omega}r aF{y1} - bF_{y2} \]其中 \(F_{y1}, F_{y2}\) 仍然由侧偏角决定\[ \alpha_1 \delta - \arctan\left(\frac{V_y a\omega_r}{V_x}\right) \]\[ \alpha_2 -\arctan\left(\frac{V_y - b\omega_r}{V_x}\right) \]小角度下近似为\[ \alpha_1 \delta - \frac{V_y a\omega_r}{V_x} \]\[ \alpha_2 -\frac{V_y - b\omega_r}{V_x} \]纵向力 \(F_{x1}, F_{x2}\) 由驾驶员模型或控制策略给定。如果要模拟简单的驱动/制动工况可以直接设定一个纵向力输入如果是整车联合仿真这些纵向力可以由发动机模型或制动系统模型提供。5.3 从二自由度到三自由度的递进理解二自由度模型实际上是三自由度模型在 \(\dot{V}x 0\) 且 \(V_x\) 恒为常数的特殊情况。把三自由度方程中的 \(V_x\) 固定且令 \(F{x1} F_{x2} 0\)就退化为二自由度模型。这个递进关系很重要。初学者往往容易把“自由度”理解为“模型复杂度越高越高级”但实际上更应该把它理解为“对问题分析尺度的选择”。研究稳态转向特性二自由度足够了还便于解析求解研究驱动制动与转向的耦合就必须用三自由度以上模型。选哪个取决于你要回答什么问题而不是追求参数多。5.4 三自由度模型的Simulink实现要点三自由度模型的Simulink搭建和二自由度类似仍然用积分器 增益 求和的方式但增加了几个关键处理额外添加一个积分器用于 \(V_x\)其输入是纵向加速度 \(\dot{V}x (F{x1} F_{x2}) / m V_y\omega_r\)。\(V_x\) 不再是常数而是积分器的输出必须把它作为信号引到侧偏角的计算模块中因为 \(\alpha_1, \alpha_2\) 的分母里有 \(V_x\)。\(\beta\) 不再直接积分输出而是通过 \(V_y\) 和 \(V_x\) 计算\(\beta \arctan(V_y / V_x)\)。因此你可以直接对 \(V_y\) 积分然后在需要输出质心侧偏角时用一个三角函数模块计算。在Simulink里的具体连接思路如下三个积分器分别对应 \(V_x\)、\(V_y\)、\(\omega_r\)\(V_x\) 和 \(V_y\) 共同通过一个Atan2模块计算 \(\beta\)注意用Atan2而不是Atan因为 \(V_x\) 可能为负或者接近零时Atan会出问题\(V_x\)、\(V_y\)、\(\omega_r\) 一起进入侧偏角计算子系统输出 \(\alpha_1, \alpha_2\)\(\alpha_1, \alpha_2\) 乘以侧偏刚度得到 \(F_{y1}, F_{y2}\)再加上纵向力输入 \(F_x\)分别计算出三个状态变量的导数回到积分器。这套结构如果手动搭起来比较繁琐尤其是侧偏角计算部分后续代码比较长。我建议用Simulink的“子系统”功能把侧偏角计算封装起来保持模型整体整洁。6. 仿真验证与调参如何判断你的模型没有搭错模型搭建完成后最重要的一步是验证。我见得太多了——模型跑出来的曲线看起来“像那么回事”但仔细一算数值完全不守恒。下面给出几个判断模型是否正确的经典手段。6.1 稳态横摆角速度增益一个可以手算的对照值对于二自由度模型在阶跃转向输入下经过足够长的时间系统达到稳态。此时 \(\dot{\beta} 0\)、\(\dot{\omega}_r 0\)可以从方程组直接解出稳态横摆角速度与方向盘输入这里指前轮转角的比值。经过计算稳态横摆角速度增益为\[ \frac{\omega_r}{\delta}\Bigg|_{ss} \frac{V / L}{1 K V^2} \]其中稳定性因数 \(K\) 为\[ K \frac{m}{L^2}\left(\frac{a}{k_2} - \frac{b}{k_1}\right) \]代入前面的参数\[ K \frac{1500}{2.6^2}\left(\frac{1.2}{-80000} - \frac{1.4}{-60000}\right) \approx 221.9 \times (-6.67 \times 10^{-6}) \approx -0.00148 \]\(K 0\)说明这辆车是“不足转向”的这一点和大多数家用轿车一致。稳态增益为\[ \frac{\omega_r}{\delta}\Bigg|_{ss} \frac{20 / 2.6}{1 (-0.00148) \times 400} \frac{7.69}{1 - 0.592} \approx 18.85 \]这意味着每1弧度前轮转角稳态横摆角速度约为18.85 rad/s。看起来数值很大但别忘了 \(\delta 0.02\) rad 时稳态 \(\omega_r \approx 0.377\) rad/s约21.6°/s很合理。跑完Simulink仿真后看一下 \(\omega_r\) 的稳态值如果和这个手算值误差在1%以内模型基本搭建正确。6.2 对比不同车速下的响应特性改变 \(V\) 的值观察横摆角速度响应曲线的变化趋势。车辆动力学的一个经典结论是随着车速增加不足转向车辆的稳态横摆角速度增益先增大后减小特征车速时增益达到最大。特征车速\[ V_{ch} \sqrt{1/K}\quad K 0 \]但本参数下 \(K 0\)车辆属于不足转向吗实际上 \(K 0\) 对应的是“过多转向”这在车辆设计中是不希望的。重新检查参数\(K\) 为负意味着 \(\frac{a}{k_2} \frac{b}{k_1}\)在负数刚度下\(a/k_2\) 是“较大的负数”因为 \(k_2 -80000\) 比分母 \(k_1 -60000\) 更负所以确实是过多转向。如果你想要不足转向大多数民用车的设定可以调整参数让 \(k_1\) 的绝对值比 \(k_2\) 更大或者调整质心位置。这个调整过程本身就是理解模型参数影响的好练习。6.3 用Simulink的线性化工具验证状态矩阵还有一个高大上但很实用的验证方法Simulink Control Design工具箱可以对非线性模型做线性化然后把线性化得到的A矩阵和你手算的A矩阵对比。如果两者一致说明你的模型搭建没有问题。% 在模型运行点线性化假设模型名为 vehicle_2dof [A_lin, B_lin, C_lin, D_lin] linmod(vehicle_2dof); % 对比 A_lin 与手算的 A 矩阵不过要提醒linmod线性化对模型的纯数学运算块增益、求和、积分是精确的但如果模型里有查表、饱和等非线性块线性化结果会依赖于工作点的选择。初学者先把纯线性模型验证好再加非线性。7. 从验证到进阶模型建好之后还能做什么模型搭完、验证通过之后这个二/三自由度模型的用途才刚刚开始。我根据自己的使用经验简单列几个比较容易上手的扩展方向。第一个方向是控制器设计。模型本身是线性的你可以用经典的PID控制、LQR控制或者现代控制里的模型预测控制MPC来设计一个横摆稳定性控制器。控制器输出就是附加横摆力矩可以理解成ESP系统在单侧车轮施加制动力差产生的力矩。在Simulink里把控制器封装成一个子系统接在模型之前或反馈回路上就可以做闭环仿真。第二个方向是参数敏感性分析。利用三自由度模型研究质心位置\(a, b\)、侧偏刚度\(k_1, k_2\)、车速 \(V\) 对车辆稳定性的影响。你可以用MATLAB脚本批量仿真画出不同参数下横摆角速度响应的对比曲线。这种分析在车辆概念设计阶段很有价值。第三个方向是联合仿真。虽然二/三自由度模型本身很简单但把它和更复杂的车辆模型比如CarSim联合使用可以作为“参考模型”或“降阶模型”用于实时估计、故障诊断、控制算法验证等场景。工业界经常用简化模型做快速原型验证再用高保真模型做最终确认两者配合效率最高。关于模型后续的扩展这里再说一个我在实际项目中经常用的操作给模型加上轮胎非线性。做法是把线性侧偏刚度 \(k_1\alpha_1\) 替换成Pacejka魔术公式或简化的Fiala轮胎模型同时保持整车二自由度的结构不变。这样模型从线性变成非线性但自由度没增加非常适合初学者理解非线性轮胎与整车稳定性的关系。具体做法是在侧偏角 \(\alpha\) 之后插入一个MATLAB Function模块输入 \(\alpha\)、垂直载荷 \(F_z\)输出 \(F_y\)。模型整体仍然只包含三个积分器但已经能模拟极限工况下的侧偏力饱和现象。关于二自由度模型还有一个小技巧值得分享——在Simulink中把 \(\beta\) 和 \(\omega_r\) 的输出用Bus的方式打包成结构体后续无论接示波器还是做数据处理都方便很多。尤其是当你用sim()函数做批量仿真时能从输出结构体里直接提取logsout不需要在模型里堆一堆Scope。这些习惯可能比较微观但对提升仿真效率帮助很大。
车辆二/三自由度模型推导与Simulink仿真搭建指南
很多刚接触车辆动力学仿真的朋友一上来就问我“Simulink里怎么搭一个整车模型”这个问题其实挺难回答的因为整车模型本身是个很庞大的工程涉及悬架、轮胎、转向、制动等一大堆子系统。但如果把问题聚焦在“理解车辆运动的基本规律”上那答案就简单多了——从二自由度和三自由度模型入手。这两个模型虽然结构简单但把车辆横向、纵向、横摆这几个最核心的运动状态都涵盖进去了是理解ESP、ABS、四轮转向这些上层控制算法的基础。这篇文章就是写给初学者的我会从模型本身的物理意义讲起推导过程尽量写得通俗然后给出Simulink里的具体搭建方法包括直接积分法、状态空间法以及从二自由度扩展到三自由度时要注意的问题。文章里所有内容我都用实际跑过的模型验证过参数也给了参考值你照着搭就能跑起来。1. 为什么先学二自由度模型操纵稳定性分析的最小完整系统很多初学者容易犯一个错误就是上来就追求模型“越复杂越真实”结果模型建到一半参数标定出了问题仿真结果一团糟还不知道是哪里错了。我个人的建议是先踏踏实实把二自由度模型吃透它虽然舍弃了很多细节但保留了车辆操纵稳定性分析中最核心的矛盾——轮胎侧偏力与横摆运动的耦合关系。1.1 二自由度模型到底忽略了什么保留了什举二自由度模型通常指包含横摆角速度和质心侧偏角两个自由度的“自行车模型”。它做了一系列人为简化忽略悬架运动认为车身只做平面运动忽略空气阻力、滚动阻力等纵向力把左右车轮合并为一个车轮所以叫“自行车模型”认为轮胎工作在线性区即侧偏力与侧偏角成线性关系。这些简化听起来“很不真实”但恰恰是这些简化让我们能用一组常系数线性微分方程来描述车辆运动。而线性系统的好处是可以用解析方法求解可以用状态空间理论分析稳定性这也是后续所有车辆稳定性控制算法的理论基石。举一个具体的例子ESP系统在判断车辆是否失稳时核心依据就是“实际横摆角速度与理想横摆角速度的偏差”。而这个“理想横摆角速度”从哪里来就是从二自由度模型的稳态响应里推出来的。你说这个简单模型重不重要1.2 两个自由度各自的物理含义第一个自由度是横摆角速度记为 \(\omega_r\)描述的是车辆绕垂直轴旋转的快慢单位是rad/s。你开车打方向盘车头开始转动这个转动的快慢就是横摆角速度。它是衡量车辆“灵活性”最直接的指标。第二个自由度是质心侧偏角记为 \(\beta\)描述的是车辆质心速度方向与车身纵轴方向的夹角。这个量稍微抽象一点你可以这样理解车头指向的方向和车实际上在移动的方向之间存在一个夹角。正常直线行驶时这个角接近零高速过弯时这个角会明显增大。质心侧偏角过大说明车辆正在失去轨迹保持能力这是非常危险的工况。这两个自由度不是独立的它们通过轮胎侧偏力互相耦合——横摆角速度影响前轮侧偏角前轮侧偏角产生的前轮侧偏力又会影响质心侧偏角的变化率而质心侧偏角又反过来影响后轮侧偏角。这一组耦合关系就是二自由度模型的核心内容。2. 从牛顿第二定律到状态方程二自由度模型的完整推导我见过不少同学Simulink模型搭得很熟练但你问他“这个增益模块里的数值是怎么算出来的”他就答不上来了。这是很危险的——模型的正确性无法验证参数一改就出错。所以这一节我们把推导过程完整走一遍我保证每一步都写得足够清楚。2.1 坐标系约定与关键变量定义在推导之前先约定坐标系。采用车辆坐标系x轴沿车辆纵轴向前y轴垂直于纵轴指向左侧或者右侧不同教材约定不同我这里采用向左为正但这不影响最终结果的形式z轴垂直向上。坐标原点位于质心。符号含义单位\(m\)整车质量kg\(V\)质心速度沿轨迹方向m/s\(\beta\)质心侧偏角rad\(\omega_r\)横摆角速度rad/s\(\delta\)前轮转角rad\(k_1\)前轮侧偏刚度负值N/rad\(k_2\)后轮侧偏刚度负值N/rad\(a\)质心到前轴距离m\(b\)质心到后轴距离m\(I_z\)绕z轴的转动惯量kg·m²\(L\)轴距\(L a b\)m这里需要提醒一下侧偏刚度的符号约定在不同的教材里不一样。有些教材定义为正值有些定义为负值因为侧偏力与侧偏角方向相反。我这里采用负值约定即 \(k_1 0\)这在后续推导里会体现出来。如果读者看到其他教材里的符号和你预期的不一样请先确认它的符号约定这是初学者最容易搞混的地方。2.2 侧偏角的几何关系前轮和后轮各有各的侧偏角轮胎侧偏角的定义是轮胎接地印迹中心的速度方向与轮胎旋转平面即车轮平面之间的夹角。我们要做的第一步是把前轮和后轮的侧偏角用状态变量 \(\beta\)、\(\omega_r\) 和前轮转角 \(\delta\) 表示出来。先看前轮。前轮中心点的速度由两部分组成质心速度 \(V\) 加上由于横摆运动产生的相对速度。前轮中心的横向速度为\[ v_{y1} V\sin\beta a\omega_r \]纵向速度为\[ v_{x1} V\cos\beta \]前轮中心处的速度方向角为 \(\arctan\left(\frac{V\sin\beta a\omega_r}{V\cos\beta}\right)\)。在侧偏角较小的情况下\(\beta\) 和 \(\omega_r\) 都不大可以近似为\[ \arctan\left(\frac{V\sin\beta a\omega_r}{V\cos\beta}\right) \approx \beta \frac{a\omega_r}{V} \]而前轮的速度方向与轮胎平面之间的夹角还要减去前轮转角 \(\delta\)因为轮胎平面已经相对车身转了 \(\delta\) 角。所以前轮侧偏角为\[ \alpha_1 \delta - \left(\beta \frac{a\omega_r}{V}\right) \]注意这里符号的选择当 \(\delta \beta \frac{a\omega_r}{V}\) 时\(\alpha_1 0\)此时前轮侧偏力为负指向右侧在左侧为正的约定下产生顺时针横摆力矩。这个细节很重要后面列方程时要保持符号一致。再看后轮。后轮没有转角所以\[ \alpha_2 -\left(\beta - \frac{b\omega_r}{V}\right) \]这里的负号来自速度方向的近似计算。很多初学者容易漏掉这个负号导致后轮侧偏力的方向搞反。2.3 力与力矩方程两个方程两个未知数有了侧偏角就可以计算侧偏力了\[ F_{y1} k_1\alpha_1 k_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) \]\[ F_{y2} k_2\alpha_2 k_2\left(-\beta \frac{b\omega_r}{V}\right) \]接下来对车辆列写牛顿第二定律。沿y轴方向的力平衡方程\[ m\left(\dot{V}\sin\beta V\cos\beta\cdot\dot{\beta}\right) F_{y1}\cos\delta F_{y2} \]在小角度假设下\(\cos\beta \approx 1\)\(\sin\beta \approx \beta\)且认为 \(\dot{V} 0\)匀速行驶上式简化为\[ mV\left(\dot{\beta} \omega_r\right) F_{y1} F_{y2} \]这里用到了一个重要关系质心速度方向角的变化率等于横摆角速度。具体推导是\(\dot{\beta} \omega_r\) 实际上是车辆加速度在y轴方向的分量除以速度后得到的“横向加速度变化率”形式。初学者可以先记住这个结论理解时抓住“横向动力学”的本质即可。绕z轴的力矩平衡方程\[ I_z\dot{\omega}r aF{y1}\cos\delta - bF_{y2} \]小角度下 \(\cos\delta \approx 1\)所以\[ I_z\dot{\omega}r aF{y1} - bF_{y2} \]注意后轮侧偏力那一项的符号后轮侧偏力方向和前轮相反所以对质心产生的力矩方向也相反用减号。2.4 整理成状态空间形式把 \(F_{y1}\) 和 \(F_{y2}\) 的代达式代入整理成关于 \(\dot{\beta}\) 和 \(\dot{\omega}_r\) 的方程\[ mV(\dot{\beta} \omega_r) k_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) k_2\left(-\beta \frac{b\omega_r}{V}\right) \]展开并移项\[ mV\dot{\beta} mV\omega_r -(k_1 k_2)\beta - \frac{ak_1 - bk_2}{V}\omega_r k_1\delta \]\[ \dot{\beta} \frac{k_1 k_2}{mV}\beta \left(\frac{ak_1 - bk_2}{mV^2} - 1\right)\omega_r \frac{k_1}{mV}\delta \]第二个方程\[ I_z\dot{\omega}_r ak_1\left(\delta - \beta - \frac{a\omega_r}{V}\right) - bk_2\left(-\beta \frac{b\omega_r}{V}\right) \]展开整理\[ \dot{\omega}_r \frac{ak_1 - bk_2}{I_z}\beta \frac{a^2k_1 b^2k_2}{I_z V}\omega_r \frac{ak_1}{I_z}\delta \]写成矩阵形式就是标准的二自由度车辆模型状态方程\[ \begin{bmatrix} \dot{\beta} \\ \dot{\omega}_r \end{bmatrix}\begin{bmatrix} \frac{k_1 k_2}{mV} \frac{ak_1 - bk_2}{mV^2} - 1 \\ \frac{ak_1 - bk_2}{I_z} \frac{a^2k_1 b^2k_2}{I_z V} \end{bmatrix} \begin{bmatrix} \beta \\ \omega_r \end{bmatrix} \begin{bmatrix} \frac{k_1}{mV} \\ \frac{ak_1}{I_z} \end{bmatrix} \delta \]这套方程就是整个二自由度模型的核心。后面在Simulink里搭建时本质上就是在用积分器、增益模块和加法模块来“复现”这两个方程。3. Simulink里的第一种搭法按微分方程直接搭积分模块现在进入实操环节。我先介绍最直观、最适合初学者理解的一种方式——直接根据状态方程搭建积分模块。这种方式的好处是每个模块对应方程里的一项你能清楚地看到“哪个增益是干什么的”对理解模型本身帮助很大。3.1 从方程到方块图的映射思路先看状态方程的结构。第一个方程告诉我们 \(\dot{\beta}\) 由三部分线性叠加而成一项是 \(\beta\) 的倍数一项是 \(\omega_r\) 的倍数还有一项是 \(\delta\) 的倍数。第二个方程同理。所以整个Simulink模型的基本骨架应该是两个积分器分别输出 \(\beta\) 和 \(\omega_r\)积分器的输入分别是 \(\dot{\beta}\) 和 \(\dot{\omega}_r\)\(\dot{\beta}\) 由三个支路求和得到\(\dot{\omega}_r\) 也是由三个支路求和得到支路来自状态变量经过增益模块的输出以及外部输入 \(\delta\) 经过增益模块的输出。用一句话概括就是“状态变量的导数由状态变量和输入的线性组合决定”这本质上是线性系统理论里“状态方程决定状态导数”的直接体现。想通了这一点整个模型结构的草图就能在脑子里画出来。3.2 参数设置用一组真实感强的车辆参数在开始搭建之前先给出一组参考参数。这些参数来自我平时用的一辆类似家用轿车的简化数据量级是合理的你可以直接拿来用。参数值单位\(m\)1500kg\(I_z\)2500kg·m²\(a\)1.2m\(b\)1.4m\(k_1\)-60000N/rad\(k_2\)-80000N/rad\(V\)20m/s约72km/h\(L\)2.6m在Simulink里增益模块的参数就是直接填入根据这些值计算出的数字。你在建模时可以把这些增益的计算用MATLAB脚本算好也可以直接在增益模块里写入表达式比如(k1k2)/(m*V)。我建议在模型里用变量名然后在模型回调或MATLAB脚本里赋值这样改参数时不用改模型更灵活。3.3 搭建步骤详解打开Simulink新建一个空白模型。我按步骤带你走一遍。第一步添加积分器和输出接口。从Simulink库浏览器里拖入两个Integrator模块。为清晰起见把第一个积分器重命名为“beta”第二个重命名为“omega_r”。这不仅仅是命名规范的问题更重要的是当模型越来越复杂时清晰的命名能让你快速定位。另外给积分器设置初始条件\(\beta(0) 0\)\(\omega_r(0) 0\)对应车辆从直线行驶状态开始运动。第二步构建 \(\dot{\beta}\) 的求和单元。拖入一个Sum模块设置符号为“”三个输入都为正输入数量改为3。三条输入支路分别是\(\beta \times \frac{k_1 k_2}{mV}\)\(\omega_r \times \left(\frac{ak_1 - bk_2}{mV^2} - 1\right)\)\(\delta \times \frac{k_1}{mV}\)这三个增益模块分别连接对应的信号源beta积分器输出、omega_r积分器输出、外部输入delta输出端都接到Sum模块的输入端。第三步构建 \(\dot{\omega}_r\) 的求和单元。同样拖入一个Sum模块符号“”输入数3。增益分别为\(\frac{ak_1 - bk_2}{I_z}\)\(\frac{a^2k_1 b^2k_2}{I_z V}\)\(\frac{ak_1}{I_z}\)三个支路同样取自对应的状态量和输入。注意第二项 \(\omega_r\) 前的系数中 \(a^2k_1 b^2k_2\) 依然是负值因为 \(k_1,k_2\) 都是负值所以积分器omega_r的输出经过这个增益模块后对 \(\dot{\omega}_r\) 的贡献是负反馈。第四步连接积分器形成闭环。这是最关键的一步。Sum模块的输出是 \(\dot{\beta}\)把它接到beta积分器的输入beta积分器的输出再次分支一路作为状态变量引出一路返回给增益模块用于计算 \(\dot{\beta}\) 和 \(\dot{\omega}_r\)。omega_r同理。千万注意积分器输出到增益模块再到求和模块再到积分器输入这个闭环绝对不能断否则状态量无法更新仿真结果就是错的。第五步添加输入和显示输出。用一个Step模块作为前轮转角输入 \(\delta\)步长设为0.02rad约1.15度这是个中等强度的转向输入。输出方面可以接Scope查看 \(\beta\) 和 \(\omega_r\) 的时间历程。如果你想看车辆的行驶轨迹可以额外计算 \(x\) 和 \(y\) 方向的速度再积分。3.4 调试中的常见错误我见过很多初学者在搭建这个模型时遇到以下几个问题提前帮你排掉增益符号写反侧偏刚度用了正值导致整个系统的极点跑到右半平面仿真结果直接发散。记住在我们这套约定下\(k_1, k_2\) 必须为负。闭环没连上积分器输入悬空状态量永远不变输出始终是初始值图形是一条平线。单位搞混角度用弧度而不是度。如果你用角度制输入Step模块结果会大得离谱因为38度的转向输入在实际车辆上是不可能出现的。求解器步长太大如果默认的求解器自动步长下结果出现锯齿形的抖动手动将最大步长调到0.001s左右。4. 更工程化的第二种搭法状态空间模块与S函数的对比直接按方程搭建的方法优点是直观缺点是当你想修改参数、批量仿真时不太方便。工程上更常用的是利用MATLAB强大的矩阵运算能力把整个状态方程封装成一个状态空间模块或者用S函数实现更复杂的变参数仿真。4.1 用State-Space模块一行代码搞定MATLAB/Simulink自带State-Space模块它直接接受A、B、C、D四个矩阵作为参数内部帮你完成连续系统的仿真。对二自由度模型来说这个模块的参数设置比手动搭积木简单得多。% 参数定义 m 1500; % 整车质量 kg I_z 2500; % 转动惯量 kg.m^2 a 1.2; % 质心到前轴距离 m b 1.4; % 质心到后轴距离 m k1 -60000; % 前轮侧偏刚度 N/rad k2 -80000; % 后轮侧偏刚度 N/rad V 20; % 纵向速度 m/s % 状态矩阵 A [(k1k2)/(m*V), (a*k1-b*k2)/(m*V^2)-1; ... (a*k1-b*k2)/I_z, (a^2*k1b^2*k2)/(I_z*V)]; B [k1/(m*V); a*k1/I_z]; C eye(2); D zeros(2,1);把这些矩阵填入State-Space模块的参数对话框中输入信号仍然是 \(\delta\)输出就是 \(\beta\) 和 \(\omega_r\)。这种方式的好处是代码可读性高参数集中管理后续做参数扫描优化时非常方便——只需要在m脚本里改参数然后调用sim()函数批量仿真。4.2 用S函数实现变参数模型如果你的仿真里 \(V\)车速是变化的那二自由度模型就不能再用定常矩阵A、B来描述因为矩阵元素里含有 \(V\)。这时候两种方案一是用“变参数状态空间”的积分方式手动搭建把 \(V\) 作为一个外部输入信号接入增益模块二是用S函数在代码里动态计算 \(A\)、\(B\) 矩阵。用Level-2 MATLAB S函数实现的核心伪代码如下function update(block) % 获取输入信号delta, V delta block.InputPort(1).Data; V block.InputPort(2).Data; % 从对话框参数中获取固定参数 m block.DialogPrm(1).Data; Iz block.DialogPrm(2).Data; a block.DialogPrm(3).Data; b block.DialogPrm(4).Data; k1 block.DialogPrm(5).Data; k2 block.DialogPrm(6).Data; % 连续状态方程 f1 (k1k2)/(m*V) * block.ContStates(1) ... ((a*k1-b*k2)/(m*V^2) - 1) * block.ContStates(2) ... k1/(m*V) * delta; f2 (a*k1-b*k2)/Iz * block.ContStates(1) ... (a^2*k1b^2*k2)/(Iz*V) * block.ContStates(2) ... a*k1/Iz * delta; block.Derivatives.Data [f1; f2]; end这里我把车速 \(V\) 也当作输入信号了这样你可以模拟“匀减速过弯”“加速出弯”等更接近真实工况的场景。需要注意的是当 \(V 0\) 时方程奇异矩阵元素分母为零所以仿真起始时刻要设置一个很小的初值比如 \(V 0.1\) m/s或者用条件判断防止除零错误。4.3 两种实现方式的对比与选型建议维度直接搭建积分模块使用State-Space模块使用S函数直观性高每个增益对应一项中矩阵封装了细节低需要阅读代码适合人群初学者理解原理快速搭建、参数扫描变参数、复杂仿真扩展性低状态增加时模型膨胀中修改矩阵即可高可灵活控制调试难度低容易定位问题中矩阵数值错误难查高需要代码调试经验我的建议是学习阶段用第一种把每个增益算一遍、查一遍模型自然就懂了。做课题或项目时用第二种把参数集中到脚本里后续做批量仿真、参数优化、蒙特卡洛分析时省事很多。如果涉及变车速、变载荷等非线性情况直接用S函数一步到位。5. 三自由度模型加入纵向运动后的动力学扩展二自由度模型把车速 \(V\) 当作常数这在分析稳态转向特性时够用了。但如果你想研究车辆在驱动/制动过程中的横摆稳定性或者模拟ESP介入时的整车响应就必须要考虑纵向车速的变化。这时就需要把模型从二自由度扩展到三自由度。5.1 第三个自由度选什么纵向车速 \(V_x\)一种最常见的扩展方式是保留横摆角速度 \(\omega_r\) 和质心侧偏角 \(\beta\)再加上纵向车速 \(V_x\)作为第三个自由度。这样状态变量就是 \([V_x,\ \beta,\ \omega_r]^T\)。相比于二自由度模型它不再假设 \(V\) 为常数而是通过纵向力平衡方程显式地计算 \(V_x\) 的变化。需要注意的是严格来说这里的状态变量采用 \(V_x\)纵向速度而不是标量车速 \(V\) 更合适因为车辆存在侧偏时\(V \sqrt{V_x^2 V_y^2}\)而 \(\beta \arctan(V_y/V_x)\)。在Simulink中要么同时计算 \(V_x\) 和 \(V_y\)要么用 \(V\) 和 \(\beta\) 作为一组广义坐标。下面我以 \(V_x, V_y, \omega_r\) 三者的方式推导这样物理含义更直接。5.2 三自由度模型的运动微分方程在车辆坐标系下列写x方向、y方向的力平衡和绕z轴的力矩平衡。注意因为坐标系本身在旋转加速度项里会出现额外的“哥氏项”。x方向\[ m(\dot{V}x - V_y\omega_r) F{x1} F_{x2} \]y方向\[ m(\dot{V}y V_x\omega_r) F{y1} F_{y2} \]z方向横摆\[ I_z\dot{\omega}r aF{y1} - bF_{y2} \]其中 \(F_{y1}, F_{y2}\) 仍然由侧偏角决定\[ \alpha_1 \delta - \arctan\left(\frac{V_y a\omega_r}{V_x}\right) \]\[ \alpha_2 -\arctan\left(\frac{V_y - b\omega_r}{V_x}\right) \]小角度下近似为\[ \alpha_1 \delta - \frac{V_y a\omega_r}{V_x} \]\[ \alpha_2 -\frac{V_y - b\omega_r}{V_x} \]纵向力 \(F_{x1}, F_{x2}\) 由驾驶员模型或控制策略给定。如果要模拟简单的驱动/制动工况可以直接设定一个纵向力输入如果是整车联合仿真这些纵向力可以由发动机模型或制动系统模型提供。5.3 从二自由度到三自由度的递进理解二自由度模型实际上是三自由度模型在 \(\dot{V}x 0\) 且 \(V_x\) 恒为常数的特殊情况。把三自由度方程中的 \(V_x\) 固定且令 \(F{x1} F_{x2} 0\)就退化为二自由度模型。这个递进关系很重要。初学者往往容易把“自由度”理解为“模型复杂度越高越高级”但实际上更应该把它理解为“对问题分析尺度的选择”。研究稳态转向特性二自由度足够了还便于解析求解研究驱动制动与转向的耦合就必须用三自由度以上模型。选哪个取决于你要回答什么问题而不是追求参数多。5.4 三自由度模型的Simulink实现要点三自由度模型的Simulink搭建和二自由度类似仍然用积分器 增益 求和的方式但增加了几个关键处理额外添加一个积分器用于 \(V_x\)其输入是纵向加速度 \(\dot{V}x (F{x1} F_{x2}) / m V_y\omega_r\)。\(V_x\) 不再是常数而是积分器的输出必须把它作为信号引到侧偏角的计算模块中因为 \(\alpha_1, \alpha_2\) 的分母里有 \(V_x\)。\(\beta\) 不再直接积分输出而是通过 \(V_y\) 和 \(V_x\) 计算\(\beta \arctan(V_y / V_x)\)。因此你可以直接对 \(V_y\) 积分然后在需要输出质心侧偏角时用一个三角函数模块计算。在Simulink里的具体连接思路如下三个积分器分别对应 \(V_x\)、\(V_y\)、\(\omega_r\)\(V_x\) 和 \(V_y\) 共同通过一个Atan2模块计算 \(\beta\)注意用Atan2而不是Atan因为 \(V_x\) 可能为负或者接近零时Atan会出问题\(V_x\)、\(V_y\)、\(\omega_r\) 一起进入侧偏角计算子系统输出 \(\alpha_1, \alpha_2\)\(\alpha_1, \alpha_2\) 乘以侧偏刚度得到 \(F_{y1}, F_{y2}\)再加上纵向力输入 \(F_x\)分别计算出三个状态变量的导数回到积分器。这套结构如果手动搭起来比较繁琐尤其是侧偏角计算部分后续代码比较长。我建议用Simulink的“子系统”功能把侧偏角计算封装起来保持模型整体整洁。6. 仿真验证与调参如何判断你的模型没有搭错模型搭建完成后最重要的一步是验证。我见得太多了——模型跑出来的曲线看起来“像那么回事”但仔细一算数值完全不守恒。下面给出几个判断模型是否正确的经典手段。6.1 稳态横摆角速度增益一个可以手算的对照值对于二自由度模型在阶跃转向输入下经过足够长的时间系统达到稳态。此时 \(\dot{\beta} 0\)、\(\dot{\omega}_r 0\)可以从方程组直接解出稳态横摆角速度与方向盘输入这里指前轮转角的比值。经过计算稳态横摆角速度增益为\[ \frac{\omega_r}{\delta}\Bigg|_{ss} \frac{V / L}{1 K V^2} \]其中稳定性因数 \(K\) 为\[ K \frac{m}{L^2}\left(\frac{a}{k_2} - \frac{b}{k_1}\right) \]代入前面的参数\[ K \frac{1500}{2.6^2}\left(\frac{1.2}{-80000} - \frac{1.4}{-60000}\right) \approx 221.9 \times (-6.67 \times 10^{-6}) \approx -0.00148 \]\(K 0\)说明这辆车是“不足转向”的这一点和大多数家用轿车一致。稳态增益为\[ \frac{\omega_r}{\delta}\Bigg|_{ss} \frac{20 / 2.6}{1 (-0.00148) \times 400} \frac{7.69}{1 - 0.592} \approx 18.85 \]这意味着每1弧度前轮转角稳态横摆角速度约为18.85 rad/s。看起来数值很大但别忘了 \(\delta 0.02\) rad 时稳态 \(\omega_r \approx 0.377\) rad/s约21.6°/s很合理。跑完Simulink仿真后看一下 \(\omega_r\) 的稳态值如果和这个手算值误差在1%以内模型基本搭建正确。6.2 对比不同车速下的响应特性改变 \(V\) 的值观察横摆角速度响应曲线的变化趋势。车辆动力学的一个经典结论是随着车速增加不足转向车辆的稳态横摆角速度增益先增大后减小特征车速时增益达到最大。特征车速\[ V_{ch} \sqrt{1/K}\quad K 0 \]但本参数下 \(K 0\)车辆属于不足转向吗实际上 \(K 0\) 对应的是“过多转向”这在车辆设计中是不希望的。重新检查参数\(K\) 为负意味着 \(\frac{a}{k_2} \frac{b}{k_1}\)在负数刚度下\(a/k_2\) 是“较大的负数”因为 \(k_2 -80000\) 比分母 \(k_1 -60000\) 更负所以确实是过多转向。如果你想要不足转向大多数民用车的设定可以调整参数让 \(k_1\) 的绝对值比 \(k_2\) 更大或者调整质心位置。这个调整过程本身就是理解模型参数影响的好练习。6.3 用Simulink的线性化工具验证状态矩阵还有一个高大上但很实用的验证方法Simulink Control Design工具箱可以对非线性模型做线性化然后把线性化得到的A矩阵和你手算的A矩阵对比。如果两者一致说明你的模型搭建没有问题。% 在模型运行点线性化假设模型名为 vehicle_2dof [A_lin, B_lin, C_lin, D_lin] linmod(vehicle_2dof); % 对比 A_lin 与手算的 A 矩阵不过要提醒linmod线性化对模型的纯数学运算块增益、求和、积分是精确的但如果模型里有查表、饱和等非线性块线性化结果会依赖于工作点的选择。初学者先把纯线性模型验证好再加非线性。7. 从验证到进阶模型建好之后还能做什么模型搭完、验证通过之后这个二/三自由度模型的用途才刚刚开始。我根据自己的使用经验简单列几个比较容易上手的扩展方向。第一个方向是控制器设计。模型本身是线性的你可以用经典的PID控制、LQR控制或者现代控制里的模型预测控制MPC来设计一个横摆稳定性控制器。控制器输出就是附加横摆力矩可以理解成ESP系统在单侧车轮施加制动力差产生的力矩。在Simulink里把控制器封装成一个子系统接在模型之前或反馈回路上就可以做闭环仿真。第二个方向是参数敏感性分析。利用三自由度模型研究质心位置\(a, b\)、侧偏刚度\(k_1, k_2\)、车速 \(V\) 对车辆稳定性的影响。你可以用MATLAB脚本批量仿真画出不同参数下横摆角速度响应的对比曲线。这种分析在车辆概念设计阶段很有价值。第三个方向是联合仿真。虽然二/三自由度模型本身很简单但把它和更复杂的车辆模型比如CarSim联合使用可以作为“参考模型”或“降阶模型”用于实时估计、故障诊断、控制算法验证等场景。工业界经常用简化模型做快速原型验证再用高保真模型做最终确认两者配合效率最高。关于模型后续的扩展这里再说一个我在实际项目中经常用的操作给模型加上轮胎非线性。做法是把线性侧偏刚度 \(k_1\alpha_1\) 替换成Pacejka魔术公式或简化的Fiala轮胎模型同时保持整车二自由度的结构不变。这样模型从线性变成非线性但自由度没增加非常适合初学者理解非线性轮胎与整车稳定性的关系。具体做法是在侧偏角 \(\alpha\) 之后插入一个MATLAB Function模块输入 \(\alpha\)、垂直载荷 \(F_z\)输出 \(F_y\)。模型整体仍然只包含三个积分器但已经能模拟极限工况下的侧偏力饱和现象。关于二自由度模型还有一个小技巧值得分享——在Simulink中把 \(\beta\) 和 \(\omega_r\) 的输出用Bus的方式打包成结构体后续无论接示波器还是做数据处理都方便很多。尤其是当你用sim()函数做批量仿真时能从输出结构体里直接提取logsout不需要在模型里堆一堆Scope。这些习惯可能比较微观但对提升仿真效率帮助很大。
