CN109579969A - 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 - Google Patents
叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 Download PDFInfo
- Publication number
- CN109579969A CN109579969A CN201811447312.1A CN201811447312A CN109579969A CN 109579969 A CN109579969 A CN 109579969A CN 201811447312 A CN201811447312 A CN 201811447312A CN 109579969 A CN109579969 A CN 109579969A
- Authority
- CN
- China
- Prior art keywords
- impeller
- transient
- vibration amplitude
- deceleration
- maximum vibration
- 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
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，从而获得特征参数μ，根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
较佳的，计算步骤中，将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
根据本发明提供的一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
较佳的，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。
较佳的，瞬态分析模块根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
较佳的，计算步骤模块将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
与现有技术相比，本发明具有如下的有益效果：
本发明通过采用简单的计算，只需要叶轮基本参数及稳态分析即可，从而解决瞬态振动最大幅值需要大量计算量的目的。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明的计算流程图；
图2为特征参数与瞬态放大因子的函数关系图；
图3为本发明实施例的叶轮结构示意图；
图4为本发明实施例叶轮在加减速瞬态工况下的真实振动时域响应图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示，本发明提供的一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，包括：
稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值。
瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子。叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ， 根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子，将获取的稳态最大振动幅值与瞬态放大因子相乘，计算得到叶轮在加减速瞬态工况下的最大振动幅值。
推导过程以单自由度系统为例，亦可推导至多自由度系统。
1)对于单自由度系统的动力学方程
其中m、c、k为系统的质量、阻尼、刚度，x为系统的振动响应，t表示时间，F表示激励力的幅值，f为用复数形式表示的激励力，对于稳态激励，f(t)＝eiωrt，对于加减速瞬态激励，其中i为虚数单位，ω为激励频率，a为加速度值，r为激励阶次。
2)其稳态最大振动幅值为其中ξ为材料阻尼比
3)其加减速瞬态工况下的最大振动幅值为其中μ为特征参数，w(z)是中间变量z的faddeeva函数中间变量z是特征参数μ的函数
4)定义瞬态放大因子λ和特征参数μ的关系为(如图2所示)
5)从而稳态最大幅值和瞬态放大因子相乘即可得到加减速瞬态工况下的最大振动幅值。
本发明所针对的叶轮包含谐调情况和失谐情况下的叶轮，即谐调情况是指，为理想情况下，各个叶轮的扇区的特性一致；失谐情况是指由于材料不均匀、加工制造误差、实际运行的磨损等因素造成的叶轮各个扇区不一致的状况。
对于如图3所示的叶轮，
1)当激励阶次r＝1时，通过ANSYS Workbench软件对叶轮进行稳态分析，获得稳态振动最大幅值为0.6406mm。
2)在当激励力的激励阶次r＝1，加速度a＝3769.91rad/s2时；本例通过ANSYSWorkbench软件对叶轮进行模态分析并结合第一步的稳态分析，获得叶轮的固有频率和阻尼比如下，叶轮在当前激励阶次下的固有频率为ωn＝3499.52rad/s，叶轮的阻尼比ξ＝0.5％。则通过公式可计算特征参数为0.2850，计算得对应的瞬态放大因子为0.5257。
3)将所得的稳态最大幅值0.6406mm和所得的瞬态放大因子0.5257相乘，即可快速获得叶轮在加减速瞬态工况下的最大振动幅值为0.3368mm。通过ANSYS Workbench软件对叶轮进行为时约2小时的瞬态分析，获得真实的加减速瞬态工况下的最大振动幅值为0.3420mm，如图4所示，相对误差为-1.56％。通过本方法可在保证计算精度的前提下，大幅提高计算效率。
在上述一种叶轮在加减速瞬态工况下最大振动幅值的获取方法的基础上，本发明还提供一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。
Claims (8)
1.一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，包括：
稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
2.根据权利要求1所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。
3.根据权利要求2所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，瞬态分析步骤中，根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
4.根据权利要求1所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，计算步骤中，将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
5.一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，包括：
稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；
瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；
计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。
6.根据权利要求5所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ωn。
7.根据权利要求6所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，瞬态分析模块根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ，根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。
8.根据权利要求5所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，计算步骤模块将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。
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 true CN109579969A (zh) | 2019-04-05 |
| CN109579969B 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 (5)
| 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 |
| US20080243329A1 (en) * | 2007-03-30 | 2008-10-02 | Jorg Hamel | Detection and Compensation of Periodic Disturbances in a Motor Vehicle Steering Device |
| 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 | 中国人民解放军战略支援部队航天工程大学 | 一种基于动态补偿技术的推力计算方法 |
-
2018
- 2018-11-29 CN CN201811447312.1A patent/CN109579969B/zh active Active
Patent Citations (5)
| 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 |
| US20080243329A1 (en) * | 2007-03-30 | 2008-10-02 | Jorg Hamel | Detection and Compensation of Periodic Disturbances in a Motor Vehicle Steering Device |
| 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 (4)
| Title |
|---|
| MATSUO T 等: ""Some transfer-function conditions for a desired maximum amplitude or exponential envelope of a closed-loop transient response"", 《IEEE会议》 * |
| 尚海 等: ""基于主成分分析法的失谐叶片加工误差源识别"", 《机电一体化》 * |
| 毕庆贞 等: ""Tntegrated post-processor for 5-axis machine tools with geomaetric errors compensation"", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS & MANUFACTURE》 * |
| 陈宗良: ""CAP1400核主泵水力设计及瞬态空化特性研究"", 《中国优秀硕士学位论文全文数据库》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109579969B (zh) | 2020-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| González et al. | Steady and unsteady radial forces for a centrifugal pump with impeller to tongue gap variation | |
| Ernst et al. | ARCH-COMP 2019 Category Report: Falsification. | |
| CN109000776B (zh) | 一种涡轮机叶片振动参数分析方法及分析装置 | |
| CN106547957B (zh) | 一种旋转环状周期结构参激弹性振动分析方法 | |
| Ashcroft et al. | On the development of a harmonic balance method for aeroelastic analysis | |
| Rahmati et al. | Nonlinear time and frequency domain methods for multirow aeromechanical analysis | |
| CN105758604B (zh) | 基于加速度等效的多轴振动试验条件裁剪方法 | |
| CN109540482B (zh) | 一种涡轮机叶片无键相同步振动参数分析方法及分析装置 | |
| Chen et al. | Extraction the unbalance features of spindle system using wavelet transform and power spectral density | |
| Sicot et al. | Time-domain harmonic balance method for aerodynamic and aeroelastic simulations of turbomachinery flows | |
| Zhou et al. | The analytical imbalance response of Jeffcott rotor during acceleration | |
| CN110504695B (zh) | 一种评估双馈风机在强迫振荡中参与度的方法 | |
| CN106840685B (zh) | 一种燃气轮机动态工况数据库样本点选取方法 | |
| CN103955143B (zh) | 一种交流永磁电机鲁棒控制器的参数整定方法 | |
| Zhai et al. | Influence of cutting parameters on force coefficients and stability in plunge milling | |
| CN109579969A (zh) | 叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统 | |
| CN112327957B (zh) | 转子低频振动多阶线谱的控制方法及装置 | |
| Novo et al. | Non linear simulation of a wave energy converter with multiple degrees of freedom using a harmonic balance method | |
| CN110829465A (zh) | 计及多调速器死区的电力系统超低频振荡风险评估方法 | |
| Reutter et al. | Comparison of Experiments, Full-annulus Calculations and Harmonic-balance-calculations of a Multi-stage Compressor | |
| CN104049541B (zh) | 一种直流电机鲁棒控制器的参数整定方法 | |
| CN105048920A (zh) | 一种考虑励磁系统调节作用的改进同步发电机抛载试验参数辨识方法 | |
| CN108491568B (zh) | 一种取向硅钢叠片损耗计算方法及系统 | |
| CN106126840A (zh) | 基于轴心轨迹的旋转机械转速计算装置及方法 | |
| Gao et al. | Resonance speed measurement of high-speed spindle using an instruction-domain-based approach |
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 |