matlab/simulink:车辆模型：四分之一车被动悬架的双质量(二自由度)simulink模型。 输入为路面不平度，输出为车轮加速度、车身加速度、车轮动载荷。 需要matlab2016b及以上版本
在汽车动力学研究中，四分之一车被动悬架的双质量（二自由度）模型是一个经典且重要的模型，它能帮助我们理解车辆悬架系统的基本动态特性。本文将详细介绍如何在Matlab/Simulink环境下搭建这一模型，并分析其输入输出特性。需注意，本文的模型搭建基于Matlab 2016b及以上版本。
一、模型原理
四分之一车被动悬架双质量模型主要考虑了车身质量和车轮质量，简化地描述了车辆悬架系统在垂直方向上的运动。它有两个自由度，分别对应车身的垂直位移和车轮的垂直位移。路面不平度作为输入，会引起车轮和车身的振动，而我们关心的输出则是车轮加速度、车身加速度以及车轮动载荷。
二、Simulink模型搭建
1. 输入模块
首先，我们需要引入路面不平度作为输入。在Simulink库中，可以使用“From Workspace”模块来实现。假设我们已经在Matlab工作区中生成了路面不平度的数据，数据格式可以是时间序列。例如，在Matlab命令行中输入以下代码生成一个简单的路面不平度信号：
t = 0:0.01:10; % 时间向量，从0到10秒，步长0.01秒 road_profile = 0.05*sin(2*pi*0.5*t); % 简单的正弦路面不平度信号然后将t和road_profile作为列向量组成一个矩阵data，并在“From Workspace”模块中设置数据输入。
2. 模型主体搭建
双质量模型的动力学方程可以通过牛顿第二定律推导得出。设车身质量为 $ms$，车轮质量为 $mu$，悬架弹簧刚度为 $ks$，悬架阻尼系数为 $cs$，轮胎刚度为 $k_t$。
车身的动力学方程为：
\[ ms \ddot{z}s = -ks (zs - zu) - cs (\dot{z}s - \dot{z}u) \]
车轮的动力学方程为：
\[ mu \ddot{z}u = ks (zs - zu) + cs (\dot{z}s - \dot{z}u) - kt (zu - r) \]
其中，$zs$ 是车身垂直位移，$zu$ 是车轮垂直位移，$r$ 是路面不平度。
在Simulink中，我们可以使用积分器模块来实现从位移到速度、加速度的转换。例如，对于车身位移 $zs$，通过两个积分器模块依次得到车身速度 $\dot{z}s$ 和车身加速度 $\ddot{z}s$。同样的方法应用于车轮位移 $zu$。
弹簧力和阻尼力可以通过相应的数学运算模块来计算。例如，弹簧力 $F{spring} = ks (zs - zu)$，可以通过减法模块得到位移差，再通过增益模块乘以弹簧刚度 $ks$ 得到。阻尼力类似计算，$F{damping} = cs (\dot{z}s - \dot{z}_u)$。
轮胎力 $F{tire} = kt (z_u - r)$，也是通过减法和增益模块实现。
3. 输出模块
车轮加速度、车身加速度以及车轮动载荷作为输出。车轮加速度可以直接从车轮位移积分得到的加速度信号获取，车身加速度同理。车轮动载荷可以通过轮胎力计算得到，即 $F{dynamic} = kt (z_u - r)$。这些输出信号可以连接到“To Workspace”模块，以便后续在Matlab中进行分析和绘图。
三、模型仿真与分析
搭建好模型后，设置合适的仿真参数，如仿真时间等。运行仿真后，我们可以在Matlab工作区中获取输出数据。例如，通过以下代码可以绘制车身加速度随时间的变化曲线：
figure; plot(tout, body_acceleration); xlabel('Time (s)'); ylabel('Body Acceleration (m/s^2)'); title('Body Acceleration in Quarter - Car Model');通过分析这些输出曲线，我们可以了解车辆在不同路面不平度输入下的振动特性，评估悬架系统的性能，比如车身加速度反映了乘坐舒适性，车轮动载荷影响轮胎的接地性等。
四分之一车被动悬架双质量（二自由度）Simulink模型为我们研究车辆悬架系统提供了一个基础且有效的工具，通过不断调整模型参数和输入信号，我们可以深入探索悬架系统的各种动态特性。
