1. 电励磁同步电机仿真概述
电励磁同步电机作为工业领域的重要动力设备，其动态过程仿真一直是电机控制领域的研究热点。这次我们要实现的Matlab仿真，将完整呈现电机从启动到运行再到制动的全过程。不同于永磁同步电机，电励磁同步电机通过调节转子励磁电流可以实现功率因数的灵活控制，这使得它在大型工业设备和电力系统中具有独特优势。
在实际工程中，我们经常需要预测电机在不同工况下的性能表现。比如在轧钢机应用中，电机需要频繁启停和变速运行；在船舶推进系统中，则更关注电机的动态响应和制动性能。通过Matlab仿真，我们可以在不接触实际设备的情况下，全面评估电机的运行特性，这对系统设计和故障诊断都大有裨益。
2. 仿真模型构建
2.1 电机数学模型建立
电励磁同步电机的动态特性可以通过派克变换后的d-q轴方程来描述。在建模时，我们需要考虑以下关键方程：
定子电压方程：
code复制v_d = R_s*i_d + dψ_d/dt - ω_r*ψ_q
v_q = R_s*i_q + dψ_q/dt + ω_r*ψ_d
转子电压方程：
code复制v_f = R_f*i_f + dψ_f/dt
磁链方程：
code复制ψ_d = L_d*i_d + L_md*i_f
ψ_q = L_q*i_q
ψ_f = L_md*i_d + L_f*i_f
电磁转矩方程：
code复制T_e = (3/2)*(P/2)*(ψ_d*i_q - ψ_q*i_d)
这些方程构成了我们仿真模型的基础。在Matlab/Simulink中，我们可以通过S-Function或者直接使用Simscape Electrical库中的同步电机模块来实现这些方程。
2.2 仿真参数设置
一个典型的10kW电励磁同步电机参数如下表所示：
| 参数名称 | 符号 | 数值 | 单位 |
|---|---|---|---|
| 额定功率 | P_N | 10 | kW |
| 定子电阻 | R_s | 0.2 | Ω |
| d轴电感 | L_d | 8.5 | mH |
| q轴电感 | L_q | 6.2 | mH |
| 励磁电阻 | R_f | 0.5 | Ω |
| 励磁电感 | L_f | 150 | mH |
| 极对数 | P | 4 | - |
| 转动惯量 | J | 0.08 | kg·m² |
在Simulink中设置这些参数时，需要注意单位的一致性。特别是电感参数，Simulink默认使用H（亨利）作为单位，而电机参数表通常给出的是mH（毫亨）。
3. 启动过程仿真实现
3.1 异步启动建模
电励磁同步电机的异步启动过程需要特别关注以下几个关键点：
- 鼠笼绕组合成电阻设置：
code复制R_b = 3*R_s + R_ext
其中R_ext为外接启动电阻，用于限制启动电流。
- 启动转矩特性：
code复制T_start = (3*V^2*R_r)/(ω_s*(R_s+R_r)^2 + (X_s+X_r)^2)
式中V为相电压，R_r为转子等效电阻，X_s和X_r分别为定转子漏抗。
在Simulink中，我们可以通过以下步骤实现异步启动模型：
- 使用Three-Phase Programmable Voltage Source模块作为电源
- 设置Synchronous Machine pu Standard模块的初始参数
- 添加Asynchronous Start Torque模块模拟鼠笼绕组效应
- 配置适当的机械负载模型
3.2 牵入同步控制
当转子转速达到95%同步转速时，需要投入励磁。这个切换点的判断逻辑可以通过以下Matlab函数实现：
matlab复制function [firing] = sync_condition(omega_act, omega_sync, threshold)
slip = (omega_sync - omega_act)/omega_sync;
if slip <= threshold
firing = 1; % 投入励磁
else
firing = 0; % 保持异步运行
end
end
励磁投入的时序控制非常关键，过早会导致启动失败，过晚则可能引起电流冲击。在实际仿真中，我们通常采用PID控制器来调节励磁电流的上升斜率：
code复制I_f(t) = K_p*e(t) + K_i∫e(t)dt + K_d*de(t)/dt
其中e(t)为转速误差信号。
4. 运行阶段仿真分析
4.1 稳态运行特性
在同步运行阶段，我们需要关注以下关键指标：
- 功率因数特性：
code复制cosφ = P/(S) = P/(√(P^2 + Q^2))
通过调节励磁电流I_f，可以改变电机的功率因数。
- V型曲线绘制：
保持负载转矩不变，改变励磁电流，记录定子电流变化，可以得到典型的V型曲线。在Simulink中，可以使用Sweep模块自动完成这个扫描过程。
4.2 动态响应测试
为评估电机的动态性能，我们可以模拟以下工况：
-
突加负载测试：
在t=1s时突然施加50%额定负载，观察转速和电流的瞬态响应。 -
励磁调节测试：
在t=2s时阶跃改变励磁电流20%，观察无功功率的变化。
这些测试可以帮助我们评估控制系统的动态性能指标，如调节时间、超调量等。
5. 能耗制动仿真实现
5.1 制动电路建模
能耗制动的关键是在定子侧接入直流电源。在Simulink中，我们可以通过以下方式实现：
- 使用Three-Phase Breaker模块切断交流电源
- 通过Controlled Voltage Source接入直流电压
- 设置适当的制动电阻值：
code复制R_br = V_dc^2 / P_br
其中P_br为预期的制动功率。
5.2 制动过程优化
为提高制动效率，我们可以采用以下策略：
-
分级制动：
根据转速分段调整直流电压大小，高速段用较高电压，低速段降低电压。 -
电流闭环控制：
code复制V_dc = K*(I_ref - I_act)
保持制动电流在安全范围内。
制动时间的理论估算：
code复制t_br = J*Δω / T_br
其中J为转动惯量，Δω为转速变化量，T_br为平均制动转矩。
6. 仿真结果分析
6.1 启动过程波形
典型的启动过程波形应包括：
- 转速从0加速到同步转速的曲线
- 定子电流的幅值变化
- 励磁电流的投入时序
- 电磁转矩的波动情况
通过FFT分析可以观察到启动过程中电流谐波成分的变化。
6.2 制动性能评估
制动性能的主要评价指标：
- 制动时间：从开始制动到完全停止的时间
- 能量损耗：制动过程中电阻消耗的总能量
- 电流峰值：制动过程中的最大电流值
这些数据可以帮助我们优化制动电阻的参数选择。
7. 实际应用中的注意事项
- 参数敏感性分析：
电机的性能对某些参数特别敏感，如：
- 定转子互感L_md：影响牵入同步能力
- 转动惯量J：影响启动和制动时间
- 电阻值：影响发热和效率
建议通过参数扫描分析这些敏感参数的影响。
- 数字仿真中的数值问题：
- 使用适当的求解器（如ode23tb）
- 设置合理的最大步长（通常为1e-4s）
- 注意代数环问题，必要时插入Memory模块
- 模型验证技巧：
- 与空载试验数据对比
- 检查功率平衡关系
- 验证稳态工作点的合理性
8. 仿真模型扩展应用
完整的电励磁同步电机模型还可以用于以下研究：
- 故障仿真：
- 定子绕组匝间短路
- 励磁绕组开路
- 轴承故障模拟
- 先进控制算法测试：
- 矢量控制实现
- 直接转矩控制
- 自适应控制策略
- 系统级仿真：
- 多电机协同控制
- 电网互动研究
- 新能源发电应用
对于想进一步研究的同学，可以尝试加入温度场模型，研究不同工况下的电机温升特性。这需要耦合电磁模型和热模型，可以使用Simscape Multibody来实现多物理场仿真。