发明名称 永磁同步电机多参数解耦在线辨识方法
摘要 本发明涉及永磁同步电机技术领域,为解决表贴式永磁同步电机多参数在线辨识时的耦合问题,实现永磁同步电机电感、定子电阻和转子磁链在线解耦辨识,为此,本发明采取的技术方案是,永磁同步电机多参数解耦在线辨识方法,包括下列步骤:(1)永磁同步电机参数辨识耦合分析;(2)解耦辨识策略:利用d轴电流注入前后的电压偏差来增加电机数学方程的阶数,以实现表贴式永磁同步电机电感、定子电阻、转子磁链多参数解耦辨识;(3)神经网络辨识器设计,针对永磁同步电机参数在线辨识问题,选取自适应神经网络结构和基于最小均方算法的权值收敛算法对电机参数进行在线辨识。本发明主要应用于永磁同步电机的设计制造。
申请公布号 CN103248306B 申请公布日期 2015.05.13
申请号 CN201310134864.8 申请日期 2013.05.24
申请人 天津大学 发明人 史婷娜;胡升;阎彦;夏长亮
分类号 H02P21/14(2006.01)I;H02P27/08(2006.01)I 主分类号 H02P21/14(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 刘国威
主权项 一种永磁同步电机多参数解耦在线辨识方法,其特征是,包括如下步骤:(1)永磁同步电机参数辨识耦合分析忽略铁心饱和;不计涡流和磁滞损耗;转子上无阻尼绕组,永磁体也无阻尼作用;反电动势为正弦;永磁同步电机在d‑q轴坐标系的电流方程通常可写成如下形式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>L</mi><mi>d</mi></msub><mfrac><msub><mi>di</mi><mi>d</mi></msub><mi>dt</mi></mfrac><mo>=</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>L</mi><mi>q</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mo>+</mo><msub><mi>u</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>q</mi></msub><mfrac><msub><mi>di</mi><mi>q</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>&psi;</mi><mi>f</mi></msub><mo>+</mo><msub><mi>u</mi><mi>q</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000011.GIF" wi="549" he="295" /></maths>式中,u<sub>d</sub>、u<sub>q</sub>、i<sub>d</sub>、i<sub>q</sub>分别为d、q轴电压、电流;L<sub>d</sub>、L<sub>q</sub>为d、q轴电感;R<sub>s</sub>为定子电阻;ω<sub>e</sub>为电气角速度;ψ<sub>f</sub>为转子磁链,采用i<sub>d</sub>=0;其中,表贴式永磁同步电机通常认为交直轴电感相等,即L<sub>d</sub>=L<sub>q</sub>=L<sub>s</sub>;(2)解耦辨识策略利用d轴电流注入前后的电压偏差来增加电机数学方程的阶数,以实现表贴式永磁同步电机电感、定子电阻、转子磁链多参数解耦辨识,短时注入d轴电流前后,永磁同步电机的电压方程可综合写成如下形式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>L</mi><mi>s</mi></msub><mfrac><mrow><msub><mi>di</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>L</mi><mi>q</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>s</mi></msub><mfrac><mrow><msub><mi>di</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>&psi;</mi><mi>f</mi></msub><mo>+</mo><msub><mi>u</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>s</mi></msub><mfrac><msub><mi>di</mi><mi>d</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>L</mi><mi>q</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>s</mi></msub><mfrac><mrow><msub><mi>di</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>L</mi><mi>d</mi></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>u</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>&psi;</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000012.GIF" wi="1131" he="559" /></maths>式中,t<sub>2</sub>表示注入瞬时d轴电流时的某一时刻,t<sub>3</sub>表示d轴瞬时电流注入结束后的某一时刻;利用上述数学方程,实时记录瞬时d轴电流注入时的d轴电压值,检测其前后变化量,结合i<sub>d</sub>=0正常运行时的电压方程,即可实现电感、定子电阻和转子磁链参数的独立辨识;(3)神经网络辨识器设计针对永磁同步电机参数在线辨识问题,选取自适应神经网络结构和基于最小均方算法的权值收敛算法对电机参数进行在线辨识,假设Z为实际系统模型的输出,α为权值收敛因子,Y为神经网络输出,则最小均方权值收敛算法的表达式为:θ<sub>i</sub>(k+1)=θ<sub>i</sub>(k)+2αU<sub>i</sub>(Z‑Y)θ<sub>i</sub>为网络权值;U<sub>i</sub>为输入信号;当系统在线稳定运行时,d、q轴电流微分项近似为0,则当系统处于d轴电流瞬时注入时的稳态时,表贴式永磁同步电机电压方程可进一步简化为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>L</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>&psi;</mi><mi>f</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>L</mi><mi>s</mi></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000021.GIF" wi="729" he="162" /></maths>当注入d轴电流结束后,i<sub>d</sub>=0,则系统稳态时的电压方程可进一步简化为下式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>L</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>&psi;</mi><mi>f</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000022.GIF" wi="487" he="175" /></maths>根据上述两式中的d轴电压方程可得注入瞬时d轴电流前后的d轴电压偏差方程为:u<sub>d</sub>(t<sub>3</sub>)‑u<sub>d</sub>(t<sub>2</sub>)=Δu<sub>d</sub>=‑R<sub>s</sub>i<sub>d</sub>(t<sub>2</sub>)利用以上三式可分别建立电感、定子电阻和转子磁链辨识网络结构,实现电感、定子电阻和转子磁链的独立辨识;神经网络辨识器设计进一步具体为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>L</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>u</mi><mi>q</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mn>3</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>&psi;</mi><mi>f</mi></msub><mo>+</mo><msub><mi>&theta;</mi><mn>2</mn></msub><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000023.GIF" wi="795" he="186" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Delta;u</mi><mi>d</mi><mo>*</mo></msubsup><mo>=</mo><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&theta;</mi><mn>3</mn></msub><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000687346960000024.GIF" wi="829" he="78" /></maths>利用上述两式建立改进型神经网络辨识结构,θ<sub>1</sub>、θ<sub>2</sub>、θ<sub>3</sub>分别为考虑了逆变器输出增益的参数辨识校正权值,其权值收敛率满足下式,<maths num="0007" id="cmaths0007"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>G</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>G</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msubsup><mi>u</mi><mi>q</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>G</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msubsup><mi>&Delta;u</mi><mi>d</mi><mo>*</mo></msubsup></mrow><mrow><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000687346960000025.GIF" wi="445" he="493" /></maths>G<sub>u</sub>为等效的电压型逆变器基波输出增益,0&lt;G<sub>u</sub>&lt;1,则永磁同步电机电感、定子电阻和转子磁链的参数迭代收敛率进一步可写成如下形式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>L</mi><mi>s</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>L</mi><mi>s</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&omega;</mi><mi>e</mi></msub><msub><mi>i</mi><mi>q</mi></msub><mo>[</mo><msubsup><mi>u</mi><mi>d</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>u</mi><mi>d</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000687346960000026.GIF" wi="835" he="87" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>R</mi><mi>s</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>R</mi><mi>s</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>&alpha;</mi><mn>2</mn></msub><mi></mi><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>[</mo><msubsup><mrow><mo>(</mo><mi>&Delta;u</mi></mrow><mi>d</mi><mo>*</mo></msubsup><mo>-</mo><msubsup><mi>&Delta;u</mi><mi>d</mi><mo>&prime;</mo></msubsup><mo>]</mo></mrow>]]></math><img file="FDA0000687346960000027.GIF" wi="813" he="91" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>&psi;</mi><mi>f</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&psi;</mi><mi>f</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>&alpha;</mi><mn>3</mn></msub><msub><mi>&omega;</mi><mi>e</mi></msub><mo>[</mo><msubsup><mi>u</mi><mi>q</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>u</mi><mi>q</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000687346960000028.GIF" wi="823" he="91" /></maths>采用上述三式的参数收敛算法即可实现表贴式永磁同步电机的电感、定子电阻和转子磁链参数的解耦在线辨识。
地址 300072 天津市南开区卫津路92号