CN116317753B - 一种矩峰同位型电机弱磁控制方法及*** - Google Patents
一种矩峰同位型电机弱磁控制方法及*** Download PDFInfo
- Publication number
- CN116317753B CN116317753B CN202310286604.6A CN202310286604A CN116317753B CN 116317753 B CN116317753 B CN 116317753B CN 202310286604 A CN202310286604 A CN 202310286604A CN 116317753 B CN116317753 B CN 116317753B
- Authority
- CN
- China
- Prior art keywords
- motor
- current
- axis
- voltage
- torque
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 101000841267 Homo sapiens Long chain 3-hydroxyacyl-CoA dehydrogenase Proteins 0.000 claims abstract description 63
- 102100029107 Long chain 3-hydroxyacyl-CoA dehydrogenase Human genes 0.000 claims abstract description 63
- JJYKJUXBWFATTE-UHFFFAOYSA-N mosher's acid Chemical compound COC(C(O)=O)(C(F)(F)F)C1=CC=CC=C1 JJYKJUXBWFATTE-UHFFFAOYSA-N 0.000 claims abstract description 63
- 230000003313 weakening effect Effects 0.000 claims abstract description 61
- 230000004907 flux Effects 0.000 claims abstract description 58
- 238000012937 correction Methods 0.000 claims description 45
- 230000003247 decreasing effect Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 230000003068 static effect Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 14
- 230000005389 magnetism Effects 0.000 description 13
- 230000000670 limiting effect Effects 0.000 description 9
- 238000004088 simulation Methods 0.000 description 5
- 230000005347 demagnetization Effects 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000011478 gradient descent method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
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/0085—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for high speeds, e.g. above nominal speed
- H02P21/0089—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for high speeds, e.g. above nominal speed using field weakening
-
- 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/18—Estimation of position or speed
-
- 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/22—Current control, e.g. using a current control loop
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明提出了一种矩峰同位型电机弱磁控制方法及***,根据测量的矩峰同位型电机的转子位置角和给定的目标转速,确定电磁转矩给定值;根据所述矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,求解所述目标函数的极值点,通过给定所述电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值;根据逆变器极限承受电压与定子电压的差值判断是否需要开启弱磁控制,开启时则修正交、直轴电流给定值;根据所述交、直轴电流给定值和获取的交、直轴电流实际值,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转,无退磁风险。
Description
技术领域
本发明属于电机驱动相关技术领域,尤其涉及一种矩峰同位型电机弱磁控制方法及***。
背景技术
本部分的陈述仅仅是提供了与本发明相关的背景技术信息,不必然构成在先技术。
近年来,具有矩峰同位特征的电机受到研究者的关注,该类电机利用非对称转子结构,使永磁转矩最大值与磁阻转矩最大值在相同的电流相位角处叠加,以达成增加输出转矩的目的。然而,该新型电机的拓扑结构会使转子主磁链矢量超前磁极中心线45°电角度,因此,按照凸极式永磁电机转子交直轴(dq)坐标定向规则(d轴定向在转子磁极中心线处,q轴超前d轴90°电角度)进行定向,会产生d轴与转子主磁链矢量互差45°的现象,该现象在传统凸极式永磁电机中不会出现。因此,矩峰同位型电机的数学模型方程与传统凸极式电机不同,需要匹配不同的控制***。
另一方面,近年来,针对传统凸极式永磁电机所设计的弱磁控制***已较为成熟。为拓宽永磁电机的调速范围,通常需要在电机运行在高速段时通入弱磁电流,减小电机产生的反电势。然而,永磁体在电机进行弱磁控制时,若弱磁电流过大,则有可能出现退磁的问题,因此,不会引起永磁体退磁问题的电机弱磁控制***具有更可靠的性能,这就要求永磁体始终处于助磁状态。然而,经研究表明,传统隐极式或凸极式电机无法兼顾永磁体助磁与最大转矩电流比(MTPA)控制,也不能兼顾永磁体助磁与最大转矩电压比(MTPV)控制,而矩峰同位型电机的出现为该问题的解决提供了条件。然而,由于矩峰同位型电机具有与传统凸极式永磁电机不同的数学模型方程,传统电机的弱磁控制***无法适用于矩峰同位型电机。
发明内容
为克服上述现有技术的不足,本发明提供了一种矩峰同位型电机弱磁控制方法,能够考虑该电机的特殊拓扑结构,拓宽该电机的调速范围,同时提供的弱磁控制***及方法无需考虑永磁体退磁的问题。
本发明都一个方面提供了一种矩峰同位型电机弱磁控制方法,包括:
根据测量的矩峰同位型电机的转子位置角和给定的目标转速,确定电磁转矩给定值;
根据所述矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,求解所述目标函数的极值点,通过给定所述电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值;
根据逆变器极限承受电压与定子电压的差值判断是否需要开启弱磁控制,开启时则修正交、直轴电流给定值;
根据所述交、直轴电流给定值和获取的交、直轴电流实际值,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转。
本发明的第二个方面提供一种矩峰同位型电机弱磁控制***,包括:角速度处理模块、转速PI控制器、MTPA模块、弱磁修正模块、三相静止坐标系-两相旋转坐标系变换器、d轴/q轴电流PI控制器、两相旋转坐标系-两相静止坐标系变换器、SVPWM模块、逆变桥;
所述转速PI控制器接收设定目标转速值和所述角速度处理模块传来的电机转子转速值作为负反馈值进行计算形成电磁转矩给定值;
所述MTPA模块根据转速PI控制器传来的电磁转矩给定值计算得到MTPA条件下电机交、直轴电流的给定值;
所述弱磁修正模块根据逆变器极限承受电压值与定子电压的差值判断是否需要开启弱磁,若开启则给所述MTPA模块发送信号,并修正MTPA模块输出的交、直轴给定电流值;
所述三相静止坐标系-两相旋转坐标系变换器利用电角度将电机输入的三相电流值变换到两相旋转坐标系下,得到交、直轴实际电流值;
所述d轴/q轴电流PI控制器接收交、直轴电流实际值和所述MTPA模块传来的电机dq轴电流的给定值进行计算得到交、直轴电压给定值;
所示两相旋转坐标系-两相静止坐标系变换器利用电角度将所述交、直轴电压给定值变换到两相静止坐标系下,得到交、直轴电压值;
所述SVPWM模块根据所述交、直轴电压值生成三相PWM信号;
所述逆变桥接收所述SVPWM模块的PWM信号生成三相电压值驱动电机运转。
以上一个或多个技术方案存在以下有益效果:
本公开提出的弱磁控制***,在输出负转矩的同时,电机电流始终处于助磁状态,因此无需考虑永磁体退磁问题。
本公开通过在MTPA模块及电流修正模块中方程求极值得到的是转矩电流比或转矩电压比的最大值,从而该模块可以实现最大转矩电流比及最大转矩电压比控制。
本公开的控制***并非只适用于矩峰同位型电机,凡是利用非对称转子结构使电机永磁转矩最大值与磁阻转矩最大值在相同或相近电流相位角处叠加的思想所设计的电机均可采用。
本公开的电流给定模块可应用至模型预测控制***,以及上述基于模型预测控制***与弱磁控制***的无传感器控制***设计。
本发明附加方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1是本发明实施例一提供的永磁体助磁状态下矩峰同位型电机的弱磁控制***控制原理图;
图2(a)是本发明中矩峰同位型电机每极拓扑结构图;
图2(b)是本发明中矩峰同位型电机每极dq轴定向位置与转子主磁链位置;
图3(a)是本发明实施例一中的电机中各转矩分量随电流相位角的变化关系;
图3(b)传统凸极式永磁电机中各转矩分量随电流相位角的变化关系;
图4是本发明实施例一中的电机各电气量时空统一矢量图;
图5是本发明实施例一中的MTPA及恒转矩方向曲线图;
图6是本发明实施例一中的MTPV方向示意图;
图7是本发明实施例一中的弱磁控制d、q轴电流给定值示意图;
图8是本发明实施例一中的仿真转速图;
图9是本发明实施例一中的仿真输出转矩图;
图10是本发明实施例一中的仿真d、q轴实际电流散点图;
图11是本发明实施例二所提供的矩峰同位型电机的弱磁控制方法应用于电机时的控制原理图;
图12(a)是本发明实施例二中的电机中各转矩随电流相位角的变化关系;
图12(b)传统凸极式永磁电机中各转矩随电流相位角的变化关系;
图13是本发明实施例二中的电机各电气量空间矢量图;
图14是本发明实施例二中的MTPA及恒转矩方向曲线图;
图15是本发明实施例二中的MTPV方向示意图;
图16是本发明实施例二中的弱磁控制d、q轴电流给定值示意图;
图17是本发明实施例二中的仿真转速图;
图18是本发明实施例二中的仿真输出转矩图;
图19是本发明实施例二中的仿真d、q轴实际电流散点图;
图20是本发明实施例二中的弱磁控制***与一般电机控制***的仿真对比图;
图21是本发明实施例二中的通过输出负转矩减速的仿真转速图;
图22是本发明实施例二中的通过输出负转矩减速的仿真输出转矩图;
图23是本发明实施例二中的通过输出负转矩减速的仿真d、q轴实际电流散点图。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。
在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
实施例一
矩峰同位型电机,是一类新型电机,通过使用非对称转子结构,将电机永磁转矩最大值与磁阻转矩最大值在相同或相近电流相位角处叠加,以达成使电机在通入相同幅值的正弦电流时,能输出更大的转矩的目的。此外,该电机每极中心线与转子主磁链矢量位置不重合,导致推导出与传统凸极式永磁电机不同的转矩方程。
如图1-图10所示,本实施例提出一种矩峰同位型电机弱磁控制方法,实现对矩峰同位型电机的控制,包括:根据测量的矩峰同位型电机的转子位置角和给定的目标转速,确定电磁转矩给定值;
根据所述矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,求解所述目标函数的极值点,通过给定所述电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值;
根据逆变器极限承受电压与定子电压的差值判断是否需要开启弱磁控制,开启时则修正交、直轴电流给定值;
根据所述交、直轴电流给定值和获取的交、直轴电流实际值,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转。
需要说明的是,本实施例所提出的控制***并非只适用于本实施例使用的矩峰同位型电机,凡是具有电机永磁转矩最大值与磁阻转矩最大值在相同或相近电流相位角处叠加特点的电机,均可采用本实施例所提出的控制***。
本实施例以矩峰同位型电机为例进行说明,如图1所示,包括:矩峰同位型电机1、直流电源2、逆变桥3、ABC-dq变换器4、q轴电流PI控制器5、d轴电流PI控制器6、dq-αβ变换器7、SVPWM模块8、光电编码器9,角速度处理模块10、转速PI控制器11、MTPA模块12、电流修正模块13、弱磁判断模块14、电流限幅模块15。
其中,光电编码器9安装在矩峰同位型电机的转子轴上,测量矩峰同位型电机的转子机械位置角θm,并分别送入角速度处理模块。
角速度处理模块10,用于根据光电编码器9测量的转子位置角θm,微分计算得到转速ωr。
转速PI控制器,用于根据角速度处理模块计算的转速ωr,以及给定的目标转速ωr *,计算得到电磁转矩给定值。
MTPA模块12,用于根据电磁转矩给定值,依托辅助函数,得到MTPA条件下电机d、q轴电流的给定值id *、iq *。
弱磁判断模块14,根据逆变器极限承受电压Usmax与定子电压Us的差值ΔU判断是否进行弱磁,即开通电流修正模块。
电流修正模块13,用于根据逆变器极限承受电压Usmax与定子电压Us的差值ΔU,计算得到电机d、q轴电流的修正值。
电流限幅模块15,用于将电机d、q轴电流的给定值id *、iq *经电流修正模块修正后的值限制在电机电流最大值内。
ABC-dq变换器4,用于利用电角度θe,将由电流互感器得到的电机输入的三相电流值变换到dq坐标系下,得到d轴与q轴实际电流值id和iq。
d轴电流PI控制器6,用于根据前述电机d轴电流的给定值id *以及d轴电流实际值id,计算得到d轴电压给定值ud *。
q轴电流PI控制器5,用于根据前述电机q轴电流的给定值iq *以及q轴电流实际值iq,计算得到q轴电压给定值uq *。
dq-αβ变换器7,用于利用电角度θe,将电压给定值ud *、uq *由d-q坐标系变换到α-β坐标系下,得到uα和uβ。
SVPWM模块8,用于基于电压给定uα和uβ得到三相的PWM信号,并送入逆变桥模块。
逆变桥模块,连接直流电压源与矩峰同位型电机,用于根据前述三相PWM信号生成三相电压值,驱动电机运转。
本实施例提供一种矩峰同位型电机弱磁控制方法,包括:
根据测量的矩峰同位型电机的转子位置角和给定的目标转速,确定电磁转矩给定值;
根据所述矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,求解所述目标函数的极值点,通过给定所述电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值;
弱磁控制根据实际电流值大小判断是否需要开启,开启时修正给定交直轴电流,而后根据所述交、直轴电流给定值和获取的交、直轴电流实际值,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转。
矩峰同位型电机每极拓扑结构如图2(a)所示,该电机每个极上有四个内置永磁体,镶嵌在转子槽内;每个极上dq轴定向位置如图2(b)所示,转子主磁链矢量的方向与位置位于dq轴的角平分线处。
矩峰同位型电机特有的拓扑结构使得永磁转矩TPM与磁阻转矩Tre最大值在相同电流相位角处叠加得到总转矩Tem,如图3(a)所示,相较于传统凸极式永磁电机(如图3(b)所示),总转矩最大值有所提升,使得电机具有更大的转矩密度。
根据磁力线分布,得出由图4所示电机空间矢量图,其中is为定子电流空间矢量,id、iq分别为is的交直轴分量,ψPM为永磁体产生的永磁磁链,永磁磁链偏移角为45°,ψ0为is产生的磁链,ψs为ψ0与ψPM的合成磁链。
则矩峰同位型电机的电磁转矩方程为:
其中,p为电机极对数,Ld为d轴电感,Lq为q轴电感,Te为电磁转矩。
如图5中虚线方向,则矩峰同位型电机恒转矩方向即电磁转矩方程的切线,可表示为:
其中,为求偏导数,
在本实施例的MTPA模块中,为求得电流和转矩的关系,即求在某个Te下的最小值,为求得在该条件下的iq与id的关系,利用拉格朗日极值定理,利用拉格朗日乘子λ引入辅助函数,得:
根据拉格朗日极值定理,要求的iq与id的关系为上述函数F的极值点,即:
舍去增根,解得:
为保证所述弱磁控制***永磁体始终处在助磁状态下,加入了以下约束条件:
id+iq>0 (6)
为直观表示上述的推导过程,绘制如图5所示的MTPA曲线图,在id-iq平面中的恒转矩曲线如虚线所示,上述问题实际上就是求解每条虚线上与原点最近的点,将这些点连接起来,就得到了MTPA曲线,也就是由拉格朗日方程公式(3)解得的根。
其中,MTPA曲线仅包含助磁区的部分。
在本实施例弱磁判断模块中,为实现自动弱磁,通过实时判断逆变器极限承受电压Usmax和定子电压Us的差值ΔU实现,具体过程如下:
其中,Udc为直流母线电压。
矩峰同位型电机的定子电压Us满足以下关系:
其中,忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,其稳态下电机d轴电压ud和q轴电压uq可表示为:
故有
其中,ωr为电机电角速度。
则ΔU为:
ΔU=Usmax-Us (11)
当ΔU<0时,进行弱磁,即开通电流修正模块。
在本实施例的电流修正模块中,需要确立电机恒转矩方向,需要确立电机电压递减方向。
其中,通过梯度下降法,确立电机电压递减方向,具体过程如下:具体为:
为了确定电压下降方向,设计代价函数S:
其中,忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,其稳态下电机d轴电压ud和q轴电压uq可表示为:
电压递减方向可由梯度下降法表示,表示S的梯度,满足:
电压递减方向即由(Ud,Uq)表示。
在电流修正模块中,需求的电机MTPV方程,需要先求得电机电压极限椭圆,推导过程如下:
忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,稳态下矩峰同位型电机d轴和q轴电压可表示为:
此时,由于逆变器两端不能无限制地承受电压,有如下的电压方程限制:
即:
其中,Umax为逆变器极限承受电压。
该电压限制方程为一椭圆,其椭圆中心为其长轴顶点为其短轴顶点为/>
电机MTPV方向为电机电压极限椭圆与电机转矩双曲线切点的连线,推导过程如下:
电机MTPV方程可由下式表示:
其中,矩峰同位型电机的电磁转矩Te方程为:
矩峰同位型电机的定子电压Us方程为:
则有,
即:
MTPV方向为式(22)的切线方向,即:
其中,
为保证所述弱磁控制***永磁体始终处在助磁状态下,本公开加入了以下约束条件:
id+iq>0 (24)
为了直观地表示上述推导过程,绘制如图6所示的MTPV曲线图,在id-iq平面中的恒转矩曲线如虚线所示,上述问题实际上就是求解每条虚线上与电机不同转速(ωr)下电机电压极限圆相切的点,将这些点连接起来,就得到了MTPV曲线。
其中,MTPV方程仅取助磁区的部分,此时其转矩输出为负,MTPV曲线仅包含助磁区的部分。
在电流修正模块中,弱磁区域由以下过程确定:
令:
其中,θ表示恒转矩方向与电压下降方向的夹角,
当cosθ>0时,电流修正模块运行在弱磁区域1(FW1),沿恒转矩方向给定修正电流,当cosθ<0时,电流修正模块运行在弱磁区域2(FW2),沿MTPV方向给定修正电流。
则电流修正模块给定的d轴和q轴弱磁电流修正值如下:
其中,α,β为修正系数。
弱磁开通后,电机电流的给定值为MTPA模块的电流给定值和电流修正模块给定的弱磁电流修正值之和。
本实施例的弱磁控制d、q轴电流给定值如图7所示,本实施例中,电机运行在-1700rpm及以下时,电机不开通弱磁,电机d、q轴给定电流仅由MTPA模块给定;当电机运行速度超过-1700rpm时,电机弱磁模块开通,此时电机的d、q轴给定电流为图7中的恒转矩段,即弱磁区域1;当电机加速至-1900rpm以上时,此时电机的d、q轴给定电流为图7中的MTPV段,即弱磁区域2。
图8、9、10给出了本实施例的仿真波形,本实施例中Udc取200V。
如图8所示,电机在0s时给定-1700rpm转速并使电机空载运行,电机根据MTPA模块恒定输出-2Nm的转速升至-1700rpm,电机不开通弱磁,电机的d、q轴给定电流处在MTPA段;1s时,给定-1900rpm转速并使电机空载运行,此时电机开通弱磁,电机恒定输出-2Nm的电磁转矩,电机的d、q轴给定电流处在恒转矩段;1.5s时,给定-2900rpm并使电机空载运行,电机开通弱磁,并根据MTPV方向单位电压输出最大转矩,电机的d、q轴给定电流处在MTPV段。上述所有的d轴电流和q轴电流均如图10中散点所示。
本发明限制d、q轴给定电流的目的是为使电机一直工作在永磁体助磁状态下,防止电机永磁体退磁。
实施例二
如图11-图23所示,本实施例提供一种矩峰同位型电机的弱磁控制方法,如图11所示,将光电编码器9安装在矩峰同位型电机的转子轴上,测量矩峰同位型电机的转子位置角θm,并进行角速度处理得到实际转速ωr,可选的,对转子位置角进行微分处理,得到所述矩峰同位型电机的实际转速。转速PI控制器11根据角速度处理模块10计算的转速ωr,以及给定的目标转速ωr *,计算得到电磁转矩给定值。
MTPA模块12根据矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,求解所述目标函数的极值点,进而通过给定电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值。
具体地,矩峰同位型电机每极拓扑结构如图2(a)所示,该电机每个极上有四个内置永磁体,分别镶嵌在转子槽内;每个极上dq轴(即交直轴)定向位置如图2(b)所示,该图中还标注了转子主磁链矢量的方向与位置,即在dq轴的角平分线处。
该电机特有的拓扑结构使得永磁转矩TPM与磁阻转矩Tre最大值在相同电流相位角处叠加得到总转矩Tem,如图12(a)所示,相较于传统凸极式永磁电机(如图12(b)所示),总转矩最大值有所提升,使得电机具有更大的转矩密度。
根据磁力线分布,得出如图13所示的电机空间矢量图,其中is为定子电流空间矢量,id、iq分别为交轴电流和直轴电流(也即is的交、直轴分量),ψPM为永磁体产生的永磁磁链,永磁磁链偏移角为45°,ψ0为is产生的磁链,ψs为ψ0与ψPM的合成磁链。
则矩峰同位型电机的电磁转矩方程为:
其中,p为电机极对数,Te为电磁转矩,Ld为直轴电感,Lq为交轴电感。
则矩峰同位型电机恒转矩方向,即电磁转矩方程的切线,可表示为:
即沿图14中虚线的切线方向。
要求MTPA(最大转矩电流比)模块中电流和转矩的关系,即求在某个Te下的最小值,为求得在该条件下的iq与id的关系,利用拉格朗日极值定理,利用拉格朗日乘子λ引入辅助函数,可得目标函数F为:
根据拉格朗日极值定理,要求的iq与id的关系为上述目标函数F的极值点,从而有:
舍去增根,解得
为了保证弱磁性能,加入了以下约束条件:
id+iq≤0 (33)
为了直观的表示上述推导过程,绘制如图14所示的MTPA曲线图,在id-iq平面中的恒转矩曲线如虚线所示,上述问题实际上就是求解每条虚线上与原点最近的点,将这些点连接起来,就得到了MTPA曲线,也就是由拉格朗日方程解得的根,从而得到符合MTPA条件的交、直轴电流给定值。
作为一可选实施方式,为了实现自动弱磁,通过实时判断逆变器极限承受电压Usmax和定子电压Us的差值ΔU实现,具体过程如下:
逆变器极限承受电压Usmax满足以下关系:
其中,Udc为直流母线电压。
根据所述矩峰同位型电机的交、直轴电流实际值确定定子电压,在具体实施中,通过电流互感器测量矩峰同位型电机的三相电流值,ABC-dq变换器4利用电角度θe,将由电流互感器得到的三相电流值变换到dq坐标系下,得到d轴与q轴实际电流值id和iq。并根据d轴与q轴实际电流值确定定子电压,矩峰同位型电机的定子电压Us满足以下关系:
其中,忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,其稳态下电机d轴电压ud和q轴电压uq可表示为:
故有
定义ΔU为:ΔU=Usmax-Us,根据定子电压与直流母线电压值之间的差值,判定是否进行弱磁控制。
具体地,当ΔU<0时,进行弱磁,即开通电流修正模块。
进一步地,通过梯度下降法,确立电机电压递减方向,具体过程如下:
为了确定电压下降方向,设计代价函数S:
其中,忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,其稳态下电机d轴电压ud和q轴电压uq可表示为:
电压递减方向可由梯度下降法表示,▽S表示S的梯度,满足:
电压递减方向即由(Ud,Uq)表示。
要求得电机MTPV方程,需要求得电机电压极限椭圆,其推导过程如下:
忽略定子电阻,由于矩峰同位型电机的永磁磁链ψPM超前d轴45°,稳态下矩峰同位型电机d轴和q轴电压可表示为:
此时,由于逆变器两端不能无限制地承受电压,有如下的电压方程限制:
即
该电压限制方程为一椭圆,其椭圆中心为其长轴顶点为其短轴顶点为/>电机MTPV方程为电机电压极限椭圆与电机转矩曲线切点的连线,其推导过程如下:电机MTPV方程可由下式表示
其中,矩峰同位型电机的电磁转矩Te方程为:
矩峰同位型电机的定子电压Us方程为:
则有,
即
MTPV方向为上式的切线方向,即
为了保证弱磁性能,加入了以下约束条件:
id+iq≤0 (50)
为了直观的表示上述推导过程,绘制如图15所示的MTPV曲线图,在id-iq平面中的恒转矩曲线如虚线所示,上述问题实际上就是求解每条虚线上与电机不同转速(ωr)下电机电压极限圆相切的点,将这些点连接起来,就得到了MTPV曲线。
其中,电流修正模块的弱磁区域由以下过程确定:
令
这里,θ表示恒转矩方向与电压下降方向的夹角;
其中,
当cosθ>0时,电流修正模块运行在弱磁区域1(FW1),沿恒转矩方向给定修正电流,当cosθ<0时,电流修正模块运行在弱磁区域2(FW2),沿MTPV方向给定修正电流;
电流修正模块给定的d轴和q轴弱磁电流修正值如下:
这里,α,β为修正系数。
弱磁开通后,电机电流的给定值为MTPA模块的电流给定值和电流修正模块给定的交、直轴弱磁修正电流之和。
可选的,在对交、直轴电流给定值进行修正之前,还通过电流限幅模块15将电机d、q轴电流给定值限制在电机最大电流内。
进一步的,由电流互感器测量得到矩峰同位型电机1的三相电流值,ABC-dq变换器4利用电角度θe,将三相电流值变换到dq坐标系下,得到d轴与q轴实际电流值id和iq;d轴电流PI控制器6根据前述电机d轴电流的给定值id *以及d轴电流实际值id,计算得到d轴电压给定值ud *;q轴电流PI控制器5根据前述电机q轴电流的给定值iq *以及q轴电流实际值iq,计算得到q轴电压给定值uq *;dq-αβ变换器7利用电角度θe,将交、直轴电压给定值由d-q坐标系变换到α-β坐标系下,得到uα和uβ;交、直轴电压给定值uα和uβ经过信号生成模块8(即SVPWM模块)生成空间矢量脉宽调制信号(即三相的PWM信号),并送入逆变桥3,逆变桥3连接有直流电压源2和矩峰同位型电机1,用于根据前述三相的PWM信号生成三相电压值,驱动电机运转。
本实施例的弱磁控制d、q轴电流给定值如图16所示,电机运行在700rpm及以下时,电机不开通弱磁,电机d、q轴给定电流仅由MTPA模块给定;当电机运行至700rpm以上时,电机弱磁模块开通,此时电机的d、q轴给定电流为图16中的恒转矩段,即弱磁区域1;当电机加速至1100rpm以上时,收到电流限幅模块的限制,电机d、q轴给定电流为图16中的极限电流圆段,当电机加速至3600rpm以上时,此时电机的d、q轴给定电流为图16中的MTPV段,即弱磁区域2。
图17、18、19给出了本实施例的仿真波形,本实施例中Udc取100V。
如图17所示,电机在0s时给定700rpm转速,电机根据MTPA模块恒定输出5Nm的转矩,转速升至700rpm,电机的d、q轴给定电流处在MTPA段;随后在0.5s带3Nm负载,此时电机根据MTPA模块恒定输出3Nm的电磁转矩,电机的d、q轴给定电流处在MTPA段;1s时,给定1100rpm并使电机空载运行,此时电机开通弱磁,电机恒定输出4Nm的电磁转矩,电机的d、q轴给定电流处在恒转矩段;1.5s时,给定3600rpm并使电机空载运行,电机开通弱磁,此时电机由于电流极限圆的限制,输出转矩不断减少,电机的d、q轴给定电流处在极限电流圆段;3s时,给定9000rpm并使电机空载运行,电机开通弱磁,并根据MTPV方向单位电压输出最大转矩,电机的d、q轴给定电流处在MTPV段。上述所有的d轴电流和q轴电流均如图19中散点所示。
仿真的结果表明,本实施例设计的弱磁控制***性能良好,能实现转速的无差跟踪,且电机的d、q轴实际电流与预期给定值相吻合。
图20给出了采用本实施例提供的弱磁控制***与不采用本实施例提供的弱磁***的电机调速范围比较图,本实施例中Udc取100V。如图20所示,当电机输出2Nm的负载时,采用本公开弱磁控制***能使电机的调速范围从1660rpm拓宽至3070rpm,体现了本实施例提供方法在具体应用中的良好性能。
同样的,该方法同样适用于输出负转矩的工况,图21、22、23给出了本实施例输出负转矩进行电机减速的仿真波形,同样的,本实施例中Udc取100V。
如图21所示,电机在0s时给定6000rpm转速,电机迅速升至给定转速;随后在3s时给定2700rpm转速,此时电机开通弱磁,电机的d、q轴给定电流处在MTPV段;5s时,给定1800rpm转速,此时电机开通弱磁,此时电机由于电流极限圆的限制,电机输出的负转矩随着转速上升而增加,电机的d、q轴给定电流处在极限电流圆段;6s时,给定1150rpm转速,电机开通弱磁,并恒定输出-2Nm的负转矩,电机的d、q轴给定电流处在恒转矩段;7s时,给定0rpm转速,电机不开通弱磁,恒定输出-2Nm的负转矩,电机的d、q轴给定电流处在MTPA段。上述所有的d轴电流和q轴电流均如图23中散点所示。仿真结果显示了,本公开弱磁控制***,同样适用于输出负转矩进行减速的工况。
根据上述求解过程,设计出控制***框图如图11所示。
控制***工作时,首先,通过光电编码器9采集电机的角度信号,用作角速度处理,角速度信号ωr与给定的目标转速ωr *作差,形成负反馈通道,差值信号经转速PI控制器11计算后得到了电机电磁转矩给定值,通过MTPA模块给定MTPA条件下的电流给定值,结合弱磁判断模块14,进入电流修正模块13中得到电机修正电流,将MTPA模块的给定电流与电流修正模块的给定电流相加通过电流限幅模块15得到电机d轴,q轴电流给定值。d轴,q轴电流给定值与经ABC-dq坐标变换后得到的d轴,q轴电流实际值做差,并经过两个电流PI控制器6、5计算得到电压给定值ud *和uq *,再经由dq-αβ坐标变换得到α-β坐标系下的电压给定uα和uβ,经过信号生成模块8后,可以得到ABC三相的PWM信号,将其输入给逆变桥3,即可生成驱动电机所需的三相电压值。
实施例二
本实施例提供一种矩峰同位型电机弱磁控制***,采用包括:
一种矩峰同位型电机弱磁控制***,包括:角速度处理模块、转速PI控制器、MTPA模块、弱磁修正模块、三相静止坐标系-两相旋转坐标系变换器、d轴/q轴电流PI控制器、两相旋转坐标系-两相静止坐标系变换器、SVPWM模块、逆变桥;
所述转速PI控制器接收设定目标转速值和所述角速度处理模块传来的电机转子转速值作为负反馈值进行计算形成电磁转矩给定值;
所述MTPA模块根据转速PI控制器传来的电磁转矩给定值计算得到MTPA条件下电机交、直轴电流的给定值;
所述弱磁修正模块根据逆变器极限承受电压值与定子电压的差值判断是否需要开启弱磁,若开启则给所述MTPA模块发送信号,并修正MTPA模块输出的交、直轴给定电流值;
所述三相静止坐标系-两相旋转坐标系变换器利用电角度将电机输入的三相电流值变换到两相旋转坐标系下,得到交、直轴实际电流值;
所述d轴/q轴电流PI控制器接收交、直轴电流实际值和所述MTPA模块传来的电机dq轴电流的给定值进行计算得到交、直轴电压给定值;
所示两相旋转坐标系-两相静止坐标系变换器利用电角度将所述交、直轴电压给定值变换到两相静止坐标系下,得到交、直轴电压值;
所述SVPWM模块根据所述交、直轴电压值生成三相PWM信号;
所述逆变桥接收所述SVPWM模块的PWM信号生成三相电压值驱动电机运转。
本领域技术人员应该明白,上述本发明的各模块或各步骤可以用通用的计算机装置来实现,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。本发明不限制于任何特定的硬件和软件的结合。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (6)
1.一种矩峰同位型电机弱磁控制方法,其特征在于,包括:
根据测量的矩峰同位型电机的转子位置角和给定的目标转速,确定电磁转矩给定值;
根据所述矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,
所述目标函数F表示为:
其中,λ为拉格朗日乘子,p为电机极对数,Ld为d轴电感,Lq为q轴电感,ψPM为永磁体产生的永磁磁链,Te为电磁转矩;
求解所述目标函数的极值点,通过给定所述电磁转矩给定值,得到符合MTPA条件的交、直轴电流给定值;
根据逆变器极限承受电压与定子电压的差值判断是否需要开启弱磁控制,开启时则修正交、直轴电流给定值;
矩峰同位型电机的定子电压Us方程为:
其中,ωr为电机电角速度;
根据所述交、直轴电流给定值和获取的交、直轴电流实际值,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转;
解所述目标函数的极值点时,为保证永磁体处于助磁条件下,加入的约束条件为:
id+iq>0
其中,iq为交轴电流、id为直轴电流;或,为了保证弱磁性能,加入约束条件为:
id+iq≤0
其中,iq为交轴电流、id为直轴电流;
通过梯度下降法确定电机电压递减方向,进而确定电机恒转矩方向;通过电机电压极限椭圆与电机转矩双曲线切点的连线确定电机的MTPV方向;
电机MTPV方程可由下式表示
MTPV方向为上式的切线方向,即
MTPV方向由(Wd,Wq)表示。
2.如权利要求1所述的矩峰同位型电机弱磁控制方法,其特征在于,根据矩峰同位型电机每极的拓扑结构,确定电磁转矩方程,所述电磁转矩方程表示为:
其中,Te为电磁转矩,p为电机极对数,ψPM为永磁体产生的永磁磁链,iq为交轴电流、id为直轴电流,Ld为直轴电感,Lq为交轴电感。
3.如权利要求1所述的矩峰同位型电机弱磁控制方法,其特征在于,根据恒转矩方向与电压下降方向的夹角的余弦即cosθ确定电流修正方向,当cosθ>0时,沿对应恒转矩方向给定修正电流;当cosθ<0时,沿所述MTPV方向给定修正电流。
4.如权利要求1所述的矩峰同位型电机弱磁控制方法,其特征在于,根据逆变器极限承受电压值与定子电压的差值计算得到电机交、直轴电流的修正值,根据所述电机交、直轴电流的修正值对MTPA条件的交、直轴电流给定值进行电流的修正。
5.如权利要求1所述的矩峰同位型电机弱磁控制方法,其特征在于,根据所述交、直轴电流给定值和获取的交、直轴电流实际值,确定电压给定值;将所述电压给定值进行坐标变换后,生成空间矢量脉宽调制信号,用以驱动所述矩峰同位型电机的运转。
6.一种矩峰同位型电机弱磁控制***,其特征在于,包括:角速度处理模块、转速PI控制器、MTPA模块、弱磁修正模块、三相静止坐标系-两相旋转坐标系变换器、d轴/q轴电流PI控制器、两相旋转坐标系-两相静止坐标系变换器、SVPWM模块、逆变桥;
所述转速PI控制器接收设定目标转速值和所述角速度处理模块传来的电机转子转速值作为负反馈值进行计算形成电磁转矩给定值;
根据矩峰同位型电机的电磁转矩方程构建辅助函数,在满足MTPA条件下利用拉格朗日乘子引入所述辅助函数得到目标函数,
所述目标函数F表示为:
其中,λ为拉格朗日乘子,p为电机极对数,Ld为d轴电感,Lq为q轴电感,ψPM为永磁体产生的永磁磁链,Te为电磁转矩;
所述MTPA模块根据转速PI控制器传来的电磁转矩给定值计算得到MTPA条件下电机交、直轴电流的给定值;
所述弱磁修正模块根据逆变器极限承受电压值与定子电压的差值判断是否需要开启弱磁,若开启则给所述MTPA模块发送信号,并修正MTPA模块输出的交、直轴给定电流值;
所述三相静止坐标系-两相旋转坐标系变换器利用电角度将电机输入的三相电流值变换到两相旋转坐标系下,得到交、直轴实际电流值;
所述d轴/q轴电流PI控制器接收交、直轴电流实际值和所述MTPA模块传来的电机dq轴电流的给定值进行计算得到交、直轴电压给定值;
矩峰同位型电机的定子电压Us方程为:
其中,ωr为电机电角速度;
所示两相旋转坐标系-两相静止坐标系变换器利用电角度将所述交、直轴电压给定值变换到两相静止坐标系下,得到交、直轴电压值;
将交、直轴电压给定值由d-q坐标系变换到α-β坐标系下,得到uα和uβ;
所述SVPWM模块根据uα和uβ生成三相PWM信号;
所述逆变桥接收所述SVPWM模块的PWM信号生成三相电压值驱动电机运转;
在所述MTPA模块中,为保证永磁体处于助磁条件下,加入的约束条件为:
id+iq>0
其中,iq为交轴电流、id为直轴电流;或,为了保证弱磁性能,加入约束条件为:
id+iq≤0
其中,iq为交轴电流、id为直轴电流;
通过梯度下降法确定电机电压递减方向,进而确定电机恒转矩方向;通过电机电压极限椭圆与电机转矩双曲线切点的连线确定电机的MTPV方向;
电机MTPV方程可由下式表示
MTPV方向为上式的切线方向,即
MTPV方向即由(Wd,Wq)表示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310286604.6A CN116317753B (zh) | 2023-03-20 | 2023-03-20 | 一种矩峰同位型电机弱磁控制方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310286604.6A CN116317753B (zh) | 2023-03-20 | 2023-03-20 | 一种矩峰同位型电机弱磁控制方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116317753A CN116317753A (zh) | 2023-06-23 |
CN116317753B true CN116317753B (zh) | 2024-04-26 |
Family
ID=86792298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310286604.6A Active CN116317753B (zh) | 2023-03-20 | 2023-03-20 | 一种矩峰同位型电机弱磁控制方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116317753B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117559857B (zh) * | 2023-11-23 | 2024-06-25 | 山东大学 | 凸极偏置型永磁同步电机的超前角弱磁控制方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651626A (zh) * | 2012-04-24 | 2012-08-29 | 北京交通大学 | 一种永磁同步电机的弱磁控制方法 |
CN106849805A (zh) * | 2017-02-09 | 2017-06-13 | 澳特卡新能源科技(上海)有限公司 | 一种电动压缩机驱动电机的弱磁控制方法 |
CN112865654A (zh) * | 2021-04-13 | 2021-05-28 | 山东大学 | 永磁聚磁式同步磁阻电机转矩最大化利用控制***及方法 |
CN113206625A (zh) * | 2021-05-31 | 2021-08-03 | 大连海事大学 | 一种用于内置式永磁同步电机最大转矩电流比控制方法 |
CN114553082A (zh) * | 2022-03-07 | 2022-05-27 | 祝涛 | 基于内埋式pmsm扭矩提升算法 |
-
2023
- 2023-03-20 CN CN202310286604.6A patent/CN116317753B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651626A (zh) * | 2012-04-24 | 2012-08-29 | 北京交通大学 | 一种永磁同步电机的弱磁控制方法 |
CN106849805A (zh) * | 2017-02-09 | 2017-06-13 | 澳特卡新能源科技(上海)有限公司 | 一种电动压缩机驱动电机的弱磁控制方法 |
CN112865654A (zh) * | 2021-04-13 | 2021-05-28 | 山东大学 | 永磁聚磁式同步磁阻电机转矩最大化利用控制***及方法 |
CN113206625A (zh) * | 2021-05-31 | 2021-08-03 | 大连海事大学 | 一种用于内置式永磁同步电机最大转矩电流比控制方法 |
CN114553082A (zh) * | 2022-03-07 | 2022-05-27 | 祝涛 | 基于内埋式pmsm扭矩提升算法 |
Non-Patent Citations (4)
Title |
---|
Control System for a Novel Dual-Rotor Permanent Magnet Synchronous Reluctance Motor Considering Torque Superposition;Gefei Zhu et al;《2022 25th International Conference on Electrical Machines and Systems (ICEMS)》;第1-5页 * |
内埋式永磁同步电机的弱磁控制策略;唐朝晖 等;《电机与控制学报》;第14卷(第5期);第68-72页 * |
内嵌式永磁同步电机变频调速***及控制策略研究;梁力波;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;第3章 * |
新型永磁聚磁式同步磁阻电机研究;刘炎;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;第5章 * |
Also Published As
Publication number | Publication date |
---|---|
CN116317753A (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9755562B2 (en) | Power convertor, controller, and method for changing carrier frequency | |
US7872435B2 (en) | Motor control apparatus | |
KR101046802B1 (ko) | 교류 회전기의 제어 장치 및 이 제어 장치를 사용한 교류회전기의 전기적 정수 측정 방법 | |
WO2016121751A1 (ja) | インバータ制御装置及びモータ駆動システム | |
CN107968609B (zh) | 永磁同步电机的弱磁控制方法和装置 | |
JP2004343963A (ja) | ブラシレスdcモータの制御装置 | |
JPH1118496A (ja) | 電気車の制御装置および制御方法 | |
CN102780433A (zh) | 一种基于电流控制的无刷直流电机瞬时转矩控制方法 | |
KR102676718B1 (ko) | 모터 회전속도 제어방법 및 제어시스템 | |
Gu et al. | Driving and braking control of PM synchronous motor based on low-resolution hall sensor for battery electric vehicle | |
CN116317753B (zh) | 一种矩峰同位型电机弱磁控制方法及*** | |
EP1681762B1 (en) | Synchronous motor driving system and method | |
CN106533294B (zh) | 一种基于线电压调制的永磁同步电机弱磁控制方法 | |
JP2000156993A (ja) | 永久磁石型同期機の制御装置及びその制御方法 | |
US9240744B2 (en) | Methods, systems and apparatus for adjusting current and/or torque commands used to control operation of an asynchronous machine | |
US9035580B2 (en) | AC motor control apparatus | |
JP2008072858A (ja) | 車両用回転電機の制御装置 | |
JP2013126284A (ja) | 電動機駆動装置 | |
JP5131051B2 (ja) | 回転機の制御装置、及び回転機の制御システム | |
CN112865654B (zh) | 永磁聚磁式同步磁阻电机转矩最大化利用控制***及方法 | |
CN102594254A (zh) | 降低永磁同步电机噪音的***及控制方法 | |
JP2009142073A (ja) | 回転機の制御装置及び制御システム | |
CN111865159B (zh) | 异步电机最大转矩电流比控制方法及*** | |
CN114553082A (zh) | 基于内埋式pmsm扭矩提升算法 | |
CN113141139A (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 |