前阵子整理了一套关于“电动汽车等速工况续驶里程仿真”的 Matlab 工程包带注释、带说明书、能直接运行拿来当学习参考正好合适。做整车能耗或者电池系统相关工作的朋友应该都清楚这类“看起来简单”的仿真包真正值钱的地方反而不是能跑出数而是能把从物理公式到代码实现之间的那层逻辑捋清楚让别人拿到以后能复现、能核对、能继续改。这篇文章就把我当时设计这套包时踩过的坑、做过的取舍、还有验证结果的方法完整拆开讲一遍。这套仿真解决的问题很聚焦一台电动汽车在水平平直道路上以恒定速度巡航不考虑坡道、侧风和驾驶循环从起始 SoC 放电到终止 SoC一共能跑多少公里。为什么先做这个窄场景因为等速工况是整个能耗仿真体系里最基础、最能被手算核验的场景也是新手入门续驶里程计算最好的切口。比起直接上 Simulink 搭整车模型把这套“阻力—功率—能量—里程”的链路先吃透后续再做综合工况会顺很多。1. 为什么等速工况值得单独做成一版仿真交付收到项目需求时第一反应往往是“等速工况太简单了不就是阻力乘以速度再除个效率吗”。但实际把交付标准拉到“可供参考学习、有说明书、全套可复现”的高度以后你会发现简单场景反而对工程纪律要求更高。代码里每一个参数都要有物理含义每一行注释都要解释“为什么这么算”说明书要把公式重新推导一遍包括单位换算这种最容易翻车的地方全都得写清楚。这套包的使用场景主要有三类一是给刚接触电动汽车能耗计算的人做入门教程从一行行代码里看懂续航是怎么算出来的二是做整车预研时快速建立“速度—续航”对照表判断一台新车的油耗水平大致落在什么区间三是作为后续动态工况仿真模型的底盘把等速这个子模块验证扎实了再去叠加加速、减速和再生制动逻辑错误才不会层层传染。我把仿真目标刻意限定得比较窄水平路、无坡道、无风、恒速巡航、电池放电至 SoC 下限即停止。窄目标带来的好处非常实际——代码可以写得很短短到每一行都能和物理公式一一对上。初学者拿到手看到的不是一个几百行互相调用的黑盒脚本而是一条能顺着注释从头走到尾的完整计算链。1.1 等速工况在整车研发里的角色很多朋友会问日常驾驶哪有一直匀速跑的等速工况太理想化了。这里要区分两个概念——综合工况循环和等速特征工况。综合工况是拿来评估一辆车贴近真实使用场景的能耗表现里面塞满了加减速、怠速、匀速段而等速特征工况是研发早期用来做快速判断的基础观测面把滚动阻力、风阻、传动效率这些因素干净地分离开。行业里做预研时通常会把 40、60、80、100、120 公里/小时这几个速度点挨个跑一遍得到一条“等速续航—车速”曲线。低速区间的续航会很高速度往上走以后空气阻力按速度平方增长轮端功率按速度三次方增长续航曲线掉得又快又陡。这条曲线本身就是一张非常有说服力的工程图表也是很多“匀速工况最大续航 XX 公里”这类宣传数据背后的物理来源。所以别嫌等速简单它一直是续航评估体系里的基础观测面。1.2 简化模型的适用边界这套包为了做到“可参考、可学习”主动做了两处简化暂时忽略电池温度对容量和内阻的影响电机效率先用一个综合常数代表。这样换来的是整个计算链条可以在一个小时内完成“读参数—跑代码—对结果”的全部动作。简化不等于随便估算。等速巡航时电机通常工作在高效率区间效率波动很小所以用常数是合理近似。但如果谁把这套模型直接拿去模拟城市综合工况误差就会明显放大。所以我特别在说明书里写了适用边界适合“中高速稳定巡航”的里程估算不代表综合续航更不代表用户实际开出来的续航。边界写清工具才不会被误用。沿着这层边界继续往下走核心就是那串能量账。2. 把续驶里程公式写成人话阻力、效率、可用能量三段式所有等速续驶里程计算都逃不开这个三段式先算阻力并换算成轮端功率再算电池可用能量最后用“能量除以功率得到时间时间乘以速度得到里程”。展开到工程层面链路是车辆阻力 → 轮端功率 → 电池侧功率 → 可用能量 → 持续放电时间 → 续驶里程每段看似都是初中物理但颗粒度怎么控制、单位怎么盯、参数怎么取才是真正拉开差距的地方。2.1 轮端阻力滚动阻力与空气阻力水平路面等速行驶车辆要克服两类阻力。滚动阻力与速度基本无关公式是F_roll f × m × gf 是滚动阻力系数沥青路面一般在 0.010 到 0.015 之间m 是整车质量包括整备质量加乘员和货物g 是重力加速度。空气阻力则随速度平方增长F_aero 0.5 × ρ × Cd × A × v²ρ 是空气密度通常取 1.225 kg/m³Cd 是风阻系数A 是车辆正面迎风面积v 是车速注意这里的单位必须用 m/s因为公式整体是国际单位制。总阻力就是两项相加F_total F_roll F_aero阻力不等于功率轮端功率等于总阻力乘以车速P_wheel F_total × v滚动阻力对应的功率只和速度的一次方成正比而空气阻力对应的功率和速度的三次方成正比。用一套典型中级车参数来感受一下滚动阻力功率从 60 公里/小时到 120 公里/小时只从大约 2.9 kW 涨到 5.9 kW空气阻力功率却从大约 1.8 kW 涨到 14 kW。所以高速段续航暴跌主凶就是风阻。速度翻一倍风阻功率翻八倍这也是“车速越高续航掉得越狠”的物理根源。2.2 从轮端到电池端效率的归并与取舍轮端机械功率要经过减速器、电机、逆变器最终从电池里取出来。严格做法是给每个环节单独建效率模型但在等速仿真这一版里我把它们合并成一个从电池端到轮端的综合效率常数记为 η_total。于是电池需要提供的功率P_battery P_wheel / η_totalη_total 的取值要特别小心。电机效率通常在 0.85 到 0.95 之间减速器在 0.95 左右逆变器在 0.97 左右三者相乘大约是 0.78 到 0.88。如果某个人在代码里写 η_total 0.95那多半是把某个环节漏算了或者把工况点想得过于理想。等速巡航点下电机效率往往能到 0.92 以上乘上减速和变流环节取 0.85 到 0.88 比较常见。最终取多少应当与自己手里的实测数据对齐不要拍脑袋。2.3 电池可用能量电压、容量和 SoC 窗口电池侧最容易出错的地方是把“标称容量”和“实际可用容量”混为一谈。假设电池额定电压为 U额定容量为 Q那么标称电量是U × Q / 1000 kWh但你不能把这块电全部用完。出于保护寿命的考虑动力电池只在某个 SoC 区间工作比如从 95% 放到 10%。于是可用能量 标称电量 × (SoC_start − SoC_end)举个例子一块 60 kWh 的电池如果写 SoC 从 1.0 放到 0计算用能量 60 kWh这在现实中基本不存在。多数乘用车实际可用窗口要按 0.85 到 0.95 折减。很多初版仿真里程多算了祸根就埋在这里。有了电池功率和可用能量持续时间t_hours 可用能量 / 电池功率里程Range v_kmh × t_hours整个链条只有三行却是整份 Matlab 代码的骨架。下一步就看怎么把这三行组织成一套能交付、能运行、能继续扩展的代码。3. Matlab 实现方案从纯计算到 SoC 递推的两级建模这套包里我准备了两级实现。第一级是纯公式计算适合快速出数、交叉核验第二级加入电池内阻和 OCV 曲线做逐秒递推更接近“仿真”这个词的真实含义。3.1 第一级工程估算法三行公式直接出结果第一级实现的核心是一个函数输入车辆参数结构体和目标车速输出续驶里程。代码压到最简后长这样function range_km cruiseRange(params, v_kmh) % cruiseRange 等速工况续驶里程快速估算 % 输入: params 车辆与电池参数结构体 % v_kmh 等速车速, km/h % 输出: range_km 续驶里程, km v_ms v_kmh / 3.6; % km/h 转 m/s F_roll params.f * params.m * params.g; % 滚动阻力, N F_aero 0.5 * params.rho * params.Cd * params.A * v_ms^2; F_total F_roll F_aero; % 总阻力, N P_wheel F_total * v_ms; % 轮端功率, W P_bat P_wheel / params.eta_total; % 电池侧功率, W E_avail params.U_bat * params.Q_Ah * (params.soc0 - params.soc_end) / 1000; % kWh t_h E_avail * 1000 / P_bat; % 小时 range_km v_kmh * t_h; end每一行都对应第 2 节的公式注释直接写物理含义。等速工况没有加减速、没有再生制动、没有换挡所以这段代码不需要任何循环一次算完。真正建议保留的参数输入区我单独放在脚本前面% 用户输入区整车参数 params.m 1500; % 整车整备质量载荷, kg params.g 9.81; % 重力加速度, m/s² params.f 0.012; % 滚动阻力系数 params.Cd 0.28; % 风阻系数 params.A 2.2; % 迎风面积, m² params.rho 1.225; % 空气密度, kg/m³ % 三电参数 params.U_bat 400; % 电池额定电压, V params.Q_Ah 150; % 电池容量, Ah对应 60 kWh params.soc0 0.95; % 起始 SoC params.soc_end 0.10; % 放电终止 SoC params.eta_total 0.88; % 电机逆变器减速器综合效率为什么把参数单独抽成结构体而不是散落成变量因为后续批量跑速度点时要循环调用函数化以后只需要在循环里改 v_kmh 一个参数。散落变量则会写出一堆重复代码还容易改一处漏一处。这里必须强调一个经验所有参与功率计算的单位必须统一。车速从 km/h 转到 m/s 是新手最容易漏的一步。漏了以后阻力数值直接偏小 3.6 倍轮端功率偏小 12.96 倍续航却会莫名胀大十几倍。很多“怎么跑都跑不出合理数”的怪现象最后查出来都是这一步的问题。3.2 第二级带 OCV 曲线和内阻的 SoC 递推仿真第一级假设电池端电压恒定这在实际中不成立。电池电压会随 SoC 下降而降低如果外部需求功率固定电压降低就意味着电流升高电流升高又进一步拉低端电压形成一种“加速放电”的非线性。所以第二级实现用逐秒递推% 等速续航递推仿真按固定时间步长更新 SoC dt 1; % 时间步长, s soc params.soc0; range_km 0; % OCV_from_SOC: 查表函数输入 SoC返回开路电压 % params.R_bat: 电池内阻, Ohm这里视为常数 while soc params.soc_end ocv OCV_from_SOC(soc); % 开路电压 P_bat P_wheel / params.eta_total; % 需求电池功率, W % 解 I P_bat / U_terminal % U_terminal ocv - I * R_bat % 得到: R_bat * I^2 - ocv * I P_bat 0 I (ocv - sqrt(ocv^2 - 4 * params.R_bat * P_bat)) / (2 * params.R_bat); soc soc - I * dt / (params.Q_Ah * 3600); % SoC 更新 range_km range_km v_ms * dt / 1000; % 里程累积, km end这里处理的是“需求功率恒定条件下端电压与电流的非线性关系”。方程来源很简单P_bat U × I而 U U_ocv − I × R_bat代入整理就是一元二次方程。求根时取较小的那个根因为大根对应过大的电流物理上不是稳定工作点。第二级比第一级多了电池内阻 R_bat 和一条 OCV-SoC 曲线。多数初版项目没有这些真实数据可以先做合理假设但必须在说明书里明确标为“假设值”并告诉用户后续拿到电芯测试数据后只需要替换查表函数就能升级模型。这样整个包的可扩展性一下就撑起来了。3.3 注释与参数化代码“可读”比“能跑”更重要“Matlab 有注释完美运行”这个承诺里含金量其实在“有注释”。脚本能跑只说明语法没错、流程通了一行行注释把“为什么这么算”讲清楚才是别人能长期复用的前提。我在这套包里坚持三个注释习惯。第一每个物理量的注释里写单位比如% 滚动阻力, N这样量纲错乱会在代码审查时一眼被揪出来第二每个公式前面写它来源的物理关系比如P_bat P_wheel / eta_total后面一定跟一句“电池功率 轮端功率 / 动力链总效率”第三所有输入参数集中在“用户输入区”不许散落在脚本中间改参数只改一处避免改漏导致结果对不上。这套规矩看起来琐碎实际调试和交付时的收益非常明显。有了代码下一步是验证。跑出来的数到底对不对凭什么说“可供参考学习”靠的就是手算与仿真互相印证。4. 校核与试运行结果和手算如何对齐数学模型是可以被“跑出来”的但能不能经得起手算的检验才是真功夫。我拿到任何仿真输出第一件事都是用笔和计算器把一组典型工况算一遍逐项对比。差在 1% 以内说明链路没问题差得远问题基本出在单位或者效率定义上。4.1 以 60 km/h 为例的纸面验算用第 3 节那组典型参数m 1500 kg、f 0.012、Cd 0.28、A 2.2 m²、ρ 1.225 kg/m³。先把车速换算成 m/s60 / 3.6 16.67 m/s滚动阻力F_roll 0.012 × 1500 × 9.81 ≈ 176.6 N空气阻力F_aero 0.5 × 1.225 × 0.28 × 2.2 × 16.67² ≈ 104.9 N总阻力约 281.5 N。轮端功率P_wheel 281.5 × 16.67 ≈ 4692 W取动力链总效率 0.88电池侧功率P_bat 4692 / 0.88 ≈ 5332 W电池 400 V、150 AhSoC 窗口 95% 到 10%E_avail 400 × 150 × 0.85 / 1000 51 kWh可持续时间t 51 / (5332/1000) ≈ 9.57 h里程Range 60 × 9.57 ≈ 574 km这个结果合理吗一台 1500 kg、60 kWh 的中级电动车在 60 km/h 等速巡航时续航接近 570 公里和同类车型宣传的“最大等速续航”处于同一量级。注意它比综合续航高不少这是等速模型的正常表现——没有加减速、没有怠速、没有空调和采暖负荷能量几乎全用在克服阻力和传动损耗上。4.2 多速度点结果对比与合理性判读套包里通常会批量跑一组速度点。用相同参数估算几个点结果是车速 km/h总阻力 N轮端功率 kW电池功率 kW续航 km60281.54.75.357480362.98.19.2445100467.713.014.8346120595.819.922.6271换算成百公里电耗100 km/h 附近约 14.8 kWh/100km120 km/h 约 18.8 kWh/100km和市面主流中级电动车的常用高速电耗区间同量级。看到这种量级对得上的数字我才放心把仿真包交出去。速度从 60 提到 120只翻了一倍续航却从 574 掉到 271跌了一半以上。这种非线性掉电行为正是风阻功率按速度三次方增长的直接体现。如果某套等速仿真跑出来“120 km/h 还能有 450 km”那基本可以断定效率取错了或者风阻系数写小了。4.3 等速模型最容易出现的几个“假错误”所谓“假错误”就是代码不报错、结果却脱离物理常识的情况。我总结最典型的三种。第一种滚动阻力系数抄成 0.12 而不是 0.012。低速时还能被风阻掩盖速度一拉高数据整体偏离真实值。第二种用了 60 kWh 名义电量却忘了乘 SoC 窗口120 km/h 续航凭空多出约 18%。第三种效率定义用反——把电机效率、逆变器效率、传动效率这三者的乘积理解成“1 除以三者之和”或者干脆漏掉其中一个环节。如果把效率从 0.88 错写为 1.14 附近电池侧功率会偏低续航偏高而且速度越高偏差越大。这类问题用肉眼在代码里很难找出来但拿手算结果一对比一分钟就能定位。所以说明书里我把“先手算、再跑仿真”列为标准动作不是可选项。验证顺利通过以后这套包就可以宣称“可参考学习”了。接下来只剩最后一件事交付物的秩序和边界。5. 交付包的内容编排与后续扩展一套完整的交付包包含四类内容源码脚本、参数输入文件、说明文档、结果示例图表。源码和参数是核心资产说明书是让这套资产能被别人使用的钥匙示例图表是校验结果的证据。5.1 说明书的写作框架说明书我建议按“五段式”组织。第一段写项目目标和适用范围明确“等速、平路、无坡、无风”这几个前提第二段写数学模型把第 2 节里的公式全部列出来符号表必须完整每个符号都要带单位第三段写代码结构用目录树说明每个文件干什么、哪个是入口脚本、哪些是函数第四段写使用步骤从打开哪个文件、改哪些参数、运行哪个脚本到输出哪张图、保存哪份结果一步步列清楚第五段写已知简化假设和后续改进方向。容易被低估的是结果示例图表。跑完主脚本后自动生成的速度—续航曲线、不同速度点的功率对比表这些既是对账的凭据也能直接用于课程报告或工作汇报。交付包面向学习参考可视化越完整参考价值越高。5.2 这套模型能改造成什么等速仿真虽然聚焦但天然是扩展的起点。我梳理过一条很顺的升级路径。第一步把恒速输入换成速度序列。把综合工况循环的数据按时间步长拆成小段每个小段内仍用等速阻力公式算瞬时功率逐段累加得到总能耗。这只需要保留第 3.2 节的递推循环把单一 v_kmh 换成速度数组即可。第二步加入再生制动。减速段不需要电池向外供能反而有一部分制动能量回馈。最简单处理是给减速段加一个回馈系数比如把减速能量乘 0.6 回充到电池。加完之后续航会比纯等速贴近真实但也必须清楚每 1% 的回馈系数假设都会影响个位数百分比的续航结论假设值必须标注。第三步引入电池温度效应。低温下可用电量下降、内阻增大续航衰减明显。每隔 10°C 对容量和内阻做一次线性修正代码主结构不用动就能把冬季续航的敏感性体现出来。做到这三步原本的等速小工具就长成了微型整车能耗仿真框架。每次新增功能时都回头用第 4 节的手算校验法跑一遍回归确保老的等速功能没有被新逻辑破坏。这就是一套能用于教学、能用于工程预研、还能横向扩展的仿真包的完整成长路线。我整理这套包时最大的体会是续驶里程仿真难的部分从来不在 Matlab 语法而在于把物理链条理清、把单位盯死、把假设说明白。只要这三件事做到位哪怕模型简化到只剩常系数和直线公式算出来的结果依然能站稳脚跟。等速工况看似最基础但它恰恰是验证这套工程纪律最合适的试验场。