1. 串联式增程混合动力系统建模概述
在新能源汽车领域，串联增程式混合动力系统因其独特的能量传递路径而备受关注。这种构型中，发动机不直接驱动车轮，而是作为发电机组的动力源，通过电能转换实现动力输出。Simulink作为多域仿真平台，能够完美呈现这种机电耦合系统的动态特性。
我去年参与开发的某商用增程式电动车项目，正是基于这种架构。实际开发中遇到最棘手的问题是动力源切换时的扭矩波动控制，这促使我深入研究了整套建模方法。下面分享的模型包含三个关键子系统：反映中国典型道路特征的工况模型、模拟人类驾驶行为的驾驶员模型，以及最核心的串联式动力总成模型。
2. 模型架构设计与核心模块解析
2.1 整车动力系统拓扑设计
串联式架构的核心特征是"发动机-发电机-电池-电机"的单一能量路径。在我们的模型中：
- 1.8L阿特金森循环发动机通过转速闭环控制驱动永磁同步发电机
- 锂离子电池组作为能量缓冲单元，采用Thevenin等效电路模型
- 驱动电机采用基于效率MAP的前馈+反馈扭矩控制
特别要注意发电机与驱动电机的功率匹配。我们通过公式P_gen = η·P_engine（η取0.82-0.85）确定发电机额定功率，再根据车辆最大需求功率确定电池放电功率补充能力。
2.2 工况模型构建要点
中国典型工况（CLTC-P）的建模需要特别注意：
matlab复制% 速度-时间曲线插值处理
time = 0:0.1:1800; % 30分钟工况
speed = interp1(StandardTime, StandardSpeed, time, 'pchip');
accel = [0 diff(speed)./diff(time)]; % 计算加速度
实际项目中我们发现原始工况数据存在加速度突变问题，通过三次样条插值和移动平均滤波处理后，车辆纵向动力学响应更符合真实情况。
2.3 驾驶员模型开发技巧
采用预瞄-跟随驾驶员模型时，关键参数设置经验：
- 预瞄距离建议取车速的1.2-1.5倍（单位：米）
- 转向延迟时间常数设为0.3-0.5秒
- 油门/制动踏板映射曲线需要实测标定
我们在江淮iEV7S实车采集的驾驶员操作数据表明，熟练驾驶员的速度跟踪误差标准差不超过0.3m/s，这个指标可作为模型验证基准。
3. 关键子系统建模与参数标定
3.1 发动机-发电机组建模
使用准静态建模方法时要注意：
matlab复制% 发动机万有特性数据插值
engine_map = load('Atkinson_1.8L.mat');
fuel_consumption = griddedInterpolant(...
engine_map.RPM, engine_map.Torque, engine_map.BSFC);
实测中发现，发动机冷启动阶段的燃油消耗量比稳态高18%-22%，需要在模型中添加温度修正系数。我们开发的动态补偿算法使燃油经济性预测误差从7.2%降至3.5%。
3.2 电池管理系统建模
锂离子电池的二阶RC模型参数辨识流程：
- 通过HPPC测试获取OCV-SOC曲线
- 用最小二乘法拟合脉冲充放电数据
- 考虑温度影响的Arrhenius方程修正
某三元锂电池包的测试数据显示，-10℃时内阻会增加65%，这在低温工况仿真中必须体现。我们开发的变参数模型在NEDC工况下SOC预测误差<1.5%。
3.3 驱动电机控制策略
永磁同步电机的效率优化控制要点：
- 基于d-q轴电流解耦的MTPA控制
- 弱磁控制转折点设为基速的1.3倍
- 考虑逆变器损耗的全局效率MAP
实测数据对比显示，采用动态效率优化策略可使城市工况能耗降低6-8%。下面是一个典型的效率查表实现：
matlab复制% 电机效率二维插值表
eff_map = griddedInterpolant(...
motor_speed, motor_torque, motor_efficiency);
current_eff = eff_map(actual_speed, command_torque);
4. 模型集成与仿真验证
4.1 系统级耦合问题处理
各子系统集成时常见问题及解决方案：
| 问题现象 | 根本原因 | 解决方法 |
|---|---|---|
| 仿真步长突变 | 变步长算法遇到刚度变化 | 固定步长+代数环检测 |
| SOC振荡 | 电池模型时间常数过小 | 增大RC网络电容值 |
| 扭矩响应延迟 | 执行器动力学未建模 | 增加一阶惯性环节 |
我们采用0.01秒固定步长配合Rosenbrock求解器，在保证精度的同时将仿真速度提升40%。
4.2 典型工况验证结果
CLTC-P工况下的关键指标对比：
| 参数 | 仿真值 | 实测值 | 误差 |
|---|---|---|---|
| 百公里电耗 | 15.2kWh | 15.8kWh | 3.8% |
| SOC变化量 | 28.5% | 27.9% | 2.1% |
| 最高电机温度 | 82℃ | 85℃ | 3.5% |
特别要注意的是，空调负载会使实际电耗增加10-15%，建议在模型中增加可选附件负载模块。
5. 工程应用中的经验技巧
5.1 模型加速方法
大型模型运行缓慢时的优化策略：
- 将查表数据转换为Simulink.LookupTable对象
- 对非关键子系统使用FMU导出
- 启用并行计算工具箱
某包含200+模块的完整车辆模型，经过优化后单次WLTC工况仿真时间从38分钟缩短至6分钟。
5.2 参数敏感性分析
采用Morris筛选法确定关键参数：
- 确定35个待分析参数及其变化范围
- 生成100组采样点进行全局仿真
- 计算各参数的μ*（平均影响度）和σ（交互作用）
分析发现电池内阻和电机峰值效率对能耗影响最为显著，这指导我们优化了测试资源配置。
5.3 硬件在环测试衔接
模型向HIL测试过渡的注意事项：
- 将连续控制器离散化时，采样时间不大于5ms
- 替换理想传感器模型为带噪声的版本
- 增加通信延迟补偿模块
我们在转鼓试验台上对比发现，经过适配的模型与实车CAN信号吻合度达到92%以上。