CN109740284A - 一种应用于动态翼型转捩判断的可变滑移窗口技术 - Google Patents
一种应用于动态翼型转捩判断的可变滑移窗口技术 Download PDFInfo
- Publication number
- CN109740284A CN109740284A CN201910054491.0A CN201910054491A CN109740284A CN 109740284 A CN109740284 A CN 109740284A CN 201910054491 A CN201910054491 A CN 201910054491A CN 109740284 A CN109740284 A CN 109740284A
- Authority
- CN
- China
- Prior art keywords
- window
- rms
- technology
- aerofoil profile
- value
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种应用于动态翼型转捩判断的可变滑移窗口技术,所述技术中脉动压力均方根值(RMS)使用式1进行计算本技术中构造一种可变的矩形窗形式,如式2所示,W=h (2)。本技术中在窗函数增大的同时,窗口宽度变窄;在窗函数减小的同时,窗口宽度变宽,本技术中窗函数与窗口宽度的变化规律由式3进行控制:m·h=m0 (3) 窗函数h由式4进行描述翼型俯仰震荡规律为:α=α0+A sin(2πft) (5)。本发明有效地解决了特征位置处迎角不确定度较大的问题。
Description
技术领域
本发明涉及数据处理方法技术领域,具体为一种应用于动态翼型转捩判断的可变滑移窗口技术。
背景技术
动态翼型的非定常气动特征在风力机叶片,直升机旋翼叶片,喷气发动机压缩机中起重要作用,所以动态翼型的非定常气动特征需要深入研究。但是,目前的翼型设计理论依然依赖于静态计算和设计准则,急需将翼型的非定常气动特性考虑在内。在翼型的动态运动过程中,人们往往希望精准捕捉翼型的非定常转捩和再次层流化特征。我们使用脉动压力方法对动态翼型的非定常特征进行了研究,并发展了“相位平均方法”和“滑移窗口技术”。“相位平均方法”不需要很高的采样率即可获得翼型的非定常特征,但需要很长的采样时间。“滑移窗口技术”解决了相位平均方法采样时间过长的不足,该技术理论上只需采集一个震荡周期的脉动压力数据,即可捕捉到翼型的非定常特征。
现有技术有如下不足:
1、“滑移窗口技术”但仍有不足,在一个完整的运动周期内,由于流动分离、再附的特征极强,导致转捩、再附特征不够明显;
2、特征位置处迎角不确定度较大。
发明内容
本发明的目的在于提供一种应用于动态翼型转捩判断的可变滑移窗口技术,以解决上述背景技术中提出的“滑移窗口技术”但仍有不足,在一个完整的运动周期内,由于流动分离、再附的特征极强,导致转捩、再附特征不够明显和特征位置处迎角不确定度较大的问题。
一种应用于动态翼型转捩判断的可变滑移窗口技术,所述技术中脉动压力均方根值(RMS)使用式1进行计算
本技术中构造一种可变的矩形窗形式,如式2所示,
W=h (2)
本技术中在窗函数增大的同时,窗口宽度变窄;在窗函数减小的同时,窗口宽度变宽,本技术中窗函数与窗口宽度的变化规律由式3进行控制:
m·h=m0 (3)
窗函数h由式4进行描述
翼型俯仰震荡规律为:
α=α0+Asin(2πft) (5)
由采样率fs引起的一个采样间隔内的角度误差Δαs为:
窗口内的角度误差即为:
Δα=2mΔαs=4mπAf/fscos(2πft) (7)
对于窗函数h,有以下性质:当RMS/RMS0值小于阈值μ0时,h随RMS/RMS0增大而缓慢增大,接近阈值μ0时,h迅速增大;当RMS/RMS0值达到阈值时,h 取最大值;RMS/RMS0值大于阈值时,h迅速减小,由于半窗口宽度不能无限大,式3和式4决定了h有最小值。
根据h函数以上的性质特点,本技术中提出了一种窗口高度变化函数:类高斯函数。
式中,c1为窗函数的放大系数,c2为窗函数的形状系数,μ0为的阈值,c3为窗函数的平移量。
进一步的,所述技术中W(t)是窗函数,是采样时刻,是采样间隔,是窗口宽度,j是窗口内数据的序号,i表示第i个窗口,p是采集到的压力,表示第i个窗口内压力平均值,代表第i个窗口内压力的均方根值(RMS)。
进一步的,所述技术中的W(t)选择为矩形窗。
进一步的,所述技术中m是半窗口宽度,h是窗函数的值,m0是基准半窗口宽度。
进一步的,所述技术中RMS是窗口内的RMS值,RMS0是参考值。
进一步的,所述技术中α0是平均迎角,f是振荡频率,A是翼型震荡的振幅。
进一步的,所述技术中Ω是翼型俯仰运动的角速度,fs是采样率,Δt是采样率引起的采样时间间隔。
与现有技术相比,本发明的有益效果是:
1、本发明中当选定阈值后,h的最大值即为c1+c3,如果,在窗函数W 的作用下,在阈值处,窗口内的RMS值会较“滑移窗口技术”的计算值大,起到了增强特征的作用,很好的解决了在一个完整的运动周期内,由于流动分离、再附的特征极强,导致转捩、再附特征不够明显的问题;
2、本发明中由于特征位置(RMS/RMS0=)处h的增大,在式3的作用下,会引起m的减小,在式7的作用下,m的减小会使角度误差减小,有效地解决了特征位置处迎角不确定度较大的问题。
附图说明
图1为本发明一种应用于动态翼型转捩判断的可变滑移窗口技术窗口高度h函数形状;
图2为本发明一种应用于动态翼型转捩判断的可变滑移窗口技术窗口高度和宽度变化曲线;
图3为本发明一种应用于动态翼型转捩判断的可变滑移窗口技术“滑移窗口技术”与“可变滑移窗口技术”处理结果对比;
图4为本发明一种应用于动态翼型转捩判断的可变滑移窗口技术“滑移窗口技术”与“可变滑移窗口技术”处理结果的迎角误差。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
请参阅图1-4,本发明提供如下技术方案:一种应用于动态翼型转捩判断的可变滑移窗口技术,首先确定类Gauss曲线的参数(选择参数m0,RMS0,c1, c2,c3和μ0,基准半窗口宽度m0:该值是可变滑移窗口技术中的不变量,m0取为m0=30~50均可,窗口高度变化函数平移量c3:c3是h函数的最小值,同时也意味着最大的半窗口宽度,本发明令最大的半窗口宽度为mmax=m0,这样使得hmin=m0/mmax=1,所以c3取值为:c3=hmin=1,窗口高度变化函数放大系数 c1:h函数最大值为c1+c3,考虑到窗口高度不能过分大,本发明取c1=3(如果窗口高度过大,式3会引起窗口宽度过小,实际应用中,窗口宽度不能小于1),窗口高度变化函数形状系数c2:该系数影响了窗口高度的变化速率,一般情况下,c2取为c2=1,RMS的参考值RMS0:该值可选择为传感器在风洞中的零状态情况下的RMS值,阈值μ0,视不同的试验条件及目的进行选择,本技术中可变滑移窗口技术使得转捩/再次层流化处的特征凸现出来,本技术中选择转捩/再次层流化处的RMS值与RMS0之比作为μ0,
使用不可变矩形窗计算基准窗口m0内的脉动信号RMS值,技术中脉动压力均方根值(RMS)使用式1进行计算
本技术中构造一种可变的矩形窗形式,如式2所示,
W=h (2)
本技术中在窗函数增大的同时,窗口宽度变窄;在窗函数减小的同时,窗口宽度变宽,本技术中窗函数与窗口宽度的变化规律由式3进行控制:
m·h=m0 (3)
计算使用式8计算新的窗口高度h窗函数h由式4进行描述
翼型俯仰震荡规律为:
α=α0+Asin(2πft) (5)
由采样率fs引起的一个采样间隔内的角度误差Δαs为:
窗口内的角度误差即为:
Δα=2mΔαs=4mπAf/fscos(2πft) (7)
对于窗函数h,有以下性质:当RMS/RMS0值小于阈值μ0时,h随RMS/RMS0增大而缓慢增大,接近阈值μ0时,h迅速增大;当RMS/RMS0值达到阈值时,h 取最大值;RMS/RMS0值大于阈值时,h迅速减小,由于半窗口宽度不能无限大,式3和式4决定了h有最小值。
根据h函数以上的性质特点,本技术中提出了一种窗口高度变化函数:类高斯函数。
式中,c1为窗函数的放大系数,c2为窗函数的形状系数,μ0为的阈值,c3为窗函数的平移量,使用新的窗函数W计算新的加权RMS值
本技术中W(t)是窗函数,是采样时刻,是采样间隔,是窗口宽度,j是窗口内数据的序号,i表示第i个窗口,p是采集到的压力,表示第i个窗口内压力平均值,代表第i个窗口内压力的均方根值(RMS);本技术中的W(t) 选择为矩形窗;本技术中m是半窗口宽度,h是窗函数的值,m0是基准半窗口宽度;本技术中RMS是窗口内的RMS值,RMS0是参考值;技术中α0是平均迎角,f是振荡频率,A是翼型震荡的振幅;技术中Ω是翼型俯仰运动的角速度,fs是采样率,Δt是采样率引起的采样时间间隔。
本技术中使用的是S809翼型的实验数据,试验雷诺数为Re=0.75×106,翼型做俯仰震荡的平均迎角为8°,振幅为10°,振荡频率为0.5Hz,采样率为 10kHz,采样时间为20s。
选择翼型沿弦向x=50mm处采集的脉动压力作为算例,计算过程中参数的选择如
Table 1所示:
Table 1 参数选择
窗口高度h函数的形状如图1所示。
图2为一个完整的周期内窗口高度和宽度的变化曲线,在转捩/再次层流化位置处,窗口高度增加,这将使该位置处的RMS值增大;同时,转捩/再次层流化处的窗口宽度变小,这将使该处迎角误差减小。
图3为滑移窗口技术与可变滑移窗口技术的处理结果对比,相比于滑移窗口技术,可变滑移窗口技术在转捩/再次层流化处的RMS值有明显增大。
本技术中使用无量纲脉冲强度(INB)来评价信号特征,如式9所示:
式中,INBi是第i时刻的无量纲脉冲强度INB值,RMSi是第i时刻的RMS 值,RMSN是第i时刻的去心邻域的RMS平均值,N是邻域大小,取为N=1000。
Table 2 特征位置处的INB
Table 2为转捩/再次层流化处两种数据处理方法的INB值对比,在两个特征位置处,相较于滑移窗口技术,可变滑移窗口技术得到的INB值有量级上的增大。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (7)
1.一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中脉动压力均方根值(RMS)使用式1进行计算
本技术中构造一种可变的矩形窗形式,如式2所示,
W=h (2)
本技术中在窗函数增大的同时,窗口宽度变窄;在窗函数减小的同时,窗口宽度变宽,本技术中窗函数与窗口宽度的变化规律由式3进行控制:
m·h=m0 (3)
窗函数h由式4进行描述
翼型俯仰震荡规律为:
α=α0+A sin(2πft) (5)
由采样率fs引起的一个采样间隔内的角度误差Δαs为:
窗口内的角度误差即为:
Δα=2mΔαs=4mπ Af/fs cos(2πft) (7)
对于窗函数h,有以下性质:当RMS/RMS0值小于阈值μ0时,h随RMS/RMS0增大而缓慢增大,接近阈值μ0时,h迅速增大;当RMS/RMS0值达到阈值时,h取最大值;RMS/RMS0值大于阈值时,h迅速减小,由于半窗口宽度不能无限大,式3和式4决定了h有最小值。
根据h函数以上的性质特点,本技术中提出了一种窗口高度变化函数:类高斯函数,
式中,c1为窗函数的放大系数,c2为窗函数的形状系数,μ0为的阈值,c3为窗函数的平移量。
2.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中W(t)是窗函数,是采样时刻,是采样间隔,是窗口宽度,j是窗口内数据的序号,i表示第i个窗口,p是采集到的压力,表示第i个窗口内压力平均值,代表第i个窗口内压力的均方根值(RMS)。
3.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中的W(t)选择为矩形窗。
4.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中m是半窗口宽度,h是窗函数的值,m0是基准半窗口宽度。
5.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中RMS是窗口内的RMS值,RMS0是参考值。
6.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中α0是平均迎角,f是振荡频率,A是翼型震荡的振幅。
7.根据权利要求1所述的一种应用于动态翼型转捩判断的可变滑移窗口技术,其特征在于:所述技术中Ω是翼型俯仰运动的角速度,fs是采样率,Δt是采样率引起的采样时间间隔。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910054491.0A CN109740284B (zh) | 2019-01-21 | 2019-01-21 | 一种应用于动态翼型转捩判断的可变滑移窗口方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910054491.0A CN109740284B (zh) | 2019-01-21 | 2019-01-21 | 一种应用于动态翼型转捩判断的可变滑移窗口方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109740284A true CN109740284A (zh) | 2019-05-10 |
CN109740284B CN109740284B (zh) | 2020-09-22 |
Family
ID=66365473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910054491.0A Active CN109740284B (zh) | 2019-01-21 | 2019-01-21 | 一种应用于动态翼型转捩判断的可变滑移窗口方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109740284B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111832159A (zh) * | 2020-06-23 | 2020-10-27 | 北京临近空间飞行器***工程研究所 | 一种基于飞行试验数据的边界层转捩阵面动态演化过程确定方法 |
CN112328963A (zh) * | 2020-09-29 | 2021-02-05 | 国创新能源汽车智慧能源装备创新中心(江苏)有限公司 | 信号的有效值计算方法、装置 |
CN112748682A (zh) * | 2020-12-07 | 2021-05-04 | 珠海格力电器股份有限公司 | 信号采样处理方法、装置、设备和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130167096A1 (en) * | 2011-12-22 | 2013-06-27 | Lsi Corporation | Location and timing window based decoupling capacitor evaulation tool and method |
CN103441762A (zh) * | 2013-09-09 | 2013-12-11 | 江南大学 | 一种基于Blackman窗三谱线插值测试ADC动态参数的方法 |
CN104368407A (zh) * | 2014-10-11 | 2015-02-25 | 哈尔滨工程大学 | 多面可视化流动加热实验装置 |
CN105005694A (zh) * | 2015-07-13 | 2015-10-28 | 天津大学 | 一种基于动态称重***的桥梁疲劳寿命频域分析方法 |
CN107908592A (zh) * | 2017-11-16 | 2018-04-13 | 中国航天空气动力技术研究院 | 一种高超声速风洞背景噪声分频分析方法 |
CN107956708A (zh) * | 2017-11-17 | 2018-04-24 | 浙江大学 | 一种基于快速谱峭度分析的泵潜在空化故障检测方法 |
CN109033525A (zh) * | 2018-06-27 | 2018-12-18 | 浙江大学 | 一种基于简化三方程转捩模型的高超声速转捩预测方法 |
-
2019
- 2019-01-21 CN CN201910054491.0A patent/CN109740284B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130167096A1 (en) * | 2011-12-22 | 2013-06-27 | Lsi Corporation | Location and timing window based decoupling capacitor evaulation tool and method |
CN103441762A (zh) * | 2013-09-09 | 2013-12-11 | 江南大学 | 一种基于Blackman窗三谱线插值测试ADC动态参数的方法 |
CN104368407A (zh) * | 2014-10-11 | 2015-02-25 | 哈尔滨工程大学 | 多面可视化流动加热实验装置 |
CN105005694A (zh) * | 2015-07-13 | 2015-10-28 | 天津大学 | 一种基于动态称重***的桥梁疲劳寿命频域分析方法 |
CN107908592A (zh) * | 2017-11-16 | 2018-04-13 | 中国航天空气动力技术研究院 | 一种高超声速风洞背景噪声分频分析方法 |
CN107956708A (zh) * | 2017-11-17 | 2018-04-24 | 浙江大学 | 一种基于快速谱峭度分析的泵潜在空化故障检测方法 |
CN109033525A (zh) * | 2018-06-27 | 2018-12-18 | 浙江大学 | 一种基于简化三方程转捩模型的高超声速转捩预测方法 |
Non-Patent Citations (1)
Title |
---|
魏斌斌等: "基于POD方法的高超声速边界层转捩判定方法", 《气体物理》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111832159A (zh) * | 2020-06-23 | 2020-10-27 | 北京临近空间飞行器***工程研究所 | 一种基于飞行试验数据的边界层转捩阵面动态演化过程确定方法 |
CN111832159B (zh) * | 2020-06-23 | 2023-08-29 | 北京临近空间飞行器***工程研究所 | 一种基于飞行试验数据的边界层转捩阵面动态演化过程确定方法 |
CN112328963A (zh) * | 2020-09-29 | 2021-02-05 | 国创新能源汽车智慧能源装备创新中心(江苏)有限公司 | 信号的有效值计算方法、装置 |
CN112748682A (zh) * | 2020-12-07 | 2021-05-04 | 珠海格力电器股份有限公司 | 信号采样处理方法、装置、设备和存储介质 |
CN112748682B (zh) * | 2020-12-07 | 2021-12-14 | 珠海格力电器股份有限公司 | 信号采样处理方法、装置、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109740284B (zh) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109740284A (zh) | 一种应用于动态翼型转捩判断的可变滑移窗口技术 | |
Dong et al. | A wind tunnel simulation of the mean velocity fields behind upright porous fences | |
Faller | Instability and transition of disturbed flow over a rotating disk | |
Vickery et al. | Lift or across-wind response of tapered stacks | |
Bandi et al. | Experimental investigation on aerodynamic characteristics of various triangular-section high-rise buildings | |
CN105719029B (zh) | 一种基于风速波动特征提取的风电功率组合预测方法 | |
Fox et al. | A profile of the four moment measures perpendicular to a shoreline, South Haven, Michigan | |
Butterfield et al. | Comparison of wind tunnel airfoil performance data with wind turbine blade data | |
CN104537198B (zh) | 一种气体流量计的数据检测方法 | |
CN108509718A (zh) | 一种基于质量守恒的远场尾流二维解析模型 | |
Holst et al. | Experimental analysis of a NACA 0021 airfoil under dynamic angle of attack variation and low Reynolds numbers | |
Radok | Deposition and erosion of snow by the wind | |
CN109682570A (zh) | 一种非对称尾部形状水翼的水力阻尼比识别方法 | |
Narasimha et al. | Turbulent bursts in the atmosphere | |
Wong et al. | Flow characteristics around two parallel adjacent square cylinders of different sizes | |
Soltani et al. | Roughness and turbulence effects on the aerodynamic efficiency of a wind turbine blade section | |
CN104678992B (zh) | 矿用主通风机稳定性与经济性实时检测方法 | |
CN106092178A (zh) | 提高测量精度的数据修正方法 | |
Emeis | How well does a power law fit to a diabatic boundary-layer wind profile | |
Wang et al. | Turbulence modelling of deep dynamic stall at low Reynolds number | |
CN103267987A (zh) | 季风涡旋定位方法 | |
Fitton et al. | Scaling of turbulence in the atmospheric surface-layer: Which anisotropy? | |
Rooney et al. | Finite aspect ratio effects on vortex shedding behind two cylinders at angles of incidence | |
Wacker | Local Wind Pressures for Rectangular Buildings in Turbulent Boundary Layers | |
Yang et al. | Aerodynamic forces on two-dimensional rectangular cylinders subjected to accelerating flow: effect of side ratio |
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 |