CN102857160A - 一种基于多线拟合的变励磁同步电机mtpa控制方法 - Google Patents
一种基于多线拟合的变励磁同步电机mtpa控制方法 Download PDFInfo
- Publication number
- CN102857160A CN102857160A CN2012103506177A CN201210350617A CN102857160A CN 102857160 A CN102857160 A CN 102857160A CN 2012103506177 A CN2012103506177 A CN 2012103506177A CN 201210350617 A CN201210350617 A CN 201210350617A CN 102857160 A CN102857160 A CN 102857160A
- Authority
- CN
- China
- Prior art keywords
- psi
- current
- excitation
- excitation magnetic
- rectangular axis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Control Of Ac Motors In General (AREA)
Abstract
本发明涉及一种基于多线拟合的变励磁同步电机MTPA控制方法,针对控制结构图中的直交轴电流解算模块,提出一种多线拟合的方法,将励磁磁链在其变化范围内进行离散化处理,然后在每一个特定的励磁磁链下,以电磁转矩为单变量进行直交轴电流解算函数的二次多项式拟合,可以得到不同励磁磁链下的多条拟合曲线。系统运行时,针对某时刻的励磁磁链,选择与其最接近的两条拟合曲线先分别进行直交轴电流的计算,然后以这两组值为基础,以此刻的励磁磁链大小为节点进行线性插值,便可得到该时刻直交轴电流值的给定值,最终完成变励磁同步电机的最大转矩电流比控制。
Description
技术领域
本发明一种变励磁同步电机最大转矩电流比控制方法,具体针对励磁电流变化的同步电机利用多线拟合方法进行最大转矩电流比(MTPA)控制的方法,属于交流电机传动技术领域。
背景技术
起动/发电一体化是未来航空交流电源系统的一个重要发展方向。目前我国飞机交流电源系统大都采用三级式无刷同步电机(原理框图见图1)作为发电机,该类发电机无起动航空发动机的功能,发动机由独立的起动机进行起动。这样的发动机-电源系统包含两套电机,使得其体积和重量较大,且系统复杂,可靠性降低。若能在原有三级式无刷同步发电机的基础上,通过控制使其运行在电动状态来完成发动机的起动,即实现起动/发电一体化,就可以省去专门的起动机,减轻机载重量和系统体积。但是,该电机主要为满足发电功能设计,在电动状态时存在如下问题:1)静止和低速状态时,励磁机输出电压较低、主发电机励磁电流较小,严重影响主发电机的带载起动能力;2)随着电机转速的升高,励磁机输出电压逐渐增大,主发电机励磁电流也逐渐增大,即主发电机励磁电流在电机起动过程中一直处于变化之中。
从电机结构上看,主发电机为电励磁凸极同步电机。为了提高主发电机在静止和低速时的带载能力,应选择最大转矩电流比控制策略。对于永磁同步电机或者励磁电流恒定的电励磁同步电机,其励磁磁链或励磁电流是保持不变的,即电磁转矩方程Te=np(Ψdiq-Ψqid)=np[Ψfiq+(Ld-Lq)idiq]中的ψf为恒值,所以针对此类电机的MTPA控制就是要实时的完成如下所示的非线性规划问题,以求得直交轴电流的参考值:
通过对上式的求解,可以获得满足MTPA控制的id、iq与Te的关系:
式(2)中函数f1和f2的精确解析式都很难确定,所以一般采用二次或三次多项式进行拟合,利用拟合多项式搭建直交轴电流解算模块,完成电机的MTPA控制。
当同步电机励磁电流随电机转速发生变化时,即励磁磁链非恒定时,由于式(1)中的ψf为变量,故由式(1)求解出的id、iq的表达式将变为:
对于式(3),函数f1和f2的精确解析式更加难于确定,且不能采用简单的二次或三次多项式进行拟合。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于多线拟合的变励磁同步电机MTPA控制方法,与传统恒励磁同步电机MTPA控制相比,变励磁同步电机MTPA控制的最大区别在于直交轴电流的MTPA解算模块。
本发明的思想在于:在变励磁同步电机MTPA控制中,直交轴电流的MTPA解算模块需要进行修改,最终的控制结构图如图2所示,其中被虚线包围的部分为传统恒励磁下同步电机的MTPA控制结构图。
技术方案
一种基于多线拟合的变励磁同步电机MTPA控制方法,其特征在于步骤如下:
步骤1:对同步电机的励磁磁链Ψf进行等距离离散化处理,得到n个离散点以Ψi表示,i=1,2…,n;
步骤2:求解每一个Ψi对应的直交轴电流非线性规划问题,其中:id为直轴电流,iq为交轴电流,Te为电磁转矩,np为电机极对数,Ld为直轴电感,Lq为交轴电感;
用二次多项式拟合方法获得近似表达式为 其中:idi表示第i个离散点处的直轴电流,iqi表示第i个离散点处的交轴电流,fid和fiq均为关于Te的二次多项式;
步骤3:当系统运行时,对某一时刻的励磁磁链Ψf,选择与之最为接近的两个励磁磁链离散点Ψi和Ψi+1,采用对应的直交轴电流函数表达式 和
步骤4:利用公式 以此刻的励磁磁链Ψf为节点进行线性插值,得到控制变励磁同步电机的d轴电流id和q轴电流iq。
有益效果
本发明提出的一种基于多线拟合的变励磁同步电机MTPA控制方法,针对控制结构图中的直交轴电流解算模块,提出一种多线拟合的方法,即针对式(3),将励磁磁链在其变化范围内进行离散化处理,然后在每一个特定的励磁磁链下,以电磁转矩为单变量进行直交轴电流解算函数的二次多项式拟合,可以得到不同励磁磁链下的多条拟合曲线。系统运行时,针对某时刻的励磁磁链,选择与其最接近的两条拟合曲线先分别进行直交轴电流的计算,然后以这两组值为基础,以此刻的励磁磁链大小为节点进行线性插值,便可得到该时刻直交轴电流值的给定值,最终完成变励磁同步电机的最大转矩电流比控制。
本发明方法具有以下有益效果:
1)采用多线拟合方法进行直交轴电流函数的拟合,解决了变励磁同步电机MTPA控制中直交轴电流函数难以获得的问题;
2)本发明方法中采用的多线拟合方法简单有效,相比于二维查表法,具有拟合程度好,所占内存小的优势。
附图说明
图1:现有技术中的三级式无刷同步起动/发电系统结构图;
图2:本发明提出的变励磁同步电机最大转矩电流比控制结构图;
图3:励磁电流随电机转速变化曲线;
图4:不同励磁磁链下直轴电流的曲线拟合
图5:不同励磁磁链下交轴电流的曲线拟合
图6:变励磁同步电机最大转矩电流比控制仿真结果——转速与转矩曲线
图7:变励磁同步电机最大转矩电流比控制仿真结果——A相电流曲线
具体实施方式
现结合实施例、附图对本发明作进一步描述:
为验证本发明方法,采用Matlab2008b-Simulink6.0进行仿真验证。仿真中电励磁同步电机参数为:极对数nP=3;定子绕组电阻Rs=10.3mΩ;直轴电感Ld=0.63mH;交轴电感Lq=0.31mH;互感Lm=6mH;给定转速nN=3000(r/min)。励磁电流随转速变化的曲线如图3所示。具体仿真条件设定为:起动负载为15Nm,之后逐渐升高至45Nm,当电机转矩达到2200(r/min)时,负载又逐渐下降,至负载为6Nm时稳定。
实施例包含的具体步骤如下:
步骤1:同步电机的励磁磁链Ψf的变化范围为0.06~0.144Wb,对励磁磁链Ψf进行等距离离散化处理,得到离散点为0.06、0.072、0.084、0.096、0.108、0.12、0.132、0.144;
步骤2求解每一个Ψi对应的直交轴电流非线性规划问题:,
针对每一个励磁磁链离散点,代入非线性规划问题 中进行求解;其中:id为直轴电流,iq为交轴电流,Te为电磁转矩,np为电机极对数,Ld为直轴电感,Lq为交轴电感;
用二次多项式拟合方法获得近似表达式为 其中:idi表示第i个离散点处的直轴电流,iqi表示第i个离散点处的交轴电流,fid和fiq均为关于Te的二次多项式;结果如下表所示,在直角坐标系中的表示如图4、5所示。
励磁磁链 | 直轴电流id拟合多项式 | 交轴电流iq拟合多项式 |
0.06 | id1=-0.00592×T2+2.6522×T-11.6664 | iq1=-0.01851×T2+4.27854×T+11.2600 |
0.072 | id2=-0.00161×T2+2.05348×T-11.404 | iq2=-0.01634×T2+4.04198×T+6.32205 |
0.084 | id3=-0.00024×T2+1.68193×T-11.758 | iq3=-0.01140×T2+3.56940×T+5.56863 |
0.096 | id4=0.00177×T2+1.2792×T-10.0375 | iq4=-0.00975×T2+3.33185×T+2.82112 |
0.108 | id5=0.00316×T2+0.94928×T-8.0756 | iq5=-0.00813×T2+3.08662×T+1.10329 |
0.12 | id6=0.00401×T2+0.6905×T-6.22845 | iq6=-0.00665×T2+2.84686×T+0.10980 |
0.132 | id7=0.0044×T2+0.49465×T-4.65968 | iq7=-0.00535×T2+2.62143×T-0.40666 |
0.144 | id8=0.0045×T2+0.35059×T-3.41219 | iq8=-0.00425×T2+2.41525×T-0.62943 |
步骤3:当系统运行时,对某一时刻的励磁磁链Ψf,选择与之最为接近的两个励磁磁链离散点Ψi和Ψi+1,采用对应的直交轴电流函数表达式 和
以步骤2中获得的不同励磁磁链下直交轴电流拟合函数为基础,编写直交轴电流MTPA解算模块。该模块的输入为励磁磁链(记为ph_f)和转矩给定值(记为T),输出为直轴电流给定值和交轴电流给定值。这里以直轴电流解算为例给出具体的实现方法,交轴电流与之类似:
if(ph_f<=0.06)
id=id1;
elseif(ph_f<=0.072)
id=id1+(id2-idl)*(ph_f-0.06)/0.012;
elseif(ph_f<=0.084)
id=id2+(id3-id2)*(ph_f-0.072)/0.012;
elseif(ph_f<=0.096)
id=id3+(id4-id3)*(ph_f-0.084)/0.012;
elseif(ph_f<=0.108)
id=id4+(id5-id4)*(ph_f-0.096)/0.012;
elseif(ph_f<=0.12)
id=id5+(id6-id5)*(ph_f-0.108)/0.012;
elseif(ph_f<=0.132)
id=id6+(id7-id6)*(ph_f-0.12)/0.012;
elseif(ph_f<=0.144)
id=id7+(id8-id7)*(ph_f-0.132)/0.012;
else
id=id8;
end
步骤4:利用公式 以此刻的励磁磁链Ψf为节点进行线性插值,得到控制变励磁同步电机的d轴电流id和q轴电流iq。
通过调节PI参数,最终获得仿真结果如图6、7所示。
以图2为基础,在Matlab中建立变励磁同步电机MTPA控制仿真模型,其中直交轴电流MTPA解算模块按照步骤(3)内容编写,其他部分均与传统恒励磁同步电机MTPA控制中的一致。
Claims (1)
1.一种基于多线拟合的变励磁同步电机MTPA控制方法,其特征在于步骤如下:
步骤1:对同步电机的励磁磁链Ψf进行等距离离散化处理,得到n个离散点以Ψi表示,i=1,2…,n;
步骤2:求解每一个Ψi对应的直交轴电流非线性规划问题,其中:id为直轴电流,iq为交轴电流,Te为电磁转矩,np为电机极对数,Ld为直轴电感,Lq为交轴电感;
用二次多项式拟合方法获得近似表达式为 其中:idi表示第i个离散点处的直轴电流,iqi表示第i个离散点处的交轴电流,fid和fiq均为关于Te的二次多项式;
步骤4:利用公式 以此刻的励磁磁链Ψf为节点进行线性插值,得到控制变励磁同步电机的d轴电流id和q轴电流iq。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012103506177A CN102857160A (zh) | 2012-09-20 | 2012-09-20 | 一种基于多线拟合的变励磁同步电机mtpa控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012103506177A CN102857160A (zh) | 2012-09-20 | 2012-09-20 | 一种基于多线拟合的变励磁同步电机mtpa控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102857160A true CN102857160A (zh) | 2013-01-02 |
Family
ID=47403419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012103506177A Pending CN102857160A (zh) | 2012-09-20 | 2012-09-20 | 一种基于多线拟合的变励磁同步电机mtpa控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102857160A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107896080A (zh) * | 2017-12-14 | 2018-04-10 | 成都雅骏新能源汽车科技股份有限公司 | 内嵌式永磁同步电机mtpa曲线拟合方法及控制系统 |
CN108063575A (zh) * | 2017-11-16 | 2018-05-22 | 中冶南方(武汉)自动化有限公司 | 自启动永磁同步电机i/f启动方法及系统 |
CN109286346A (zh) * | 2017-07-20 | 2019-01-29 | 南通杰夫电气有限公司 | 基于一元一次函数拟合的永磁同步电机mtpa控制方法 |
CN109768751A (zh) * | 2018-12-24 | 2019-05-17 | 北方工业大学 | 一种用于永磁同步电机的最大转矩电流比控制方法 |
CN111115127A (zh) * | 2018-10-30 | 2020-05-08 | 豪夫迈·罗氏有限公司 | 估计驱动系统的操作状态的方法和驱动系统 |
CN111948537A (zh) * | 2020-08-11 | 2020-11-17 | 臻驱科技(上海)有限公司 | 凸极式永磁同步电机最大转矩磁链比工作点的标定方法 |
CN112003508A (zh) * | 2020-09-18 | 2020-11-27 | 蔚然(南京)动力科技有限公司 | 电机无位置传感器控制方法、装置 |
CN112740537A (zh) * | 2020-12-24 | 2021-04-30 | 深圳市英威腾电气股份有限公司 | 永磁同步电机的mtpa控制方法、装置、系统及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020180398A1 (en) * | 2001-06-01 | 2002-12-05 | Hitachi, Ltd. | Driving equipment for synchronous motors |
US6936991B2 (en) * | 2002-06-03 | 2005-08-30 | Ballard Power Systems Corporation | Method and apparatus for motor control |
CN102223133A (zh) * | 2011-06-02 | 2011-10-19 | 西北工业大学 | 一种凸极式永磁同步电机的最大转矩控制方法 |
-
2012
- 2012-09-20 CN CN2012103506177A patent/CN102857160A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020180398A1 (en) * | 2001-06-01 | 2002-12-05 | Hitachi, Ltd. | Driving equipment for synchronous motors |
US6936991B2 (en) * | 2002-06-03 | 2005-08-30 | Ballard Power Systems Corporation | Method and apparatus for motor control |
CN102223133A (zh) * | 2011-06-02 | 2011-10-19 | 西北工业大学 | 一种凸极式永磁同步电机的最大转矩控制方法 |
Non-Patent Citations (4)
Title |
---|
XING SHAOBANG ET AL: "Efficiency Improvement on PMSM Twelve Sectors DTC System", 《CONTROL AND DECISION CONFERENCE (CCDC), 2011 CHINESE》 * |
刘国林等: "基于MTPA的永磁同步电动机矢量控制系统", 《变频器世界》 * |
李军等: "基于分段曲线拟合的IPMSM最大转矩电流比控制研究", 《四川大学学报》 * |
李庆扬等: "《数值分析》", 31 July 2006, 武汉市:华中科技大学出版社 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109286346A (zh) * | 2017-07-20 | 2019-01-29 | 南通杰夫电气有限公司 | 基于一元一次函数拟合的永磁同步电机mtpa控制方法 |
CN108063575A (zh) * | 2017-11-16 | 2018-05-22 | 中冶南方(武汉)自动化有限公司 | 自启动永磁同步电机i/f启动方法及系统 |
CN108063575B (zh) * | 2017-11-16 | 2020-06-09 | 中冶南方(武汉)自动化有限公司 | 自启动永磁同步电机i/f启动方法及系统 |
CN107896080A (zh) * | 2017-12-14 | 2018-04-10 | 成都雅骏新能源汽车科技股份有限公司 | 内嵌式永磁同步电机mtpa曲线拟合方法及控制系统 |
CN107896080B (zh) * | 2017-12-14 | 2019-11-19 | 成都雅骏新能源汽车科技股份有限公司 | 内嵌式永磁同步电机mtpa曲线拟合方法及控制系统 |
CN111115127A (zh) * | 2018-10-30 | 2020-05-08 | 豪夫迈·罗氏有限公司 | 估计驱动系统的操作状态的方法和驱动系统 |
CN109768751A (zh) * | 2018-12-24 | 2019-05-17 | 北方工业大学 | 一种用于永磁同步电机的最大转矩电流比控制方法 |
CN111948537A (zh) * | 2020-08-11 | 2020-11-17 | 臻驱科技(上海)有限公司 | 凸极式永磁同步电机最大转矩磁链比工作点的标定方法 |
CN111948537B (zh) * | 2020-08-11 | 2022-12-02 | 臻驱科技(上海)有限公司 | 凸极式永磁同步电机最大转矩磁链比工作点的标定方法 |
CN112003508A (zh) * | 2020-09-18 | 2020-11-27 | 蔚然(南京)动力科技有限公司 | 电机无位置传感器控制方法、装置 |
CN112740537A (zh) * | 2020-12-24 | 2021-04-30 | 深圳市英威腾电气股份有限公司 | 永磁同步电机的mtpa控制方法、装置、系统及设备 |
CN112740537B (zh) * | 2020-12-24 | 2023-02-17 | 深圳市英威腾电气股份有限公司 | 永磁同步电机的mtpa控制方法、装置、系统及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102857160A (zh) | 一种基于多线拟合的变励磁同步电机mtpa控制方法 | |
CN102710208B (zh) | 基于温度变化的永磁同步电机自适应控制方法 | |
CN101969289B (zh) | 改善电机在场削弱区域工作期间的扭矩线性的增益调节 | |
CN102223133B (zh) | 一种凸极式永磁同步电机的最大转矩控制方法 | |
CN103378788B (zh) | 变频空调用压缩机的驱动方法和装置 | |
CN102857159A (zh) | 基于二元二次函数拟合的变励磁同步电机mtpa控制方法 | |
CN101957884B (zh) | 混合励磁同步电机发电系统建模方法 | |
CN103762922A (zh) | 一种交流伺服弱磁调速的方法 | |
CN110212710A (zh) | 一种车用永磁同步电机设计方法 | |
CN202535303U (zh) | 变频空调用压缩机的驱动装置 | |
CN103595324B (zh) | 一种混合励磁电机弱磁控制方法 | |
CN103427442B (zh) | 可变速抽水蓄能交流励磁发电电动机的同步并网控制方法 | |
CN103401506A (zh) | 一种电动车用隐极式混合励磁电机直接转矩控制方法 | |
CN102904518B (zh) | 一种同步发电机q轴参数在线辨识方法 | |
Bonisławski et al. | Unconventional control system of hybrid excited synchronous machine | |
CN105071735A (zh) | 基于t-1简化模型的异步电机节能控制方法 | |
CN104135203B (zh) | 基于带隙变步长法的异步电机弱磁控制方法 | |
Golovko et al. | MATHEMATICAL SIMULATION OF AUTONOMOUS WIND ELECTRIC INSTALLATION WITH MAGNETOELECTRIC GENERATOR. | |
CN103546083B (zh) | 用于调节用来控制异步电机的操作的电流和/或转矩命令的方法、系统和装置 | |
Li et al. | Analytical derivation of efficiency map of an induction machine for electric vehicle applications | |
CN105978434B (zh) | 一种永磁同步电机最大转矩电流比控制方法及控制系统 | |
CN103904975B (zh) | 一种凸极式同步电动机的控制方法 | |
CN105353220A (zh) | 三相异步电动机电学性能参数的辨识方法 | |
CN104022716A (zh) | 系数拟合的变励磁同步电机最大转矩电流比控制方法 | |
CN103684173A (zh) | 一种双馈风力发电机的变流器控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130102 |