EPS与二自由度转向模型的Simulink实现与工程实践

1. 项目概述：当EPS遇上二自由度转向模型

在车辆动力学仿真领域，EPS（电动助力转向系统）与二自由度转向模型的结合一直是底盘开发的核心课题。最近我在搭建这个经典模型的Simulink S函数实现时，发现现有文献大多停留在理论推导，缺少从零开始的完整工程实践指南。本文将分享如何用MATLAB 2025最新语法构建这个联合仿真模型，过程中会特别说明EPS助力特性曲线拟合、轮胎侧偏刚度参数化处理等工程细节。

这个模型的价值在于：它既能用于EPS控制策略的快速验证（比如不同车速下的助力曲线设计），又能评估整车转向特性（不足/过多转向趋势）。对于从事底盘电控的工程师，掌握这套建模方法可以节省约40%的算法开发时间。下面我会从模型原理拆解开始，逐步展示S函数的具体实现。

2. 模型原理与关键参数

2.1 EPS系统建模要点

EPS模型的核心是建立方向盘转矩到转向电机电流的映射关系。实测数据显示，典型的EPS特性曲线包含三个关键区段：

• 死区（±0.5Nm）：忽略微小转矩以避免误触发
• 线性助力区（0.5-4Nm）：助力增益随车速增大而减小
• 饱和区（>4Nm）：限幅保护电机

在MATLAB中，我采用分段多项式拟合实现这个非线性特性：

matlab复制function current = eps_curve(torque, v)
    % 死区处理
    if abs(torque) < 0.5
        current = 0;
    % 线性区
    elseif abs(torque) <=4
        current = torque * (12 - 0.2*v); 
    % 饱和区
    else
        current = sign(torque) * (48 - 0.8*v);
    end
end

2.2 二自由度转向模型推导

整车模型采用经典的自行车模型，核心方程包括：

code复制m*v*(dβ/dt + r) = Fyf + Fyr
Iz*dr/dt = a*Fyf - b*Fyr

其中轮胎侧向力采用线性假设：

code复制Fyf = -Cf*αf = -Cf*(δ - β - a*r/v)
Fyr = -Cr*αr = -Cr*(-β + b*r/v)

关键提示：实际仿真中需添加转向系统刚度系数Ks（通常取800-1200N·m/rad），否则会导致前轮转角δ与方向盘输入θ的线性关系失真。

3. S函数实现详解

3.1 模块接口设计

在Level-2 MATLAB S函数中，需要明确定义：

matlab复制function sys=mdlDerivatives(t,x,u)
    % 输入u: [方向盘转矩θ, 车速v]
    % 状态x: [横摆角速度r, 侧偏角β, 前轮转角δ]
    % 输出y: [横摆角速度r, 侧向加速度ay]
    
    % EPS电流计算
    i_motor = eps_curve(u(1), u(2)); 
    
    % 转向系统动力学
    delta_dot = (i_motor*Kt - Ks*x(3))/Js;
    
    % 整车动力学
    [beta_dot, r_dot] = vehicle_dynamics(x(2),x(1),x(3),u(2));
    
    sys = [r_dot; beta_dot; delta_dot]; 
end

3.2 参数初始化技巧

在mdlInitializeSizes中设置可调参数时，建议采用结构体组织：

matlab复制function params = init_params()
    params.m = 1850;   % 整车质量(kg)
    params.Iz = 3500;  % 横摆惯量(kg·m²)
    params.Cf = 80000; % 前轮侧偏刚度(N/rad)
    params.Cr = 95000; % 后轮侧偏刚度(N/rad)
    params.Kt = 0.12;  % 电机转矩常数(N·m/A)
end

避坑指南：轮胎刚度单位易混淆，有些文献用N/deg需要转换。建议在注释中明确标注单位。

4. 仿真验证与问题排查

4.1 典型测试场景设置

建议分三个阶段验证：

1. 静态验证（v=0）：检查方向盘转矩到前轮转角的传动比
2. 阶跃转向输入：观察横摆角速度响应超调量
3. 正弦扫频测试：分析相位滞后与幅值衰减

4.2 常见异常与解决方案

5. 模型应用扩展

5.1 与ADAS系统联合仿真

将模型输出接入车道保持控制器：

matlab复制function delta_cmd = lka_control(y, ref)
    % y: [r, ay]
    % ref: 期望路径曲率
    Kp = 0.8; Ki = 0.2;
    persistent integral_err;
    
    err = ref - y(1)/v;
    integral_err = integral_err + err*Ts;
    
    delta_cmd = Kp*err + Ki*integral_err;
end

5.2 参数辨识实践

基于实车数据反求轮胎刚度：

1. 采集蛇形试验的r/ay信号
2. 构建优化问题：

matlab复制function err = objective(C)
    sim_out = run_model(C);
    err = norm(sim_out.r - measured_r);
end
C_opt = fminsearch(@objective, [80000,95000]);

6. 工程经验沉淀

在多个量产项目实践中，我总结了几个关键认知：

1. 冬季测试时轮胎刚度会降低约15%，需在模型中加入温度补偿系数
2. EPS死区设置过小会导致高速工况"抢方向盘"现象
3. 横摆角速度传感器噪声会显著影响控制性能，建议在模型输出端添加2Hz低通滤波

模型最终精度验证标准：在80km/h阶跃转向工况下，横摆角速度稳态误差应小于5%，峰值响应时间误差在0.1秒内。达到这个水平后，模型即可用于大多数控制策略开发场景。