发明名称 一种永磁同步电机准无差拍模型预测磁链控制方法
摘要 本发明公开了一种永磁同步电机准无差拍模型预测磁链控制方法,在有限控制框架内,首先计算下一时刻的定子磁链矢量参考值,再根据参考值按照无差拍控制思想,获得目标电压矢量;通过目标电压矢量的位置角判断目标电压矢量所在扇区,利用扇区选择3个有效电压矢量;然后,根据3个有效电压矢量预测下一时刻的定子磁链矢量;最后,通过优化价值函数获得逆变器最优开关状态,逆变器根据最优开关状态输出电压给永磁同步电机。本方法的价值函数中不含有权值计算,在进行价值函数优化时,只需针对3个有效电压矢量进行优化,减少了算法运算时间。
申请公布号 CN105763120A 申请公布日期 2016.07.13
申请号 CN201610188463.4 申请日期 2016.03.29
申请人 东南大学 发明人 花为;黄文涛;程明
分类号 H02P21/14(2016.01)I;H02P21/12(2016.01)I 主分类号 H02P21/14(2016.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 黄成萍
主权项 一种永磁同步电机准无差拍模型预测磁链控制方法,其特征在于:包括如下步骤:(1)计算参考转矩T<sub>e</sub><sup>*</sup>(k):将参考速度ω<sup>*</sup>(k)与实际反馈速度ω(k)的差值e(k)输入PI控制器(1),根据公式(1.1)计算参考转矩T<sub>e</sub><sup>*</sup>(k);<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><msub><mi>T</mi><mi>e</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>K</mi><mi>P</mi></msub><mo>+</mo><mfrac><msub><mi>K</mi><mi>I</mi></msub><mi>s</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1.1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952427950000011.GIF" wi="646" he="135" /></maths>其中:K<sub>P</sub>和K<sub>I</sub>分别为PI控制器(1)的比例增益和积分增益;(2)计算(k+1)时刻定子磁链矢量参考ψ<sub>s</sub><sup>*</sup>(k+1):将定子磁链幅值参考|ψ<sub>s</sub><sup>*</sup>(k)|和参考转矩T<sub>e</sub><sup>*</sup>(k)输入转矩角计算模块(2‑1),根据公式(2.1)计算转矩角参考δ<sup>*</sup>(k);将转子位置角θ<sub>r</sub>(k)和实际反馈速度ω(k)输入转子位置角预测模块(2‑2),根据公式(2.2)预测(k+1)时刻转子位置角θ<sub>r</sub>(k+1);将转矩角参考δ<sup>*</sup>(k)和(k+1)时刻转子位置角θ<sub>r</sub>(k+1)相加,得到(k+1)时刻定子磁链位置角θ<sub>s</sub>(k+1);将定子磁链幅值参考|ψ<sub>s</sub><sup>*</sup>(k)|和(k+1)时刻定子磁链位置角θ<sub>s</sub>(k+1)输入定子磁链矢量参考计算模块(2‑3),根据公式(2.4)计算(k+1)时刻定子磁链矢量参考ψ<sub>s</sub><sup>*</sup>(k+1);<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>&delta;</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>arcsin</mi><mfrac><mrow><mn>2</mn><msub><mi>L</mi><mi>s</mi></msub><msup><msub><mi>T</mi><mi>e</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mn>3</mn><msub><mi>P</mi><mi>r</mi></msub><mo>|</mo><msup><msub><mi>&psi;</mi><mi>s</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><msup><msub><mi>&psi;</mi><mi>f</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2.1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952427950000012.GIF" wi="882" he="151" /></maths>θ<sub>r</sub>(k+1)=θ<sub>r</sub>(k)+ω(k)T<sub>s</sub>  (2.2)θ<sub>s</sub>(k+1)=δ<sup>*</sup>(k)+θ<sub>r</sub>(k+1)  (2.3)ψ<sub>s</sub><sup>*</sup>(k+1)=|ψ<sub>s</sub><sup>*</sup>(k)|exp(jθ<sub>s</sub>(k+1))  (2.4)其中:L<sub>s</sub>为永磁同步电机同步电感,P<sub>r</sub>为永磁同步电机极对数,|ψ<sub>f</sub><sup>*</sup>(k)|为永磁体磁链幅值,T<sub>s</sub>为预测控制的采样时间;(3)有效电压矢量选择:将(k+1)时刻定子磁链矢量参考ψ<sub>s</sub><sup>*</sup>(k+1)输入目标矢量计算模块(3‑1),根据公式(3.1)计算目标电压矢量u<sub>obj</sub>(k+1);将目标电压矢量u<sub>obj</sub>(k+1)输入目标电压矢量角计算模块(3‑2),根据公式(3.2)计算目标电压矢量与α轴夹角θ<sub>u</sub>(k+1);将目标电压矢量与α轴夹角θ<sub>u</sub>(k+1)输入扇区判断模块(3‑3),根据目标电压矢量与α轴夹角θ<sub>u</sub>(k+1)判断目标电压矢量u<sub>obj</sub>(k+1)的扇区号N(k+1);将扇区号N(k+1)输入选择模块(3‑4),获得三个基本电压矢量u<sub>i</sub>(k+1),矢量编号<img file="FDA0000952427950000023.GIF" wi="454" he="54" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>u</mi><mrow><mi>o</mi><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><msub><mi>&psi;</mi><mi>s</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&psi;</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><msub><mi>T</mi><mi>s</mi></msub></mfrac><mo>+</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3.1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952427950000021.GIF" wi="974" he="139" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>u</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>arc</mi><mi> </mi><mi>t</mi><mi>a</mi><mi>n</mi><mfrac><mrow><msub><mi>u</mi><mrow><mi>o</mi><mi>b</mi><mi>j</mi><mi>&beta;</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>u</mi><mrow><mi>o</mi><mi>b</mi><mi>j</mi><mi>&alpha;</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3.2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952427950000022.GIF" wi="742" he="151" /></maths>其中:ψ<sub>s</sub>(k)为k时刻定子磁链矢量,R<sub>s</sub>为定子电阻,i<sub>s</sub>(k)为定子电流,u<sub>objβ</sub>(k+1)和u<sub>objα</sub>(k+1)分别为u<sub>obj</sub>(k+1)的α轴和β轴分量;(4)计算(k+1)时刻定子磁链矢量预测值ψ<sub>s</sub>(k+1):将定子电流i<sub>s</sub>(k)和三个基本电压矢量u<sub>i</sub>(k+1)输入定子磁链矢量预测模块(4),根据公式(4.1)计算(k+1)时刻定子磁链矢量预测值ψ<sub>s</sub>(k+1);ψ<sub>s</sub>(k+1)=ψ<sub>s</sub>(k)+T<sub>s</sub>(u<sub>i</sub>(k+1)‑R<sub>s</sub>i<sub>s</sub>(k))  (4.1)(5)选择逆变器最优开关状态:将(k+1)时刻定子磁链矢量参考ψ<sub>s</sub><sup>*</sup>(k+1)和(k+1)时刻定子磁链矢量预测值ψ<sub>s</sub>(k+1)输入优化模块(5),根据公式(5.1)计算价值函数g<sub>i</sub>,价值函数g<sub>i</sub>取最小值时,其对应的基本电压矢量确定为最优基本电压矢量u<sub>opt</sub>,根据最优基本电压矢量u<sub>opt</sub>获得最优开关状态S<sub>a,b,c</sub>;g<sub>i</sub>=|ψ<sub>s</sub><sup>*</sup>(k+1)‑ψ<sub>s</sub>(k+1)|  (5.1)(6)输出最优电压:逆变器将最优开关状态S<sub>a,b,c</sub>转化为最优电压输送给永磁同步电机。
地址 211189 江苏省南京市江宁区东南大学路2号