CN116464673A - 一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 - Google Patents
一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 Download PDFInfo
- Publication number
- CN116464673A CN116464673A CN202310419350.0A CN202310419350A CN116464673A CN 116464673 A CN116464673 A CN 116464673A CN 202310419350 A CN202310419350 A CN 202310419350A CN 116464673 A CN116464673 A CN 116464673A
- Authority
- CN
- China
- Prior art keywords
- rotor
- unbalance
- signal
- molecular pump
- identification model
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000000725 suspension Substances 0.000 title claims abstract description 26
- 238000012937 correction Methods 0.000 claims abstract description 58
- 230000004044 response Effects 0.000 claims abstract description 25
- 230000001360 synchronised effect Effects 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 50
- 238000006073 displacement reaction Methods 0.000 claims description 30
- 238000005457 optimization Methods 0.000 claims description 18
- 230000005672 electromagnetic field Effects 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 9
- 238000005339 levitation Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
- F04D29/661—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps
- F04D29/662—Balancing of rotors
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D19/00—Axial-flow pumps
- F04D19/02—Multi-stage pumps
- F04D19/04—Multi-stage pumps specially adapted to the production of a high vacuum, e.g. molecular pumps
- F04D19/048—Multi-stage pumps specially adapted to the production of a high vacuum, e.g. molecular pumps comprising magnetic bearings
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/05—Shafts or bearings, or assemblies thereof, specially adapted for elastic fluid pumps
- F04D29/056—Bearings
- F04D29/058—Bearings magnetic; electromagnetic
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
- F04D29/661—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps
- F04D29/666—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps by means of rotor construction or layout, e.g. unequal distribution of blades or vanes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/30—Compensating imbalance
- G01M1/32—Compensating imbalance by adding material to the body to be tested, e.g. by correcting-weights
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2260/00—Function
- F05D2260/81—Modelling or simulation
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Testing Of Balance (AREA)
Abstract
本发明涉及一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,包括:基于运行中的磁悬浮分子泵构建转子的动力学模型,根据所述动力学模型建立不平衡量识别模型;对所述不平衡量识别模型进行校正优化,获取优化后的不平衡量识别模型;基于所述优化后的不平衡量识别模型识别校正质量,将所述校正质量加入转子的运行中,完成磁悬浮分子泵在线动平衡。本发明能够准确反映转子常规控制模式下的不平衡量,并且有效解决计算过程中同步响应混杂干扰与模型参数误差的问题,从而实现转子不平衡量的准确识别。
Description
技术领域
本发明涉及主动磁轴承控制技术领域,特别是涉及一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法。
背景技术
分子泵是利用高速旋转的涡轮转子,在多级涡轮结构配合下,获得定向速度,将能量传输给气体分子,进而压缩驱向排气口的一种真空泵,广泛应用于空间科学、材料科学、电真空工业、微电子与半导体工业等领域。同机械轴承分子泵相比,磁悬浮分子泵在使用寿命、无油润滑、噪声环保等方面具有明显优势,是分子泵技术发展的必然趋势。
由于材质不均匀以及加工装配中存在的误差,甚至部分转子自身设计上的非对称性等多种因素,导致转子的惯性主轴和几何轴不能完全重合,进而不可避免地引起转子不平衡。即使是已经完成动平衡的分子泵系统,在经过长久运行后,由于转子磨损或叶片腐蚀等问题的影响,残余不平衡质量逐渐积累,造成动平衡效果变差。由于磁悬浮分子泵运行在高速旋转的状态下,逐渐积累的转子质量不平衡若不能被及时检测,便会激发显著的同步振动,造成机械振动,甚至造成转子与接地轴承发生碰撞从而引发严重事故。且过大的系统振动会导致功放长期工作在饱和区,使得转子悬浮精度差,不利于系统的稳定控制。因此,残余不平衡量是磁悬浮旋转机械产品健康检测的一个重要指标,实现转子在实际工作转速状态下的高效不平衡量检测与高精度在线动平衡是分子泵系统应用与维护的重要前提。
转子动平衡方法主要分为离线动平衡和在线动平衡两类。无论是离线还是在线动平衡,基本思路都是通过物理增重或去重,实现转子惯性主轴与几何轴重合。相比于传统的离线动平衡,在线动平衡使用转子系统自身具备的传感器与算法实现,有着无需拆机、动平衡效果更好等优势,因此有着广泛的应用前景。传统在线动平衡方法主要有模态平衡法、影响系数法等,但这些方法依赖多次添加试重来得到理想结果,动平衡效率很低。
与专利CN106153256A相比,本发明的创新点在于:考虑磁悬浮分子泵实际运行条件下的控制方法,构建基于常规控制的不平衡量识别模型。且无需添加试重,仅通过添加旋转电磁场的方式便实现了不平衡量的高精度识别。与专利CN103604563A相比,本发明引入了刚度系数校正优化环节,提高了不平衡量的识别精度,且仅需要添加一次同频旋转电磁场,提高了动平衡效率。
发明内容
本发明的目的是提供一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,不需要对转子运行状态进行限制,无需添加试重,仅需加入一次旋转电磁场实现刚度系数的校正优化,实现转子实际运行条件下,一次启车就能够辨识出不平衡量,具有不平衡量识别精度高,操作效率高和不影响分子泵系统稳定性等优点。
为实现上述目的,本发明提供了如下方案:
一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,包括:
基于运行中的磁悬浮分子泵构建转子的动力学模型,根据所述动力学模型建立不平衡量识别模型;
对所述不平衡量识别模型进行校正优化,获取优化后的不平衡量识别模型;
基于所述优化后的不平衡量识别模型识别校正质量,将所述校正质量加入转子的运行中,完成磁悬浮分子泵在线动平衡。
进一步地,建立所述不平衡量识别模型包括:
基于转子常规控制模式下的运行状态构建包含转子不平衡扰动力的不平衡量识别模型,通过对所述包含转子不平衡扰动力的不平衡量识别模型中的转子不平衡扰动力进行消除,生成所述不平衡量识别模型。
进一步地,所述不平衡扰动力包括:转子静不平衡扰动力和转子动不平衡扰动力,其中,所述转子静不平衡扰动力通过轴承控制消除,所述转子动不平衡扰动力通过双平面校正法消除。
进一步地,所述不平衡量识别模型为:
M=TV
其中,M为校正质量矩阵,T为转换系数矩阵,V为同步响应。
进一步地,对所述不平衡量识别模型进行校正优化,获取优化后的不平衡量识别模型包括:
利用正余弦序列相关运算对所述磁悬浮分子泵的转子进行信号同步响应提取,获取第一转速同频信号;在所述第一转速同频信号中加入旋转电磁矢量场,获取第二转速同频信号;
基于所述第一转速同频信号和所述第二转速同频信号对所述不平衡量识别模型进行刚度系数校正优化,生成所述优化后的不平衡量识别模型。
进一步地,利用正余弦序列相关运算对所述磁悬浮分子泵进行信号同步响应提取包括:
获取转子径向位移信号与轴承控制电流信号,其中,所述转子径向位移信号与所述轴承控制电流信号均包括转速同频信号与其他低高频干扰信号;将所述转子径向位移信号和所述轴承控制电流信号分别通过正余弦序列相关运算,消除所述其他低高频干扰信号,提取第一转速同频信号。
进一步地,将所述转子径向位移信号和所述轴承控制电流信号分别通过正余弦序列相关运算,消除所述其他低高频干扰信号,提取第一转速同频信号包括:
根据周期信号的分解原理,传感器采集的所述转子径向位移信号和所述轴承控制电流信号可以分解为不同幅值初始相位的各频率信号之和,与同频正余弦序列进行求均值:
其中,uin(t)为输入信号序列;ωi为不同频率值,包括转频ω0与其他频率ωi(i≠0);Ai表示频率ωi信号的幅值;θi为频率Ti信号的初始相位;所述第一转速同频信号为uout(t)=A0sin(ω0t+θ0)。
进一步地,生成所述优化后的不平衡量识别模型包括:
通过对比所述第一转速同频信号和所述第二转速同频信号的变化,计算校正优化的转换系数矩阵,基于所述校正优化的转换系数矩阵对所述不平衡量识别模型的刚度系数进行校正优化,生成所述优化后的不平衡量识别模型。
进一步地,计算所述校正优化的转换系数矩阵的数学模型为:
T2(V1-Vv)=Mv
其中,T2为优化的转换系数矩阵,V1与Vv分别为第一转速同频信号矩阵和第二转速同频信号矩阵,Mv为旋转电磁场的等效离心质量矩阵。
进一步地,对所述不平衡量识别模型的刚度系数进行校正优化的计算模型为:
K2=-[A1(Vv-V1),A2(Vv-V1)][B1(Vv-V1)+C1,B2(Vv-V1)+C2]-1
其中,K2为刚度系数矩阵,A1,B1,C1,A2,B2,C2为过程运算矩阵。
本发明的有益效果为:
本发明基于转子常规控制模式下的运行状态进行动力学建模,构建符合转子实际运行条件下的不平衡量识别数学模型,能够准确反映转子实际运行时的不平衡量;本发明在兼顾参数校正精度的基础上,提出一种无需添加试重与转子停机的参数校正方法,设计旋转电磁场,通过对比添加前后转子同步响应矩阵,计算出系统刚度系数的校正优化值,无需转子多次启停,能够降低停机成本且得到较高的参数校正精度,有效解决了计算过程中同步响应混杂干扰与模型参数误差的问题,实现转子不平衡量的准确识别。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法流程图;
图2为本发明实施例的主动磁悬浮转子结构及模型坐标定义图;
图3为本发明实施例的转子不平衡示意图,其中,图(a)为转子静不平衡示意图,图(b)为转子动不平衡示意图;
图4为本发明实施例的转子结构示意图;
图5为本发明实施例的同频提取器结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本实施例提供了一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,如图1所示,包括:
S1.基于常规控制模式下的转子运动状态构建不平衡量识别模型。
基于刚性转子动力学,构建包含转子不平衡扰动力的不平衡量识别模型,其中,转子不平衡扰动力包括转子静不平衡扰动力和动不平衡扰动力,由于惯性主轴与旋转轴不重合产生离心力与离心力矩,利用力的分解原理,如图2所示,将其分解为惯性主轴-几何轴和几何轴-旋转轴偏移产生的离心力和离心力矩,几何轴-旋转轴偏移产生的离心力和离心力矩可以通过轴承控制消除,惯性主轴-几何轴偏移产生的离心力和离心力矩可以通过双平面校正法的校正质量离心力和离心力矩消除;通过双平面校正法得到反映至两配重盘的校正质量,其中不平衡量识别包含利用双平面校正法反映到A、B两配重盘的校正质量,利用校正质量离心力/力矩与静不平衡/动不平衡扰动力抵消;基于内嵌于磁悬浮轴承线圈内侧的转子径向位移传感器和线圈内的电流传感器,获得第一信号。转子结构如图4所示。
S1具体流程如下:
定义几何轴g、惯性主轴w和旋转轴r三个轴,质心由C表示,其垂直于旋转轴r的径向平面定义为U,显然质心C为惯性主轴与平面U的交点。而几何轴g和旋转轴w与平面U的交点分别记为几何中心O和Π。
如图3(a)所示,以旋转中心П为原点构建ПXY坐标系,定义O在惯性固定坐标系ΠXY中的坐标为(x,y)。由于Y坐标相值超前X方向相位π/2,因此可利用复系数法,将Π到O的向量可以表示为R=x+y×j。其中,j为虚数单位。同样的,O到C的距离,即质心偏心距el可表示为:
其中,el与分别代表由于转子质心与几何中心的初始相对距离与角度,其由转子自身静平衡质量特性决定,不随转子转速变化;Ω为转子旋转角速度。
因此,在ΠXY坐标系中,转子质心位置可表示为:
C=R+el
关于质心C的运动方程的特征为:
其中,m代表转子质量,是从∏到C的位移矢量的二阶导数,f是作用在转子上的电磁合力。
结合以上方程得到:
根据欧拉方程,惯性主轴ω相对于几何轴g的倾斜运动为:
其中,egω和分别代表rαβ坐标系下几何轴与惯性主轴偏角,即动态不平衡的初始角模量与初始相位角。
如图3(b)所示,在rαβ坐标系中,转子惯性主轴角度可表示为:
w=γ+egω
利用刚体动力学的欧拉方程,惯性主轴的倾斜运动可通过其相对旋转轴的角位移矢量w可表达为:
其中,Jr和Jz分别代表了转子的横向和极转动惯量,w是从旋转轴r到惯性主轴ω的角位移矢量,P代表作用在转子上的电磁合力矩。
得到具有动态不平衡的等效转动运动方程为:
其中,与/>分别代表旋转轴不平行于几何轴产生的离心力矩和陀螺力矩,JregωΩ2与JzegωΩ2则表示惯性主轴偏移几何轴产生的力矩。
根据以上方程,解得惯性主轴与几何轴的距离与偏角为:
其中,Js满足Js=Jr-Jz。
不平衡量识别磁悬浮电磁力可表示为:
其中,Fa,Fb分别代表转子在AB端所受的轴承合力;ia,ib分别为AB端轴承线圈控制电流;sa,sb分别为转子在AB端轴承线圈处的径向位移;kia,kib分别为AB端电磁力-电流刚度;ksa,ksb分别为AB端电磁力-位移刚度。
径向轴承力计算方程为:
其中,Fax,Fbx,Fay,Fby分别代表四组轴承线圈的径向磁轴承作用力;iax,iay,ibx和iby是代表相应线圈绕组中的电流;sax,say,sbx和sby分别代表转子在AB端线圈轴承平面上相对磁轴承几何中心的转子径向位移。
利用复系数法,各物理量可表达为:
由此,转子所受电磁力产生的合力F与合力矩P的计算方程为:
其中,Lma,Lmb分别代表AB端轴承线圈Ba,Bb距离转子几何中心O的垂直距离。
通过电磁合力与扭矩与几何轴-旋转轴偏移产生的离心力和扭矩构建方程为:
其中,R代表几何轴-旋转轴的偏离量,w代表几何轴-旋转轴的倾角,Jr和Jz代表转子的横向和极转动惯量。
如图4所示,由于轴承空间位置限制,位移传感器无法安装在线圈轴承Ba,Bb处,需要安装在轴承两侧Sa,Sb处实现线圈轴承处位移间接测量。位移传感器测量转子在Sa,Sb平面上的径向位移记为xax,xay,xbx和xby。
其中,Lsa,Lsb分别代表AB端传感器安装平面Sa,Sb距离几何中心O的垂直距离,而Lca,Lcb分别代表AB端配重盘Pa,Pb距离几何中心O的距离,ra和rb分别代表AB端配重盘上校正质量安装位置相对旋转轴心的径向距离。
利用传感器采集的位移信息,可以得到几何轴与旋转轴的径向偏移R与角度偏移γ:
由此得到轴承线圈处转子径向位移计算公式为:
因此得到磁轴承电磁合力与合力矩的计算方程为:
其中,ks=ksa+ksb,Γ=ksbLmb-ksaLma,Λ=Lma 2ksa+Lmb 2ksb。
双平面校正法计算的校正质量可表示为:
其中,mcax,mcay为平衡盘A端所需校正质量在x,y方向的分量,mcbx,mcby为平衡盘B端所需校正质量在x,y方向的分量,j为虚数单元,mca,mcb是平衡盘A、B端的校正质量的复系数表达。
转子不平衡的消除是通过在转子配重盘Pa,pb添加校正质量实现,分别添加特定大小和角度的校正质量mca和mcb,其产生的离心力可表示为:
其中,fca,fcb是校正质量产生的离心力,ra,rb是校正质量安装位置相对转轴的距离,Ω是转子平衡转速。
转子动平衡是通过校正质量离心力,抵消惯性主轴与几何轴偏移造成的转子不平衡,其表现为转子整体所受的作用力和力矩为零,可以表示为:
其中,m为转子质量,Jr为转子极惯性矩,el为几何轴与惯性主轴的偏离矩,egw为几何轴与惯性主轴倾角,Lca,Lcb分别为AB端配重面与质心的距离。
理论上,受到不平衡质量的影响,刚性转子在平衡转速下仅做同步响应。因此按照圆周运动的动力学规律,得到广义位移和倾角的假设等式为:
利用这样的假设,可以直接利用磁轴承系统中的线圈电流和转子位移响应来计算校正质量,而不必求解转子质心的加速度,角速度和角加速度等物理量,大大简化了不平衡求解过程。得到不平衡质量的计算方程为:
为直观表达传感器测量数据与不平衡质量关系,利用矩阵形式表达不平衡质量计算式为:
M=TV
其中,M=[mca,mcb]T代表不平衡质量,V=[xa,xb,ia,ib]T代表转子同步响应,表示校正质量变换系数矩阵。
S2.将转子加速至目标转速,进行同步响应提取,获取第一转速同频信号。
构建利用正余弦序列相关运算实现信号同步响应提取,其中,利用正余弦序列相关运算实现信号同步响应提取包括转子径向位移信号与轴承控制电流信号,转子径向位移信号与轴承控制电流信号均包括转速同频信号与其他低高频干扰信号;将转子径向位移信号和轴承控制电流信号分别通过相关运算,消除低高频干扰信号,从而提取出转速同频信号,即为所述同步响应。同频提取器结构框图如图5所示。
S2具体流程如下:
考虑到输入信号uin(t)中混杂有转速同频量ω0的同步响应与其他频率ωi(i≠0)的干扰信号,按照周期信号的分解原理,可将去零偏后的输入信号分解为:
其中,uin(t)为输入信号序列,ωi为不同频率值,包括转频ω0与其他频率ωi(i≠0),Ai代表频率ωi信号的幅值,θi代表信号初始相位,转速同频量可表示为uout(t)=A0sin(ω0t+θ0)。
为求算A0与θ0,可将同频量分解为正余弦分量:
uout(t)=usinsin(ω0t)+ucoscos(ω0t)
在已构建正余弦基准序列的基础上,只需计算正余弦分量幅值usin与ucos即可。以usin为例,对0~T时间内信号进行运算:
由于转子径向位移与线圈控制电流有限,当信号采集时间足够长时,可以认为T>>Ai。将信号分为转速同频量序列A0sin(ω0t+θ0)sin(ω0t)与干扰频率信号Aisin(ωit+θi)sin(ω0t)。按照三角函数积分规律,各频率信号处理后幅值可分别简化为:
因此可以计算得到usin=A0cos(θ0)。同理可得ucos=A0sin(θ0),从而得到输出信号uout(t)=A0sin(ω0t+θ0),由此实现了同频信号的提取。
S3.在目标转速运行状态下的转子中加入旋转电磁矢量场,获取第二转速同频信号,进行刚度系数和不平衡量识别模型的优化,识别校正质量。
构建基于旋转电磁场的系统刚度系数校正优化流程,其中,基于旋转电磁场的系统刚度系数校正优化包括电流-力刚度系数和位移-力刚度系数;获取实际转子径向位移信号和轴承控制电流信号,基于正余弦序列相关运算实现信号同步响应提取,添加旋转电磁矢量场,对比前后转子同步响应,实现系统刚度系数的数值校正优化。具体流程如下:
添加的旋转电磁场可表示为:
其中,U1x,U1y,U2x,U2y分别代表AB端轴承线圈产生旋转电压相对于霍尔基准角位置在x,y方向上的分量。
为尽可能降低电磁场对转子运行稳定性的影响,设置旋转电压产生的电磁合力为0,即有:
假设轴承线圈的电阻电感值分别为R,L,则旋转电磁场控制电压Uv与其在轴承线圈中的同频电流Iv满足:
Uv=RIv+Liv=(R+LΩj)Iv
为简化表达,用Iv=[Iva;Ivb]代表添加的同频电流,由此产生的AB端轴承同频电磁力fva=kiaIva,fvb=kibIvb。其中fva,fvb分别代表添加在AB端的同频电磁力,Iva,Ivb分别表示需要施加的同频电流。
根据轴承电磁力与质量离心力等效规律,可以得出旋转电磁场等效的离心质量Mv=[Mva;Mvb]的表达式为:
旋转电磁场产生等效离心质量按照实虚部展开写为:
Mvr=[real(Mv),imag(Mv)]。
校正质量变换系数矩阵优化校正过程的在线动平衡流程为:首先将转子加速到平衡转速,获得普通控制模式下的同步响应V1,利用设计参数计算的初始变换系数矩阵T1计算不平衡质量记为M1:
M1=T1V1在控制程序中添加同频控制电压,对比前后转子同步响应的变化,计算校正优化的转换系数矩阵的原理方程为:
T2(V1-Vv)=Mv
其中,T2为优化的转换系数矩阵,V1与Vv是旋转电磁场添加前后的转子同步响应矩阵,Mv为旋转电磁场的等效离心质量矩阵。
具体的,保持转子运行在平衡转速下,添加对应等效离心质量Mvr的旋转电磁场Uv,得到新状态下的同步响应Vv。将同步响应矩阵V1,Vv按实虚部展开为
V1r=[real(V1),imag(V1)],Vvr=[real(Vv),imag(Vv)]。
结合前后两次平衡状态下的转子同步响应,得到状态转换矩阵的计算原理方程为:
T2(V1r-Vvr)=Mvr
为利用矩阵运算实现刚度系数的直接求解,定义刚度系数矩阵K=[ksa,ksb,kia,kib]。
将校正质量变换系数矩阵拆分为两个行矩阵:
T=[Tup;Tdown]
其中,Tup=[Ca,Da,Ea,Fa],Tdown=[Cb,Db,Eb,Fb]。
将刚度系数矩阵提取出来有:
其中,A1,A2分别为1×4矩阵,B1,B2为4×4方阵,其表达式分别为:
类似的,Mvr中的刚度系数矩阵提取出,可得方程:
其中,C1,C2的表达式分别为:
校正质量变换矩阵校正优化原理方程可拆分为上下两个行矩阵,利用矩阵计算规律,分别得到两个变换系数行矩阵的校正方程为:
结合以上方程可以得到:
其中,B1(Vv-V1)+C1与B2(Vv-V1)+C2均为4行2列矩阵。
为简化运算过程,将两项等式重新结合,分别为:
4×4方阵KB=[B1(Vv-V1)+C1,B2(Vv-V1)+C2];
1×4矩阵KA=[A1(Vv-V1),A2(Vv-V1)]。
由此,得到刚度系数矩阵的计算方程为:
K2=-[A1(Vv-V1),A2(Vv-V1)][B1(Vv-V1)+C1,B2(Vv-V1)+C2]-1
其中,K2为刚度系数矩阵,A1,B1,C1为过程运算矩阵。
简化为:
通过计算得到转子平衡转速条件下的实际刚度系数数值,从而得到校正后质量变换系数矩阵T2,从而得到实际校正质量计算方程为:
M2=T2V1。
S4.停止转子运行,将获取的校正质量添加至转子中,判断残余不平衡量是否满足平衡精度等级,若满足,则结束平衡;若不满足,则再次进行步骤S2~S4,直到残余不平衡量满足平衡精度等级为止。
本发明基于转子常规控制模式下的运行状态进行动力学建模,构建符合转子实际运行条件下的不平衡量识别数学模型,能够准确反映转子实际运行时的不平衡量;本发明在兼顾参数校正精度的基础上,提出一种无需添加试重与转子停机的参数校正方法,设计旋转电磁场,通过对比添加前后转子同步响应矩阵,计算出系统刚度系数的校正优化值,无需转子多次启停,能够降低停机成本且得到较高的参数校正精度,有效解决了计算过程中同步响应混杂干扰与模型参数误差的问题,实现转子不平衡量的准确识别。
以上所述的实施例仅是对本发明优选方式进行的描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (10)
1.一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,包括:
基于运行中的磁悬浮分子泵构建转子的动力学模型,根据所述动力学模型建立不平衡量识别模型;
对所述不平衡量识别模型进行校正优化,获取优化后的不平衡量识别模型;
基于所述优化后的不平衡量识别模型识别校正质量,将所述校正质量加入转子的运行中,完成磁悬浮分子泵在线动平衡。
2.根据权利要求1所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,建立所述不平衡量识别模型包括:
基于转子常规控制模式下的运行状态构建包含转子不平衡扰动力的不平衡量识别模型,通过对所述包含转子不平衡扰动力的不平衡量识别模型中的转子不平衡扰动力进行消除,生成所述不平衡量识别模型。
3.根据权利要求2所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,所述不平衡扰动力包括:转子静不平衡扰动力和转子动不平衡扰动力,其中,所述转子静不平衡扰动力通过轴承控制消除,所述转子动不平衡扰动力通过双平面校正法消除。
4.根据权利要求3所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,所述不平衡量识别模型为:
M=TV
其中,M为校正质量矩阵,T为转换系数矩阵,V为同步响应。
5.根据权利要求1所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,对所述不平衡量识别模型进行校正优化,获取优化后的不平衡量识别模型包括:
利用正余弦序列相关运算对所述磁悬浮分子泵的转子进行信号同步响应提取,获取第一转速同频信号;在所述第一转速同频信号中加入旋转电磁矢量场,获取第二转速同频信号;
基于所述第一转速同频信号和所述第二转速同频信号对所述不平衡量识别模型进行刚度系数校正优化,生成所述优化后的不平衡量识别模型。
6.根据权利要求5所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,利用正余弦序列相关运算对所述磁悬浮分子泵进行信号同步响应提取包括:
获取转子径向位移信号与轴承控制电流信号,其中,所述转子径向位移信号与所述轴承控制电流信号均包括转速同频信号与其他低高频干扰信号;将所述转子径向位移信号和所述轴承控制电流信号分别通过正余弦序列相关运算,消除所述其他低高频干扰信号,提取第一转速同频信号。
7.根据权利要求6所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,将所述转子径向位移信号和所述轴承控制电流信号分别通过正余弦序列相关运算,消除所述其他低高频干扰信号,提取第一转速同频信号包括:
根据周期信号的分解原理,传感器采集的所述转子径向位移信号和所述轴承控制电流信号可以分解为不同幅值初始相位的各频率信号之和,与同频正余弦序列进行求均值:
其中,uin(t)为输入信号序列;ωi为不同频率值,包括转频ω0与其他频率ωi(i≠0);Ai表示频率ωi信号的幅值;θi为频率ωi信号的初始相位;所述第一转速同频信号为uout(t)=A0sin(ω0t+θ0)。
8.根据权利要求5所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,生成所述优化后的不平衡量识别模型包括:
通过对比所述第一转速同频信号和所述第二转速同频信号的变化,计算校正优化的转换系数矩阵,基于所述校正优化的转换系数矩阵对所述不平衡量识别模型的刚度系数进行校正优化,生成所述优化后的不平衡量识别模型。
9.根据权利要求8所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,计算所述校正优化的转换系数矩阵的数学模型为:
T2(V1-Vv)=Mv
其中,T2为优化的转换系数矩阵,V1与Vv分别为第一转速同频信号矩阵和第二转速同频信号矩阵,Mv为旋转电磁场的等效离心质量矩阵。
10.根据权利要求8所述的基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法,其特征在于,对所述不平衡量识别模型的刚度系数进行校正优化的计算模型为:
K2=-[A1(Vv-V1),A2(Vv-V1)][B1(Vv-V1)+C1,B2(Vv-V1)+C2]-1
其中,K2为刚度系数矩阵,A1,B1,C1,A2,B2,C2为过程运算矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310419350.0A CN116464673A (zh) | 2023-04-19 | 2023-04-19 | 一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310419350.0A CN116464673A (zh) | 2023-04-19 | 2023-04-19 | 一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116464673A true CN116464673A (zh) | 2023-07-21 |
Family
ID=87178442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310419350.0A Pending CN116464673A (zh) | 2023-04-19 | 2023-04-19 | 一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116464673A (zh) |
-
2023
- 2023-04-19 CN CN202310419350.0A patent/CN116464673A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108716471B (zh) | 一种磁悬浮分子泵转子极小位移主动控制方法 | |
CN101639395B (zh) | 一种改进的高速主轴全息动平衡方法 | |
CN107133387B (zh) | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 | |
CN104331565B (zh) | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 | |
CN111177943B (zh) | 一种转子不平衡系数变角度迭代搜寻的不平衡补偿方法 | |
CN110261083B (zh) | 一种基于视觉的磁悬浮转子振动力抑制效果测量方法 | |
CN106289776B (zh) | 一种柔性磁悬浮轴承转子刚度阻尼辨识方法 | |
CN110145541A (zh) | 一种基于相位稳定的磁悬浮轴承转子不平衡运动控制方法 | |
CN105478245B (zh) | 基于主轴振动检测的双自由度精密离心机副轴动不平衡量辨识方法 | |
Zheng et al. | Rotor balancing for magnetically levitated TMPs integrated with vibration self-sensing of magnetic bearings | |
CN102778335A (zh) | 一种基于等效初相矢的各向异性转子全息动平衡方法 | |
CN112953344B (zh) | 一种无轴承异步电机转子不平衡振动补偿控制方法 | |
KR101977761B1 (ko) | 회전기기의 불평형 감시방법 및 그 장치 | |
CN116464673A (zh) | 一种基于旋转电磁矢量场的磁悬浮分子泵在线动平衡方法 | |
CN109687799B (zh) | 基于假定旋转坐标的磁轴承转子分层振动补偿方法 | |
WO2023216451A1 (zh) | 一种超重力离心机动平衡长期监测方法 | |
CN109847952B (zh) | 一种基于驱动电流的双轴精密离心机回转台动平衡方法 | |
CN112307580A (zh) | 一种计入动态运行刚度的高精度在线智能动平衡方法 | |
CN103217284B (zh) | 一种旋转机械模态阻尼系数的测量方法 | |
Yu et al. | Unbalanced variation after assembly and double-speed influence coefficient method in the threshing drum | |
CN113188719B (zh) | 一种检控定子固连式磁悬浮旋转关节在线动平衡系统 | |
CN113565874A (zh) | 一种基于变步长最小均方差的磁悬浮传感器干扰抑制方法 | |
CN1103443C (zh) | 非对称转子的全息动平衡方法 | |
CN115169028A (zh) | 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 | |
Yin et al. | Self-Identification and Balancing for AMB-Rotor System of Turbomolecular Pumps |
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 |