永磁同步直线电机建模仿真及优化教学 从入门到精通永磁同步直线电机高速精密绘图仪笔尖能在纸上跑出米每秒级速度却连发丝粗细的误差都没有晶圆台托着指甲盖大的芯片在光刻机里微米级挪位卡得死死的这些“直来直去还准到离谱”的家伙核心大多都是永磁同步直线电机——PMLSM今天咱们就从工科生最挠头的“模型怎么搭”“搭完怎么跑”“跑出来怎么改”三个点唠从只会搜论文摘要到能对着自己的仿真结果调参数调结构。第一站入门篇别信“模型都一样”先搭个“最朴素的能用就行”的PMLSM模型很多新手上来就啃有限元分析FEA的黑箱子跑三天两夜出来个波形还不对参数的锅何苦呢PMLSM本质就是把旋转的永磁同步电机PMSM“剪一刀展开成直线”所以先用MATLAB/Simulink里的PMSM模块“改一改直线的味道”搭个控制环全的集中参数模型先爽爽高速定位的感觉再说。第一步先摸清楚“剪一刀展开”到底改了啥参数磁链的变化率旋转的是“角速度ω2πn/60”展开成直线就是“线速度v2πnτ/60”τ是极距简单说就是相邻两个磁铁南极北极的直线距离电磁推力公式旋转的是“电磁转矩Te”展开成直线就是把旋转的半径一般是气隙半径rg换成线速度对应的系数——其实就是Te/rg但集中参数模型里更常用的推导是Fem (3π/pτ)ψfiq (3π/pτ)(Ld - Lq)idiqp是极对数ψf是永磁体磁链幅值LdLq是d-q轴电感对和PMSM的转矩公式几乎一模一样只是前面多了个(3π/pτ)的转换因子还要加个“直线电机特有的麻烦”项不对入门先别加比如齿槽力、端部力这些非线性先放一放等控制环跑通再喂进去第二步Simulink里的“直线版PMSM”速成别自己从零画电感模块、坐标变换模块太费头发MATLAB/Simulink里自带Motor Control BlocksetMCB或者Simscape Electrical Specialized Power Systems都有现成的PMSM子系统咱们只需要改几个参数就行% 先打开一个Simulink空白模型然后用代码帮咱们调下PMSM的参数 % 假设咱们搭的是一台超入门级的小推力PMLSM极对数p2极距τ20mm永磁体磁链ψ_f0.08WbL_dL_q0.001H忽略凸极效应入门更简单定子电阻R_s0.5Ω额定推力F_rated10N open_system(new_system(PMLSM_Basic)); add_block(mcb/PMSM, PMLSM_Basic/PMLSM); set_param(PMLSM_Basic/PMLSM, Ld, 0.001); set_param(PMLSM_Basic/PMLSM, Lq, 0.001); set_param(PMLSM_Basic/PMLSM, Rs, 0.5); set_param(PMLSM_Basic/PMLSM, FluxPM, 0.08); set_param(PMLSM_Basic/PMLSM, PolePairs, 2); % 最重要的把“机械转动惯量J”换成直线的“质量M”把“粘滞阻尼系数B”换成直线的“粘滞阻尼系数B_l” % 转动惯量和质量的转换哦对Simulink的PMSM模块如果是“机械端口选择为旋转机械”才需要直接加Simscape的“Rotational-Translational Interface”更简单 add_block(simscape/Foundation Library/Mechanical/Mechanical Sensors and Transducers/Rotational-Translational Interface, PMLSM_Basic/Rot2Trans); set_param(PMLSM_Basic/Rot2Trans, LeadScrewPitch, 4*0.02); % LeadScrewPitch其实就是2πτp不对入门咱们直接简化τ*2p让线速度v角速度ω*(τ*p)因为展开后的线速度和角速度的关系可以这么凑极对数p转一圈2π弧度走2p个极距也就是2pτ的距离所以没错vω*pτ那LeadScrewPitch就填pτ0.04m然后把剩下的FOC控制环Id0控制入门推力只靠q轴电流最直观搭起来转速环PI→电流环PID→SVPWM→逆变器→直线电机→Rot2Trans接口出来拉质量块。第三步跑个小测试爽爽咱们给位置指令加个1m/s的梯形速度曲线目标位置5m质量块设成0.1kg小质量更容易看出效果粘滞阻尼设成0.01N/(m/s)。跑出来的波形应该是什么样速度环先冲再稳推力和q轴电流成正比因为Id0而且忽略凸极和非线性位置曲线刚好在5m的时候停下误差大概在±0.001m如果不是±0.001m那就调一下位置环PI的参数比如把P调大一点别超调太多就行。第二站进阶篇喂进“真实的直线电机麻烦”仿真才不是纸上谈兵入门的那个模型跑出来太“完美”了真实的PMLSM有齿槽力、端部力这些推力波动还有摩擦力、负载扰动这些外部干扰如果不加这些优化出来的参数拿到实物上根本不能用。咱们今天就加两个最典型的齿槽力和库仑粘滞摩擦力。第一步齿槽力怎么加别硬算有限元的结果先用经验公式凑齿槽力是因为定子铁芯有齿有槽永磁体路过齿的时候吸力大路过槽的时候吸力小反复横跳产生的。入门级的齿槽力经验公式很简单% 齿槽力经验公式F_cog F_cog0*sin(2πp x / τ) % F_cog0是齿槽力幅值入门级的小推力PMLSM大概是额定推力的5%也就是0.5N % 先写个MATLAB Function放到Simulink里 function F_cog cogging_force(x) p 2; tau 0.02; F_cog0 0.5; F_cog F_cog0 * sin(2*pi*p*x / tau); end这个公式对不对大概对如果要更准的可以先做个简单的二维有限元分析Ansys Maxwell二维瞬态或者Harmonic把齿槽力随位置变化的曲线导出来再用MATLAB的“Curve Fitting Toolbox”拟合一下比如加个二次谐波三次谐波之类的但入门进阶先用这个正弦曲线就行。第二步库仑粘滞摩擦力怎么加Simulink里有现成的模块吗有Simscape Electrical Specialized Power Systems里的“Mechanical Load”模块或者Simscape里的“Translational Friction”模块都可以直接选“库仑粘滞”模型。咱们还是用Simscape的吧更真实% 把入门篇里的质量块换成Simscape的“Translational Mass”模块再加上“Translational Friction”模块 % 库仑摩擦力F_coul设成0.2N粘滞阻尼系数B_l还是0.01N/(m/s) add_block(simscape/Foundation Library/Mechanical/Translational Elements/Translational Friction, PMLSM_Basic/Friction); set_param(PMLSM_Basic/Friction, CoulombFriction, 0.2); set_param(PMLSM_Basic/Friction, ViscousFriction, 0.01);第三步喂进麻烦后的波形对比咱们再跑一遍刚才的1m/s梯形速度曲线对比一下加麻烦前后的推力和位置波形推力波形加之前是完美的梯形因为速度环是梯形电流环PI跟踪加之后推力波形有了小波动而且启动和停止的时候还有库仑摩擦力的尖峰位置波形加之前几乎没有误差加之后启动的时候有一点点滞后因为库仑摩擦力阻碍运动停止的时候可能会有一点点抖动因为推力波动和库仑摩擦力的换向如果要把这个抖动消掉可以试试滑模控制SMC或者模型预测控制MPC但这是下一站优化篇的内容了咱们先记下来这个问题。第三站精通篇调参数调结构让仿真结果无限接近实物需求现在咱们的模型“有血有肉”了接下来就是优化——要么优化控制参数PI、SMC、MPC的系数要么优化结构参数极距τ、永磁体厚度h、齿槽宽比bs/bt。咱们今天就分别优化一个控制参数用遗传算法GA调SMC的滑模面系数和趋近律系数结构参数用Ansys Maxwell二维瞬态Isight集成调τ和h。第一步控制参数优化——遗传算法调SMC的滑模系数Id0的PI控制对付不了推力波动和库仑摩擦力咱们先把控制环换成SMC然后用GA调系数。SMC的核心是滑模面s和趋近律咱们选最常用的线性滑模面和指数趋近律% 位置环的线性滑模面s c1*(x_ref - x) c2*(v_ref - v) % 指数趋近律ds/dt -ε*sign(s) - k*s % GA的目标函数是位置跟踪误差的积分绝对值IAEJ ∫|x_ref - x|dt % 约束条件是c10, c20, ε0, k0 % 先写个SMC的控制逻辑MATLAB Function再写个GA的优化脚本 % 优化脚本大概长这样简化版别直接复制粘贴要结合自己的模型 fitnessfcn (c) PMLSM_SMC_IAE(c); % c是[c1,c2,ε,k] nvars 4; lb [1,1,0.1,0.1]; % 下界 ub [100,100,10,10]; % 上界 options optimoptions(ga,Display,iter,MaxGenerations,50); [x,fval,exitflag,output] ga(fitnessfcn,nvars,[],[],[],[],lb,ub,[],options);把优化出来的c1,c2,ε,k放到SMC的控制逻辑里再跑一遍刚才的梯形速度曲线IAE会比PI控制小很多大概能小50%以上位置误差也会从±0.001m降到±0.0001m左右爽不爽第二步结构参数优化——Maxwell二维Isight集成调τ和h刚才的优化是“改算法改参数”现在咱们改“硬件结构”——目标是减小齿槽力幅值Fcog0因为齿槽力是推力波动的主要来源之一。咱们选两个最容易改的结构参数极距τ和永磁体厚度h目标函数是Fcog0最小约束条件是额定推力F_em≥10N总成本可以简化为永磁体体积VhτL*2p最小L是定子铁芯长度尽量小。怎么把Maxwell和Isight集成起来其实很简单Isight里有个“COMSOL Multiphysics”或者“Ansys Maxwell”的插件咱们先在Maxwell里做个二维瞬态模型设置好τ和h为设计变量设置好Fcog0和Fem为输出变量然后把这个模型导入Isight加上遗传算法或者粒子群算法PSO作为优化算法设置好约束条件和目标函数点击“Run”就行。永磁同步直线电机建模仿真及优化教学 从入门到精通永磁同步直线电机优化出来的结果大概是什么样比如τ可能从原来的20mm变成19.5mm刚好和齿槽宽凑成半整数齿槽配合齿槽力会大幅减小h可能从原来的5mm变成4mm减小永磁体体积但保证额定推力F_cog0会从原来的0.5N降到0.1N左右爽不爽最后唠两句从入门到精通PMLSM的建模仿真优化关键是先搭简单模型再慢慢加真实的麻烦最后调参数调结构。别一开始就啃FEA的黑箱子别一开始就用复杂的算法先用MATLAB/Simulink的集中参数模型爽爽等有了感觉再往上加东西。还有仿真只是辅助最终还是要做实物验证——因为仿真里的很多参数都是假设的比如定子电阻Rs会随温度变化永磁体磁链ψf会随温度和磁场变化这些都是仿真里很难完全考虑到的。好了今天的唠嗑就到这里有什么问题欢迎在评论区留言