CN107196543A - 共直流母线开绕组异步电机系统零序环流抑制方法 - Google Patents

共直流母线开绕组异步电机系统零序环流抑制方法 Download PDF

Info

Publication number
CN107196543A
CN107196543A CN201710541205.4A CN201710541205A CN107196543A CN 107196543 A CN107196543 A CN 107196543A CN 201710541205 A CN201710541205 A CN 201710541205A CN 107196543 A CN107196543 A CN 107196543A
Authority
CN
China
Prior art keywords
msub
mrow
mtd
mtr
inverter
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
Application number
CN201710541205.4A
Other languages
English (en)
Other versions
CN107196543B (zh
Inventor
杨淑英
姚乐
张兴
李浩源
谢震
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201710541205.4A priority Critical patent/CN107196543B/zh
Publication of CN107196543A publication Critical patent/CN107196543A/zh
Application granted granted Critical
Publication of CN107196543B publication Critical patent/CN107196543B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/493—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode the static converters being arranged for operation in parallel
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M1/00—Details of apparatus for conversion
    • H02M1/12—Arrangements for reducing harmonics from ac input or output
    • 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
    • H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
    • H02P27/08—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
    • H02P27/085—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation wherein the PWM mode is adapted on the running conditions of the motor, e.g. the switching frequency
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M1/00—Details of apparatus for conversion
    • H02M1/12—Arrangements for reducing harmonics from ac input or output
    • H02M1/123—Suppression of common mode voltage or current

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Inverter Devices (AREA)
  • External Artificial Organs (AREA)

Abstract

本发明要解决的技术问题为针对现有技术中单纯从调制策略环节抑制零序环流，忽略逆变器非线性、反电动势零序分量，导致共直流母线开绕组异步电机定子电流畸变明显、零序环流明显；基于此本发明提出了一种零序电压可控的空间矢量PWM调制方案，通过对共直流母线开绕组异步电机系统零序电压的主动控制实现对零序环流的抑制，同时避免使用零矢量，减小开绕组电机系统的共模电压。

Description

共直流母线开绕组异步电机系统零序环流抑制方法
技术领域
本发明属于电气工程领域的逆变器控制技术，具体涉及一种共直流母线开绕组异步电机系统零序环流抑制方法。
背景技术
相比于传统的单逆变器驱动异步电机调速系统，共直流母线开绕组异步电机系统具有容错性高、拓扑结构灵活、直流侧电压等级低等优点。基于此，共直流母线开绕组异步电机系统越来越多地受到学界的关注。
虽然共直流母线开绕组异步电机具有良好的应用前景，但其作为一种较新的电机拓扑结构，还面临着诸多有待解决的问题，其中之一便是系统的零序环流问题。系统的零序环流由系统的零序电压引起，零序电压的来源主要有，逆变器调制策略本身产生的零序电压；逆变器非线性因素导致的零序电压；开绕组异步电机发电动势中的零序分量等。
目前，针对开绕组异步电机零序环流问题有学术论文对此做了深入的理论分析，如2014年IEEE文献“Comparative Evaluation of SVPWM Strategies for a DualInverter fed Open-End Winding Induction Motor Drive with a Single DC PowerSupply”(“共直流母线双逆变器供电开绕组异步电机驱动的SVPWM(Space Vector PulseWidth Modulation)策略比较”——2014年第四十届IEEE工业电子学会年会)所述，其中涉及到以下几种调制方案：
1)120°解耦SVPWM策略；
2)180°解耦SVPWM策略；
3)SAZE(The Sample Averaged Zero-sequence Elimination)SVPWM策略；
4)双电源供电策略；
2016年IEEE文献“Angular Modulation of Dual-Inverter Fed Open-End Motorfor Electrical Vehicle Applications”(“电动汽车用双逆变器供电开绕组电机角度调制”——2016年4月IEEE电力电子协会期刊，第31卷，第4期)，该文章指出调制策略中使用零矢量会使系统存在较大的共模电压。
通过上述文章分析以及相关结果可以发现，方案1)、2)、3)中开绕组异步电机定子电流存在明显的畸变，零序环流分量明显。这4种方案存在着明显的不足，如下：
(1)方案1)、2)、3)未考虑到逆变器非线性导致的零序电压分量；
(2)方案1)、2)、3)未考虑到开绕组异步电机反电动势中的零序分量产生的零序电压；
(3)方案4)采用两个独立直流源供电，切断了零序环流的回路，但是增加了系统的硬件体积与成本。
(4)4种方案中均使用零矢量，导致开绕组电机系统存在较大的共模电压。
发明内容
本发明要解决的技术问题为针对现有技术中单纯从调制策略环节抑制零序环流，忽略逆变器非线性、反电动势零序分量，导致共直流母线开绕组异步电机定子电流畸变明显、零序环流明显；基于此本发明提出了一种零序电压可控的空间矢量PWM调制方案，通过对共直流母线开绕组异步电机系统零序电压的主动控制实现对零序环流的抑制，同时避免使用零矢量，减小开绕组电机系统的共模电压。
为解决本发明的技术问题，所采用的技术方案如下。
本发明提供了一种共直流母线开绕组异步电机系统零序环流抑制方法，本抑制方法涉及的共直流母线开绕组异步电机拓扑结构具体描述如下：直流母线电压为Udc；直流侧设置有一个稳压电容C，电容C的正极连接第一逆变器和第二逆变器正极，电容C负极连接第一逆变器和第二逆变器负极；第一逆变器第一相桥臂包括两只开关管S11、S12以及两只二极管D11、D12，S11与S12串联，S11发射极连接S12集电极，D11正极连接S11发射极，负极连接S11集电极，D12正极连接S12发射极，负极连接S12集电极；第一逆变器第二相桥臂包括两只开关管S13、S14以及两只二极管D13、D14，S13与S14串联，S13发射极连接S14集电极；D13正极连接S13发射极，负极连接S13集电极，D14正极连接S14发射极，负极连接S14集电极；第一逆变器第三相桥臂包括两只开关管S15、S16以及两只二极管D15、D16，S15与S16串联，S15发射极连接S16集电极，D15正极连接S15发射极，负极连接S15集电极，D16正极连接S16发射极，负极连接S16集电极；第二逆变器第一相桥臂包括两只开关管S21、S22以及两只二极管D21、D22，S21与S22串联，S21发射极连接S22集电极，D21正极连接S21发射极，负极连接S21集电极，D22正极连接S22发射极，负极连接S22集电极；第二逆变器第二相桥臂包括两只开关管S23、S24以及两只二极管D23、D24，S23与S24串联，S23发射极连接S24集电极，D23正极连接S23发射极，负极连接S23集电极，D24正极连接S24发射极，负极连接S24集电极；第二逆变器第三相桥臂包括两只开关管S25、S26以及两只二极管D25、D26，S25与S26串联，S25发射极连接S26集电极，D25正极连接S25发射极，负极连接S25集电极，D26正极连接S26发射极，负极连接S26集电极；开绕组异步电机三相定子绕组分别为A、B、C，A相两端连接第一逆变器第一相桥臂中点a1和第二逆变器第一相桥臂中点a2，B相两端连接第一逆变器第二相桥臂中点b1和第二逆变器第二相桥臂中点b2，C相两端连接第一逆变器第三相桥臂中点c1和第二逆变器第三相桥臂中点c2；
本抑制方法包括以下步骤：
步骤1，采集开绕组异步电机定子上三相电流iA、iB、iC和转子电角速度ωr；
步骤2，将步骤1采集得到的开绕组异步电机定子上三相电流iA、iB、iC按照下式(1)、(2)变换到同步旋转dq坐标系下：
其中iα、iβ为两相静止坐标系αβ中开绕组异步电机定子α轴电流和β轴电流，i0为开绕组异步电机系统的零序环流，id、iq为旋转坐标系dq坐标系中开绕组异步电机定子d轴电流和q轴电流、θ为转子定向角度；
步骤3，根据步骤2获得的id、iq、i0计算调节器输入值：
其中分别为设定的开绕组异步电机定子d轴电流，q轴电流以及零序环流，idin为定子d轴电流调节器输入值，iqin为定子q轴电流调节器输入值，i0in为零序环流调节器输入值；
步骤4，设设定参考电压矢量为OU，计算设定参考电压矢量OU在αβ坐标系下的分量uα、uβ：
其中udout为定子d轴电流调节器输出值，uqout定子q轴电流调节器输出值，θ为转子定向角度，uα、uβ为设定参考电压矢量OU在αβ坐标系下的分量；
步骤5，将设定参考电压矢量OU在αβ坐标系下的分量uα、uβ按照下式变换到abc坐标系下，得到ua、ub、uc：
其中OU为设定参考电压矢量，ua、ub、uc为设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量；
步骤6，按照下式对设定参考电压矢量OU进行180°解耦：
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，U0为零序环流调节器输出值，u0为第一逆变器零轴电压目标值，u0'为第二逆变器零轴电压目标值；
步骤7，根据下表对电压矢量OU1、OU2所在αβ坐标系中扇区位置进行判断：
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，ua、ub、uc分别设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量，扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系
中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域；
步骤8，根据步骤7中获得的第一逆变器和第二逆变器的合成参考电压矢量OU1、OU2所在αβ坐标系中扇区位置信息，按照下表选择第一逆变器和第二逆变器参与合成的电压空间矢量；
其中扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域；U1、U2、U3、U4、U5、U6为第一逆变器产生的电压空间矢量，U1’、U2’、U3’、U4’、U5’、U6’为第二逆变器产生的电压空间矢量；
第一逆变器产生的电压空间矢量对应的开关状态如下表所示：
其中sa1、sb1、sc1分别表示第一逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示上开关管导通，下开关管关断；等于0时表示上开关管关断，下开关管导通；
第二逆变器产生的电压空间矢量对应的开关状态如下表：
其中sa2、sb2、sc2分别第二逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示桥臂上开关管导通，下开关管关断；等于0时表示桥臂上开关管关断，下开关管导通；
步骤9，计算第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量的作用时间；
第一逆变器根据步骤8的选择结果，在OU1所在扇区内选择电压空间矢量U1、U2、U3、U4，根据伏秒平衡原理建立下式：
求解得到t1、t2、t3、t4；
第二逆变器根据步骤8的选择结果，在OU2所在扇区内选择电压空间矢量U1'、U2'、U3'、U4'，根据伏秒平衡原理建立下式：
求解得到t1'、t2'、t3'、t4'；
其中，Ts为一个开关周期的时间，u0为第一逆变器零轴电压目标值，U1、U2、U3、U4为步骤8中选择的第一逆变器产生的电压空间矢量，U10、U20、U30、U40分别为U1、U2、U3、U4在零轴上的分量，U1α、U2α、U3α、U4α分别为U1、U2、U3、U4在α轴上的分量，U1β、U2β、U3β、U4β分别为U1、U2、U3、U4在β轴上的分量，t1、t2、t3、t4分别为U1、U2、U3、U4的作用时间，U1'、U2'、U3'、U4'为步骤4中选择的第二逆变器产生的电压空间矢量，u0'为第二逆变器零轴电压目标值，U10'、U20'、U30'、U40'分别为U1'、U2'、U3'、U4'在零轴上的分量，U'1α、U'2α、U'3α、U'4α分别为U1'、U2'、U3'、U4'在α轴上的分量，U'1β、U'2β、U'3β、U'4β分别为U1'、U2'、U3'、U4'在β轴上的分量。t1'、t2'、t3'、t4'分别为U1'、U2'、U3'、U4'的作用时间；
步骤10，根据步骤9计算得到的第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量U1、U2、U3、U4和U1'、U2'、U3'、U4'的作用时间画出第一逆变器和第二逆变器的开关波形；
步骤11，根据步骤10得到的开关波形得到第一逆变器和第二逆变器中开关管的翻转时刻；
步骤12，根据步骤11获取的第一逆变器和第二逆变器中开关管的翻转时刻更新比较寄存器的值，从而由数字信号处理器产生两个逆变器的开关信号，完成两个逆变器的电压矢量调制信号的生成。
优选地，其特征在于，步骤1中开绕组异步电机三相定子电流iA、iB、iC是利用电流传感器实时采样得到。
优选地，其特征在于，步骤1中转子电角速度的ωr采集步骤如下：
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N；
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr，其计算公式为：
在公式(9)中，M为光电式旋转编码器旋转一周所产生的脉冲数，P为开绕组异步电机极对数，T为预设采样周期。
优选地，其特征在于，步骤2中转子定向角度θ的获得方式为：
θ＝∫ωrdt (10)
同时限制θ的范围为0°-360°；其中，ωr为转子电角速度，dt为时间的微分。
优选地，其特征在于，步骤3中d轴电流调节器、q轴电流调节器、零序环流调节器均为比例积分调节器。
附图说明
图1是共直流母线开绕组异步电机拓扑结构。
图2是本发明的控制框图。
图3是第一逆变器输出的电压空间矢量图。
图4是第二逆变器输出的电压空间矢量图。
图5是第一逆变器扇区I的矢量合成原理图
图6是第一逆变器在扇区Ⅰ的开关波形图。
图7是u0设置为定值0时共直流母线开绕组异步电机定子单相电流变化示意图。
图8是u0设置为定值0时共直流母线开绕组异步电机零序环流变化示意图。
图9是本发明调制方式下共直流母线开绕组异步电机定子单相电流变化示意图。
图10是本发明调制方式下共直流母线开绕组异步电机零序环流变化示意图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整的描述。显然所描述的实施例仅是本发明实施例的一部分，基于本发明的实施例，本领域的技术人员在没有做出创造性劳动的前提下获得的其它实施例，都属于本专利的保护范围。
本发明提出了一种零序电压可控的空间矢量PWM调制方案，通过对共直流母线开绕组异步电机系统零序电压的主动控制实现对零序环流的抑制，同时避免使用零矢量，减小开绕组电机系统的共模电压。
图1是本发明涉及的共直流母线开绕组异步电机拓扑结构。描述如下：直流母线电压为Udc；直流侧设置有一个稳压电容C，电容C的正极连接第一逆变器和第二逆变器正极，电容C负极连接第一逆变器和第二逆变器负极；第一逆变器第一相桥臂包括两只开关管S11、S12以及两只二极管D11、D12，S11与S12串联，S11发射极连接S12集电极，D11正极连接S11发射极，负极连接S11集电极，D12正极连接S12发射极，负极连接S12集电极；第一逆变器第二相桥臂包括两只开关管S13、S14以及两只二极管D13、D14，S13与S14串联，S13发射极连接S14集电极；D13正极连接S13发射极，负极连接S13集电极，D14正极连接S14发射极，负极连接S14集电极；第一逆变器第三相桥臂包括两只开关管S15、S16以及两只二极管D15、D16，S15与S16串联，S15发射极连接S16集电极，D15正极连接S15发射极，负极连接S15集电极，D16正极连接S16发射极，负极连接S16集电极；第二逆变器第一相桥臂包括两只开关管S21、S22以及两只二极管D21、D22，S21与S22串联，S21发射极连接S22集电极，D21正极连接S21发射极，负极连接S21集电极，D22正极连接S22发射极，负极连接S22集电极；第二逆变器第二相桥臂包括两只开关管S23、S24以及两只二极管D23、D24，S23与S24串联，S23发射极连接S24集电极，D23正极连接S23发射极，负极连接S23集电极，D24正极连接S24发射极，负极连接S24集电极；第二逆变器第三相桥臂包括两只开关管S25、S26以及两只二极管D25、D26，S25与S26串联，S25发射极连接S26集电极，D25正极连接S25发射极，负极连接S25集电极，D26正极连接S26发射极，负极连接S26集电极；开绕组异步电机三相定子绕组分别为A、B、C，A相两端连接第一逆变器第一相桥臂中点a1和第二逆变器第一相桥臂中点a2，B相两端连接第一逆变器第二相桥臂中点b1和第二逆变器第二相桥臂中点b2，C相两端连接第一逆变器第三相桥臂中点c1和第二逆变器第三相桥臂中点c2。
图2是本发明的控制框图。根据图2，本发明包括以下步骤：
步骤1，采集开绕组异步电机定子上三相电流iA、iB、iC和转子电角速度ωr。
所述三相定子电流iA、iB、iC是利用电流传感器实时采样得到；
所述转子电角速度的ωr采集步骤如下：
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N；
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr，其计算公式为：
在公式(9)中，M为光电式旋转编码器旋转一周所产生的脉冲数，P为开绕组异步电机极对数，T为预设采样周期。
步骤2，将步骤1采集得到的开绕组异步电机定子上三相电流iA、iB、iC按照下式(1)、(2)变换到同步旋转dq坐标系下：
其中iα、iβ为两相静止坐标系αβ中开绕组异步电机定子α轴电流和β轴电流，i0为开绕组异步电机系统的零序环流，id、iq为旋转坐标系dq坐标系中开绕组异步电机定子d轴电流和q轴电流、θ为转子定向角度。
转子定向角度θ的获得方式为：
θ＝∫ωrdt (10)
同时限制θ的范围为0°-360°；其中，ωr为转子电角速度，dt为时间的微分。
步骤3，根据步骤2获得的id、iq、i0计算调节器输入值：
其中分别为设定的开绕组异步电机定子d轴电流，q轴电流以及零序环流，idin为定子d轴电流调节器输入值，iqin为定子q轴电流调节器输入值，i0in为零序环流调节器输入值。在实施例中，所述d轴电流调节器、q轴电流调节器、零序环流调节器均为比例积分调节器。
步骤4，设设定参考电压矢量为OU，计算设定参考电压矢量OU在αβ坐标系下的分量uα、uβ：
其中udout为定子d轴电流调节器输出值，uqout定子q轴电流调节器输出值，θ为转子定向角度，uα、uβ为设定参考电压矢量OU在αβ坐标系下的分量。
步骤5，将设定参考电压矢量OU在αβ坐标系下的分量uα、uβ按照下式变换到abc坐标系下，得到ua、ub、uc：
其中OU为设定参考电压矢量，ua、ub、uc为设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量。
步骤6，按照下式对设定参考电压矢量OU进行180°解耦：
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，U0为零序环流调节器输出值，u0为第一逆变器零轴电压目标值，u0'为第二逆变器零轴电压目标值。
步骤7，根据下表对电压矢量OU1、OU2所在αβ坐标系中扇区位置进行判断：
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，ua、ub、uc分别设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量，扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域，例如，当ua≥0，ub<0，uc<0时，可判定参考电压矢量OU1所在αβ坐标系中扇区位置为扇区Ⅰ，参考电压矢量OU2所在αβ坐标系中扇区位置为扇区Ⅳ。
步骤8，根据步骤7中获得的第一逆变器和第二逆变器的合成参考电压矢量OU1、OU2所在αβ坐标系中扇区位置信息，按照下表选择第一逆变器和第二逆变器参与合成的电压空间矢量。
其中扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域；U1、U2、U3、U4、U5、U6为第一逆变器产生的电压空间矢量，U1’、U2’、U3’、U4’、U5’、U6’为第二逆变器产生的电压空间矢量，例如，第一逆变器的合成参考电压矢量OU1所在αβ坐标系中扇区位置为扇区Ⅰ，其选择参与合成的电压空间矢量为U6、U1、U2、U3。第二逆变器的合成参考电压矢量OU2所在αβ坐标系中扇区位置为扇区Ⅰ，其选择参与合成的电压空间矢量为U6’、U1’、U2’、U3’。
第一逆变器产生的电压空间矢量对应的开关状态如下表所示：
其中sa1、sb1、sc1分别表示第一逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示上开关管导通，下开关管关断；等于0时表示上开关管关断，下开关管导通，例如，当电压矢量为U1时，开关状态sa1sb1sc1为100，即此时sa1为1，sb1为0，sc1也为0。
第二逆变器产生的电压空间矢量对应的开关状态如下表：
其中sa2、sb2、sc2分别第二逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示桥臂上开关管导通，下开关管关断；等于0时表示桥臂上开关管关断，下开关管导通，例如，当电压矢量为U1’时，开关状态sa2sb2sc2为100，即此时sa2为1，sb2为0，sc2也为0。
图3是第一逆变器输出的电压空间矢量图。图4是第二逆变器输出的电压空间矢量图。图5是第一逆变器扇区I的矢量合成原理图。
步骤9，计算第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量的作用时间；
第一逆变器根据步骤8的选择结果，在OU1所在扇区内选择电压空间矢量U1、U2、U3、U4，根据伏秒平衡原理建立下式：
求解得到t1、t2、t3、t4。
第二逆变器根据步骤8的选择结果，在OU2所在扇区内选择电压空间矢量U1'、U2'、U3'、U4'，根据伏秒平衡原理建立下式：
求解得到t1'、t2'、t3'、t4'。
其中，Ts为一个开关周期的时间，u0为第一逆变器零轴电压目标值，U1、U2、U3、U4为步骤8中选择的第一逆变器产生的电压空间矢量，U10、U20、U30、U40分别为U1、U2、U3、U4在零轴上的分量，U1α、U2α、U3α、U4α分别为U1、U2、U3、U4在α轴上的分量，U1β、U2β、U3β、U4β分别为U1、U2、U3、U4在β轴上的分量，t1、t2、t3、t4分别为U1、U2、U3、U4的作用时间，U1'、U2'、U3'、U4'为步骤4中选择的第二逆变器产生的电压空间矢量，u0'为第二逆变器零轴电压目标值，U10'、U20'、U30'、U40'分别为U1'、U2'、U3'、U4'在零轴上的分量，U'1α、U'2α、U'3α、U'4α分别为U1'、U2'、U3'、U4'在α轴上的分量，U'1β、U'2β、U'3β、U'4β分别为U1'、U2'、U3'、U4'在β轴上的分量。t1'、t2'、t3'、t4'分别为U1'、U2'、U3'、U4'的作用时间。
步骤10，根据步骤9计算得到的第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量U1、U2、U3、U4和U1'、U2'、U3'、U4'的作用时间画出第一逆变器和第二逆变器的开关波形。
步骤11，根据步骤10得到的开关波形得到第一逆变器和第二逆变器中开关管的翻转时刻。
步骤12，根据步骤11获取的第一逆变器和第二逆变器中开关管的翻转时刻更新比较寄存器的值，从而由数字信号处理器产生两个逆变器的开关信号，完成两个逆变器的电压矢量调制信号的生成。
图6是第一逆变器在扇区Ⅰ的开关波形图。由图6可见，第一逆变器在第一扇区选择矢量U6、U1、U2、U3，其作用时间分别为T6、T1、T2、T3,从图中可以看出开关管的反转时刻，由此设置数字处理器比较寄存器的值分别为从而由数字信号处理器产生两个逆变器的开关信号，完成两个逆变器的空间矢量调制信号的生成。
图7是u0设置为定值0时共直流母线开绕组异步电机定子单相电流变化示意图。图8是u0设置为定值0时共直流母线开绕组异步电机零序环流变化示意图。图9是本发明调制方式下共直流母线开绕组异步电机定子单相电流变化示意图。图10是本发明调制方式下共直流母线开绕组异步电机零序环流变化示意图。由此，可以看出本发明对于共直流母线开绕组异步电机系统零序环流起到了很好的抑制作用。

Claims (5)

1.一种共直流母线开绕组异步电机系统零序环流抑制方法，本抑制方法涉及的共直流母线开绕组异步电机拓扑结构具体描述如下：直流母线电压为Udc；直流侧设置有一个稳压电容C，电容C的正极连接第一逆变器和第二逆变器正极，电容C负极连接第一逆变器和第二逆变器负极；第一逆变器第一相桥臂包括两只开关管S11、S12以及两只二极管D11、D12，S11与S12串联，S11发射极连接S12集电极，D11正极连接S11发射极，负极连接S11集电极，D12正极连接S12发射极，负极连接S12集电极；第一逆变器第二相桥臂包括两只开关管S13、S14以及两只二极管D13、D14，S13与S14串联，S13发射极连接S14集电极；D13正极连接S13发射极，负极连接S13集电极，D14正极连接S14发射极，负极连接S14集电极；第一逆变器第三相桥臂包括两只开关管S15、S16以及两只二极管D15、D16，S15与S16串联，S15发射极连接S16集电极，D15正极连接S15发射极，负极连接S15集电极，D16正极连接S16发射极，负极连接S16集电极；第二逆变器第一相桥臂包括两只开关管S21、S22以及两只二极管D21、D22，S21与S22串联，S21发射极连接S22集电极，D21正极连接S21发射极，负极连接S21集电极，D22正极连接S22发射极，负极连接S22集电极；第二逆变器第二相桥臂包括两只开关管S23、S24以及两只二极管D23、D24，S23与S24串联，S23发射极连接S24集电极，D23正极连接S23发射极，负极连接S23集电极，D24正极连接S24发射极，负极连接S24集电极；第二逆变器第三相桥臂包括两只开关管S25、S26以及两只二极管D25、D26，S25与S26串联，S25发射极连接S26集电极，D25正极连接S25发射极，负极连接S25集电极，D26正极连接S26发射极，负极连接S26集电极；开绕组异步电机三相定子绕组分别为A、B、C，A相两端连接第一逆变器第一相桥臂中点a1和第二逆变器第一相桥臂中点a2，B相两端连接第一逆变器第二相桥臂中点b1和第二逆变器第二相桥臂中点b2，C相两端连接第一逆变器第三相桥臂中点c1和第二逆变器第三相桥臂中点c2；
本抑制方法的特征在于，包括以下步骤：
步骤1，采集开绕组异步电机定子上三相电流iA、iB、iC和转子电角速度ωr；
步骤2，将步骤1采集得到的开绕组异步电机定子上三相电流iA、iB、iC按照下式(1)、(2)变换到同步旋转dq坐标系下：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&amp;alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&amp;beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>A</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>B</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>C</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>q</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>&amp;alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>&amp;beta;</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中iα、iβ为两相静止坐标系αβ中开绕组异步电机定子α轴电流和β轴电流，i0为开绕组异步电机系统的零序环流，id、iq为旋转坐标系dq坐标系中开绕组异步电机定子d轴电流和q轴电流、θ为转子定向角度；
步骤3，根据步骤2获得的id、iq、i0计算调节器输入值：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mrow> <mi>d</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>i</mi> <mi>d</mi> <mo>*</mo> </msubsup> <mo>-</mo> <msub> <mi>i</mi> <mi>d</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>i</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>-</mo> <msub> <mi>i</mi> <mi>q</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mrow> <mn>0</mn> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>i</mi> <mn>0</mn> <mo>*</mo> </msubsup> <mo>-</mo> <msub> <mi>i</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> 1
其中分别为设定的开绕组异步电机定子d轴电流，q轴电流以及零序环流，idin为定子d轴电流调节器输入值，iqin为定子q轴电流调节器输入值，i0in为零序环流调节器输入值；
步骤4，设设定参考电压矢量为OU，计算设定参考电压矢量OU在αβ坐标系下的分量uα、uβ：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>d</mi> <mi>o</mi> <mi>u</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>q</mi> <mi>o</mi> <mi>u</mi> <mi>t</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中udout为定子d轴电流调节器输出值，uqout定子q轴电流调节器输出值，θ为转子定向角度，uα、uβ为设定参考电压矢量OU在αβ坐标系下的分量；
步骤5，将设定参考电压矢量OU在αβ坐标系下的分量uα、uβ按照下式变换到abc坐标系下，得到ua、ub、uc：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mi>a</mi> </msub> <mo>=</mo> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mi>b</mi> </msub> <mo>=</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mi>c</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中OU为设定参考电压矢量，ua、ub、uc为设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量；
步骤6，按照下式对设定参考电压矢量OU进行180°解耦：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>OU</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>O</mi> <mi>U</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>OU</mi> <mn>2</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>O</mi> <mi>U</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>U</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>u</mi> <mn>0</mn> <mo>,</mo> </msubsup> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>U</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，U0为零序环流调节器输出值，u0为第一逆变器零轴电压目标值，u0'为第二逆变器零轴电压目标值；
步骤7，根据下表对电压矢量OU1、OU2所在αβ坐标系中扇区位置进行判断：
OU1所在扇区 Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ OU2所在扇区 Ⅳ Ⅴ Ⅵ Ⅰ Ⅱ Ⅲ ua ua≥0 ua≥0 ua<0 ua<0 ua<0 ua≥0 ub ub<0 ub≥0 ub≥0 ub≥0 ub<0 ub<0 uc uc<0 uc<0 uc<0 uc≥0 uc≥0 uc≥0
其中OU1为第一逆变器合成参考电压矢量，OU2为第二逆变器合成参考电压矢量，ua、ub、uc分别设定参考电压矢量OU在abc坐标系中a轴、b轴和c轴上的分量，扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域；
步骤8，根据步骤7中获得的第一逆变器和第二逆变器的合成参考电压矢量OU1、OU2所在αβ坐标系中扇区位置信息，按照下表选择第一逆变器和第二逆变器参与合成的电压空间矢量；
其中扇区Ⅰ表示αβ坐标系中0°-60°区域，α轴所在位置为0°，β轴所在位置为90°，扇区Ⅱ表示αβ坐标系中60°-120°区域，扇区Ⅲ表示αβ坐标系中120°-180°区域，扇区Ⅳ表示αβ坐标系中180°-240°区域，扇区Ⅴ表示αβ坐标系中240°-300°区域，扇区Ⅵ表示αβ坐标系中300°-360°区域；U1、U2、U3、U4、U5、U6为第一逆变器产生的电压空间矢量，U1’、U2’、U3’、U4’、U5’、U6’为第二逆变器产生的电压空间矢量；
第一逆变器产生的电压空间矢量对应的开关状态如下表所示：
其中sa1、sb1、sc1分别表示第一逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示上开关管导通，下开关管关断；等于0时表示上开关管关断，下开关管导通；
第二逆变器产生的电压空间矢量对应的开关状态如下表：
其中sa2、sb2、sc2分别第二逆变器第一桥臂、第二桥臂、第三桥臂的开关状态，等于1时表示桥臂上开关管导通，下开关管关断；等于0时表示桥臂上开关管关断，下开关管导通；
步骤9，计算第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量的作用时间；
第一逆变器根据步骤8的选择结果，在OU1所在扇区内选择电压空间矢量U1、U2、U3、U4，根据伏秒平衡原理建立下式：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <mi>U</mi> <mn>10</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <mi>U</mi> <mn>20</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <mi>U</mi> <mn>30</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <mi>U</mi> <mn>40</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>1</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>2</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>3</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>4</mn> <mi>&amp;alpha;</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>1</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>2</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>3</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <mi>U</mi> <mrow> <mn>4</mn> <mi>&amp;beta;</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> 3
求解得到t1、t2、t3、t4；
第二逆变器根据步骤8的选择结果，在OU2所在扇区内选择电压空间矢量U1'、U2'、U3'、U4'，根据伏秒平衡原理建立下式：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msup> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>,</mo> </msup> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mn>10</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mn>20</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mn>30</mn> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mn>40</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msup> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mo>,</mo> </msup> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>1</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>2</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>3</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>4</mn> <mi>&amp;alpha;</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <msup> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mo>,</mo> </msup> <mo>=</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>1</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>2</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>3</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>3</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mn>4</mn> </msub> <msub> <msup> <mi>U</mi> <mo>,</mo> </msup> <mrow> <mn>4</mn> <mi>&amp;beta;</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>=</mo> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>,</mo> </msup> <mo>+</mo> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>,</mo> </msup> <mo>+</mo> <msup> <msub> <mi>t</mi> <mn>3</mn> </msub> <mo>,</mo> </msup> <mo>+</mo> <msup> <msub> <mi>t</mi> <mn>4</mn> </msub> <mo>,</mo> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
求解得到t1'、t2'、t3'、t4'；
其中，Ts为一个开关周期的时间，u0为第一逆变器零轴电压目标值，U1、U2、U3、U4为步骤8中选择的第一逆变器产生的电压空间矢量，U10、U20、U30、U40分别为U1、U2、U3、U4在零轴上的分量，U1α、U2α、U3α、U4α分别为U1、U2、U3、U4在α轴上的分量，U1β、U2β、U3β、U4β分别为U1、U2、U3、U4在β轴上的分量，t1、t2、t3、t4分别为U1、U2、U3、U4的作用时间，U1'、U2'、U3'、U4'为步骤4中选择的第二逆变器产生的电压空间矢量，u0'为第二逆变器零轴电压目标值，U10'、U20'、U30'、U40'分别为U1'、U2'、U3'、U4'在零轴上的分量，U'1α、U'2α、U'3α、U'4α分别为U1'、U2'、U3'、U4'在α轴上的分量，U'1β、U'2β、U'3β、U'4β分别为U1'、U2'、U3'、U4'在β轴上的分量。t1'、t2'、t3'、t4'分别为U1'、U2'、U3'、U4'的作用时间；
步骤10，根据步骤9计算得到的第一逆变器和第二逆变器每个扇区选择的相应的4个电压空间矢量U1、U2、U3、U4和U1'、U2'、U3'、U4'的作用时间画出第一逆变器和第二逆变器的开关波形；
步骤11，根据步骤10得到的开关波形得到第一逆变器和第二逆变器中开关管的翻转时刻；
步骤12，根据步骤11获取的第一逆变器和第二逆变器中开关管的翻转时刻更新比较寄存器的值，从而由数字信号处理器产生两个逆变器的开关信号，完成两个逆变器的电压矢量调制信号的生成。
2.根据权利要求1所述的一种共直流母线开绕组异步电机系统零序环流抑制方法，其特征是步骤1中开绕组异步电机三相定子电流iA、iB、iC是利用电流传感器实时采样得到。
3.根据权利要求1所述的一种共直流母线开绕组异步电机系统零序环流抑制方法，其特征是步骤1中转子电角速度的ωr采集步骤如下：
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N；
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr，其计算公式为：
<mrow> <msub> <mi>T&amp;omega;</mi> <mi>r</mi> </msub> <mo>=</mo> <mfrac> <mi>N</mi> <mrow> <mi>M</mi> <mo>&amp;times;</mo> <mi>T</mi> </mrow> </mfrac> <mo>&amp;times;</mo> <mi>P</mi> <mo>&amp;times;</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
在公式(9)中，M为光电式旋转编码器旋转一周所产生的脉冲数，P为开绕组异步电机极对数，T为预设采样周期。
4.根据权利要求1所述的一种共直流母线开绕组异步电机系统零序环流抑制方法，其特征是步骤2中转子定向角度θ的获得方式为：
θ＝∫ωrdt (10)
同时限制θ的范围为0°-360°；其中，ωr为转子电角速度，dt为时间的微分。
5.根据权利要求1所述的一种共直流母线开绕组异步电机系统零序环流抑制方法，其特征是步骤3中d轴电流调节器、q轴电流调节器、零序环流调节器均为比例积分调节器。
CN201710541205.4A 2017-07-05 2017-07-05 共直流母线开绕组异步电机系统零序环流抑制方法 Active CN107196543B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710541205.4A CN107196543B (zh) 2017-07-05 2017-07-05 共直流母线开绕组异步电机系统零序环流抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710541205.4A CN107196543B (zh) 2017-07-05 2017-07-05 共直流母线开绕组异步电机系统零序环流抑制方法

Publications (2)

Publication Number Publication Date
CN107196543A true CN107196543A (zh) 2017-09-22
CN107196543B CN107196543B (zh) 2018-07-27

Family

ID=59880773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710541205.4A Active CN107196543B (zh) 2017-07-05 2017-07-05 共直流母线开绕组异步电机系统零序环流抑制方法

Country Status (1)

Country Link
CN (1) CN107196543B (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107994828A (zh) * 2017-12-29 2018-05-04 武汉英弗耐斯电子科技有限公司 一种永磁同步电机交直轴电感测量方法
CN108418489A (zh) * 2018-03-16 2018-08-17 江苏大学 抑制零序电流的绕组开放式永磁同步电机矢量控制方法
CN108988718A (zh) * 2018-08-15 2018-12-11 华中科技大学 抑制零序电流和共模电压的交流电机驱动系统及方法
CN109861622A (zh) * 2019-04-08 2019-06-07 浙江大学 共直流母线开绕组永磁同步电机系统的低开关频率控制策略
CN109889128A (zh) * 2019-01-31 2019-06-14 华中科技大学 一种双三相绕组交流电机共模噪声抑制方法和系统
CN110112946A (zh) * 2018-01-31 2019-08-09 北京金风科创风电设备有限公司 控制变流器的机侧逆变器的交流输出端电压的方法和装置
CN110932533A (zh) * 2019-12-06 2020-03-27 合肥工业大学 共中线开绕组电机控制变流器拓扑高频共模电压抑制方法
CN111146801A (zh) * 2019-12-06 2020-05-12 合肥工业大学 共直流母线双逆变器光伏发电系统的零序电流抑制方法
CN111181464A (zh) * 2020-02-13 2020-05-19 东南大学 一种基于零序调磁的记忆电机调磁控制方法及系统
CN111224597A (zh) * 2020-02-19 2020-06-02 东华大学 一种开绕组永磁同步电机的零序电流抑制方法
CN111245335A (zh) * 2018-11-28 2020-06-05 广东威灵汽车部件有限公司 逆变器的控制装置、方法、开绕组电机的控制系统
CN111342694A (zh) * 2018-11-28 2020-06-26 广东威灵汽车部件有限公司 逆变器的控制装置、方法和开绕组电机的控制系统
CN112510981A (zh) * 2020-12-08 2021-03-16 合肥工业大学 一种电磁兼容友好型开绕组拓扑结构
CN113364326A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种瞬时误差最小的3d-svpwm调制策略过调制方法及系统
CN113364330A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种输出谐波最小化3d-svpwm调制策略过调制方法及系统
CN113364332A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-azspwm调制策略快速过调制方法及系统
CN113364328A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-svpwm调制策略过调制输出最大化方法及系统
CN113364329A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-azspwm调制策略过调制输出最大化方法及系统
CN113364325A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-svpwm调制策略快速过调制方法及系统
CN113364331A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种瞬时误差最小的3d-azspwm调制策略过调制方法及系统
CN113364327A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种输出谐波最小化3d-azspwm调制策略过调制方法及系统
CN116155067A (zh) * 2023-04-18 2023-05-23 湖南大学 基于零序环流抑制的开关频率优化控制模型、方法及装置
CN116827210A (zh) * 2023-08-24 2023-09-29 四川大学 一种开绕组电机三维空间矢量调制方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001061283A (ja) * 1999-08-19 2001-03-06 Yaskawa Electric Corp ３相中性点クランプ式ｐｗｍインバータ装置
US20040032755A1 (en) * 2002-06-13 2004-02-19 Riku Pollanen Method in connection with converter bridges
US20050046508A1 (en) * 2003-08-28 2005-03-03 Vacon Oyj Pulse-width modulation method for a frequency converter
CN101291130A (zh) * 2008-06-13 2008-10-22 株洲南车时代电气股份有限公司 一种基于空间矢量的同步调制方法及系统
CN102624276A (zh) * 2012-04-25 2012-08-01 杭州电子科技大学 一种新颖的交流伺服逆变器死区效应补偿的方法
CN102684529A (zh) * 2012-05-30 2012-09-19 无锡市大元广盛电气有限公司 一种空间电压矢量脉宽调制方法
CN104467005A (zh) * 2014-01-02 2015-03-25 艾伏新能源科技(上海)股份有限公司 T型三电平三相四桥臂光伏并网发电系统及其控制方法
CN104539200A (zh) * 2014-12-03 2015-04-22 广东美的制冷设备有限公司 电机的转速控制器、电机的转速控制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001061283A (ja) * 1999-08-19 2001-03-06 Yaskawa Electric Corp ３相中性点クランプ式ｐｗｍインバータ装置
US20040032755A1 (en) * 2002-06-13 2004-02-19 Riku Pollanen Method in connection with converter bridges
US20050046508A1 (en) * 2003-08-28 2005-03-03 Vacon Oyj Pulse-width modulation method for a frequency converter
CN101291130A (zh) * 2008-06-13 2008-10-22 株洲南车时代电气股份有限公司 一种基于空间矢量的同步调制方法及系统
CN102624276A (zh) * 2012-04-25 2012-08-01 杭州电子科技大学 一种新颖的交流伺服逆变器死区效应补偿的方法
CN102684529A (zh) * 2012-05-30 2012-09-19 无锡市大元广盛电气有限公司 一种空间电压矢量脉宽调制方法
CN104467005A (zh) * 2014-01-02 2015-03-25 艾伏新能源科技(上海)股份有限公司 T型三电平三相四桥臂光伏并网发电系统及其控制方法
CN104539200A (zh) * 2014-12-03 2015-04-22 广东美的制冷设备有限公司 电机的转速控制器、电机的转速控制方法

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107994828A (zh) * 2017-12-29 2018-05-04 武汉英弗耐斯电子科技有限公司 一种永磁同步电机交直轴电感测量方法
CN110112946A (zh) * 2018-01-31 2019-08-09 北京金风科创风电设备有限公司 控制变流器的机侧逆变器的交流输出端电压的方法和装置
CN110112946B (zh) * 2018-01-31 2020-11-10 北京金风科创风电设备有限公司 控制变流器的机侧逆变器的交流输出端电压的方法和装置
CN108418489B (zh) * 2018-03-16 2020-06-09 江苏大学 抑制零序电流的绕组开放式永磁同步电机矢量控制方法
CN108418489A (zh) * 2018-03-16 2018-08-17 江苏大学 抑制零序电流的绕组开放式永磁同步电机矢量控制方法
CN108988718A (zh) * 2018-08-15 2018-12-11 华中科技大学 抑制零序电流和共模电压的交流电机驱动系统及方法
CN108988718B (zh) * 2018-08-15 2020-10-02 华中科技大学 抑制零序电流和共模电压的方法
CN111342694A (zh) * 2018-11-28 2020-06-26 广东威灵汽车部件有限公司 逆变器的控制装置、方法和开绕组电机的控制系统
CN111245335A (zh) * 2018-11-28 2020-06-05 广东威灵汽车部件有限公司 逆变器的控制装置、方法、开绕组电机的控制系统
CN109889128A (zh) * 2019-01-31 2019-06-14 华中科技大学 一种双三相绕组交流电机共模噪声抑制方法和系统
CN109861622A (zh) * 2019-04-08 2019-06-07 浙江大学 共直流母线开绕组永磁同步电机系统的低开关频率控制策略
CN110932533A (zh) * 2019-12-06 2020-03-27 合肥工业大学 共中线开绕组电机控制变流器拓扑高频共模电压抑制方法
CN110932533B (zh) * 2019-12-06 2021-08-10 合肥工业大学 共中线开绕组电机控制变流器拓扑高频共模电压抑制方法
CN111146801B (zh) * 2019-12-06 2022-10-04 合肥工业大学 共直流母线双逆变器光伏发电系统的零序电流抑制方法
CN111146801A (zh) * 2019-12-06 2020-05-12 合肥工业大学 共直流母线双逆变器光伏发电系统的零序电流抑制方法
CN111181464B (zh) * 2020-02-13 2021-06-25 东南大学 一种基于零序调磁的记忆电机调磁控制方法及系统
CN111181464A (zh) * 2020-02-13 2020-05-19 东南大学 一种基于零序调磁的记忆电机调磁控制方法及系统
CN111224597A (zh) * 2020-02-19 2020-06-02 东华大学 一种开绕组永磁同步电机的零序电流抑制方法
CN112510981A (zh) * 2020-12-08 2021-03-16 合肥工业大学 一种电磁兼容友好型开绕组拓扑结构
CN113364331A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种瞬时误差最小的3d-azspwm调制策略过调制方法及系统
CN113364332A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-azspwm调制策略快速过调制方法及系统
CN113364328A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-svpwm调制策略过调制输出最大化方法及系统
CN113364329A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-azspwm调制策略过调制输出最大化方法及系统
CN113364325A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种3d-svpwm调制策略快速过调制方法及系统
CN113364330A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种输出谐波最小化3d-svpwm调制策略过调制方法及系统
CN113364327A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种输出谐波最小化3d-azspwm调制策略过调制方法及系统
CN113364328B (zh) * 2021-07-01 2022-09-16 合肥工业大学 一种3d-svpwm调制策略过调制输出最大化方法及系统
CN113364326A (zh) * 2021-07-01 2021-09-07 合肥工业大学 一种瞬时误差最小的3d-svpwm调制策略过调制方法及系统
CN113364329B (zh) * 2021-07-01 2022-12-02 合肥工业大学 一种3d-azspwm调制策略过调制输出最大化方法及系统
CN116155067A (zh) * 2023-04-18 2023-05-23 湖南大学 基于零序环流抑制的开关频率优化控制模型、方法及装置
CN116827210A (zh) * 2023-08-24 2023-09-29 四川大学 一种开绕组电机三维空间矢量调制方法
CN116827210B (zh) * 2023-08-24 2023-11-21 四川大学 一种开绕组电机三维空间矢量调制方法

Also Published As

Publication number Publication date
CN107196543B (zh) 2018-07-27

Similar Documents

Publication Publication Date Title
CN107196543A (zh) 共直流母线开绕组异步电机系统零序环流抑制方法
CN101204003B (zh) 电力变换控制装置以及电力变换控制方法
CN107565868B (zh) 一种五相永磁同步电机开路故障下的容错控制系统及方法
CN103780191B (zh) 开绕组永磁同步电机串联补偿矢量控制方法
CN103997267B (zh) 一种开绕组永磁同步电机的串联补偿直接转矩控制方法
Poddar et al. Sensorless field-oriented control for double-inverter-fed wound-rotor induction motor drive
CN103138671B (zh) 一种永磁同步电机逆变器死区效应补偿方法及系统
CN110022106B (zh) 一种基于高频信号注入的永磁同步电机无位置传感器控制方法
CN107005194A (zh) 多绕组电动机驱动控制装置
CN107231111B (zh) 五相永磁同步电机高负载能力的直接转矩控制方法
CN101286724A (zh) 一种永磁同步电机直接转矩控制方法
CN110299882B (zh) 混合供电型开绕组永磁同步电机三矢量模型预测控制方法
CN103441726A (zh) 基于比例谐振调节器的双三相永磁电机矢量控制方法
CN107689760A (zh) 基于矩阵变换器的永磁电机无位置矢量控制系统与方法
CN101145751A (zh) 交流励磁双馈异步风力发电机转子电流无延时控制方法
CN103401503A (zh) 一种在谐波平面在线辨识双三相电机参数的方法及装置
CN102204082A (zh) 永磁体同步电动机的控制装置
CN105610369A (zh) 一种基于滑模观测器的异步电机磁链观测方法
CN107196571A (zh) 一种双电机串联预测型直接转矩控制方法
CN111431453A (zh) 基于逆变器单相故障的开绕组永磁同步电机容错控制方法
CN107204727A (zh) 一种小容量薄膜电容永磁同步电机直轴电流给定控制方法
CN104753419B (zh) 一种基于电压切割法的开绕组永磁同步电机系统控制方法
CN109981017B (zh) 一种开路故障下的两单元同相位永磁同步电机容错控制系统及方法
CN108054962A (zh) 一种无刷直流电机无位置传感器控制方法
CN109347392B (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