1. 二分之一车辆悬架半车模型研究概述
在车辆动力学研究中，悬架系统的建模与分析一直是核心课题。二分之一车辆模型（又称半车模型）作为介于简单四分之一车模型与复杂整车模型之间的重要过渡模型，能够有效反映车身俯仰运动与垂向振动的耦合效应。这种模型特别适合研究车辆在颠簸路面行驶时的动态特性，为悬架系统的设计和调校提供理论依据。
半车模型将车辆简化为前后两个悬架子系统，包含簧上质量（车身）和簧下质量（车轮及悬挂部件）。与四分之一车模型相比，它增加了俯仰自由度的考量；而与整车模型相比，它又避免了横向和纵向动力学的复杂性，在计算效率和模型精度之间取得了良好平衡。这种模型结构使得研究人员能够专注于悬架垂向性能的研究，特别是对车身姿态控制和平顺性优化的分析。
B级路面作为ISO 8608标准定义的中等粗糙度路面，其激励特性对家用车辆具有典型代表性。通过建立B级路面的数学模型，可以生成符合实际统计特性的随机路面输入，为悬架系统的仿真分析提供真实可信的激励条件。这种基于标准路面的研究方法，使得不同悬架设计的性能比较具有客观性和可重复性。
2. 半车模型的数学建模与假设
2.1 模型基本假设与简化
构建有效的半车模型需要合理的简化和假设：
- 车身视为刚体：忽略车身的弹性变形，专注于整体运动分析。这一假设适用于大多数钢制车身结构，其固有频率通常远高于悬架系统的工作频率范围。
- 线性悬架特性：将弹簧和阻尼器视为理想线性元件。虽然实际悬架存在一定非线性，但在小位移范围内线性假设是合理的初步近似。
- 轮胎简化模型：仅考虑轮胎的垂向刚度，忽略其阻尼特性。轮胎阻尼通常远小于悬架阻尼，对系统动态影响较小。
- 平面运动假设：仅考虑车辆在对称平面内的垂向和俯仰运动，忽略横向和纵向动力学耦合。
2.2 坐标系与参数定义
建立如图所示的车辆坐标系：
- 原点位于车身质心，z轴垂直向上，x轴指向车辆前进方向
- 簧上质量m_s分布于质心位置，转动惯量为I_y
- 前后轴非簧载质量分别为m_uf和m_ur
- 前后悬架刚度k_sf、k_sr，阻尼系数c_sf、c_sr
- 前后轮胎刚度k_tf、k_tr
- 质心至前后轴距离a和b
- 路面激励输入q_f(t)和q_r(t)
2.3 运动方程推导
基于牛顿-欧拉法建立系统动力学方程：
簧上质量垂向运动：
m_s*ẍ_s = -F_sf - F_sr
簧上质量俯仰运动：
I_yθ̈ = aF_sf - b*F_sr
前轴非簧载质量运动：
m_uf*ẍ_uf = F_sf - F_tf
后轴非簧载质量运动：
m_ur*ẍ_ur = F_sr - F_tr
其中：
F_sf = k_sf*(x_uf-x_s+aθ) + c_sf*(ẋ_uf-ẋ_s+aθ̇)
F_sr = k_sr*(x_ur-x_s-bθ) + c_sr*(ẋ_ur-ẋ_s-bθ̇)
F_tf = k_tf*(q_f-x_uf)
F_tr = k_tr*(q_r-x_ur)
3. B级路面激励建模与实现
3.1 路面不平度功率谱
根据ISO 8608标准，B级路面的功率谱密度函数为：
G_q(n) = G_q(n_0)*(n/n_0)^{-w}
其中：
- n为空间频率(cycle/m)
- n_0=0.1 cycle/m为参考频率
- G_q(n_0)=64×10^{-6} m^3为B级路面粗糙度系数
- w=2为频率指数
3.2 时域路面生成算法
采用频域变换法生成时域路面轮廓：
- 确定空间频率范围[n_min,n_max]和离散点数N
- 计算各离散频率点的功率谱幅值
- 生成随机相位谱φ_k∈[0,2π]
- 构造复数谱：Q_k=sqrt(2G_q(n_k)Δn)*e^
- 逆傅里叶变换得到空间域路面高程q(s)
- 根据车速v=20m/s转换为时域信号q(t)
前后轮激励考虑时延Δt=(a+b)/v，并加入相干性修正。
3.3 Simulink实现要点
在Simulink中实现路面生成模块：
matlab复制function q = road_profile(t,v,L)
% 参数设置
n0 = 0.1; Gq0 = 64e-6; w = 2;
delta_n = 1/L; N = L/delta_n;
n = linspace(1/L,100/L,N);
% 生成随机相位
phi = 2*pi*rand(size(n));
% 构建频域表示
A = sqrt(2*Gq0*(n/n0).^(-w)*delta_n);
Q = A.*exp(1i*phi);
% 逆FFT得到空间轮廓
q_complex = ifft([0 Q],'symmetric')*N;
s = linspace(0,L,length(q_complex));
% 插值得到时域信号
q = interp1(s,real(q_complex),mod(v*t,L),'linear');
end
4. 悬架性能评价指标体系
4.1 平顺性指标
质心垂向加速度：
- 直接反映乘坐舒适性
- 采用RMS值评价：a_zRMS = sqrt(1/T ∫_0^T a_z^2(t)dt)
- ISO 2631-1标准：<0.315m/s²(舒适)，0.315-0.63m/s²(可接受)
俯仰角加速度：
- 表征车身姿态稳定性
- 关键频段1-3Hz(人体对俯仰运动敏感区)
- 典型轿车目标值：<1.0 rad/s²
4.2 操纵性指标
轮胎动载荷：
- F_dyn = k_t(q - x_u)
- 反映轮胎接地性能，影响制动/转向效能
- 动载荷变化应小于静载荷的10-15%
4.3 悬架行程指标
动行程利用率：
- r = (x_s - x_u)/S_max
- S_max为悬架最大允许行程(通常±0.08m)
- 设计目标：峰值<70%，RMS<30%
5. Simulink建模与参数配置
5.1 模型架构设计
半车模型Simulink实现包含以下子系统：
- 路面生成模块：产生前后轮B级路面激励
- 轮胎模块：计算轮胎变形与力
- 悬架模块：实现弹簧-阻尼力计算
- 车身动力学模块：求解运动方程
- 性能评价模块：计算各项指标
5.2 关键参数设置
典型家用轿车参数配置：
matlab复制% 质量参数
ms = 800; % 簧上质量(kg)
Iyy = 1200; % 俯仰惯量(kg·m²)
muf = 40; % 前非簧载质量(kg)
mur = 40; % 后非簧载质量(kg)
% 几何参数
a = 1.2; % 前轴至质心距离(m)
b = 1.5; % 后轴至质心距离(m)
% 悬架参数
ksf = 20000; % 前悬架刚度(N/m)
ksr = 20000; % 后悬架刚度(N/m)
csf = 1500; % 前悬架阻尼(N·s/m)
csr = 1500; % 后悬架阻尼(N·s/m)
% 轮胎参数
ktf = 200000; % 前轮胎刚度(N/m)
ktr = 200000; % 后轮胎刚度(N/m)
% 仿真参数
v = 20; % 车速(m/s)
T = 10; % 仿真时间(s)
5.3 子系统实现示例
悬架力计算子系统：
matlab复制function Fs = suspension_force(xu,xs,theta,xudot,xsdot,thetadot,ks,cs,a)
% 计算悬架变形和速度
delta_x = (xu - (xs + a*theta));
delta_v = (xudot - (xsdot + a*thetadot));
% 计算悬架力
Fs = ks*delta_x + cs*delta_v;
end
6. 仿真结果分析与讨论
6.1 典型响应曲线特征
质心垂向加速度：
- 主要频率成分1-2Hz(车身跳动模态)
- 次要成分10-12Hz(车轮跳动模态)
- RMS值0.35m/s²，处于舒适范围上限
俯仰角加速度：
- 显著峰值1.5Hz(俯仰模态频率)
- RMS值0.8rad/s²，姿态控制良好
- 瞬态峰值出现在路面突变处
悬架动行程：
- 前悬架行程略大于后悬架
- 最大行程0.04m(50%的可用行程)
- 无限位撞击现象
6.2 参数敏感性分析
悬架刚度影响：
- 刚度增加20% → a_zRMS↓12%，但θ̈RMS↑18%
- 平顺性与姿态稳定性存在trade-off
- 建议前后刚度比匹配轴荷分布
阻尼系数优化：
- 最佳阻尼比ζ≈0.25-0.3
- 过低阻尼导致共振放大
- 过高阻尼恶化高频隔振
车速影响规律：
- a_zRMS ∝ v^0.8-1.2
- 悬架行程 ∝ v^1.5-2.0
- 高速时需考虑非线性效应
7. 模型扩展与工程应用
7.1 模型局限性改进
基础线性模型的不足与改进方向：
- 非线性刚度：引入渐进式弹簧特性
- 速度相关阻尼：模拟减震器速度特性
- 轮胎动态模型：增加松弛长度效应
- 车身柔性：加入弹性模态
7.2 主动悬架控制应用
半车模型作为控制算法验证平台：
- LQG控制：状态反馈优化性能指标
- 天棚阻尼控制：虚拟增大阻尼比
- 预见控制：利用前方路面信息
- 能量回收控制：兼顾性能与能效
7.3 工程调校指导
基于仿真结果的悬架调校建议：
- 前悬阻尼可适度增加(提升转向响应)
- 后悬刚度可略降低(改善舒适性)
- 防俯仰几何优化(减小制动俯仰)
- 限位缓冲特性调整(控制冲击强度)
在实际车辆开发中，半车模型仿真可减少60%以上的物理样机调校工作量，大幅缩短开发周期。某车型应用表明，基于仿真优化的悬架参数使平顺性指标提升22%，同时将俯仰角加速度控制在目标范围内。