PHEV仿真模型开发：比亚迪唐DM混动系统建模实践

1. 项目概述：PHEV仿真模型开发背景与价值

在新能源汽车研发领域，建模仿真技术已成为缩短开发周期、降低测试成本的核心手段。以比亚迪唐DM为代表的P2-P4三擎四驱混动系统，凭借其复杂的动力耦合机制和能量管理策略，对仿真建模提出了极高要求。本项目基于Matlab/Simulink平台，完整复现了该混动架构的整车控制逻辑与动态特性。

这个仿真模型的价值主要体现在三个方面：首先，通过模块化建模可独立验证各子系统（如HCU控制策略、变速箱换挡逻辑）的功能正确性；其次，支持在虚拟环境中进行极端工况测试（如电机失效时的扭矩补偿策略）；最后，为实车控制参数标定提供先验数据，减少台架试验次数。根据实测对比，该模型在NEDC工况下的燃油消耗率仿真误差小于3.2%，动力电池SOC预测偏差控制在5%以内。

提示：完整模型包含27个核心子系统模块，需Matlab R2020b及以上版本运行，仿真步长建议设置为1ms以保证数值稳定性

2. 模型架构设计与关键组件解析

2.1 三擎四驱动力拓扑结构

比亚迪唐DM采用的P2-P4混联架构具有典型代表性：

• P2位置：ISG电机（30kW）通过离合器与发动机曲轴直连，实现启停、扭矩辅助和发电功能
• P4位置：后桥双驱动电机（2×180kW）独立驱动，构成电子四驱系统
• 动力耦合点：AMT变速箱输出轴通过分动箱将动力分配至前后轴

在Simulink中，我们使用Simscape Driveline库构建该拓扑结构。特别需要注意各旋转部件的惯性参数设置，例如：

matlab复制% 发动机转动惯量设置
engine.Inertia = 0.25; % kg·m²
% ISG电机转子惯量
ISG_rotor.Inertia = 0.018; 

2.2 整车HCU控制架构

混合动力控制单元(HCU)采用分层式设计：

1. 上层策略层

   • 驾驶意图识别（基于加速踏板开度变化率）
   • 工作模式决策（EV/HEV/Charge等）
   • 全局能量分配算法

2. 下层执行层

   • 发动机扭矩请求仲裁
   • 电机扭矩分配逻辑
   • 变速箱换挡时序控制

在模型中，我们使用Stateflow实现模式切换的状态机。典型状态转移条件如下表所示：

3. 核心子系统建模细节

3.1 发动机动态模型

采用均值有效压力(MEP)建模方法，关键参数包括：

• 万有特性曲线（基于台架试验数据拟合）
• 涡轮迟滞效应（一阶惯性环节模拟）
• 燃油消耗率MAP图

matlab复制function Torque = EngineModel(RPM, Throttle)
    % 基于查表法的扭矩计算
    persistent TorqueMap;
    if isempty(TorqueMap)
        load('EngineData.mat','TorqueMap'); 
    end
    Torque = interp2(TorqueMap, RPM, Throttle);
end

注意：发动机冷启动工况需额外考虑水温对摩擦扭矩的影响，建议添加温度修正系数

3.2 驱动电机效率优化

永磁同步电机采用场定向控制(FOC)策略，在模型中需要：

1. 建立dq轴坐标系下的电压方程
2. 实现MTPA（最大转矩电流比）算法
3. 考虑弱磁控制区特性

效率优化效果对比如下：

3.3 电池管理系统建模

锂离子电池组采用二阶RC等效电路模型，需重点关注：

• SOC估算（安时积分+EKF修正）
• 温度场分布模型
• 健康状态(SOH)衰减算法

电池参数辨识流程：

1. 脉冲放电实验获取OCV-SOC曲线
2. 最小二乘法拟合R0、Rp、Cp参数
3. 循环老化试验建立容量衰减模型

4. 整车集成与仿真验证

4.1 模型接口规范

为保证各子系统协同工作，必须统一信号接口：

• 物理量单位制（Nm, rad/s, kW等）
• 总线命名规则（如HCU_CAN1、BMS_CAN2）
• 采样时间同步机制

典型信号定义示例：

matlab复制% 发动机控制总线定义
EngineCmd = Simulink.Bus;
EngineCmd.Elements(1).Name = 'TorqueReq';
EngineCmd.Elements(1).DataType = 'single';
EngineCmd.Elements(2).Name = 'RPM_Target';
EngineCmd.Elements(2).DataType = 'uint16';

4.2 典型工况测试结果

在WLTC循环工况下，模型输出与实车数据对比如下：

4.3 硬件在环(HIL)测试部署

将模型部署到dSPACE SCALEXIO系统时需注意：

1. 固定步长离散化处理连续模块
2. 优化函数调用层次减少延迟
3. 关键信号添加Scope监控节点

5. 工程经验与问题排查

5.1 常见数值稳定性问题

• 代数环问题：在扭矩分配逻辑中出现的代数环可通过引入单位延迟(z^-1)解决
• 过零检测失效：离合器结合过程需调整Simulink的过零检测容差
• 刚性系统振荡：传动系模型使用ode23tb求解器表现更优

5.2 参数敏感性分析

通过Morris筛选法识别出对油耗影响最大的三个参数：

1. 发动机怠速转速（±50rpm导致油耗变化2.1%）
2. 电机峰值效率点（1%效率变化影响续航1.8km）
3. 电池内阻温度系数（10%偏差导致SOC估算误差3.5%）

5.3 模型加速技巧

• 对查表模块启用"Lookup Table Optimizer"
• 将MATLAB Function块转换为C-MEX S函数
• 使用并行计算工具箱加速参数扫描

在实际项目中，我们通过代码生成技术将核心控制算法部署到TC297芯片，相比解释执行速度提升17倍。这个过程中发现Simulink的ERT代码生成模板需要针对Autosar架构进行以下调整：

1. 将全局变量封装到RTE接口
2. 添加BSW服务调用钩子函数
3. 配置多核任务调度映射