CN111464099B - 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法 - Google Patents

单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法 Download PDF

Info

Publication number
CN111464099B
CN111464099B CN202010406190.2A CN202010406190A CN111464099B CN 111464099 B CN111464099 B CN 111464099B CN 202010406190 A CN202010406190 A CN 202010406190A CN 111464099 B CN111464099 B CN 111464099B
Authority
CN
China
Prior art keywords
plane
torque
suspension
voltage vector
current
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
Application number
CN202010406190.2A
Other languages
English (en)
Other versions
CN111464099A (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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN202010406190.2A priority Critical patent/CN111464099B/zh
Publication of CN111464099A publication Critical patent/CN111464099A/zh
Application granted granted Critical
Publication of CN111464099B publication Critical patent/CN111464099B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/141Flux estimation
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P27/00Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04Arrangements 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/06Arrangements 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/08Arrangements 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/12Arrangements 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 pulsing by guiding the flux vector, current vector or voltage vector on a circle or a closed curve, e.g. for direct torque control

Abstract

本发明提出一种单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,利用空间矢量调制方法,在双机电能量转换平面上输出满足要求的电压矢量,降低转矩及悬浮脉动。其分别在转矩控制平面和悬浮力控制平面实现二者的解耦控制,有效地增强了转子悬浮运行的平稳性;利用逆变器电压矢量直接实现电磁转矩和悬浮力的控制,加速了电磁转矩和悬浮力控制的动态响应;采用空间矢量调制方式,有效的降低了电磁转矩和悬浮力脉动。

Description

单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法
技术领域
本发明属于无轴承磁通切换电机领域,尤其涉及一种单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法。
背景技术
单绕组无轴承磁通切换电机利用空间对称绕组对气隙磁场进行相反方向调制,从而产生满足转子悬浮运行悬浮力。为了实现转子悬浮状态旋转运行,六相绕组中需要同时流过控制转子切向旋转的转矩电流分量和控制转子径向悬浮的悬浮电流分量,即利用一套绕组同时实现转子转矩和悬浮力控制。该方法有利于电机输出转矩能力的充分发挥,但同时也使得转矩和悬浮力控制相互耦合。
目前实现该种电机的转矩及悬浮力控制一般采用电流闭环控制,需要静止坐标系和旋转坐标系之间的坐标转换,控制算法计算量大,且逆变器输出电压通过电机电感转化为电流控制转矩和悬浮力,属于间接控制,存在滞后。降低了转矩及悬浮力的动态响应,降低了转子悬浮运行的稳定性。
单绕组无轴承磁通切换电机属于多相电机,可利用六相逆变器输出电压矢量分别在转矩及悬浮平面直接实现该种电机转矩和悬浮力的控制,能够解决电流闭环控制存在间接性、滞后性等问题。利用多相电机电压矢量的丰富性,基于开关表可实现多平面的控制要求。但基于开关表的实现方案并不能很好的满足控制要求,转矩及悬浮力脉动均较大。
发明内容
本发明针对六相单绕组无轴承磁通切换电机,提出一种单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法。该方法利用空间矢量调制方法,在双机电能量转换平面上输出满足要求的电压矢量,降低转矩及悬浮脉动。
本发明具体采用以下技术方案:
一种单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于,包括以下步骤:
步骤S1:利用电流传感器及AD转换通道,检测出六相绕组电流iA~iF;利用转子位置角传感器及检测通道,检测出转子位置角θr;利用转子X和Y方向径向位移传感器及检测通道,检测出转子X和Y方向径向位移x、y;
步骤S2:将所述六相绕组电流iA~iF输入T6坐标变换计算模块,输出转矩平面电流iαT、iβT,悬浮平面电流iαS,iβS以及零序平面电流i01,i02
Figure GDA0002839891070000021
其中T6矩阵为:
Figure GDA0002839891070000022
步骤S3:将所述转矩平面电流iαT、iβT,转子位置角θr输入磁链模型计算模块,得到转矩平面磁链ψαTβT和转矩平面磁链幅值|ψsT|:
Figure GDA0002839891070000023
Figure GDA0002839891070000024
Figure GDA0002839891070000025
其中LT是转矩平面电感值,ψfαT,ψfβT是永磁磁链矢量ψf在转矩平面静止坐标系上的投影,|ψf|为永磁磁链矢量幅值;这些电感值和永磁磁链幅值可通过查阅电机设计手册或通过电机有限元分析软件获得。
步骤S4:将电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|、转矩平面磁链幅值|ψST|输入转矩角变化量计算模块,得到转矩角变化量Δδ:
Figure GDA0002839891070000031
其中Pn为电机极对数;
步骤S5:将所述转矩角变化量Δδ、转矩平面磁链ψαT、ψβT、转矩平面磁链幅值给定值
Figure GDA0002839891070000032
转矩平面磁链幅值|ψST|、同步角速度ωr输入转矩平面磁链增量计算模块,得到转矩平面磁链增量ΔψαT、ΔψβT
Figure GDA0002839891070000033
其中Ts为控制周期;
步骤S6:将所述悬浮平面电流iαS,iβS、悬浮力相位差
Figure GDA0002839891070000034
悬浮力系数K、转子位置角θr输入悬浮力计算模块,输出XY坐标系下的实际悬浮力值Fx和Fy
Figure GDA0002839891070000035
步骤S7:将X方向径向位移给定值x*减去X方向径向位移x,Y方向径向位移给定值y*减去Y方向径向位移y;所得差值分别输入X方向径向位移PI控制器和Y方向径向位移PI控制器,输出期望悬浮力值
Figure GDA0002839891070000036
和Fy *
Figure GDA0002839891070000037
其中kp为PI控制器的比例放大系数,kI为PI控制器的积分放大系数,t为时间;
步骤S8:将X方向的期望悬浮力值
Figure GDA0002839891070000038
减去实际悬浮力值Fx,Y方向的期望悬浮力值
Figure GDA0002839891070000039
减去实际悬浮力值Fy;输出悬浮力差值ΔFx和ΔFy
Figure GDA0002839891070000041
步骤S9:将所述悬浮力差值ΔFx和ΔFy、悬浮力相位差
Figure GDA0002839891070000042
悬浮力系数K、转子位置角θr输入悬浮平面磁链增量计算模块,得到悬浮平面磁链增量ΔψαS和ΔψβS
Figure GDA0002839891070000043
其中,LS为电机悬浮平面电感值;
步骤S10:将所述转矩平面磁链增量ΔψαT、ΔψβT和悬浮平面磁链增量ΔψαS、ΔψβS,转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS输入期望电压矢量计算模块,得到转矩平面期望电压矢量
Figure GDA0002839891070000044
转矩平面期望电压矢量幅值
Figure GDA0002839891070000045
转矩平面期望电压矢量角度θT,悬浮平面期望电压矢量
Figure GDA0002839891070000046
悬浮平面期望电压矢量幅值
Figure GDA0002839891070000047
悬浮平面期望电压矢量角度θS
Figure GDA0002839891070000048
Figure GDA0002839891070000049
其中Rs为定子电阻值;
步骤S11:根据六相单绕组无轴承磁通切换电机的基本电压矢量分布,分别将转矩平面和悬浮平面以60°为单位各划分为6个扇区;将所述转矩平面期望电压矢量角度θT输入转矩平面扇区判断模块,得到转矩平面期望电压矢量扇区号kT;将悬浮平面期望电压矢量角度θS输入悬浮平面扇区判断模块,得到悬浮平面期望电压矢量扇区号kS
步骤S12:根据所述转矩平面期望电压矢量扇区号kT,查转矩平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKT1和uKT2;根据悬浮平面期望电压矢量扇区号kS,查悬浮平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKS1和uKS2
步骤S13:将所述转矩平面期望电压矢量幅值
Figure GDA0002839891070000051
悬浮平面期望电压矢量幅值
Figure GDA0002839891070000052
控制周期Ts输入转矩平面作用时间和悬浮平面作用时间计算模块,得到转矩平面作用时间TKT和悬浮平面作用时间TKS
Figure GDA0002839891070000053
Figure GDA0002839891070000054
步骤S14:将所述转矩平面期望电压矢量扇区号kT、转矩平面期望电压矢量角度θT、转矩平面作用时间TKT、转矩平面期望电压矢量幅值
Figure GDA0002839891070000055
直流母线电压UDC输入转矩平面电压矢量作用时间计算模块,得到转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2
Figure GDA0002839891070000056
Figure GDA0002839891070000057
步骤S15:将所述悬浮平面期望电压矢量扇区号kS、悬浮平面期望电压矢量角度θS、悬浮平面作用时间TKS、悬浮平面期望电压矢量幅值
Figure GDA0002839891070000061
直流母线电压UDC输入悬浮平面电压矢量作用时间计算,得到悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2
Figure GDA0002839891070000062
Figure GDA0002839891070000063
步骤S16:将转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2,悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2分别输入转矩平面零矢量作用时间计算模块、悬浮平面零矢量作用时间计算模块,得到转矩平面零矢量作用时间tKT0,悬浮平面零矢量作用时间tKS0
tKT0=TKT-tKT1-tKT2
tKS0=TKS-tKS1-tKS2
步骤S17:由冲量等效原则,将所述转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2、悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2、以及转矩平面零矢量作用时间tKT0和悬浮平面零矢量作用时间tKS0转化为连接电机的逆变桥臂上桥臂导通时间TA~TF,用以控制无轴承磁通切换电机转矩和悬浮平面,实现电机悬浮旋转。
优选地,步骤S4中所述电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|计算方法如下:
步骤S41:将电磁转矩给定值
Figure GDA0002839891070000064
减去电机电磁转矩Te,得到电磁转矩误差ΔTe
Figure GDA0002839891070000071
步骤S42:将转矩平面磁链ψαT、ψβT输入位置角计算模块,输出转矩平面磁链矢量ψST=ψαT+jψβT在转矩平面所处的位置角θψST,再减去转子位置角θr,得到转矩角δ:
Figure GDA0002839891070000072
δ=θψSTr
步骤S43:将转矩平面磁链幅值给定值
Figure GDA0002839891070000073
减去转矩平面磁链幅值|ψST|得到转矩平面磁链幅值增量Δ|ψST|:
Figure GDA0002839891070000074
优选地,步骤S6中所述悬浮力相位差
Figure GDA0002839891070000075
悬浮力系数K计算方法如下:
步骤S61:将所述转矩平面电流iαT,iβT和转子位置角θr输入旋转变换模块,输出转矩平面qT轴电流iqT
iqT=-sinθr×iαT+cosθr×iβT
步骤S62;将所述转矩平面qT轴电流iqT输入悬浮力相位差计算模块,输出悬浮力相位差
Figure GDA0002839891070000076
Figure GDA0002839891070000077
其中,kqT为单位正向悬浮电流与单位qT轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM为单位正向悬浮电流与永磁体产生的悬浮力基波幅值;
步骤S63:将所述转矩平面qT轴电流iqT输入悬浮力系数计算模块,输出悬浮力系数K:
Figure GDA0002839891070000081
优选地,步骤S11中所述转矩平面期望电压矢量扇区号kT和悬浮平面期望电压矢量扇区号kS具体计算方法如下:
Figure GDA0002839891070000082
Figure GDA0002839891070000083
其中,函数int()表示取整数含义,θT∈[0,2π),θS∈[0,2π)。
优选地,实现该方法的硬件系统包括:分别连接无轴承磁通切换电机的六相逆变器、电机转子位置角检测电路和XY径向位移检测电路;六相绕组电流采样电路设置在所述无轴承磁通切换电机和六相逆变器之间;所述电机转子位置角检测电路经QEP测速模块连接中央控制器;所述六相逆变器经隔离驱动模块连接中央控制器,经滤波电容和整流器连接直流电源;直流母线电压采样电路连接六相逆变器;所述中央控制器分别连接上位机、可编程ASIC器件、故障检测模块和AD调理模块;所述故障检测模块连接可编程ASIC器件;所述故障检测模块和AD调理模块采集电压、电流、位移传感器信号。
优选地,所述六相逆变器的直流母线电压采用直流电源;所述六相逆变器的开关管采用IGBT或MOSFET,中央控制器采用DSP或单片机;所述六相绕组电流采样电路采用霍尔电流传感器与运算放大器构成,或采用绕组串功率电阻后接差分运算放大器构成;所述直流母线电压采样电路采用霍尔电压传感器与运算放大器构成,或采用并联电阻分压后接由运算放大器构成的电压跟随器构成;所述电机转子位置角检测电路采用旋转编码器后接电平转换电路构成,或采用旋转变压器后接解码电路构成;所述XY径向位移检测电路采用电涡流传感器后接运算放大器构成,或采用线性光耦后接运算放大器构成;所述可编程ASIC器件采用CPLD或FPGA;所述故障检测模块和AD调理模块均由采用运算放大器构成;将电压、电流信号经过AD调理电路及故障检测模块分别输入中央控制器的A/D采样模块和可编程ASIC器件;转子位移信号经过AD调理电路进行电平转换后输入中央控制器的A/D采样模块;所述位置角检测电路输出的脉冲信号通过QEP测速模块输入中央控制器;所述中央控制器执行步骤S1-步骤S17,输出逆变桥臂开关信号,经由隔离驱动模块控制六相逆变器中的功率开关管的开关动作。
与现有技术相比,本发明及其优选方案有以下有益效果:
(1)分别在转矩控制平面和悬浮力控制平面实现二者的解耦控制,有效地增强了转子悬浮运行的平稳性;
(2)利用逆变器电压矢量直接实现电磁转矩和悬浮力的控制,加速了电磁转矩和悬浮力控制的动态响应;
(3)采用空间矢量调制方式,有效的降低了电磁转矩和悬浮力脉动。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
图1是本发明实施例单绕组无轴承磁通切换电机低转矩及悬浮力直接控制系统框图;
图2是本发明实施例SVPWM调制模块框图;
图3是本发明实施例六相单绕组无轴承磁通切换电机横截面示意图;
图4是本发明实施例六相单绕组无轴承磁通切换电机转矩平面基本电压矢量分布图;
图5是本发明实施例六相单绕组无轴承磁通切换电机悬浮平面基本电压矢量分布图;
图6是本发明实施例驱动系统硬件结构示意图;
图7是本发明实施例转矩平面坐标系定义示意图;
图8是本发明实施例悬浮平面坐标系定义示意图;
图9是本发明实施例转矩平面磁链增量示意图;
图10是本发明实施例悬浮平面磁链增量示意图;
图11是本发明实施例六相单绕组无轴承磁通切换电机完整电压矢量分布图(转矩平面电压矢量uαT+juβT分布);
图12是本发明实施例六相单绕组无轴承磁通切换电机完整电压矢量分布图(悬浮平面电压矢量uαS+juβS分布);
图13是本发明实施例电压矢量合成图。
具体实施方式
为让本专利的特征和优点能更明显易懂,下文特举实施例,作详细说明如下:
本实施例所提的单绕组无轴承磁通切换电机低转矩及悬浮力直接控制系统框图如图1所示。由六相逆变器、无轴承磁通切换电机、磁链模型、转矩模型、位置角计算模块、转矩角变化量计算模块、转矩平面磁链增量计算模块、悬浮力系数计算模块、悬浮力相位差计算模块、悬浮力计算模块、悬浮平面磁链增量计算模块、X方向径向位移PI控制器、Y方向径向位移PI控制器、期望电压矢量计算模块、SVPWM调制模块等环节构成。
其中,电机的转子位置角θr通过转子位置角传感器及检测通道获得。电机转子X和Y方向径向位移x、y通过转子X和Y方向径向位移传感器及检测通道获得。定子六相绕组电流iA~iF通过电流传感器及检测通道获得。将自然坐标系下的定子六相绕组电流iA~iF送入T6坐标变换模块,输出转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS;把转矩平面电流iαT、iβT送入磁链模型计算模块,输出转矩平面磁链ψαT、ψβT及其对应的幅值|ψST|;转矩平面磁链幅值给定值
Figure GDA0002839891070000101
减去转矩平面磁链幅值|ψST|得到转矩平面磁链幅值增量Δ|ψST|。把转矩平面磁链ψαT、ψβT送入位置角计算模块,输出转矩平面磁链矢量ψST=ψαT+jψβT在转矩平面所处的位置角θψST,再减去转子位置角θr,得到转矩角δ。把转矩平面磁链ψαT、ψβT,转矩平面电流iαT、iβT送入转矩模型计算模块,得到电机电磁转矩Te。电磁转矩给定值
Figure GDA0002839891070000102
减去电机电磁转矩Te,得到电磁转矩误差ΔTe。将电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|、转矩平面磁链幅值|ψST|送入转矩角变化量计算模块,得到转矩角变化量Δδ;把转矩角变化量Δδ、转矩平面磁链ψαT、ψβT,转矩平面磁链幅值给定值
Figure GDA0002839891070000103
转矩平面磁链幅值|ψST|送入转矩平面磁链增量计算模块,得到转矩平面磁链增量ΔψαT、ΔψβT;将X方向径向位移给定值x*减去X方向径向位移x,得到的差值送入X方向径向位移PI控制器,输出X方向期望悬浮力值
Figure GDA00028398910700001114
将Y方向径向位移给定值y*减去Y方向径向位移y,得到的差值送入Y方向径向位移PI控制器,输出Y方向期望悬浮力值
Figure GDA0002839891070000112
将悬浮平面电流iαT、iβT和转子位置角θr送入旋转变换模块,输出转矩平面旋转坐标系下qT轴电流iqT,再将该值送入悬浮力相位差
Figure GDA0002839891070000113
计算模块,输出悬浮力相位差
Figure GDA0002839891070000114
将转矩平面旋转坐标系下qT轴电流iqT送入悬浮力系数K计算模块,输出悬浮力系数K。把转子位置角θr,悬浮力相位差
Figure GDA0002839891070000115
悬浮力系数K以及悬浮平面电流iαS、iβS送入悬浮力计算模块,输出XY坐标系下的悬浮力Fx和Fy
其中XY坐标系定义见图3,X轴线与A1线圈轴线重合。将期望悬浮力值
Figure GDA0002839891070000116
Figure GDA0002839891070000117
分别减去实际悬浮力值Fx和Fy,得到的差值送入悬浮平面磁链增量计算模块,输出悬浮平面磁链增量ΔψαS、ΔψβS;将转矩平面磁链增量ΔψαT、ΔψβT和悬浮平面磁链增量ΔψαS、ΔψβS送入期望电压矢量计算模块,并经SVPWM调制模块输出逆变器各上桥臂导通时间TA~TF以控制电机。
其中期望电压矢量计算模块和SVPWM调制模块框图如图2所示。将转矩平面磁链增量ΔψαT、ΔψβT和悬浮平面磁链增量ΔψαS、ΔψβS,转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS,控制周期Ts,定子电阻值Rs送入期望电压矢量计算模块。得到转矩平面期望电压矢量
Figure GDA0002839891070000118
转矩平面期望电压矢量幅值
Figure GDA0002839891070000119
转矩平面期望电压矢量角度θT,悬浮平面期望电压矢量
Figure GDA00028398910700001110
悬浮平面期望电压矢量幅值
Figure GDA00028398910700001111
悬浮平面期望电压矢量角度θS;将转矩平面期望电压矢量幅值
Figure GDA00028398910700001112
悬浮平面期望电压矢量幅值
Figure GDA00028398910700001113
和控制周期Ts送入转矩平面作用时间和悬浮平面作用时间计算模块,得到转矩平面作用时间TKT和悬浮平面作用时间TKS。将转矩平面期望电压矢量角度θT送入转矩平面扇区判断,得到转矩平面期望电压矢量扇区号kT,其中转矩平面电压矢量扇区划分如图4所示,将悬浮平面期望电压矢量角度θS送入悬浮平面扇区判断,得到悬浮平面期望电压矢量扇区号kS,其中悬浮平面电压矢量扇区划分如图5所示;将转矩平面期望电压矢量扇区号kT、转矩平面期望电压矢量角度θT、转矩平面作用时间TKT、转矩平面期望电压矢量幅值
Figure GDA0002839891070000121
直流母线电压UDC送入转矩平面电压矢量作用时间计算模块,得到转矩平面电压矢量作用时间tKT1和tKT2。将悬浮平面期望电压矢量扇区号kS、悬浮平面期望电压矢量角度θS、悬浮平面作用时间TKS、悬浮平面期望电压矢量幅值
Figure GDA0002839891070000122
直流母线电压UDC送入悬浮平面电压矢量作用时间计算模块,得到悬浮平面电压矢量作用时间tKS1和tKS2;将转矩平面电压矢量作用时间tKT1和tKT2送入转矩平面零矢量作用时间计算模块,得到转矩平面零矢量作用时间tKT0。将悬浮平面电压矢量作用时间tKS1和tKS2送入悬浮平面零矢量作用时间计算模块,得到悬浮平面零矢量作用时间tKS0;根据转矩平面期望电压矢量扇区号kT,查转矩平面基本电压矢量分布图,如图4所示,选择转矩平面基本电压矢量uKT1和uKT2。根据悬浮平面期望电压矢量扇区号kS,查悬浮平面基本电压矢量分布图,如图5所示,选择悬浮平面基本电压矢量uKS1和uKS2;将转矩平面基本电压矢量uKT1和uKT2、转矩平面零矢量作用时间tKT0、转矩平面电压矢量作用时间tKT1和tKT2、悬浮平面基本电压矢量uKS1和uKS2、悬浮平面零矢量作用时间tKS0、悬浮平面电压矢量作用时间tKS1和tKS2送入逆变桥臂占空比合成,根据伏秒乘积等效原则,输出逆变器各上桥臂导通时间TA~TF,用以控制无轴承磁通切换电机转矩和悬浮平面,实现电机悬浮旋转。
电机结构如图3所示,电机12个U型铁芯,每个U型铁芯之间夹一个沿切向充磁的永磁体,充磁方向交替相反,转子有10个齿。定子每相绕组串联绕在空间相互垂直的定子齿上,构成6相对称绕组。其中A相和D相绕组空间对称、B相和E相绕组空间对称、C相和F相绕组空间对称,六相绕组轴线空间上互差60°机械角。若需要调节电机转速,则利用转速闭环控制输出转矩给定值
Figure GDA0002839891070000131
本发明的实施例驱动系统硬件结构如图6所示。包括:整流器、滤波电容、直流母线电压采样电路、六相逆变器、无轴承磁通切换电机、六相绕组电流采样电路、电机转子位置角检测电路、XY径向位移检测电路、隔离驱动、中央控制器、上位机等。其中六相逆变器直流母线电压也可以采用合适的直流电源提供。六相逆变器的开关管可采用IGBT或MOSFET,中央控制器采用DSP或单片机。电流采样电路采用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成。采用霍尔方案可以有效实现控制回路与主回路的电气隔离,采用绕组串功率电阻方案可以降低驱动系统成本。直流母线电压采样电路采用霍尔电压传感器与运算放大器相结合方式构成,也可以采用并联电阻分压后接由运算放大器构成的电压跟随器相结合方式构成。转子位置角检测电路可以采用旋转编码器后接电平转换电路构成,也可以采用旋转变压器后接解码电路构成,其中前者成本较低,但位置角采样精度受编码器线数限制,而后者成本较高,但位置角采样精度较高。XY径向位移检测电路采用电涡流传感器后接运算放大器相结合方式构成,也可以采用线性光耦后接运算放大器相结合方式构成。可编程ASIC器件可以采用CPLD或FPGA。故障检测模块和AD调理模块均可采用运算放大器构成。将电压、电流信号经过AD调理电路及故障检测模块分别送到中央控制器的A/D采样模块和可编程ASIC器件。转子位移信号经过AD调理电路进行电平转换后送至中央控制器的A/D采样模块。位置角检测电路输出的脉冲信号通过QEP测速模块送至中央控制器。根据取得的信号和本发明的直接控制方法,输出逆变桥臂开关信号,经由隔离驱动去控制逆变器中的功率开关管的开关动作。
图7、图8为无轴承磁通切换电机的坐标系定义。其中图7为转矩平面坐标系定义,图8为悬浮平面坐标系定义。
如图7所示,A~F为电机六相绕组轴线。根据A相中耦合永磁磁链正最大位置判断得A相绕组轴线处于图3中转子齿中心线逆时针超前A1线圈轴线9°机械角处。六相绕组轴线空间上互差60°机械角。αTβT为转矩平面静止坐标系,dTqT为旋转坐标系。其中αT轴与A相绕组轴线重合。dT轴与永磁体在六相绕组轴线中耦合合成磁链矢量ψf方向重合,与αT轴沿逆时针方向夹角表示为转子位置角θr。转矩平面磁链矢量ψST=ψαT+jψβT。转矩平面电流矢量iST=iαT+jiβT。转矩角δ为转矩平面磁链矢量ψST和永磁磁链矢量ψf的夹角。
悬浮平面坐标系定义如图8所示。XY为水平-垂直直角坐标系,X轴与图二中的A1线圈轴线方向一致,与A相绕组主线相差9°机械角。αSβS为悬浮平面静止坐标系,dSqS为旋转坐标系。αS轴和dS轴之间的夹角为
Figure GDA0002839891070000141
其中
Figure GDA0002839891070000142
为悬浮力相位差。悬浮平面磁链矢量ψSS=ψαS+jψβS。悬浮平面电流矢量iSS=iαS+jiβS。虚拟磁链矢量
Figure GDA0002839891070000143
其中K为悬浮力系数。悬浮角γ为悬浮平面磁链矢量ψSS和虚拟磁链矢量ψ的夹角。
利用T6恒功率矩阵,由六相绕组电流iA~iF得到转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS,零序平面电流i01、i02
Figure GDA0002839891070000144
其中,T6恒功率矩阵为:
Figure GDA0002839891070000145
转矩平面磁链ψαTβT为:
Figure GDA0002839891070000151
其中LT是转矩平面电感值,ψfαT,ψfβT是永磁磁链矢量ψf在转矩平面静止坐标系上的投影
Figure GDA0002839891070000152
其中θr为转子位置角。|ψf|为永磁磁链矢量幅值。转矩平面电感值LT和永磁磁链矢量幅值|ψf|可通过电机有限元分析软件获得。
电机电磁转矩为:
Te=pnαTiβTβTiαT) (4)
其中,pn为电机极数。
利用转子位置角θr,将转矩平面电流iαT、iβT,转矩平面磁链ψαT、ψβT,电磁转矩Te旋转变换至转矩平面旋转坐标系dTqT。得到旋转坐标系dTqT下的电流idT、idT,磁链ψdT、ψqT和电磁转矩Te
Figure GDA0002839891070000153
Figure GDA0002839891070000154
其中Ld为电机转矩平面dT轴电感,Lq为电机转矩平面qT轴电感
Te=pndTiqTqTidT) (7)
由图7可以得到转矩平面磁链矢量ψST和旋转坐标系dTqT下的磁链ψdT、ψqT的关系:
ψST=ψdT+jψdT (8)
Figure GDA0002839891070000155
其中,|ψST|为转矩平面磁链幅值,转矩角δ为转矩平面磁链矢量ψST与dT轴之间的夹角。
将式(9)代入(7)中,得到变形后的电磁转矩:
Te=PnST|(iqTcosδ-idTsinδ) (10)
再将式(6)和式(9)代入式(10)中,用磁链代替电流得到:
Figure GDA0002839891070000161
由直接转矩控制理论可知,保持转矩控制平面定子磁链幅值|ψST|恒定,控制转矩角δ,即可控制电磁转矩的大小。
根据电机产生的悬浮力原理可得XY坐标系下的悬浮力Fx、Fy方程:
Figure GDA0002839891070000162
K为悬浮力系数:
Figure GDA0002839891070000163
Figure GDA0002839891070000164
为悬浮力相位差:
Figure GDA0002839891070000165
其中,kqT为单位正向悬浮电流与单位qT轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM为单位正向悬浮电流与永磁体产生的悬浮力基波幅值。二者皆可通过电机有限元分析软件获得。
将(12)式展开可得到:
Figure GDA0002839891070000166
由前文对悬浮平面坐标系定义可知虚拟磁链矢量ψ,悬浮平面电流矢量iS,悬浮平面磁链矢量ψSS为:
Figure GDA0002839891070000171
iS=iαS+jiβS (17)
ψSS=L2·iS (18)
其中,L2为悬浮平面电感值。
将(16)(17)(18)代入(15)式可得:
Figure GDA0002839891070000172
通过控制悬浮平面磁链幅值|ψSS|大小以及悬浮角γ即可控制悬浮力的大小和方向。
为了实现低转矩及悬浮力直接控制,采用基于空间矢量调制的直接控制策略。其基本思路是根据转矩及悬浮力控制的需要,求出在一个周期内电机实际需要的期望电压矢量,利用平行四边形矢量相加原则,将两个及以上的基本电压矢量各工作一定时间来合成所需要的一定大小和方向的电压矢量,精确地补偿转矩和磁链增量。
上述思路的第一个关键是空间期望电压矢量的计算模块。图9为转矩平面磁链增量示意图。ψST为当前转矩平面磁链矢量,ΔψST为转矩平面磁链增量矢量,
Figure GDA0002839891070000173
为期望转矩平面磁链矢量。
从图中可以看出:
Figure GDA0002839891070000174
进一步可得到在转矩平面磁链增量ΔψαT、ΔψβT为:
Figure GDA0002839891070000181
Figure GDA0002839891070000182
带入(21)式可得:
Figure GDA0002839891070000183
式中,ωr为同步角速度,Ts为控制周期,转矩平面磁链幅值给定值
Figure GDA0002839891070000184
为定值。转矩角变化量Δδ可通过电磁转矩求得。
Figure GDA0002839891070000185
实际控制中,控制方式属于离散控制,控制周期为Ts,则将式(23)离散化可得:
Figure GDA0002839891070000186
由式(24)可求得转矩角变化量Δδ为:
Figure GDA0002839891070000191
当转矩平面产生式(22)所述的磁链增量时,系统所加转矩平面期望电压矢量
Figure GDA0002839891070000192
为:
Figure GDA0002839891070000193
图10为悬浮平面磁链增量示意图。其中F为当前悬浮力矢量,F*为期望悬浮力矢量,悬浮力矢量差ΔF=F*-F。ψSS为当前悬浮平面磁链矢量,
Figure GDA0002839891070000194
为期望悬浮平面磁链矢量。悬浮磁链矢量差
Figure GDA0002839891070000195
当前悬浮平面磁链矢量ψSS与αS轴夹角为γ+u,其中
Figure GDA0002839891070000196
当前悬浮力矢量F在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000197
其中,
Figure GDA0002839891070000198
期望悬浮力矢量F*在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000199
悬浮力矢量差ΔF在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000201
当前悬浮磁链矢量ψSS在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000202
期望悬浮磁链矢量
Figure GDA0002839891070000203
在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000204
悬浮磁链矢量差ΔψSS在悬浮平面αSβS坐标系下的分量为:
Figure GDA0002839891070000205
将式(29)代入式(32)可得:
Figure GDA0002839891070000206
式中,
Figure GDA0002839891070000207
当悬浮平面产生式(33)所述的磁链增量时,系统所加悬浮平面期望电压矢量
Figure GDA0002839891070000208
为:
Figure GDA0002839891070000211
在求出转矩和悬浮平面期望电压矢量后,基于空间矢量调制的直接控制策略的第二个关键是生成六相逆变桥臂的开关信号,用以驱动六相逆变桥臂,实现电机的悬浮旋转。
在αTβTαSβS0102轴系中,六相相电压如下式:
Figure GDA0002839891070000212
则可计算得到:
Figure GDA0002839891070000213
根据式(36)可得到转矩平面和悬浮平面完整电压矢量分布如图11、图12所示。
为在一个控制周期TS内实现对两个平面磁链的精确控制,在一个周期内进一步分时段控制不同的平面。每个平面的作用时间由期望电压矢量幅值按比例求取:
Figure GDA0002839891070000221
在转矩平面控制时间TKT内,从六相单绕组无轴承磁通切换电机的64个电压矢量中选取2个及以上电压矢量合成转矩平面期望电压矢量
Figure GDA0002839891070000222
这几个电压矢量在悬浮平面的合成矢量为零矢量。同理,在悬浮平面控制时间TKS内,选取数个电压矢量合成悬浮平面期望电压矢量
Figure GDA0002839891070000223
在悬浮平面的合成矢量为零矢量。基于上述原则,从六相单绕组无轴承磁通切换电机的完整电压矢量分布图中,选取符合要求的电压矢量构成图4、图5所示的基本电压矢量分布图。从图4、图5中可以看出,所选取的电压矢量在转矩平面内为有效电压矢量时,在悬浮平面内为零矢量;在悬浮平面内为有效电压矢量时,在转矩平面内为零矢量。
确定完工作平面作用时间以及基本电压矢量后,进一步的需要确定各电压矢量的具体作用时间。现以转矩平面第一扇区为例作具体说明。当转矩平面期望电压矢量
Figure GDA0002839891070000224
位于第一扇区时,如图13所示的电压矢量合成图。
由伏秒平衡原理可知
Figure GDA0002839891070000225
其中,uKT1和uKT2分别对应期望电压矢量所在扇区两侧的基本电压矢量。由图12所示的几何关系可得:
Figure GDA0002839891070000231
进一步推广至转矩平面其他扇区可得:
Figure GDA0002839891070000232
同理可得合成悬浮平面期望电压矢量的两个有效电压矢量工作时间:
Figure GDA0002839891070000233
转矩平面和悬浮平面零电压矢量工作时间为:
Figure GDA0002839891070000234
根据冲量等效原则,可将两个平面内的各电压矢量作用时间直接转化为各相导通PWM波总的作用时间,进而控制电机实现转子的悬浮旋转。
本实施例具体的工作过程包括如下步骤:
(1)利用电流传感器及AD转换通道,检测出六相绕组电流iA~iF;利用转子位置角传感器及检测通道,检测出转子位置角θr;利用转子X和Y方向径向位移传感器及检测通道,检测出转子X和Y方向径向位移x、y。
(2)把六相绕组电流iA~iF送入T6坐标变换计算模块,输出转矩平面电流iαT、iβT,悬浮平面电流iαS,iβS以及零序平面电流i01,i02
Figure GDA0002839891070000241
其中T6矩阵为:
Figure GDA0002839891070000242
(3)把转矩平面电流iαT、iβT,转子位置角θr送入磁链模型计算模块,得到转矩平面磁链ψαTβT和转矩平面磁链幅值|ψsT|:
Figure GDA0002839891070000243
Figure GDA0002839891070000244
Figure GDA0002839891070000245
其中LT是转矩平面电感值,ψfαT,ψfβT是永磁磁链矢量ψf在转矩平面静止坐标系上的投影,|ψf|为永磁磁链矢量幅值。这些电感值和永磁磁链幅值可通过查阅电机设计手册或通过电机有限元分析软件获得。
(4)把电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|、转矩平面磁链幅值|ψST|送入转矩角变化量计算模块,得到转矩角变化量Δδ:
Figure GDA0002839891070000246
其中Pn为电机极对数。
(5)把转矩角变化量Δδ、转矩平面磁链ψαT、ψβT、转矩平面磁链幅值给定值
Figure GDA0002839891070000251
转矩平面磁链幅值|ψST|、同步角速度ωr送入转矩平面磁链增量计算模块,得到转矩平面磁链增量ΔψαT、ΔψβT
Figure GDA0002839891070000252
其中Ts为控制周期。
(6)把悬浮平面电流iαS,iβS、悬浮力相位差
Figure GDA0002839891070000253
悬浮力系数K、转子位置角θr送入悬浮力计算模块,输出XY坐标系下的实际悬浮力值Fx和Fy
Figure GDA0002839891070000254
(7)把X方向径向位移给定值x*减去X方向径向位移x,Y方向径向位移给定值y*减去Y方向径向位移y。所得差值分别送入X方向径向位移PI控制器和Y方向径向位移PI控制器,输出期望悬浮力值
Figure GDA0002839891070000255
和Fy *
Figure GDA0002839891070000256
其中kp为PI控制器的比例放大系数,kI为PI控制器的积分放大系数
(8)把X方向的期望悬浮力值
Figure GDA0002839891070000257
减去实际悬浮力值Fx,Y方向的期望悬浮力值
Figure GDA0002839891070000258
减去实际悬浮力值Fy。输出悬浮力差值ΔFx和ΔFy
Figure GDA0002839891070000259
(9)把悬浮力差值ΔFx和ΔFy、悬浮力相位差
Figure GDA00028398910700002510
悬浮力系数K、转子位置角θr送入悬浮平面磁链增量计算模块,得到悬浮平面磁链增量ΔψαS和ΔψβS
Figure GDA0002839891070000261
其中,LS为电机悬浮平面电感值,该值可通过电机有限元分析软件获得。
(10)将转矩平面磁链增量ΔψαT、ΔψβT和悬浮平面磁链增量ΔψαS、ΔψβS,转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS送入期望电压矢量计算模块。得到转矩平面期望电压矢量
Figure GDA0002839891070000262
转矩平面期望电压矢量幅值
Figure GDA0002839891070000263
转矩平面期望电压矢量角度θT,悬浮平面期望电压矢量
Figure GDA0002839891070000264
悬浮平面期望电压矢量幅值
Figure GDA0002839891070000265
悬浮平面期望电压矢量角度θS
Figure GDA0002839891070000266
Figure GDA0002839891070000267
其中Ts为控制周期,Rs为定子电阻值。
(11)根据六相单绕组无轴承磁通切换电机基本电压矢量分布,分别将转矩平面和悬浮平面以60°为单位各划分为6个扇区,如图4、图5所示。将转矩平面期望电压矢量角度θT送入转矩平面扇区判断,得到转矩平面期望电压矢量扇区号kT;将悬浮平面期望电压矢量角度θS送入悬浮平面扇区判断,得到悬浮平面期望电压矢量扇区号kS
(12)根据转矩平面期望电压矢量扇区号kT,查转矩平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKT1和uKT2;根据悬浮平面期望电压矢量扇区号kS,查悬浮平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKS1和uKS2
(13)把转矩平面期望电压矢量幅值
Figure GDA0002839891070000271
悬浮平面期望电压矢量幅值
Figure GDA0002839891070000272
控制周期Ts送入转矩平面作用时间和悬浮平面作用时间计算模块,得到转矩平面作用时间TKT和悬浮平面作用时间TKS
Figure GDA0002839891070000273
Figure GDA0002839891070000274
(14)将转矩平面期望电压矢量扇区号kT、转矩平面期望电压矢量角度θT、转矩平面作用时间TKT、转矩平面期望电压矢量幅值
Figure GDA0002839891070000275
直流母线电压UDC送入转矩平面电压矢量作用时间计算模块,得到转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2
Figure GDA0002839891070000276
Figure GDA0002839891070000277
(15)将悬浮平面期望电压矢量扇区号kS、悬浮平面期望电压矢量角度θS、悬浮平面作用时间TKS、悬浮平面期望电压矢量幅值
Figure GDA0002839891070000278
直流母线电压UDC送入悬浮平面电压矢量作用时间计算模块,得到悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2
Figure GDA0002839891070000281
Figure GDA0002839891070000282
(16)转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2,悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2分别送入转矩平面零矢量作用时间、悬浮平面零矢量作用时间计算模块,得到转矩平面零矢量作用时间tKT0,悬浮平面零矢量作用时间tKS0
tKT0=TKT-tKT1-tKT2
tKS0=TKS-tKS1-tKS2
(17)由冲量等效原则,将转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2、悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2、以及转矩平面零矢量作用时间tKT0和悬浮平面零矢量作用时间tKS0转化为逆变桥臂上桥臂导通时间TA~TF,用以控制无轴承磁通切换电机转矩和悬浮平面,实现电机悬浮旋转。
上述步骤(4)中所需的电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|计算方法如下:
(4.1)电磁转矩给定值
Figure GDA0002839891070000283
减去电机电磁转矩Te,得到电磁转矩误差ΔTe
Figure GDA0002839891070000284
(4.2)把转矩平面磁链ψαT、ψβT送入位置角计算模块,输出转矩平面磁链矢量ψST=ψαT+jψβT在转矩平面所处的位置角θψST,再减去转子位置角θr,得到转矩角δ:
Figure GDA0002839891070000285
δ=θψSTr
(4.3)把转矩平面磁链幅值给定值
Figure GDA0002839891070000291
减去转矩平面磁链幅值|ψST|得到转矩平面磁链幅值增量Δ|ψST|:
Figure GDA0002839891070000292
上述步骤(6)中所需的悬浮力相位差
Figure GDA0002839891070000293
悬浮力系数K计算方法如下:
(6.1)把转矩平面电流iαT,iβT和转子位置角θr送入旋转变换模块,输出转矩平面qT轴电流iqT
iqT=-sinθr×iαT+cosθr×iβT
(6.2)把转矩平面qT轴电流iqT送入悬浮力相位差计算模块,输出悬浮力相位差
Figure GDA0002839891070000294
Figure GDA0002839891070000295
其中,kqT为单位正向悬浮电流与单位qT轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM为单位正向悬浮电流与永磁体产生的悬浮力基波幅值。二者皆可通过电机有限元分析软件获得
(6.3)把转矩平面qT轴电流iqT送入悬浮力系数计算模块,输出悬浮力系数K:
Figure GDA0002839891070000296
上述步骤(9)中所需的悬浮力相位差
Figure GDA0002839891070000297
悬浮力系数K计算方法与步骤(6.1)(6.2)(6.3)相同。
上述步骤(11)中转矩平面期望电压矢量扇区号kT和悬浮平面期望电压矢量扇区号kS具体计算方法如下:
Figure GDA0002839891070000301
Figure GDA0002839891070000302
其中,函数int()表示取整数含义,θT∈[0,2π),θS∈[0,2π)。
本专利不局限于上述最佳实施方式,任何人在本专利的启示下都可以得出其它各种形式的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本专利的涵盖范围。

Claims (6)

1.一种单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于,包括以下步骤:
步骤S1:利用电流传感器及AD转换通道,检测出六相绕组电流iA~iF;利用转子位置角传感器及检测通道,检测出转子位置角θr;利用转子X和Y方向径向位移传感器及检测通道,检测出转子X和Y方向径向位移x、y;
步骤S2:将所述六相绕组电流iA~iF输入T6坐标变换计算模块,输出转矩平面电流iαT、iβT,悬浮平面电流iαS,iβS以及零序平面电流i01,i02
Figure FDA0002839891060000011
其中T6矩阵为:
Figure FDA0002839891060000012
步骤S3:将所述转矩平面电流iαT、iβT,转子位置角θr输入磁链模型计算模块,得到转矩平面磁链ψαTβT和转矩平面磁链幅值|ψsT|:
Figure FDA0002839891060000013
Figure FDA0002839891060000014
Figure FDA0002839891060000015
其中LT是转矩平面电感值,ψfαT,ψfβT是永磁磁链矢量ψf在转矩平面静止坐标系上的投影,|ψf|为永磁磁链矢量幅值;
步骤S4:将电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|、转矩平面磁链幅值|ψST|输入转矩角变化量计算模块,得到转矩角变化量Δδ:
Figure FDA0002839891060000021
其中Pn为电机极对数;
步骤S5:将所述转矩角变化量Δδ、转矩平面磁链ψαT、ψβT、转矩平面磁链幅值给定值
Figure FDA0002839891060000022
转矩平面磁链幅值|ψST|、同步角速度ωr输入转矩平面磁链增量计算模块,得到转矩平面磁链增量ΔψαT、ΔψβT
Figure FDA0002839891060000023
其中Ts为控制周期;
步骤S6:将所述悬浮平面电流iαS,iβS、悬浮力相位差
Figure FDA0002839891060000024
悬浮力系数K、转子位置角θr输入悬浮力计算模块,输出XY坐标系下的实际悬浮力值Fx和Fy
Figure FDA0002839891060000025
步骤S7:将X方向径向位移给定值x*减去X方向径向位移x,Y方向径向位移给定值y*减去Y方向径向位移y;所得差值分别输入X方向径向位移PI控制器和Y方向径向位移PI控制器,输出期望悬浮力值
Figure FDA0002839891060000026
Figure FDA0002839891060000027
Figure FDA0002839891060000028
其中kp为PI控制器的比例放大系数,kI为PI控制器的积分放大系数,t为时间;
步骤S8:将X方向的期望悬浮力值
Figure FDA0002839891060000031
减去实际悬浮力值Fx,Y方向的期望悬浮力值
Figure FDA0002839891060000032
减去实际悬浮力值Fy;输出悬浮力差值ΔFx和ΔFy
Figure FDA0002839891060000033
步骤S9:将所述悬浮力差值ΔFx和ΔFy、悬浮力相位差
Figure FDA0002839891060000034
悬浮力系数K、转子位置角θr输入悬浮平面磁链增量计算模块,得到悬浮平面磁链增量ΔψαS和ΔψβS
Figure FDA0002839891060000035
其中,LS为电机悬浮平面电感值;
步骤S10:将所述转矩平面磁链增量ΔψαT、ΔψβT和悬浮平面磁链增量ΔψαS、ΔψβS,转矩平面电流iαT、iβT,悬浮平面电流iαS、iβS输入期望电压矢量计算模块,得到转矩平面期望电压矢量
Figure FDA0002839891060000036
转矩平面期望电压矢量幅值
Figure FDA0002839891060000037
转矩平面期望电压矢量角度θT,悬浮平面期望电压矢量
Figure FDA0002839891060000038
悬浮平面期望电压矢量幅值
Figure FDA0002839891060000039
悬浮平面期望电压矢量角度θS
Figure FDA00028398910600000310
Figure FDA0002839891060000041
其中Rs为定子电阻值;
步骤S11:根据六相单绕组无轴承磁通切换电机的基本电压矢量分布,分别将转矩平面和悬浮平面以60°为单位各划分为6个扇区;将所述转矩平面期望电压矢量角度θT输入转矩平面扇区判断模块,得到转矩平面期望电压矢量扇区号kT;将悬浮平面期望电压矢量角度θS输入悬浮平面扇区判断模块,得到悬浮平面期望电压矢量扇区号kS
步骤S12:根据所述转矩平面期望电压矢量扇区号kT,查转矩平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKT1和uKT2;根据悬浮平面期望电压矢量扇区号kS,查悬浮平面基本电压矢量分布图,选择其所在扇区相邻的两个矢量uKS1和uKS2
步骤S13:将所述转矩平面期望电压矢量幅值
Figure FDA0002839891060000042
悬浮平面期望电压矢量幅值
Figure FDA0002839891060000043
控制周期Ts输入转矩平面作用时间和悬浮平面作用时间计算模块,得到转矩平面作用时间TKT和悬浮平面作用时间TKS
Figure FDA0002839891060000044
Figure FDA0002839891060000045
步骤S14:将所述转矩平面期望电压矢量扇区号kT、转矩平面期望电压矢量角度θT、转矩平面作用时间TKT、转矩平面期望电压矢量幅值
Figure FDA0002839891060000046
直流母线电压UDC输入转矩平面电压矢量作用时间计算模块,得到转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2
Figure FDA0002839891060000051
Figure FDA0002839891060000052
步骤S15:将所述悬浮平面期望电压矢量扇区号kS、悬浮平面期望电压矢量角度θS、悬浮平面作用时间TKS、悬浮平面期望电压矢量幅值
Figure FDA0002839891060000053
直流母线电压UDC输入悬浮平面电压矢量作用时间计算,得到悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2
Figure FDA0002839891060000054
Figure FDA0002839891060000055
步骤S16:将转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2,悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2分别输入转矩平面零矢量作用时间计算模块、悬浮平面零矢量作用时间计算模块,得到转矩平面零矢量作用时间tKT0,悬浮平面零矢量作用时间tKS0
tKT0=TKT-tKT1-tKT2
tKS0=TKS-tKS1-tKS2
步骤S17:由冲量等效原则,将所述转矩平面基本电压矢量uKT1和uKT2的作用时间tKT1和tKT2、悬浮平面基本电压矢量uKS1和uKS2作用时间tKS1和tKS2、以及转矩平面零矢量作用时间tKT0和悬浮平面零矢量作用时间tKS0转化为连接电机的逆变桥臂上桥臂导通时间TA~TF,用以控制无轴承磁通切换电机转矩和悬浮平面,实现电机悬浮旋转。
2.根据权利要求1所述的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于:
步骤S4中所述电磁转矩误差ΔTe、转矩角δ、转矩平面磁链幅值增量Δ|ψST|计算方法如下:
步骤S41:将电磁转矩给定值
Figure FDA0002839891060000061
减去电机电磁转矩Te,得到电磁转矩误差ΔTe
Figure FDA0002839891060000062
步骤S42:将转矩平面磁链ψαT、ψβT输入位置角计算模块,输出转矩平面磁链矢量ψST=ψαT+jψβT在转矩平面所处的位置角
Figure FDA0002839891060000068
再减去转子位置角θr,得到转矩角δ:
Figure FDA0002839891060000063
Figure FDA0002839891060000064
步骤S43:将转矩平面磁链幅值给定值
Figure FDA0002839891060000065
减去转矩平面磁链幅值|ψST|得到转矩平面磁链幅值增量Δ|ψST|:
Figure FDA0002839891060000066
3.根据权利要求1所述的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于:
步骤S6中所述悬浮力相位差
Figure FDA0002839891060000067
悬浮力系数K计算方法如下:
步骤S61:将所述转矩平面电流iαT,iβT和转子位置角θr输入旋转变换模块,输出转矩平面qT轴电流iqT
iqT=-sinθr×iαT+cosθr×iβT
步骤S62;将所述转矩平面qT轴电流iqT输入悬浮力相位差计算模块,输出悬浮力相位差
Figure FDA0002839891060000071
Figure FDA0002839891060000072
其中,kqT为单位正向悬浮电流与单位qT轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM为单位正向悬浮电流与永磁体产生的悬浮力基波幅值;
步骤S63:将所述转矩平面qT轴电流iqT输入悬浮力系数计算模块,输出悬浮力系数K:
Figure FDA0002839891060000073
4.根据权利要求1所述的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于:
步骤S11中所述转矩平面期望电压矢量扇区号kT和悬浮平面期望电压矢量扇区号kS具体计算方法如下:
Figure FDA0002839891060000074
Figure FDA0002839891060000075
其中,函数int()表示取整数含义,θT∈[0,2π),θS∈[0,2π)。
5.根据权利要求1所述的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于:实现该方法的硬件系统包括:分别连接无轴承磁通切换电机的六相逆变器、电机转子位置角检测电路和XY径向位移检测电路;六相绕组电流采样电路设置在所述无轴承磁通切换电机和六相逆变器之间;所述电机转子位置角检测电路经QEP测速模块连接中央控制器;所述六相逆变器经隔离驱动模块连接中央控制器,经滤波电容和整流器连接直流电源;直流母线电压采样电路连接六相逆变器;所述中央控制器分别连接上位机、可编程ASIC器件、故障检测模块和AD调理模块;所述故障检测模块连接可编程ASIC器件;所述故障检测模块和AD调理模块采集电压、电流、位移传感器信号。
6.根据权利要求5所述的单绕组无轴承磁通切换电机低转矩及悬浮力的直接控制方法,其特征在于:所述六相逆变器的直流母线电压采用直流电源;所述六相逆变器的开关管采用IGBT或MOSFET,中央控制器采用DSP或单片机;所述六相绕组电流采样电路采用霍尔电流传感器与运算放大器构成,或采用绕组串功率电阻后接差分运算放大器构成;所述直流母线电压采样电路采用霍尔电压传感器与运算放大器构成,或采用并联电阻分压后接由运算放大器构成的电压跟随器构成;所述电机转子位置角检测电路采用旋转编码器后接电平转换电路构成,或采用旋转变压器后接解码电路构成;所述XY径向位移检测电路采用电涡流传感器后接运算放大器构成,或采用线性光耦后接运算放大器构成;所述可编程ASIC器件采用CPLD或FPGA;所述故障检测模块和AD调理模块均由运算放大器构成;将电压、电流信号经过AD调理电路及故障检测模块分别输入中央控制器的A/D采样模块和可编程ASIC器件;转子位移信号经过AD调理电路进行电平转换后输入中央控制器的A/D采样模块;所述位置角检测电路输出的脉冲信号通过QEP测速模块输入中央控制器;所述中央控制器执行步骤S1-步骤S17,输出逆变桥臂开关信号,经由隔离驱动模块控制六相逆变器中的功率开关管的开关动作。
CN202010406190.2A 2020-05-14 2020-05-14 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法 Active CN111464099B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010406190.2A CN111464099B (zh) 2020-05-14 2020-05-14 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010406190.2A CN111464099B (zh) 2020-05-14 2020-05-14 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法

Publications (2)

Publication Number Publication Date
CN111464099A CN111464099A (zh) 2020-07-28
CN111464099B true CN111464099B (zh) 2021-03-02

Family

ID=71680431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010406190.2A Active CN111464099B (zh) 2020-05-14 2020-05-14 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法

Country Status (1)

Country Link
CN (1) CN111464099B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111682820A (zh) * 2020-08-10 2020-09-18 福州大学 一种单绕组无轴承磁通切换电机直接磁链控制方法及系统
CN112671269B (zh) * 2021-01-07 2022-05-13 中国石油大学(华东) 电磁斥力悬浮装置的悬浮力和转矩的协调控制方法
CN113839592B (zh) * 2021-09-18 2024-02-13 福州大学 时间最优无轴承磁通切换电机转矩及悬浮力预测控制方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4198866B2 (ja) * 2000-05-25 2008-12-17 株式会社荏原製作所 ベアリングレスモータの制御装置
CN101599737B (zh) * 2009-04-22 2011-05-18 南京航空航天大学 空间矢量调制的永磁开关磁链发电机电压控制方法
US10177627B2 (en) * 2015-08-06 2019-01-08 Massachusetts Institute Of Technology Homopolar, flux-biased hysteresis bearingless motor
CN105406784B (zh) * 2015-12-14 2017-11-17 江苏大学 单绕组无轴承电机转矩和悬浮力直接控制器及构造方法

Also Published As

Publication number Publication date
CN111464099A (zh) 2020-07-28

Similar Documents

Publication Publication Date Title
CN111464099B (zh) 单绕组无轴承磁通切换电机低转矩及悬浮力的控制方法
CN111193462B (zh) 单绕组无轴承磁通切换电机转子悬浮旋转直接控制方法
CN111682820A (zh) 一种单绕组无轴承磁通切换电机直接磁链控制方法及系统
CN101615876B (zh) 一种隐极式永磁同步电机的调速控制系统和方法
CN107005194A (zh) 多绕组电动机驱动控制装置
CN105406784B (zh) 单绕组无轴承电机转矩和悬浮力直接控制器及构造方法
CN110034719B (zh) 一种六相单绕组无轴承磁通切换电机缺一相容错控制方法
CN108199639B (zh) 六相和三相双绕组悬浮互补型无轴承磁通电机驱动方法
Hiware et al. Indirect field oriented control for induction motor
CN108512473B (zh) 三相四开关永磁同步电机调速系统的直接转矩控制方法
CN108199640B (zh) 缺相容错型六相和三相双绕组悬浮无轴承磁通电机驱动方法
CN104201965B (zh) 定子永磁型无轴承同步电机转子悬浮控制方法
CN109936320A (zh) 一种基于占空比调制的双电机串联直接转矩控制方法
CN110061676B (zh) 一种基于磁链观测器的无轴承永磁同步电机控制器
CN112968642B (zh) 六相单绕组无轴承磁通切换电机线性自抗扰控制方法
CN205509912U (zh) 单绕组无轴承电机转矩和悬浮力直接控制器
CN106130429A (zh) 无轴承永磁同步电机预测控制器及构造方法
Toliyat et al. Position-sensorless control of surface-mount permanent-magnet AC (PMAC) motors at low speeds
CN113328669A (zh) 高频旋转电流注入无轴承磁通切换电机转子偏心观测方法
Grabner et al. Feedback control of a novel bearingless torque motor using an extended FOC method for PMSMs
Wu et al. Research on Direct Torque Control Based on RZVSVPWM of PMSM
CN109981012B (zh) 六相单绕组无轴承磁通切换电机转子悬浮运行控制方法
CN111585476B (zh) 一种无权重系数的双电机驱动系统预测转矩控制方法
CN113839592B (zh) 时间最优无轴承磁通切换电机转矩及悬浮力预测控制方法
CN112994560B (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