CN109579969B - 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 - Google Patents

叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 Download PDF

Info

Publication number
CN109579969B
CN109579969B CN201811447312.1A CN201811447312A CN109579969B CN 109579969 B CN109579969 B CN 109579969B CN 201811447312 A CN201811447312 A CN 201811447312A CN 109579969 B CN109579969 B CN 109579969B
Authority
CN
China
Prior art keywords
transient
impeller
acceleration
vibration amplitude
deceleration
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
CN201811447312.1A
Other languages
English (en)
Other versions
CN109579969A (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.)
Shanghai Moxiang Machinery Technology Co ltd
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201811447312.1A priority Critical patent/CN109579969B/zh
Publication of CN109579969A publication Critical patent/CN109579969A/zh
Application granted granted Critical
Publication of CN109579969B publication Critical patent/CN109579969B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明提供了一种叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统，包括：稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。本发明通过采用简单的计算，只需要叶轮基本参数及稳态分析即可，从而解决瞬态振动最大幅值需要大量计算量的目的。

Description

叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统
技术领域
本发明涉及计算、测量技术领域，具体地，涉及叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统。
背景技术
目前计算叶轮在加减速瞬态工况下最大振动幅值的方法为杜哈美(Duhamel)积分方法或其他时域下的数值积分方法，该类方法将时间轴划分成一系列的离散点，在每一步应用瞬态力的变化，累计每一瞬态力下的响应，从而获得加减速瞬态工况下的振动情况，最后获得其振动最大幅值。
因此，现有技术中的这种时域下数值积分方法，其虽能获得准确的加减速瞬态工况下的最大振动幅值，但由于计算量巨大，最大振动幅值的获取速度较慢。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统。
根据本发明提供的一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，包括：
稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
较佳的，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。
较佳的，瞬态分析步骤中，根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，
Figure BDA0001885997390000011
根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
较佳的，计算步骤中，将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
根据本发明提供的一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
较佳的，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。
较佳的，瞬态分析模块根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，
Figure BDA0001885997390000021
根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
较佳的，计算步骤模块将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
与现有技术相比，本发明具有如下的有益效果：
本发明通过采用简单的计算，只需要叶轮基本参数及稳态分析即可，从而解决瞬态振动最大幅值需要大量计算量的目的。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明的计算流程图；
图2为特征参数与瞬态放大因子的函数关系图；
图3为本发明实施例的叶轮结构示意图；
图4为本发明实施例叶轮在加减速瞬态工况下的真实振动时域响应图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示，本发明提供的一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，包括：
稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值。
瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子。叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，
Figure BDA0001885997390000031
Figure BDA0001885997390000032
根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子，将获取的稳态最大振动幅值与瞬态放大因子相乘，计算得到叶轮在加减速瞬态工况下的最大振动幅值。
推导过程以单自由度系统为例，亦可推导至多自由度系统。
1)对于单自由度系统的动力学方程
Figure BDA0001885997390000033
其中m、c、k为系统的质量、阻尼、刚度，x为系统的振动响应，t表示时间，F表示激励力的幅值，f为用复数形式表示的激励力，对于稳态激励，f(t)＝eiωrt，对于加减速瞬态激励，
Figure BDA0001885997390000034
其中i为虚数单位，ω为激励频率，a为加速度值，r为激励阶次。
2)其稳态最大振动幅值为
Figure BDA0001885997390000035
其中ξ为材料阻尼比
3)其加减速瞬态工况下的最大振动幅值为
Figure BDA0001885997390000036
其中μ为特征参数，
Figure BDA0001885997390000037
w(z)是中间变量z的faddeeva函数
Figure BDA0001885997390000038
中间变量z是特征参数μ的函数
Figure BDA0001885997390000039
4)定义瞬态放大因子λ和特征参数μ的关系为
Figure BDA00018859973900000310
(如图2所示)
5)从而稳态最大幅值和瞬态放大因子相乘即可得到加减速瞬态工况下的最大振动幅值。
本发明所针对的叶轮包含谐调情况和失谐情况下的叶轮，即谐调情况是指，为理想情况下，各个叶轮的扇区的特性一致；失谐情况是指由于材料不均匀、加工制造误差、实际运行的磨损等因素造成的叶轮各个扇区不一致的状况。
对于如图3所示的叶轮，
1)当激励阶次r＝1时，通过ANSYS Workbench软件对叶轮进行稳态分析，获得稳态振动最大幅值为0.6406mm。
2)在当激励力的激励阶次r＝1，加速度a＝3769.91rad/s2时；本例通过ANSYSWorkbench软件对叶轮进行模态分析并结合第一步的稳态分析，获得叶轮的固有频率和阻尼比如下，叶轮在当前激励阶次下的固有频率为ωn＝3499.52rad/s，叶轮的阻尼比ξ＝0.5％。则通过公式
Figure BDA0001885997390000041
可计算特征参数为0.2850，计算得对应的瞬态放大因子为0.5257。
3)将所得的稳态最大幅值0.6406mm和所得的瞬态放大因子0.5257相乘，即可快速获得叶轮在加减速瞬态工况下的最大振动幅值为0.3368mm。通过ANSYS Workbench软件对叶轮进行为时约2小时的瞬态分析，获得真实的加减速瞬态工况下的最大振动幅值为0.3420mm，如图4所示，相对误差为-1.56％。通过本方法可在保证计算精度的前提下，大幅提高计算效率。
在上述一种叶轮在加减速瞬态工况下最大振动幅值的获取方法的基础上，本发明还提供一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (2)

1.一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，包括：
稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值；
叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn；
瞬态分析步骤中，根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，
Figure FDA0002321951950000011
根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子；
计算步骤中，将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
2.一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值；
叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn；
瞬态分析模块根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，
Figure FDA0002321951950000012
根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子；
计算步骤模块将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
CN201811447312.1A 2018-11-29 2018-11-29 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 Active CN109579969B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811447312.1A CN109579969B (zh) 2018-11-29 2018-11-29 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811447312.1A CN109579969B (zh) 2018-11-29 2018-11-29 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统

Publications (2)

Publication Number Publication Date
CN109579969A CN109579969A (zh) 2019-04-05
CN109579969B true CN109579969B (zh) 2020-05-22

Family

ID=65924059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811447312.1A Active CN109579969B (zh) 2018-11-29 2018-11-29 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统

Country Status (1)

Country Link
CN (1) CN109579969B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4222339C2 (de) * 1992-07-08 1994-04-28 First Order System Verfahren zur vollständigen Zustandsbeobachtung an einer Anlage, die als mechanisches System 1. Ordnung beschreibbar ist
JP2010112769A (ja) * 2008-11-05 2010-05-20 Kyushu Institute Of Technology 時刻歴応答シミュレーション方法、そのプログラム及び装置
CN105745724A (zh) * 2013-12-10 2016-07-06 施耐德电器工业公司 电气装置的短路设备、以及用于消灭可能出现在这种装置中的电弧的系统
CN108829946A (zh) * 2018-05-29 2018-11-16 中国人民解放军战略支援部队航天工程大学 一种基于动态补偿技术的推力计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE502007001846D1 (de) * 2007-03-30 2009-12-10 Ford Global Tech Llc Verfahren zur Detektion periodischer Störungen in der Lenkeinrichtung eines Kraftfahrzeuges sowie Ve

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4222339C2 (de) * 1992-07-08 1994-04-28 First Order System Verfahren zur vollständigen Zustandsbeobachtung an einer Anlage, die als mechanisches System 1. Ordnung beschreibbar ist
JP2010112769A (ja) * 2008-11-05 2010-05-20 Kyushu Institute Of Technology 時刻歴応答シミュレーション方法、そのプログラム及び装置
CN105745724A (zh) * 2013-12-10 2016-07-06 施耐德电器工业公司 电气装置的短路设备、以及用于消灭可能出现在这种装置中的电弧的系统
CN108829946A (zh) * 2018-05-29 2018-11-16 中国人民解放军战略支援部队航天工程大学 一种基于动态补偿技术的推力计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"CAP1400核主泵水力设计及瞬态空化特性研究";陈宗良;《中国优秀硕士学位论文全文数据库》;20161231;全文 *
"Some transfer-function conditions for a desired maximum amplitude or exponential envelope of a closed-loop transient response";Matsuo T 等;《IEEE会议》;19941231;全文 *
"Tntegrated post-processor for 5-axis machine tools with geomaetric errors compensation";毕庆贞 等;《International Journal of Machine Tools & Manufacture》;20151231;第94卷;65-73页 *

Also Published As

Publication number Publication date
CN109579969A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
Lehfuss et al. Comparison of multiple power amplification types for power hardware-in-the-loop applications
CN108646574B (zh) 一种基于幂次吸引律的离散重复控制器及其控制方法
CN104639003B (zh) 一种交流伺服系统的转动惯量辨识方法
CN104573172A (zh) 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置
CN106546400B (zh) 一种多输入多输出非高斯随机振动试验系统及试验方法
CN108983615A (zh) 基于反双曲正弦吸引律的离散双周期重复控制器
CN111523182A (zh) 一种基于流固耦合的叶片振动响应分析方法
CN109117504A (zh) 一种双向功能梯度曲壳振动分析方法
CN115422654A (zh) 基于cfd/csd技术的跨/超声速飞行器高效高精度颤振时域分析方法
Li et al. Iterative learning control with initial rectifying action for nonlinear continuous systems
CN109579969B (zh) 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统
Alhawwary et al. DNS and LES of the flow over the T106C turbine using the high-order FR/CPR method
CN111684706A (zh) 伺服系统阻力特性获取方法、控制方法和系统及存储装置
CN105512352B (zh) 一种面向气动载荷的分析方法
CN111707528A (zh) 一种基于累积损伤等效的动力学环境时频转换方法
CN110702313B (zh) 基于变转速起动高精度识别柔性转子不平衡激振力的方法
CN110657933B (zh) 一种地震模拟振动台迭代控制方法
Ma et al. Fractional order model reference adaptive control for a hydraulic driven flight motion simulator
Zheng et al. A fractional order controller design based on Bode’s ideal transfer function and bode’s ideal cut-off ideas
Malwatkar et al. PID controllers tuning for improved performance of unstable processes
Narwariya et al. Parametric study on Harmonic Analysis of anti-symmetric laminated composite Plate
Sheng et al. Identification for active vibration control of flexible structure based on Prony algorithm
Bretz Progress in the Migration of Flight Test Analysis Routines to Python
Brezina et al. Modeling and control of system with pump and pipeline by pole placement method
Bonhage et al. Transient amplification of maximum vibration amplitudes

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
TR01 Transfer of patent right

Effective date of registration: 20231007

Address after: 200030 No. 1954, Huashan Road, Shanghai, Xuhui District

Patentee after: Bi Qingzhen

Address before: 200240 No. 800, Dongchuan Road, Shanghai, Minhang District

Patentee before: SHANGHAI JIAO TONG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231120

Address after: 201109 Building 1, No. 600, Jianchuan Road, Minhang District, Shanghai

Patentee after: Shanghai Moxiang Machinery Technology Co.,Ltd.

Address before: 200030 No. 1954, Huashan Road, Shanghai, Xuhui District

Patentee before: Bi Qingzhen

TR01 Transfer of patent right