确定模型水轮机转轮叶片头部空化的声学方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103557933103557933A(43)申请公布日2014.02.05(21)申请号201310589186.4(22)申请日2013.11.20(71)申请人哈尔滨电机厂有限责任公司地址150040黑龙江省哈尔滨市香坊区三大动力路99号哈尔滨电机厂有限责任公司技术管理部(72)发明人赵越(51)Int.Cl.G01H17/00(2006.01)G01M15/00(2006.01)权权利要求书3页利要求书3页说明书6页说明书6页附图1页附图1页(54)发明名称确定模型水轮机转轮叶片头部空化的声学方法(57)摘要确定模型水轮机转轮叶片头部空化的声学方法涉及模型水轮机转轮叶片头部空化判别领域。本发明通过分析Δf1／Δf2值随空化系数变化的趋势来判定模型水轮机转轮叶片头部处是否发生空化。在获取模型水轮机运行转轮的声波信号后，建立了水声信号高低频能量分布趋势的数学模型，利用发生模型水轮机转轮叶片头部空化时，表征特征值Δf1／Δf2随空化系数变化趋势的一次分段函数的分界点恰好为模型水轮机转轮叶片头部空化初始发生点这一特性，通过确定该一次分段函数的分界点的空化系数来确定模型水轮机转轮叶片头部空化的发生。CN103557933ACN103579ACN103557933A权利要求书1/3页1.一种确定模型水轮机转轮叶片头部处空化的声学方法，其特征是：该方法包括如下步骤：1）启动模型水轮机试验系统；2）调整模型水轮机运行工况，使模型水轮机转轮叶片头部处于未空化状态；3）保持模型水轮机运行工况稳定，对水声信号进行采集；4）建立模型水轮机水声信号截取后的时间序列：x(n)＝s(n)W(n)式中：x(n)为模型水轮机水声信号截取后的时间序列；s(n)为模型水轮机水声信号采样后的时间序列；W(n)为窗函数；5）计算模型水轮机水声信号的频谱：式中：Sxx为模型水轮机水声信号的功率谱；Δt为模型水轮机水声信号的采样时间间隔；N为模型水轮机水声信号的采样数；X(k)为对x(n)进行傅立叶变换后的频谱函数，k＝0,1,2,3,4,…，N；X*(k)为X(k)的共轭频谱函数，6）利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势，具体方法如下：在水轮机水声信号的功率谱中，能量的变化趋势是连续的，且低频区域呈现出随着水声频率的升高能量值陡降的趋势，而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多，低频区域与高频区域之间存在着一个明显的交界点，将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合，则该一次分段函数可表示为：式中：是一次分段函数；xm是一次分段函数的交点；2CN103557933A权利要求书2/3页k0、k1是一次分段函数中一次项的系数；b0、b1是一次分段函数中常数项；7）利用一次分段函数来表示水声功率谱中低频能量的分布趋势，具体方法如下：利用项6所述方法，在水声功率谱(0，f0)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述，其中，f0为利用项6所述方法计算出的水轮机水声信号功率谱中高低频区域能量的交界点所对应的频率，该一次分段函数的交点所对应的频率记为f1；8）计算Δf1按下式计算Δf1：Δf1＝f0-f19）利用一次分段函数来表示水声功率谱中高频能量的分布趋势，具体方法如下：利用项6所述方法，在水声功率谱(f0,f2)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述，该一次分段函数的交点所对应的频率记为f2；10、计算Δf2按下式计算Δf2：Δf2＝f2-f011、采集不同空化系数下的水声信号：逐步降低模型水轮机的空化系数，在不同的空化系数下重复步骤3至步骤10，直至水轮机转轮叶片头部处出现明显的空化现象；12、确定模型水轮机转轮叶片头部处空化：具体判定方法如下:将Δf1／Δf2随空化系数σ的变化趋势采用最小二乘法拟合的一次分段函数y′来拟合，则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮叶片头部处空化初始发生的位置，该一次分段函数y′可表示为：式中：y′是一次分段函数；x′m是一次分段函数y′的交点；k′0、k′1是一次分段函数y′中一次项的系数；b0′、b1′是一次分段函数y′中常数项；令任一组测量Vi=yi-y′i数据(xi，yi)对一次分段函数y′的偏差为Vi，则有Vi=yi-y′i,即：假设x′i≤x′m时有n1个测量数据，x′i＞x′m时有n2个测量数据，即n1+n2＝N，则偏差的平方和Qi为：3CN103557933A权利要求书3/3页系数的变化趋势偏差的平方和Qi，其中最小Qi值所