发明名称 一种交流伺服系统速度环控制器参数自整定方法
摘要 本发明公开了一种交流伺服系统速度环控制器参数自整定方法,该方法首先从电流的反馈信号<img file="201410529147x100004dest_path_image001.GIF" wi="12" he="27" />和转速的反馈信号<img file="2819dest_path_image002.GIF" wi="14" he="18" />中,根据基于模型参考自适应的辨识机理辨识出当前系统的转动惯量<img file="dest_path_image003.GIF" wi="14" he="27" />;然后根据实际的动态要求以及执行器的输出性能确定性能指标<img file="dest_path_image005.GIF" wi="18" he="18" />;最后根据当前系统转动惯量的估计值以及性能指标计算出最优控制器的参数。该方法实现简单,参数调节较少,普适性强,对工况(转动惯量)变化较大的情况,有良好的适应性,能够满足高性能永磁同步电机交流伺服领域应用需求,适用于工况(转动惯量)变化剧烈的情况。
申请公布号 CN104242770A 申请公布日期 2014.12.24
申请号 CN201410529147.X 申请日期 2014.10.09
申请人 南京科远自动化集团股份有限公司;东南大学 发明人 李世华;颜赟达;吴斌;王军晓;王会明;刘国耀;沈德明;朱广斌;彭帅
分类号 H02P21/00(2006.01)I 主分类号 H02P21/00(2006.01)I
代理机构 南京汇盛专利商标事务所(普通合伙) 32238 代理人 陈扬
主权项 一种交流伺服系统速度环控制器参数自整定方法,其特征在于:该方法首先通过基于模型参考自适应辨识机理对永磁同步电机调速系统的转动惯量进行辨识,然后将辨识出来的惯量值代入系统的模型并建立性能指标,最后根据贝尔曼动态规划求解最优控制器;具体步骤如下:1)给定一个速度指令信号,根据下式,从q轴电流反馈信号i<sub>q</sub>以及转速信号ω中辨识出电机当前工况下的转动惯量;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>&omega;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mn>2</mn><mi>&omega;</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><mover><mi>b</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>T</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000583382820000011.GIF" wi="951" he="81" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>b</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>b</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>&beta;</mi><mfrac><mrow><msub><mi>&Delta;T</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><msubsup><mi>&beta;&Delta;T</mi><mi>e</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000583382820000012.GIF" wi="1410" he="155" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>J</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mrow><mover><mi>b</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000583382820000013.GIF" wi="264" he="140" /></maths>其中,ΔT<sub>e</sub>(k)=K<sub>t</sub>[i<sub>q</sub>(k)‑i<sub>q</sub>(k‑1)],<img file="FDA0000583382820000014.GIF" wi="458" he="73" />ω(k)是第k个采样周期的转速值,<img file="FDA0000583382820000015.GIF" wi="116" he="74" />是第k个采样周期的转速值的估计值,i<sub>q</sub>(k)是第k个采样周期的电流值,β为正常数;2)采用<img file="FDA0000583382820000016.GIF" wi="130" he="85" />的控制策略,得到速度环控制器;把电流环传递函数近似为常数1,系统的动力学方程可以表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>t</mi></msub><msubsup><mi>i</mi><mi>q</mi><mo>*</mo></msubsup></mrow><mi>J</mi></mfrac><mo>-</mo><mfrac><mi>B&omega;</mi><mi>J</mi></mfrac><mo>-</mo><mfrac><msub><mi>T</mi><mi>L</mi></msub><mi>J</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000583382820000017.GIF" wi="1140" he="149" /></maths>其中,ω为转速信号,<img file="FDA0000583382820000018.GIF" wi="46" he="85" />为q轴电流给定信号,J为转动惯量,B为粘滞摩擦系数,<img file="FDA0000583382820000019.GIF" wi="288" he="129" />n<sub>p</sub>为磁极对数,ψ<sub>f</sub>为转子磁链,系统的离散化模型可以表示为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>BT</mi><mi>s</mi></msub><mi>J</mi></mfrac><mo>)</mo></mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mi>t</mi></msub><msub><mi>T</mi><mi>s</mi></msub></mrow><mi>J</mi></mfrac><msubsup><mi>i</mi><mi>q</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>J</mi></mfrac><msub><mi>T</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005833828200000110.GIF" wi="1519" he="148" /></maths>其中,T<sub>s</sub>为采样时间,系统的状态空间表达式为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>ax</mi><mo>+</mo><mi>bu</mi></mrow>]]></math><img file="FDA00005833828200000115.GIF" wi="263" he="74" /></maths>                           (23)y=cx其中,x为状态,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>x</mi><mi>T</mi></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mi>&omega;</mi></mtd><mtd><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA00005833828200000112.GIF" wi="300" he="83" /></maths>y为输出量,y=ω;u为控制量,<img file="FDA00005833828200000113.GIF" wi="147" he="86" /><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><mi>B</mi><mi>J</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00005833828200000114.GIF" wi="319" he="209" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><msub><mi>K</mi><mi>t</mi></msub><mi>J</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000583382820000021.GIF" wi="231" he="218" /></maths>c=(1 0);将步骤1)中得到的惯量的辨识值,代入(23)式,得到当前工况下的系统模型;建立性能指标H,见下式;<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mo>&infin;</mo></msubsup><mrow><mo>(</mo><msup><mi>Qe</mi><mn>2</mn></msup><mo>+</mo><mi>R</mi><msubsup><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>q</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>)</mo></mrow><mi>dt</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000583382820000022.GIF" wi="1278" he="132" /></maths>其中,e=ω<sub>r</sub>‑ω,ω<sub>r</sub>为转速的给定信号,Q、R为正常数;根据贝尔曼动态规划求解出此时的最优控制量,见下式:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>i</mi><mi>q</mi><mo>*</mo></msubsup><mo>=</mo><msub><mi>m</mi><mn>1</mn></msub><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mi>&omega;dt</mi><mo>+</mo><msub><mi>m</mi><mn>2</mn></msub><mi>&omega;</mi><mo>+</mo><mi>n</mi><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>&omega;</mi><mi>r</mi></msub><mi>dt</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000583382820000023.GIF" wi="1317" he="109" /></maths>其中,P为满足黎卡提代数方程Pa+a<sup>T</sup>P‑PbR<sup>‑1</sup>b<sup>T</sup>P+c<sup>T</sup>Qc=0的对称正定常数矩阵,m=‑R<sup>‑1</sup>B<sup>T</sup>P=(m<sub>1</sub> m<sub>2</sub>),n=R<sup>‑1</sup>b<sup>T</sup>(PbR<sup>‑1</sup>b<sup>T</sup>‑a<sup>T</sup>)<sup>‑1</sup>c<sup>T</sup>Q。
地址 211100 江苏省南京市江宁区秣陵街道清水亭东路1266号