专利摘要

专利摘要

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

权利要求

1.一种叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，包括：

稳态分析步骤：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；

瞬态分析步骤：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；

计算步骤：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。

2.根据权利要求1所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ω_n。

3.根据权利要求2所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，瞬态分析步骤中，根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ， 根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。

4.根据权利要求1所述的叶轮在加减速瞬态工况下最大振动幅值的获取方法，其特征在于，计算步骤中，将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。

5.一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，包括：

稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；

瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；

计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。

6.根据权利要求5所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，叶轮的基本参数包括叶轮材料阻尼比ξ和叶轮固有频率ω_n。

7.根据权利要求6所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，瞬态分析模块根据加减速瞬态工况下的激励力的加减速值a和激励阶次r，从而获得特征参数μ， 根据特征参数和瞬态放大因子之间的一一对应的函数关系，计算得到瞬态放大因子。

8.根据权利要求5所述的叶轮在加减速瞬态工况下最大振动幅值的获取系统，其特征在于，计算步骤模块将获取的稳态最大振动幅值与瞬态放大因子相乘，得到叶轮在加减速瞬态工况下的最大振动幅值。

说明书

技术领域

本发明涉及计算、测量技术领域，具体地，涉及叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统。

背景技术

目前计算叶轮在加减速瞬态工况下最大振动幅值的方法为杜哈美(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)＝e^iω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/s²时；本例通过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％。通过本方法可在保证计算精度的前提下，大幅提高计算效率。

在上述一种叶轮在加减速瞬态工况下最大振动幅值的获取方法的基础上，本发明还提供一种叶轮在加减速瞬态工况下最大振动幅值的获取系统，包括：

稳态分析模块：根据叶轮在稳态工况下的振动情况获取稳态最大振动幅值；

瞬态分析模块：对叶轮在加减速瞬态工况下的振动情况以及叶轮的基本参数获取瞬态放大因子；

计算模块：根据获取的稳态最大振动幅值与瞬态放大因子计算得到叶轮在加减速瞬态工况下的最大振动幅值。

本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

叶轮在加减速瞬态工况下最大振动幅值的获取方法及系统专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。