简介基于Matlab的三相异步电动机起动、调速和制动特性仿真设计报告面向电气工程、自动化及电机控制方向的学生和工程技术人员用于理解异步电动机动态特性并掌握Simulink仿真建模方法。资源为1份完整docx文档整体约1.47MB内容涵盖电动机基础理论、直接起动/星三角起动/变频器起动、变频调速/转子电阻调速/能量回馈调速以及反接制动/能耗制动/再生制动等典型工况的原理说明、仿真模型搭建、参数设置与结果分析并附有任务书、摘要和设计流程便于对照学习或直接参考。已有525人学习适合作为课程设计、毕业设计或自学三相异步电动机仿真与运行特性的系统资料。1. 三相异步电动机特性仿真在 Matlab 里到底仿什么一台三相异步电动机直接起动时电流能冲到额定值的 5 到 7 倍转速却要等一两秒才爬到位正常运行后想调速要么改变极对数、调定子电压要么改变供电频率停机时如果直接断电靠摩擦慢慢停很多产线根本等不起。这三个问题分别对应起动、调速、制动特性是电机学里最容易讲清楚、也最容易在参数上做错的地方。用 Matlab 把这三种工况放进同一套电机参数里仿真可以一次性对比起动冲击电流、调速范围和制动时间不用再靠手算等效电路去猜波形。本文适合正在做电机课程设计、准备设备选型或者刚开始接触变频器和软起动器的工程师照着下面的 dq 模型和 Simulink 步骤可以自行复现。2. 写状态方程还是拖模块三相异步电动机仿真模型的两种搭法2.1 为什么用 dq 模型而不是直接解 abc 三相模型直接对 abc 三相写出定转子电压方程电感矩阵里会出现转子位置角的正余弦项。转子一转这些项全部随时间变化数值积分每一步都要重新组装带时变电感的矩阵而且在三相电压不平衡或者暂态电流很大的时候步长被逼得很小。dq 变换把三相绕组等效到随磁场同步旋转的两相坐标上时变耦合项被消掉电感矩阵变成常数矩阵方程退化成一组常系数微分方程ode45 就能稳定求解。常见的做法是采用等幅值 Park 变换。变换前后的功率数值不一致但电流幅值和实际相电流峰值直接对应画起动电流冲击时读起来很直观。相应地电磁转矩表达式里会出现 1.5 这个系数对应等幅值坐标变换的功率修正。稳态特性计算仍然推荐用 T 型等效电路磁路饱和、集肤效应在等效电路里更容易手动补偿动态过程、起动暂态、制动过程则必须用 dq 状态方程。2.2 三相异步电动机的状态方程怎么列电压、磁链与转矩鼠笼型转子的电压方程右侧为零只有定子侧有输入电压。额定运行点附近假设磁路线性、忽略铁耗那么微分方程系统由五阶状态组成定子 dq 轴磁链、转子 dq 轴磁链外加机械转速。表 1 给出了 abc 模型与 dq 模型的差别方便在一开始就决定自己做仿真用哪一种。对比项abc 三相模型dq 旋转坐标模型电感矩阵含转子位置角时变常数矩阵状态个数7 阶6 个绕组电流加转速5 阶4 个磁链加转速数值积分步长敏感容易发散ode45 默认容差即可稳定适用场景三相不平衡、定子匝间短路等起动、调速、制动等对称工况dq 坐标下的定子电压方程、磁链方程与运动方程按下式联立。定子侧输入为 u_ds 和 u_qs转子侧短路即 u_dr u_qr 0ψds Ls·ids Lm·idr ψqs Ls·iqs Lm·iqr ψdr Lr·idr Lm·ids ψqr Lr·iqr Lm·iqs电磁转矩在等幅值 Park 变换下写为Te 1.5·p·(ψds·iqs - ψqs·ids)运动方程为J·dωm/dt Te - Tl。注意这里的 ωm 是机械角速度单位 rad/s而磁链微分方程里出现的电角速度是 ωe转差电角速度是 ωe 减去 p 倍的机械角速度极对数 p 把两者桥接起来。初学最容易漏的就是 1.5 这个系数或者极对数漏掉之后转矩数值整体偏小起动时间会被拉长。2.3 用 M 脚本把三相异步电动机模型写成可调函数下面这段代码是全文仿真的核心所有起动、调速、制动仿真都复用它。函数输入是时间 t 和状态向量 x输出是五个状态变量的导数符合 ode45 对目标函数的要求。function dx motor_dq(t, x) % 三相异步电动机 dq 模型鼠笼转子短路 % x [psid_s; psiq_s; psid_r; psiq_r; omega_m] % ---- 电机参数3kW 四级电机数量级---- Rs 0.435; % 定子电阻欧姆 Rr 0.816; % 转子电阻折算值欧姆 Lls 0.01; Llr 0.01; % 定转子漏感亨 Lm 0.4; % 励磁电感亨 p 2; % 极对数 J 0.089; % 转动惯量kg·m^2 Tl 0; % 负载转矩空载起动设 0 Ls Lls Lm; Lr Llr Lm; den Ls * Lr - Lm^2; % ---- 由磁链反算电流 ---- i_ds (Lr * x(1) - Lm * x(3)) / den; i_qs (Lr * x(2) - Lm * x(4)) / den; i_dr (Ls * x(3) - Lm * x(1)) / den; i_qr (Ls * x(4) - Lm * x(2)) / den; % ---- 定子电压直接起动三相平衡正弦电压的 dq 表示 ---- f 50; % 供电频率 Hz we 2*pi*f; % 同步电角速度 U_phase 220; % 定子相电压有效值 u_ds U_phase * sqrt(2); % d 轴恒定电压 u_qs 0; % ---- 转子磁链微分鼠笼转子短路 ---- dpsi_dr -Rr * i_dr (we - p*x(5)) * x(4); dpsi_qr -Rr * i_qr - (we - p*x(5)) * x(3); % ---- 定子磁链微分 ---- dpsi_ds u_ds - Rs * i_ds we * x(2); dpsi_qs u_qs - Rs * i_qs - we * x(1); % ---- 转矩与运动方程 ---- Te 1.5 * p * (x(1) * i_qs - x(2) * i_ds); domega (Te - Tl) / J; dx [dpsi_ds; dpsi_qs; dpsi_dr; dpsi_qr; domega]; end代码里采用恒定的 d 轴电压分量来等效三相正弦供电省去了在每一个积分步计算三相瞬时值的开销。参数里 Rs、Rr 和漏感直接决定了起动电流和最大转矩漏感偏大起动电流偏小但起动转矩也会掉励磁电感 Lm 主要影响空载电流。改参数时优先保持 Rs、Rr 数量级不变只调 Lm 和漏感来匹配你手里的电机铭牌。初始状态全为零对应转速零、磁场零的断电状态直接仿真就是完整的直接起动过程。需要在 Simulink 里对照验证时可以在 Matlab 命令窗口执行powerlib打开 Simscape Electrical 库里面 Asynchronous Machine SI Units 模块的参数页与上面 M 脚本的参数一一对应。安装 Matlab 时勾选 Simulink 和 Simscape Electrical 两个组件机器上没装的话仿真会报“块未定义”之类的错误。3. 起动特性仿真直接起动、Y-Δ 起动与软起动的暂态对比3.1 起动电流为什么是额定电流的 5 到 7 倍起动瞬间转速为零转差率 s 1转子等效电阻 Rr/s 就等于 Rr 本身数值很小这时的定子电流只受定转子电阻和漏抗限制。用上一节的参数估算堵转阻抗约 6 欧姆量级起动电流峰值接近 50 安培而额定相电流峰值只有 9 安培左右。这就是阶跃电压激励下的典型响应电流瞬间冲到峰值随着转速上升、转差率减小转子等效电阻 Rr/s 逐渐变大电流回落整个过程是一条“尖峰后衰减”的曲线。3.2 直接起动仿真看电流尖峰和转速爬升将motor_dq函数保存后在另一个脚本里调用 ode45并计算定子电流幅值和机械转速% 直接起动仿真 t_end 1.5; % 仿真时长 1.5 秒 [t, y] ode45(motor_dq, [0 t_end], zeros(5,1)); % 反算定子电流 dq 分量等幅值变换下幅值即相电流峰值 Rs 0.435; Rr 0.816; Lls 0.01; Llr 0.01; Lm 0.4; Ls Lls Lm; Lr Llr Lm; den Ls * Lr - Lm^2; i_ds (Lr * y(:,1) - Lm * y(:,3)) / den; i_qs (Lr * y(:,2) - Lm * y(:,4)) / den; i_peak sqrt(i_ds.^2 i_qs.^2); % 相电流峰值 n y(:,5) * 60 / (2*pi); % 机械转速 r/min figure; yyaxis left; plot(t, i_peak, LineWidth, 1.2); ylabel(定子电流峰值 (A)); yyaxis right; plot(t, n, LineWidth, 1.2); ylabel(转速 (r/min)); xlabel(时间 (s));从图上可以读到两个关键特征。第一个是起动瞬间电流峰值大约是额定峰值的 5 到 7 倍这个倍数由漏抗决定与负载轻重基本无关第二个是电流从峰值开始回落的时间点和转速接近同步转速的时间点吻合说明电流下降不是电压变化引起的而是反电动势随转速升高建立起来的结果。如果观察到电流长时间不回落优先检查负荷转矩 Tl 是否设得太大或者转子参数 Rr 数量级不对。3.3 Y-Δ 起动怎么仿真降压到额定电压的 1/√3Y-Δ 起动只适用于正常运行时三角形接法的电机。起动时把定子绕组改成星形每相绕组两端电压降为三角形接法时的 1/√3转矩与电压平方成正比起动转矩也随之降为原来的三分之一。仿真做法很简单把motor_dq函数里的u_ds从220*sqrt(2)改成220*sqrt(2)/sqrt(3)其余参数不变重跑一遍直接起动脚本即可。表 2 是三种起动方式在同一台电机上的对比结果数值来自上面这组参数在空载起动下的稳态读取。起动方式起动电流倍数起动转矩倍数适用场景直接起动57 倍额定电流12 倍额定转矩小容量、电网容量充裕Y-Δ 起动1/3相对直接起动1/3正常 Δ 接、轻载起动软起动电压斜坡24 倍可调按电压平方下降大惯量、对机械冲击敏感Y-Δ 起动节省的是电网冲击代价是起动转矩同步缩水。如果负载要求的起动转矩超过额定转矩的一半Y-Δ 起动可能起不来。判断方法是在仿真图里叠加负载转矩曲线看转速曲线是否在某处停滞停滞就说明电磁转矩与负载转矩的交点低于同步转速需要换更高起动转矩的方案。3.4 软起动用电压斜坡限制电流冲击软起动不是斩波调压而是在起动时间窗内让定子电压从零线性上升。对应到 M 脚本只需要把电压赋值处改成随 t 变化的斜坡% 软起动0.8 秒内电压从 0 线性升到额定 t_ramp 0.8; u_ds 220 * sqrt(2) * min(t / t_ramp, 1);这个写法避免了在状态方程里引入离散逻辑ode45 可以平滑跨过斜坡段。电压上升速率越慢电流峰值越小但起动时间越长。需要注意软起动不能提高起动转矩因为电压降低的同时转矩按平方关系下降风机、水泵这类轻载起动场合最合适重载场合应该考虑转子串电阻或者直接上变频起动。4. 调速特性仿真调压调速与变频调速的稳态和动态边界4.1 调速的可能性空间从转速公式看三个调节变量三相异步电动机的同步转速由n1 60·f1 / p决定实际转速是n n1·(1-s)。所以调速有三个切入路径改变极对数 p改变供电频率 f1或者改变转差率 s。变极调速是有级的一般只有 2 到 3 挡改变转差率最常用的是调定子电压变频调速则同时影响同步转速和转差率是当前最主流的方案。仿真时只需要修改motor_dq函数里的频率和电压赋值同一个状态方程可以覆盖这三种思路。4.2 调压调速转差率的变化为什么非常有限调压调速的原理是基于转矩与电压平方成正比电压降低机械特性曲线整体下移与负载转矩的交点向转差率增大的方向移动。先看稳态特性用 T 型等效电路画不同电压下的机械特性簇% 调压调速机械特性计算 f1 50; p 2; ns 60*f1/p; % 同步转速 R1 0.435; R2 0.816; X1 1.0; X2 1.0; omega_s 2*pi*f1/p; % 同步机械角速度 s linspace(0.001, 0.5, 1000); U_list [0.5 0.75 1.0] * 220; % 不同电压倍数 figure; hold on; for U U_list T 3 * U^2 * R2 ./ s ./ ... ( (R1 R2./s).^2 (X1X2).^2 ) / omega_s; plot((1-s)*ns, T, LineWidth, 1.2); end % 加一条负载转矩曲线泵类负载 T ∝ n^2 n (1-s)*ns; plot(n, 20*(n/ns).^2, k--, LineWidth, 1.2); xlabel(转速 (r/min)); ylabel(转矩 (N·m)); legend(0.5U,0.75U,1.0U,负载);从图里能看到一个关键边界在转差率很小的工作段特性曲线很硬电压从额定降到一半转速几乎没变只有进入转差率较大的非线性段转速变化才明显但这时转子铜耗增大电机发热严重。因此调压调速只适合风机泵类平方转矩负载以及对调速精度要求不高的场合。仿真时要注意把负载从恒转矩改成Tl k·n²的形式否则会得到“电压调不下去”的错误结论。4.3 变频调速V/f 恒定曲线与机械特性簇变频调速在基频以下要保持磁通近似恒定通常按U1/f1 常数给定电压。频率降低后漏抗也随之降低用等效电路画曲线时定子漏抗必须同步修改否则低频段会算出比实际偏大的转矩。% 变频调速机械特性簇基频以下恒转矩 p 2; R1 0.435; R2 0.816; Lls 0.01; Llr 0.01; % 漏感亨 s linspace(0.001, 0.3, 1000); f_list [10 25 40 50]; % 不同频率 figure; hold on; for f1 f_list U1 220 * f1 / 50; % V/f 恒定 X1 2*pi*f1*Lls; X2 2*pi*f1*Llr; % 漏抗随频率变化 omega_s 2*pi*f1/p; ns_f 60*f1/p; T 3 * U1^2 * R2 ./ s ./ ... ( (R1 R2./s).^2 (X1X2).^2 ) / omega_s; plot((1-s)*ns_f, T, LineWidth, 1.2); end xlabel(转速 (r/min)); ylabel(转矩 (N·m)); legend(10Hz,25Hz,40Hz,50Hz);低频时定子电阻压降占比明显增大实际工程会在 V/f 曲线上叠加一个起始电压补偿也就是常说的 boost。仿真中可以在U1的计算式里加一个固定偏置电压比如U1 220*f1/50 15然后对比加与不加的起动转矩差异。做闭环调速时转速外环可以用 Matlab 的 PID Tuner 自动整定在 Simulink 里搭一个带转速反馈的 V/f 控制框双击 PID Controller 模块打开 Tune 按钮根据阶跃响应目标直接生成 Kp、Ki 参数比自己试凑快得多。4.4 调速结果的三个验证指标仿真完成后不要只看转速曲线好不好看要量化验证调速性能。静差率按下式计算静差率 (n0 - nN) / n0 · 100%其中 n0 是同一频率下的空载转速nN 是额定负载转速静差率越小特性越硬。调速范围D nmax / nmin前提是最低转速 nmin 仍满足静差率要求而不是零转速也算进去。第三个指标是转矩能力基频以下恒转矩调速基频以上电压达到额定后只能弱磁升速转矩随转速上升而下降。把这三项数据从仿真曲线里读出来与电机铭牌和负载需求对照才能判断一个调速方案是否真的可用。5. 制动特性仿真三种制动的波形检查方法和参数收口5.1 能耗制动、反接制动与回馈制动的本质区别从转差率看三种制动最清楚能耗制动是将定子从电网断开并通入直流产生静止磁场转子导体切割静止磁场产生制动转矩转速为零时制动转矩自然消失不会反转反接制动是任意对调两相电源线让旋转磁场反向旋转转差率大于 1制动转矩大但转速过零时必须切断电源否则电机反向起动回馈制动发生在实际转速高于同步转速时转差率为负电机处于发电状态。三者的适用场景完全不同能耗制动用于快速平稳停车反接制动用于要求制动转矩最大的场合回馈制动用于变频器减速过程或位能性负载下放。5.2 能耗制动的直流电流参数选择能耗制动时接通定子的直流电流一般取电机空载电流的 2 到 4 倍太小制动力不足太大会烧绕组。仿真时可以在三相电流源里把频率设为 0或者直接在 M 脚本里把 u_qs 设为 0、u_ds 改为直流电压值同时去掉交流频率项。制动时间可以从波形上读取从制动投入瞬间到转速降到额定转速的 5% 以下所用的时间就是工程上常见的“停车时间”。5.3 用一个小脚本从仿真数据里读出制动时间Simulink 中做反接制动只要把三相电压源的 B、C 相交换即可转速过零切断用 Relay 模块检测转速符号变化。仿真完成后把 Scope 数据保存到工作区用下面的代码计算制动时间out sim(brake_sim.slx); t out.tout; n out.yout{1}.Values.Data; % 转速信号 t_brake_start 2.0; % 制动投入时刻 idx find(t t_brake_start n 0.05*max(n), 1); tb t(idx) - t_brake_start; fprintf(制动时间: %.3f s\n, tb);改变直流电流的大小或者反接电压的幅值再对比 tb 的变化就能看出制动转矩与制动时间之间的取舍。实际调试时常用set_param一次性修改多个参数并循环仿真把不同参数下的转速波形叠加在同一张图上比一次一次手动改参数再截图要高效得多。本文还有配套的精品资源点击获取
简介基于Matlab的三相异步电动机起动、调速和制动特性仿真设计报告面向电气工程、自动化及电机控制方向的学生和工程技术人员用于理解异步电动机动态特性并掌握Simulink仿真建模方法。资源为1份完整docx文档整体约1.47MB内容涵盖电动机基础理论、直接起动/星三角起动/变频器起动、变频调速/转子电阻调速/能量回馈调速以及反接制动/能耗制动/再生制动等典型工况的原理说明、仿真模型搭建、参数设置与结果分析并附有任务书、摘要和设计流程便于对照学习或直接参考。已有525人学习适合作为课程设计、毕业设计或自学三相异步电动机仿真与运行特性的系统资料。1. 三相异步电动机特性仿真在 Matlab 里到底仿什么一台三相异步电动机直接起动时电流能冲到额定值的 5 到 7 倍转速却要等一两秒才爬到位正常运行后想调速要么改变极对数、调定子电压要么改变供电频率停机时如果直接断电靠摩擦慢慢停很多产线根本等不起。这三个问题分别对应起动、调速、制动特性是电机学里最容易讲清楚、也最容易在参数上做错的地方。用 Matlab 把这三种工况放进同一套电机参数里仿真可以一次性对比起动冲击电流、调速范围和制动时间不用再靠手算等效电路去猜波形。本文适合正在做电机课程设计、准备设备选型或者刚开始接触变频器和软起动器的工程师照着下面的 dq 模型和 Simulink 步骤可以自行复现。2. 写状态方程还是拖模块三相异步电动机仿真模型的两种搭法2.1 为什么用 dq 模型而不是直接解 abc 三相模型直接对 abc 三相写出定转子电压方程电感矩阵里会出现转子位置角的正余弦项。转子一转这些项全部随时间变化数值积分每一步都要重新组装带时变电感的矩阵而且在三相电压不平衡或者暂态电流很大的时候步长被逼得很小。dq 变换把三相绕组等效到随磁场同步旋转的两相坐标上时变耦合项被消掉电感矩阵变成常数矩阵方程退化成一组常系数微分方程ode45 就能稳定求解。常见的做法是采用等幅值 Park 变换。变换前后的功率数值不一致但电流幅值和实际相电流峰值直接对应画起动电流冲击时读起来很直观。相应地电磁转矩表达式里会出现 1.5 这个系数对应等幅值坐标变换的功率修正。稳态特性计算仍然推荐用 T 型等效电路磁路饱和、集肤效应在等效电路里更容易手动补偿动态过程、起动暂态、制动过程则必须用 dq 状态方程。2.2 三相异步电动机的状态方程怎么列电压、磁链与转矩鼠笼型转子的电压方程右侧为零只有定子侧有输入电压。额定运行点附近假设磁路线性、忽略铁耗那么微分方程系统由五阶状态组成定子 dq 轴磁链、转子 dq 轴磁链外加机械转速。表 1 给出了 abc 模型与 dq 模型的差别方便在一开始就决定自己做仿真用哪一种。对比项abc 三相模型dq 旋转坐标模型电感矩阵含转子位置角时变常数矩阵状态个数7 阶6 个绕组电流加转速5 阶4 个磁链加转速数值积分步长敏感容易发散ode45 默认容差即可稳定适用场景三相不平衡、定子匝间短路等起动、调速、制动等对称工况dq 坐标下的定子电压方程、磁链方程与运动方程按下式联立。定子侧输入为 u_ds 和 u_qs转子侧短路即 u_dr u_qr 0ψds Ls·ids Lm·idr ψqs Ls·iqs Lm·iqr ψdr Lr·idr Lm·ids ψqr Lr·iqr Lm·iqs电磁转矩在等幅值 Park 变换下写为Te 1.5·p·(ψds·iqs - ψqs·ids)运动方程为J·dωm/dt Te - Tl。注意这里的 ωm 是机械角速度单位 rad/s而磁链微分方程里出现的电角速度是 ωe转差电角速度是 ωe 减去 p 倍的机械角速度极对数 p 把两者桥接起来。初学最容易漏的就是 1.5 这个系数或者极对数漏掉之后转矩数值整体偏小起动时间会被拉长。2.3 用 M 脚本把三相异步电动机模型写成可调函数下面这段代码是全文仿真的核心所有起动、调速、制动仿真都复用它。函数输入是时间 t 和状态向量 x输出是五个状态变量的导数符合 ode45 对目标函数的要求。function dx motor_dq(t, x) % 三相异步电动机 dq 模型鼠笼转子短路 % x [psid_s; psiq_s; psid_r; psiq_r; omega_m] % ---- 电机参数3kW 四级电机数量级---- Rs 0.435; % 定子电阻欧姆 Rr 0.816; % 转子电阻折算值欧姆 Lls 0.01; Llr 0.01; % 定转子漏感亨 Lm 0.4; % 励磁电感亨 p 2; % 极对数 J 0.089; % 转动惯量kg·m^2 Tl 0; % 负载转矩空载起动设 0 Ls Lls Lm; Lr Llr Lm; den Ls * Lr - Lm^2; % ---- 由磁链反算电流 ---- i_ds (Lr * x(1) - Lm * x(3)) / den; i_qs (Lr * x(2) - Lm * x(4)) / den; i_dr (Ls * x(3) - Lm * x(1)) / den; i_qr (Ls * x(4) - Lm * x(2)) / den; % ---- 定子电压直接起动三相平衡正弦电压的 dq 表示 ---- f 50; % 供电频率 Hz we 2*pi*f; % 同步电角速度 U_phase 220; % 定子相电压有效值 u_ds U_phase * sqrt(2); % d 轴恒定电压 u_qs 0; % ---- 转子磁链微分鼠笼转子短路 ---- dpsi_dr -Rr * i_dr (we - p*x(5)) * x(4); dpsi_qr -Rr * i_qr - (we - p*x(5)) * x(3); % ---- 定子磁链微分 ---- dpsi_ds u_ds - Rs * i_ds we * x(2); dpsi_qs u_qs - Rs * i_qs - we * x(1); % ---- 转矩与运动方程 ---- Te 1.5 * p * (x(1) * i_qs - x(2) * i_ds); domega (Te - Tl) / J; dx [dpsi_ds; dpsi_qs; dpsi_dr; dpsi_qr; domega]; end代码里采用恒定的 d 轴电压分量来等效三相正弦供电省去了在每一个积分步计算三相瞬时值的开销。参数里 Rs、Rr 和漏感直接决定了起动电流和最大转矩漏感偏大起动电流偏小但起动转矩也会掉励磁电感 Lm 主要影响空载电流。改参数时优先保持 Rs、Rr 数量级不变只调 Lm 和漏感来匹配你手里的电机铭牌。初始状态全为零对应转速零、磁场零的断电状态直接仿真就是完整的直接起动过程。需要在 Simulink 里对照验证时可以在 Matlab 命令窗口执行powerlib打开 Simscape Electrical 库里面 Asynchronous Machine SI Units 模块的参数页与上面 M 脚本的参数一一对应。安装 Matlab 时勾选 Simulink 和 Simscape Electrical 两个组件机器上没装的话仿真会报“块未定义”之类的错误。3. 起动特性仿真直接起动、Y-Δ 起动与软起动的暂态对比3.1 起动电流为什么是额定电流的 5 到 7 倍起动瞬间转速为零转差率 s 1转子等效电阻 Rr/s 就等于 Rr 本身数值很小这时的定子电流只受定转子电阻和漏抗限制。用上一节的参数估算堵转阻抗约 6 欧姆量级起动电流峰值接近 50 安培而额定相电流峰值只有 9 安培左右。这就是阶跃电压激励下的典型响应电流瞬间冲到峰值随着转速上升、转差率减小转子等效电阻 Rr/s 逐渐变大电流回落整个过程是一条“尖峰后衰减”的曲线。3.2 直接起动仿真看电流尖峰和转速爬升将motor_dq函数保存后在另一个脚本里调用 ode45并计算定子电流幅值和机械转速% 直接起动仿真 t_end 1.5; % 仿真时长 1.5 秒 [t, y] ode45(motor_dq, [0 t_end], zeros(5,1)); % 反算定子电流 dq 分量等幅值变换下幅值即相电流峰值 Rs 0.435; Rr 0.816; Lls 0.01; Llr 0.01; Lm 0.4; Ls Lls Lm; Lr Llr Lm; den Ls * Lr - Lm^2; i_ds (Lr * y(:,1) - Lm * y(:,3)) / den; i_qs (Lr * y(:,2) - Lm * y(:,4)) / den; i_peak sqrt(i_ds.^2 i_qs.^2); % 相电流峰值 n y(:,5) * 60 / (2*pi); % 机械转速 r/min figure; yyaxis left; plot(t, i_peak, LineWidth, 1.2); ylabel(定子电流峰值 (A)); yyaxis right; plot(t, n, LineWidth, 1.2); ylabel(转速 (r/min)); xlabel(时间 (s));从图上可以读到两个关键特征。第一个是起动瞬间电流峰值大约是额定峰值的 5 到 7 倍这个倍数由漏抗决定与负载轻重基本无关第二个是电流从峰值开始回落的时间点和转速接近同步转速的时间点吻合说明电流下降不是电压变化引起的而是反电动势随转速升高建立起来的结果。如果观察到电流长时间不回落优先检查负荷转矩 Tl 是否设得太大或者转子参数 Rr 数量级不对。3.3 Y-Δ 起动怎么仿真降压到额定电压的 1/√3Y-Δ 起动只适用于正常运行时三角形接法的电机。起动时把定子绕组改成星形每相绕组两端电压降为三角形接法时的 1/√3转矩与电压平方成正比起动转矩也随之降为原来的三分之一。仿真做法很简单把motor_dq函数里的u_ds从220*sqrt(2)改成220*sqrt(2)/sqrt(3)其余参数不变重跑一遍直接起动脚本即可。表 2 是三种起动方式在同一台电机上的对比结果数值来自上面这组参数在空载起动下的稳态读取。起动方式起动电流倍数起动转矩倍数适用场景直接起动57 倍额定电流12 倍额定转矩小容量、电网容量充裕Y-Δ 起动1/3相对直接起动1/3正常 Δ 接、轻载起动软起动电压斜坡24 倍可调按电压平方下降大惯量、对机械冲击敏感Y-Δ 起动节省的是电网冲击代价是起动转矩同步缩水。如果负载要求的起动转矩超过额定转矩的一半Y-Δ 起动可能起不来。判断方法是在仿真图里叠加负载转矩曲线看转速曲线是否在某处停滞停滞就说明电磁转矩与负载转矩的交点低于同步转速需要换更高起动转矩的方案。3.4 软起动用电压斜坡限制电流冲击软起动不是斩波调压而是在起动时间窗内让定子电压从零线性上升。对应到 M 脚本只需要把电压赋值处改成随 t 变化的斜坡% 软起动0.8 秒内电压从 0 线性升到额定 t_ramp 0.8; u_ds 220 * sqrt(2) * min(t / t_ramp, 1);这个写法避免了在状态方程里引入离散逻辑ode45 可以平滑跨过斜坡段。电压上升速率越慢电流峰值越小但起动时间越长。需要注意软起动不能提高起动转矩因为电压降低的同时转矩按平方关系下降风机、水泵这类轻载起动场合最合适重载场合应该考虑转子串电阻或者直接上变频起动。4. 调速特性仿真调压调速与变频调速的稳态和动态边界4.1 调速的可能性空间从转速公式看三个调节变量三相异步电动机的同步转速由n1 60·f1 / p决定实际转速是n n1·(1-s)。所以调速有三个切入路径改变极对数 p改变供电频率 f1或者改变转差率 s。变极调速是有级的一般只有 2 到 3 挡改变转差率最常用的是调定子电压变频调速则同时影响同步转速和转差率是当前最主流的方案。仿真时只需要修改motor_dq函数里的频率和电压赋值同一个状态方程可以覆盖这三种思路。4.2 调压调速转差率的变化为什么非常有限调压调速的原理是基于转矩与电压平方成正比电压降低机械特性曲线整体下移与负载转矩的交点向转差率增大的方向移动。先看稳态特性用 T 型等效电路画不同电压下的机械特性簇% 调压调速机械特性计算 f1 50; p 2; ns 60*f1/p; % 同步转速 R1 0.435; R2 0.816; X1 1.0; X2 1.0; omega_s 2*pi*f1/p; % 同步机械角速度 s linspace(0.001, 0.5, 1000); U_list [0.5 0.75 1.0] * 220; % 不同电压倍数 figure; hold on; for U U_list T 3 * U^2 * R2 ./ s ./ ... ( (R1 R2./s).^2 (X1X2).^2 ) / omega_s; plot((1-s)*ns, T, LineWidth, 1.2); end % 加一条负载转矩曲线泵类负载 T ∝ n^2 n (1-s)*ns; plot(n, 20*(n/ns).^2, k--, LineWidth, 1.2); xlabel(转速 (r/min)); ylabel(转矩 (N·m)); legend(0.5U,0.75U,1.0U,负载);从图里能看到一个关键边界在转差率很小的工作段特性曲线很硬电压从额定降到一半转速几乎没变只有进入转差率较大的非线性段转速变化才明显但这时转子铜耗增大电机发热严重。因此调压调速只适合风机泵类平方转矩负载以及对调速精度要求不高的场合。仿真时要注意把负载从恒转矩改成Tl k·n²的形式否则会得到“电压调不下去”的错误结论。4.3 变频调速V/f 恒定曲线与机械特性簇变频调速在基频以下要保持磁通近似恒定通常按U1/f1 常数给定电压。频率降低后漏抗也随之降低用等效电路画曲线时定子漏抗必须同步修改否则低频段会算出比实际偏大的转矩。% 变频调速机械特性簇基频以下恒转矩 p 2; R1 0.435; R2 0.816; Lls 0.01; Llr 0.01; % 漏感亨 s linspace(0.001, 0.3, 1000); f_list [10 25 40 50]; % 不同频率 figure; hold on; for f1 f_list U1 220 * f1 / 50; % V/f 恒定 X1 2*pi*f1*Lls; X2 2*pi*f1*Llr; % 漏抗随频率变化 omega_s 2*pi*f1/p; ns_f 60*f1/p; T 3 * U1^2 * R2 ./ s ./ ... ( (R1 R2./s).^2 (X1X2).^2 ) / omega_s; plot((1-s)*ns_f, T, LineWidth, 1.2); end xlabel(转速 (r/min)); ylabel(转矩 (N·m)); legend(10Hz,25Hz,40Hz,50Hz);低频时定子电阻压降占比明显增大实际工程会在 V/f 曲线上叠加一个起始电压补偿也就是常说的 boost。仿真中可以在U1的计算式里加一个固定偏置电压比如U1 220*f1/50 15然后对比加与不加的起动转矩差异。做闭环调速时转速外环可以用 Matlab 的 PID Tuner 自动整定在 Simulink 里搭一个带转速反馈的 V/f 控制框双击 PID Controller 模块打开 Tune 按钮根据阶跃响应目标直接生成 Kp、Ki 参数比自己试凑快得多。4.4 调速结果的三个验证指标仿真完成后不要只看转速曲线好不好看要量化验证调速性能。静差率按下式计算静差率 (n0 - nN) / n0 · 100%其中 n0 是同一频率下的空载转速nN 是额定负载转速静差率越小特性越硬。调速范围D nmax / nmin前提是最低转速 nmin 仍满足静差率要求而不是零转速也算进去。第三个指标是转矩能力基频以下恒转矩调速基频以上电压达到额定后只能弱磁升速转矩随转速上升而下降。把这三项数据从仿真曲线里读出来与电机铭牌和负载需求对照才能判断一个调速方案是否真的可用。5. 制动特性仿真三种制动的波形检查方法和参数收口5.1 能耗制动、反接制动与回馈制动的本质区别从转差率看三种制动最清楚能耗制动是将定子从电网断开并通入直流产生静止磁场转子导体切割静止磁场产生制动转矩转速为零时制动转矩自然消失不会反转反接制动是任意对调两相电源线让旋转磁场反向旋转转差率大于 1制动转矩大但转速过零时必须切断电源否则电机反向起动回馈制动发生在实际转速高于同步转速时转差率为负电机处于发电状态。三者的适用场景完全不同能耗制动用于快速平稳停车反接制动用于要求制动转矩最大的场合回馈制动用于变频器减速过程或位能性负载下放。5.2 能耗制动的直流电流参数选择能耗制动时接通定子的直流电流一般取电机空载电流的 2 到 4 倍太小制动力不足太大会烧绕组。仿真时可以在三相电流源里把频率设为 0或者直接在 M 脚本里把 u_qs 设为 0、u_ds 改为直流电压值同时去掉交流频率项。制动时间可以从波形上读取从制动投入瞬间到转速降到额定转速的 5% 以下所用的时间就是工程上常见的“停车时间”。5.3 用一个小脚本从仿真数据里读出制动时间Simulink 中做反接制动只要把三相电压源的 B、C 相交换即可转速过零切断用 Relay 模块检测转速符号变化。仿真完成后把 Scope 数据保存到工作区用下面的代码计算制动时间out sim(brake_sim.slx); t out.tout; n out.yout{1}.Values.Data; % 转速信号 t_brake_start 2.0; % 制动投入时刻 idx find(t t_brake_start n 0.05*max(n), 1); tb t(idx) - t_brake_start; fprintf(制动时间: %.3f s\n, tb);改变直流电流的大小或者反接电压的幅值再对比 tb 的变化就能看出制动转矩与制动时间之间的取舍。实际调试时常用set_param一次性修改多个参数并循环仿真把不同参数下的转速波形叠加在同一张图上比一次一次手动改参数再截图要高效得多。本文还有配套的精品资源点击获取
