前几天有个刚入行的朋友来问我车辆二自由度模型的方程他推导起来没问题一打开 Simulink 就不知道模块该怎么摆。这个场景我太熟了几年前我第一次用 Simulink 搭二自由度车辆模型时也是对着状态方程发了半天呆。其实这类模型是所有整车操纵动力学仿真的起点二自由度模型解决“这辆车在给定方向盘输入下怎么横向运动”的核心问题三自由度模型再叠加一个侧倾运动更贴近真实车辆的动态表现。这篇文章就面向刚开始接触车辆仿真的初学者把核心方程、Simulink 建模思路、参数取值方法、常见错误一次讲清楚。看完你不仅能搭出一个能跑的模型还能用手算结果去校核仿真曲线知道自己搭出来的东西到底对不对。1. 二自由度模型到底在描述什么从整车到“自行车模型”的简化思路1.1 “自行车模型”是怎么来的为什么它能成立真实车辆在空间里有六个运动自由度沿 x、y、z 三个方向的平动加上绕这三根轴的转动。如果全都要建模悬架、转向、轮胎、空气动力学全都得进来复杂度瞬间爆炸。工程上研究操纵稳定性时重点关注的是车辆在水平面内的运动也就是侧向平动和横摆转动这两个自由度。于是就有了经典的二自由度模型也叫“自行车模型”。这个模型的简化思路非常直接把整车看成一个刚性车身前后轴各用一个等效车轮代替前轮转角直接作为系统输入忽略转向系统内部结构。这样一来整个车辆就简化成“前一个轮、后一个轮”的自行车结构。两个状态量分别是侧向速度 v沿 y 轴单位 m/s和横摆角速度 r绕 z 轴单位 rad/s。为什么这种简化能成立核心前提是在正常的稳态或准稳态转向工况下侧向加速度通常不超过 0.4g此时轮胎工作在线性侧偏区悬架对横向运动的影响可以折算到等效侧偏刚度里。换句话说二自由度模型不是“忽略”了悬架和轮胎的复杂性而是把它们的净效果合并成了几个等效参数。就像研究飞机飞行动力学时先用刚体方程抓住主要矛盾具体的气动弹性细节留到后续迭代再考虑。1.2 状态方程推导写出来才知道 Simulink 里该放什么模型方程看起来吓人其实推导过程很直接。先用 ISO 车辆坐标系约定x 轴向前y 轴向左z 轴向上横摆角速度 r 绕 z 轴逆时针为正。轮胎侧偏特性写成前轴侧偏力 F_yf -C_f · α_f后轴侧偏力 F_yr -C_r · α_r其中 C_f、C_r 是前后轴等效侧偏刚度正值单位 N/radα_f、α_r 是前后轮侧偏角。侧偏角根据运动学关系得到α_f δ - (v a·r) / uα_r -(v - b·r) / u这里 δ 是前轮转角a 是质心到前轴距离b 是质心到后轴距离u 是纵向车速。注意符号比较绕前轮侧偏角是车轮朝向与实际运动方向的夹角包含了前轮转角输入 δ后轮没有转向输入所以只有运动项。然后写牛顿方程。侧向力平衡m · (v_dot u·r) F_yf F_yr横摆力矩平衡I_z · r_dot a·F_yf - b·F_yr把侧偏力和侧偏角代进去整理成状态空间形式[v_dot] [ -(CfCr)/(m*u) -(a*Cf-b*Cr)/(m*u) - u ] [v] [ Cf/m ] [ ] [ ] [ ] [ ] δ [r_dot] [ -(a*Cf-b*Cr)/(Iz*u) -(a^2*Cfb^2*Cr)/(Iz*u)] [r] [ a*Cf/Iz ]看到这个矩阵你就知道 Simulink 里该放什么了两个积分器分别对 v_dot 和 r_dot 积分得到 v 和 r再把 v、r 通过对应增益反馈回去加上输入 δ 的前馈增益就是完整的闭环信号流。矩阵里每一项都有物理含义A(1,1) 是侧向速度对侧向加速度的阻尼项数值越大侧向速度收敛越快。A(1,2) 里的那个“-u”是横摆运动对侧向运动的耦合项这是二自由度模型最核心的耦合关系。车速越高横摆角速度对侧向速度的影响越强。A(2,1) 是侧向运动对横摆力矩的反馈它的大小直接决定了车辆的转向特性。A(2,2) 是横摆运动的阻尼项决定了横摆角速度自身的收敛速度。1.3 一个关键点恒定车速假设与线性轮胎假设二自由度模型里纵向车速 u 不是状态变量而是模型参数。仿真时要在 workspace 里把它设成一个固定值或者用阶跃信号在起步阶段快速拉到一个恒定值。这意味着这个模型不适用于急加速、急刹车工况它只适合分析匀速行驶时的转向响应。线性轮胎假设同样重要。侧偏刚度 C 是侧偏角 α 在小角度范围内的斜率一旦侧偏角超过约 5~8 度轮胎侧偏力进入饱和区线性模型的误差会越来越大。所以用二自由度模型做仿真前轮转角输入不要给太大0.05 rad约 2.9 度以内是比较安全的范围。要是你非给 0.5 rad那结果看起来“发散”其实是正常的因为模型假设已经失效了。2. Simulink 里搭二自由度模型的完整过程2.1 先理清信号流输入、状态、输出的关系搭建之前先画一条信号链前轮转角 δ 是输入状态是 v 和 r输出可以根据需要扩展比如侧向加速度 ay、侧偏角 α_f、α_r、转向半径 R 等。侧向加速度的计算公式是ay v_dot u·r注意这里用了 v_dot而不是直接用 v 乘以某个系数。因为侧向加速度包含两项车身坐标系下的横向加速度 v_dot以及因横摆运动产生的向心加速度 u·r。在 Simulink 里v_dot 就是第一个积分器的输入端信号直接把它引出来用就行不需要自己再算一遍微分这也避免了不必要的代数环问题。2.2 模块搭建五步走从积分器到示波器下面按步骤操作以一个全新的 Simulink 模型为例新建模型拖入两个 Integrator 模块分别命名为“v”和“r”。积分器的初始值都设为 0代表车辆初始时刻处于直线行驶状态。拖入两个 Sum 模块或 Add 模块分别计算 v_dot 和 r_dot。将 v 积分器的输出接入第一个 Sum 的一路输入r 积分器的输出接入第二个 Sum 的一路输入。拖入六个 Gain 模块按状态方程的反馈路径连接。具体是v 反馈到 v_dot 的增益是 A(1,1)r 反馈到 v_dot 的增益是 A(1,2)v 反馈到 r_dot 的增益是 A(2,1)r 反馈到 r_dot 的增益是 A(2,2)。输入 δ 到 v_dot 的增益是 B(1)输入 δ 到 r_dot 的增益是 B(2)。反馈路径连好之后把每个增益模块的数值填成 workspace 变量名比如 g11、g12、g21、g22、b1、b2方便后续批量改参数。这是实操中的一个重要技巧不要在 Gain 模块里直接填数字否则每次改车速都要手动改六个模块容易漏改出错。添加输出从 v_dot 和 r 信号引出接入一个 Sum 模块得到 ay再接 Scope 或 To Workspace 模块保存数据。侧偏角 α_f、α_r 也同理用 Gain 和 Sum 组合算出来方便后续与理论值对比。2.3 参数代入与仿真准备一个具体的数值算例光说抽象跑个具体算例才有感觉。取一组典型中级轿车参数参数数值单位整车质量 m1500kg质心到前轴距离 a1.2m质心到后轴距离 b1.4m前轴等效侧偏刚度 C_f80000N/rad后轴等效侧偏刚度 C_r100000N/rad横摆转动惯量 I_z2600kg·m²纵向车速 u30m/s先手工计算状态矩阵元素m·u 45000a·C_f - b·C_r 1.2 × 80000 - 1.4 × 100000 96000 - 140000 -44000于是A(1,1) -(80000 100000) / 45000 -4.0A(1,2) -(-44000) / 45000 - 30 0.978 - 30 -29.02A(2,1) -(-44000) / (2600 × 30) 0.564A(2,2) -(1.44 × 80000 1.96 × 100000) / (2600 × 30) -(115200 196000) / 78000 -3.99B(1) 80000 / 1500 53.33B(2) 1.2 × 80000 / 2600 36.92把六个数值填到对应 Gain 模块里建议直接填变量 g11 -4.0、g12 -29.02、g21 0.564、g22 -3.99、b1 53.33、b2 36.92。输入给阶跃信号幅值 δ 0.05 rad仿真时间 10 秒求解器用 ode45步长选 auto。跑完之后怎么验证结果用手算稳态横摆角速度增益来校验。稳定性因数K_us m·(b·C_r - a·C_f) / (L² · C_f · C_r)代入L 2.6 mb·C_r - a·C_f 140000 - 96000 44000所以K_us 1500 × 44000 / (2.6² × 80000 × 100000) ≈ 0.00122 s²/m稳态横摆角速度增益r_ss / δ u / (L · (1 K_us · u²)) 30 / (2.6 × (1 0.00122 × 900)) 5.50因此稳态 r_ss 5.50 × 0.05 0.275 rad/s。如果仿真曲线收敛到 0.275 附近说明模型搭对了如果差很多先回头查符号和单位。3. 从二自由度到三自由度加一个侧倾自由度到底改变了什么3.1 侧倾运动方程与耦合项二自由度模型忽略了车身侧倾。实际转弯时离心力让车身绕侧倾轴转动侧倾角 φ 会改变轮胎的载荷分布和车轮外倾角反过来又影响侧偏特性。三自由度模型就是在二自由度基础上把车身侧倾运动加进去状态向量变成 [v, r, φ_dot, φ]。车身侧倾力矩平衡方程可以写成I_x · φ_ddot C_φ · φ_dot K_φ · φ m_s · h · (v_dot u·r)这里 I_x 是簧载质量绕侧倾轴的转动惯量C_φ 是悬架侧倾阻尼K_φ 是悬架侧倾刚度m_s 是簧载质量h 是簧载质心到侧倾轴的垂直距离也就是侧倾力臂。等号右边的 m_s · h · ay 就是侧向惯性力产生的侧倾力矩。注意到 φ 方程里出现了 ay v_dot u·r而 v_dot 本身又来自二自由度模型的状态方程所以三个自由度之间存在耦合。在 Simulink 里实现时直接从 v 积分器的输入端把 v_dot 信号引出来做一次增益计算再接进 φ_ddot 的求和模块这个信号流是很自然的不会产生代数环。3.2 侧倾对等效侧偏刚度和侧偏角的影响增加侧倾自由度之后轮胎侧偏模型需要修正。主要有两个渠道第一是侧倾转向效应。车身侧倾时悬架导向机构会使车轮产生附加转角等效为对侧偏角的修正α_f δ - (v a·r) / u E_f · φα_r -(v - b·r) / u E_r · φE_f、E_r 是侧倾转向系数单位 rad/rad。前悬架一般设计成侧倾时前轮略微向转弯内侧偏转E_f 取正值后悬架的侧倾转向系数对车辆转向特性影响非常敏感很多高性能车就是靠调这个参数来平衡不足转向和过度转向。第二是侧倾外倾效应和载荷转移。车辆侧倾时外侧车轮载荷增加、内侧车轮载荷减少而轮胎侧偏刚度随垂向载荷的变化是非线性的载荷增加时侧偏刚度增大但增量递减载荷减少时侧偏刚度减小。综合效果是左右轮胎加起来的总等效侧偏刚度比没有载荷转移时小。这就是为什么侧倾会“软化”轮胎。对初学者来说先理解趋势就行三自由度模型里等效侧偏刚度实际上变小了所以相同输入下车辆响应会更“活”。3.3 二三自由度模型仿真结果的差异在哪同样一组车辆参数把三自由度模型搭出来跟二自由度对比跑一遍最直观的差异出现在瞬态阶段。假设 K_φ 60000 N·m/radm_s 1300 kgh 0.45 m稳态侧向加速度 ay 5 m/s²那么稳态侧倾力矩约 2925 N·m稳态侧倾角约 0.049 rad也就是 2.8 度。这个角度看着不大但乘以侧倾转向系数之后对侧偏角的修正可能达到 0.01~0.02 rad足以让横摆角速度响应出现可观测的变化。仿真曲线上你会看到三自由度模型的横摆角速度响应通常比二自由度模型有更大的超调量振荡衰减稍慢稳态值可能略有差异。侧倾角的响应曲线则表现为先快速增大、再逐渐收敛到稳态。如果侧倾刚度偏小侧倾角偏大车辆在瞬态阶段的“侧倾惯性”会让横摆响应显得更迟钝。这也是很多底盘调校工程师说的“侧倾响应跟横摆响应要匹配”的由来。4. 参数怎么取仿真结果可信度的关键4.1 车辆质量与几何参数整车质量 m 最好用实测值没有条件就用车辆公告参数或同级别车的公开数据。真正需要细分的是簧载质量 m_s 和非簧载质量 m_u因为侧倾方程里用到的 m_s 是簧载质量二自由度方程里的 m 是整车质量。初学阶段可以简单取 m_s 0.85~0.9 倍的整车质量。a、b 的确定比想象中容易。知道前后轴静载荷 F_zf、F_zr 就能算a F_zr · L / (m · g)b F_zf · L / (m · g)也就是说前轴载荷比例越大质心越靠前b 越大。乘用车前后轴载荷比例大约是 5545 到 6040。横摆转动惯量 I_z 最准确的获取方式是台架测量没有条件可以用经验公式估算。一个常用的范围是 I_z ≈ m · a · b 乘以一个系数系数通常在 0.8~1.2 之间也可以用 I_z ≈ m · L² / 12 到 m · L² / 3 粗估。对同一组参数做敏感性分析时会发现I_z 对横摆响应速度影响明显但对稳态增益影响不大所以初学时取个合理值就能跑通流程。下面给一组可直接使用的中级轿车参考参数参数数值单位整车质量 m1500kg簧载质量 m_s1300kg质心到前轴距离 a1.2m质心到后轴距离 b1.4m轴距 L2.6m横摆转动惯量 I_z2600kg·m²侧倾转动惯量 I_x600kg·m²侧倾力臂 h0.45m前悬架侧倾刚度 K_φf40000N·m/rad后悬架侧倾刚度 K_φr30000N·m/rad悬架侧倾阻尼 C_φ4000N·m·s/rad4.2 侧偏刚度取值的三种来源侧偏刚度是整个模型里最难精确取值、也是对结果影响最大的参数。它的定义是侧偏角在零度附近的侧偏力斜率轿车子午线轮胎的单胎侧偏刚度大约在 30000~150000 N/rad 之间载荷越大侧偏刚度越大胎压越高侧偏刚度越大。最好的数据来源是轮胎试验台的实测数据但对于个人学习场景通常拿不到。退而求其次有三种可行方案第一查同级别车型的公开论文或技术报告很多车辆动力学论文会列出所用参数找一辆轴距、质量相近的车把侧偏刚度按载荷比例折算一下也能用。第二用经验公式粗估。近似认为轮胎峰值侧偏力在侧偏角 8 度约 0.14 rad左右出现那么 C_α ≈ F_z / 0.14。比如前轴载荷 8500 N单胎侧偏刚度 ≈ 4250 / 0.14 ≈ 30000 N/rad前轴总侧偏刚度 C_f ≈ 60000 N/rad。这个办法很粗糙只适合建模起步阶段。第三用 Magic Formula 轮胎模型参数表取简化线性区间的斜率。以后想上非线性模型时同一套参数还能直接复用。这里要特别强调二自由度模型对 C_f 与 C_r 的比值非常敏感因为稳定性取决于 a·C_f 与 b·C_r 的相对大小。如果你的模型仿真结果对参数稍一变化就剧烈波动先检查一下是不是 C_f 和 C_r 给反了或者单胎刚度忘了乘以 2每个轴有两轮C_f 是左右轮刚度之和。4.3 参数扫描与敏感性分析模型搭好之后强烈建议做一个参数扫描这能帮你快速理解每个参数的作用。用 MATLAB 脚本批量调用 Simulink 模型把车速 u 从 20 变到 40 m/s每隔 5 m/s 跑一次看稳态横摆角速度增益和特征值怎么变。脚本结构很经典% 定义参数 m 1500; a 1.2; b 1.4; Iz 2600; Cf 80000; Cr 100000; L a b; u_list 20:5:40; for i 1:length(u_list) u u_list(i); % 计算状态矩阵 A11 -(Cf Cr) / (m * u); A12 -(a * Cf - b * Cr) / (m * u) - u; A21 -(a * Cf - b * Cr) / (Iz * u); A22 -(a^2 * Cf b^2 * Cr) / (Iz * u); B1 Cf / m; B2 a * Cf / Iz; % 赋值给模型中的变量 assignin(base, g11, A11); assignin(base, g12, A12); assignin(base, g21, A21); assignin(base, g22, A22); assignin(base, b1, B1); assignin(base, b2, B2); % 仿真 simOut sim(vehicle_2dof, StopTime, 10); % 提取横摆角速度稳态值 r_ss(i) simOut.r.Data(end); end % 与理论值对比 K_us m * (b * Cr - a * Cf) / (L^2 * Cf * Cr); r_theory u_list ./ (L * (1 K_us * u_list.^2)) * 0.05; plot(u_list, r_ss, o, u_list, r_theory, -);跑完你会发现不足转向车辆的稳态增益随车速先上升后下降在特征车速处达到峰值同一次扫描还能画出根轨迹看特征值虚部怎么变化理解阻尼特性随车速变差的趋势。5. 初学者的常见错误与排查实录5.1 常见错误速查表这些年见过不少新人在这个模型上翻车问题高度集中直接列成表错误现象可能原因解决方法仿真结果直接发散数值爆炸状态方程符号搞反单位不统一步长过大重新推导矩阵核对每一项符号统一用 m/s、rad、N/rad稳态值对不上理论计算B 矩阵算错侧偏刚度单位写成 N/deg用手算稳态增益做静态校核输出曲线高频抖动存在代数环信号采样率过低加 Memory 模块或改用 State-Space 模块仿真特别慢系统呈现刚性特征步长设置不合理改用 ode15s 等刚性求解器报 Algebraic Loop 警告信号在同一个步长内形成无延迟循环检查是否把积分器输出以外的信号用于反馈低速正常、高速结果完全不对漏了 A(1,2) 里的 -u 耦合项检查状态矩阵 A(1,2) 是否包含 -u输入转角信号看着正常结果却饱和角度单位混用把角度当弧度输入0.05 弧度约 2.9 度0.05 度约 0.00087 弧度5.2 代数环与数值求解问题代数环是 Simulink 初学者最容易撞上的警告。它的本质是某个信号的输出直接或间接反馈到自己的输入而且中间没有任何延迟模块Simulink 不得不在每个仿真步长里迭代求解这个环。二自由度模型里如果不小心把 v_dot 信号先算出来、再用它去修正输入就有可能形成代数环。解决办法有两个一是在反馈路径上插入 Memory 模块相当于引入一个步长的延迟打破迭代环二是改用 State-Space 模块把 A、B、C、D 矩阵一次性填进去内部处理方式不会有代数环。但要注意Memory 模块会引入一个采样周期的延迟如果仿真步长较大这个延迟会影响瞬态精度。更干净的思路是重构信号流确保反馈信号都是从积分器输出端引出来的不要从积分器输入端反推。求解器选择上也有些讲究。二自由度模型一般是光滑的常微分方程ode45 默认配置就能跑得很好。如果加了侧倾方程之后出现高频振荡、仿真步长骤减可以切换到 ode15s 这类刚性求解器速度会有明显改善。5.3 结果“不稳定”时先查模型还是先查参数仿真曲线发散时第一反应不要急着调参数先做三件事第一静态校核。用手算公式算一遍稳态增益和特征值看理论值是否稳定。特征值实部为负是稳定的必要条件如果实部为正说明模型参数本身就不稳定不是数值求解的锅。第二单位校核。把模型里所有模块的增益列一张表逐一确认单位。最经典的错误是把车速写成 30 km/h 而不是 30 m/s或者把侧偏刚度写成单胎值忘了乘 2。第三信号方向校核。把反馈增益的正负号和方程里的符号对一遍。我自己的经验是90% 的“莫名发散”最后都出在单位换算和漏了负号上真正模型原理有问题的情况反而不多。6. 仿真结果怎么解读从开环响应看车辆稳定性6.1 阶跃转向输入下的瞬态与稳态响应给模型一个阶跃前轮转角输入观察横摆角速度 r 的时域响应这是整车操纵稳定性分析最经典的测试。响应曲线通常包含三个阶段上升段、可能的超调振荡段、稳态段。继续用前面的参数算特征值。A 矩阵的迹是 -7.99行列式约 32.33特征值约 -3.995 ± 4.046j。这说明系统是一个稳定的欠阻尼二阶系统无阻尼自然频率约 5.7 rad/s约 0.9 Hz阻尼比约 0.70。按照二阶系统公式估算超调量约 4.6%上升时间约 0.3 秒稳态值 0.275 rad/s。把这些理论值和仿真曲线对照看能很直观地理解“欠阻尼”意味着什么。如果换一组参数让阻尼比变小比如增大车速 u特征值虚部会比实部下降得更慢振荡就会更明显。这也是为什么高速变线时车辆动态感觉更“晃”的数学解释。6.2 特征值与稳定性阈值的可视化分析把车速当成变量每个车速都算一遍特征值画在复平面上就是根轨迹。对不足转向车辆特征值实部始终为负但随着车速增加特征值逐渐向虚轴靠近说明稳定性裕度在下降。对过度转向车辆存在一个临界车速超过之后特征值实部变正系统失稳。临界车速的公式是u_crit sqrt(C_f · C_r · L² / (m · (b·C_r - a·C_f)))代入前面的参数b·C_r - a·C_f 44000分母 1500 × 44000所以临界车速超过 100 m/s说明这台车在正常车速范围内都是不足转向特性不会出现固有失稳。反过来如果某辆车的参数使 b·C_r - a·C_f 为负临界车速可能在 30~40 m/s 左右那这辆车高速行驶就非常危险这也是汽车工程里严格杜绝过度转向特性的原因。6.3 从增益差异理解不足转向与过度转向判定转向特性的另一个直观方法是看稳态横摆角速度增益随车速的变化趋势不足转向K_us 0增益先随车速增大到达特征车速后开始下降曲线有个明显的峰值。中性转向K_us 0增益与车速成正比直线上升。过度转向K_us 0增益随车速增大而加速上升在临界车速处趋于无穷。用当前参数算一下 20、30、40 m/s 三个车速的增益车速 m/s稳态横摆角速度增益趋势205.17上升段305.50接近峰值405.21下降段峰值出现在特征车速 u_ch sqrt(1 / K_us) ≈ 28.6 m/s。这就是为什么有些车在 100 km/h 左右变道时感觉特别灵活再往上反而感觉“迟钝”了一点的原因。这种用一个小模型就能解释的真实驾驶感受也是我觉得二自由度模型特别值得花时间吃透的地方。最后说一点个人体会。我建议初学者不要嫌二自由度模型太简单先在纸上把状态方程一步步推导一遍再用积分器和增益模块搭出来最后用 State-Space 模块交叉验证一遍。这个过程走下来你对状态空间表达、信号流方向、参数灵敏度、稳定性判据的理解都会比直接拖一个大而全的模型扎实得多。后续要扩展方向也很清楚把线性轮胎换成魔术公式进入非线性区把恒定车速假设放开做纵侧向耦合或者接 Carsim 做联合仿真。但不管往哪个方向走二自由度模型打底的理解能力都会让你少踩很多坑。
前几天有个刚入行的朋友来问我车辆二自由度模型的方程他推导起来没问题一打开 Simulink 就不知道模块该怎么摆。这个场景我太熟了几年前我第一次用 Simulink 搭二自由度车辆模型时也是对着状态方程发了半天呆。其实这类模型是所有整车操纵动力学仿真的起点二自由度模型解决“这辆车在给定方向盘输入下怎么横向运动”的核心问题三自由度模型再叠加一个侧倾运动更贴近真实车辆的动态表现。这篇文章就面向刚开始接触车辆仿真的初学者把核心方程、Simulink 建模思路、参数取值方法、常见错误一次讲清楚。看完你不仅能搭出一个能跑的模型还能用手算结果去校核仿真曲线知道自己搭出来的东西到底对不对。1. 二自由度模型到底在描述什么从整车到“自行车模型”的简化思路1.1 “自行车模型”是怎么来的为什么它能成立真实车辆在空间里有六个运动自由度沿 x、y、z 三个方向的平动加上绕这三根轴的转动。如果全都要建模悬架、转向、轮胎、空气动力学全都得进来复杂度瞬间爆炸。工程上研究操纵稳定性时重点关注的是车辆在水平面内的运动也就是侧向平动和横摆转动这两个自由度。于是就有了经典的二自由度模型也叫“自行车模型”。这个模型的简化思路非常直接把整车看成一个刚性车身前后轴各用一个等效车轮代替前轮转角直接作为系统输入忽略转向系统内部结构。这样一来整个车辆就简化成“前一个轮、后一个轮”的自行车结构。两个状态量分别是侧向速度 v沿 y 轴单位 m/s和横摆角速度 r绕 z 轴单位 rad/s。为什么这种简化能成立核心前提是在正常的稳态或准稳态转向工况下侧向加速度通常不超过 0.4g此时轮胎工作在线性侧偏区悬架对横向运动的影响可以折算到等效侧偏刚度里。换句话说二自由度模型不是“忽略”了悬架和轮胎的复杂性而是把它们的净效果合并成了几个等效参数。就像研究飞机飞行动力学时先用刚体方程抓住主要矛盾具体的气动弹性细节留到后续迭代再考虑。1.2 状态方程推导写出来才知道 Simulink 里该放什么模型方程看起来吓人其实推导过程很直接。先用 ISO 车辆坐标系约定x 轴向前y 轴向左z 轴向上横摆角速度 r 绕 z 轴逆时针为正。轮胎侧偏特性写成前轴侧偏力 F_yf -C_f · α_f后轴侧偏力 F_yr -C_r · α_r其中 C_f、C_r 是前后轴等效侧偏刚度正值单位 N/radα_f、α_r 是前后轮侧偏角。侧偏角根据运动学关系得到α_f δ - (v a·r) / uα_r -(v - b·r) / u这里 δ 是前轮转角a 是质心到前轴距离b 是质心到后轴距离u 是纵向车速。注意符号比较绕前轮侧偏角是车轮朝向与实际运动方向的夹角包含了前轮转角输入 δ后轮没有转向输入所以只有运动项。然后写牛顿方程。侧向力平衡m · (v_dot u·r) F_yf F_yr横摆力矩平衡I_z · r_dot a·F_yf - b·F_yr把侧偏力和侧偏角代进去整理成状态空间形式[v_dot] [ -(CfCr)/(m*u) -(a*Cf-b*Cr)/(m*u) - u ] [v] [ Cf/m ] [ ] [ ] [ ] [ ] δ [r_dot] [ -(a*Cf-b*Cr)/(Iz*u) -(a^2*Cfb^2*Cr)/(Iz*u)] [r] [ a*Cf/Iz ]看到这个矩阵你就知道 Simulink 里该放什么了两个积分器分别对 v_dot 和 r_dot 积分得到 v 和 r再把 v、r 通过对应增益反馈回去加上输入 δ 的前馈增益就是完整的闭环信号流。矩阵里每一项都有物理含义A(1,1) 是侧向速度对侧向加速度的阻尼项数值越大侧向速度收敛越快。A(1,2) 里的那个“-u”是横摆运动对侧向运动的耦合项这是二自由度模型最核心的耦合关系。车速越高横摆角速度对侧向速度的影响越强。A(2,1) 是侧向运动对横摆力矩的反馈它的大小直接决定了车辆的转向特性。A(2,2) 是横摆运动的阻尼项决定了横摆角速度自身的收敛速度。1.3 一个关键点恒定车速假设与线性轮胎假设二自由度模型里纵向车速 u 不是状态变量而是模型参数。仿真时要在 workspace 里把它设成一个固定值或者用阶跃信号在起步阶段快速拉到一个恒定值。这意味着这个模型不适用于急加速、急刹车工况它只适合分析匀速行驶时的转向响应。线性轮胎假设同样重要。侧偏刚度 C 是侧偏角 α 在小角度范围内的斜率一旦侧偏角超过约 5~8 度轮胎侧偏力进入饱和区线性模型的误差会越来越大。所以用二自由度模型做仿真前轮转角输入不要给太大0.05 rad约 2.9 度以内是比较安全的范围。要是你非给 0.5 rad那结果看起来“发散”其实是正常的因为模型假设已经失效了。2. Simulink 里搭二自由度模型的完整过程2.1 先理清信号流输入、状态、输出的关系搭建之前先画一条信号链前轮转角 δ 是输入状态是 v 和 r输出可以根据需要扩展比如侧向加速度 ay、侧偏角 α_f、α_r、转向半径 R 等。侧向加速度的计算公式是ay v_dot u·r注意这里用了 v_dot而不是直接用 v 乘以某个系数。因为侧向加速度包含两项车身坐标系下的横向加速度 v_dot以及因横摆运动产生的向心加速度 u·r。在 Simulink 里v_dot 就是第一个积分器的输入端信号直接把它引出来用就行不需要自己再算一遍微分这也避免了不必要的代数环问题。2.2 模块搭建五步走从积分器到示波器下面按步骤操作以一个全新的 Simulink 模型为例新建模型拖入两个 Integrator 模块分别命名为“v”和“r”。积分器的初始值都设为 0代表车辆初始时刻处于直线行驶状态。拖入两个 Sum 模块或 Add 模块分别计算 v_dot 和 r_dot。将 v 积分器的输出接入第一个 Sum 的一路输入r 积分器的输出接入第二个 Sum 的一路输入。拖入六个 Gain 模块按状态方程的反馈路径连接。具体是v 反馈到 v_dot 的增益是 A(1,1)r 反馈到 v_dot 的增益是 A(1,2)v 反馈到 r_dot 的增益是 A(2,1)r 反馈到 r_dot 的增益是 A(2,2)。输入 δ 到 v_dot 的增益是 B(1)输入 δ 到 r_dot 的增益是 B(2)。反馈路径连好之后把每个增益模块的数值填成 workspace 变量名比如 g11、g12、g21、g22、b1、b2方便后续批量改参数。这是实操中的一个重要技巧不要在 Gain 模块里直接填数字否则每次改车速都要手动改六个模块容易漏改出错。添加输出从 v_dot 和 r 信号引出接入一个 Sum 模块得到 ay再接 Scope 或 To Workspace 模块保存数据。侧偏角 α_f、α_r 也同理用 Gain 和 Sum 组合算出来方便后续与理论值对比。2.3 参数代入与仿真准备一个具体的数值算例光说抽象跑个具体算例才有感觉。取一组典型中级轿车参数参数数值单位整车质量 m1500kg质心到前轴距离 a1.2m质心到后轴距离 b1.4m前轴等效侧偏刚度 C_f80000N/rad后轴等效侧偏刚度 C_r100000N/rad横摆转动惯量 I_z2600kg·m²纵向车速 u30m/s先手工计算状态矩阵元素m·u 45000a·C_f - b·C_r 1.2 × 80000 - 1.4 × 100000 96000 - 140000 -44000于是A(1,1) -(80000 100000) / 45000 -4.0A(1,2) -(-44000) / 45000 - 30 0.978 - 30 -29.02A(2,1) -(-44000) / (2600 × 30) 0.564A(2,2) -(1.44 × 80000 1.96 × 100000) / (2600 × 30) -(115200 196000) / 78000 -3.99B(1) 80000 / 1500 53.33B(2) 1.2 × 80000 / 2600 36.92把六个数值填到对应 Gain 模块里建议直接填变量 g11 -4.0、g12 -29.02、g21 0.564、g22 -3.99、b1 53.33、b2 36.92。输入给阶跃信号幅值 δ 0.05 rad仿真时间 10 秒求解器用 ode45步长选 auto。跑完之后怎么验证结果用手算稳态横摆角速度增益来校验。稳定性因数K_us m·(b·C_r - a·C_f) / (L² · C_f · C_r)代入L 2.6 mb·C_r - a·C_f 140000 - 96000 44000所以K_us 1500 × 44000 / (2.6² × 80000 × 100000) ≈ 0.00122 s²/m稳态横摆角速度增益r_ss / δ u / (L · (1 K_us · u²)) 30 / (2.6 × (1 0.00122 × 900)) 5.50因此稳态 r_ss 5.50 × 0.05 0.275 rad/s。如果仿真曲线收敛到 0.275 附近说明模型搭对了如果差很多先回头查符号和单位。3. 从二自由度到三自由度加一个侧倾自由度到底改变了什么3.1 侧倾运动方程与耦合项二自由度模型忽略了车身侧倾。实际转弯时离心力让车身绕侧倾轴转动侧倾角 φ 会改变轮胎的载荷分布和车轮外倾角反过来又影响侧偏特性。三自由度模型就是在二自由度基础上把车身侧倾运动加进去状态向量变成 [v, r, φ_dot, φ]。车身侧倾力矩平衡方程可以写成I_x · φ_ddot C_φ · φ_dot K_φ · φ m_s · h · (v_dot u·r)这里 I_x 是簧载质量绕侧倾轴的转动惯量C_φ 是悬架侧倾阻尼K_φ 是悬架侧倾刚度m_s 是簧载质量h 是簧载质心到侧倾轴的垂直距离也就是侧倾力臂。等号右边的 m_s · h · ay 就是侧向惯性力产生的侧倾力矩。注意到 φ 方程里出现了 ay v_dot u·r而 v_dot 本身又来自二自由度模型的状态方程所以三个自由度之间存在耦合。在 Simulink 里实现时直接从 v 积分器的输入端把 v_dot 信号引出来做一次增益计算再接进 φ_ddot 的求和模块这个信号流是很自然的不会产生代数环。3.2 侧倾对等效侧偏刚度和侧偏角的影响增加侧倾自由度之后轮胎侧偏模型需要修正。主要有两个渠道第一是侧倾转向效应。车身侧倾时悬架导向机构会使车轮产生附加转角等效为对侧偏角的修正α_f δ - (v a·r) / u E_f · φα_r -(v - b·r) / u E_r · φE_f、E_r 是侧倾转向系数单位 rad/rad。前悬架一般设计成侧倾时前轮略微向转弯内侧偏转E_f 取正值后悬架的侧倾转向系数对车辆转向特性影响非常敏感很多高性能车就是靠调这个参数来平衡不足转向和过度转向。第二是侧倾外倾效应和载荷转移。车辆侧倾时外侧车轮载荷增加、内侧车轮载荷减少而轮胎侧偏刚度随垂向载荷的变化是非线性的载荷增加时侧偏刚度增大但增量递减载荷减少时侧偏刚度减小。综合效果是左右轮胎加起来的总等效侧偏刚度比没有载荷转移时小。这就是为什么侧倾会“软化”轮胎。对初学者来说先理解趋势就行三自由度模型里等效侧偏刚度实际上变小了所以相同输入下车辆响应会更“活”。3.3 二三自由度模型仿真结果的差异在哪同样一组车辆参数把三自由度模型搭出来跟二自由度对比跑一遍最直观的差异出现在瞬态阶段。假设 K_φ 60000 N·m/radm_s 1300 kgh 0.45 m稳态侧向加速度 ay 5 m/s²那么稳态侧倾力矩约 2925 N·m稳态侧倾角约 0.049 rad也就是 2.8 度。这个角度看着不大但乘以侧倾转向系数之后对侧偏角的修正可能达到 0.01~0.02 rad足以让横摆角速度响应出现可观测的变化。仿真曲线上你会看到三自由度模型的横摆角速度响应通常比二自由度模型有更大的超调量振荡衰减稍慢稳态值可能略有差异。侧倾角的响应曲线则表现为先快速增大、再逐渐收敛到稳态。如果侧倾刚度偏小侧倾角偏大车辆在瞬态阶段的“侧倾惯性”会让横摆响应显得更迟钝。这也是很多底盘调校工程师说的“侧倾响应跟横摆响应要匹配”的由来。4. 参数怎么取仿真结果可信度的关键4.1 车辆质量与几何参数整车质量 m 最好用实测值没有条件就用车辆公告参数或同级别车的公开数据。真正需要细分的是簧载质量 m_s 和非簧载质量 m_u因为侧倾方程里用到的 m_s 是簧载质量二自由度方程里的 m 是整车质量。初学阶段可以简单取 m_s 0.85~0.9 倍的整车质量。a、b 的确定比想象中容易。知道前后轴静载荷 F_zf、F_zr 就能算a F_zr · L / (m · g)b F_zf · L / (m · g)也就是说前轴载荷比例越大质心越靠前b 越大。乘用车前后轴载荷比例大约是 5545 到 6040。横摆转动惯量 I_z 最准确的获取方式是台架测量没有条件可以用经验公式估算。一个常用的范围是 I_z ≈ m · a · b 乘以一个系数系数通常在 0.8~1.2 之间也可以用 I_z ≈ m · L² / 12 到 m · L² / 3 粗估。对同一组参数做敏感性分析时会发现I_z 对横摆响应速度影响明显但对稳态增益影响不大所以初学时取个合理值就能跑通流程。下面给一组可直接使用的中级轿车参考参数参数数值单位整车质量 m1500kg簧载质量 m_s1300kg质心到前轴距离 a1.2m质心到后轴距离 b1.4m轴距 L2.6m横摆转动惯量 I_z2600kg·m²侧倾转动惯量 I_x600kg·m²侧倾力臂 h0.45m前悬架侧倾刚度 K_φf40000N·m/rad后悬架侧倾刚度 K_φr30000N·m/rad悬架侧倾阻尼 C_φ4000N·m·s/rad4.2 侧偏刚度取值的三种来源侧偏刚度是整个模型里最难精确取值、也是对结果影响最大的参数。它的定义是侧偏角在零度附近的侧偏力斜率轿车子午线轮胎的单胎侧偏刚度大约在 30000~150000 N/rad 之间载荷越大侧偏刚度越大胎压越高侧偏刚度越大。最好的数据来源是轮胎试验台的实测数据但对于个人学习场景通常拿不到。退而求其次有三种可行方案第一查同级别车型的公开论文或技术报告很多车辆动力学论文会列出所用参数找一辆轴距、质量相近的车把侧偏刚度按载荷比例折算一下也能用。第二用经验公式粗估。近似认为轮胎峰值侧偏力在侧偏角 8 度约 0.14 rad左右出现那么 C_α ≈ F_z / 0.14。比如前轴载荷 8500 N单胎侧偏刚度 ≈ 4250 / 0.14 ≈ 30000 N/rad前轴总侧偏刚度 C_f ≈ 60000 N/rad。这个办法很粗糙只适合建模起步阶段。第三用 Magic Formula 轮胎模型参数表取简化线性区间的斜率。以后想上非线性模型时同一套参数还能直接复用。这里要特别强调二自由度模型对 C_f 与 C_r 的比值非常敏感因为稳定性取决于 a·C_f 与 b·C_r 的相对大小。如果你的模型仿真结果对参数稍一变化就剧烈波动先检查一下是不是 C_f 和 C_r 给反了或者单胎刚度忘了乘以 2每个轴有两轮C_f 是左右轮刚度之和。4.3 参数扫描与敏感性分析模型搭好之后强烈建议做一个参数扫描这能帮你快速理解每个参数的作用。用 MATLAB 脚本批量调用 Simulink 模型把车速 u 从 20 变到 40 m/s每隔 5 m/s 跑一次看稳态横摆角速度增益和特征值怎么变。脚本结构很经典% 定义参数 m 1500; a 1.2; b 1.4; Iz 2600; Cf 80000; Cr 100000; L a b; u_list 20:5:40; for i 1:length(u_list) u u_list(i); % 计算状态矩阵 A11 -(Cf Cr) / (m * u); A12 -(a * Cf - b * Cr) / (m * u) - u; A21 -(a * Cf - b * Cr) / (Iz * u); A22 -(a^2 * Cf b^2 * Cr) / (Iz * u); B1 Cf / m; B2 a * Cf / Iz; % 赋值给模型中的变量 assignin(base, g11, A11); assignin(base, g12, A12); assignin(base, g21, A21); assignin(base, g22, A22); assignin(base, b1, B1); assignin(base, b2, B2); % 仿真 simOut sim(vehicle_2dof, StopTime, 10); % 提取横摆角速度稳态值 r_ss(i) simOut.r.Data(end); end % 与理论值对比 K_us m * (b * Cr - a * Cf) / (L^2 * Cf * Cr); r_theory u_list ./ (L * (1 K_us * u_list.^2)) * 0.05; plot(u_list, r_ss, o, u_list, r_theory, -);跑完你会发现不足转向车辆的稳态增益随车速先上升后下降在特征车速处达到峰值同一次扫描还能画出根轨迹看特征值虚部怎么变化理解阻尼特性随车速变差的趋势。5. 初学者的常见错误与排查实录5.1 常见错误速查表这些年见过不少新人在这个模型上翻车问题高度集中直接列成表错误现象可能原因解决方法仿真结果直接发散数值爆炸状态方程符号搞反单位不统一步长过大重新推导矩阵核对每一项符号统一用 m/s、rad、N/rad稳态值对不上理论计算B 矩阵算错侧偏刚度单位写成 N/deg用手算稳态增益做静态校核输出曲线高频抖动存在代数环信号采样率过低加 Memory 模块或改用 State-Space 模块仿真特别慢系统呈现刚性特征步长设置不合理改用 ode15s 等刚性求解器报 Algebraic Loop 警告信号在同一个步长内形成无延迟循环检查是否把积分器输出以外的信号用于反馈低速正常、高速结果完全不对漏了 A(1,2) 里的 -u 耦合项检查状态矩阵 A(1,2) 是否包含 -u输入转角信号看着正常结果却饱和角度单位混用把角度当弧度输入0.05 弧度约 2.9 度0.05 度约 0.00087 弧度5.2 代数环与数值求解问题代数环是 Simulink 初学者最容易撞上的警告。它的本质是某个信号的输出直接或间接反馈到自己的输入而且中间没有任何延迟模块Simulink 不得不在每个仿真步长里迭代求解这个环。二自由度模型里如果不小心把 v_dot 信号先算出来、再用它去修正输入就有可能形成代数环。解决办法有两个一是在反馈路径上插入 Memory 模块相当于引入一个步长的延迟打破迭代环二是改用 State-Space 模块把 A、B、C、D 矩阵一次性填进去内部处理方式不会有代数环。但要注意Memory 模块会引入一个采样周期的延迟如果仿真步长较大这个延迟会影响瞬态精度。更干净的思路是重构信号流确保反馈信号都是从积分器输出端引出来的不要从积分器输入端反推。求解器选择上也有些讲究。二自由度模型一般是光滑的常微分方程ode45 默认配置就能跑得很好。如果加了侧倾方程之后出现高频振荡、仿真步长骤减可以切换到 ode15s 这类刚性求解器速度会有明显改善。5.3 结果“不稳定”时先查模型还是先查参数仿真曲线发散时第一反应不要急着调参数先做三件事第一静态校核。用手算公式算一遍稳态增益和特征值看理论值是否稳定。特征值实部为负是稳定的必要条件如果实部为正说明模型参数本身就不稳定不是数值求解的锅。第二单位校核。把模型里所有模块的增益列一张表逐一确认单位。最经典的错误是把车速写成 30 km/h 而不是 30 m/s或者把侧偏刚度写成单胎值忘了乘 2。第三信号方向校核。把反馈增益的正负号和方程里的符号对一遍。我自己的经验是90% 的“莫名发散”最后都出在单位换算和漏了负号上真正模型原理有问题的情况反而不多。6. 仿真结果怎么解读从开环响应看车辆稳定性6.1 阶跃转向输入下的瞬态与稳态响应给模型一个阶跃前轮转角输入观察横摆角速度 r 的时域响应这是整车操纵稳定性分析最经典的测试。响应曲线通常包含三个阶段上升段、可能的超调振荡段、稳态段。继续用前面的参数算特征值。A 矩阵的迹是 -7.99行列式约 32.33特征值约 -3.995 ± 4.046j。这说明系统是一个稳定的欠阻尼二阶系统无阻尼自然频率约 5.7 rad/s约 0.9 Hz阻尼比约 0.70。按照二阶系统公式估算超调量约 4.6%上升时间约 0.3 秒稳态值 0.275 rad/s。把这些理论值和仿真曲线对照看能很直观地理解“欠阻尼”意味着什么。如果换一组参数让阻尼比变小比如增大车速 u特征值虚部会比实部下降得更慢振荡就会更明显。这也是为什么高速变线时车辆动态感觉更“晃”的数学解释。6.2 特征值与稳定性阈值的可视化分析把车速当成变量每个车速都算一遍特征值画在复平面上就是根轨迹。对不足转向车辆特征值实部始终为负但随着车速增加特征值逐渐向虚轴靠近说明稳定性裕度在下降。对过度转向车辆存在一个临界车速超过之后特征值实部变正系统失稳。临界车速的公式是u_crit sqrt(C_f · C_r · L² / (m · (b·C_r - a·C_f)))代入前面的参数b·C_r - a·C_f 44000分母 1500 × 44000所以临界车速超过 100 m/s说明这台车在正常车速范围内都是不足转向特性不会出现固有失稳。反过来如果某辆车的参数使 b·C_r - a·C_f 为负临界车速可能在 30~40 m/s 左右那这辆车高速行驶就非常危险这也是汽车工程里严格杜绝过度转向特性的原因。6.3 从增益差异理解不足转向与过度转向判定转向特性的另一个直观方法是看稳态横摆角速度增益随车速的变化趋势不足转向K_us 0增益先随车速增大到达特征车速后开始下降曲线有个明显的峰值。中性转向K_us 0增益与车速成正比直线上升。过度转向K_us 0增益随车速增大而加速上升在临界车速处趋于无穷。用当前参数算一下 20、30、40 m/s 三个车速的增益车速 m/s稳态横摆角速度增益趋势205.17上升段305.50接近峰值405.21下降段峰值出现在特征车速 u_ch sqrt(1 / K_us) ≈ 28.6 m/s。这就是为什么有些车在 100 km/h 左右变道时感觉特别灵活再往上反而感觉“迟钝”了一点的原因。这种用一个小模型就能解释的真实驾驶感受也是我觉得二自由度模型特别值得花时间吃透的地方。最后说一点个人体会。我建议初学者不要嫌二自由度模型太简单先在纸上把状态方程一步步推导一遍再用积分器和增益模块搭出来最后用 State-Space 模块交叉验证一遍。这个过程走下来你对状态空间表达、信号流方向、参数灵敏度、稳定性判据的理解都会比直接拖一个大而全的模型扎实得多。后续要扩展方向也很清楚把线性轮胎换成魔术公式进入非线性区把恒定车速假设放开做纵侧向耦合或者接 Carsim 做联合仿真。但不管往哪个方向走二自由度模型打底的理解能力都会让你少踩很多坑。
