CN101183241A - 一种设计磁悬浮高速转子系统进动交叉参数的方法 - Google Patents
一种设计磁悬浮高速转子系统进动交叉参数的方法 Download PDFInfo
- Publication number
- CN101183241A CN101183241A CNA2007101797228A CN200710179722A CN101183241A CN 101183241 A CN101183241 A CN 101183241A CN A2007101797228 A CNA2007101797228 A CN A2007101797228A CN 200710179722 A CN200710179722 A CN 200710179722A CN 101183241 A CN101183241 A CN 101183241A
- Authority
- CN
- China
- Prior art keywords
- precession
- rotor
- high speed
- designing
- magnetic levitation
- 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.)
- Granted
Links
Images
Landscapes
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
一种设计磁悬浮高速转子系统进动交叉参数的方法,建立磁悬浮闭环转子系统复系数动力学模型,绘制最高速下复系数开环传递函数的负频Nyquist曲线并计算进动交叉的低通截止频率,然后根据进动模态的相角裕度要求确定最高速下的校正目标,搜索设计频率和确定最高速下的进动交叉增益,进而确定整个转速范围内的进动交叉增益,确保转子进动在整个转速范围内的相角稳定裕度,实现了进动交叉反馈的鲁棒稳定设计。本发明基于负频率Nyquist曲线,提出磁悬浮高速转子系统由相角裕度设计进动交叉参数的方法,相对通常的多变量系统设计方法如状态空间分析方法等,不仅非常直观,而且具有良好的鲁棒性,因而更适合应用于离心机、高精度数控车床、透平机、储能飞轮,以及磁悬浮飞轮和磁悬浮控制力矩陀螺等实际的磁悬浮高速转子系统中。
Description
技术领域
本发明涉及一种设计磁悬浮高速转子系统进动交叉参数的方法,可以用于磁悬浮转子系统进动交叉参数的设计。
背景技术
相对于传统的机械滚珠轴承,磁轴承具有无接触、刚度和阻尼主动可控等突出优点,因而没有摩擦和磨损,也无需润滑、允许转子高速旋转、振动小、支承精度高,特别适合于超净环境设备、高转速设备和要求低振动,高精度,长寿命的航天设备。目前,磁悬浮高速转子系统已经在离心机、高精度数控车床、透平机、储能飞轮,以及磁悬浮飞轮和磁悬浮控制力矩陀螺等民用和航天设备中得到日益广泛的应用,并有取代机械轴承的趋势。
但是,高速转子尤其是大惯量比的扁平高速转子具有强烈的陀螺效应,使转子在静态悬浮时解耦的两自由度转动在有转速时发生耦合产生进动和章动,且随转速升高磁悬浮转子趋于不稳定。交叉反馈是一种抑制陀螺效应,提高磁悬浮高速转子进动和章动稳定性的有效方法,对于提升磁悬浮转子的最高稳定转速具有显著效果,其中的关键是要设计合适的交叉参数(主要包括滤波器截止频率和交叉增益)。由于多变量控制理论固有的机理不明晰的缺点,采用多变量理论进行交叉参数设计不能保证稳定裕度,不利于实际应用,使得交叉参数的设计成为一个难点,实际应用时主要通过实验方法确定交叉参数,缺乏理论依据。
目前,双频Bode图开始用于磁悬浮高速转子的稳定性分析,引申的双频Nyquist图已经用于章动交叉参数的设计,并取得良好的设计效果。但是,进动交叉参数设计与章动交叉设计的存在显著区别:(1)章动交叉参数的设计s=jω的正频率域进行,而进动设计是在s=-jω的负频率域进行;(2)章动频率随转速变化而大范围变化,而进动频率在整个转速范围内变化较小;(3)章动交叉参数设计时可以直接选择单位增益临界频率作为设计频率,但由于LPF在负频率域中为超前环节,直接选择单位增益临界频率作为设计频率通常没有设计解。因此,磁悬浮高速转子章动交叉参数的设计步骤不能直接应用于进动交叉参数设计。
发明内容
本发明的技术解决问题:针对进动的特殊性,提供一种基于负频Nyquist曲线的进动交叉参数设计方法,解决了磁悬浮高速转子进动交叉参数设计问题,确保磁悬浮高速转子系统的进动稳定裕度。
本发明的技术解决方案是:建立磁悬浮闭环转子系统复系数动力学模型,绘制最高速下复系数开环传递函数的负频Nyquist曲线并计算进动交叉的低通截止频率,然后根据进动模态的相角裕度要求确定最高速下的校正目标、搜索设计频率和确定最高速下的进动交叉增益,进而确定整个转速范围内的进动交叉增益,确保转子进动在整个转速范围内的相角稳定裕度。
附图说明
图1为本发明的磁悬浮闭环转子系统及坐标系定义;
图2为本发明的控制器框图;
图3为本发明进动交叉校正前后复系数开环传递函数的负频Nyquist曲线(最高速下)。
图4为本发明的流程图;
具体实施方式
以一种磁悬浮转子系统为例说明本发明方法的具体实施方式,设计目标为Frmax=400Hz下进动模态相角裕度达到γp=20°。
磁悬浮闭环转子系统及坐标系定义如图1所示,该系统由位移传感器、控制器、功放、电磁铁和转子构成,oxyz为转子坐标系,o点位于转子质心,x和y轴沿转子径向且只跟随转子的径向转动而不跟随自转,z轴沿转子轴向。α、β为转子径向绕x,y轴转动的角位移,Ω为转子的自转转速。图中只画出转子A、B两端y方向的电磁铁,x方向与之类似。磁悬浮转子系统每个自由度都由位移传感器测量该自由度上的转子位移,如果转子不在给定的零位上,则误差信号通过控制器运算后,由功放输出相应的控制电流,驱动磁轴承电磁铁产生适当的磁力吸引转子回到给定位置上。根据陀螺技术方程建立磁悬浮转子系统的径向转动运动动力学微分方程模型:
上式中α、β表示转子径向相对定子绕x,y轴转动的角位移,Jx=Jy和Jz分别为转子径向和轴向的转动惯量,H=JzΩ为转子角动量,Ω=2πFr为转子转速,pdx和pdy为转子径向的扰动力矩,ki和kh为磁轴承的位移刚度和电流刚度,ks为磁轴承位移传感器灵敏度,lm和ls分别为磁轴承和位移传感器到转子中心的距离,gc、gcr和gw分别为PID控制器、交叉控制器和功放的输入-输出变换算子,即有 L表示拉氏变换,s为算子,gc(s)、gcr(s)和gw(s)分别为PID控制器、进动交叉控制器和功放的传递函数。
本发明方法具体实现时采用的控制器的框图如图2所示,包括PID控制器和进动交叉控制器两部分,其中进动交叉控制器又由LPF环节和增益环节kr11串连构成。转子位移信号输入控制器后,一方面经过PID控制器计算得到一部分控制量,另一方面,x自由度和y自由度的转子位移信号各自求差后输入到进动交叉控制器,计算后分别与y自由度和x自由度PID控制器输出求和或求差,作为总的控制器输出信号馈送到功放。
图3为本发明采用进动交叉校正前后复系数开环传递函数的负频Nyquist曲线(最高速下)。该图坐标系横轴和纵轴分别为实轴和虚轴,曲线上每个点的横坐标和纵坐标分别对应gOL(-jω)值的实部和虚部,ω在(0,+∞)之间取值。图中带点实线代表校正前曲线,带点虚线代表校正后的曲线,频率范围2~700Hz,每Hz一个点,箭头表示频率升高方向。O点为复平面原点,虚线圆表示单位圆。穿越点p1和p2分别为校正前后Nyquist曲线与单位圆的交点,p11为进动交叉校正设计频率点。校正前p1点的频率、相位和相角裕度分别为fp1=12Hz、φp1=176°、γp1=-4°,进动交叉校正设计频率即p11点频率为fp11=20Hz,校正后的p2点频率、相位和相角裕度分别为fp2=20Hz、φp1=160°、γp1=20°。
本发明进动交叉参数设计的流程图如图4所示,具体步骤如下:
(1)建立磁悬浮闭环转子系统复系数动力学模型:由图1和图2根据陀螺技术方程建立磁悬浮转子系统的径向转动运动的动力学微分方程:
用Jrr统一表示Jx=Jy,且令Φ=a+jβ,pd=pdx+jpdy,其中j为虚数单位,将微分方程组的第一式乘以j再加到第二式,再做拉氏变换得到:
令 gceff(s)=2lmlskiksgw(s)[gc(s)+gcr(s)],则系统的复系数动力学模型(这里采用开环传递函数)为:
gOL(s)=goeff(s)gceff(s)
(2)绘制最高速时复系数开环传递函数的负频Nyquist曲线:详见图3的带点虚线。
(3)低通截止频率选取:静态悬浮时的主极点频率为2π×80rad/s,选择进动交叉LPF截止频率等于静态悬浮时主极点频率,即ωr1=2π×80rad/s。
(4)确定最高速下的校正目标:按照设计指标要求γp=20°,则校正目标的相角为180°-γp=160°,即图3中的p2点。
(5)在最高速下搜索设计频率ωp11(=2πfp11):可由下式解出ωp11
其中gcrLPF为交叉低通LPF算子,p2点为p11点期望的校正目标,也就是校正后Nyquist曲线与单位圆的交点。在校正前的负频Nyquist曲线(图3)中查出fp1=12Hz,在ω>ωp1=2π×12rad/s范围内搜索p11点,得到fp11=20Hz,即ωp11=2π×20rad/s,对应校正量=0.3832。
(6)确定最高速下的进动交叉增益:根据
确定最高速下的进动交叉增益为krl1d=0.3048。
(7)确定整个转速范围内的交叉增益:在0<Fr<Frmax范围内,进动交叉增益与转速成正比,即krh1=krh1dFr/Frmax=0.3048×Fr/400,其中Fr和Frmax分别为转子转速和转子最高转速。
本发明中,磁轴承控制器可以是模拟或数字的;磁轴承功放可以是模拟功放或数字功放,可以是线性功放或开关功放,开关功放又可以是滞环(Hysteresis)型、采样保持(Sample/Hold)型、脉宽调制(PWM)型、最小脉宽(MPW)型和三电平型,功放输出功率管的配置形式可以是半桥或全桥;磁轴承电磁铁可以是采用电磁偏置或永磁偏置方式;位移传感器可以是电涡流型传感器或电感型传感器。
Claims (8)
1.一种设计磁悬浮高速转子系统进动交叉参数的方法,其特征在于:建立磁悬浮闭环转子系统复系数动力学模型,绘制最高速下复系数开环传递函数的负频Nyquist曲线并计算进动交叉的低通截止频率,然后根据进动模态的相角裕度要求确定最高速下的校正目标、搜索设计频率和确定最高速下的进动交叉增益,进而确定整个转速范围内的进动交叉增益,确保转子进动在整个转速范围内的相角稳定裕度。
2.根据权利要求1所述的设计磁悬浮高速转子系统进动交叉参数的方法,其特征在于:所述的建立磁悬浮闭环转子系统复系数动力学模型的步骤包括:
(1)根据陀螺技术方程建立磁悬浮转子系统的径向转动运动动力学微分方程模型:
上式中α、β表示转子径向相对定子绕x,y轴转动的角位移,Jx=Jy和Jz分别为转子径向和轴向的转动惯量,H=JzΩ为转子角动量,Ω=2πFr为转子转速,pdx和pdy为转子径向的扰动力矩,ki和kh为磁轴承的位移刚度和电流刚度,ks为磁轴承位移传感器灵敏度,lm和ls分别为磁轴承和位移传感器到转子中心的距离,gc、gcr和gw分别为PID控制器、进动交叉控制器和功放的输入-输出变换算子,即有 L表示拉氏变换,s为算子,gc(s)、gcr(s)和gw(s)分别为PID控制器、进动交叉控制器和功放的传递函数;
(2)采用复数变换将步骤(1)中的实系数两变量方程转化为复系数单变量形式,得到复系数单变量等效系统。用Jrr统一表示Jx=Jy,且令φ=α+jβ,pd=pdx+jpdy,其中j为虚数单位,将微分方程组的第一式乘以j再加到第二式,再做拉氏变换得到:
令 gceff(s)=2lmlskiksgw(s)[gc(s)+gcr(s)],则系统的复系数开环传递函数为:
gOL(s)=goeff(s)gceff(s)。
3.根据权利要求1所述的设计磁悬浮高速转子系统进动交叉参数的方法,其特征在于:所述的最高速时复系数开环传递函数的负频Nyquist曲线,其坐标系横轴和纵轴分别为实轴和虚轴,曲线上每个点的横坐标和纵坐标分别对应gOL(-jω)值的实部和虚部,ω在(0,+∞)之间取值。
4.根据权利要求1所述的设计磁悬浮高速转子系统进动交叉参数的方法,其特征在于:所述的低通截止频率选择为静态悬浮时主极点频率的0.9~1倍。
8.根据权利要求1所述的设计磁悬浮高速转子系统进动交叉参数的方法,其特征在于:所述整个转速范围内的交叉增益在0<Fr<Frmax可以按下述公式确定:
krh1=krh1dFr/Frmax
其中Fr和Frmax分别为转子转速和转子最高转速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710179722A CN100587633C (zh) | 2007-12-17 | 2007-12-17 | 一种设计磁悬浮高速转子系统进动交叉参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710179722A CN100587633C (zh) | 2007-12-17 | 2007-12-17 | 一种设计磁悬浮高速转子系统进动交叉参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101183241A true CN101183241A (zh) | 2008-05-21 |
CN100587633C CN100587633C (zh) | 2010-02-03 |
Family
ID=39448548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710179722A Expired - Fee Related CN100587633C (zh) | 2007-12-17 | 2007-12-17 | 一种设计磁悬浮高速转子系统进动交叉参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100587633C (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101571720B (zh) * | 2009-06-03 | 2011-02-02 | 北京航空航天大学 | 一种磁悬浮飞轮三自由度动量交换控制方法 |
CN104166345A (zh) * | 2014-08-01 | 2014-11-26 | 中国人民解放军装备学院 | 一种磁悬浮控制力矩陀螺转子系统解耦和扰动抑制方法 |
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
CN106200380A (zh) * | 2016-07-15 | 2016-12-07 | 浙江工业大学 | 基于非线性观测的磁悬浮系统跟踪控制方法 |
CN107797450A (zh) * | 2017-10-11 | 2018-03-13 | 北京航空航天大学 | 一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法 |
CN107807533A (zh) * | 2017-12-04 | 2018-03-16 | 南京航空航天大学 | 基于交叉解耦陷波方法的同频振动力矩抑制控制方法 |
CN108828951A (zh) * | 2018-07-20 | 2018-11-16 | 南京航空航天大学 | 磁悬浮轴承多模型鲁棒无扰动切换控制方法 |
CN111289275A (zh) * | 2019-12-31 | 2020-06-16 | 北京化工大学 | 旋转机械稳定性辨识方法、装置、计算机设备及存储介质 |
-
2007
- 2007-12-17 CN CN200710179722A patent/CN100587633C/zh not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101571720B (zh) * | 2009-06-03 | 2011-02-02 | 北京航空航天大学 | 一种磁悬浮飞轮三自由度动量交换控制方法 |
CN104166345A (zh) * | 2014-08-01 | 2014-11-26 | 中国人民解放军装备学院 | 一种磁悬浮控制力矩陀螺转子系统解耦和扰动抑制方法 |
CN104166345B (zh) * | 2014-08-01 | 2019-05-07 | 中国人民解放军装备学院 | 一种磁悬浮控制力矩陀螺转子系统解耦和扰动抑制方法 |
CN104331565A (zh) * | 2014-11-10 | 2015-02-04 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
CN104331565B (zh) * | 2014-11-10 | 2017-07-07 | 河海大学常州校区 | 轴类磁悬浮刚性转子系统的动力学建模方法及控制方法 |
CN106200380A (zh) * | 2016-07-15 | 2016-12-07 | 浙江工业大学 | 基于非线性观测的磁悬浮系统跟踪控制方法 |
CN107797450A (zh) * | 2017-10-11 | 2018-03-13 | 北京航空航天大学 | 一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法 |
CN107807533A (zh) * | 2017-12-04 | 2018-03-16 | 南京航空航天大学 | 基于交叉解耦陷波方法的同频振动力矩抑制控制方法 |
CN108828951A (zh) * | 2018-07-20 | 2018-11-16 | 南京航空航天大学 | 磁悬浮轴承多模型鲁棒无扰动切换控制方法 |
CN111289275A (zh) * | 2019-12-31 | 2020-06-16 | 北京化工大学 | 旋转机械稳定性辨识方法、装置、计算机设备及存储介质 |
CN111289275B (zh) * | 2019-12-31 | 2021-06-01 | 北京化工大学 | 旋转机械稳定性辨识方法、装置、计算机设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN100587633C (zh) | 2010-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100587633C (zh) | 一种设计磁悬浮高速转子系统进动交叉参数的方法 | |
CN100462885C (zh) | 一种磁悬浮反作用飞轮开环高精度不平衡振动控制系统 | |
CN108716471B (zh) | 一种磁悬浮分子泵转子极小位移主动控制方法 | |
Lei et al. | Control of flexible rotor systems with active magnetic bearings | |
CN102506860B (zh) | 一种基于加速度反馈与前馈的惯性稳定装置及其控制方法 | |
CN104950919B (zh) | 一种设计磁悬浮转子自平衡系统自适应滤波器稳定性参数的方法 | |
CN101247097B (zh) | 一种设计磁悬浮扁平高速转子系统陷波器参数的方法 | |
WO2022041298A1 (zh) | 磁悬浮转子系统及其极微振动控制方法和控制装置 | |
CN108897218B (zh) | 一种基于混合奇次重复控制器的磁悬浮转子奇次谐波振动抑制方法 | |
Liu et al. | Research on automatic balance control of active magnetic bearing‐rigid rotor system | |
Fang et al. | Adaptive complete suppression of imbalance vibration in AMB systems using gain phase modifier | |
CN110145541B (zh) | 一种基于相位稳定的磁悬浮轴承转子不平衡运动控制方法 | |
CN102323825B (zh) | 一种航天器机动时dgmscmg系统的力矩补偿控制方法 | |
Zheng et al. | Tracking compensation control for nutation mode of high-speed rotors with strong gyroscopic effects | |
CN108008745A (zh) | 基于相位滞后补偿的磁悬浮飞轮同频振动力抑制方法及控制系统 | |
CN100354869C (zh) | 一种判定磁悬浮转子系统径向转动稳定性的方法 | |
CN108710286A (zh) | 基于交叉解耦陷波器的同频振动力矩分层控制方法 | |
Huang et al. | A review of active magnetic bearing control technology | |
Tang et al. | Active vibration control of the flexible rotor to pass the first bending critical speed in high energy density magnetically suspended motor | |
Tang et al. | High-speed rotor’s mechanical design and stable suspension based on inertia-ratio for gyroscopic effect suppression | |
Su et al. | The precise control of a double gimbal MSCMG based on modal separation and feedback linearization | |
Peng et al. | High precision synchronous vibration suppression for a MSFW subject to phase lag influence | |
CN114322971B (zh) | 基于双二阶广义积分器的磁悬浮转子同频振动力抑制方法 | |
Tang et al. | Active vibration control of the flexible rotor in high energy density magnetically suspended motor with mode separation method | |
CN102102704A (zh) | 五自由度交流主动磁轴承α阶逆系统解耦控制器构造方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100203 Termination date: 20181217 |