Simulink整车建模必学：汽车行驶阻力计算子系统搭建与封装详解

作为一个长期做整车动力性仿真的人我有一个很深的体会很多初学者在 Simulink 里搭建整车模型时第一步往往不是纠结于发动机模型也不是纠结于换挡策略而是卡在了最基础、最枯燥的“行驶阻力计算”上。你可能觉得行驶阻力不就是套公式吗滚动阻力、空气阻力、坡度阻力、加速阻力四个一加不就完事了但真正用 Simulink 搭建子系统时你会发现单位换算怎么处理速度信号从哪里来坡度要不要做正余弦分解为什么我算出来的阻力数值和手算差了好几个数量级这些都是非常现实的问题。本文是《基于 Simulink 的汽车动力性能建模与仿真》系列的第 005 讲我们专门来拆解计算汽车行驶阻力子系统的完整建模过程。我会从物理公式出发逐步搭建 Simulink 模型再做子系统封装最后给出一套可复用的参数化模板。无论你是正在做课程设计、毕业设计还是在做整车性能仿真项目这篇文章都能帮你绕开常见的坑。1. 为什么要单独做一个“行驶阻力子系统”在很多整车动力学的 Simulink 模型中阻力计算往往被直接塞进一个巨大的“整车纵向动力学”模块里。这样做的后果是后期想单独分析某一种阻力对动力性的影响时你会发现模块内部一团乱麻信号线交叉纵横根本无从下手。1.1 行驶阻力在整车仿真中的位置我们先看一张简化的整车纵向动力学信号流驾驶员模型 → 踏板信号 → 动力源模型(发动机/电机) → 传动系模型 → 驱动力 ↓ 车速信号 ← 整车动力学模型 ← 合力(驱动力 - 行驶阻力) ↑ 行驶阻力计算子系统(本文建模对象)可以看到行驶阻力计算子系统处在整车动力学模型的反馈回路中。它接收车速、加速度、坡度等信号输出总行驶阻力。这个阻力会与驱动力做差再通过整车质量计算加速度最后积分得到车速。如果你搭建的是纯前向仿真模型这个子系统几乎是必做的一环。即使你使用的是后向仿真根据工况需求车速反推驱动力阻力计算模块同样不可缺少。1.2 四种行驶阻力的物理意义汽车在水平道路上等速行驶时必须克服来自地面的滚动阻力和来自空气的空气阻力。当汽车在坡道上上坡行驶时还需要克服坡度阻力。当汽车加速行驶时还需要克服加速阻力。这四种阻力的物理本质不同阻力类型产生的根本原因与什么因素有关滚动阻力 Ff轮胎与地面接触区的弹性迟滞损失车重、轮胎类型、路面条件、车速空气阻力 Fw车身表面与空气摩擦、前后压力差风阻系数、迎风面积、车速平方关系坡度阻力 Fi重力沿坡道方向的分量车重、坡度角/坡度值加速阻力 Fj汽车加速时惯性力车重、旋转质量换算系数、加速度在建模之前必须明白这四种阻力不是简单的“有没有”的关系而是在不同工况下的组合。例如汽车在平直路面匀速行驶时坡度阻力和加速阻力都可以视为 0在坡道起步时加速阻力和坡度阻力同时存在。子系统的建模目标就是输入必要的工况信号车速、加速度、坡度等在内部按照物理公式计算四种阻力并输出总阻力。2. Simulink 建模环境与准备工作在开始建模之前先确认一下环境。本文的示例以 MATLAB R2023b 之后的常见界面为例但实际上建模思路在 R2018b 到最新版本中都是通用的。2.1 软件版本与工具包搭建这个子系统只需要 MATLAB 和 Simulink 基础模块库不需要额外的 Simscape、Vehicle Dynamics Blockset 等专业工具箱。也就是说即使你只有最基础的 Simulink 许可证也可以完成本讲的全部内容。如果你有 Vehicle Dynamics Blockset当然可以直接调用现成的阻力模型块但那样反而不利于理解底层物理关系。本讲我们坚持“从零搭建”让你彻底搞清楚每个模块在干什么。2.2 示例车辆参数准备为了让模型可运行我们需要准备一组示例车辆参数。这里以某典型乘用车为参考参数如下% 车辆参数定义脚本vehicle_params.m % 建议在搭建模型前先运行将参数写入工作区 m 1500; % 整车质量单位 kg g 9.8; % 重力加速度单位 m/s^2 f 0.015; % 滚动阻力系数良好沥青路面 CD 0.32; % 空气阻力系数 A 2.2; % 迎风面积单位 m^2 rho 1.2258; % 空气密度单位 kg/m^3 delta 1.06; % 旋转质量换算系数请注意这里的参数不是固定的你在实际应用中应该根据具体车型修改。滚动阻力系数 f 与路面种类、车速和轮胎气压有关一般在 0.012 到 0.018 之间空气阻力系数 CD 与车身造型强相关轿车的 CD 值通常在 0.28 到 0.40 之间。2.3 建模的整体思路本讲的核心思路是“先搭建单个阻力计算模块再封装为子系统”。具体步骤如下先明确每个阻力公式的输入输出。在 Simulink 中分别搭建四个阻力计算模块。用 Add/Sum 模块合成总阻力。打包成子系统并用 Mask 实现参数化。连接 Signal Builder 或 Step 信号做仿真验证。下面我们进入核心的建模环节。3. 行驶阻力公式的核心理论拆解在动鼠标之前先把公式彻底搞懂。这里有一个很关键的细节国内教材和工程中车速单位常常混用有的用 km/h公里每小时有的用 m/s米每秒。这直接决定了空气阻力公式的系数。3.1 滚动阻力公式滚动阻力的经典公式为Ff W * f * cos(alpha)其中W 为车轮法向载荷在水平路面上 W m * gf 为滚动阻力系数alpha 为坡度角。在坡度不大时例如道路坡度小于 10%cos(alpha) ≈ 1因此很多工程模型中直接简化为Ff m * g * f在本文的建模中我们为了兼顾后续做坡道工况的扩展性会在模型中保留 cos 项。如果你确定始终在平路工况下仿真也可以去掉 cos 项使模型更简洁。3.2 空气阻力公式空气阻力的经典公式为Fw 0.5 * rho * CD * A * v^2但这个公式中车速 v 的单位必须是 m/s且这个阻力是相对风速产生的。国内教材中更常见的写法是Fw CD * A * v_a^2 / 21.15其中 v_a 的单位是 km/h。为什么会有 21.15 这个数字因为它已经将空气密度、单位换算等常数合并进去了。推导如下Fw 0.5 * rho * CD * A * (v_a / 3.6)^2 (0.5 * 1.2258 / 12.96) * CD * A * v_a^2 0.0473 * CD * A * v_a^2 ≈ CD * A * v_a^2 / 21.15所以在建模时必须先明确你模型中的车速信号单位是什么再决定用哪个版本的公式。如果混用结果会差 3.6 的平方倍也就是约 13 倍这是一个非常容易踩的坑。在本讲中我们采用 SI 单位制车速统一使用 m/s因此使用第一种公式Fw 0.5 * rho * CD * A * v^2。3.3 坡度阻力公式坡度阻力本质上是重力沿坡道方向的分量Fi m * g * sin(alpha)在实际工程中道路坡度 i 通常定义为坡道垂直高度与水平距离之比即 i tan(alpha)。当坡度角较小时例如 i 0.2tan(alpha) ≈ sin(alpha) ≈ alpha此时可以用Fi m * g * i这个简化公式在工程上非常常用。但如果你做的是大坡度工况仿真比如越野车爬陡坡就必须使用 sin(alpha)。本讲的建模中我会同时把两种形式讲清楚默认使用坡度值 i作为输入因为实际道路坡度数据往往以百分比形式给出并在模块内部完成 Fi m * g * i 的计算。3.4 加速阻力公式加速阻力包括平动质量惯性力和旋转质量惯性力两部分。工程上的简化公式为Fj delta * m * a其中delta 为旋转质量换算系数通常大于 1a 为车辆纵向加速度单位 m/s^2。delta 的物理含义是发动机飞轮、离合器、变速器齿轮、传动轴、车轮等旋转部件在加速时需要额外消耗一部分扭矩等效于整车质量增大了 delta 倍。对于普通乘用车delta 一般在 1.05 到 1.35 之间具体数值与挡位和传动比有关。工程上如果做粗略估算可以取常数 1.1 左右如果做精确仿真delta 可以根据挡位查表得到。3.5 总行驶阻力合成整车行驶阻力为F_total Ff Fw Fi Fj注意这里的合成是纯代数和。因为四种阻力都作用在车辆前进的反方向上方向一致。这看起来非常简单但在 Simulink 建模时你需要考虑单位制、信号维度、初始值等问题。下面我们进入实际操作。4. 行驶阻力子系统的完整建模实战现在我们开始搭建模型。本文采用的方案是先分别搭建四个阻力计算模块再用一个求和模块输出总阻力。为了便于后续复用我们会把整个阻力计算做成一个带 Mask 封装的可配置子系统。4.1 创建新的 Simulink 模型在 MATLAB 命令窗口输入simulink或者直接点击 MATLAB 主页上的 Simulink 图标打开 Simulink 开始页选择“Blank Model”创建一个空模型。保存模型文件名为Vehicle_Dynamics_Demo.slx然后在模型窗口中打开模块库浏览器快捷键 CtrlShiftL我们需要用到的模块主要来自这几个库Simulink → Sources信号源Simulink → Math Operations数学运算Simulink → Continuous积分、微分Simulink → Sinks信号输出Simulink → Signal Routing信号路由Simulink → Commonly Used Blocks常用模块4.2 搭建滚动阻力模块滚动阻力模型相对最简单。我们用 Constant 模块作为重力载荷输入Gain 模块乘以滚动阻力系数。具体步骤从 Simulink → Sources 中拖入一个Constant模块。从 Simulink → Math Operations 中拖入一个Gain模块。双击 Constant将数值设置为m*g注意这里可以直接写表达式Simulink 会在运行时计算。双击 Gain将增益设置为f。模块连接如下Constant(m*g) → Gain(× f) → 滚动阻力 Ff这里的 Constant 模块参数表达式中m 和 g 都来自工作区。如果你已经运行了前文中的 vehicle_params.m 脚本Simulink 会自动读取这些变量。为了后续子系统的参数化我们实际上不需要在 Constant 模块里写死 m 和 g。更好的做法是让 Constant 模块直接代表整车重力 W然后再乘以 f。这样后续用 Mask 封装时只需暴露 m、g、f 三个参数即可。为了让模型更清晰我给模块命名如下Constant 重命名为重力载荷 WGain 重命名为滚动阻力系数 f输出信号线标注为Ff4.3 搭建空气阻力模块空气阻力的核心是平方运算。拖入一个Inport模块作为车速输入命名为v。拖入一个Math Function模块函数选择square平方。拖入两个Gain模块。拖入一个Product模块。具体连接逻辑v → Math Function(square) → Gain(0.5*rho) → Gain(CD*A) → Fw也可以将两个 Gain 合并为一个v → Math Function(square) → Gain(0.5*rho*CD*A) → Fw我建议使用两个 Gain 分开相乘这样在查看模型时物理意义更明确第一个 Gain 是动压头 0.5rho第二个 Gain 是空气阻力系数与迎风面积的乘积 CDA。需要注意如果车速单位是 m/s那么空气阻力的单位是 N牛顿这个量级通常在几百牛顿到上千牛顿。如果你发现计算结果数值异常大或异常小首先检查车速信号单位。4.4 搭建坡度阻力模块坡度阻力的建模有两种方式完整版用 sin和简化版用坡度值 i。本讲采用简化版因为大多数道路坡度数据都以 i 的形式给出。拖入一个Inport模块作为坡度输入命名为i。拖入一个Constant模块数值设置为m*g。拖入一个Product模块。连接逻辑i → Product → Fi m*g → Product 的另一个输入如果你希望完整版可以在 Product 之前加一个 Trigonometric Functionsin模块将坡度角 alpha 换算为 sin(alpha)。这个我们留作扩展不展开。这里有一个工程细节道路坡度 i 是一个无量纲量例如 5% 的坡度记为 0.05。在 Simulink 中输入时可以直接用 Constant 或 Inport 提供。4.5 搭建加速阻力模块加速阻力模块稍微特殊因为它需要加速度信号 a。在整车动力学模型中加速度通常不是直接给定的而是由驱动力减阻力的合力除以质量得到。但在构建独立阻力模块时我们可以把加速度当作一个输入信号。拖入一个Inport模块命名为a。拖入一个Constant模块数值为delta*m。拖入一个Product模块。连接逻辑a → Product → Fj delta*m → Product 的另一个输入这里再提醒一次delta 是旋转质量换算系数通常大于 1。如果你发现加速阻力偏小可以先查一下 delta 是否取值太小。4.6 合成总阻力输出现在我们有了四个阻力信号 Ff、Fw、Fi、Fj。用 Add/Sum 模块将它们合成总阻力。拖入一个Add模块在 Math Operations 中。将 Add 的输入端口数量改为 4双击 Add设置 Inputs 为。分别将四个阻力信号连接到 Add 的输入。在 Add 输出端连接一个Outport模块命名为F_total。连接完成后整个模型内部结构应该如下Ff ──┐ Fw ──┼── Add() ── F_total Fi ──┤ Fj ──┘到这里一个最基本的计算行驶阻力子系统的内部逻辑已经搭好了。目前这个模型还只是散落的模块下面我们把它封装成一个完整的子系统。5. 子系统封装与参数化5.1 使用 Subsystem 模块封装你当然可以用鼠标框选全部模块然后右键选择“Create Subsystem from Selection”Simulink 会自动为你创建子系统。但更推荐的方式是提前规划好输入输出然后手动搭建 Subsystem。手动封装步骤如下从 Simulink → Ports Subsystems 中拖入一个Subsystem模块。双击打开子系统内部。在内部删除默认的 Inport/Outport按上述 4.2 到 4.6 小节的步骤重新搭建。将子系统的输入端口分别命名为v、a、i。将子系统的输出端口命名为F_total。此时子系统外部看起来就是一个三输入一输出的模块v ──┐ a ──┼── [阻力计算子系统] ── F_total i ──┘5.2 使用 Mask 实现参数化现在子系统内部还引用了工作区的变量 m、g、f、CD、A、rho、delta。为了让这个子系统在不同的车辆参数之间快速切换我们使用 Mask 编辑器进行封装。具体操作右键点击 Subsystem 模块选择Mask → Create System Mask。在 Mask 编辑器中选择Parameters Dialog标签页。添加 7 个参数分别对应 m、g、f、CD、A、rho、delta。在Initialization标签页中将变量的赋值代码写为m m; g g; f f; CD CD; A A; rho rho; delta delta;注意这里的左侧是内部变量名右侧是 Mask 参数名。如果你的 Mask 参数名与内部变量名相同那么这一行也不需要显式写Simulink mask 会自动把参数传递给内部的 Constant/Gain 模块表达式。为了让用户方便填写建议在 Mask 中给每个参数都加上描述文本。例如参数名描述默认值m整车质量 (kg)1500g重力加速度 (m/s^2)9.8f滚动阻力系数0.015CD空气阻力系数0.32A迎风面积 (m^2)2.2rho空气密度 (kg/m^3)1.2258delta旋转质量换算系数1.06封装之后你只需要在 Mask 参数对话框中修改数值整个阻力计算子系统的内部参数就会同步更新不需要再打开子系统一层层去改。5.3 添加额外测量输出可选扩展在工程调试时我们往往不仅需要总阻力还想观察某一种阻力分量的大小。为此可以在子系统内部添加一个Mux或Bus Creator将四个阻力分量聚合后一起输出。具体做法在子系统的输出侧添加一个Mux模块。将四路阻力信号连接到 Mux 的输入端。Mux 输出连接到另一个 Outport命名为F_components。此时子系统外部就有了两个输出端口F_total和F_components。这样可以同时观察总阻力和各分量对调试和验证非常有帮助。6. 仿真验证与结果分析模型搭好了参数也封装好了下面我们来验证模型的正确性。6.1 创建测试信号在子系统外部拖入几个信号源使用Step模块模拟车速信号 v 的变化。使用Constant模块给加速度 a 和坡度 i 一个固定值。或者更直观地使用Signal Builder模块创建多段工况信号模拟加速、匀速、爬坡等工况。为了方便对比我们先用简单的常量信号做手算验证。6.2 手算与仿真结果对比假设我们设置v 20 m/s约 72 km/ha 0.5 m/s^2i 0水平路面。基于前面的参数手算如下滚动阻力 Ff m * g * f 1500 * 9.8 * 0.015 220.5 N 空气阻力 Fw 0.5 * rho * CD * A * v^2 0.5 * 1.2258 * 0.32 * 2.2 * 400 172.5 N 坡度阻力 Fi m * g * i 0 N 加速阻力 Fj delta * m * a 1.06 * 1500 * 0.5 795 N 总阻力 F_total 220.5 172.5 0 795 1188 N在 Simulink 中你把上述信号输入后仿真结束时用 Scope 或 Display 模块观察得到的 F_total 应该稳定在 1188 N 左右。如果结果偏差很大说明模型内部参数或单位出了问题。下表是不同车速下的仿真验证参考车速 v (m/s)空气阻力 Fw (N)总阻力 F_total不含加速(N)1043.1263.620172.5393.030388.0608.540689.9910.4注意空气阻力随车速呈平方关系增长这就是为什么高速时空气阻力成为最主要阻力来源。6.3 使用 Scope 观察动态变化为了更直观地观察阻力变化过程建议在模型中加入Scope。将车速信号 v 用Ramp模块模拟线性增长例如从 0 开始以 1 m/s² 的斜率增加到 50。将加速度 a 用Step模块设置为 1。将坡度 i 设为 0。运行仿真 50 秒观察 Scope 中 F_total 的变化曲线。你会看到空气阻力随时间呈二次曲线增长加速阻力恒定滚动阻力恒定总阻力在低车速时以滚动阻力为主高车速时空气阻力快速占主导。7. 常见问题与排查思路在实际建模过程中读者经常遇到一些典型问题这里集中说明。问题现象常见原因解决思路仿真结果比手算大很多或小很多单位不统一比如车速用了 m/s 却套用了 km/h 公式统一使用 SI 单位制并在模块标注中写明单位加速阻力模块输出始终为 0加速度输入信号未正确连接或 a 是常量 0检查 Inport 是否连接到正确的信号源子系统 Mask 修改参数后模型无变化Mask 参数没有传给内部模块或表达式写成了常量检查 Mask Initialization 代码确认变量映射模型中出现代数环加速度信号由总阻力反馈计算而加速阻力又依赖加速度使用 Memory 模块或延迟模块打断代数环或从上游直接提供加速度空气阻力计算值偏小CD 或 A 的单位/数值不对rho 忘乘检查公式 Fw0.5rhoCDAv^2注意 rho 的数量级约 1.2仿真步长过小导致仿真缓慢阻力模型中有高频动态或 Derivative 模块使用固定步长求解器或避免使用 Derivative改用滤波后的加速度信号7.1 关于 Derivative 模块的警告很多初学者喜欢用 Derivative 模块对车速求导得到加速度。但在 Simulink 中Derivative 模块对噪声非常敏感容易导致仿真结果振荡甚至无法收敛。推荐做法是在整车模型中加速度信号应该由动力学方程直接计算得到即a (F_drive - F_resistance) / (delta * m)而不是对速度求导。如果只能获得速度信号可以使用一阶低通滤波器对速度做平滑处理后再求导但要注意滤波延迟。7.2 坡度输入的处理细节如果你使用的是坡度角 alpha单位 rad而不是坡度值 i需要增加一个 sin 模块。很多工程数据给出的是百分比坡度例如 5% 表示 i 0.05此时 Fi m * g * i。两者不要混用。8. 工程实践建议与后续扩展到这里我们已经完成了一个完整的汽车行驶阻力计算子系统。最后再给你几点工程实践建议以及后续学习的方向。8.1 命名规范与模型管理在实际项目中模型里的信号线、模块名、变量名都要有清晰规范。我建议信号线命名使用“物理量_单位”的格式例如v_mps、F_total_N。车辆参数统一用一个初始化脚本管理不要在 Constant 模块里写死数值。每个子系统都要有注释说明输入输出和公式来源。模型文件命名建议包含版本号例如Vehicle_Dynamics_Demo_v1_0.slx。8.2 从阻力子模型到完整整车动力学模型完成阻力模型之后下一步可以搭建完整的整车纵向动力学模型。整车纵向运动方程如下delta * m * dv/dt F_drive - F_total其中 F_drive 来自发动机/电机扭矩经过传动系传递到驱动轮的驱动力。你可以在 Simulink 中搭建如下结构发动机扭矩 → 变速器 → 主减速器 → 驱动力 F_drive ↓ [阻力计算子系统] → F_total → 求和 → 加速度 → 积分 → 车速这就是一个最基本的整车前向仿真模型。基于这个模型你可以继续研究最高车速仿真最大爬坡度仿真原地起步加速时间仿真NEDC/WLTC 循环工况仿真8.3 进阶方向当你已经完全掌握阻力子系统的搭建后可以考虑以下进阶方向将阻力计算从静态公式升级为动态查表模型滚动阻力系数 f 可以做成随车速变化的查表空气阻力系数 CD 可以做成随高度/风速变化的动态参数。与 Simscape 结合将机械传动系统用 Simscape 搭建阻力计算作为外部力输入可以模拟更真实的物理连接。使用 MATLAB App Designer 搭建 GUI将 Simulink 模型嵌入到 GUI 中实现参数输入和仿真结果可视化。代码生成如果你的目标是嵌入式实时仿真可以学习 Simulink Coder将阻力计算模型生成 C 代码。以上内容就是本讲关于“计算汽车行驶阻力子系统建模”的完整讲解。从物理公式到 Simulink 实现从单个模块到参数化封装再到仿真验证和常见问题排查每一步都是可以直接复用的。如果你在搭建过程中遇到卡住的地方优先检查单位制和信号连接这两个是新手最常见的问题根源。你可以在自己的 MATLAB 环境中动手搭建一遍然后把车辆参数换成你手头的车型数据验证一下不同工况下的阻力分布特性。如果有任何疑问欢迎在评论区交流讨论。