CN112968642B - 六相单绕组无轴承磁通切换电机线性自抗扰控制方法 - Google Patents

六相单绕组无轴承磁通切换电机线性自抗扰控制方法 Download PDF

Info

Publication number
CN112968642B
CN112968642B CN202110178182.1A CN202110178182A CN112968642B CN 112968642 B CN112968642 B CN 112968642B CN 202110178182 A CN202110178182 A CN 202110178182A CN 112968642 B CN112968642 B CN 112968642B
Authority
CN
China
Prior art keywords
rotor
current
suspension
axis
phase
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
CN202110178182.1A
Other languages
English (en)
Other versions
CN112968642A (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 CN202110178182.1A priority Critical patent/CN112968642B/zh
Publication of CN112968642A publication Critical patent/CN112968642A/zh
Application granted granted Critical
Publication of CN112968642B publication Critical patent/CN112968642B/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/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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/05Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for damping motor oscillations, e.g. for reducing hunting
    • 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/22Current control, e.g. using a current control loop
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Connection Of Motors, Electrical Generators, Mechanical Devices, And The Like (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

本发明提出六相单绕组无轴承磁通切换电机线性自抗扰控制方法,电机通过定子绕组的定子电流对定子与转子之间的气隙磁场进行调制;定子绕组为对称的六相绕组,以同方向的定子电流悬浮分量来产生磁悬浮力,使电机的转子径向悬浮;所述转子在悬浮状态时,其旋转过程的转矩与磁悬浮力之间强耦合;所述自抗扰控制方法中,先推导出转子位移与定子电流之间的关系,并构造成利于控制器设计的标准模型,设计含模型扰动信息的线性扩张状态观测器及线性反馈控制律来控制定子电流,从而提升转子启动悬浮过程的稳态性能和动态性能,降低转子悬浮时的径向位移和转速脉动;本发明能解决电机转子径向位移脉动和转速脉动较大从而影响电机稳定运行的高性能控制问题。

Description

六相单绕组无轴承磁通切换电机线性自抗扰控制方法
技术领域
本发明涉及无轴承磁通切换电机设计技术领域,尤其是六相单绕组无轴承磁通切换电机线性自抗扰控制方法。
背景技术
定子永磁型磁通切换电机永磁体仅放置于定子侧,转子仅由硅钢片堆叠而成,无永磁体和绕组,可以有效避免永磁体因温升导致的退磁风险,具有转矩密度高、工作效率高、转子运行鲁棒性强、适宜高速运行等优点,但由于电机转子采用机械轴承支撑,机械摩擦限制了转子转速的提高,并且轴承润滑会导致污染等问题。
为了克服机械轴承支撑带来的不利问题,将无轴承技术引入该电机中,以磁悬浮方式使转子悬浮,从而构成无轴承磁通切换电机,该方法有利于电机输出转矩能力的充分发挥。
由于无轴承磁通切换电机在转子悬浮、旋转过程中转矩和悬浮力存在强耦合,同时由于转子质量分布不均、电机制造工艺、定子和转子间的不平衡磁拉力等因素,使得转子存在严重的径向位移脉动和转速脉动,目前尚未有很好的控制方法。电机长期工作在较大的振动工况下,易出现机械故障。
发明内容
本发明提出六相单绕组无轴承磁通切换电机线性自抗扰控制方法,能解决电机转子径向位移脉动和转速脉动较大从而影响电机稳定运行的高性能控制问题。
本发明采用以下技术方案。
六相单绕组无轴承磁通切换电机线性自抗扰控制方法,所述电机通过定子绕组的定子电流对定子与转子之间的气隙磁场进行调制;所述定子绕组为对称的六相绕组,以同方向的定子电流悬浮分量来产生磁悬浮力,使电机的转子径向悬浮;所述转子在悬浮状态时,其旋转过程的转矩与磁悬浮力之间强耦合;
所述自抗扰控制方法中,先推导出转子位移与定子电流之间的关系,并构造成利于控制器设计的标准模型,设计含模型扰动信息的线性扩张状态观测器及线性反馈控制律来控制定子电流,从而提升转子启动悬浮过程的稳态性能和动态性能,降低转子悬浮时的径向位移和转速脉动。
所述自抗扰控制方法通过逆变器来控制定子电流,在推导转子位移与定子电流之间的关系,设计线性扩张状态观测器、线性反馈控制律时,采用以下方法;
步骤S1、获取六相定子电流iA~iF,转子电角度θr,转子机械角度θΩ,转子电角频率ωr,转子转速n,转子在X方向的径向位移x和Y方向的径向位移y,将转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure BDA0002940649830000023
定义为状态变量zx2和Y方向的偏心速度
Figure BDA0002940649830000024
定义为状态变量zy2
Figure BDA0002940649830000021
其中,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量,Tx为转子在X方向的总扰动,Ty为转子在Y方向的总扰动;
步骤S2、利用步骤S1得到的转子转速n,将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure BDA0002940649830000022
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
步骤S3、利用步骤S2得到的zn1、zn2和bk,计算交轴转矩电流给定iq *
Figure BDA0002940649830000031
其中,n*为转速给定,kpn为转速控制比例系数;
步骤S4、根据直流母线电压UDC、开关管信号SA~SF,计算α轴悬浮平面电压uαS和β轴悬浮平面电压uβS
Figure BDA0002940649830000032
步骤S5、利用步骤S1得到的iA~iF,步骤S4得到的uαS、uβS计算α轴悬浮平面感应电势EαS和β轴悬浮平面感应电势EβS
Figure BDA0002940649830000033
步骤S6、利用步骤S1得到的θΩ、ωr、x和y计算转子在X方向的已知扰动微分hx和转子在Y方向的已知扰动微分hy
Figure BDA0002940649830000034
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,iq'为交轴转矩电流的微分,a1,a2,a3,a4,h1和h2为扰动补偿系数;
步骤S7、计算转子悬浮控制增益b:
Figure BDA0002940649830000041
其中,lrt为转子质心长度,k1为偏心系数,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流;
步骤S8:将Tx为转子在X方向的总扰动定义为状态变量zx3,将Ty为转子在Y方向的总扰动定义为状态变量zy3,设Tx的微分为hx,Ty的微分为hy
步骤S9:利用步骤1的x、y、zx1、zx2、zy1、zy2,步骤S6的hx、hy,步骤S7的b,得到悬浮线性扩张状态观测器为:
Figure BDA0002940649830000042
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure BDA0002940649830000043
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
步骤S10:利用步骤S9的
Figure BDA0002940649830000044
步骤S7的b计算ix *和iy *
Figure BDA0002940649830000051
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
步骤S11:利用步骤S1得到的θr,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure BDA0002940649830000052
其中,id *=0;
步骤S12:利用步骤S1的θr,步骤S10的ix *和iy *,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure BDA0002940649830000053
其中,
Figure BDA0002940649830000054
步骤S13:利用步骤S11的iαT *、iβT *,步骤S12的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure BDA0002940649830000055
其中,零序电流给定io1 *=0,io2 *=0;
步骤S14:利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制。
步骤S6中,扰动补偿系数可通过有限元仿真数据计算得到,具体包括以下步骤:
步骤S61:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中径向位移耦合电感Le',dSqS悬浮平面磁链中永磁偏差磁链Fme,dSqS悬浮平面磁链中永磁偏差磁链Fme',径向位移耦合磁链Fle,转子在X方向单边磁拉力函数fx,转子在Y方向的单边磁拉力函数fy,dS轴悬浮平面磁链ψx和qS轴悬浮平面磁链ψy计算扰动补偿系数:
Figure BDA0002940649830000061
其中,
Figure BDA0002940649830000062
步骤S7中,偏心系数k1可通过有限元仿真数据计算得到,具体包括以下步骤:
步骤S71:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中永磁偏差磁链Fme,绕组自电感L0,计算偏心系数k1
Figure BDA0002940649830000071
步骤S14中,六相逆变桥臂开关状态量SA~SF计算可以采用电流滞环控制策略得到:
Figure BDA0002940649830000072
时,Si=1(i=A~F)  公式十七;
当时
Figure BDA0002940649830000073
Si=0(i=A~F)  公式十八;
式中,ε为允许设定的电流控制误差。
所述电机的转子包括10个齿结构;电机的定子包括10个U型铁心;每个U型铁心之间夹有一个沿切向充磁的永磁体,且各永磁体的充磁方向交替相反;电机定子的每相绕组串联绕在空间相互垂直的定子齿上,构成6相对称绕组。
所述6相对称绕组分别为A相绕组、B相绕组、C相绕组、D相绕组、E相绕组、F相绕组;
在电机的转矩控制坐标系下,其中A相和D相绕组空间对称、B相和E相绕组空间对称、C相和F相绕组空间对称;六相绕组轴线空间上互差60°机械角,绕组中同时流过控制电机切向旋转的六相对称转矩电流iAT~iFT和控制转子径向悬浮的六相对称悬浮电流iAS~iFS,同时为保证电机转子产生的悬浮力正比于悬浮电流,令空间对称绕组中流过的悬浮电流大小方向均相等,即iA=iAT+iAS,iB=iBT+iBS,iC=iCT+iCS,iD=iDT+iDS,iE=iET+iES,iF=iFT+iFS;其中iAT=-iDT,iET=-iBT,iCT=-iFT,iAS=iDS,iES=iBS,iCS=iFS
所述控制方法为内模控制方法,所述电机的驱动系统硬件结构包括整流电路、滤波电容、直流母线电压采集电路、六相逆变器、无轴承磁通切换电机、六相绕组电流采集电路、电机转子位置角采集电路、转子径向x和y位移采集电路、隔离驱动、中央控制器、人机接口;
所述六相逆变器功率管采用IGBT或MOSFET,中央控制器采用DSP或单片机;
所述六相绕组电流采集电路采用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成;
所述直流母线电压采集电路采用霍尔电压传感器与运算放大器相结合方式构成,或是采用并联电阻分压后接由运算放大器构成的电压跟随器相结合方式构成;
所述转子位置角检测电路采用旋转编码器后接电平转换电路构成,或是采用旋转变压器后接解码电路构成;
所述转子径向x和y位移采集电路采用电涡流传感器后接运算放大器相结合方式构成,或是采用线性光耦后接运算放大器相结合方式构成;
当电机工作时,六相绕组电流采集电路、直流母线电压采集电路、转子径向x和y位移采集电路把采集数据以弱电信号方式输出至中央控制器的A/D转换模块,转子位置角检测电路把采集数据以脉冲信号方式输出至中央控制器的QEP模块,中央控制器根据收到的采集数据和内模控制方法,向逆变器输出逆变桥臂开关信号,经由隔离驱动去控制逆变器中的功率开关管的开关动作。
当电机工作时,绕组中同时流过控制电机转子切向旋转的六相对称转矩电流iAT~iFT,其中iAT=-iDT,iET=-iBT,iCT=-iFT;利用恒功率矩阵将电机A-F自然坐标系定子电流变化到静止直角坐标系αTβT,转矩电流在αTβT坐标系的投影为iαT、iβT,再将转矩电流由静止直角坐标系αTβT变化到dq旋转坐标系中,转矩电流在dq坐标系的投影id、iq;其中,T6恒功率矩阵为:
Figure BDA0002940649830000081
令空间对称绕组中流过的悬浮电流iAS~iFS大小方向均相等,即iAS=iDS,iES=iBS,iCS=iFS;A相绕组轴线处于转子齿中心线逆时针超前A1线圈9°机械角,根据xy方向磁拉力分析可知转子将受到在空间近似45°方向悬浮力,由此建立悬浮控制坐标系定义。
在建立悬浮控制坐标系定义时,定义转矩控制坐标系的XY坐标系为水平-垂直直角坐标系,使其X轴与A相绕组的A1线圈轴线重合,与A相绕组轴线相差9°机械角,采用转矩控制坐标系定义近似的定义来建立静止直角坐标系αSβS,旋转坐标系dSqS;dS与A相绕组之间的夹角为
Figure BDA0002940649830000091
其中,θr为转子位置角,角度参数
Figure BDA0002940649830000092
kq表示单位正向悬浮电流与单位q轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM表示单位正向悬浮电流与永磁体产生的悬浮力基波幅值。悬浮电流在αSβS坐标系的投影为iαS、iβS;悬浮电流在dSqS坐标系的投影ix、iy
在定义转子旋转的陀螺效应坐标系时,定义静止坐标系x1-y1-z1与转子同步旋转的坐标系i-j-k,其中k轴与转子旋转轴重合,转子底部由调心轴承固定。定义转子绕i、j、k三个轴旋转的转动惯量Ii、Ij、Ik。其中由于转子结构对称,故Ii=Ij
在上述坐标系下,设转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure BDA0002940649830000094
定义为状态变量zx2和Y方向的偏心速度
Figure BDA0002940649830000093
定义为状态变量zy2,扩张状态变量zx3表示转子在X方向的总扰动,扩张状态变量zy3表示转子在Y方向的总扰动,则悬浮系统状态方程可以进一步表示为:
Figure BDA0002940649830000101
其中,b为转子悬浮控制增益,hx为转子在X方向的总扰动Tx的微分,hy为转子在Y方向的总扰动Ty的微分,Y为悬浮状态输出;
所述标准模型包括转子的转子动力学模型、气隙磁场的气隙磁链模型、六相单绕组无轴承磁通切换电机模型标准型;
所述转子的转子动力学模型为
Figure BDA0002940649830000102
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,Ik为k轴的转动惯量,fx为转子在X方向的单边磁拉力,fy为转子在Y方向的单边磁拉力,lrt为转子质心长度,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量;
所述定子与转子之间的气隙磁场的气隙磁链模型为
Figure BDA0002940649830000103
其中
Figure BDA0002940649830000111
其中,k1、k2、kψx和kψy为偏心系数,Le为dSqS悬浮平面磁链中iq电流电感系数,Fme为dSqS悬浮平面磁链中永磁偏差磁链,L0为绕组自电感,ψx为dS轴悬浮平面磁链,ψy为qS轴悬浮平面磁链;
将公式B2代入公式B3;推导出六相单绕组无轴承磁通切换电机模型标准型:
Figure BDA0002940649830000112
其中,fux为转子在X方向的未建模动态,fuy为转子在Y方向的未建模动态,
Figure BDA0002940649830000113
Figure BDA0002940649830000114
由公式B6计算hx与hy为:
Figure BDA0002940649830000115
其中
Figure BDA0002940649830000116
其中,Le'为dSqS悬浮平面磁链中径向位移耦合电感,Fme'为dSqS悬浮平面磁链中永磁偏差磁链,Fle为径向位移耦合;
针对公式B1建立悬浮线性扩张状态观测器;
Figure BDA0002940649830000121
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure BDA0002940649830000122
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
根据公式B9得到的悬浮线性扩张状态观测器的输出
Figure BDA0002940649830000123
计算ix *和iy *
Figure BDA0002940649830000124
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure BDA0002940649830000125
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
由公式B11转速线性扩张状态观测器的输出
Figure BDA0002940649830000126
计算交轴转矩电流给定iq *
Figure BDA0002940649830000131
其中,n*为转速给定,kpn为转速控制比例系数;
根据公式B12,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure BDA0002940649830000132
根据公式B10,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure BDA0002940649830000133
其中
Figure BDA0002940649830000134
根据公式B13的iαT *、iβT *,式(14)的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure BDA0002940649830000135
其中,零序电流给定io1 *=0,io2 *=0。
利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制。
与现有技术相比,本发明有以下有益效果:
1、本发明得到了传统方法中需要在线辨识或经验整定的控制器增益的数学表达式,采用了含有已知模型扰动信息的线性扩张状态观测器及线性反馈控制律,有效增强了转子悬浮、旋转控制中电机抗干扰能力,降低了转子径向位移脉动和转速脉动,提高了转子悬浮运行的稳定性。
2、本发明在线性自抗扰控制方法设计中,考虑了转子悬浮运行的转子动力学模型和磁链模型,允许转子稳定悬浮运行的转速范围更广阔。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
附图1为本发明实施例的方法原理示意图;
附图2为本发明实施例的六相单绕组无轴承磁通切换电机横截面示意图;
附图3为本发明实施例的驱动系统硬件结构示意图;
附图4为本发明实施例的转矩控制坐标系定义示意图;
附图5为本发明实施例的悬浮控制坐标系定义示意图;
附图6为本发明实施例的转子动态偏心运动坐标系定义示意图。
具体实施方式
如图所示,六相单绕组无轴承磁通切换电机线性自抗扰控制方法,所述电机通过定子绕组的定子电流对定子与转子之间的气隙磁场进行调制;所述定子绕组为对称的六相绕组,以同方向的定子电流悬浮分量来产生磁悬浮力,使电机的转子径向悬浮;所述转子在悬浮状态时,其旋转过程的转矩与磁悬浮力之间强耦合;
所述自抗扰控制方法中,先推导出转子位移与定子电流之间的关系,并构造成利于控制器设计的标准模型,设计含模型扰动信息的线性扩张状态观测器及线性反馈控制律来控制定子电流,从而提升转子启动悬浮过程的稳态性能和动态性能,降低转子悬浮时的径向位移和转速脉动。
所述自抗扰控制方法通过逆变器来控制定子电流,在推导转子位移与定子电流之间的关系,设计线性扩张状态观测器、线性反馈控制律时,采用以下方法;
步骤S1、获取六相定子电流iA~iF,转子电角度θr,转子机械角度θΩ,转子电角频率ωr,转子转速n,转子在X方向的径向位移x和Y方向的径向位移y,将转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure BDA0002940649830000151
定义为状态变量zx2和Y方向的偏心速度
Figure BDA0002940649830000152
定义为状态变量zy2
Figure BDA0002940649830000153
其中,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量,Tx为转子在X方向的总扰动,Ty为转子在Y方向的总扰动;
步骤S2、利用步骤S1得到的转子转速n,将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure BDA0002940649830000154
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
步骤S3、利用步骤S2得到的zn1、zn2和bk,计算交轴转矩电流给定iq *
Figure BDA0002940649830000155
其中,n*为转速给定,kpn为转速控制比例系数;
步骤S4、根据直流母线电压UDC、开关管信号SA~SF,计算α轴悬浮平面电压uαS和β轴悬浮平面电压uβS
Figure BDA0002940649830000161
步骤S5、利用步骤S1得到的iA~iF,步骤S4得到的uαS、uβS计算α轴悬浮平面感应电势EαS和β轴悬浮平面感应电势EβS
Figure BDA0002940649830000162
步骤S6、利用步骤S1得到的θΩ、ωr、x和y计算转子在X方向的已知扰动微分hx和转子在Y方向的已知扰动微分hy
Figure BDA0002940649830000163
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,iq'为交轴转矩电流的微分,a1,a2,a3,a4,h1和h2为扰动补偿系数;
步骤S7、计算转子悬浮控制增益b:
Figure BDA0002940649830000164
其中,lrt为转子质心长度,k1为偏心系数,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流;
步骤S8:将Tx为转子在X方向的总扰动定义为状态变量zx3,将Ty为转子在Y方向的总扰动定义为状态变量zy3,设Tx的微分为hx,Ty的微分为hy
步骤S9:利用步骤1的x、y、zx1、zx2、zy1、zy2,步骤S6的hx、hy,步骤S7的b,得到悬浮线性扩张状态观测器为:
Figure BDA0002940649830000171
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure BDA0002940649830000172
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
步骤S10:利用步骤S9的
Figure BDA0002940649830000173
步骤S7的b计算ix *和iy *
Figure BDA0002940649830000174
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
步骤S11:利用步骤S1得到的θr,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure BDA0002940649830000181
其中,id *=0;
步骤S12:利用步骤S1的θr,步骤S10的ix *和iy *,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure BDA0002940649830000182
其中,
Figure BDA0002940649830000183
步骤S13:利用步骤S11的iαT *、iβT *,步骤S12的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure BDA0002940649830000184
其中,零序电流给定io1 *=0,io2 *=0;
步骤S14:利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制。
步骤S6中,扰动补偿系数可通过有限元仿真数据计算得到,具体包括以下步骤:
步骤S61:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中径向位移耦合电感Le',dSqS悬浮平面磁链中永磁偏差磁链Fme,dSqS悬浮平面磁链中永磁偏差磁链Fme',径向位移耦合磁链Fle,转子在X方向单边磁拉力函数fx,转子在Y方向的单边磁拉力函数fy,dS轴悬浮平面磁链ψx和qS轴悬浮平面磁链ψy计算扰动补偿系数:
Figure BDA0002940649830000191
其中,
Figure BDA0002940649830000192
步骤S7中,偏心系数k1可通过有限元仿真数据计算得到,具体包括以下步骤:
步骤S71:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中永磁偏差磁链Fme,绕组自电感L0,计算偏心系数k1
Figure BDA0002940649830000193
步骤S14中,六相逆变桥臂开关状态量SA~SF计算可以采用电流滞环控制策略得到:
Figure BDA0002940649830000194
时,Si=1(i=A~F)  公式十七;
当时
Figure BDA0002940649830000201
Si=0(i=A~F)  公式十八;
式中,ε为允许设定的电流控制误差。
如图2所示,所述电机的转子包括10个齿结构;电机的定子包括10个U型铁心;每个U型铁心之间夹有一个沿切向充磁的永磁体,且各永磁体的充磁方向交替相反;电机定子的每相绕组串联绕在空间相互垂直的定子齿上,构成6相对称绕组。
所述6相对称绕组分别为A相绕组、B相绕组、C相绕组、D相绕组、E相绕组、F相绕组;
本例中,图4为转矩控制坐标系定义;在电机的转矩控制坐标系下,其中A相和D相绕组空间对称、B相和E相绕组空间对称、C相和F相绕组空间对称;六相绕组轴线空间上互差60°机械角,绕组中同时流过控制电机切向旋转的六相对称转矩电流iAT~iFT和控制转子径向悬浮的六相对称悬浮电流iAS~iFS,同时为保证电机转子产生的悬浮力正比于悬浮电流,令空间对称绕组中流过的悬浮电流大小方向均相等,即iA=iAT+iAS,iB=iBT+iBS,iC=iCT+iCS,iD=iDT+iDS,iE=iET+iES,iF=iFT+iFS;其中iAT=-iDT,iET=-iBT,iCT=-iFT,iAS=iDS,iES=iBS,iCS=iFS
所述控制方法为内模控制方法,如图3所示,所述电机的驱动系统硬件结构包括整流电路、滤波电容、直流母线电压采集电路、六相逆变器、无轴承磁通切换电机、六相绕组电流采集电路、电机转子位置角采集电路、转子径向x和y位移采集电路、隔离驱动、中央控制器、人机接口;
所述六相逆变器功率管采用IGBT或MOSFET,中央控制器采用DSP或单片机;
所述六相绕组电流采集电路采用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成;
所述直流母线电压采集电路采用霍尔电压传感器与运算放大器相结合方式构成,或是采用并联电阻分压后接由运算放大器构成的电压跟随器相结合方式构成;
所述转子位置角检测电路采用旋转编码器后接电平转换电路构成,或是采用旋转变压器后接解码电路构成;
所述转子径向x和y位移采集电路采用电涡流传感器后接运算放大器相结合方式构成,或是采用线性光耦后接运算放大器相结合方式构成;
当电机工作时,六相绕组电流采集电路、直流母线电压采集电路、转子径向x和y位移采集电路把采集数据以弱电信号方式输出至中央控制器的A/D转换模块,转子位置角检测电路把采集数据以脉冲信号方式输出至中央控制器的QEP模块,中央控制器根据收到的采集数据和内模控制方法,向逆变器输出逆变桥臂开关信号,经由隔离驱动去控制逆变器中的功率开关管的开关动作。
当电机工作时,绕组中同时流过控制电机转子切向旋转的六相对称转矩电流iAT~iFT,其中iAT=-iDT,iET=-iBT,iCT=-iFT;利用恒功率矩阵将电机A-F自然坐标系定子电流变化到静止直角坐标系αTβT,转矩电流在αTβT坐标系的投影为iαT、iβT,再将转矩电流由静止直角坐标系αTβT变化到dq旋转坐标系中,转矩电流在dq坐标系的投影id、iq;其中,T6恒功率矩阵为:
Figure BDA0002940649830000211
令空间对称绕组中流过的悬浮电流iAS~iFS大小方向均相等,即iAS=iDS,iES=iBS,iCS=iFS;由图1可知,A相绕组轴线处于转子齿中心线逆时针超前A1线圈9°机械角,根据xy方向磁拉力分析可知转子将受到在空间近似45°方向悬浮力,由此建立悬浮控制坐标系定义。
如图5所示,在建立悬浮控制坐标系定义时,定义转矩控制坐标系的XY坐标系为水平-垂直直角坐标系,使其X轴与A相绕组的A1线圈轴线重合,与A相绕组轴线相差9°机械角,采用转矩控制坐标系定义近似的定义来建立静止直角坐标系αSβS,旋转坐标系dSqS;dS与A相绕组之间的夹角为
Figure BDA0002940649830000221
其中,θr为转子位置角,角度参数
Figure BDA0002940649830000222
kq表示单位正向悬浮电流与单位q轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM表示单位正向悬浮电流与永磁体产生的悬浮力基波幅值。悬浮电流在αSβS坐标系的投影为iαS、iβS;悬浮电流在dSqS坐标系的投影ix、iy
如图6为转子旋转的陀螺效应坐标系定义,在定义转子旋转的陀螺效应坐标系时,定义静止坐标系x1-y1-z1与转子同步旋转的坐标系i-j-k,其中k轴与转子旋转轴重合,转子底部由调心轴承固定。定义转子绕i、j、k三个轴旋转的转动惯量Ii、Ij、Ik。其中由于转子结构对称,故Ii=Ij
在上述坐标系下,设转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure BDA0002940649830000223
定义为状态变量zx2和Y方向的偏心速度
Figure BDA0002940649830000224
定义为状态变量zy2,扩张状态变量zx3表示转子在X方向的总扰动,扩张状态变量zy3表示转子在Y方向的总扰动,则悬浮系统状态方程可以进一步表示为:
Figure BDA0002940649830000231
其中,b为转子悬浮控制增益,hx为转子在X方向的总扰动Tx的微分,hy为转子在Y方向的总扰动Ty的微分,Y为悬浮状态输出;
所述标准模型包括转子的转子动力学模型、气隙磁场的气隙磁链模型、六相单绕组无轴承磁通切换电机模型标准型;
所述转子的转子动力学模型为
Figure BDA0002940649830000232
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,Ik为k轴的转动惯量,fx为转子在X方向的单边磁拉力,fy为转子在Y方向的单边磁拉力,lrt为转子质心长度,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量;
所述定子与转子之间的气隙磁场的气隙磁链模型为
Figure BDA0002940649830000233
其中
Figure BDA0002940649830000241
其中,k1、k2、kψx和kψy为偏心系数,Le为dSqS悬浮平面磁链中iq电流电感系数,Fme为dSqS悬浮平面磁链中永磁偏差磁链,L0为绕组自电感,ψx为dS轴悬浮平面磁链,ψy为qS轴悬浮平面磁链;
将公式B2代入公式B3;推导出六相单绕组无轴承磁通切换电机模型标准型:
Figure BDA0002940649830000242
其中,fux为转子在X方向的未建模动态,fuy为转子在Y方向的未建模动态,
Figure BDA0002940649830000243
Figure BDA0002940649830000244
由公式B6计算hx与hy为:
Figure BDA0002940649830000245
其中
Figure BDA0002940649830000246
其中,Le'为dSqS悬浮平面磁链中径向位移耦合电感,Fme'为dSqS悬浮平面磁链中永磁偏差磁链,Fle为径向位移耦合;
针对公式B1建立悬浮线性扩张状态观测器;
Figure BDA0002940649830000251
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure BDA0002940649830000252
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
根据公式B9得到的悬浮线性扩张状态观测器的输出
Figure BDA0002940649830000253
计算ix *和iy *
Figure BDA0002940649830000254
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure BDA0002940649830000255
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
由公式B11转速线性扩张状态观测器的输出
Figure BDA0002940649830000256
计算交轴转矩电流给定iq *
Figure BDA0002940649830000261
其中,n*为转速给定,kpn为转速控制比例系数;
根据公式B12,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure BDA0002940649830000262
根据公式B10,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure BDA0002940649830000263
其中
Figure BDA0002940649830000264
根据公式B13的iαT *、iβT *,式(14)的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure BDA0002940649830000265
其中,零序电流给定io1 *=0,io2 *=0。
利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。

Claims (7)

1.六相单绕组无轴承磁通切换电机线性自抗扰控制方法,所述电机通过定子绕组的定子电流对定子与转子之间的气隙磁场进行调制;所述定子绕组为对称的六相绕组,以同方向的定子电流悬浮分量来产生磁悬浮力,使电机的转子径向悬浮;所述转子在悬浮状态时,其旋转过程的转矩与磁悬浮力之间强耦合;其特征在于:所述自抗扰控制方法中,先推导出转子位移与定子电流之间的关系,并构造成利于控制器设计的标准模型,设计含模型扰动信息的线性扩张状态观测器及线性反馈控制律来控制定子电流,从而提升转子启动悬浮过程的稳态性能和动态性能,降低转子悬浮时的径向位移和转速脉动;
所述自抗扰控制方法通过逆变器来控制定子电流,在推导转子位移与定子电流之间的关系,设计线性扩张状态观测器、线性反馈控制律时,采用以下方法;
步骤S1、获取六相定子电流iA~iF,转子电角度θr,转子机械角度θΩ,转子电角频率ωr,转子转速n,转子在X方向的径向位移x和Y方向的径向位移y,将转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure FDA0004009634650000011
定义为状态变量zx2和Y方向的偏心速度
Figure FDA0004009634650000012
定义为状态变量zy2
Figure FDA0004009634650000013
其中,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量,Tx为转子在X方向的总扰动,Ty为转子在Y方向的总扰动;
步骤S2、利用步骤S1得到的转子转速n,将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure FDA0004009634650000021
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
步骤S3、利用步骤S2得到的zn1、zn2和bk,计算交轴转矩电流给定iq *
Figure FDA0004009634650000022
其中,n*为转速给定,kpn为转速控制比例系数;
步骤S4、根据直流母线电压UDC、开关管信号SA~SF,计算α轴悬浮平面电压uαS和β轴悬浮平面电压uβS
Figure FDA0004009634650000023
步骤S5、利用步骤S1得到的iA~iF,步骤S4得到的uαS、uβS计算α轴悬浮平面感应电势EαS和β轴悬浮平面感应电势EβS
Figure FDA0004009634650000024
步骤S6、利用步骤S1得到的θΩ、ωr、x和y计算转子在X方向的已知扰动微分hx和转子在Y方向的已知扰动微分hy
Figure FDA0004009634650000031
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,
Figure FDA0004009634650000032
为交轴转矩电流的微分,a1,a2,a3,a4,h1和h2为扰动补偿系数;
步骤S7、计算转子悬浮控制增益b:
Figure FDA0004009634650000033
其中,lrt为转子质心长度,k1为偏心系数,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流;
步骤S8:将Tx为转子在X方向的总扰动定义为状态变量zx3,将Ty为转子在Y方向的总扰动定义为状态变量zy3,设Tx的微分为hx,Ty的微分为hy
步骤S9:利用步骤1的x、y、zx1、zx2、zy1、zy2,步骤S6的hx、hy,步骤S7的b,得到悬浮线性扩张状态观测器为:
Figure FDA0004009634650000034
Figure FDA0004009634650000041
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure FDA0004009634650000042
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
步骤S10:利用步骤S9的
Figure FDA0004009634650000043
步骤S7的b计算ix *和iy *
Figure FDA0004009634650000044
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
步骤S11:利用步骤S1得到的θr,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure FDA0004009634650000045
其中,id *=0;
步骤S12:利用步骤S1的θr,步骤S10的ix *和iy *,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure FDA0004009634650000046
其中,
Figure FDA0004009634650000047
步骤S13:利用步骤S11的iαT *、iβT *,步骤S12的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure FDA0004009634650000051
其中,零序电流给定io1 *=0,io2 *=0;
步骤S14:利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制;
步骤S6中,扰动补偿系数可通过有限元仿真数据计算得到,具体包括以下步骤:
步骤S61:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中径向位移耦合电感Le',dSqS悬浮平面磁链中永磁偏差磁链Fme,dSqS悬浮平面磁链中永磁偏差磁链Fme',径向位移耦合磁链Fle,转子在X方向单边磁拉力函数fx,转子在Y方向的单边磁拉力函数fy,dS轴悬浮平面磁链ψx和qS轴悬浮平面磁链ψy计算扰动补偿系数:
Figure FDA0004009634650000052
其中,
Figure FDA0004009634650000053
步骤S7中,偏心系数k1可通过有限元仿真数据计算得到,具体包括以下步骤:步骤S71:得到dSqS悬浮平面磁链中iq电流电感Le系数,dSqS悬浮平面磁链中永磁偏差磁链Fme,绕组自电感L0,计算偏心系数k1
Figure FDA0004009634650000061
2.根据权利要求1所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:步骤S14中,六相逆变桥臂开关状态量SA~SF计算可以采用电流滞环控制策略得到:
Figure FDA0004009634650000062
时,Si=1(i=A~F)    公式十七;
当时
Figure FDA0004009634650000063
Si=0(i=A~F)    公式十八;
式中,ε为允许设定的电流控制误差。
3.根据权利要求1所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:所述电机的转子包括10个齿结构;电机的定子包括10个U型铁心;每个U型铁心之间夹有一个沿切向充磁的永磁体,且各永磁体的充磁方向交替相反;电机定子的每相绕组串联绕在空间相互垂直的定子齿上,构成6相对称绕组。
4.根据权利要求3所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:所述6相对称绕组分别为A相绕组、B相绕组、C相绕组、D相绕组、E相绕组、F相绕组;
在电机的转矩控制坐标系下,其中A相和D相绕组空间对称、B相和E相绕组空间对称、C相和F相绕组空间对称;六相绕组轴线空间上互差60°机械角,绕组中同时流过控制电机切向旋转的六相对称转矩电流iAT~iFT和控制转子径向悬浮的六相对称悬浮电流iAS~iFS,同时为保证电机转子产生的悬浮力正比于悬浮电流,令空间对称绕组中流过的悬浮电流大小方向均相等,即iA=iAT+iAS,iB=iBT+iBS,iC=iCT+iCS,iD=iDT+iDS,iE=iET+iES,iF=iFT+iFS;其中iAT=-iDT,iET=-iBT,iCT=-iFT,iAS=iDS,iES=iBS,iCS=iFS
5.根据权利要求4所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:所述控制方法为内模控制方法,所述电机的驱动系统硬件结构包括整流电路、滤波电容、直流母线电压采集电路、六相逆变器、无轴承磁通切换电机、六相绕组电流采集电路、电机转子位置角采集电路、转子径向x和y位移采集电路、隔离驱动、中央控制器、人机接口;
所述六相逆变器功率管采用IGBT或MOSFET,中央控制器采用DSP或单片机;
所述六相绕组电流采集电路采用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成;
所述直流母线电压采集电路采用霍尔电压传感器与运算放大器相结合方式构成,或是采用并联电阻分压后接由运算放大器构成的电压跟随器相结合方式构成;
所述转子位置角检测电路采用旋转编码器后接电平转换电路构成,或是采用旋转变压器后接解码电路构成;
所述转子径向x和y位移采集电路采用电涡流传感器后接运算放大器相结合方式构成,或是采用线性光耦后接运算放大器相结合方式构成;
当电机工作时,六相绕组电流采集电路、直流母线电压采集电路、转子径向x和y位移采集电路把采集数据以弱电信号方式输出至中央控制器的A/D转换模块,转子位置角检测电路把采集数据以脉冲信号方式输出至中央控制器的QEP模块,中央控制器根据收到的采集数据和内模控制方法,向逆变器输出逆变桥臂开关信号,经由隔离驱动去控制逆变器中的功率开关管的开关动作。
6.根据权利要求5所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:当电机工作时,绕组中同时流过控制电机转子切向旋转的六相对称转矩电流iAT~iFT,其中iAT=-iDT,iET=-iBT,iCT=-iFT;利用恒功率矩阵将电机A-F自然坐标系定子电流变化到静止直角坐标系αTβT,转矩电流在αTβT坐标系的投影为iαT、iβT,再将转矩电流由静止直角坐标系αTβT变化到dq旋转坐标系中,转矩电流在dq坐标系的投影id、iq;其中,T6恒功率矩阵为:
Figure FDA0004009634650000081
令空间对称绕组中流过的悬浮电流iAS~iFS大小方向均相等,即iAS=iDS,iES=iBS,iCS=iFS;A相绕组轴线处于转子齿中心线逆时针超前A1线圈9°机械角,根据xy方向磁拉力分析可知转子将受到在空间近似45°方向悬浮力,由此建立悬浮控制坐标系定义。
7.根据权利要求5所述的六相单绕组无轴承磁通切换电机线性自抗扰控制方法,其特征在于:在建立悬浮控制坐标系定义时,定义转矩控制坐标系的XY坐标系为水平-垂直直角坐标系,使其X轴与A相绕组的A1线圈轴线重合,与A相绕组轴线相差9°机械角,采用转矩控制坐标系定义近似的定义来建立静止直角坐标系αSβS,旋转坐标系dSqS;dS与A相绕组之间的夹角为
Figure FDA0004009634650000082
其中,θr为转子位置角,角度参数
Figure FDA0004009634650000083
kq表示单位正向悬浮电流与单位q轴电流气隙磁场相互作用产生的悬浮力基波幅值,kPM表示单位正向悬浮电流与永磁体产生的悬浮力基波幅值;悬浮电流在αSβS坐标系的投影为iαS、iβS;悬浮电流在dSqS坐标系的投影ix、iy
在定义转子旋转的陀螺效应坐标系时,定义静止坐标系x1-y1-z1与转子同步旋转的坐标系i-j-k,其中k轴与转子旋转轴重合,转子底部由调心轴承固定;定义转子绕i、j、k三个轴旋转的转动惯量Ii、Ij、Ik;其中由于转子结构对称,故Ii=Ij
在上述坐标系下,设转子的径向位移x定义为状态变量zx1,转子的径向位移y定义为状态变量zy1,转子在X方向的偏心速度
Figure FDA0004009634650000091
定义为状态变量zx2和Y方向的偏心速度
Figure FDA0004009634650000092
定义为状态变量zy2,扩张状态变量zx3表示转子在X方向的总扰动,扩张状态变量zy3表示转子在Y方向的总扰动,则悬浮系统状态方程可以进一步表示为:
Figure FDA0004009634650000093
其中,b为转子悬浮控制增益,hx为转子在X方向的总扰动Tx的微分,hy为转子在Y方向的总扰动Ty的微分,Y为悬浮状态输出;
所述标准模型包括转子的转子动力学模型、气隙磁场的气隙磁链模型、六相单绕组无轴承磁通切换电机模型标准型;
所述转子的转子动力学模型为
Figure FDA0004009634650000094
其中,ga为重力加速度,m为转子质量,h为转子重心高度,Ii为i轴的转动惯量,Ik为k轴的转动惯量,fx为转子在X方向的单边磁拉力,fy为转子在Y方向的单边磁拉力,lrt为转子质心长度,kPM为单位悬浮电流分量与永磁磁动势相互作用产生的悬浮力基波分量幅值,kq是单位悬浮电流分量与单位iq电流相互作用产生的悬浮力基波分量幅值,iq为交轴转矩电流,ix为转子在x方向悬浮的dS轴悬浮电流分量,iy为转子在y方向悬浮的qS轴悬浮电流分量;
所述定子与转子之间的气隙磁场的气隙磁链模型为
Figure FDA0004009634650000101
其中
Figure FDA0004009634650000102
其中,k1、k2、kψx和kψy为偏心系数,Le为dSqS悬浮平面磁链中iq电流电感系数,Fme为dSqS悬浮平面磁链中永磁偏差磁链,L0为绕组自电感,ψx为dS轴悬浮平面磁链,ψy为qS轴悬浮平面磁链;
将公式B2代入公式B3;推导出六相单绕组无轴承磁通切换电机模型标准型:
Figure FDA0004009634650000103
其中,fux为转子在X方向的未建模动态,fuy为转子在Y方向的未建模动态,
Figure FDA0004009634650000104
Figure FDA0004009634650000105
由公式B6计算hx与hy为:
Figure FDA0004009634650000106
其中
Figure FDA0004009634650000111
其中,Le'为dSqS悬浮平面磁链中径向位移耦合电感,Fme'为dSqS悬浮平面磁链中永磁偏差磁链,Fle为径向位移耦合;
针对公式B1建立悬浮线性扩张状态观测器;
Figure FDA0004009634650000112
其中,β1、β2和β3为悬浮线性扩张状态观测器带宽,
Figure FDA0004009634650000113
为悬浮线性扩张状态观测器的输出,ix *为转子在x方向悬浮的dS轴悬浮电流分量给定,iy *为转子在y方向悬浮的qS轴悬浮电流分量给定;
根据公式B9得到的悬浮线性扩张状态观测器的输出
Figure FDA0004009634650000114
计算ix *和iy *
Figure FDA0004009634650000115
其中,kp为悬浮控制比例系数,kd为悬浮控制微分系数;
将转子转速n定义为状态变量zn1,将转速总扰动定义为状态变量zn2,则转速线性扩张状态观测器为:
Figure FDA0004009634650000121
其中,λ1和λ2为转速线性扩张状态观测器带宽,bk为转速控制增益,并且转速控制增益bk是k轴转动惯量Jk的倒数,iq *为交轴转矩电流给定;
由公式B11转速线性扩张状态观测器的输出
Figure FDA0004009634650000122
计算交轴转矩电流给定iq *
Figure FDA0004009634650000123
其中,n*为转速给定,kpn为转速控制比例系数;
根据公式B12,通过旋转变换计算α轴转矩平面电流分量给定iαT *和β轴转矩平面电流给定iβT *
Figure FDA0004009634650000124
根据公式B10,通过旋转变换计算α轴悬浮平面电流分量给定iαS *和β轴悬浮平面电流给定iβS *
Figure FDA0004009634650000125
其中
Figure FDA0004009634650000126
根据公式B13的iαT *、iβT *,式(14)的iαS *、iβS *计算六相定子电流给定iA *~iF *
Figure FDA0004009634650000131
其中,零序电流给定io1 *=0,io2 *=0;
利用六相定子电流给定iA *~iF *与定子六相绕组电流iA~iF,得到控制六相逆变桥臂开关状态量SA~SF,其中当对应相桥臂上管导通,下管关断时Si=1,i=A~F;当对应相桥臂上管关断,下管导通时Si=0,i=A~F;上管和下管互补导通,在SA~SF控制作用下,逆变器输出满足要求的定子电流,实现转子径向悬浮控制。
CN202110178182.1A 2021-02-09 2021-02-09 六相单绕组无轴承磁通切换电机线性自抗扰控制方法 Active CN112968642B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110178182.1A CN112968642B (zh) 2021-02-09 2021-02-09 六相单绕组无轴承磁通切换电机线性自抗扰控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110178182.1A CN112968642B (zh) 2021-02-09 2021-02-09 六相单绕组无轴承磁通切换电机线性自抗扰控制方法

Publications (2)

Publication Number Publication Date
CN112968642A CN112968642A (zh) 2021-06-15
CN112968642B true CN112968642B (zh) 2023-04-07

Family

ID=76284522

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110178182.1A Active CN112968642B (zh) 2021-02-09 2021-02-09 六相单绕组无轴承磁通切换电机线性自抗扰控制方法

Country Status (1)

Country Link
CN (1) CN112968642B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115085609B (zh) * 2022-07-27 2022-11-04 南京工程学院 一种单绕组磁悬浮电机滑模系统控制方法
CN115326440B (zh) * 2022-10-14 2022-12-30 西南交通大学 一种永磁电动悬浮实验装置及实验方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105048913A (zh) * 2015-07-31 2015-11-11 河南科技大学 基于电流补偿的无轴承异步电机不平衡振动控制系统
CN108716471A (zh) * 2018-04-03 2018-10-30 北京航空航天大学 一种磁悬浮分子泵转子极小位移主动控制方法
CN109981012A (zh) * 2019-04-25 2019-07-05 福州大学 六相单绕组无轴承磁通切换电机转子悬浮运行内模控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105048913A (zh) * 2015-07-31 2015-11-11 河南科技大学 基于电流补偿的无轴承异步电机不平衡振动控制系统
CN108716471A (zh) * 2018-04-03 2018-10-30 北京航空航天大学 一种磁悬浮分子泵转子极小位移主动控制方法
CN109981012A (zh) * 2019-04-25 2019-07-05 福州大学 六相单绕组无轴承磁通切换电机转子悬浮运行内模控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Active Disturbance Rejection Control for Suspension System of Bearingless Motor with Rectifier Circuit;Chen Li等;《2013 Fifth International Conference on Measuring Technology and Mechatronics Automation》;20130404;第143-146页 *
基于自抗扰控制器的开关磁阻电动机转矩控制系统;司利云等;《微特电机》;20071128(第11期);第35-39页 *
无轴承永磁薄片电机径向扰动前馈补偿控制;丁强等;《电气传动》;20141231(第12期);第58-61页 *

Also Published As

Publication number Publication date
CN112968642A (zh) 2021-06-15

Similar Documents

Publication Publication Date Title
Oshima et al. Characteristics of a permanent magnet type bearingless motor
Huang et al. Analysis and control of multiphase permanent-magnet bearingless motor with a single set of half-coiled winding
CN112968642B (zh) 六相单绕组无轴承磁通切换电机线性自抗扰控制方法
CN110034719B (zh) 一种六相单绕组无轴承磁通切换电机缺一相容错控制方法
CN108199639B (zh) 六相和三相双绕组悬浮互补型无轴承磁通电机驱动方法
CN111193462B (zh) 单绕组无轴承磁通切换电机转子悬浮旋转直接控制方法
CN112019109B (zh) 双绕组无轴承磁通切换电机转子悬浮快速精准控制方法
CN108199640B (zh) 缺相容错型六相和三相双绕组悬浮无轴承磁通电机驱动方法
CN104201965B (zh) 定子永磁型无轴承同步电机转子悬浮控制方法
CN111682820A (zh) 一种单绕组无轴承磁通切换电机直接磁链控制方法及系统
Chen et al. Radial displacement sensorless control of bearingless flux-switching permanent magnet machines based on difference of symmetric-winding flux linkages
CN113328669B (zh) 高频旋转电流注入无轴承磁通切换电机转子偏心观测方法
CN108054973A (zh) 降低悬浮电流的单绕组定子永磁型磁通切换电机驱动方法
CN107040172A (zh) 一种无轴承磁通切换电机转子径向偏移观测方法
Yu et al. Reduction of cross-coupling for bearingless doubly salient electromagnetic motor using nonsynchronous rotating coordinate transformation
CN109067280B (zh) 一种基于空间对称绕组磁链和的无轴承磁通切换电机转子径向位移观测方法
CN108258956B (zh) 一种单绕组宽转子齿结构无轴承开关磁阻发电机控制方法
CN109981012B (zh) 六相单绕组无轴承磁通切换电机转子悬浮运行控制方法
CN113162315B (zh) 一种四自由度磁悬浮开关磁阻电机及共励磁驱动方法
Pei et al. A novel flux-reversal bearingless slice motor controlled by rotor angle independent suspension current
CN105071733B (zh) 无轴承转子径向位移分离控制与综合补偿方法
Ye et al. Development of bearingless induction motors and key technologies
CN107317537B (zh) 六相逆变器供电的无轴承磁通切换电机驱动方法及系统
CN112701984A (zh) 单绕组无轴承磁通切换电机缺相邻两相转子悬浮控制方法
CN113437917B (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