电励磁同步电机EESM和永磁同步电机相比转子磁场完全可调因此在需要宽范围弱磁、频繁起动制动或者对磁场控制精度要求高的场合它一直是工程上绕不开的方案。但“磁场可调”这四个字放在Matlab/Simulink仿真里意味着三件事都会变复杂启动过程怎么把转子平滑地带起来运行过程怎么把磁场和转矩解耦控制停机时怎么把转子里的动能可控地损耗掉。我最近把“启动运行能耗制动”三个过程放在同一个仿真模型里完整跑了一遍从异步起动到同步运行再到断电制动全程用Matlab/Simulink搭了一套可复用的框架。这篇文章就把模型结构、参数整定、切换逻辑和实际调试中踩过的坑全部写出来给正在做同步电机仿真尤其是准备做瞬态过程研究或者控制器开发的朋友做个参考。这套仿真的价值不在于某一个阶段的波形多漂亮而在于它把一台电励磁同步电机从零转速到额定转速、带载稳定运行、再到制动停机的完整生命周期拉通了一遍。你拿到手里可以直接改参数也能把它当成后续做无速度传感器、弱磁控制、故障诊断的基座模型。我尽量把每一步选择背后的理由说清楚而不是只给你一堆结论。1. 电励磁同步电机的数学模型三阶段仿真的地基不管是启动、运行还是制动Simulink模型的核心都不是那几个封装好的电机模块而是电机本体和负载的数学模型。电励磁同步电机比异步电机多了一条独立的转子励磁回路比永磁同步电机多了一个可以控制的磁场源所以dq轴方程里必须把励磁绕组的动态和它与定子d轴磁链的耦合单独写出来。1.1 从物理结构到dq轴方程定子侧是三相交流绕组转子侧有励磁绕组和阻尼绕组。在这里我先忽略阻尼绕组的影响只保留励磁绕组因为变频启动方案下阻尼绕组作用有限而能耗制动阶段阻尼绕组主要影响初始瞬态对整体趋势影响不大。在转子磁场定向的dq旋转坐标系下定子电压方程写成$$u_{sd} R_s i_{sd} \frac{d\psi_{sd}}{dt} - \omega_e \psi_{sq}$$$$u_{sq} R_s i_{sq} \frac{d\psi_{sq}}{dt} \omega_e \psi_{sd}$$磁链方程$$\psi_{sd} L_d i_{sd} M_{sf} i_f$$$$\psi_{sq} L_q i_{sq}$$$$\psi_f L_f i_f M_{sf} i_{sd}$$励磁绕组的电压方程$$u_f R_f i_f \frac{d\psi_f}{dt}$$电磁转矩$$T_e \frac{3}{2} p \left( \psi_{sd} i_{sq} - \psi_{sq} i_{sd} \right)$$再加上运动方程$$J \frac{d\omega_m}{dt} T_e - T_L - B \omega_m$$这套方程是后面所有Simulink模块的数学依据。这里最值得注意的是$\psi_{sd}$里出现了$M_{sf} i_f$这一项也就是说转子励磁电流会直接改变定子d轴的磁链水平。你在做id0控制时转矩表达式简化为$T_e \frac{3}{2} p M_{sf} i_f i_{sq}$磁场和转矩的控制通道就彻底分开了。1.2 励磁绕组动态为什么不能忽略很多入门仿真喜欢把励磁电流当成一个恒定值直接给进去这在稳态运行分析里问题不大但是在启动和制动过程里会出大问题。励磁绕组是一个大电感小电阻回路时间常数通常有几十毫秒甚至上百毫秒。启动瞬间如果励磁电流突变转子磁链不能突变它和定子磁链之间的耦合会产生剧烈的瞬态振荡表现出来就是转矩波动、电流尖峰。我在这套模型里把励磁回路完整建了出来而不是用常数代替。Simulink里可以直接用受控电流源加电感电阻来搭但更稳妥的做法是直接用基本的Simscape电力电子元件或者自己写一个S-Function实现上述方程。我的做法是分模块搭建dq轴的磁链积分模块、坐标变换模块、转矩计算模块、运动方程模块励磁绕组单独一个子系统。这样每个环节都能单独看波形定位问题也方便。1.3 仿真参数选用及标幺问题下面是我用的一组10kW机组估算参数不是某个具体型号但量级是真实的。仿真时建议先用这组参数把框架跑通再替换成你手头电机的实际参数。参数数值说明额定功率10 kW输出机械功率额定电压380 V线电压有效值额定频率50 Hz同步频率极对数p2额定转速1500 rpm定子电阻Rs0.2 Ω75度热态值定子d轴电感Ld10 mH含互感和漏感定子q轴电感Lq8 mH凸极效应励磁电阻Rf12 Ω励磁绕组热态电阻励磁电感Lf0.25 H励磁时间常数约20ms定子励磁互感Msf0.3 H决定磁场耦合强度转动惯量J0.4 kg·m²电机加负载折算值阻尼系数B0.005 N·m·s/rad机械损耗等效额定励磁电流3 A对应额定气隙磁场关于标幺制我个人建议在模型内部用国际单位制只在控制器里换算。标幺值虽然看起来优雅但调试时一旦某个基值给错实际物理量就对不上。直接用SI单位每个变量是什么量级一目了然。等模型稳定了你想切标幺再切。2. Simulink整体架构与三阶段切换逻辑三阶段仿真最容易犯的错误是把启动、运行、制动做成三个互相独立的仿真模型然后用“结果粘贴”的方式讲故事。那样做出来的模型没有连续物理过程切换瞬间的表现全是假的。我的做法是把三个阶段放在同一个Simulink模型里通过一个状态机切换电气拓扑和控制策略。2.1 顶层模型拆分模型顶层分成四个部分电机本体模块dq方程积分、坐标反变换、三相输出变频器与定子切换模块逆变桥、交流接触器模型、制动电阻励磁回路模块可控直流源、励磁斜坡给定控制与切换逻辑模块Stateflow状态机、PI调节器、V/f生成器、制动电阻投切信号其中“定子切换模块”是这套模型和普通电机模型最大的区别。运行阶段定子接在逆变器上制动阶段定子要通过接触器切到三相对称制动电阻上。Simulink里我用三个理想开关元件模拟接触器一组开关接到逆变器一组开关接到制动电阻利用互锁逻辑保证切换瞬间不会出现两路同时接通的情况。2.2 阶段切换的判断条件与Stateflow实现Stateflow状态机里我建了三个状态STARTUP、RUN、BRAKE。切换条件不是简单地按时间点切而是按系统实际状态判断。STARTUP到RUN电机转速达到额定转速的95%且持续0.2秒。这个滞回和持续时间条件能避免转速波动导致的误切换。如果在启动过程中负载过重导致一直升不到95%转速模型就停留在STARTUP状态并给出报警这样对复现真实故障场景也有帮助。RUN到BRAKE收到外部制动指令或者运行计时器到达预设时间。我模型里用的是“运行2秒后自动进入制动”这样跑完一次仿真就能把三个阶段波型全部看到。BRAKE到停机完成转速低于10rpm此时惯性已经很小继续制动意义不大制动接触器断开仿真结束。用Stateflow而不是普通逻辑门的好处是切换过程的优先级、抖动抑制、可读性都更好。你不想引入Stateflow的话也可以用Matlab Function加上升沿触发器实现但阶段多了以后代码会变得很乱。2.3 三阶段电气拓扑的变化下面表格把三个阶段的关键差异列出来方便对照理解阶段定子连接励磁回路控制策略目标启动逆变器输出励磁给定斜坡上升V/f控制带低频补偿平滑升速到额定转速运行逆变器输出励磁电流恒定3Aid0双闭环矢量控制抗负载扰动转速稳定制动三相对称制动电阻励磁保持3A无逆变器控制电阻耗能快速可控消耗动能能耗制动阶段逆变器必须完全封锁这一点很多第一次仿真的人容易忽略。如果制动瞬间逆变器还在输出那就变成了再生制动或者反接制动能量去向完全不同。能耗制动的核心就是“定子脱离电源、转子继续励磁、定子接电阻”让电机工作在发电机状态把机械能变成电能耗在电阻上。3. 启动阶段斜坡升频升压的控制实现电励磁同步电机的启动工程上就两种主流路线。一种是异步启动转子带阻尼绕组接电网后靠异步转矩把转子拉起来接近同步速再投励牵入同步另一种是变频启动逆变器从低频低压开始逐步升频升压全程都在同步电机状态下加速。我建模时采用的是变频启动原因很直接异步启动的冲击电流和冲击转矩太大放在仿真里对数值求解很不友好而且后续做带载启动、S曲线加减速这类优化时变频启动的灵活性强太多。3.1 启动阶段为什么用V/f控制而不是矢量控制启动阶段转速很低反电动势很小电流环所需要的解耦电压基本建立不起来。矢量控制在低速段需要精确的转子位置和磁链观测一旦磁链观测初始值不对电流环就会发散。V/f控制没有电流闭环只保证电压矢量幅值和频率按比例变化低速段的鲁棒性反而好。我采用的V/f曲线是最低频率2Hz对应输出电压约50V然后按7.6V/Hz比例上升到50Hz时达到380V。这个“最低频率下的电压补偿”是为了克服定子电阻压降。如果不加补偿低频段定子电阻上的压降占比很大实际励磁电流建立不起来电机在2-5Hz区间会“爬行”转速半天不动电流却很大。3.2 Simulink里的具体搭建方式用Matlab Function生成频率斜坡和电压幅值然后经过积分得到角度再做三相正弦波生成function [freq, vref, theta] vf_generator(t, t_ramp, f_base, v_base) % t_ramp: 升频时间, 我设成1秒 % f_base: 50Hz, v_base: 380V if t 0.01 freq 2; elseif t t_ramp freq 2 (f_base - 2) * (t - 0.01) / (t_ramp - 0.01); else freq f_base; end vref v_base * (freq / f_base); % 简化的V/f比实际可加IR补偿 theta 2 * pi * integral_freq; % 角度由频率积分得到 v_a vref * sqrt(2/3) * sin(theta); v_b vref * sqrt(2/3) * sin(theta - 2*pi/3); v_c vref * sqrt(2/3) * sin(theta 2*pi/3); end注意这里的角度不能用freq直接乘时间必须对频率做积分。你直接用梯形公式累加就行Simulink里的Integrator模块比手写积分更稳。3.3 启动阶段仿真波形与判断标准我仿真的结果是0.4秒时转速到达500rpm0.8秒时到达1200rpm1秒左右跟上50Hz同步转速约1420rpm的机械转速加转差修正后进入同步。定子电流在整个启动过程中被限制在额定电流的1.6倍以内没有出现明显尖峰。这里有一个关键指标启动过程中瞬态转矩平均值应当始终为正且不能长时间低于负载转矩否则转速会卡在某个平台。如果你仿真时发现启动阶段转速上升缓慢或者电流振荡先看励磁电流是不是已经稳定建立。再强调一次励磁回路时间常数约20毫秒启动时我给的是斜坡励磁指令从0开始用100毫秒上升至3A这样避免了磁场突变引发的转矩振荡。启动结束后励磁电流稳定在3A整个启动过程才算是健康状态。4. 运行阶段id0矢量控制与转速闭环启动完成、转速进入同步范围后控制策略就要从开环的V/f切换到闭环矢量控制。这里的切换不是简单的“换一个控制器”而是要对电流状态做连续性处理不然切换瞬间电流会跳变。4.1 为什么要选id0控制电励磁同步电机在dq坐标系下转矩表达式是$$T_e \frac{3}{2} p \left( \psi_{sd} i_{sq} - \psi_{sq} i_{sd} \right)$$如果让id0励磁磁链全部由转子励磁电流建立转矩变成$$T_e \frac{3}{2} p M_{sf} i_f i_{sq}$$这就实现了磁场和转矩的完全解耦。iq控制转矩if控制磁场互不干扰。它和永磁同步电机的id0控制逻辑一模一样只是磁场的来源从永磁体换成了励磁绕组。这样做的好处是算法简单、物理意义清晰适合作为第一版控制器。缺点是需要励磁电流足够大对磁阻转矩完全没有利用不过对于追求稳定性的基础仿真这点损失可以接受。4.2 双闭环调节器参数整定外环是转速环内环是dq轴电流环。电流环带宽我选500-1000rad/s转速环带宽选50-100rad/s。参数整定用带宽法公式很简单$$Kp_{i} L \times \omega_{bw}, \quad Ki_{i} R_s \times \omega_{bw}$$$$Kp_{n} \frac{J \times \omega_{bn}}{1.5 p M_{sf} i_f}, \quad Ki_{n} \frac{Kp_{n}}{T_n}$$其中$T_n$是转速环积分时间常数我取0.1秒。按这个公式可以得到一组初值然后微调。实际整定时你会发现电流环的d轴和q轴电感不一样所以d轴和q轴的Kp要分开算别图省事用一个值。参数计算值实际采用值d轴电流环Kp10 mH × 800 rad/s 87.5d轴电流环Ki0.2 Ω × 800 160150q轴电流环Kp8 mH × 800 6.46q轴电流环Ki0.2 Ω × 800 160150转速环Kp0.4 × 60 / (1.5×2×0.3×3) ≈ 8.99转速环Ki9 / 0.1 9085上表是初值放到模型里还要根据仿真波形调。判断标准是电流阶跃响应无超调或者超调小于5%转速阶跃响应超调小于10%稳态转速误差小于0.5%。4.3 运行阶段带载与突加负载表现运行阶段我让电机在1.2秒时带上了额定负载9.5N·m考虑机械损耗后约等于额定输出。稳态时转速稳定在1500rpmid稳定在0附近iq稳定在对应负载电流约9A。1.5秒时我又做了一次突加50%负载的测试转速瞬间跌落了大约35rpm电流环和转速环在150毫秒内把转速拉回稳态。这说明双闭环参数基本是合理的。这个阶段如果想要更真实可以把负载模型从常数改成风机类负载或者恒功率负载。我建议至少留一个负载接口后面做负载突变扰动实验时非常有用。另外运行阶段的励磁电流控制也很关键我没有做弱磁但如果你的工况要求高速弱磁就需要增加一个励磁电流调节环当转速超过额定值后按1/转速的比例减小if保证反电动势不过压。这是EESM相对永磁电机最大的优势你可以在这个模型基础上直接扩展。5. 能耗制动阶段制动电阻选型与转矩衰减曲线能耗制动是这套模型里我个人觉得最有工程价值的部分。普通永磁同步电机做能耗制动只能靠短接定子绕组而电励磁同步电机可以在切断定子电源后继续控制转子励磁从而维持制动过程中的磁场强度制动的可控制性明显更好。5.1 能耗制动的物理过程从运行状态收到制动指令之后逆变器三相驱动全部封锁接触器把定子绕组切换到一组三相对称电阻上。这时电机转子仍有励磁电流由于转子还在转动定子绕组中会感应出三相电动势感应电流流过制动电阻产生制动转矩。随着转速下降感应电动势和电流同步衰减制动转矩也持续减小直到转速接近零。制动的本质是转子储存的动能 定子电阻和制动电阻上的焦耳热 励磁回路输入的能量 摩擦损耗。其中绝大部分动能最终消耗在制动电阻上。5.2 制动电阻的近似计算方法制动电阻阻值直接影响制动时间和制动电流峰值。阻值越小制动功率越大制动时间越短但定子电流峰值也越高可能超过电机允许的最大电流。我用的近似计算逻辑是设定一个期望的最大制动功率P_brake比如取额定功率的2倍即20kW。制动开始瞬间定子感应线电压约等于额定380V三相总制动功率近似为$$P_{brake} \approx \frac{3 U_{ph}^2}{R_b}$$其中$U_{ph}$是相电压有效值约220V。反解$$R_b \frac{3 \times 220^2}{20000} \approx 7.3 \Omega$$实际取标称值8Ω。初始制动电流峰值约等于相电压幅值除以单相回路电阻$$I_{peak} \approx \frac{220 \times \sqrt{2}}{8} \approx 39 A$$这个值超过了额定电流但因为制动过程持续时间短且制动电阻承担了大部分功率电机不会过热。如果想限制在额定电流附近就把电阻加大到15Ω制动时间会拉长到1秒以上。这是一个典型的“时间vs电流”折中。制动时间的粗略估算忽略摩擦和励磁损耗动能$E_k 0.5 J \omega^2$平均制动功率约10kW则制动时间约为$$t_{brake} \approx \frac{0.5 \times 0.4 \times (157 rad/s)^2}{10000} \approx 0.49 s$$仿真结果基本符合实测约0.55秒到0.6秒多出来的时间是因为制动转矩随转速下降不是恒定的。5.3 仿真结果衰减时间与制动深度我模型里的制动阶段波形是这样的制动指令发出后约20ms接触器完成切换定子电流在初始时刻到达约34A的峰值然后随转速线性下降转速从1500rpm在0.55秒内跌到接近0制动电流的包络基本是一条近似线性的衰减曲线和理论一致。整个过程励磁电流保持不变这是电励磁同步电机能耗制动的标志性特征。这里特别提醒一点制动电阻的功率等级选择。仿真里电阻不发热但实际工程里制动电阻一般选额定功率为制动平均功率的10%-20%因为它是短时工作制。如果你直接把仿真用的电阻拖到实物里长时间反复制动电阻会过热烧毁。6. 几个在实际仿真中真正折磨人的问题这部分是我个人把这套模型从零搭起来时耗费时间最多的地方。网上很多论文只展示最后的漂亮波形不会告诉你中间踩过多少坑。我挑四个最典型的说。6.1 励磁环节一上来就发散第一次跑模型我直接把励磁电压按稳态值36V一步给进去结果定子电流前10毫秒内出现几万安培的尖峰Simulink直接报数值错误。原因在上面说过励磁电流突变转子磁链瞬态变化经互感耦合到定子侧产生巨大的瞬态电动势而电流环还没介入定子侧就是纯电感回路电流必然冲高。解决方法是励磁斜坡给定从0开始用100ms斜坡上升到3A或者通过调节励磁电压环的输出限幅来实现。这在实际工程里也有对应——投励过程从来不是直接突加额定励磁电压而是要控制励磁电流上升速率。6.2 阶段切换瞬间定子电流冲击从启动的V/f控制切到运行的矢量控制时如果两个控制器的输出电压矢量和相位不一致切换的下一秒定子上就会叠加出一个很大的电压差产生电流冲击。我在驱动器设计里经常会遇到“无缝切换”这个词仿真也一样需要。我的处理方法是在切换指令发出后让控制器先用前一个阶段的电压矢量作为初始值然后通过一个一阶低通滤波器对电压指令做50ms的过渡让逆变器输出电压缓慢过渡到矢量控制的输出。这样电流没有突变切换波形几乎看不到痕迹。6.3 制动阶段求解器卡死或者波形振荡能耗制动阶段定子回路由大电感和小电阻串联时间常数比较小同时励磁回路时间常数又比较大系统里同时存在快变和慢变两部分数值刚性很强。用默认的ode45跑步长会变得极小仿真速度慢到无法接受有时候还会发散。我用的是ode23t变步长求解器相对容差设为1e-4最大步长限制为1毫秒。这样既保证了制动瞬间的精度又能把2.5秒的完整三阶段仿真实时性跑得很快。如果你的模型里加了开关器件和接触器动作Simulink里尽量用Simscape Electrical的阻抗模型而不是理想开关模型数值冲击会明显减少。6.4 无阻尼绕组导致制动结束阶段转子“爬行”能耗制动到转速很低时制动转矩随着反电动势按比例衰减最后可能还剩下几个rpm的转速长时间停不下来。我的模型里如果不做任何处理转速会在8rpm左右停留很久仿真的最后阶段看起来很拖沓。我加了一个判断逻辑当转速低于10rpm时封锁制动接触器并投入一个很小的直流制动电流或者直接结束仿真流程。在实际工业设备里这个低速爬行也是能耗制动的常见短板通常要靠机械抱闸或直流制动补一刀。你在仿真里提前把这个逻辑做进去后面做工程验证时会有参考。做完这套三阶段仿真最大的感触是电励磁同步电机的每一个阶段单独拿出来都不算难难的是三个阶段的平滑衔接和参数协调。启动阶段的V/f参数会影响切入矢量控制时的电流冲击制动电阻阻值不只影响制动时间还会反过来影响接触器切换瞬间的过电压。我那组参数最后能跑出平稳的完整曲线是在启动斜坡时间、励磁斜坡时间、制动电阻几十个组合里试出来的。建议你也别迷信某一组“最优参数”而是多用几个工况去压这个模型把每段过渡的物理过程看过一遍这套仿真才算真正属于你了。
