发明名称 一种基于自抗扰和锁相环技术的永磁同步电机系统无传感器速度控制方法
摘要 一种基于自抗扰和锁相环技术的永磁同步电机系统无传感器速度控制方法,包括:建立永磁同步电机系统的数学模型,初始化系统状态及控制参数;设计扩张状态观测器,估计电机的反电动势;基于三角函数,设计锁相环系统,通过估计的反电动势提取电机转子的速度和位置信号,并基于自抗扰设计电机的速度控制器。本发明能够有效解决永磁同步电机系统的复杂不稳定问题,简化电机硬件设计,实现永磁同步电机系统的无传感器稳定控制。
申请公布号 CN105827160A 申请公布日期 2016.08.03
申请号 CN201610158167.X 申请日期 2016.03.18
申请人 浙江工业大学 发明人 陈强;董方;陶亮;郭一军
分类号 H02P6/08(2016.01)I;H02P6/182(2016.01)I 主分类号 H02P6/08(2016.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种基于自抗扰和锁相环技术的永磁同步电机无传感器速度控制方法,其特征在于:包括以下步骤:步骤1,建立永磁同步电机系统数学模型,初始化系统状态及控制参数,过程如下:1.1,永磁同步电机的数学模型表示如下:<maths num="0001"><math><![CDATA[<mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>&alpha;</mi></msub><mo>=</mo><mo>-</mo><mfrac><mi>R</mi><mi>L</mi></mfrac><msub><mi>i</mi><mi>&alpha;</mi></msub><mo>-</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><msub><mi>e</mi><mi>&alpha;</mi></msub><mo>+</mo><mfrac><msub><mi>u</mi><mi>&alpha;</mi></msub><mi>L</mi></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>&beta;</mi></msub><mo>=</mo><mo>-</mo><mfrac><mi>R</mi><mi>L</mi></mfrac><msub><mi>i</mi><mi>&beta;</mi></msub><mo>-</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><msub><mi>e</mi><mi>&beta;</mi></msub><mo>+</mo><mfrac><msub><mi>u</mi><mi>&beta;</mi></msub><mi>L</mi></mfrac></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mrow>]]></math><img file="FDA0000944462510000011.GIF" wi="1150" he="270" /></maths>其中,i<sub>α</sub>,i<sub>β</sub>分别为α轴和β轴上的定子电流;u<sub>α</sub>,u<sub>β</sub>分别为α轴和β轴上的定子电压;R为定子的等效电阻;e<sub>α</sub>,e<sub>β</sub>分别为α‑β轴上的反电动势,表示如下:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>e</mi><mi>&alpha;</mi></msub><mo>=</mo><mo>-</mo><msub><mi>l</mi><mi>e</mi></msub><mi>&omega;</mi><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>e</mi><mi>&beta;</mi></msub><mo>=</mo><msub><mi>l</mi><mi>e</mi></msub><mi>&omega;</mi><mi>cos</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000012.GIF" wi="1123" he="165" /></maths>其中,l<sub>e</sub>为反电动势系数;θ和ω分别是转子的速度和位置信号;1.2,式(1)被转化为:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>f</mi><mo>(</mo><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo><mo>+</mo><mi>b</mi><mo>.</mo><msub><mi>u</mi><mi>o</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000013.GIF" wi="1124" he="62" /></maths>其中x<sub>1</sub>为[i<sub>α</sub> i<sub>β</sub>]<sup>T</sup>,x<sub>2</sub>表示<img file="FDA0000944462510000014.GIF" wi="217" he="87" />u<sub>o</sub>为[u<sub>α</sub> u<sub>β</sub>]<sup>T</sup>,b是控制器的增益并满足b=1/L;d(t)为系统总的扰动;根据式(3),设计以下的反电动势观测器来估计反电动势e<sub>α</sub>和e<sub>β</sub>;<maths num="0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>e</mi><mn>1</mn></msub><mo>=</mo><msub><mi>z</mi><mn>01</mn></msub><mo>-</mo><msub><mi>i</mi><mi>&alpha;</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>01</mn></msub><mo>=</mo><mi>b</mi><mo>.</mo><mrow><mo>(</mo><mrow><msub><mi>u</mi><mi>&alpha;</mi></msub><mo>-</mo><msub><mi>z</mi><mn>01</mn></msub><mo>.</mo><mi>R</mi><mo>-</mo><msub><mi>z</mi><mn>02</mn></msub></mrow><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>01</mn></msub><mo>.</mo><msub><mi>e</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>02</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&beta;</mi><mn>02</mn></msub><mo>.</mo><msub><mi>e</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>e</mi><mn>2</mn></msub><mo>=</mo><msub><mi>z</mi><mn>11</mn></msub><mo>-</mo><msub><mi>i</mi><mi>&beta;</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>11</mn></msub><mo>=</mo><mi>b</mi><mo>.</mo><mrow><mo>(</mo><mrow><msub><mi>u</mi><mi>&beta;</mi></msub><mo>-</mo><msub><mi>z</mi><mn>11</mn></msub><mo>.</mo><mi>R</mi><mo>-</mo><msub><mi>z</mi><mn>12</mn></msub></mrow><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>11</mn></msub><mo>.</mo><msub><mi>e</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>12</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&beta;</mi><mn>12</mn></msub><mo>.</mo><msub><mi>e</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000015.GIF" wi="1180" he="463" /></maths>其中,z<sub>01</sub>和z<sub>11</sub>分别是i<sub>α</sub>和i<sub>β</sub>的观测值;e<sub>1</sub>和e<sub>2</sub>分别是估计误差;z<sub>02</sub>和z<sub>12</sub>都是反电动势的估计值;β<sub>01</sub>,β<sub>02</sub>都是正常量;1.3,根据式(2),获得:<maths num="0005"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>&alpha;</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>e</mi><mi>&beta;</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>l</mi><mi>e</mi><mn>2</mn></msubsup><msup><mi>&omega;</mi><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000021.GIF" wi="996" he="70" /></maths>并且有:<maths num="0006"><math><![CDATA[<mrow><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>-</mo><mover><mi>&theta;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mi>c</mi><mi>o</mi><mi>s</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>-</mo><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mi>s</mi><mi>i</mi><mi>n</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000022.GIF" wi="1172" he="71" /></maths>其中,<img file="FDA0000944462510000023.GIF" wi="35" he="62" />为转子位置θ的估计值;根据式(5)和式(6),得到:<maths num="0007"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>e</mi><mi>&alpha;</mi></msub><mo>.</mo><mi>cos</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>+</mo><msub><mi>e</mi><mi>&beta;</mi></msub><mo>.</mo><mi>sin</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>=</mo><mo>-</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><mi>&omega;</mi><mo>.</mo><mi>sin</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>.</mo><mi>cos</mi><mover><mi>&theta;</mi><mo>^</mo></mover><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><mi>&omega;</mi><mo>.</mo><mi>cos</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>.</mo><mi>sin</mi><mover><mi>&theta;</mi><mo>^</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><mi>&omega;</mi><mo>.</mo><mi>sin</mi><mrow><mo>(</mo><mrow><mi>&theta;</mi><mo>-</mo><mover><mi>&theta;</mi><mo>^</mo></mover></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><mi>&omega;</mi><mo>.</mo><mi>sin</mi><mrow><mo>(</mo><mover><mi>&theta;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000024.GIF" wi="1422" he="254" /></maths>其中,<img file="FDA0000944462510000025.GIF" wi="38" he="58" />是估计误差,当<img file="FDA0000944462510000026.GIF" wi="37" he="58" />变得很小的时候,得到<img file="FDA0000944462510000027.GIF" wi="235" he="70" />经过锁相环的调节之后,最终得到<img file="FDA0000944462510000028.GIF" wi="158" he="62" />因此<img file="FDA0000944462510000029.GIF" wi="115" he="63" />成立;1.4,锁相环的传递函数和闭环误差传递函数如下所示:<maths num="0008"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mi>p</mi></msub><mo>.</mo><mi>s</mi><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mi>I</mi></msub></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mi>p</mi></msub><mo>.</mo><mi>s</mi><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mi>I</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009444625100000210.GIF" wi="1166" he="143" /></maths><maths num="0009"><math><![CDATA[<mrow><msub><mi>G</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><mi>e</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>e</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><msup><mi>s</mi><mn>2</mn></msup><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mrow><mi>p</mi><mo>.</mo></mrow></msub><mo>.</mo><mi>s</mi><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><mo>.</mo><msub><mi>k</mi><mi>I</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009444625100000211.GIF" wi="1222" he="142" /></maths>由于转子的位置信号是阶跃信号,因此锁相环的稳态误差传递函数为:<maths num="0010"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>e</mi><mrow><mo>(</mo><mi>&infin;</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>lim</mi><mrow><mi>s</mi><mo>&RightArrow;</mo><mn>0</mn></mrow></munder><mo>.</mo><mi>&Delta;</mi><mi>e</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>lim</mi><mrow><mi>s</mi><mo>&RightArrow;</mo><mn>0</mn></mrow></munder><mfrac><msup><mi>s</mi><mn>2</mn></msup><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><msub><mi>k</mi><mi>p</mi></msub><mi>s</mi><mo>+</mo><msub><mi>l</mi><mi>e</mi></msub><msub><mi>k</mi><mi>p</mi></msub></mrow></mfrac><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009444625100000212.GIF" wi="1269" he="142" /></maths>步骤2,基于自抗扰的速度控制器设计,过程如下:2.1,电机的运动学模型建立如下:<maths num="0011"><math><![CDATA[<mrow><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mrow><mn>3</mn><msub><mi>n</mi><mi>p</mi></msub><msub><mi>&psi;</mi><mi>f</mi></msub></mrow><mrow><mn>2</mn><mi>J</mi></mrow></mfrac><msub><mi>i</mi><mi>q</mi></msub><mo>-</mo><mfrac><mi>B</mi><mi>J</mi></mfrac><mi>&omega;</mi><mo>-</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><msubsup><mi>bi</mi><mi>q</mi><mo>*</mo></msubsup><mo>+</mo><mi>a</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009444625100000213.GIF" wi="1286" he="126" /></maths>其中,a(t)为系统的总扰动;系统的状态方程为:<maths num="0012"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>b</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>a</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><mi>&omega;</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009444625100000214.GIF" wi="1118" he="142" /></maths>2.2,线性扩张状态观测器被设计为:<maths num="0013"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>e</mi><mn>0</mn></msub><mo>=</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&beta;</mi><mn>21</mn></msub><mo>.</mo><msub><mi>e</mi><mn>0</mn></msub><mo>+</mo><mi>b</mi><mo>.</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&beta;</mi><mn>22</mn></msub><mo>.</mo><msub><mi>e</mi><mn>0</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000031.GIF" wi="1118" he="223" /></maths>其中,z<sub>1</sub>为ω(t)的观测值;z<sub>2</sub>为系统总扰动的估计值;β<sub>21</sub>和β<sub>22</sub>都是正常量;2.3,非线性误差反馈控制器设计如下:u(t)=(k(ω<sub>m</sub>‑ω)‑z<sub>2</sub>)/b    (14)其中,u(t)为电机的速度控制器输入信号;ω<sub>m</sub>为电机速度的参考信号;b为电机速度的参考信号;2.4,为了简化线性扩张状态观测器,u被设计为:<maths num="0014"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mfrac><mrow><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>+</mo><msub><mi>u</mi><mn>0</mn></msub></mrow><mi>b</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000032.GIF" wi="1077" he="110" /></maths>忽略系统扰动的影响,式(12)被简化为:<maths num="0015"><math><![CDATA[<mrow><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo><mo>+</mo><msub><mi>u</mi><mn>0</mn></msub><mo>&ap;</mo><msub><mi>u</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000033.GIF" wi="1166" he="63" /></maths>根据式(12)和式(16),系统的闭环传递函数为:<maths num="0016"><math><![CDATA[<mrow><msub><mi>G</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>k</mi><mrow><mi>s</mi><mo>+</mo><mi>k</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944462510000034.GIF" wi="1094" he="119" /></maths>通过选择k,速度控制系统能够保持稳定。
地址 310014 浙江省杭州市下城区朝晖六区潮王路18号浙江工业大学