双三相PMSM缺相容错控制策略与Simulink实现-嵌云网-嵌入式AI开发资源站

双三相PMSM缺相容错控制策略与Simulink实现

科学声音

1. 双三相PMSM缺相容错控制概述

双三相永磁同步电机（Dual Three-Phase PMSM）因其高功率密度、低转矩脉动和容错能力强的特点，在航空航天、电动汽车等高可靠性应用领域备受关注。与传统三相电机相比，双三相结构通过两组三相绕组（通常相位差30°）实现了功率分流和冗余控制。当发生单相开路故障时，系统仍能通过剩余五相绕组继续运行，但需要特殊的容错控制策略来维持性能。

缺相运行的核心挑战在于：

电流空间矢量分布不对称导致的谐波问题
转矩输出能力下降（最大输出转矩约为正常状态的69.4%）
定子铜耗增加和效率降低
控制系统需要实时重构

Simulink作为电机控制系统设计的标准工具，提供了从算法开发到硬件实现的完整仿真环境。其模块化设计特别适合实现复杂的容错控制算法，包括：

故障检测与诊断模块
坐标变换与解耦控制模块
虚拟矢量生成模块
在线参数辨识模块

2. 缺相故障建模与特性分析

2.1 双三相PMSM健康状态模型

在自然坐标系下，双三相PMSM的电压方程可表示为：

code复制[V]6×1 = [R]6×6[I]6×1 + d/dt[L]6×6[I]6×1 + ωψm[E]6×1

其中R为电阻矩阵，L为电感矩阵，ψm为永磁体磁链，E为反电动势系数矩阵。通过矢量空间解耦变换（VSD），可将系统分解为：

α-β子空间：产生有效转矩
z1-z2子空间：引起谐波损耗
零序子空间：中性点隔离时为零

2.2 单相缺相故障特性

当F相发生开路故障时，系统表现出：

电流特性：

剩余五相电流幅值不均衡
α-β子空间电流轨迹畸变
z1-z2子空间出现谐波电流（THD可达141.8%）

转矩特性：

平均转矩下降约30.6%
出现二次谐波转矩脉动（达直流分量的40.76%）

损耗特性：

铜耗直流量增加38.09%
效率降低11.85%

实验数据表明，传统控制策略在缺相情况下会导致：

相电流THD从5.6%恶化至141.8%
转矩脉动从1N·m增至1.7N·m

3. 混合容错控制策略设计

3.1 控制策略分区原理

根据输出转矩范围采用不同控制方式：

低转矩区（0-0.542 p.u.）：最小损耗（ML）策略
高转矩区（0.542-0.694 p.u.）：最大转矩（MT）策略
过渡区采用混合控制策略

关键参数说明：p.u.表示标幺值，以额定转矩Te为基准，1 p.u.=Te

3.2 ML控制策略实现

目标函数：

code复制J_ML = Σ(i_k^2), k∈{A,B,C,D,E}

最小化各相电流幅值平方和，降低铜耗

实现步骤：

通过拉格朗日乘数法求解约束优化问题
设置z1-z2子空间参考电流iz1=iz2=0
得到各相电流分配系数：

code复制i_k = a_k·Iα + b_k·Iβ

其中系数a_k、b_k通过对称分量法确定

3.3 MT控制策略实现

目标函数：

code复制J_MT = max(|i_k|)

限制最大相电流幅值，提高转矩输出能力

实现步骤：

使用Matlab优化工具箱求解非线性约束问题
获得各相电流最优分配方案
计算z1-z2子空间参考电流

3.4 混合控制策略

定义混合系数λ∈[0,1]：

code复制T = 0.542λ + 0.694(1-λ)

当λ=0.5时，输出转矩0.618 p.u.，此时：

ML策略贡献50%电流分量
MT策略贡献50%电流分量

动态调整规则：

检测当前总输出转矩T
计算λ=(0.694-T)/(0.694-0.542)
按比例混合两种策略的输出

4. Simulink实现细节

4.1 系统架构设计

完整仿真模型包含：

电机本体模块：

实现双三相PMSM的六相模型
集成故障注入接口

故障检测模块：

基于电流频谱分析的实时诊断
故障相定位逻辑

混合控制器：

转矩需求分区逻辑
ML/MT策略切换机制
电流参考生成

空间矢量调制：

虚拟矢量重构算法
12扇区SVPWM实现

4.2 关键模块参数设置

坐标变换模块：

matlab复制% 健康状态变换矩阵
T6s = (2/3)*[1, -0.5, -0.5, sqrt(3)/2, -sqrt(3)/2, 0;
             0, sqrt(3)/2, -sqrt(3)/2, 0.5, 0.5, -1;
             1, -0.5, -0.5, -sqrt(3)/2, sqrt(3)/2, 0;
             0, -sqrt(3)/2, sqrt(3)/2, 0.5, 0.5, -1;
             1, 1, 1, 0, 0, 0];
             
% 缺相状态修正矩阵（F相开路）
T5s = T6s(:,[1,2,3,4,5]); % 移除F相列

电流分配优化模块：

matlab复制function [a_k, b_k] = currentAllocation(phase)
    % ML策略系数表
    persistent ml_coeff = [
        0.897, 0.441;   % A相
        0.276, 0.851;   % B相
        -0.723, 0.851;  % C相
        -0.897, 0.441;  % D相
        -0.276, -0.851]; % E相
    
    % MT策略系数表
    persistent mt_coeff = [
        0.945, 0.328;
        0.145, 0.989;
        -0.800, 0.600;
        -0.945, 0.328;
        -0.145, -0.989];
    
    switch control_mode
        case 'ML'
            a_k = ml_coeff(phase,1);
            b_k = ml_coeff(phase,2);
        case 'MT'
            a_k = mt_coeff(phase,1);
            b_k = mt_coeff(phase,2);
    end
end

4.3 虚拟矢量重构技术

等幅虚拟矢量：

幅值统一为0.277Udc
零矢量占比约47.4%
适合低速轻载（THD=7.8%）

最大幅值虚拟矢量：

消除零矢量
最大幅值达0.527Udc
适合高速重载（THD=15.7%）

实现示例：

matlab复制% 虚拟矢量合成（以V1为例）
V_virtual = 0.0681*V16 + 0.4593*V18 + 0.4228*V26 + 0.0498*V27;

5. 仿真结果分析

5.1 稳态性能对比

指标	健康状态	传统容错	混合策略
转矩脉动	1.0 N·m	1.7 N·m	1.2 N·m
相电流THD	5.6%	141.8%	7.8%
最大输出转矩	1.0 p.u.	0.542 p.u.	0.694 p.u.
铜耗增加	-	42.3%	19.8%

5.2 动态响应测试

加速测试（100→300rpm）：

转速上升时间：0.15s
超调量：<5%
转矩跟踪误差：±0.3N·m

负载突变（3→7N·m）：

转矩恢复时间：0.08s
转速波动：<2rpm
电流调节时间：0.05s

5.3 关键波形展示

混合控制下相电流：

各相幅值：A(4.2A), B(3.8A), C(4.5A), D(4.1A), E(3.9A)
相位差：72°±3°（理想72°）

转矩-铜耗曲线：

0.5 p.u.转矩时，铜耗降低23.5%
0.6 p.u.转矩时，定子损耗减少19%

6. 工程实现建议

参数敏感性处理：

电感变化±20%时，需在线更新变换矩阵
电阻温漂每10°C需重新标定

实时性优化：

电流环控制周期≤100μs
速度环控制周期≤1ms
故障检测延时<2ms

实验平台搭建要点：

采用Infineon XE164微控制器
配置12位ADC采样（采样率≥50kHz）
死区时间设置为2μs
母线电压检测精度±0.5%

实际调试中发现，采用如下配置可改善性能：

PWM开关频率：10kHz
电流滤波截止频率：2kHz
速度观测器带宽：100Hz

在TI LAUNCHXL-F28379D开发板上实测表明：

算法执行时间：85μs
代码占用Flash：56KB
RAM消耗：12KB

7. 扩展应用与改进方向

多故障模式扩展：

可适配双相短路故障
支持相邻两相开路情况

先进控制算法融合：

模型预测控制（MPC）
滑模观测器（SMO）
自适应模糊PID

硬件在环测试：

使用Speedgoat实时目标机
仿真步长设置为50μs
添加±5%参数扰动测试鲁棒性

一个实用的改进技巧是：在DSP中预存多种故障模式的变换矩阵，通过查表法快速切换控制策略，可将故障响应时间缩短至1ms以内。