发明名称 一种感应电机最大转矩电流比控制方法
摘要 本发明公开了一种感应电机最大转矩电流比控制方法,属于感应电机控制的技术领域。所述控制方法在感应电机处于轻载运行的工况下,将PI和预测函数控制相结合,得到一种可应用于多输入多输出系统的新的控制方法,把预测函数控制器输出的控制量作为逆变器电路的输入信号,来改变逆变器的输出电压,从而控制电机转速,实现对感应电机转速的跟踪控制。控制精度高、跟踪速度较快、稳态误差小,抗干扰能力强,既能保证感应电机具有较好的稳定性与动态性能,又能提高电机动态时的运行效率。
申请公布号 CN104135205B 申请公布日期 2016.07.13
申请号 CN201410320017.5 申请日期 2014.07.08
申请人 南京信息工程大学 发明人 郭伟;夏友亮;王汉杰;郁振波;韩丹丹
分类号 H02P21/22(2016.01)I 主分类号 H02P21/22(2016.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 杨楠
主权项 一种感应电机最大转矩电流比控制方法,包括速度环和电流环,所述速度环为PI控制,其特征在于:所述电流环为多变量PI预测函数控制,所述多变量PI预测函数控制的步骤如下:步骤1、初始化感应电机控制系统的各项参数:定子电阻R<sub>s</sub>,转子电阻R<sub>r</sub>,定子自感L<sub>s</sub>,转子自感L<sub>r</sub>,互感L<sub>m</sub>,转动惯量J,额定频率f,转子极对数n<sub>p</sub>,额定转速n<sub>N</sub>,电机同步旋转速度ω<sub>a</sub>,建立基于转子磁场定向的感应电机动态数学模型,确定上述模型的参考轨迹,并将基于转子磁场定向的感应电机动态数学模型转化为状态空间方程,得出系数矩阵A<sub>m</sub>、B<sub>m</sub>、C<sub>m</sub>;所述基于转子磁场定向的感应电机动态数学模型的表达式为:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mrow><mi>s</mi><mi>d</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mrow><mi>s</mi><mi>q</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><mo>/</mo><msub><mi>&sigma;L</mi><mi>s</mi></msub></mrow></mtd><mtd><msub><mi>&omega;</mi><mi>a</mi></msub></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>a</mi></msub><mo>/</mo><mi>&sigma;</mi></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><mo>/</mo><msub><mi>&sigma;L</mi><mi>s</mi></msub></mrow></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>i</mi><mrow><mi>s</mi><mi>d</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mrow><mi>s</mi><mi>q</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mn>1</mn><mo>/</mo><msub><mi>&sigma;L</mi><mi>s</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mn>1</mn><mo>/</mo><msub><mi>&sigma;L</mi><mi>s</mi></msub></mrow></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>u</mi><mrow><mi>s</mi><mi>d</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>s</mi><mi>q</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000960605360000011.GIF" wi="1125" he="150" /></maths><img file="FDA0000960605360000012.GIF" wi="349" he="151" />其中,σ=1‑L<sub>m</sub><sup>2</sup>/L<sub>s</sub>L<sub>r</sub>,i<sub>sd</sub>为定子电流直轴分量,i<sub>sq</sub>为定子电流交轴分量,u<sub>sd</sub>为定子电压直轴分量,u<sub>sq</sub>为定子电压交轴分量;所述状态空间方程的表达式为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>X</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mi>m</mi></msub><msub><mi>X</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mi>m</mi></msub><mi>U</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Y</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mi>m</mi></msub><msub><mi>X</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000960605360000013.GIF" wi="732" he="150" /></maths>其中,Y<sub>m</sub>(k)为k时刻模型预测输出向量,X<sub>m</sub>(k)、X<sub>m</sub>(k‑1)分别为k、k‑1时刻模型状态值向量,U(k‑1)为k‑1时刻控制输入向量;步骤2、根据下式计算控制量向量U(k):U(k)=(L<sub>a</sub>+L<sub>b</sub>+L<sub>c</sub>)F<sub>n</sub>(0)<sup>T</sup>,其中:<maths num="0003"><math><![CDATA[<mrow><msub><mi>L</mi><mi>a</mi></msub><mo>=</mo><mo>-</mo><mo>&lsqb;</mo><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><msub><mi>QGF</mi><mi>n</mi></msub><mo>+</mo><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msub><mi>RF</mi><mi>n</mi></msub><mo>&rsqb;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>K</mi><mi>P</mi></msub><mo>+</mo><msub><mi>K</mi><mi>i</mi></msub><mo>)</mo></mrow><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><mi>Q</mi><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000014.GIF" wi="1110" he="71" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>L</mi><mi>b</mi></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><msub><mi>QGF</mi><mi>n</mi></msub><mo>+</mo><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msub><mi>RF</mi><mi>n</mi></msub><mo>&rsqb;</mo><mo>&CenterDot;</mo><mn>2</mn><msub><mi>K</mi><mi>P</mi></msub><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><msup><mi>Qq</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000015.GIF" wi="1038" he="71" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>L</mi><mi>c</mi></msub><mo>=</mo><mo>-</mo><mo>&lsqb;</mo><msub><mi>K</mi><mi>i</mi></msub><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><msub><mi>QGF</mi><mi>n</mi></msub><mo>+</mo><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msub><mi>RF</mi><mi>n</mi></msub><mo>&rsqb;</mo><mo>&CenterDot;</mo><msub><mi>K</mi><mi>P</mi></msub><msubsup><mi>F</mi><mi>n</mi><mi>T</mi></msubsup><msup><mi>G</mi><mi>T</mi></msup><msup><mi>Qq</mi><mrow><mo>-</mo><mn>2</mn></mrow></msup><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000016.GIF" wi="1041" he="71" /></maths>F<sub>n</sub>=[f<sub>1</sub>(i) f<sub>2</sub>(i) … f<sub>J</sub>(i)],i=1,2,…,P‑1,则F<sub>n</sub>(0)=[f<sub>1</sub>(0) f<sub>2</sub>(0) … f<sub>J</sub>(0)],<maths num="0006"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><msup><mfenced open = "(" close = ")"><mtable><mtr><mtd><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup></mtd><mtd><msubsup><mi>G</mi><mn>2</mn><mi>T</mi></msubsup></mtd><mtd><mo>...</mo></mtd><mtd><msubsup><mi>G</mi><mi>p</mi><mi>T</mi></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo><msub><mi>G</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>m</mi></msub><msubsup><mi>A</mi><mi>m</mi><mrow><mi>P</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>m</mi></msub><mo>+</mo><msub><mi>C</mi><mi>m</mi></msub><msubsup><mi>A</mi><mi>m</mi><mrow><mi>P</mi><mo>-</mo><mn>2</mn></mrow></msubsup><msub><mi>B</mi><mi>m</mi></msub><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>C</mi><mi>m</mi></msub><msub><mi>B</mi><mi>m</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000017.GIF" wi="1494" he="87" /></maths><maths num="0007"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><msup><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><msubsup><mi>D</mi><mn>1</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>D</mi><mn>2</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msubsup><mi>D</mi><mi>P</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000018.GIF" wi="730" he="86" /></maths><maths num="0008"><math><![CDATA[<mrow><msub><mi>D</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mi>m</mi></msub><msubsup><mi>A</mi><mi>m</mi><mi>i</mi></msubsup><msub><mi>X</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Y</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>C</mi><mi>m</mi></msub><msub><mi>X</mi><mi>m</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><mi>c</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>+</mo><msup><mi>&alpha;</mi><mi>i</mi></msup><mo>&lsqb;</mo><mi>c</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Y</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>P</mi><mo>,</mo></mrow>]]></math><img file="FDA0000960605360000019.GIF" wi="1580" he="79" /></maths>其中,c(k)=[c<sub>1</sub>(k) c<sub>2</sub>(k) … c<sub>N</sub>(k)]<sup>T</sup>,<img file="FDA0000960605360000021.GIF" wi="1302" he="231" />U(k)是多变量PI预测函数在第k个时刻的控制量输出向量;K<sub>P</sub>、K<sub>i</sub>分别为广义比例项、积分项系数矩阵;f<sub>j</sub>(i)为基函数在t=(k+i)T<sub>s</sub>时的取值,T<sub>s</sub>为采样周期,F<sub>n</sub>为基函数的值构成的向量,J为基函数的阶数,j为基函数个数的索引,从1到J之间整数;Q和R分别表示误差加权矩阵和控制加权矩阵;q<sup>‑1</sup>和q<sup>‑2</sup>为延时算子;Y<sub>P</sub>(k)为k时刻感应电机输出的定子电压直轴分量和定子电压交轴分量组成的向量;c(k)为k时刻感应电机定子电流直轴分量和定子电流交轴分量的参考值组成的向量;X<sub>m</sub>(k)为感应电机的模型状态向量;T<sub>r</sub>是参考轨迹的期望响应时间;α<sup>i</sup>为由<img file="FDA0000960605360000022.GIF" wi="54" he="62" />组成的第i时刻参考轨迹衰减因子矩阵;P为预测步长;i为第i步预测时刻;N为输出变量的个数;步骤3、将步骤2计算所得的控制量向量U(k)作为多变量PI预测函数控制的输入,经过反Park变换及空间矢量规则,产生空间矢量信号,将该信号施加到感应电机的功率驱动电路上,控制电压源型逆变器开关的开通与关断,使电压源型逆变器输出不同幅值和频率的三相交流电,控制感应电机的转速。
地址 215101 江苏省苏州市吴中区木渎镇中山东路70号吴中科技创业园2号楼2310室