CN111510043B - 基于倍角坐标系的旋转高频信号注入***及位置提取方法 - Google Patents
基于倍角坐标系的旋转高频信号注入***及位置提取方法 Download PDFInfo
- Publication number
- CN111510043B CN111510043B CN202010240784.0A CN202010240784A CN111510043B CN 111510043 B CN111510043 B CN 111510043B CN 202010240784 A CN202010240784 A CN 202010240784A CN 111510043 B CN111510043 B CN 111510043B
- Authority
- CN
- China
- Prior art keywords
- axis
- module
- current
- voltage
- frequency
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 63
- 238000002347 injection Methods 0.000 title claims abstract description 42
- 239000007924 injection Substances 0.000 title claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 82
- 238000006243 chemical reaction Methods 0.000 claims abstract description 63
- 238000004458 analytical method Methods 0.000 claims abstract description 50
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 230000001360 synchronised effect Effects 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims description 66
- 150000001875 compounds Chemical class 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 101100129590 Schizosaccharomyces pombe (strain 972 / ATCC 24843) mcp5 gene Proteins 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000003137 locomotive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
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/22—Current control, e.g. using a current control loop
-
- 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
-
- 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/24—Vector control not involving the use of rotor position or rotor speed sensors
-
- 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
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
- H02P25/022—Synchronous motors
- H02P25/024—Synchronous motors controlled by supply frequency
- H02P25/026—Synchronous motors controlled by supply frequency thereby detecting the rotor position
-
- 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
-
- 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
- H02P2205/00—Indexing scheme relating to controlling arrangements characterised by the control loops
- H02P2205/01—Current loop, i.e. comparison of the motor current with a current reference
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开了一种基于倍角坐标系的旋转高频信号注入***及位置提取方法,属于无位置传感器控制技术领域。本发明的注入***包括:电源电路、整流器、逆变器、永磁同步电机、电机负载模块、电机电流采集模块、电机电压采集模块、Clark变换模块、Park变换模块、倍角坐标系分析计算模块、响应电流信号分析模块、转子位置提取模块、速度/角度计算模块、转速误差模块、d轴电流误差模块、d轴电流PI模块、q轴电流误差模块、高频信号注入模块、Park逆变换模块等。本发明能实现电机在低速情况下的无位置运行,获得预估转子位置信号,在旋转高频信号注入法的基础上,降低了低通滤波器的使用,简化了信号提取过程,在响应速度以及信号提取精度上具有优越性。
Description
技术领域
本发明涉及一种基于倍角坐标系的旋转高频信号注入***及位置提取方法,属于无位置传感器控制技术领域。
背景技术
永磁同步电机以其体积小、效率高、功率密度高和易控制、使用周期长和可靠性强等特点在伺服驱动领域具有主流优势,在要求高控制精度和高可靠性的场合如舰船推进、数控机床、机车牵引、电动汽车和家用电器等许多领域获得极为广泛的应用,成为各国学者的研究热点。
在电动工具领域,其体积要求及工作环境较为复杂,而传统的位置式传感器难以满足其加工精度以及工作要求,所以无位置技术显得尤为重要。为了满足低速下转子位置跟踪稳定、带载能力强等条件,选用了旋转高频信号注入法作为方案。
旋转高频信号注入法由于需要额外注入高频信号,同时转子位置信号中夹杂直流分量,基频以及高频信号波,因此,如何准确提取所需信号频率显得尤为重要。传统的信号提取方法是通过在电机两相静止坐标系下分别注入频率相同的高频余弦信号和正弦信号,形成与注入信号频率一致的圆形旋转磁场。检测该激励信号产生的电流响应,通过设计特定的信号处理,最终获得转子位置与转速信息,实现无位置传感器控制。该方法引入带通、低通等滤波器,经过调制输出信号幅值。由于提取方法存在带通、低通滤波器,引入了幅值误差以及相位偏移,且造成了软件延时,***相对复杂等问题。同时,提取效果完全取决于滤波器的处理效果,因而提取效果不佳。
发明内容
为了克服现有技术的不足,本发明提出了一种基于倍角坐标系的旋转高频信号注入***及位置提取方法,通过建立倍角坐标系并通过简单三角函数计算来实现高频信号幅值的提取方法。
本发明为解决其技术问题采用如下技术方案:
一种基于倍角坐标系的旋转高频信号注入***,包括:电源电路、整流器、逆变器、永磁同步电机、电机负载模块、电机电流采集模块、电机电压采集模块、第一Clark变换模块、第二Clark变换模块、第一Park变换模块、第二Park变换模块、倍角坐标系分析计算模块、响应电流信号分析模块、转子位置提取模块、速度/角度计算模块、转速误差模块、转速PI模块、d轴电流误差模块、d轴电流PI模块、q轴电流误差模块、q轴电流PI模块、高频信号注入模块、Park逆变换模块、空间矢量脉宽调制模块和低通滤波器模块;其中:所述电源电路、整流器、逆变器、永磁同步电机和电机负载模块顺序连接,所述d轴电流误差模块、d轴电流PI模块、Park逆变换模块和空间矢量脉宽调制模块顺序连接,所述转速误差模块、转速PI模块、q轴电流误差模块、q轴电流PI模块和Park逆变换模块顺序连接,所述电机电流采集模块、第一Clark变换模块、第一Park变换模块、低通滤波器模块和q轴电流误差模块顺序连接,所述电机电压采集模块、第二Clark变换模块、倍角坐标系分析计算模块、响应电流信号分析模块和转子位置提取模块顺序连接,所述第二Park变换模块、转子位置提取模块、速度/角度计算模块顺序连接,所述速度/角度计算模块分别与转速误差模块、Park逆变换模块、第一Park变换模块、第二Park变换模块连接,所述高频信号分别与Park逆变换模块、第二Park变换模块连接;
所述电源电路为单相交流电源,用于为整流器提供单相交流电;
所述整流器为单相不控整流器,用于将单相输入交流电整流为直流电,并向逆变器供电;
所述逆变器为三相电压源型逆变器,用于接收空间矢量脉宽调制模块的电压脉冲,并根据电压脉冲控制永磁同步电机;
所述电机负载模块为外部负载,用于永磁同步电机加载/卸载;
所述电机电流采集模块用于采集永磁同步电机三相电流,并发送至第一Clark变换模块;
所述电机电压采集模块用于采集永磁同步电机三相电压,并发送至第二Clark变换模块;
所述第一Clark变换模块和第二Clark变换模块用于将永磁同步电机三相电流转换为α、β轴电流,电机三相电压转换为α、β轴电压,并发送至第一Park变换模块和倍角坐标系分析计算模块;
所述第一Park变换模块和第二Park变换模块用于将α、β轴电流转换为d轴估计电流、q轴估计电流并发送至低通滤波器模块,将轴高频注入电压转换为d轴实际高频电压、q轴实际高频电压并发送至转子位置提取模块;
所述低通滤波器模块用于将d轴估计电流、q轴估计电流滤除高频电流成分,得到d轴基波电流、q轴基波电流,并发送至d轴电流误差模块、q轴电流误差模块;
所述倍角坐标系分析计算模块用于将α、β轴电流和电压进行分析计算,得到倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流,并发送至响应电流信号分析模块;
所述响应电流信号分析模块用于将γ、δ轴高频响应电流与α、β轴高频响应电流进行分析与计算,得到包含转子位置信息的混合信号,并发送至转子位置提取模块;
所述转子位置提取模块用于将包含转子位置信息的混合信号处理后得到转子位置正余弦信号,并发送至速度/角度计算模块;
所述速度/角度计算模块用于将转子位置正余弦信息通过相应计算,得到预估转子转速、角度及角度误差信息,并发送至转速误差模块、Park逆变换模块、第一Park变换模块和第二Park变换模块;
所述转速误差模块用于将给定转速与所述速度/角度计算模块计算得到的预估转速对比,得到转速误差,并将其发送至转速PI模块;
所述转速PI模块用于将转速误差进行PI调节,得到q轴给定电流,并发送至q轴电流误差模块;
所述d轴电流误差模块用于将d轴给定电流与d轴基波电流对比,得到d轴电流误差,并将其发送至d轴电流PI模块;
所述q轴电流误差模块用于将q轴给定电流与q轴基波电流对比,得到q轴电流误差,并将其发送至q轴电流PI模块;
所述d轴电流PI模块用于将d轴电流误差进行PI调节,得到d轴估计电压,并发送至电压Park逆变换模块;
所述q轴电流PI模块用于将q轴电流误差进行PI调节,得到q轴估计电压,并发送至电压Park逆变换模块;
所述高频信号注入模块用于给电机预估d、q轴注入高频信号,发送到Park逆变换模块及第二Park变换模块;
所述Park逆变换模块用于将d、q轴电压转换为α轴电压、β轴电压,并发送至空间矢量脉宽调制模块;
所述空间矢量脉宽调制模块用于根据α、β轴电压、母线电压计算得到电压脉冲,并发送至逆变器。
基于倍角坐标系的旋转高频信号注入***的位置提取方法,包括以下步骤:
(1)实时采集单相交流输入电压、直流母线电压幅值及相位,实时采集永磁同步电机三相电流及三相电压,并利用信号提取方法来提取转子位置;
(2)对所述的永磁同步电机abc三相电流和电压进行Clark变换,得到α、β轴电流和电压,对所述α、β轴电流进行Park变换,得到d、q轴估计电流,对所述轴高频注入电压进行Park变换,得到d、q轴实际高频电压;
(3)对所述的α、β轴电流和电压按照倍角坐标系信号分析方法进行计算,得到倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流;
(4)对所述的γ、δ轴高频响应电流与α、β轴高频响应电流按照响应电流信号分析方法进行计算,得到包含转子位置信息的混合信号;
(5)对所述的包含转子位置信息的混合信号按照转子位置提取原理进行转子位置提取,得到转子位置角正余弦信号,通过相应计算获得预估转速;
(6)计算给定电机转速与所述预估转速的误差,并对转速误差进行PI调节;得到q轴给定电流;
(7)计算d、q轴给定电流与d、q轴基波电流误差,并对电流误差进行PI调节,得到d、q轴电压估计值;
(8)给电机预估d、q轴注入高频信号,与d、q轴电压估计值叠加进行Park逆变换,得到α、β轴电压;
(9)根据α、β轴电压、母线电压,对逆变器进行SVPWM调制,并通过逆变器控制电机。
步骤(1)所述信号提取的方法包含如下步骤:
三相电流通过第一Clark变换模块得到α轴实际电流、β轴实际电流,三相电压通过第二Clark变换模块得到α轴实际电压、β轴实际电压,通过倍角坐标系分析计算模块求取倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流,并通过响应电流信号分析模块进一步求得包含转子位置信息的混合信号,再通过转子位置提取模块求出转子位置信息角的正余弦信号,进一步通过速度/角度计算模块获得预估转子转速、角度及角度误差信息。
步骤(3)所述倍角坐标系定义如下:
通过将α和β轴逆时针旋转2θe来建立新的测量旋转坐标系,记为γ、δ坐标系,其中γ轴超前同步旋转坐标轴d轴θe度,δ轴超前同步旋转坐标轴q轴θe度;假设α、β轴高频响应电流为iαh、iβh,γ、δ轴高频响应电流为iγh、iδh,则两者关系表示为式(1):
式中,θe为磁极位置电角度,即同步旋转坐标d、q轴与静止坐标轴α、β轴之间的角度。
步骤(3)所述倍角坐标系分析方法过程如下:
注入的高频电压信号与经过d、q轴电流PI模块的d、q轴电压相叠加,经过Park逆变换以及SVPWM调制,并通过逆变器控制电机;
通过电机电压电流采集模块得到三相电压和电流,并进行Clark变换,得到α、β轴电压和电流,对所述α、β轴电压和电流通过带通滤波器滤波得到高频电压信号uαh、uβh和高频电流信号iαh、iβh;
计算永磁同步电机定子电压的基本方程,如式(3)表示:
式中,uα、uβ为α、β轴定子电压,iα、iβ为α、β轴定子电流,Rs为定子绕组电阻,ωe为电机电角速度,λf为转子励磁磁链,Ls为电感矩阵,如式(4)表示:
计算理想低速***下的注入高频电压信号在α、β轴产生的电压如式(5)所示:
对上式进行积分和微分变换并展开,得到等式(6):
进一步的,通过对高频电压积分,结合高频电流计算来提取高频注入电压在倍角坐标系下的高频电流响应iγh、iδh,如式(7):
步骤(4)所述响应电流分析方法过程如下:
将iγh、iδh进行Park逆变换,得到同步旋转d、q轴响应电流,同时将iαh、iβh进行Park变化,也得到同步旋转d、q轴响应电流;理想情况下,变换后的电流相等,记为式(8):
将上式展开,得到以下方程:
两边同时除以cosθe,得到以下方程:
因此角度θe的正切值由下式得到:
其中:inum、iden代表求取第一种转子正切值的分子与分母,i'num、i'den代表求取第二种转子正切值的分子与分母;
通过结合倍角参考坐标系,求取相应倍角坐标系响应电流,实现对转子位置角进行提取;
进一步的,对计算转子位置正切值的分子分母进行分析计算,得到下列表达式:
分析分子分母信号组成形式及波形,发现电流信号中含有基波、高频信号波,且为乘积形式表示;在实际信号中,通过将inum、iden及i'num、i'den送入转子位置提取模块,结合Park变换后的实际d、q轴高频电压udh、uqh即求取转子位置角度的正余弦信息。
步骤(5)所述转子位置提取过程如下:
对inum、iden及i'num、i'den进行分析,并进一步去除幅值系数后,重新记为inum1、iden1及i'num1、i'den1,则由下式表示:
同时,对注入高频电压信号进行Park变换,得到同步旋转d、q轴系下的高频电压:
同样对其进一步去除幅值系数后,重新记为udh1、uqh1,则由下式表示:
根据简单三角函数数学关系进行角度正余弦信号提取,方法为:
本发明的有益效果如下:
1、通过倍角坐标系对高频响应电流进行提取,在简单数学计算后,简化位置信号公式表达形式,使用转子位置提取模块提取转子位置的正余弦信号,有效避免了高频项对解调的耦合效果,简化了***。
2、相较于以往的滤波器调制的信号提取方法,减少了滤波器的使用,同时信号提取更快,***更加简便,且不需要引入额外的解调信号。因此,以提取转子位置信号更加优越为前提,实现电机旋转高频信号注入法在低速下带载运行,优化了转子位置估计的准确性和鲁棒性,增强无位置控制***的稳定性及实用性,控制简单有效。
3、简化信号提取环节,根据理论推导得到包含转子位置角的混合信号形式,在固定注入频率下提取转子位置角正余弦信息,达到准确、快速的信号提取目的,同时增强***鲁棒性及实用性。
4、去除解调***所需的低通滤波器的设计过程,从而避免了低通滤波器引入的延时、相位偏移及幅值损耗等问题,避免了调节的复杂性,简化了调节过程,有效增加了信号提取的快速性和准确性。另一方面,准确和快速的提取信号使得信号中含有的转子位置信号更加精确,转子位置估计更加精确、快速,误差减小。
附图说明
图1为本发明一种基于倍角坐标系的旋转高频信号注入法信号提取***结构框图。
图2为本发明倍角参考坐标系示意图。
图3为本发明倍角坐标系分析模块结构框图。
图4为本发明响应电流分析模块结构框图。
图5为本发明转子位置提取模块结构框图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例仅用于解释本发明,并非用于限定本发明的范围。
本发明公开了一种基于倍角坐标系的旋转高频信号注入***及位置提取方法,图1为本发明的***结构框图,包括:电源电路、整流器、逆变器、永磁同步电机、电机负载模块、电机电流采集模块、电机电压采集模块、Clark(abc轴-αβ轴)变换模块、Park(αβ轴-dq轴)变换模块、低通滤波器模块、倍角坐标系分析计算模块、响应电流信号分析模块、转子位置提取模块、速度/角度计算模块、转速误差模块、转速PI(比例积分)模块、d轴电流误差模块、d轴电流PI(比例积分)模块、q轴电流误差模块、q轴电流PI(比例积分)模块、高频信号注入模块、Park逆变换(dq轴-αβ轴)模块、脉宽调制模块、和低通滤波器模块。
所述电源电路为单相交流电源,用于为整流器提供单相交流电。
所述整流器为单相不控整流器,用于将单相输入交流电整流为直流电,并向逆变器供电。
所述逆变器为三相电压源型逆变器,用于接收脉宽调制模块的电压脉冲,并根据电压脉冲控制电机。
所述电机负载模块为外部负载,用于电机加载/卸载。
所述电机电流采集模块用于采集电机三相电流,并发送至第一Clark变换模块。
所述电机电压采集模块用于采集电机三相电压,并发送至第二Clark变换模块。
所述第一Clark变换模块和第二Clark变换模块用于将电机三相电流转换为α、β轴电流,电机三相电压转换为α、β轴电压,并发送至倍角坐标系分析计算模块和第一Park变换模块。
所述第一Park变换模块和第二Park变换模块用于将α、β轴电流转换为d轴估计电流、q轴估计电流并发送至低通滤波器模块,将轴高频注入电压转换为d轴实际高频电压、q轴实际高频电压并发送至转子位置提取模块。
所述低通滤波器模块用于将d轴估计电流、q轴估计电流滤除高频电流成分,得到d轴基波电流、q轴基波电流,并发送至d轴电流误差模块、q轴电流误差模块。
所述倍角坐标系分析计算模块用于将α、β轴电流和电压进行分析计算,得到倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流,并发送至响应电流信号分析模块。
所述响应电流信号分析模块用于将γ、δ轴高频响应电流与α、β轴高频响应电流进行分析与计算,得到包含转子位置信息的混合信号,并发送至转子位置提取模块。
所述转子位置提取模块用于将包含转子位置信息的混合信号处理后得到转子位置正余弦信号,并发送至速度/角度计算模块。
所述速度/角度计算模块用于将转子位置正余弦信息通过相应计算,得到预估转子转速、角度及角度误差信息,并发送至转速误差模块、Park逆变换模块、第一Park变换模块和第二Park变换模块。
所述转速误差模块用于将给定转速与所述速度/角度计算模块计算得到的预估转速对比,得到转速误差,并将其发送至转速PI模块。
所述转速PI模块用于将转速误差进行PI调节,得到q轴给定电流,并发送至q轴电流误差模块。
所述d轴电流误差模块用于将d轴给定电流与d轴基波电流对比,得到d轴电流误差,并将其发送至d轴电流PI模块。
所述q轴电流误差模块用于将q轴给定电流与q轴基波电流对比,得到q轴电流误差,并将其发送至q轴电流PI模块。
所述d轴电流PI模块用于将d轴电流误差进行PI调节,得到d轴估计电压,并发送至电压Park逆变换模块。
所述q轴电流PI模块用于将q轴电流误差进行PI调节,得到q轴估计电压,并发送至电压Park逆变换模块。
所述高频信号注入模块用于给电机预估d、q轴注入高频信号,发送到Park逆变换模块及第二Park变换模块。
所述Park逆变换模块用于将d、q轴电压转换为α轴电压、β轴电压,并发送至脉宽调制模块。
所述脉宽调制模块为空间矢量脉宽调制,用于根据α、β轴电压、母线电压计算得到电压脉冲,并发送至逆变器。
一种基于倍角坐标系的旋转高频信号注入***的位置提取方法,包括以下步骤:在预估d、q轴上注入高频信号,引起***的电流电压响应,实时采集电机三相电流Ia、Ib、Ic;对电机abc三相电流进行Clark变换,得到α、β轴电流iα、iβ;实时采集电机三相电压Ua、Ub、Uc;对电机abc三相电压进行Clark变换,得到α、β轴电压uα、uβ;对α、β轴电流进行Park变换,得到d、q轴估计电流通过低通滤波器得到d、q基波电流id_f、iq_f;通过倍角坐标系分析计算模块求取倍角坐标系γ、δ轴高频响应电流iγh、iδh及α、β轴高频响应电流iαh、iβh;并通过响应电流分析模块进一步求得包含转子位置信息的混合信号inum、iden及i'num、i'den;对轴高频注入电压经Park变换,得到d、q轴实际高频电压udh、uqh;结合inum、iden及i'num、i'den,通过转子位置提取模块进一步求取转子位置正余弦信号sinθe、cosθe,最后通过速度/角度计算模块得到预估转子转速角度及角度误差Δθe,计算给定电机转速n*与所述预估电机转速的误差,并对转速误差进行PI调节,得到q轴给定电流iq*;计算d、q轴给定电流与估计d、q轴电流误差,并对电流误差进行PI调节,得到d、q轴电压估计值与轴高频注入电压叠加得到d、q轴计算电压ud*、uq*;对d、q轴计算电压进行Park逆变换,得到α、β轴电压uα、uβ;根据αβ轴电压、母线电压,对逆变器进行SVPWM调制(空间矢量脉宽调制),并通过逆变器控制电机。所示倍角坐标系信号分析模块,包括倍角坐标系分析计算模块、响应电流分析模块、转子位置提取模块。
如图2所示为倍角参考坐标系示意图。
对倍角参考坐标系的定义说明如下:
为方便分析,通过将α和β轴逆时针旋转2θe来建立新的测量旋转坐标系,记为γ、δ坐标系,其中γ轴超前同步旋转坐标轴d轴θe度,δ轴超前同步旋转坐标轴q轴θe度。
假设α、β轴高频响应电流为iαh、iβh,γ、δ轴高频响应电流为iγh、iδh,则两者关系可以表示为式(1):
式中,θe为磁极位置电角度,即同步旋转坐标d、q轴与静止坐标轴α、β轴之间的角度。
如图3所示为倍角坐标系分析模块结构框图。
图中,iγh、iδh为所求倍角坐标系高频响应电流,iα、iβ为三相电流经Clark变换后的α、β轴电流,uα、uβ为三相电压经Clark变换后的α、β轴电压,经带通滤波器后得到α、β轴高频响应电流和电压iαh、iβh和uαh、uβh。对uαh、uβh积分后进行适当变换,同时结合iαh、iβh可得到倍角坐标系高频响应电流信号iγh、iδh。具体原理如下:
注入的高频电压信号与经过d、q轴电流PI模块的d、q轴电压相叠加,经过Park逆变换以及SVPWM调制,并通过逆变器控制电机;
通过电机电压电流采集模块得到三相电压和电流,并进行Clark变换,得到α、β轴电压和电流,对所述α、β轴电压和电流通过带通滤波器滤波得到高频电压信号uαh、uβh和高频电流信号iαh、iβh;
计算永磁同步电机定子电压的基本方程,如式(3)表示:
式中,uα、uβ为α、β轴定子电压,iα、iβ为α、β轴定子电流,Rs为定子绕组电阻,ωe为电机电角速度,λf为转子励磁磁链,Ls为电感矩阵,如式(4)表示:
计算理想低速***下的注入高频电压信号在α、β轴产生的电压如式(5)所示:
其中:uαh、uβh分别为注入高频电压信号在α、β轴产生的电压响应信号;
对上式进行积分和微分变换并展开,可得到等式(6):
进一步的,可以通过对高频电压积分,结合高频电流计算来提取高频注入电压在倍角坐标系下的高频电流响应iγh、iδh,如式(7):
如图4所示为响应电流分析模块结构框图。
图中,inum、iden分别为计算转子位置角度正切值的分子及分母部分,i'num、i'den分别为等价的计算转子位置角度正切值的分子及分母部分,且均由iαh、iβh及iγh、iδh进行简单计算得出,具体计算原理如下:
对于旋转高频电压注入引起的α、β轴高频响应电流iαh、iβh,会在γ、δ轴得到高频响应电流iγh、iδh。将iγh、iδh进行Park逆变换,得到同步旋转d、q轴响应电流,同时将iαh、iβh进行Park变化,也得到同步旋转d、q轴响应电流。理想情况下,变换后的电流相等,记为式(8):
将上式展开,可得到以下方程:
两边同时除以cosθe,可得到以下方程:
因此角度θe的正切值可以由下式得到:
可知,通过结合倍角参考坐标系,求取相应倍角坐标系响应电流,可以实现对转子位置角进行提取。
进一步的,对计算转子位置正切值的分子分母进行分析计算,得到下列表达式:
分析分子分母信号组成形式及波形,发现电流信号中含有基波、高频信号波,且为乘积形式表示。因此在实际信号中,通过将inum、iden及i'num、i'den送入转子位置提取模块,结合Park变换后的实际d、q轴高频电压udh、uqh即可求取转子位置角度的正余弦信息。
如图5所示为转子位置提取模块整体结构框图。
图中,对inum、iden及i'num、i'den进行分析,并进一步去除幅值系数后,重新记为inum1、iden1及i'num1、i'den1,则可由下式表示:
同时,对注入高频电压信号进行Park变换,可以得到同步旋转d、q轴系下的高频电压:
同样对其进一步去除幅值系数后,重新记为udh1、uqh1,则可由下式表示:
根据简单三角函数数学关系进行角度正余弦信号提取,方法为:
Claims (6)
1.一种基于倍角坐标系的旋转高频信号注入***,其特征在于,包括:电源电路、整流器、逆变器、永磁同步电机、电机负载模块、电机电流采集模块、电机电压采集模块、第一Clark变换模块、第二Clark变换模块、第一Park变换模块、第二Park变换模块、倍角坐标系分析计算模块、响应电流信号分析模块、转子位置提取模块、速度与角度计算模块、转速误差模块、转速PI模块、d轴电流误差模块、d轴电流PI模块、q轴电流误差模块、q轴电流PI模块、高频信号注入模块、Park逆变换模块、空间矢量脉宽调制模块和低通滤波器模块;其中:所述电源电路、整流器、逆变器、永磁同步电机和电机负载模块顺序连接,所述d轴电流误差模块、d轴电流PI模块、Park逆变换模块和空间矢量脉宽调制模块顺序连接,所述转速误差模块、转速PI模块、q轴电流误差模块、q轴电流PI模块和Park逆变换模块顺序连接,所述电机电流采集模块、第一Clark变换模块、第一Park变换模块、低通滤波器模块和q轴电流误差模块顺序连接,所述电机电压采集模块、第二Clark变换模块、倍角坐标系分析计算模块、响应电流信号分析模块和转子位置提取模块顺序连接,所述第二Park变换模块、转子位置提取模块、速度与角度计算模块顺序连接,所述速度与角度计算模块分别与转速误差模块、Park逆变换模块、第一Park变换模块、第二Park变换模块连接,所述高频信号注入模块分别与Park逆变换模块、第二Park变换模块连接;
所述电源电路为单相交流电源,用于为整流器提供单相交流电;
所述整流器为单相不控整流器,用于将单相输入交流电整流为直流电,并向逆变器供电;
所述逆变器为三相电压源型逆变器,用于接收空间矢量脉宽调制模块的电压脉冲,并根据电压脉冲控制永磁同步电机;
所述电机负载模块为外部负载,用于永磁同步电机加载/卸载;
所述电机电流采集模块用于采集永磁同步电机三相电流,并发送至第一Clark变换模块;
所述电机电压采集模块用于采集永磁同步电机三相电压,并发送至第二Clark变换模块;
所述第一Clark变换模块用于将永磁同步电机三相电流转换为α、β轴电流,并发送至第一Park变换模块和倍角坐标系分析计算模块;所述第二Clark变换模块用于将电机三相电压转换为α、β轴电压,并发送至倍角坐标系分析计算模块;
所述第一Park变换模块用于将α、β轴电流转换为d轴估计电流、q轴估计电流并发送至低通滤波器模块;所述第二Park变换模块用于将预估d、q轴高频注入电压转换为d轴实际高频电压、q轴实际高频电压并发送至转子位置提取模块;
所述低通滤波器模块用于将d轴估计电流、q轴估计电流滤除高频电流成分,得到d轴基波电流、q轴基波电流,并发送至d轴电流误差模块、q轴电流误差模块;
所述倍角坐标系分析计算模块用于将α、β轴电流和电压进行分析计算,得到倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流,并发送至响应电流信号分析模块;所述倍角坐标系定义如下:
通过将α和β轴逆时针旋转2θe来建立新的测量旋转坐标系,记为γ、δ坐标系,即为倍角坐标系,其中γ轴超前同步旋转坐标轴d轴θe度,δ轴超前同步旋转坐标轴q轴θe度;
假设α、β轴高频响应电流为iαh、iβh,γ、δ轴高频响应电流为iγh、iδh,则两者关系表示为式(1):
式中,θe为磁极位置电角度,即同步旋转坐标d、q轴与静止坐标轴α、β轴之间的角度;
所述响应电流信号分析模块用于将γ、δ轴高频响应电流与α、β轴高频响应电流进行分析与计算,得到包含转子位置信息的混合信号,并发送至转子位置提取模块;
所述转子位置提取模块用于将包含转子位置信息的混合信号处理后得到转子位置正余弦信号,并发送至速度与角度计算模块;
所述速度与角度计算模块用于将转子位置正余弦信息通过相应计算,得到预估转子转速、角度及角度误差信息,并将预估转子转速发送至转速误差模块,将角度发送至Park逆变换模块和第一Park变换模块,将角度误差发送至第二Park变换模块;
所述转速误差模块用于将给定转速与所述速度与角度计算模块计算得到的预估转子转速对比,得到转速误差,并将其发送至转速PI模块;
所述转速PI模块用于将转速误差进行PI调节,得到q轴给定电流,并发送至q轴电流误差模块;
所述d轴电流误差模块用于将d轴给定电流与d轴基波电流对比,得到d轴电流误差,并将其发送至d轴电流PI模块;
所述q轴电流误差模块用于将q轴给定电流与q轴基波电流对比,得到q轴电流误差,并将其发送至q轴电流PI模块;
所述d轴电流PI模块用于将d轴电流误差进行PI调节,得到d轴估计电压,并发送至Park逆变换模块;
所述q轴电流PI模块用于将q轴电流误差进行PI调节,得到q轴估计电压,并发送至Park逆变换模块;
所述高频信号注入模块用于给电机预估d、q轴注入高频信号,发送到Park逆变换模块及第二Park变换模块;
所述Park逆变换模块用于将d、q轴电压转换为α轴电压、β轴电压,并发送至空间矢量脉宽调制模块;
所述空间矢量脉宽调制模块用于根据α、β轴电压、母线电压计算得到电压脉冲,并发送至逆变器。
2.根据权利要求1所述的基于倍角坐标系的旋转高频信号注入***的位置提取方法,其特征在于,包括以下步骤:
(1)实时采集单相交流输入电压、直流母线电压幅值及相位,实时采集永磁同步电机三相电流及三相电压,并利用信号提取方法来提取转子位置;
(2)对所述的永磁同步电机abc三相电流和电压进行Clark变换,得到α、β轴电流和电压,对所述α、β轴电流进行Park变换,得到d、q轴估计电流,对所述预估d、q轴高频注入电压进行Park变换,得到d、q轴实际高频电压;
(3)对所述的α、β轴电流和电压按照倍角坐标系信号分析方法进行计算,得到倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流;
(4)对所述的γ、δ轴高频响应电流与α、β轴高频响应电流按照响应电流信号分析方法进行计算,得到包含转子位置信息的混合信号;
(5)对所述的包含转子位置信息的混合信号按照转子位置提取原理进行转子位置提取,得到转子位置角正余弦信号,通过相应计算获得预估转子转速;
(6)计算给定电机转速与所述预估转子转速的误差,并对转速误差进行PI调节,得到q轴给定电流;
(7)计算d、q轴给定电流与d、q轴基波电流误差,并对电流误差进行PI调节,得到d、q轴电压估计值;
(8)给电机预估d、q轴注入高频信号,与d、q轴电压估计值叠加进行Park逆变换,得到α、β轴电压;
(9)根据α、β轴电压、母线电压,对逆变器进行SVPWM调制,并通过逆变器控制电机。
3.根据权利要求2所述的基于倍角坐标系的旋转高频信号注入***的位置提取方法,其特征在于,步骤(1)所述信号提取的方法包含如下步骤:
三相电流通过第一Clark变换模块得到α轴实际电流、β轴实际电流,三相电压通过第二Clark变换模块得到α轴实际电压、β轴实际电压,通过倍角坐标系分析计算模块求取倍角坐标系γ、δ轴高频响应电流及α、β轴高频响应电流,并通过响应电流信号分析模块进一步求得包含转子位置信息的混合信号,再通过转子位置提取模块求出转子位置信息角的正余弦信号,进一步通过速度与角度计算模块获得预估转子转速、角度及角度误差信息。
4.根据权利要求2所述的基于倍角坐标系的旋转高频信号注入***的位置提取方法,其特征在于,步骤(3)所述倍角坐标系信号分析方法过程如下:
假设***在预估d、q轴上注入高频电压信号,如式(2):
注入的高频电压信号与经过d、q轴电流PI模块的d、q轴电压估计值相叠加,经过Park逆变换以及SVPWM调制,并通过逆变器控制电机;
通过电机电压电流采集模块得到三相电压和电流,并进行Clark变换,得到α、β轴电压和电流,对所述α、β轴电压和电流通过带通滤波器滤波得到高频电压信号uαh、uβh和高频电流信号iαh、iβh;
计算永磁同步电机定子电压的基本方程,如式(3)表示:
式中,uα、uβ为α、β轴定子电压,iα、iβ为α、β轴定子电流,Rs为定子绕组电阻,ωe为电机电角速度,λf为转子励磁磁链,Ls为电感矩阵,如式(4)表示:
对上式进行积分和微分变换并展开,得到等式(6):
进一步的,通过对高频电压积分,结合高频电流计算来提取高频注入电压在倍角坐标系下的高频电流响应iγh、iδh,如式(7):
5.根据权利要求4所述的基于倍角坐标系的旋转高频信号注入***的位置提取方法,其特征在于,步骤(4)所述响应电流信号分析方法过程如下:
将iγh、iδh进行Park逆变换,得到同步旋转d、q轴响应电流,同时将iαh、iβh进行Park变换,也得到同步旋转d、q轴响应电流;理想情况下,变换后的电流相等,记为式(8):
将上式展开,得到以下方程:
两边同时除以cosθe,得到以下方程:
因此角度θe的正切值由下式得到:
其中:inum、iden代表求取第一种转子正切值的分子与分母,i'num、i'den代表求取第二种转子正切值的分子与分母;
通过结合倍角坐标系,求取相应倍角坐标系响应电流,实现对转子位置角进行提取;
进一步的,对计算转子位置正切值的分子分母进行分析计算,得到下列表达式:
分析分子分母信号组成形式及波形,发现电流信号中含有基波、高频信号波,且为乘积形式表示;在实际信号中,通过将inum、iden及i'num、i'den送入转子位置提取模块,结合Park变换后的d、q轴实际高频电压udh、uqh求取转子位置角度的正余弦信息。
6.根据权利要求5所述的基于倍角坐标系的旋转高频信号注入***的位置提取方法,其特征在于,步骤(5)所述转子位置提取过程如下:
对inum、iden及i'num、i'den进行分析,并进一步去除幅值系数后,重新记为inum1、iden1及i'num1、i'den1,则由下式表示:
同时,对注入高频电压信号进行Park变换,得到同步旋转d、q轴系下的高频电压:
进一步去除幅值系数后,重新记为udh1、uqh1,则由下式表示:
根据三角函数数学关系进行角度正余弦信号提取,方法为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010240784.0A CN111510043B (zh) | 2020-03-31 | 2020-03-31 | 基于倍角坐标系的旋转高频信号注入***及位置提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010240784.0A CN111510043B (zh) | 2020-03-31 | 2020-03-31 | 基于倍角坐标系的旋转高频信号注入***及位置提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111510043A CN111510043A (zh) | 2020-08-07 |
CN111510043B true CN111510043B (zh) | 2022-04-22 |
Family
ID=71863897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010240784.0A Active CN111510043B (zh) | 2020-03-31 | 2020-03-31 | 基于倍角坐标系的旋转高频信号注入***及位置提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111510043B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113708691B (zh) * | 2021-09-01 | 2023-12-22 | 广东汇天航空航天科技有限公司 | 转子运行数据估测方法、计算设备、电机控制方法及*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110474587A (zh) * | 2019-09-03 | 2019-11-19 | 哈尔滨理工大学 | 无源控制高频信号注入的永磁同步电机无位置传感器控制***及控制方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006230169A (ja) * | 2005-02-21 | 2006-08-31 | Toshiba Corp | 同期機の制御装置 |
JP5561516B2 (ja) * | 2009-07-06 | 2014-07-30 | 株式会社ジェイテクト | モータ制御装置および車両用操舵装置 |
JP5652664B2 (ja) * | 2011-10-21 | 2015-01-14 | アイシン・エィ・ダブリュ株式会社 | 回転電機制御装置 |
CN103929109B (zh) * | 2014-05-12 | 2016-06-29 | 哈尔滨工业大学 | 一种大功率内置式永磁同步电机无位置传感器控制***及控制方法 |
JP6390337B2 (ja) * | 2014-10-21 | 2018-09-19 | 株式会社デンソー | 回転電機の制御装置 |
CN109495047B (zh) * | 2018-12-28 | 2021-06-04 | 东北大学 | 一种基于高频信号注入的永磁同步电机无传感器控制方法 |
CN109889117B (zh) * | 2019-04-04 | 2020-09-01 | 合肥工业大学 | 基于旋转高频注入法的ipmsm位置观测方法、***及驱动*** |
-
2020
- 2020-03-31 CN CN202010240784.0A patent/CN111510043B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110474587A (zh) * | 2019-09-03 | 2019-11-19 | 哈尔滨理工大学 | 无源控制高频信号注入的永磁同步电机无位置传感器控制***及控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111510043A (zh) | 2020-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101964624B (zh) | 永磁同步电机的无传感器控制*** | |
CN108288936B (zh) | 一种永磁直线电机低速无位置传感器控制方法 | |
CN110022106B (zh) | 一种基于高频信号注入的永磁同步电机无位置传感器控制方法 | |
CN107241042B (zh) | 基于并联epll的脉振高频信号注入法信号提取***及策略 | |
CN108900127B (zh) | 考虑交叉耦合效应的ipmsm低速段无位置传感器控制方法 | |
CN107508521B (zh) | 一种永磁同步电机的无速度传感器控制方法和*** | |
CN114598206B (zh) | 永磁同步电机宽速域转子位置观测器设计方法 | |
CN111245330A (zh) | 基于过采样的永磁电机方波注入无位置传感器控制方法 | |
CN109039204B (zh) | 车用永磁同步电机无位置传感器模型预测控制***及方法 | |
CN109412484B (zh) | 一种三相感应电机无传感器矢量控制的转子磁链估算方法 | |
CN106655952A (zh) | 一种检测永磁同步电机转子初始位置的电流包络线法 | |
CN106026831B (zh) | 一种无位置传感器控制下的电流解耦方法 | |
CN108512475B (zh) | 一种基于脉振高频电压注入法的转子位置检测方法 | |
CN111371365A (zh) | 基于三次谐波空间高频方波信号注入的五相ipmsm无位置传感器控制方法 | |
CN111510043B (zh) | 基于倍角坐标系的旋转高频信号注入***及位置提取方法 | |
CN113489409B (zh) | 基于旋转综合矢量的永磁电机高频注入位置估计方法 | |
CN113972874B (zh) | 一种永磁同步电机无位置传感器控制的位置误差补偿方法 | |
CN109546904B (zh) | 双三相永磁同步电机的转子位置检测方法 | |
US11863107B2 (en) | Device and method for controlling rotary electric machine | |
CN110649847A (zh) | 一种pmslm在低速段的无位置传感器控制方法 | |
CN114337402B (zh) | 一种同步磁阻电机电感参数辨识及转子位置混合估计方法 | |
CN113676106B (zh) | 用于六相永磁电机的双套绕组高频注入的无位置方法 | |
CN113381655B (zh) | 消除单电流传感器控制的电机***中采样延迟误差的方法 | |
CN113489410B (zh) | 周期互补高频方波注入的无传感器控制方法 | |
CN114726273A (zh) | 一种基于改进滑模观测器的pcb盘式永磁同步电机无感控制***及方法 |
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 |