本发明涉及流体机械设计领域，特指涉及一种泵流动诱导振动性能综合评价方法。
背景技术：
泵作用一种通用机械，在众多工业领域中发挥着不可取代的重要作用。随着社会的发展，对泵运行时的稳定性提出了新的要求，较低的振动水平不但可以节省能量提高性能，而且对保证泵的寿命也至关重要。
泵的振动有两种：一是机械振动引起的系统振动，二是流动诱导引起的振动。前者主要是受设计和制造影响，现已通过主动控制技术等得到了较好的解决，而后者主要是由于泵内部的非定常流动所引起的，其作用机理还在不断的研究中。压力脉动作为泵内非定常流动特征的具体表现形式，同时也是引发泵流动诱导振动的主要因素。但是目前对于压力脉动的分析方法大多只关注主频处(最大幅值对应的频率)的幅值分析，对于次主频处(通常是叶频，也可能是轴频、动静干涉频率等)的幅值考虑的较少。虽然这些方法能够在一定程度上反应泵的振动水平，但当次主频等频率处的幅值较大时，这种分析方法就显得不够全面。因此，急需提出一种能够综合评价泵流动诱导振动性能的方法。但是迄今为止，尚未见泵流动诱导振动综合评价方法的相关报道。
技术实现要素：
针对现有技术中存在不足，本发明提供了一种泵流动诱导振动性能综合评价方法，较为全面地反映泵的振动性能。
本发明是通过以下技术手段实现上述技术目的的。
一种泵流动诱导振动性能综合评价方法，包含以下步骤：
步骤一：选取叶轮出口处作为监测点，采用试验测试或数值计算得到监测点的压力脉动数据，并计算出无量纲化的压力脉动系数时域变化函数；
步骤二：对压力脉动系数时域变化函数进行快速傅里叶变化以获得频域变化函数，进行全频域搜索，按幅值从大到小选取前三个频率为计算频率点；
步骤三：采用层次分析法确定计算频率点的幅值在泵流动诱导振动整体评价中的权重因子，通过对各个计算频率点的幅值和与之对应的权重因子的计算获得泵三阶综合振动压强，根据所述泵三阶综合振动压强评判泵流动诱导振动性能，所述泵三阶综合振动压强越大，则泵流动诱导振动性能越差，泵三阶综合振动压强越小，则泵流动诱导振动性能越好。
优选地，所述步骤一中，利用试验测试获取数据时，选取叶轮出口处作为监测点，当泵稳定运行后方可进行数据采样，采样频率fs选定为1000f1，f1为泵的轴频，持续时间t＝2t，t为待测泵的旋转周期；
利用数值计算获取压力脉动数据时，将正确的定常计算结果设置为非定常计算的初始条件，并设置非定常计算时间步长
优选地，所述步骤一中，获取压力脉动数据后，除去压力脉动数据中的无效数据并与时间信息相匹配，获得压力脉动时域变化函数fp(t)，采用excel或origin等数据处理软件将得到的压力脉动时域变化函数fp(t)变化为压力脉动系数cp的时域变化函数fc(t)，实现所选变量的无量纲化，其中
优选地，所述步骤二具体如下：
(1)通过origin或matlab等数据分析软件对得到的压力脉动系数时域变化函数fc(t)进行快速傅里叶变化以获取相应的频域变化函数；
(2)采用matlab或excel或origin等数据分析软件对步骤(1)获取的频域变化函数进行全局搜索，将各频率点的数据按照幅值由大到小的顺序依次排列，选取幅值排在前三位的频率点作为计算频率点。
优选地，所述步骤三具体如下：
s1：根据各计算频率点及其幅值的相互关系构建中间判断层矩阵a，具体如下：
三个所述计算频率点的数据按照幅值由大到小依次排列为(f1,a1)、(f2,a2)、(f3,a3)，定义元素aij(i<j，i和j的取值为1、2或3)的取值为最接近bij的一个正整数，
aij的含义为：关于计算频率点幅值对于泵流动诱导振动的重要性，第i个计算频率点幅值为第j个计算频率点幅值的aij倍，例如：a12＝3，其含义为：第一个计算频率幅值对泵整体振动的重要性为第二个计算频率的3倍。以aij构建中间判断层矩阵
s2：使用规范列平均法对中间判断层a进行计算，获得各个计算频率点幅值的权重因子；
首先对中间判断层矩阵a进行归一化处理，得到矩阵
s3：计算得出泵三阶综合振动压强
优选地，所述数据分析软件为matlab软件或origin软件。
本发明的有益效果：
1)中间过程采用压力脉动系数作为无量纲变量，适用于不同种类泵在不同工况下的振动性能评估，具有广泛的使用前景。
2)利用层次分析法计算出不同频率点幅值的权重因子，能较为全面地表达出各频率点幅值对泵流动诱导振动性能的重要程度。
3)该评价方法涉及了频域上幅值最大的三个频率点，计算出的泵三阶综合振动压强可较全面地、简洁地反映泵流动诱导振动性能。
附图说明
图1为本发明所述一种泵流动诱导振动性能综合评价方法的流程图。
图2为实施例中泵出口处的压力脉动系数时域变化图。
图3为实施例中泵出口处的压力脉动系数频域变化图。
具体实施方式
下面结合附图以及具体实施例对本发明作进一步的说明，但本发明的保护范围并不限于此。
本实施例采用比转数为27的五级离心泵，其设计流量q＝5m3/h，扬程h＝36m，转速n＝2900r/min，叶轮叶片数z＝8，导叶数为6。
如图1所示，一种泵流动诱导振动性能综合评价方法，包含以下步骤：
步骤一：选取叶轮出口处作为监测点，采用试验测试或数值计算得到监测点的压力脉动数据，并计算出无量纲化压力脉动系数时域变化函数，具体如下：
1)选取叶轮出口处作为监测点，采用试验测试或数值计算得到监测点的压力脉动数据，本实施例利用数值计算获取压力脉动数据，将正确的定常计算结果设置为非定常计算的初始条件，并设置非定常计算时间步长
2)除去压力脉动数据中的无效数据并与时间信息相匹配，获得压力脉动时域变化函数fp(t)，采用excel软件软件将得到的压力脉动时域变化函数fp(t)变化为压力脉动系数cp的时域变化函数fc(t)，如图2所示，实现无量纲化，其中
步骤二：对压力脉动系数时域变化函数进行快速傅里叶变化以获得频域变化函数，如图3所示，进行全频域搜索，按幅值从大到小选取前三个频率为计算频率点，具体如下：
(1)通过matlab软件对得到的压力脉动系数时域变化函数fc(t)进行快速傅里叶变化以获取相应的频域变化函数；
(2)采用origin软件对步骤(1)获取的频域变化函数进行全局搜索，将各频率点的数据按照幅值由大到小的顺序依次排列，选取幅值排在前三位的频率点作为计算频率点。
步骤三：采用层次分析法确定计算频率点的幅值在所述综合评价方法中的权重因子，通过对各个计算频率点的幅值和其对应的权重因子进行计算获得泵三阶综合振动压强，根据所述泵三阶综合振动压强评判泵流动诱导振动性能，所述泵三阶综合振动压强越大，泵流动诱导振动性能越差，具体如下：
s1：根据各计算频率点及其幅值的相互关系构建中间判断层矩阵a，具体如下：
三个所述计算频率点的数据按照幅值由大到小依次排列为(f1,a1)、(f2,a2)、(f3,a3)，(f1,a1)、(f2,a2)、(f3,a3)，具体数值为(386.667,0.0246)、(48.333,0.0041)、(870,0.0024)。
定义元素aij(i<j，i和j的取值为1、2或3)的取值为最接近bij的一个正整数，
s2：使用规范列平均法对中间判断层a进行计算，获得各个计算频率点幅值的权重因子；
首先对中间判断层矩阵a进行归一化处理，得到矩阵
s3：计算得出泵三阶综合振动压强
l＝(0.7862*0.0246+0.1396*0.0041+0.0742*0.0024)*93200＝1872.4775，根据此数值可以进行被选泵和其他泵之间的流动诱导振动性能综合评估与比较。l数值小的泵，其流动诱导振动性能优于l数值大的泵。
所述实施例为本发明的优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。