发明名称 一种永磁同步电机的H无穷转速估计方法
摘要 永磁同步电机的H无穷转速估计方法,包括电机的矢量控制系统,所述矢量控制系统反馈通道上设有一个H无穷转速观测器,通过粒子群算法实现H无穷转速观测器参数的自主寻优,获得转速估计值。本发明优化了传统转速观测器易受系统参数变化的影响、算法复杂难于实现和抗干扰能力弱等缺点,依据H无穷控制理论实现了永磁同步电机的转速估计,实现了粒子群算法H无穷速度观测器与机械速度传感器的完美替换,从而降低了系统成本,且保证了系统的稳定性,有效改善了永磁同步电机的动态性能。
申请公布号 CN103647493B 申请公布日期 2015.12.02
申请号 CN201310567325.3 申请日期 2013.11.14
申请人 西安理工大学 发明人 钟彦儒;朱群;尹忠刚
分类号 H02P21/14(2006.01)I;H02P21/13(2006.01)I;H02P25/02(2006.01)I 主分类号 H02P21/14(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 李娜
主权项 永磁同步电机的H无穷转速估计方法,其特征在于,包括电机的矢量控制系统,所述矢量控制系统反馈通道上设有一个H无穷转速观测器,通过粒子群算法实现H无穷转速观测器参数的离线寻优,获得转速估计值,具体包括如下步骤:步骤1,通过电机的电磁转矩和机械运动方程推导得到永磁同步电机数学模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>A</mi><mi>x</mi><mo>+</mo><msub><mi>B</mi><mn>1</mn></msub><mi>w</mi><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mi>v</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><msub><mi>C</mi><mn>2</mn></msub><mi>x</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>z</mi><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><mi>x</mi><mo>+</mo><msub><mi>D</mi><mn>12</mn></msub><mi>v</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000797988680000011.GIF" wi="520" he="259" /></maths>式中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&omega;</mi><mrow><mi>r</mi><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>&theta;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&epsiv;</mi><mi>&theta;</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><msub><mi>B</mi><mi>m</mi></msub><mo>/</mo><msub><mi>J</mi><mi>m</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000797988680000012.GIF" wi="758" he="295" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>B</mi><mn>1</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>J</mi><mi>m</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>B</mi><mn>2</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>K</mi><mi>t</mi></msub><mo>/</mo><msub><mi>J</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo><mi>w</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>T</mi><mi>L</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>&theta;</mi><mrow><mi>r</mi><mi>m</mi></mrow><mo>*</mo></msubsup></mtd></mtr><mtr><mtd><msubsup><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mrow><mi>r</mi><mi>m</mi></mrow><mo>*</mo></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>C</mi><mn>2</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>C</mi><mn>1</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>q</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>q</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>q</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000797988680000013.GIF" wi="1871" he="307" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>D</mi><mn>12</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>r</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000797988680000014.GIF" wi="243" he="326" /></maths>其中,矩阵A、B<sub>1</sub>、B<sub>2</sub>为系数矩阵,e<sub>θ</sub>是θ<sub>rm</sub><sup>*</sup>和θ<sub>rm</sub>之间的位置误差,θ<sub>rm</sub><sup>*</sup>为给定参考量,ε<sub>θ</sub>是位置误差的积分,v为控制输入信号,x为状态变量,<img file="FDA0000797988680000015.GIF" wi="48" he="58" />为状态变量x的一阶导数,w是干扰输入信号,变量z为评价信号,C<sub>1</sub>和D<sub>12</sub>为加权矩阵,q<sub>i</sub>&gt;0(i=1,2,3)和r&gt;0为加权系数;C<sub>1</sub>用来调节干扰抑制效果,D<sub>12</sub>决定控制输入的大小;步骤2,根据步骤1中永磁同步电机数学模型搭建H无穷转速观测器结构,反馈控制率由下式得到,也就是H无穷转速观测器的给定表达式:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>v</mi><mo>*</mo></msup><mo>=</mo><mi>K</mi><mi>x</mi><mo>=</mo><mo>-</mo><mrow><mo>&lsqb;</mo><mrow><msup><mrow><mo>(</mo><msup><msub><mi>D</mi><mn>12</mn></msub><mi>T</mi></msup><msub><mi>D</mi><mn>12</mn></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><msub><mi>D</mi><mn>12</mn></msub><mi>T</mi></msup><msub><mi>C</mi><mn>1</mn></msub><mo>+</mo><msup><mrow><mo>(</mo><msup><msub><mi>D</mi><mn>12</mn></msub><mi>T</mi></msup><msub><mi>D</mi><mn>12</mn></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mi>B</mi><mn>2</mn><mi>T</mi></msubsup><mi>P</mi></mrow><mo>&rsqb;</mo></mrow><mi>x</mi><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><msup><mi>r</mi><mn>2</mn></msup></mfrac><mo>&lsqb;</mo><mtable><mtr><mtd><mrow><msub><mi>K</mi><mi>t</mi></msub><mo>/</mo><msub><mi>J</mi><mi>m</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable><mo>&rsqb;</mo><mi>P</mi><mi>x</mi></mrow>]]></math><img file="FDA0000797988680000021.GIF" wi="1951" he="170" /></maths>其中,K为反馈控制率,D<sub>12</sub><sup>T</sup>为D<sub>12</sub>的转置,矩阵P为Riccati方程的解,γ为给定的上界,K<sub>t</sub>为转矩常数,J<sub>m</sub>为转子的转动惯量;步骤3,步骤2中的给定表达式中含有C<sub>1</sub>和D<sub>12</sub>当中4个可调参数,将粒子群算法引入H无穷转速观测器中参数的求解,利用粒子群算法对这4个可调参数进行离线寻优,其步骤为:a.粒子的位置:x=(q<sub>1</sub>,q<sub>2</sub>,q<sub>3</sub>,r)<sup>T</sup>;b.选取适应度函数:<img file="FDA0000797988680000022.GIF" wi="451" he="125" />式中,e(t)为给定转速与H无穷转速观测器估计转速的差值;c.保留每个粒子的历史最优位置,以及所有粒子的全局最优位置;d.根据历史最优和全局最优更新粒子的速度和位置;每个粒子的速度和位置更新后,对更新后的位置中的每个变量进行取最近整数的操作;e.如果达到预设的迭代次数,则执行下一步;否则返回b;f.输出gbest和pbest,算法运行结束;g.输出gbest对应的H无穷观测器,该观测器即为所求;步骤4,将给定转速ω<sub>rm</sub><sup>*</sup>与估计转速ω<sub>rm</sub>作差得到的误差信号Δω<sub>rm</sub>送入PI调节器,得到系统转矩值,然后采用最大转矩电流比控制器,得到d、q轴电流的期望给定值;步骤5,q轴电流的期望给定值与输出电流作差,得到电流误差Δi<sub>q</sub>,同理得到Δi<sub>d</sub>,经过PI调节器输出u<sub>sd</sub><sup>*</sup>、u<sub>sq</sub><sup>*</sup>,u<sub>sd</sub><sup>*</sup>、u<sub>sq</sub><sup>*</sup>再经park反变换输出u<sub>sα</sub><sup>*</sup>、u<sub>sβ</sub><sup>*</sup>,最后通过空间矢量脉宽调制模块输出六路PWM信号供给逆变器工作,逆变器输出将直流母线电压V<sub>dc</sub>以PWM波的形式将电压施加到永磁同步电机上。
地址 710048 陕西省西安市金花南路5号