做整车底盘相关研究的朋友，应该都有同感：悬架这个题目，看着简单，真正弄明白却要绕不少弯。它既涉及振动理论、路面激励，又牵扯舒适性与操稳等多目标权衡。我在日常项目里，一提到车辆悬架建模仿真，首选几乎都是Matlab/Simulink——这套环境从建模、仿真到数据分析的链路很顺，调试周期短，团队协作时也容易统一平台。很多应届生或刚从试验岗转仿真岗的工程师，第一篇仿真报告就是从1/4车悬架模型开始的。这篇文章把我这些年用Simulink搭悬架模型的思路、参数、调试经验整理出来，给想入门车辆动力学仿真、或者已经搭过但总觉得结果不对的朋友做个参考。

1. 悬架仿真这件事，难在哪、怎么切

1.1 为什么悬架是整车仿真的“第一课”

很多人会问：整车动力学那么多内容，转向、制动、轮胎、空气动力学都比悬架“看着酷”，为什么偏偏悬架是入门的首选？我的理解是，悬架模型恰好踩在“物理清晰”和“结果可验证”的平衡点上。

悬架系统的核心功能有三条：支撑车身、衰减路面冲击、维持车轮与路面附着。这三条分别对应三个可量化的评价指标——车身加速度、悬架动挠度、轮胎动载荷。也就是说，你搭出来的模型到底行不行，不需要上试车场，直接看这三个指标的曲线就能判断。这种“模型-指标-试验”闭环的学习路径，在整车领域里是最容易建立直觉的。

从模型复杂度上看，悬架又有天然的层级结构：1/4车模型（quarter-car）两自由度，适合理解振动本质；半车（half-car）模型四自由度，加入俯仰运动；整车模型则把左右轮、前后轴都考虑进去，适合做车辆动力学联合仿真。这种由简到繁的升级路径，正好对应新手到熟手的成长过程。

我个人建议，第一次做悬架仿真，不要一上来就搭整车。先在1/4车模型上把方程、模块、评价指标全部吃透，再往半车和整车扩展，效率反而高得多。

1.2 建模仿真的完整闭环：从物理简图到结果判读

车辆动力学仿真经常出现“模型搭起来了，但说不清结果对不对”的情况。究其原因，是跳过了中间的物理建模和量纲检查环节。我的习惯是严格按照下面六步走：

1. 画物理简图，标清楚各质量、刚度、阻尼、路面输入方向；
2. 对簧上质量和簧下质量分别列牛顿方程；
3. 整理成仿真可用的标准形式（要么解出最高阶导数，要么写成状态空间）；
4. 在Simulink中搭模块，先做开环检查，不做任何控制；
5. 设计工况（正弦扫频、随机路面、脉冲输入），设置求解器参数；
6. 用评价指标判读结果，反推模型和参数是否合理。

这套流程看着朴素，但能挡住80%的返工。很多新手在第二步和第三步偷懒，直接在Simulink里“边想边画”，结果符号错了、单位错了，仿真出一堆无穷大，回头查起来反而更费时间。

2. Simulink搭建悬架模型的三种思路

2.1 从微分方程到积分器：最直观也最容易上手

1/4车两自由度模型的运动方程是这个形式的，簧上质量方程：

m_s * z_s'' + c_s * (z_s' - z_u') + k_s * (z_s - z_u) = 0

簧下质量方程：

m_u * z_u'' + c_s * (z_u' - z_s') + k_s * (z_u - z_s) + k_t * (z_u - z_r) = 0

其中z_s是车身位移，z_u是车轮位移，z_r是路面输入，m_s是簧上质量，m_u是簧下质量，k_s是悬架刚度，c_s是减振器阻尼系数，k_t是轮胎等效刚度。

在Simulink里搭积分器模型，核心思想是把最高阶导数解出来，然后用积分器串联得到速度和位移。对第一个方程做变换：

z_s'' = (-c_s * (z_s' - z_u') - k_s * (z_s - z_u)) / m_s

搭法我展开讲一下：

• 拖入两个Integrator模块，前一个输出z_s'（车身速度），后一个输出z_s（车身位移）；
• 在第一个积分器输入端，用Sum模块做负反馈组合：把z_s' - z_u'乘以-c_s/m_s、z_s - z_u乘以-k_s/m_s加到一起；
• z_u和z_u'来自车轮侧的子模型，不要自己猜数值，用Goto/From或信号线拉过去。

车轮侧同理，解出z_u''：

z_u'' = (-c_s * (z_u' - z_s') - k_s * (z_u - z_s) - k_t * (z_u - z_r)) / m_u

这里需要两个积分器生成z_u'和z_u，再从系统外部引入路面输入z_r。

这套搭法最大的好处是模块结构和物理一一对应，调试时哪个信号不对就能立刻定位。缺点也很明显：信号线数量多，模型放大后比较乱；改动参数需要进Gain模块一个个改，不方便做批量扫参。如果是临时验证物理概念，我推荐这种方式。

注意：积分器初始条件默认是0，仿真开始前如果悬架处于静平衡位置，z_s和z_u的初值都保持默认0即可。但如果你要加入车身初始位移，必须想清楚是相对静平衡还是绝对位移，否则第一个时间步就会产生很大的瞬态响应。

2.2 状态空间法：参数化调参和接入控制的利器

如果要做控制器设计或者批量参数扫描，我强烈建议改用状态空间法。把系统的四个状态定义为：

x = [z_s, z_s', z_u, z_u']'

输入u = z_r，状态空间方程写成：

x' = A * x + B * u

A矩阵和B矩阵分别为：

A = [0 1 0 0; -k_s/m_s, -c_s/m_s, k_s/m_s, c_s/m_s; 0 0 0 1; k_s/m_u, c_s/m_u, -(k_s+k_t)/m_u, -c_s/m_u]

B = [0; 0; 0; k_t/m_u]

仿真时在Simulink中直接拖入State-Space模块，把A、B按上述矩阵填入，再接一个常量输入或路面信号即可。如果需要输出悬架动挠度、轮胎动载荷，可以扩展C矩阵：

y = [z_s; z_s'; z_u; z_u'; z_s - z_u; k_t*(z_u - z_r)]

对应的C和D矩阵：

C = [1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 1; 1 0 -1 0; 0 0 k_t 0]

D = [0; 0; 0; 0; 0; -k_t]

这套做法的优势在后续特别明显。比如你想在整车控制器里加一个天棚阻尼（skyhook）控制律，只需要把状态向量引出来，在状态空间模块外部做反馈，再由控制器输出一个控制力叠加到A矩阵对应的力项上。此时物理结构没变，只是状态反馈发生了变化，完全不用重新搭积分器模型。

另外，在MATLAB脚本里可以用ss(A,B,C,D)直接建立悬架模型，配合bode、step等函数在频域里快速查看响应，再导入Simulink做时域验证。可以说状态空间法是“建模、分析、验证”三者通用的语言。

2.3 传递函数法：适合频域分析，但别滥用

第三种思路是把轮胎到车身的传递函数直接写出来。忽略簧下质量动态，只考虑单自由度简化时，车身位移对路面输入的传递函数近似为：

H(s) = (c_s * s + k_s) / (m_s * s^2 + c_s * s + k_s)

这个式子写进Transfer Fcn模块，一个模块就能完成悬架模型，特别适合做频率响应分析和隔振性能校核。还可以用MATLAB的bode函数快速绘制幅频特性，定量观察车身共振峰和隔振起始频率。

但我不建议在复杂工况和多输入场景下用传递函数法。原因有两点：第一，它天然是线性时不变模型，没法处理非线性阻尼、限位块、变刚度弹簧；第二，路面激励往往是多点多输入，传递函数法对单一输出点方便，但要同时得到动挠度和轮胎动载荷就需要搭多个传递函数，模型维护成本会上升。我的建议是，传递函数法只用来做“快速理解和频域验证”，正式项目里的主模型优先考虑状态空间或积分器结构。

2.4 悬架模型参数怎么定才靠谱

建模时参数不准，仿真结果自然没有参考价值。这里给一组典型的乘用车1/4车模型参数，也是我平时教学和预研时常用的基准值：

拿到参数后，可以快速估算两个关键频率。车身固有频率约为：

f_s = (1 / (2π)) * sqrt(k_s / m_s) ≈ 1.35 Hz

车轮系统固有频率约为：

f_u = (1 / (2π)) * sqrt((k_s + k_t) / m_u) ≈ 11.1 Hz

这两个频率决定了悬架的隔振区间。车身共振峰在1.3 Hz附近，一般来说需要靠阻尼压住；11 Hz附近是车轮共振区，如果阻尼设计不好，轮跳明显，轮胎动载荷就会很大。仿真结果里如果峰值频率和这两个理论值偏差过大，多半是参数填错了，首先要检查单位。

我踩过最典型的坑：把质量直接填成重力称重得到的“公斤数”当然没问题，但有人会把重量单位搞混，比如把重力N当作质量kg填进去，结果系统固有频率直接偏高近三倍。Simulink里所有物理量一律用SI制，养成习惯能省掉大量排错时间。

3. 路面激励与仿真配置实操

3.1 正弦输入：最稳的基础测试

验证模型第一个可以用的工况就是正弦激励。在Simulink中用Sine Wave模块，幅值设置0.02 m，频率从0.5 Hz开始逐渐增加到25 Hz，看车身加速度和悬架动挠度的变化。

更高效的扫频方式是用Chirp信号（Simulink里叫Chirp Signal模块），设置起始频率0.5 Hz，终止频率25 Hz，扫描时间20 s。这样一次仿真就能看出系统共振位置。如果你的模型是正确的，车身加速度幅值应该在1.3 Hz附近出现第一个峰，在11 Hz附近出现第二个峰。

我建议把Chirp扫频作为模型的“冒烟测试”，任何一次对模型的改动之后，先跑一遍扫频，确认共振频率没有大幅漂移，再去做随机路面和脉冲工况。这一步看起来简单，但对排查参数改动引入的错误非常有效。

3.2 随机路面：用白噪声生成真实不平度

车辆实际行驶的路面是随机过程，工程上常用路面功率谱密度描述。ISO 8608标准将路面按粗糙程度分为A到H级，C级路面代表一般公路，D级代表破损较严重的路面。如果想了解路面不平度的影响，可以用简化的白噪声成型滤波法生成时域信号。

简化思路是：路面不平度可以看作白噪声经过一个低通滤波器，滤波器截止频率与车速和路面波长有关。在Simulink中可以这样搭：

• 用Band-Limited White Noise模块，Noise power设为G0*v，其中G0是路面谱参考值，v是车速；
• Sample time设为0.001 s；
• 后面串联一个一阶低通滤波器，传递函数为1/(s/w_c + 1)，其中w_c = 2piv*n0，n0取0.1 1/m；
• 滤波器输出即为近似的路面位移时域序列。

以C级路面、车速20 m/s为例，G0取256e-6 m^3，算得截止频率约12.6 Hz。噪声的幅值需要根据实际道路粗糙度做标定，更严格的办法是用FFT合成路面高程序列，再通过From Workspace模块输入Simulink。

我也提供一个简单的MATLAB脚本，用于生成路面高程序列：

% 生成C级随机路面的简化时域序列 v = 20; % 车速 m/s dt = 0.001; % 时间步长 s t = (0:dt:20)'; % 20秒时长 n0 = 0.1; % 参考空间频率 G0 = 256e-6; % C级路面谱参考值 w_c = 2*pi*v*n0; % 截止频率 A = sqrt(G0*v/(2*pi)); rng(2024); zr = zeros(size(t)); for k = 2:length(t) zr(k) = exp(-w_c*dt)*zr(k-1) + A*sqrt(2*w_c*dt)*randn; end

脚本并不追求严格的高保真路面谱，但足够用来对比不同悬架参数在随机激励下的响应差别。跑完把zr变量保存到工作区，在Simulink里用From Workspace模块导入，注意设置采样时间与脚本一致。

3.3 脉冲输入：模拟减速带和路肩

脉冲工况用来考察悬架在瞬态冲击下的表现。最简单的做法是用Signal Builder模块画一个单脉冲，比如持续0.05 s、幅值0.05 m的梯形波，相当于以一定车速经过减速带时的轮心位移输入。

我实测下来，脉冲响应对阻尼参数最敏感。阻尼小时，车身会出现多次衰减振荡；阻尼大时，首峰降低但高频振动可能增大。通过脉冲响应的峰值时间和收敛速度，可以直观判断阻尼设置是否合理，这也是我调试减振器参数时最常用的快速手段。

脉冲信号也可以由MATLAB Function模块生成。比如用高斯脉冲近似：

z_r = 0.05 * exp(-((t - 0.5).^2) / 0.01)

仿真时间设2 s就足够观察完整响应。注意高斯脉冲的宽度参数要根据车速和障碍物宽度调整，否则等效冲击速度不对。

3.4 求解器设置与步长选择

Simulink默认的ode45变步长求解器对大多数悬架模型都够用，但有几个细节会影响结果质量。

第一，路面高频成分越丰富，最大步长必须越小。建议把Max step size设置成最高关注频率对应周期的1/20。如果仿真关注到25 Hz，周期是0.04 s，最大步长设为0.002 s比较稳妥。

第二，如果模型里有高频轮胎动力学，或者后续要接控制器，推荐用固定步长求解器ode4（四阶龙格-库塔），步长取0.001 s或更小。固定步长便于对比不同组参数的结果，也不会因为变步长自适应产生时间轴不一致的问题。

第三，对脉冲输入，建议开启Zero Crossing Detection选项，避免信号突变处步长过大导致峰值失真。Simulink里默认是开启的，但如果手动改了求解器配置，记得检查一下。

我常用的配置如下：

这套配置在大多数入门项目中都不会出错，跑出来的曲线也足够平滑。

4. 结果怎么读：三个关键评价指标

4.1 车身加速度：乘坐舒适性的硬指标

车身垂向加速度是评价乘坐舒适性最直接的指标。常用的统计量是加速度均方根值RMS，ISO 2631标准给出了不同频段的人体敏感度权重，但工程预研阶段直接看整体RMS已经能反映趋势。

在Simulink中把车身加速度信号送入RMS模块，或者导出到工作区后用rms()函数计算。我用这个指标做过一个阻尼参数扫描实验：保持其他参数不变，把阻尼系数从1000变到2500 N·s/m，加速度RMS呈现先降后升的趋势，最低点在1500附近，说明这个阻尼值确实落在合理区间。

看曲线时还要注意，车身加速度高频成分往往来自簧下质量振动，不能只盯着峰值，更别在峰值上叠加高频毛刺。如果曲线出现明显抖动，先检查步长设置，再检查路面激励滤波是否到位。

4.2 悬架动挠度：检查是否顶到限位块

悬架动挠度定义为悬架行程的变化量，即z_s - z_u。它反映悬架在冲击工况下有没有压缩到底、拉伸到顶。设计悬架限位行程通常在±100 mm左右，动挠度如果频繁触及限位，说明弹簧太软或阻尼不足。

判读动挠度时，我习惯看两个值：一是峰值，判断是否超出物理限位；二是均方根值，判断悬架行程的平均利用率。在随机路面工况下，动挠度RMS大致控制在行程上限的1/3以内比较合理，否则悬架会频繁打底，舒适性和零部件寿命都会受影响。

动挠度和车身加速度之间存在典型矛盾：偏软的弹簧能降低车身加速度，但会增大动挠度；偏硬的弹簧反过来。仿真调参时如果只优化一个指标，往往会牺牲另一个，这也是悬架设计被称为“折中艺术”的原因。

4.3 轮胎动载荷：车轮抓地的底线

轮胎动载荷表达式为k_t*(z_u - z_r)，它衡量车轮与路面之间的接触力波动。用动载荷除以静态载荷，得到无量纲的相对动载。相对动载超过1，意味着轮胎在某个时刻离地了——这时候车辆失去侧向力和驱动力，非常危险。

我在做随机路面工况时，会在Scope里同时看轮胎动载荷和静态载荷线。静态载荷等于m_s+m_u的重力，即约395*9.8≈3871 N。如果动载荷峰值接近或超过这条线，说明轮胎“飘”了。车轮共振频率附近的激励最容易造成这种现象，所以路面输入包含11 Hz附近能量时，要特别关注轮胎动载荷。

从控制角度讲，悬架阻尼和刚度对轮胎动载荷的影响比较复杂，不能简单认为“阻尼越大抓地越好”。阻尼过大会“锁死”悬架，反而把路面冲击直接传给轮胎，增大动载荷波动。多组参数对比时，我会把三个指标放到一张表里综合评估。

4.4 多指标联合评估：一张表看透参数变化

把三个指标放在一起，结论会清晰很多。下面是一组典型扫描结果示例：

这组结果表明，1500 N·s/m附近是舒适性较好且动载不高的均衡点，而2500时轮胎动载荷明显恶化。简单说，仿真结果不能只看一列，三项指标联合起来才是完整评价。

5. 我在调试中踩过的坑

5.1 符号写反，模型瞬间发散

用积分器建悬架模型时最常见的错误是反馈符号搞反。我见过不少朋友把Sum模块的加减号接反，结果悬架变成正反馈，Simulink跑出一片无穷大。判断方法很简单：给0.1 m/s的初始速度，观察车身位移是否在0附近振荡衰减。如果发散，立刻检查反馈环路中所有和阻尼、刚度相关的符号。

一个排查技巧是把Kp项先设成0，只保留阻尼项，确认响应是收敛的再逐步把刚度加回来。这样能把错误定位到具体环节。

5.2 代数环带来的“仿真卡死”

在Simulink中如果信号回路里没有积分或延迟单元，就可能形成代数环，求解器会疯狂迭代甚至报错。我的悬架模型里出现过一次：轮胎力信号直接反馈给轮胎位移输入，中间没有经过积分器，形成瞬时循环。

解决办法通常是三类：在环路中插入Memory模块或小惯性环节打破代数环；调整信号流向，确保反馈信号取自积分器输出；对传递函数法模型，用较小的惯性时间常数近似。代数环最大的风险是结果看起来也能跑完，但其实引入了额外延迟，影响精度。

5.3 单位混乱：最隐蔽的低级错误

又一次模型跑出诡异结果，我查了一下午，最后发现是参数单位搞错了——弹簧刚度输入的是N/mm而不是N/m，差了1000倍。悬架模型的单位问题集中在三处：刚度、阻尼、路面输入。刚度如果写成N/mm，模型“软”得离谱；阻尼写成N·s/mm，系统变成过阻尼；路面输入用mm而位置信号用m，直接造成幅值错乱。

我的做法是在参数脚本里把单位写进注释，并用MATLAB的变量命名强制区分，比如选择ks_unit = 'N/m'。跑模型前先做一次量纲检查：刚度除以质量，开根号后频率应该落在0.5~2 Hz附近，不符合就说明单位有误。

5.4 Bus Selector 选不到信号

用Bus对象做信号打包后，你也可能遇到Bus Selector里“没有可选信号”的问题。这通常是因为总线信号没有正确命名，或者连接线没有形成真正的Bus信号。解决办法是在创建总线时给每个信号用BusCreator并显式命名，不要依赖默认的“1”“2”信号名。此外，如果信号经过Goto/From传递，要保证Bus对象定义在全模型可见范围，否则Bus Selector无法识别。

5.5 想更真实？从1/4车走向Carsim联合仿真

入门阶段把1/4车模型做好后，如果想更进一步，Carsim和Simulink联合仿真是非常自然的下一步。Carsim提供整车轮辋、转向、车身等完整参数，Simulink里可以保留自己写的悬架控制算法或自定义悬架模型，通过S-Function接口实现双向数据交换。用Carsim做整车动力学背景，用Simulink做控制器设计，是工程实践中很常见的技术组合。

但我还是建议先把1/4车模型吃透再上联合仿真。Carsim带来的整车自由度非常多，如果连悬架三个基本指标都没建立直觉，面对整车结果很容易迷失方向。先学会走路，再学跑，这个顺序不会错。

最后说两句

这么多年搭悬架模型，我最大的体会是：仿真不是把模型“跑通”就结束了，而是要用仿真去理解悬架系统的物理本质。一个看似普通的1/4车模型，背后是固有频率、阻尼比、路面激励谱、多目标权衡这些底层概念。把这些概念在Simulink里亲手验证一遍，比看十篇理论文章都有用。你在参数扫描时会发现，弹簧、阻尼的调整方向和直觉完全一致；你也会发现，那些教科书里的折中曲线，真的不是纸上谈兵。我到现在做整车项目，遇到复杂悬架问题时，还是会先回到这个最小模型上做快速预判，它就像一个需要长期打磨但永远带着的瑞士军刀。希望这篇内容能帮你把这把刀打磨顺手。