CN107612436B - 基于永磁电机谐波反电势的转子位置计算方法 - Google Patents
基于永磁电机谐波反电势的转子位置计算方法 Download PDFInfo
- Publication number
- CN107612436B CN107612436B CN201710851734.4A CN201710851734A CN107612436B CN 107612436 B CN107612436 B CN 107612436B CN 201710851734 A CN201710851734 A CN 201710851734A CN 107612436 B CN107612436 B CN 107612436B
- Authority
- CN
- China
- Prior art keywords
- rotor
- harmonic
- formula
- sine
- steps
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 21
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000010276 construction Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 11
- 230000004907 flux Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开一种基于永磁电机谐波反电势的转子位置计算方法。首先对5、7次谐波反电势进行标幺化处理;然后,通过对标幺化谐波反电势进行和差化积处理,得到含有转子位置正余弦作为因数的变量;接着,利用积化和差、三倍角公式构造解调参考波,对提取得到的谐波反电势进行变换,得到幅值恒定的转子位置的正弦与余弦值,从而解调出转子位置的正余弦值;最后,利用反正切运算,通过正弦与余弦值得到转子位置信息,从而得到转子位置。本发明可以规避基于基波反电势计算谐波转子位置计算方法中对电机参数敏感性,可以得到鲁棒性更强的无传感器算法,便于及时获悉转子实时所处位置并对转子位置进行调整,从而保证永磁电机处于高效的工作状态。
Description
技术领域
本发明公开一种转子位置计算方法,特别涉及一种基于永磁电机谐波反电势的转子位置的计算方法,属于永磁电机无传感器算法领域。
背景技术
近20年,永磁电机的无传感器算法得到了广泛的应用。在很多应用场合下,无传感器控制策略是必不可少的。有的场合受限于空间或重量的限制,不能安装旋转位置传感器。而诸如风电等大型永磁电机应用的场合,更是无法在其轴上安装旋转位置编码器。不过,无传感器控制技术应用的最大优势在于降低成本。在一般永磁电机的应用场合,旋转位置传感器往往是设备成本中比重很大的一部分,因此,省去此类传感器无疑可以大大提高永磁电机驱动***的经济性。
通过电机电压、电流及反电势等变量计算永磁电机转子位置的方法是永磁电机无位置传感器算法的基础。传统的无传感器方法均是基于基波反电势或基波磁链来计算转子位置角度。一般使用在固定坐标系下或者观测旋转坐标系下得到的转子磁链或者反电势的正交分量,通过反正切函数直接得到转子位置角度。但是由于这种方法基于基波模型,而电机为了正常工作必须有基波电流,所以不可避免地会受到电机定子参数,如定子电阻与电枢电感的影响。这就对转子位置提取的精度产生了负面影响。
发明内容
为了克服现有技术的不足,本发明提供了一种基于永磁电机谐波反电势的转子位置计算方法,可以不必通过基波反电势或基波磁链来获得转子位置信息。
一种基于永磁电机谐波反电势的转子位置计算方法,步骤如下:
(1)对谐波反电势进行标幺化处理;
(2)通过对标幺化谐波反电势进行和差化积得到含有转子位置正余弦作为因数的变量;
(3)构造解调参考波;
(4)从步骤(2)所得变量中,使用步骤(3)构造的解调参考波解调出幅值恒定的转子位置正余弦信息;
(5)利用转子位置正余弦信息得到转子位置。
所述步骤(1)包括如下步骤:
(1A)对如式(1)所示的α-β坐标系下的5次谐波反电势7次谐波反电势进行标幺化处理,得到如式(2)、式(3)的标幺制5、7次谐波反电势:
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小。
所述步骤(2)包括如下步骤:
(2A)使用和差化积公式,对5、7次谐波反电势,即式(2)、式(3)做如下运算:
所述步骤(3)包括如下步骤:
(3A)需构造更多与sin(6θr)、cos(6θr)相关的变量作为解调参考波,以解调出sin(6θr)与cos(6θr),对式(2)、式(3)使用积化和差公式,得:
(3B)联立步骤(3A)中得到的式(8)至式(11),得:
(3C)由三倍角公式:
将式(12)、式(13)分别代入式(14)、式(15)可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr) (16)
sin(6θr)=3sin(2θr)-4sin3(2θr) (17)。
所述步骤(4)包括如下步骤:
(4A)由式(16)、式(17)以及式(4)—(7)解调出转子位置的正余弦值
所述步骤(5)包括如下步骤:
(5A)使用反正切函数通过转子位置正余弦值得到转子位置:
本发明的有益效果:
本发明实现了一种基于谐波反电势的转子位置计算方法,可以规避基于基波反电势计算谐波转子位置计算方法中对电机参数敏感性。而对于谐波反电势的提取,则可以采用将对应次谐波抑制至0的方式,从而提高反电势提取的鲁棒性。此方法为永磁电机无传感器算法提供了新的转子位置计算方式,在与对参数鲁棒的谐波反电势提取方法相结合后可以得到鲁棒性更强的无传感器算法。
附图说明
图1为本发明基于永磁电机谐波反电势的转子位置计算方法的一种流程图。
图2为5、7次谐波反电势标幺化后在α轴方向的分量sin(5θr)、sin(7θr)的曲线。
图3为通过提出的计算方法计算求得的转子所处位置θr随时间t的变化曲线与设定的参考转子位置随时间t的变化曲线。
具体实施方式
无传感器算法中的转子位置信息可以不通过基波反电势来提取,而通过谐波反电势来提取。由于谐波反电势对于电机出力的作用很小,所以我们发现若通过将谐波电流抑制为0的方法提取谐波反电势,则获得的转子信息可以是对参数不敏感的。
下面结合附图和实例对本发明作进一步的阐述。
本发明基于永磁电机谐波反电势的转子位置计算方法的一种流程图如图1所示。
首先,为得到幅值相等的5次谐波反电势与7次谐波反电势,以进行后续的变换,需对谐波反电势进行标幺化处理。如图2所示为5、7次谐波反电势经过标幺化处理后在α轴方向的分量sin(5θr)、sin(7θr)。
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小。
考虑到5θr=6θr-θr、7θr=6θr+θr,因此使用和差化积公式能够得到与目标θr有关的条件,故对5、7次谐波反电势做如下运算:
此时得到的有关基波量的条件由于系数部分为变量sin(6θr)、cos(6θr)而非常数,所以无法直接得到θr。因此必须得到6θr的正弦与余弦的信息,从而将sinθr、cosθr的系数变换为常数。为得到6θr的正弦与余弦的信息,首先考虑得到2θr:
经过简单的消元计算可得:
为了得到sin(6θr)、cos(6θr)以确定前文已得的有关基波量的条件的系数部分将sin(2θr)、cos(2θr)代入三倍角公式:
可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr)
sin(6θr)=3sin(2θr)-4sin3(2θr)。
此时,sin(6θr)、cos(6θr)的信息已经被构造出来,故可以将sin(6θr)、cos(6θr)作为参考波与上文得到的条件相乘。构造如下:
在已知了sinθr与cosθr后,可方便地通过反正切变换得到转子位置θr的值:
为了能够更直观、快速地计算求得的转子所处位置θr,在线计算时应选择调用atan2函数而非普通的的反正切函数atan。这是因为atan2函数具有解出四象限角度的能力,也即其输出范围是-π至π之间,而非像atan一样在-π/2至π/2之间。
另外,由于此种方法涉及很多三角函数运算,因此在在线计算的时候需调用CPU的快速计算函数库,如TI公司DSP所使用的fastRTS_fbu32.h库,其三角函数与反三角函数的计算才用高精度的查表法,而非一般的math.h库在计算三角函数时采用的数值计算方式,因而可以极大地提高本算法的计算速度。
实测转子所处位置θr与参考转子位置如图3所示,由于实测转子所处位置θr处于-π至π之间,而参考转子位置处于0至2π之间,因此印证了通过所提方法计算求得的转子所处位置θr与参考转子位置是吻合的,并由此证明了所提方法的有效性。
Claims (1)
1.一种基于永磁电机谐波反电势的转子位置计算方法,其特征在于,步骤如下:
(1)对谐波反电势进行标幺化处理;
(2)通过对标幺化谐波反电势进行和差化积得到含有转子位置正余弦作为因数的变量;
(3)构造解调参考波;
(4)从步骤(2)所得变量中,使用步骤(3)构造的解调参考波解调出幅值恒定的转子位置正余弦信息;
(5)利用转子位置正余弦信息得到转子位置;
所述步骤(1)包括如下步骤:
(1A)对如式(1)所示的α-β坐标系下的5次谐波反电势7次谐波反电势进行标幺化处理,得到如式(2)、式(3)的标幺制5、7次谐波反电势:
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小;
所述步骤(2)包括如下步骤:
(2A)使用和差化积公式,对5、7次谐波反电势,即式(2)、式(3)做如下运算:
所述步骤(3)包括如下步骤:
(3A)需构造更多与sin(6θr)、cos(6θr)相关的变量作为解调参考波,以解调出sin(6θr)与cos(6θr),对式(2)、式(3)使用积化和差公式,得:
(3B)联立步骤(3A)中得到的式(8)至式(11),得:
(3C)由三倍角公式:
将式(12)、式(13)分别代入式(14)、式(15)可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr) (16)
sin(6θr)=3sin(2θr)-4sin3(2θr) (17);
所述步骤(4)包括如下步骤:
(4A)由式(16)、式(17)以及式(4)—(7)解调出转子位置的正余弦值
所述步骤(5)包括如下步骤:
(5A)使用反正切函数通过转子位置正余弦值得到转子位置:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851734.4A CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851734.4A CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107612436A CN107612436A (zh) | 2018-01-19 |
CN107612436B true CN107612436B (zh) | 2019-11-22 |
Family
ID=61061201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710851734.4A Expired - Fee Related CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107612436B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110380733A (zh) * | 2019-05-21 | 2019-10-25 | 日本电产凯宇汽车电器(江苏)有限公司 | 一种汽车刹车***控制器用旋转变压器解码方法 |
EP3754079B1 (en) * | 2019-06-18 | 2022-09-14 | Memminger-IRO GmbH | A yarn delivery device and a method for delivering yarn to a textile machine |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
CN106841901A (zh) * | 2017-03-09 | 2017-06-13 | 浙江大学 | 一种变频器驱动内嵌式永磁同步电机定子匝间短路故障诊断方法 |
CN106877768A (zh) * | 2017-02-20 | 2017-06-20 | 杭州富生电器有限公司 | 多相永磁电机转子位置辨识方法 |
CN107070318A (zh) * | 2017-05-09 | 2017-08-18 | 天津大学 | 一种永磁无刷直流电机换相转矩脉动的谐波分析法 |
US9739815B2 (en) * | 2012-07-20 | 2017-08-22 | Brose Fahrzeugteile Gmbh & Co. Kommanditsellschaft, Wurzburg | Method for determining the rotor position of an electrically-commuted multi-phase direct current motor |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101694216B1 (ko) * | 2012-12-10 | 2017-01-09 | 한국전자통신연구원 | 전동기 구동 모듈 및 브러쉬리스 전동기 시스템 |
-
2017
- 2017-09-19 CN CN201710851734.4A patent/CN107612436B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9739815B2 (en) * | 2012-07-20 | 2017-08-22 | Brose Fahrzeugteile Gmbh & Co. Kommanditsellschaft, Wurzburg | Method for determining the rotor position of an electrically-commuted multi-phase direct current motor |
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
CN106877768A (zh) * | 2017-02-20 | 2017-06-20 | 杭州富生电器有限公司 | 多相永磁电机转子位置辨识方法 |
CN106841901A (zh) * | 2017-03-09 | 2017-06-13 | 浙江大学 | 一种变频器驱动内嵌式永磁同步电机定子匝间短路故障诊断方法 |
CN107070318A (zh) * | 2017-05-09 | 2017-08-18 | 天津大学 | 一种永磁无刷直流电机换相转矩脉动的谐波分析法 |
Also Published As
Publication number | Publication date |
---|---|
CN107612436A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3571698B2 (ja) | モータの回転速度制御装置 | |
CN1301586C (zh) | 电动机的磁极位置检测方法及磁极位置检测装置和使用该装置的电动机控制装置 | |
CN101615825B (zh) | 一种永磁同步电机旋转变压器标定及调零***及其运作方法 | |
JP5015307B2 (ja) | 線形ホールセンサーを利用した直交座標形速度具現方法 | |
ES2810876T3 (es) | Procedimiento para la identificación de la anisotropía magnética de una máquina eléctrica de campo giratorio | |
CN108196213A (zh) | 一种旋转变压器的零位角测试装置、方法及*** | |
TW201722060A (zh) | 用於具有脈衝寬度調變驅動之多相位永磁馬達的線性霍爾效應感測器 | |
JPWO2013136612A1 (ja) | 角度検出装置 | |
JP5660381B2 (ja) | 回転角検出装置 | |
CN101561263A (zh) | 一种永磁同步电机转子位置检测方法及其使用的旋转变压器 | |
CN107769636B (zh) | 一种无位置传感器永磁同步电机转子位置检测方法 | |
CN107612436B (zh) | 基于永磁电机谐波反电势的转子位置计算方法 | |
CN104426449B (zh) | 电动机控制装置以及涡轮分子泵 | |
CN108900127A (zh) | 考虑交叉耦合效应的ipmsm低速段无位置传感器控制方法 | |
CN103078572A (zh) | 一种永磁同步电机高精度转子位置估计方法 | |
CN207780217U (zh) | 一种旋转变压器的零位角测试装置 | |
CN103701395A (zh) | 一种基于正反序列谐波注入的电机转子初位估计方法 | |
US20200373863A1 (en) | Slotless synchronous permanent magnet motor | |
CN107834914A (zh) | 基于霍尔的无刷电机转子角度计算方法 | |
CN111769779A (zh) | 基于改进型Luenberger观测器的PMSM直接转矩控制方法 | |
KR20180019920A (ko) | 전동기의 회전자 각도 추정 장치 및 방법 | |
CN108322121A (zh) | 一种无人机电机转子初始位置检测方法、装置、设备及存储介质 | |
Wang et al. | Position estimation of outer rotor PMSM using linear Hall effect sensors and neural networks | |
CN110086399B (zh) | 一种永磁同步电机转子位置复合检测及起动运行方法 | |
Reddy et al. | A low cost sense coil based position sensing system for SRM implemented in a SoC FPGA |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191122 |