CN115169028A - 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 - Google Patents
一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 Download PDFInfo
- Publication number
- CN115169028A CN115169028A CN202210708009.2A CN202210708009A CN115169028A CN 115169028 A CN115169028 A CN 115169028A CN 202210708009 A CN202210708009 A CN 202210708009A CN 115169028 A CN115169028 A CN 115169028A
- Authority
- CN
- China
- Prior art keywords
- rotor
- notch filter
- adaptive notch
- active magnetic
- bearing system
- 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
- 239000000725 suspension Substances 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004146 energy storage Methods 0.000 title abstract description 16
- 230000003044 adaptive effect Effects 0.000 claims abstract description 31
- 238000006073 displacement reaction Methods 0.000 claims abstract description 26
- 238000012546 transfer Methods 0.000 claims abstract description 21
- 230000014509 gene expression Effects 0.000 claims abstract description 15
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000005339 levitation Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16C—SHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
- F16C32/00—Bearings not otherwise provided for
- F16C32/04—Bearings not otherwise provided for using magnetic or electric supporting means
- F16C32/0406—Magnetic bearings
- F16C32/044—Active magnetic bearings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/16—Mechanical energy storage, e.g. flywheels or pressurised fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
本发明涉及一种飞轮储能主动磁悬浮轴承系统不平衡控制方法,在转子周边设置位移传感器,将位移传感器通过带频率估计单元的自适应陷波器与控制器相连,包括以下步骤:建立含不平衡质量的主动磁悬浮轴承系统模型,并推导不平衡力在x、y方向的分量表达式;确定所述自适应陷波器的传递函数,并通过所述自适应陷波器提取所述转子位移信号中的同频振动信号;通过自适应陷波器的输出信号估计所述转子的转速,并确定估计转速的表达式;将经过带频率估计单元的自适应陷波器后得到的不平衡同频量反馈回主动磁悬浮轴承系统中,以消除同频电流,实现转子的自动平衡。本发明在无需外加转速传感器情况下实现了磁悬浮飞轮的自动平衡。
Description
技术领域
本发明涉及磁悬浮转子振动抑制领域,特别是涉及一种飞轮储能主动磁悬浮轴承系统不平衡控制方法。
背景技术
飞轮储能主要是通过高速旋转来储存能量的,如何减小轴承机械摩擦损耗是保证飞轮储能良好性能的关键问题。由于主动磁轴承具有无机械摩擦、无需润滑、电磁力可控等优点,因此在高速主轴、高速电机以及飞轮储能等高速旋转机械领域应用越来越广泛。飞轮转子材质不均、零件加工装配误差等原因会造成转子质量不平衡,该偏心质量引起的激振力会使转子产生与转速同频的振动,这种振动如果不加以消除或减小,磁悬浮转子的控制精度将会受到严重影响,甚至可能使系统损毁。由于主动磁轴承系统工作条件的改变和校正精度的限制,应用传统的现场动平衡技术后转子仍然存在残余不平衡量。因此利用主动磁轴承刚度和阻尼可调的特点,选取合理的算法抑制该不平衡振动是十分必要的。
目前针对不平衡振动的抑制算法已有诸多学者研究,抑制方法主要有两种,第一种是不平衡补偿,主要思想是通过提取转子位移信号中的同频振动信号,计算产生一个与不平衡力大小相等方向相反的力,两者相互抵消,意在实现“零位移控制”,提高转子旋转精度的同时约束转子绕几何轴旋转,目前常用的控制方法主要有影响系数法、基于模型辨识的方法等,但该方法高速下很容易导致功放电路饱和,且反作用力会传递到基座上产生振动,因此该方法适用于对转子位移精度要求高的场合。第二种是自动平衡,主要思想是直接抵消或过滤掉位移传感器输出信号中的同频成分,使控制器对同频振动信号不予响应,使转子绕惯性轴旋转,意在实现“零电流或零轴承力控制”,目前常用的方法有LMS算法、广义陷波器等,这种方法可以避免功放电路饱和,因此它是高速转子振动抑制的最佳选择。
然而,上述方法都需要有转子的角速度反馈才能实现,传统方案中通常由霍尔转速传感器、光电编码器等获得。然而,角速度传感器在许多应用中都存在一些限制,包括机械尺寸有限、可靠性低以及成本限制等问题,且一旦角速度传感器存在较大的误差或发生故障时,控制系统将无法获得准确的速度信号,因此需要对转子的实际速度进行估计。针对磁悬浮轴承转子的角速度估计问题,Lee等提出角位置微分来估计转子转速,但该方法很容易引进噪声。张染等将转子位移信号中的高倍频分量通过Hilbert谱分析得到转速信号,该方法虽然精度有所提高,但倍频分量在低速时不稳定。刘虎等根据已有的位移信号和转速信号构建BP神经网络测速模型,但该方法过度依赖转子的初始位移和转速数据,且计算量较大。
发明内容
本发明提供一种飞轮储能主动磁悬浮轴承系统不平衡控制方法,能够克服现有磁悬浮飞轮储能系统不平衡振动抑制算法的不足。
本发明解决其技术问题所采用的技术方案是:提供一种主动磁悬浮轴承系统不平衡控制方法,在转子周边设置位移传感器,将所述位移传感器通过带频率估计单元的自适应陷波器与控制器相连,包括以下步骤:
(1)建立含不平衡质量的主动磁悬浮轴承系统模型,并推导不平衡力在x、y方向的分量表达式;
(2)确定所述自适应陷波器的传递函数,并通过所述自适应陷波器提取所述转子位移信号中的同频振动信号;
(3)通过自适应陷波器的输出信号估计所述转子的转速,并确定估计转速的表达式;
(4)将经过带频率估计单元的自适应陷波器后得到的不平衡同频量反馈回主动磁悬浮轴承系统中,以消除同频电流,实现转子的自动平衡。
所述步骤(1)中主动磁悬浮轴承系统模型为:其中,xc和yc为转子的几何中心C在固定坐标系O-xy中的坐标,xg和yg为转子的质心G在固定坐标系O-xy 中的坐标,e为转子的偏心距,φ为转子的旋转角度。
所述步骤(2)中所述自适应陷波器的两相传递函数为:令s=jω,当所述自适应陷波器的增益ζ不为零时,则有:其中,D(s)和Q(s)分别表示自适应陷波器的两相传递函数,X1(s)和X2(s)分别为所述自适应陷波器的两相输出信号,s为拉普拉斯算子,Ω为所述自适应陷波器的中心频率,ω表示实际转子转速,ΔΩ为转速的一个极小值。
所述步骤(2)中当估计转速为所述自适应陷波器的中心频率Ω时,输出信号X1(s)将完全跟踪输入的转子的位移信号v,输出信号X2(s)为与输出信号X1(s)幅值绝对值相等的正交信号,此时与转速同频的不平衡振动量被提取出来。
所述主动磁悬浮轴承系统从不平衡扰动dx到控制电流ic的传递函数表达式为:其中, Q(s)=s2+2ζΩs+Ω2,Gs(s)、Gp(s)、C(s)、A(s)分别是位移传感器、转子、控制器和功率放大器的传递函数。
有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明无需外加转速传感器,依靠简单的频率估计单元就能监测飞轮转速,从而实现磁悬浮飞轮的自动平衡,整个方法计算量低、参数设计明确简单。
附图说明
图1是飞轮储能主动磁悬浮轴承系统结构的截面图;
图2是飞轮储能主动磁悬浮轴承系统中不平衡质量示意图;
图3是本发明实施方式的控制框图;
图4是本发明实施方式中自适应陷波器的框图;
图5是本发明实施方式中频率估计单元的框图;
图6是本发明实施方式中控制器的框图;
图7是信噪比10dB条件下补偿前后转子上端的位移轨迹和频率估计轨迹图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明的实施方式涉及一种主动磁悬浮轴承系统不平衡控制方法,该方法可以应用于飞轮储能主动磁悬浮轴承系统中,其结构如图1所示,包括:磁悬浮飞轮转子4,设置在所述磁悬浮飞轮转子4侧边的径向位移传感器1和轴向位移传感器5。径向位移传感器1 和轴向位移传感器5均与磁悬浮控制部分3相连,磁悬浮控制部分3的输出端连接功率放大器2,攻率放大器2向磁悬浮飞轮转子4施加放大后的控制电流。图3所示的是本实施方式中不平衡控制方法的控制框图。本实施方式中使用的磁悬浮控制部分3包括相互连接的带频率估计单元的自适应陷波器和控制器。在进行不平衡控制时,具体包括以下步骤:
步骤1,建立含不平衡质量的飞轮储能主动磁悬浮轴承系统模型,并推导不平衡力在 x、y方向的分量表达式。如图2所示,在飞轮储能主动磁悬浮轴承系统中,由于质量不平衡的存在,转子几何中心C和质心G不重合。转子几何中心C在固定坐标系O-xy中的坐标为(xc,yc),质心G在固定坐标系O-xy中的坐标为(xg,yg),Ω为转子的转速,由此可得转子质心和形心间坐标关系和不平衡力在x、y方向上的分量方程为:
步骤2,确定所述自适应陷波器的传递函数,并通过所述自适应陷波器提取所述飞轮转子位移信号中的同频振动信号。
图4是本实施方式中的自适应陷波器的框图,其两相传递函数为:令s=jω,当所述自适应陷波器的增益ζ不为零时,则有:其中,D(s)和Q(s)分别表示自适应陷波器的两相传递函数,X1(s)和X2(s)分别为所述自适应陷波器的两相输出信号,s为拉普拉斯算子,Ω为所述自适应陷波器的中心频率,ω表示实际转子转速,ΔΩ为转速的一个极小值。
当估计转速为所述自适应陷波器的中心频率Ω时,输出信号X1(s)将完全跟踪输入的飞轮转子的位移信号v,输出信号X2(s)为与输出信号X1(s)幅值绝对值相等的正交信号,此时与转速同频的不平衡振动量被提取出来。
步骤3,通过自适应陷波器的输出信号估计所述飞轮转子的转速,并确定估计转速的表达式。本步骤中通过频率估计单元估计飞轮转子的转速,该频率估计单元的结构如图5所示,其通过估计所述飞轮转子的转速,其中,x1和x2分别表示输出信号X1(s)的值和输出信号X2(s)的值。
从监测转子位移信号v到误差信号e’的传递函数(即带频率估计单元的自适应陷波器的传递函数)表达式为:从不平衡扰动dx到控制电流ic的传递函数表达式为:其中, Q(s)=s2+2ζΩs+Ω2,Gs(s)、Gp(s)、C(s)、A(s)分别是位移传感器、刚性转子、控制器和功率放大器的传递函数。
步骤4,将经过带频率估计单元的自适应陷波器后得到的不平衡同频量反馈回飞轮储能主动磁悬浮轴承系统中,以消除同频电流,实现转子的自动平衡。本步骤中,频率估计单元正确辨识转子频率时,不平衡扰动dx到控制电流ic的传递函数为零,将经过带频率估计单元的自适应陷波器提取得到的同频振动信号反馈回原主动磁轴承飞轮储能系统中去,用以消除同频电流,实现转子的自动平衡。
图7是信噪比10dB条件下补偿前后转子上端的位移轨迹和频率估计轨迹,通过对比可知,采用本实施方式的方法后能够较好的抑制该不平衡振动,且本实施方式的方法无需外加转速传感器,依靠简单的频率估计单元就能监测飞轮转速,从而实现磁悬浮飞轮的自动平衡,整个方法计算量低、参数设计明确简单。
Claims (9)
1.一种主动磁悬浮轴承系统不平衡控制方法,其特征在于,在转子周边设置位移传感器,将所述位移传感器通过带频率估计单元的自适应陷波器与控制器相连,包括以下步骤:
(1)建立含不平衡质量的主动磁悬浮轴承系统模型,并推导不平衡力在x、y方向的分量表达式;
(2)确定所述自适应陷波器的传递函数,并通过所述自适应陷波器提取所述转子位移信号中的同频振动信号;
(3)通过自适应陷波器的输出信号估计所述转子的转速,并确定估计转速的表达式;
(4)将经过带频率估计单元的自适应陷波器后得到的不平衡同频量反馈回主动磁悬浮轴承系统中,以消除同频电流,实现转子的自动平衡。
5.根据权利要求4所述的主动磁悬浮轴承系统不平衡控制方法,其特征在于,所述步骤(2)中当估计转速为所述自适应陷波器的中心频率Ω时,输出信号X1(s)将完全跟踪输入的转子的位移信号v,输出信号X2(s)为与输出信号X1(s)幅值绝对值相等的正交信号,此时与转速同频的不平衡振动量被提取出来。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210708009.2A CN115169028A (zh) | 2022-06-22 | 2022-06-22 | 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210708009.2A CN115169028A (zh) | 2022-06-22 | 2022-06-22 | 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115169028A true CN115169028A (zh) | 2022-10-11 |
Family
ID=83487708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210708009.2A Pending CN115169028A (zh) | 2022-06-22 | 2022-06-22 | 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115169028A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115857362A (zh) * | 2023-03-01 | 2023-03-28 | 坎德拉(深圳)新能源科技有限公司 | 储能飞轮转子的同频振动抑制方法、磁轴承控制器 |
-
2022
- 2022-06-22 CN CN202210708009.2A patent/CN115169028A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115857362A (zh) * | 2023-03-01 | 2023-03-28 | 坎德拉(深圳)新能源科技有限公司 | 储能飞轮转子的同频振动抑制方法、磁轴承控制器 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108716471B (zh) | 一种磁悬浮分子泵转子极小位移主动控制方法 | |
CN110145541B (zh) | 一种基于相位稳定的磁悬浮轴承转子不平衡运动控制方法 | |
US5084643A (en) | Virtual rotor balancing in magnetic bearings | |
CN107133387B (zh) | 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法 | |
JPH0572177B2 (zh) | ||
CN102707626B (zh) | 一种基于位移刚度力超前补偿的自动平衡磁悬浮转子系统 | |
CN110762120A (zh) | 一种基于磁悬浮轴承转子系统的高回转精度控制方法 | |
CN113067523B (zh) | 一种基于角域陷波滤波的磁悬浮电机振动抑制方法 | |
CN111752153B (zh) | 一种基于1.5阶混合重复控制器的谐波电流抑制方法 | |
CN111650975B (zh) | 一种基于多阶重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN115169028A (zh) | 一种飞轮储能主动磁悬浮轴承系统不平衡控制方法 | |
CN107727088A (zh) | 一种基于非线性自适应控制的全主动磁轴承系统惯性轴辨识方法 | |
CN115001336A (zh) | 一种全转速磁悬浮高速电机转子振动控制方法 | |
CN114322971B (zh) | 基于双二阶广义积分器的磁悬浮转子同频振动力抑制方法 | |
CN115750592A (zh) | 一种主动磁悬浮轴承转子的解耦和不平衡振动控制方法 | |
CN116696944A (zh) | 一种磁悬浮轴承的不平衡振动抑制方法及系统 | |
CN109687799B (zh) | 基于假定旋转坐标的磁轴承转子分层振动补偿方法 | |
Hutterer et al. | Unbalance compensation of a magnetically levitated rotor for the whole operating range | |
CN113282115B (zh) | 相移广义积分扩张状态观测器、振动抑制方法及系统 | |
KR100940367B1 (ko) | 자기 베어링을 사용하는 회전체에서 회전체 불균형에 의한회전위치오차 보상장치 및 그 방법 | |
CN115126775B (zh) | 一种磁悬浮旋转机械转速估计及不平衡振动抑制方法 | |
CN113623239B (zh) | 一种不平衡磁拉力控制方法、装置、系统、设备和介质 | |
Cui et al. | Harmonic Vibration Suppression of AMB Rotor Based on VSS Repetitive Controller Considering Mass Imbalance and Sensor Runout | |
Marino et al. | Tracking control for sensorless induction motors with uncertain load torque and resistances | |
CN114114919B (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 |