CN108631672A - 计及最优占空比调制的永磁同步电机预测磁链控制方法 - Google Patents
计及最优占空比调制的永磁同步电机预测磁链控制方法 Download PDFInfo
- Publication number
- CN108631672A CN108631672A CN201810424591.3A CN201810424591A CN108631672A CN 108631672 A CN108631672 A CN 108631672A CN 201810424591 A CN201810424591 A CN 201810424591A CN 108631672 A CN108631672 A CN 108631672A
- Authority
- CN
- China
- Prior art keywords
- vector
- flux linkage
- duty ratio
- moment
- formula
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/141—Flux estimation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种计及最优占空比调制的永磁同步电机预测磁链控制方法，通过将磁链矢量的有效值最小化原理引入到基本电压矢量的占空比计算中，并结合电压矢量筛选方法，可在降低处理器计算时间的同时有效抑制转矩和磁链脉动，提高系统的可靠性和稳态性能。首先将转矩增量转化为负载角增量，将控制目标统一为定子磁链矢量；然后根据无差拍思想进行基本电压矢量筛选，在减小备选矢量的前提下可以快速计算出有效电压矢量的作用时间；进而通过磁链矢量的有效值最小化原理计算出电压矢量的占空比，并根据磁链矢量的综合误差确定所选电压矢量和矢量占空比，有效抑制转矩和磁链脉动。
Description
技术领域
本发明涉及一种电机预测磁链控制方法，特别是一种计及最优占空比调制的永磁同步电机预测磁链控制方法。
背景技术
有限状态集模型预测控制是一种基于系统动态模型的计算，可根据不同设计要求设计价值函数，由于概念简单，动态响应快，易于考虑系统的约束以及具有优秀的多变量控制能力等优点，近些年吸引了国内外大量学者对其在电力传动领域的应用近些了研究。传统的模型预测控制目标函数包含转矩和磁链不同量纲的控制变量，需要繁复的权重系数调试；专注于转矩增量控制而对定子磁链的增量作近似处理，导致较大磁链控制偏差；每个控制周期施加一个电压矢量，转矩和磁链脉动较大，稳定性能相对较差。因此，研究一种能够实现算法复杂度减小、转矩脉动降低、磁链控制精度提高和权重系数消除的预测控制算法有着广阔的发展前景。
发明内容
本发明所要解决的技术问题是提供一种计及最优占空比调制的永磁同步电机预测磁链控制方法，用于电机驱动系统的转矩和磁链脉动抑制。
为解决上述技术问题，本发明所采用的技术方案是：
一种计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于包含以下步骤：
步骤一：在每一个控制周期中，利用转速PI控制器实时计算转矩参考值Te *；
步骤二：通过PARK坐标变换器实时计算逆变器输出d/q轴电流id/iq，进而将获得的d/q轴电流输入转矩估计模块实时计算转矩Te，然后利用转矩PI控制器实时计算负载角增量Δδsf；
步骤三：基于最大转矩电流比原理，通过磁链幅值参考值估计模块计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)；
步骤四：将获得的d/q轴电流输入磁链估计模块实时获得磁链矢量ψsd/ψsq，进而将获得的磁链矢量输入负载角估计模块实时计算负载角δsf，然后结合负载角增量Δδsf获得(k+1)时刻的负载角参考值
步骤五：将获得的(k+1)时刻磁链矢量参考幅值ψs *(k+1)和负载角参考值δsf(k+1)输入定子磁链矢量参考值计算模块得到(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)；
步骤六：将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)；
步骤七：将获得的参考电压矢量输入电压矢量筛选模块，进而根据扇区选取两个基本非零矢量ui(i＝1,2)；
步骤八：将两个基本非零矢量ui(i＝1,2)、磁链矢量ψs和零矢量u0输入占空比计算模块得到非零矢量对应的占空比di(i＝1,2)；
步骤九：将(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、基本电压矢量ui(i＝1,2)、占空比di(i＝1,2)和磁链矢量ψs输入最小化目标函数模块得到所选电压矢量uopt和占空比dopt，进而将开关状态送给逆变器得到三相电驱动永磁同步电机。
进一步地，所述步骤一具体为
将参考转速n*与实际转速n的差值en输入转速PI控制器，根据公式(1.1)获得参考转矩Te *：
其中，KPn和KIn分别为转速PI控制器的比例增益和积分增益，s为拉普拉斯积分项。
进一步地，所述步骤二中电机负载角增量Δδsf获取方法为
2.1将采集得到的逆变器输出三相电流和电机转子磁链位置角输入PARK坐标变换器，根据公式(2.1)进行坐标变换获得d/q轴电流id/iq；
其中，ia、ib和ic为逆变器输出三相电流，θr为转子磁链位置角；
2.2将获得的d/q轴电流id/iq输入转矩估计模块，根据公式(2.2)实时计算转矩Te；
其中，P为极对数，ψf为永磁体磁链幅值，Ld和Lq分别为直轴和交轴电感；
2.3将获得的参考转矩Te *与实时转矩Te的差值eT输入转矩PI控制器，根据公式(2.3)获得负载角增量Δδsf；
其中，KPT和KIT分别为转速PI控制器的比例增益和积分增益。
进一步地，所述步骤三中(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)获取方法为
将获得的参考转矩Te *输入磁链幅值参考值估计模块，根据公式(3.1)计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)：
进一步地，所述步骤四中计算(k+1)时刻负载角δsf(k+1)的方法具体为
4.1将获得的d/q轴电流输入磁链估计模块，根据公式(4.1)计算磁链矢量ψsd/ψsq；
4.2将获得的磁链矢量输入负载角估计模块，根据公式(4.2)实时计算电机负载角δsf，然后根据公式(4.3)将实时负载角增量与实时负载角相加得到(k+1)时刻的负载角
进一步地，所述步骤五中计算(k+1)时刻参考磁链矢量ψsd *(k+1)/ψsq *(k+1)方法为
将获得的定子磁链矢量参考幅值ψs *(k+1)和负载角δsf(k+1)输入定子磁链矢量参考值计算模块，根据公式(5.1)计算(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)：
进一步地，所述步骤六中计算(k+1)时刻参考电压矢量ud *(k+1)/uq *(k+1)的方法为
将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块，根据公式(6.1)计算得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)：
其中，Rs为定子电阻，Ts为采样时间，ωe为转子电角速度。
进一步地，所述步骤七中选取两个基本非零矢量ui(i＝1,2)的方法为
根据公式(7.1)将获得的参考电压矢量经过反PARK变换得到uα *(k+1)/uβ *(k+1)；进而根据公式(7.2)获得参考电压矢量在αβ轴的参考角度θs *；然后将θs *按[0,π/3)、[π/3,2π/3)、[2π/3,π)、[π,4π/3)、[4π/3,5π/3)、[5π/3,2π)划分为6个扇区，分别编号I、II、III、IV、V、VI，选取得到基本非零矢量ui(i＝1,2)：
进一步地，所述步骤八中计算两个基本非零矢量对应的占空比di(i＝1,2)的方法为
将选取得到基本非零矢量ui(i＝1,2)、零矢量u0、磁链矢量ψsd/ψsq输入占空比计算模块，根据公式(8.1)计算非零矢量和零矢量分别作用时磁链沿d/q轴的变化率，进而计算公式(8.2)的有效值最小值，计算得到的占空比di(i＝1,2)如公式(8.3)所示：
其中，S11/S12分别为非零矢量/零矢量沿d轴的磁链变化率，S21/S22分别为非零矢量/零矢量沿q轴的磁链变化率，uid(i＝1,2)为选取的非零矢量沿d轴的分量，uiq(i＝1,2)为选取的非零矢量沿d轴的分量；
其中，S1为非零矢量和零矢量沿d轴的磁链变化率之差，S1＝S11-S12，S2为非零矢量和零矢量沿q轴的磁链变化率之差，S2＝S21-S22，Er为磁链矢量有效值与磁链矢量给定值之差，ψsd e为定子磁链沿d轴的初始误差，ψsd e＝ψsd-ψsd *(k+1)，ψsq e为定子磁链沿q轴的初始误差，ψsq e＝ψsq-ψsq *(k+1)；
其中，C＝S1Ts(2ψsd e+S12Ts)+S2Ts(2ψsq e+S22Ts)，A＝S1(-2S1-S12)+S2(-2S2-S22)。
进一步地，所述步骤九中将(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、基本电压矢量ui(i＝1,2)、占空比di(i＝1,2)和磁链矢量ψs输入最小化目标函数模块，选取使得公式(8.2)值取最小时的电压矢量uopt和占空比dopt，进而将开关状态送给逆变器得到三相电驱动永磁同步电机。
本发明与现有技术相比，具有以下优点和效果：
1、采用永磁同步电机定子磁链矢量为控制目标，消除了传统模型预测转矩控制中的权值。
2、根据无差拍思想计算理想电压矢量所在扇区，通过判断扇区快速筛选电压矢量，减小算法的计算次数，使得数字处理器的运算负担降低。
3、电压矢量和占空比同时选择以最小化价值函数，提高了稳态性能，降低了开关频率。
4、根据有效值最小化原理计算磁链矢量的综合误差，有效抑制转矩和磁链脉动。
附图说明
图1是本发明的计及最优占空比调制的永磁同步电机预测磁链控制方法流程图。
图2是本发明的计及最优占空比调制的永磁同步电机预测磁链控制方法定子磁链d轴分量变化率。
图3是本发明的计及最优占空比调制的永磁同步电机预测磁链控制方法定子磁链q轴分量变化率。
图4是本发明的计及最优占空比调制的永磁同步电机预测磁链控制方法无差拍控制矢量图。
图5是本发明的计及最优占空比调制的永磁同步电机预测磁链控制方法矢量图。
具体实施方式
下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发明的解释而本发明并不局限于以下实施例。
如图1所示，一种计及最优占空比调制的永磁同步电机预测磁链控制方法，包括如下步骤：
步骤1：在每一个控制周期中，利用转速PI控制器1实时计算转矩参考值Te *。
将参考转速n*与实际转速n的差值en输入转速PI控制器1，根据公式(1.1)获得参考转矩Te *：
其中，KPn和KIn分别为转速PI控制器1的比例增益和积分增益，s为拉普拉斯积分项。
步骤2：通过PARK坐标变换器2实时计算逆变器输出d/q轴电流id/iq，进而将获得的d/q轴电流输入转矩估计模块3实时计算转矩Te，然后利用转矩PI控制器4实时计算负载角增量Δδsf。
电机负载角增量Δδsf获取方法包括如下步骤：
1)将采集得到的逆变器输出三相电流和电机转子磁链位置角输入PARK坐标变换器2，根据公式(2.1)进行坐标变换获得d/q轴电流id/iq；
其中，ia、ib和ic为逆变器输出三相电流，θr为转子磁链位置角。
2)将获得的d/q轴电流id/iq输入转矩估计模块3，根据公式(2.2)实时计算转矩Te；
其中，P为极对数，ψf为永磁体磁链幅值，Ld和Lq分别为直轴和交轴电感。
3)将获得的参考转矩Te *与实时转矩Te的差值eT输入转矩PI控制器4，根据公式(2.3)获得负载角增量Δδsf；
其中，KPT和KIT分别为转速PI控制器1的比例增益和积分增益。
步骤3：基于最大转矩电流比原理，通过磁链幅值参考值估计模块5计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)。
(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)获取方法为：将获得的参考转矩Te *输入磁链幅值参考值估计模块5，根据公式(3.1)计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)：
步骤4：将获得的d/q轴电流输入磁链估计模块6实时获得磁链矢量ψsd/ψsq，进而将获得的磁链矢量输入负载角估计模块7实时计算负载角δsf，然后结合负载角增量Δδsf获得(k+1)时刻的负载角参考值
(k+1)时刻负载角δsf(k+1)的方法包括如下步骤：
1)将获得的d/q轴电流输入磁链估计模块6，根据公式(4.1)计算磁链矢量ψsd/ψsq；
2)将获得的磁链矢量输入负载角估计模块7，根据公式(4.2)实时计算电机负载角δsf，然后根据公式(4.3)将实时负载角增量与实时负载角相加得到(k+1)时刻的负载角
步骤5：将获得的(k+1)时刻磁链矢量参考幅值ψs *(k+1)和负载角参考值δsf(k+1)输入定子磁链矢量参考值计算模块8得到(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)。
计算(k+1)时刻参考磁链矢量ψsd *(k+1)/ψsq *(k+1)方法为：将获得的定子磁链矢量参考幅值ψs *(k+1)和负载角δsf(k+1)输入定子磁链矢量参考值计算模块8，根据公式(5.1)计算(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)：
步骤6：将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块9得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)。
计算(k+1)时刻参考电压矢量ud *(k+1)/uq *(k+1)的方法为：将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块9，根据公式(6.1)计算得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)：
其中，Rs为定子电阻，Ts为采样时间，ωe为转子电角速度。
步骤7：将获得的参考电压矢量输入电压矢量筛选模块10，进而根据扇区选取两个基本非零矢量ui(i＝1,2)。
选取两个基本非零矢量ui(i＝1,2)的方法为：根据公式(7.1)将获得的参考电压矢量经过反PARK变换得到uα *(k+1)/uβ *(k+1)；进而根据公式(7.2)获得参考电压矢量在αβ轴的参考角度θs *；然后将θs *按[0,π/3)、[π/3,2π/3)、[2π/3,π)、[π,4π/3)、[4π/3,5π/3)、[5π/3,2π)划分为6个扇区，分别编号I、II、III、IV、V、VI，表1所示为电压矢量筛选表，根据表1选取得到基本非零矢量ui(i＝1,2)：
表1电压矢量筛选表
步骤8：将两个基本非零矢量ui(i＝1,2)、磁链矢量ψs和零矢量u0输入占空比计算模块11得到非零矢量对应的占空比di(i＝1,2)。
计算两个基本非零矢量对应的占空比di(i＝1,2)的方法为：将选取得到基本非零矢量ui(i＝1,2)、零矢量u0、磁链矢量ψsd/ψsq输入占空比计算模块，根据公式(8.1)计算非零矢量和零矢量分别作用时磁链沿d/q轴的变化率，进而计算公式(8.2)的有效值最小值，计算得到的占空比di(i＝1,2)如公式(8.3)所示：
其中，S11/S12分别为非零矢量/零矢量沿d轴的磁链变化率，S21/S22分别为非零矢量/零矢量沿q轴的磁链变化率，uid(i＝1,2)为选取的非零矢量沿d轴的分量，uiq(i＝1,2)为选取的非零矢量沿d轴的分量。
其中，S1为非零矢量和零矢量沿d轴的磁链变化率之差，S1＝S11-S12，S2为非零矢量和零矢量沿q轴的磁链变化率之差，S2＝S21-S22，Er为磁链矢量有效值与磁链矢量给定值之差，ψsd e为定子磁链沿d轴的初始误差，ψsd e＝ψsd-ψsd *(k+1)，ψsq e为定子磁链沿q轴的初始误差，ψsq e＝ψsq-ψsq *(k+1)。
其中，C＝S1Ts(2ψsd e+S12Ts)+S2Ts(2ψsq e+S22Ts)，A＝S1(-2S1-S12)+S2(-2S2-S22)。
步骤9：将(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、基本电压矢量ui(i＝1,2)、占空比di(i＝1,2)和磁链矢量ψs输入最小化目标函数模块12得到所选电压矢量uopt和占空比dopt，进而将开关状态送给逆变器13得到三相电驱动永磁同步电机14。
工作原理：本发明方法通过将磁链矢量的有效值最小化原理引入到基本电压矢量的占空比计算中，并结合电压矢量筛选方法，可在降低处理器计算时间的同时有效抑制转矩和磁链脉动，提高系统的可靠性和稳态性能。首先将转矩增量转化为负载角增量，将控制目标统一为定子磁链矢量，避免权值系数整定问题；然后根据无差拍思想进行基本电压矢量筛选，在减小备选矢量的前提下可以快速计算出有效电压矢量的作用时间，降低系统采样频率；进而通过磁链矢量的有效值最小化原理计算出电压矢量的占空比，并根据磁链矢量的综合误差确定所选电压矢量和矢量占空比，有效抑制转矩和磁链脉动。
本说明书中所描述的以上内容仅仅是对本发明所作的举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种修改或补充或采用类似的方式替代，只要不偏离本发明说明书的内容或者超越本权利要求书所定义的范围，均应属于本发明的保护范围。
Claims (10)
1.一种计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于包含以下步骤：
步骤一：在每一个控制周期中，利用转速PI控制器实时计算转矩参考值Te *；
步骤二：通过PARK坐标变换器实时计算逆变器输出d/q轴电流id/iq，进而将获得的d/q轴电流输入转矩估计模块实时计算转矩Te，然后利用转矩PI控制器实时计算负载角增量Δδsf；
步骤三：基于最大转矩电流比原理，通过磁链幅值参考值估计模块计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)；
步骤四：将获得的d/q轴电流输入磁链估计模块实时获得磁链矢量ψsd/ψsq，进而将获得的磁链矢量输入负载角估计模块实时计算负载角δsf，然后结合负载角增量Δδsf获得(k+1)时刻的负载角参考值δsf *(k+1)；
步骤五：将获得的(k+1)时刻磁链矢量参考幅值ψs *(k+1)和负载角参考值δsf(k+1)输入定子磁链矢量参考值计算模块得到(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)；
步骤六：将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)；
步骤七：将获得的参考电压矢量输入电压矢量筛选模块，进而根据扇区选取两个基本非零矢量ui(i＝1,2)；
步骤八：将两个基本非零矢量ui(i＝1,2)、磁链矢量ψs和零矢量u0输入占空比计算模块得到非零矢量对应的占空比di(i＝1,2)；
步骤九：将(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、基本电压矢量ui(i＝1,2)、占空比di(i＝1,2)和磁链矢量ψs输入最小化目标函数模块得到所选电压矢量uopt和占空比dopt，进而将开关状态送给逆变器得到三相电驱动永磁同步电机。
2.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤一具体为
将参考转速n*与实际转速n的差值en输入转速PI控制器，根据公式(1.1)获得参考转矩Te *：
其中，KPn和KIn分别为转速PI控制器的比例增益和积分增益，s为拉普拉斯积分项。
3.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤二中电机负载角增量Δδsf获取方法为
2.1将采集得到的逆变器输出三相电流和电机转子磁链位置角输入PARK坐标变换器，根据公式(2.1)进行坐标变换获得d/q轴电流id/iq；
其中，ia、ib和ic为逆变器输出三相电流，θr为转子磁链位置角；
2.2将获得的d/q轴电流id/iq输入转矩估计模块，根据公式(2.2)实时计算转矩Te；
其中，P为极对数，ψf为永磁体磁链幅值，Ld和Lq分别为直轴和交轴电感；
2.3将获得的参考转矩Te *与实时转矩Te的差值eT输入转矩PI控制器，根据公式(2.3)获得负载角增量Δδsf；
其中，KPT和KIT分别为转速PI控制器的比例增益和积分增益。
4.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤三中(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)获取方法为
将获得的参考转矩Te *输入磁链幅值参考值估计模块，根据公式(3.1)计算(k+1)时刻的定子磁链矢量参考幅值ψs *(k+1)：
5.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤四中计算(k+1)时刻负载角δsf(k+1)的方法具体为
4.1将获得的d/q轴电流输入磁链估计模块，根据公式(4.1)计算磁链矢量ψsd/ψsq；
4.2将获得的磁链矢量输入负载角估计模块，根据公式(4.2)实时计算电机负载角δsf，然后根据公式(4.3)将实时负载角增量与实时负载角相加得到(k+1)时刻的负载角δsf *(k+1)；
δs * f(k+1)＝δsf+Δδsf(4.3)。
6.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤五中计算(k+1)时刻参考磁链矢量ψsd *(k+1)/ψsq *(k+1)方法为
将获得的定子磁链矢量参考幅值ψs *(k+1)和负载角δsf(k+1)输入定子磁链矢量参考值计算模块，根据公式(5.1)计算(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)：
7.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤六中计算(k+1)时刻参考电压矢量ud *(k+1)/uq *(k+1)的方法为
将获得的(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、磁链矢量ψsd/ψsq和d/q轴电流id/iq输入无差拍计算模块，根据公式(6.1)计算得到(k+1)时刻的参考电压矢量ud *(k+1)/uq *(k+1)：
其中，Rs为定子电阻，Ts为采样时间，ωe为转子电角速度。
8.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤七中选取两个基本非零矢量ui(i＝1,2)的方法为
根据公式(7.1)将获得的参考电压矢量经过反PARK变换得到uα *(k+1)/uβ *(k+1)；进而根据公式(7.2)获得参考电压矢量在αβ轴的参考角度θs *；然后将θs *按[0,π/3)、[π/3,2π/3)、[2π/3,π)、[π,4π/3)、[4π/3,5π/3)、[5π/3,2π)划分为6个扇区，分别编号I、II、III、IV、V、VI，选取得到基本非零矢量ui(i＝1,2)：
9.按照权利要求1所述的计及最优占空比调制的永磁同步电机
预测磁链控制方法，其特征在于：所述步骤八中计算两个基本非零矢量对应的占空比di(i＝1,2)的方法为
将选取得到基本非零矢量ui(i＝1,2)、零矢量u0、磁链矢量ψsd/ψsq输入占空比计算模块，根据公式(8.1)计算非零矢量和零矢量分别作用时磁链沿d/q轴的变化率，进而计算公式(8.2)的有效值最小值，计算得到的占空比di(i＝1,2)如公式(8.3)所示：
其中，S11/S12分别为非零矢量/零矢量沿d轴的磁链变化率，S21/S22分别为非零矢量/零矢量沿q轴的磁链变化率，uid(i＝1,2)为选取的非零矢量沿d轴的分量，uiq(i＝1,2)为选取的非零矢量沿d轴的分量；
其中，S1为非零矢量和零矢量沿d轴的磁链变化率之差，S1＝S11-S12，S2为非零矢量和零矢量沿q轴的磁链变化率之差，S2＝S21-S22，Er为磁链矢量有效值与磁链矢量给定值之差，ψsd e为定子磁链沿d轴的初始误差，ψsd e＝ψsd-ψsd *(k+1)，ψsq e为定子磁链沿q轴的初始误差，ψsq e＝ψsq-ψsq *(k+1)；
其中，C＝S1Ts(2ψsd e+S12Ts)+S2Ts(2ψsq e+S22Ts)，A＝S1(-2S1-S12)+S2(-2S2-S22)。
10.按照权利要求1所述的计及最优占空比调制的永磁同步电机预测磁链控制方法，其特征在于：所述步骤九中将(k+1)时刻的参考磁链矢量ψsd *(k+1)/ψsq *(k+1)、基本电压矢量ui(i＝1,2)、占空比di(i＝1,2)和磁链矢量ψs输入最小化目标函数模块，选取使得公式(8.2)值取最小时的电压矢量uopt和占空比dopt，进而将开关状态送给逆变器得到三相电驱动永磁同步电机。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810424591.3A CN108631672B (zh) | 2018-05-07 | 2018-05-07 | 计及最优占空比调制的永磁同步电机预测磁链控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810424591.3A CN108631672B (zh) | 2018-05-07 | 2018-05-07 | 计及最优占空比调制的永磁同步电机预测磁链控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108631672A true CN108631672A (zh) | 2018-10-09 |
| CN108631672B CN108631672B (zh) | 2020-05-12 |
Family
ID=63695541
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810424591.3A Active CN108631672B (zh) | 2018-05-07 | 2018-05-07 | 计及最优占空比调制的永磁同步电机预测磁链控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108631672B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109600085A (zh) * | 2018-12-05 | 2019-04-09 | 天津工业大学 | 基于可变控制集的永磁同步电机直接预测占空比控制方法 |
| CN111431460A (zh) * | 2020-04-25 | 2020-07-17 | 南通大学 | 一种永磁同步电机无传感器模型预测磁链控制方法 |
| CN111726046A (zh) * | 2020-07-28 | 2020-09-29 | 南通大学 | 一种计及占空比优化的非对称六相pmsm模型预测磁链控制方法 |
| CN111800056A (zh) * | 2020-07-21 | 2020-10-20 | 中国石油大学(华东) | 一种基于新型开关表的永磁同步电机三矢量模型预测转矩控制方法 |
| CN113067515A (zh) * | 2021-04-13 | 2021-07-02 | 南通大学 | 一种计及占空比约束的永磁同步电机三矢量模型预测磁链控制方法 |
| CN113315385A (zh) * | 2021-05-20 | 2021-08-27 | 山东大学 | 海洋直驱永磁同步风电变流器动态级联预测控制方法 |
| CN113489398A (zh) * | 2021-06-03 | 2021-10-08 | 西安理工大学 | 一种内置式永磁同步电机无位置传感器参数误差补偿策略 |
| CN114079412A (zh) * | 2021-11-19 | 2022-02-22 | 天津大学 | 一种基于相电压占空比计算的电机预测控制方法 |
| CN115549540A (zh) * | 2022-10-28 | 2022-12-30 | 南京航空航天大学 | 基于坐标变换抑制转矩脉动的电励磁双凸极电机控制方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105763120A (zh) * | 2016-03-29 | 2016-07-13 | 东南大学 | 一种永磁同步电机准无差拍模型预测磁链控制方法 |
| CN106936356A (zh) * | 2017-04-24 | 2017-07-07 | 东南大学盐城新能源汽车研究院 | 矢量筛选和占空比结合的电机模型预测控制系统及方法 |
| CN107017810A (zh) * | 2017-04-24 | 2017-08-04 | 东南大学盐城新能源汽车研究院 | 永磁同步电机无权值模型预测转矩控制系统及方法 |
-
2018
- 2018-05-07 CN CN201810424591.3A patent/CN108631672B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105763120A (zh) * | 2016-03-29 | 2016-07-13 | 东南大学 | 一种永磁同步电机准无差拍模型预测磁链控制方法 |
| CN106936356A (zh) * | 2017-04-24 | 2017-07-07 | 东南大学盐城新能源汽车研究院 | 矢量筛选和占空比结合的电机模型预测控制系统及方法 |
| CN107017810A (zh) * | 2017-04-24 | 2017-08-04 | 东南大学盐城新能源汽车研究院 | 永磁同步电机无权值模型预测转矩控制系统及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 牛峰等: "基于占空比调制的永磁同步电机直接转矩控制", 《电工技术学报》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109600085B (zh) * | 2018-12-05 | 2020-08-14 | 天津工业大学 | 基于可变控制集的永磁同步电机直接预测占空比控制方法 |
| CN109600085A (zh) * | 2018-12-05 | 2019-04-09 | 天津工业大学 | 基于可变控制集的永磁同步电机直接预测占空比控制方法 |
| CN111431460A (zh) * | 2020-04-25 | 2020-07-17 | 南通大学 | 一种永磁同步电机无传感器模型预测磁链控制方法 |
| CN111431460B (zh) * | 2020-04-25 | 2023-10-31 | 南通大学 | 一种永磁同步电机无传感器模型预测磁链控制方法 |
| CN111800056A (zh) * | 2020-07-21 | 2020-10-20 | 中国石油大学(华东) | 一种基于新型开关表的永磁同步电机三矢量模型预测转矩控制方法 |
| CN111726046B (zh) * | 2020-07-28 | 2022-03-04 | 南通大学 | 一种计及占空比优化的非对称六相pmsm模型预测磁链控制方法 |
| CN111726046A (zh) * | 2020-07-28 | 2020-09-29 | 南通大学 | 一种计及占空比优化的非对称六相pmsm模型预测磁链控制方法 |
| CN113067515A (zh) * | 2021-04-13 | 2021-07-02 | 南通大学 | 一种计及占空比约束的永磁同步电机三矢量模型预测磁链控制方法 |
| CN113067515B (zh) * | 2021-04-13 | 2024-05-17 | 南通大学 | 一种计及占空比约束的永磁同步电机三矢量模型预测磁链控制方法 |
| CN113315385B (zh) * | 2021-05-20 | 2023-02-28 | 山东大学 | 海洋直驱永磁同步风电变流器动态级联预测控制方法 |
| CN113315385A (zh) * | 2021-05-20 | 2021-08-27 | 山东大学 | 海洋直驱永磁同步风电变流器动态级联预测控制方法 |
| CN113489398B (zh) * | 2021-06-03 | 2022-04-01 | 西安理工大学 | 一种内置式永磁同步电机无位置传感器参数误差补偿策略 |
| CN113489398A (zh) * | 2021-06-03 | 2021-10-08 | 西安理工大学 | 一种内置式永磁同步电机无位置传感器参数误差补偿策略 |
| CN114079412A (zh) * | 2021-11-19 | 2022-02-22 | 天津大学 | 一种基于相电压占空比计算的电机预测控制方法 |
| CN114079412B (zh) * | 2021-11-19 | 2023-04-18 | 天津大学 | 一种基于相电压占空比计算的电机预测控制方法 |
| CN115549540A (zh) * | 2022-10-28 | 2022-12-30 | 南京航空航天大学 | 基于坐标变换抑制转矩脉动的电励磁双凸极电机控制方法 |
| CN115549540B (zh) * | 2022-10-28 | 2023-05-05 | 南京航空航天大学 | 基于坐标变换抑制转矩脉动的电励磁双凸极电机控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108631672B (zh) | 2020-05-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108631672A (zh) | 计及最优占空比调制的永磁同步电机预测磁链控制方法 | |
| CN107017810A (zh) | 永磁同步电机无权值模型预测转矩控制系统及方法 | |
| CN103401503B (zh) | 一种在谐波平面在线辨识双三相电机参数的方法及装置 | |
| CN108736778A (zh) | 一种永磁同步电机双矢量预测磁链控制方法 | |
| CN207166388U (zh) | 矢量筛选和占空比结合的电机模型预测控制系统 | |
| CN106533300B (zh) | 一种基于速度环模糊控制和高频注入法的无传感器控制系统 | |
| CN109274303A (zh) | 一种基于扩展控制集的异步电机模型预测控制方法 | |
| CN106788054B (zh) | 一种基于旋转高频注入法的无速度传感器控制方法 | |
| CN110336501A (zh) | 一种内嵌式永磁同步电机模型预测控制方法 | |
| CN109039189A (zh) | 基于几何法的永磁同步电机两矢量预测控制系统及方法 | |
| CN107612446A (zh) | 一种内置式永磁同步电机模型预测转矩控制方法 | |
| CN109525158A (zh) | 空调压缩机无差拍电流预测控制方法和系统 | |
| CN105680752B (zh) | 一种永磁同步电机饱和电感参数的辨识方法及系统 | |
| CN109617490A (zh) | 一种五相永磁电机直接转矩控制下的最大转矩电流比控制方法 | |
| CN105871281A (zh) | 一种改进的永磁同步电机模型预测控制算法 | |
| CN106849809A (zh) | 一种基于磁链误差矢量法的svm‑dtc电机控制算法 | |
| US11522480B2 (en) | SPMSM sensorless composite control method with dual sliding-mode observers | |
| CN109039204A (zh) | 车用永磁同步电机无位置传感器模型预测控制系统及方法 | |
| CN112803861B (zh) | 一种永磁同步电机三矢量模型预测控制的无零矢量算法 | |
| CN113037170B (zh) | 电机控制方法、装置及终端设备 | |
| CN110112974A (zh) | 电机控制方法、控制器、存储介质及电机驱动系统 | |
| CN105591575B (zh) | 一种隐极式永磁同步电机直接特征控制系统及控制方法 | |
| CN110530083A (zh) | 一种压缩机电机控制方法、装置及空调器 | |
| CN206640528U (zh) | 永磁同步电机无权值模型预测转矩控制系统 | |
| CN109245640A (zh) | 一种异步电机模型预测控制方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210630 Address after: 226019 No.205, building 6, Nantong University, No.9, Siyuan Road, Nantong City, Jiangsu Province Patentee after: Center for technology transfer, Nantong University Address before: 226000 9 Siyuan Road, Chongchuan District, Nantong City, Jiangsu Province Patentee before: NANTONG University |
|
| TR01 | Transfer of patent right | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20181009 Assignee: NANTONG SHENGBANG ROLL MAKING CO.,LTD. Assignor: Center for technology transfer, Nantong University Contract record no.: X2021980016487 Denomination of invention: Predictive flux linkage control method of permanent magnet synchronous motor considering optimal duty cycle modulation Granted publication date: 20200512 License type: Common License Record date: 20211228 |
|
| EE01 | Entry into force of recordation of patent licensing contract |