CN113700739A - 用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 - Google Patents
用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 Download PDFInfo
- Publication number
- CN113700739A CN113700739A CN202110947289.8A CN202110947289A CN113700739A CN 113700739 A CN113700739 A CN 113700739A CN 202110947289 A CN202110947289 A CN 202110947289A CN 113700739 A CN113700739 A CN 113700739A
- Authority
- CN
- China
- Prior art keywords
- temperature
- rotating speed
- frequency
- wave trap
- modal
- 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
Classifications
-
- 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
- F16C32/0444—Details of devices to control the actuation of the electromagnets
- F16C32/0451—Details of controllers, i.e. the units determining the power to be supplied, e.g. comparing elements, feedback arrangements with P.I.D. control
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- 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)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Discrete Mathematics (AREA)
- Mechanical Engineering (AREA)
- Electromagnetism (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
Description
技术领域
本发明涉及一种用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组。
背景技术
目前,磁悬浮控制技术在飞轮储能以及航空航天等领域应用广泛,与机械轴承相比磁轴承有更为显著的优点,高速转子系统已经在离心机,高精度数控机床、透平机、储能飞轮、以及磁悬浮飞轮得到广泛应用。但转子的额定转速较高时,磁轴承需要更高的控制带宽,强陀螺效应和非线性因素会使转子出现震荡扰动,影响磁悬浮控制系统的正常运作。磁悬浮飞轮转子的扰动根频域特征可分为同频分量、倍频分量以及模态分量三类,实验测试结果表明,在飞轮转速变化,倍频分量与机械固有频率重合共振的情况下,倍频分量更加突出,此时会产生自激振荡,经过陷波器处理后虽能抑制此状况,但随转速和温度的变化,模态频率也会随着发生变化,此时的陷波器不能有效的将移动变化的模态频率滤除,自激振荡的出现会导致转子在高转速下失稳。同时,经实验发现随着转速升高导致细长型转子变短的微小形变使模态频率变大,额定转速下长时间运行下温度升高使细长型转子变长,使固有频率在原来基础上变小,原来的滤除模态的陷波器就失去作用。
发明内容
为解决此问题,本发明提出了一种用于磁悬浮飞轮的变化模态频率计算方法,应用于磁悬浮系统,计算出不同转速和温度下的模态频率。
本发明所采用的技术方案为:一种用于磁悬浮飞轮的的变化模态频率计算方法,如下,用示波器对AD采样芯片采样的电压信号进行FFT计算,得到静态下飞轮转子的模态频率为f0;当飞轮转子在升速过程中,同时测得变化的模态频率:在常温20°下,转速分别达到5000r/min和15000r/min时,模态频率分别为f1、f2,额定转速15000r/min下长时间运行温度80°时,模态频率为f3,
设转速为n,温度为t时,
在温度20°、转速达到额定转速15000r/min之前的模态频率fa为: f0+f转速;
当转速n达到额定转速15000r/min时,温度为t时的模态频率fb为: f0+f转速+f温度。
本发明的另一目的是提供采用如上所述的变化模态频率计算方法得出的一种自适应陷波器组,将飞轮转子在不同温度和转速条件下的变化模态频率滤除,抑制高速飞轮转子出现扰动震荡,使其保持稳定。
本发明所采用的技术方案为:采用如上方法的一种自适应陷波器组包括多个陷波器,分别滤除:在温度20°时,飞轮转子,静态时的模态频率,同时在温度20°时,随转速提高,转速8000r/min 时的模态频率、转速12000r/min时的模态频率、到达额定转速 15000r/min时的模态频率;保持额定转速15000r/min时,在温度 50°时的模态频率和在温度70°时的模态频率。
进一步的,本发明还具有如下技术特征:所述的自适应陷波器组包括6个陷波器0-5,陷波器0采用切比雪夫Ⅱ陷波器,在温度20°时,陷波器0滤除飞轮转子静态时的模态频率,同时在温度 20°时,随转速提高,陷波器1滤除转速8000r/min时的模态频率、陷波器2滤除转速12000r/min时的模态频率,陷波器3滤除到达额定转速15000r/min时的模态频率;保持额定转速15000r/min时,陷波器4滤除在温度50°时的模态频率,陷波器5滤除在温度70°时的模态频率。
进一步的,提供采用如上所述的自适应陷波器组得出的一种分段滤除变化模态频率的方法,如下:飞轮转子静态悬浮启用陷波器0,判断转速到达8000r/min时,启用陷波器1,没有达到时用陷波器0,陷波器1启用后判断转速是否达到12000r/min,达到时启用陷波器2,没有达到时返回上一个判断条件,陷波器2启用后判断转速是否达到15000r/min,达到时启用陷波器3,没有达到时返回上一个判断条件,此时已达到额定转速15000r/min,再判断温度条件:温度在大于40°小于60°时用陷波器4,没有达到这一温度条件返回上一判断条件,温度大于60°小于80°时用陷波器5,没有达到时返回上一个判断条件。
本发明的优点及有益效果:本发明提高了磁悬浮飞轮系统的运行的安全性,有效消除了系统的自激振荡。
附图说明
图1为本发明控制系统结构简图。
图2为本发明程序流程图。
具体实施方式
下面根据说明书附图举例对本发明做进一步的说明:
实施例1
本实施例公开的一种磁悬浮控制系统,由DSP芯片和FPGA芯片共同作用,DSP控制程序中包含PID算法、陷波器滤波算法和交叉算法以及陀螺效应补偿,DSP芯片将调理后的输出信号发送给FPGA 芯片,FPGA芯片产生PWM信号通过驱动功率桥得到期望电流。
如图1所示,DSP1将经过处理的输出信号发送给FPGA2, FPGA2输出pwm信号在经过功率桥3得到期望电流,作用到飞轮转子4上,传感器5从飞轮采集的转速和温度信号在DSP1程序中判断启用不同的陷波器,根据AD采样芯片6的电流信号和位置电压信号反馈到DSP1形成闭环实现飞轮悬浮。
本实施例以某型号飞轮为实验样本,转子的结构型式为G50号钢,分析计算得到一阶模态为635HZ,设计用切比雪夫Ⅱ陷波器将静态一阶模态滤除,在控制系统中解决模态变化的问题可用分段陷波滤波器将移动的模态滤除,实现磁悬浮功能正常运作,为验证提出设计方法的有效性,通过实验验证。根据实验数据采样到几组数据:
在常温20°下,转速分别达到5000r/min和15000r/min时,模态频率分别为690HZ、880HZ,额定转速下15000r/min长时间运行直到温度达到40°和80°时,模态频率又变为760HZ和580HZ,那么根据样机的实验数据分析,设转速为n,温度为t,
在常温20°下转速达到额定转速之前的模态频率fa为:
转速达到额定转速,长时间运行之后的模态频率fb为:
根据上述的变化模态频率的计算公式,设计出6个陷波器。陷波器0是滤除静态下的模态频率635HZ,陷波器1是滤除 8000r/min时的模态频率690HZ,陷波器2是滤除12000r/min下的模态频率770HZ,陷波器3是滤除15000r/min下的模态频率880HZ,陷波器4是滤除额定转速15000r/min时,温度到50°的模态频率 730HZ,陷波器5是滤除额定转速下温度到70°的模态频率630HZ。
采用如上的自适应陷波器组的得出的一种分段滤除变化模态频率的方法,如下,飞轮转子静态悬浮启用陷波器0,判断转速到达8000r/min时,启用陷波器1,没有达到时用陷波器0,陷波器1 启用后判断转速是否达到12000r/min,达到时启用陷波器2,没有达到时返回上一个判断条件,陷波器2启用后判断转速是否达到 15000r/min,达到时启用陷波器3,没有达到时返回上一个判断条件,此时已达到额定转速,以此类推再以同样的方式判断温度条件:温度在大于40°小于60°时用陷波器4,没有达到这一温度条件返回上一判断条件,温度大于60°小于80°时用陷波器5,没有达到时返回上一个判断条件。
Claims (4)
1.一种用于磁悬浮飞轮的的变化模态频率计算方法,其特征在于,方法如下,用示波器对AD采样芯片采样的电压信号进行FFT计算,得到静态下飞轮转子的模态频率为f0;当飞轮转子在升速过程中,同时测得变化的模态频率:在常温20°下,转速分别达到5000r/min和15000r/min时,模态频率分别为f1、f2,额定转速15000r/min下长时间运行温度80°时,模态频率为f3,
设转速为n,温度为t时,
不同转速和温度下的模态频率计算如下:
在温度20°、转速达到额定转速15000r/min之前的模态频率fa为:f0+f转速;
当转速n达到额定转速15000r/min时,温度为t时的模态频率fb为:f0+f转速+f温度。
2.根据权利要求1所述的一种用于磁悬浮飞轮的的变化模态频率计算方法得出的一种自适应陷波器组,其特征在于,包括多个陷波器,分别滤除:在温度20°时,飞轮转子,静态时的模态频率,同时在温度20°时,随转速提高,转速8000r/min时的模态频率、转速12000r/min时的模态频率、到达额定转速15000r/min时的模态频率;保持额定转速15000r/min时,在温度50°时的模态频率和在温度70°时的模态频率。
3.根据权利要求2所述的一种自适应陷波器组,其特征在于:包括6个陷波器0-5,在温度20°时,陷波器0采用切比雪夫Ⅱ陷波器滤除静态时的模态频率,同时在温度20°时,随转速提高,陷波器1滤除转速8000r/min时的模态频率、陷波器2滤除转速12000r/min时的模态频率,陷波器3滤除到达额定转速15000r/min时的模态频率;保持额定转速15000r/min时,陷波器4滤除在温度50°时的模态频率,陷波器5滤除在温度70°时的模态频率。
4.采用根据权利要求3所述的一种自适应陷波器组得出的一种分段滤除变化模态频率的方法,其特征在于,方法如下:飞轮转子静态悬浮启用陷波器0,判断转速到达8000r/min时,启用陷波器1,没有达到时用陷波器0,陷波器1启用后判断转速是否达到12000r/min,达到时启用陷波器2,没有达到时返回上一个判断条件,陷波器2启用后判断转速是否达到15000r/min,达到时启用陷波器3,没有达到时返回上一个判断条件,此时已达到额定转速15000r/min,再判断温度条件:温度在大于40°小于60°时用陷波器4,没有达到这一温度条件返回上一判断条件,温度大于60°小于80°时用陷波器5,没有达到时返回上一个判断条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110947289.8A CN113700739B (zh) | 2021-08-18 | 2021-08-18 | 用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110947289.8A CN113700739B (zh) | 2021-08-18 | 2021-08-18 | 用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113700739A true CN113700739A (zh) | 2021-11-26 |
CN113700739B CN113700739B (zh) | 2023-02-03 |
Family
ID=78653299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110947289.8A Active CN113700739B (zh) | 2021-08-18 | 2021-08-18 | 用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113700739B (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1086628A (en) * | 1965-04-02 | 1967-10-11 | Westwind Turbines Ltd | Improvements in or relating to the measurement of mass flow rate of gases |
JPH05231428A (ja) * | 1991-06-13 | 1993-09-07 | Hitachi Ltd | 磁気軸受の制御方法および制御装置 |
JPH0612827U (ja) * | 1992-07-21 | 1994-02-18 | セイコー精機株式会社 | 磁気軸受装置 |
US5833173A (en) * | 1996-02-13 | 1998-11-10 | The Boeing Company | Aircraft frequency adaptive modal suppression system |
DE102009038775A1 (de) * | 2009-08-25 | 2011-03-03 | Siemens Aktiengesellschaft | Regeleinrichtung und Verfahren zur Regelung eines Magnetlagers |
CN102011799A (zh) * | 2010-12-01 | 2011-04-13 | 北京奇峰聚能科技有限公司 | 一种高可靠储能飞轮磁轴承数字控制系统 |
CN102169046A (zh) * | 2010-12-22 | 2011-08-31 | 北京航空航天大学 | 一种磁悬浮机电设备弹性模态在线测试系统 |
CN102606504A (zh) * | 2012-03-29 | 2012-07-25 | 北京中科科仪股份有限公司 | 一种温度和转速检测机构 |
US20170114757A1 (en) * | 2015-10-21 | 2017-04-27 | Andrew Currie | Rotating, self-excited, asymmetric radio frequency resonant cavity turbine for energy storage and power production |
US20170175564A1 (en) * | 2015-12-16 | 2017-06-22 | Daniel Schlak | Flywheel with Inner Turbine, Intermediate Compressor, and Outer Array of Magnets |
CN107271126A (zh) * | 2017-05-31 | 2017-10-20 | 北京航空航天大学 | 一种基于磁悬浮技术的转子模态频率及模态振型在线测试装置 |
CN107797450A (zh) * | 2017-10-11 | 2018-03-13 | 北京航空航天大学 | 一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法 |
CN108984815A (zh) * | 2018-04-20 | 2018-12-11 | 江苏大学 | 基于工况的车载飞轮电池径向悬浮支承系统的控制方法 |
CN109038936A (zh) * | 2018-09-29 | 2018-12-18 | 李平 | 储能姿控两用同心反转双飞轮机电装置 |
CN110345158A (zh) * | 2019-07-22 | 2019-10-18 | 珠海格力电器股份有限公司 | 轴承控制方法、装置、存储介质及处理器 |
US10629234B1 (en) * | 2019-06-27 | 2020-04-21 | Western Digital Technologies, Inc. | Data storage device computing notch filters for multiple temperatures |
-
2021
- 2021-08-18 CN CN202110947289.8A patent/CN113700739B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1086628A (en) * | 1965-04-02 | 1967-10-11 | Westwind Turbines Ltd | Improvements in or relating to the measurement of mass flow rate of gases |
JPH05231428A (ja) * | 1991-06-13 | 1993-09-07 | Hitachi Ltd | 磁気軸受の制御方法および制御装置 |
JPH0612827U (ja) * | 1992-07-21 | 1994-02-18 | セイコー精機株式会社 | 磁気軸受装置 |
US5833173A (en) * | 1996-02-13 | 1998-11-10 | The Boeing Company | Aircraft frequency adaptive modal suppression system |
DE102009038775A1 (de) * | 2009-08-25 | 2011-03-03 | Siemens Aktiengesellschaft | Regeleinrichtung und Verfahren zur Regelung eines Magnetlagers |
CN102011799A (zh) * | 2010-12-01 | 2011-04-13 | 北京奇峰聚能科技有限公司 | 一种高可靠储能飞轮磁轴承数字控制系统 |
CN102169046A (zh) * | 2010-12-22 | 2011-08-31 | 北京航空航天大学 | 一种磁悬浮机电设备弹性模态在线测试系统 |
CN102606504A (zh) * | 2012-03-29 | 2012-07-25 | 北京中科科仪股份有限公司 | 一种温度和转速检测机构 |
US20170114757A1 (en) * | 2015-10-21 | 2017-04-27 | Andrew Currie | Rotating, self-excited, asymmetric radio frequency resonant cavity turbine for energy storage and power production |
US20170175564A1 (en) * | 2015-12-16 | 2017-06-22 | Daniel Schlak | Flywheel with Inner Turbine, Intermediate Compressor, and Outer Array of Magnets |
CN107271126A (zh) * | 2017-05-31 | 2017-10-20 | 北京航空航天大学 | 一种基于磁悬浮技术的转子模态频率及模态振型在线测试装置 |
CN107797450A (zh) * | 2017-10-11 | 2018-03-13 | 北京航空航天大学 | 一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法 |
CN108984815A (zh) * | 2018-04-20 | 2018-12-11 | 江苏大学 | 基于工况的车载飞轮电池径向悬浮支承系统的控制方法 |
CN109038936A (zh) * | 2018-09-29 | 2018-12-18 | 李平 | 储能姿控两用同心反转双飞轮机电装置 |
US10629234B1 (en) * | 2019-06-27 | 2020-04-21 | Western Digital Technologies, Inc. | Data storage device computing notch filters for multiple temperatures |
CN110345158A (zh) * | 2019-07-22 | 2019-10-18 | 珠海格力电器股份有限公司 | 轴承控制方法、装置、存储介质及处理器 |
Non-Patent Citations (6)
Title |
---|
崔培玲等: "主被动磁悬浮转子的不平衡振动自适应控制", 《光学精密工程》 * |
崔培玲等: "基于相移陷波器的磁轴承不平衡振动全频自适应控制", 《振动与冲击》 * |
徐敬勃等: "某磁浮飞轮转子系统有限元分析", 《计算机辅助工程》 * |
苏振中等: "磁轴承用全桥功率放大器死区效应分析与补偿算法设计", 《中国电机工程学报》 * |
郑世强等: "磁悬浮分子泵高速转子章动相位裕度跟踪补偿控制", 《机械工程学报》 * |
韩东江等: "气膜约束对轴系固有频率影响的试验", 《航空动力学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113700739B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Active balancing control of AMB-rotor systems using a phase-shift notch filter connected in parallel mode | |
CN107070342A (zh) | 一种带负载状态观测器的永磁同步电机控制系统 | |
CN104647132B (zh) | 一种基于磁悬浮轴承电主轴的铣削颤振主动控制方法 | |
CN111682815B (zh) | 一种含高频干扰重构的风力机舱偏航控制方法 | |
CN104214045A (zh) | 双馈式变速变桨风力发电机组的独立变桨距控制方法 | |
CN104300863A (zh) | 一种变负载永磁同步电机调速的自适应滑模控制方法 | |
CN106681152A (zh) | 一种利用pso优化陷波器参数进行伺服共振抑制的方法 | |
CN113765453B (zh) | 宽-窄极特征的磁悬浮开关磁阻电机悬浮控制系统 | |
CN111478365B (zh) | 一种直驱风电机组虚拟同步机控制参数的优化方法及系统 | |
CN114263564B (zh) | 考虑不变桨风速范围的风电机组有功功率控制方法及系统 | |
CN113700739B (zh) | 用于磁悬浮飞轮的变化模态频率计算方法及自适应陷波器组 | |
Tang et al. | Active vibration control of the flexible rotor to pass the first bending critical speed in high energy density magnetically suspended motor | |
CN110212574B (zh) | 考虑虚拟惯量的风电控制参数协调设置方法 | |
CN114884143A (zh) | 基于转子动能调节的风电机组输出功率虚拟滤波控制方法 | |
Kaneko | Study on transient vibration of mistuned bladed disk passing through resonance | |
Chen et al. | Indirect model reference adaptive control for a class of linear fractional order systems | |
Shi et al. | Model predictive control of direct-driven surface-mounted permanent magnet synchronous generator based on active disturbance rejection controller | |
Wei et al. | A simplified analysis method and suppression of the modalities of a magnetic levitation turbo rotor system | |
CN112486019B (zh) | 不确定风力发电机系统的最大功率跟踪模糊控制方法 | |
Mohamed et al. | New hybrid sensorless speed of a non-salient pole PMSG coupled to wind turbine using a modified switching algorithm | |
CN115288929A (zh) | 一种计及载荷计算的风电机组偏航控制方法及系统 | |
Sun et al. | Unbalance vibration compensation of magnetic bearing systems based on beetle antennae search algorithm | |
CN113326647A (zh) | 一种水轮机轴系转子动力学模态计算方法 | |
CN114243762A (zh) | 一种风机并网的分析与控制方法 | |
Meingast et al. | Investigation of rotor-casing interactions in the centrifugal compressor of a helicopter engine |
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 |