首页 / 专利库 / 物理 / 拉梅常数 / 一种利用SH波检测旋翼结冰冰层厚度的方法

一种利用SH波检测旋翼结冰层厚度的方法

阅读:472发布:2020-07-08

专利汇可以提供一种利用SH波检测旋翼结冰层厚度的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种利用SH波检测旋翼结 冰 冰层厚度的方法,利用SH波的频散特性,得出旋翼上下表面 覆盖 有不同厚度冰层时较厚层的冰层厚度,包括:建立物理模型及控制方程;根据控制方程,利用数值分析 软件 分析控制方程;计算得出一阶模态相速度只与冰层厚层的厚度值有关的区间;利用上述区间,根据SH波一阶模态频散特性确定冰层厚层的厚度值;根据上下冰层表面位移关系确定冰层厚薄分布情况。本发明能够实现利用SH波对旋翼上下表面覆盖不同厚度冰层时较厚冰层厚度值的检测,为旋翼结冰 传感器 的设计提供理论依据。,下面是一种利用SH波检测旋翼结冰层厚度的方法专利的具体信息内容。

1.一种利用SH波检测旋翼结冰层厚度的方法,其特征在于,包括步骤如下:
1)建立物理模型及控制方程;
2)利用控制方程计算得出一阶模态相速度与冰层厚层的厚度值相关的区间;
3)利用上述得到的区间,根据SH波频率与上述得到的一阶模态相速度值确定冰层厚层的厚度值;
4)根据上下冰层表面关于旋翼对称处点的一阶模态位移幅值大小关系确定冰层厚层与冰层薄层在旋翼上下表面的分布情况。
2.根据权利要求1所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,上述步骤1)具体包括:根据旋翼上下表面覆盖的冰层建立三层板物理模型;根据SH波的波动方程与边界条件得到关于上下冰层层厚、频率和波速的控制方程。
3.根据权利要求1所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,上述步骤1)具体包括:旋翼结冰时,上下表面均出现结冰现象,即上下表面将有冰层覆盖,出现一个表面的冰层厚度大于或等于另一个表面的冰层厚度;据此建立物理模型为三层板结构,以板厚中点为o点,板长方向为x1轴,板厚方向为x2轴,垂直于板向外方向为x3轴,其中,中间层为旋翼,上下层为冰层,冰层与旋翼刚性连接,冰层表面为自由表面;
SH波的波动方程为:
其中,u3为x3方向位移,cTi(i=1,2,3)为层中的体剪切波波速: μi为拉梅常数,ρi为各材料密度;由SH波的波动方程可知,各层中x3方向的位移表示为:
u3i=[Ai sin(qix2)+Bi cos(qix2)]exp[i(kx1-ωt)]
其中, k为波数:k=ω/c,c为沿x方向行波速度,ω为声波频率;
Ai、Bi为未知常数;
利用位移和应的关系 表示出应力,再利用边界条件得到关于Ai、Bi为未知数的齐次线性方程组:
若要求此线性方程组有非零解,则其系数行列式为零,即|amn|6×6=0,其中amn元素为线性方程组的系数项;由于旋翼和冰层的材料参数已知,旋翼厚度为已知,可知控制方程:f(h1,h2,f,c)=0只与上下表面冰层厚度h1、h2、频率f和波速c有关。
4.根据权利要求3所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,对上述的控制方程f(h1,h2,f,c)=0的分析通过固定其中任意两个变量,利用MATLAB计算得到剩余两个变量的关系曲线,通过所得到的关系曲线来判断变量之间的关系。
5.根据权利要求3所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,上述步骤2)具体包括:首先根据控制方程f(h1,h2,f,c)=0,利用MATLAB计算求出不同h1、h2值下一阶模态频散关系曲线,从频散关系曲线中得到,当SH波波长大于旋翼厚度时,一阶模态相速度受到上下层冰层的共同影响或者出现上下层冰层厚度变化不影响一阶模态相速度的情况,对于这两种情况都无法用来判断冰层厚度;因此,用于检测冰层厚度的SH波波长在选择时应小于旋翼厚度;
当SH波波长小于旋翼厚度时,根据控制方程f(h1,h2,f,c)=0,固定频率f,得到h1取不同值时,h2与一阶模态相速度的关系曲线;从h1和h2对一阶模态相速度值c的影响得出,max(h1,h2)值对一阶模态相速度值的影响明显占优,而min(h1,h2)值对一阶模态相速度值影响微弱,并且在max(h1,h2)值固定时,随着频率的增加,min(h1,h2)值对一阶模态相速度值影响减小;
对于控制方程f(h1,h2,f,c)=0,若max(h1,h2)值固定,由对称性,假设固定h1,每一个频率f下,都能求出与一阶模态相速度的关系曲线,即f(h1,h2,f,c)=0;当频率值为f0时,一阶模态相速度中最大值与最小值之差δ等于某一预设小值ε0时,此时忽略冰层薄层对一阶模态相速度的影响,并且由于频率越大,δ越小,使得在大于频率值f0时,一阶模态相速度认为只受到冰层厚层的厚度值影响;当冰层厚层的厚度值连续取值时,绘制出一条曲线,其横坐标为上述频率值,纵坐标为冰层薄层对一阶模态相速度影响忽略时的一阶模态相速度;并且处于该曲线下方的频率与一阶模态相速度点只与冰层厚层的厚度对应,而与冰层薄层的厚度无关;
根据f(h1,h2,f,c)=0,得出c=g(h1,h2,f),已知当f≥f0时,覆盖冰层薄层的厚度值对一阶模态相速度影响忽略,则c=g(h1,f),该函数即一阶模态相速度只与冰层厚层的厚度值有关的区间。
6.根据权利要求3所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,上述步骤4)具体包括:由控制方程f(h1,h2,f,c)=0得到的位移的表达式u=f1(h1,h2,f);利用MATLAB求出在每一个频率下两个厚度值取任意值时的位移情况;
根据u=f1(h1,h2,f),通过对h1、h2在各频率下不同的取值,得到不同h1、h2取值时上下冰层表面同一x1坐标处点的一阶模态位移幅值的大小关系;该幅值大小关系对应于厚度关系,即该位移幅值大的那层覆盖冰层厚度大,而位移幅值小的那层覆盖冰层厚度小,若位移幅值相等,则上下表面覆盖冰层厚度也相等。
7.根据权利要求1所述的利用SH波检测旋翼结冰冰层厚度的方法,其特征在于,上述步骤
3)具体包括:根据检测时SH波的中心频率与所得到的一阶模态相速度值,若属于上述区间,则满足此时一阶模态相速度与冰层厚层的厚度的对应关系,即
若不属于上述区间,则此时一阶模态相速度受到厚层厚度与薄层厚度共同影响,无法准确判断覆盖冰层厚度,应当增加中心频率,使得频率与一阶模态相速度值处于上述区间内;若所得到的一阶模态相速度等于冰层中体剪切波波速时,可知此时频率过高,应降低检测时SH波的中心频率。

说明书全文

一种利用SH波检测旋翼结冰层厚度的方法

技术领域

[0001] 本发明属于无损检测技术领域,特别涉及一种利用SH波检测旋翼结冰冰层厚度的方法。

背景技术

[0002] 自1970年以来,当直升机被允许在仪表飞行规则下进行飞行后,直升机飞行遇到了越来越多的层结冰问题。结冰对直升机性能的影响主要体现在旋翼结冰问题上,旋翼结冰后会对飞行产生不利的影响,如旋翼损伤甚至机毁人亡,因此对于冰层厚度的检测显得尤为重要。目前的技术仍无法避免飞机旋翼结冰问题,但是通过检测飞机旋翼结冰后冰层的厚度,尤其是较厚的冰层厚度,在结冰初期进行预警,或冰层达到一定厚度时发出警告,让飞行员知道旋翼冰层对飞行的影响,从而做出操作判断,避免发生重大事故是非常有意义的。
[0003] 结冰检测传感器的发展,研究者已经提出了多种检测方法:光学法、热学法、电学法、波导法等,由于声表面波检测方法具有灵敏度高,可无线无源等优点,成为较为理想的飞机结冰检测的方法。但目前研究仅限于对旋翼单侧覆盖有冰层的厚度检测,显然与旋翼上下表面均会出现结冰现象且结冰厚度并不一定相等的实际工况不相符,并且无法对达到一定厚度的冰层进行预警,导致无法满足对旋翼结冰状况进行准确评估。

发明内容

[0004] 针对于上述现有技术的不足,本发明的目的在于提供一种利用SH波检测旋翼结冰冰层厚度的方法,以解决现有技术研究仅限于对旋翼单侧覆盖有冰层的厚度检测,与旋翼上下表面均会出现结冰现象且结冰厚度并不一定相等的实际工况不相符,无法对达到一定厚度的冰层进行预警,导致无法满足对旋翼结冰状况进行准确评估的问题。
[0005] 为达到上述目的,本发明的一种利用SH波检测旋翼结冰冰层厚度的方法,包括步骤如下:
[0006] 1)建立物理模型及控制方程;
[0007] 2)利用控制方程计算得出一阶模态相速度与冰层厚层的厚度值相关的区间;
[0008] 3)利用上述得到的区间,根据SH波频率与上述得到的一阶模态相速度值确定冰层厚层的厚度值;
[0009] 4)根据上下冰层表面关于旋翼对称处点的一阶模态位移幅值大小关系确定冰层厚层与冰层薄层在旋翼上下表面的分布情况。
[0010] 优选地,上述步骤1)具体包括:根据旋翼上下表面覆盖的冰层建立三层板物理模型;根据SH波的波动方程与边界条件得到关于上下冰层层厚、频率和波速的控制方程。
[0011] 优选地,上述步骤1)具体包括:旋翼结冰时,上下表面均出现结冰现象,即上下表面将有冰层覆盖,出现一个表面的冰层厚度大于或等于另一个表面的冰层厚度;据此建立物理模型为三层板结构,以板厚中点为o点,板长方向为x1轴,板厚方向为x2轴,垂直于板向外方向为x3轴,其中,中间层为旋翼,上下层为冰层,冰层与旋翼刚性连接,冰层表面为自由表面;
[0012] SH波的波动方程为:
[0013]
[0014] 其中,u3为x3方向位移,cTi(i=1,2,3)为层中的体剪切波波速: μi为拉梅常数,ρi为各材料密度;由SH波的波动方程可知,各层中x3方向的位移表示为:
[0015] u3i=[Ai sin(qix2)+Bi cos(qix2)]exp[i(kx1-ωt)]
[0016] 其中, k为波数:k=ω/c,c为沿x方向行波速度,ω为声波频率;Ai、Bi为未知常数;
[0017] 利用位移和应的关系 表示出应力,再利用边界条件得到关于Ai、Bi为未知数的齐次线性方程组:
[0018]
[0019] 若要求此线性方程组有非零解,则其系数行列式为零,即|amn|6×6=0,其中amn元素为线性方程组的系数项;由于旋翼和冰层的材料参数已知,旋翼厚度为已知,可知控制方程:f(h1,h2,f,c)=0只与上下表面冰层厚度h1、h2、频率f和波速c有关。
[0020] 优选地,对上述的控制方程f(h1,h2,f,c)=0的分析通过固定其中任意两个变量,利用MATLAB计算得到剩余两个变量的关系曲线,通过所得到的关系曲线来判断变量之间的关系。
[0021] 优选地,上述步骤2)具体包括:首先根据控制方程f(h1,h2,f,c)=0,利用MATLAB计算求出不同h1、h2值下一阶模态频散关系曲线,从频散关系曲线中得到,当SH波波长大于旋翼厚度时,一阶模态相速度受到上下层冰层的共同影响或者出现上下层冰层厚度变化不影响一阶模态相速度的情况,对于这两种情况都无法用来判断冰层厚度;因此,用于检测冰层厚度的SH波波长在选择时应小于旋翼厚度;
[0022] 当SH波波长小于旋翼厚度时,根据控制方程f(h1,h2,f,c)=0,固定频率f,得到h1取不同值时,h2与一阶模态相速度的关系曲线;从h1和h2对一阶模态相速度值c的影响得出,max(h1,h2)值对一阶模态相速度值的影响明显占优,而min(h1,h2)值对一阶模态相速度值影响微弱,并且在max(h1,h2)值固定时,随着频率的增加,min(h1,h2)值对一阶模态相速度值影响减小;
[0023] 对于控制方程f(h1,h2,f,c)=0,若max(h1,h2)值固定,由对称性,假设固定h1,每一个频率f下,都能求出与一阶模态相速度的关系曲线,即f(h1,h2,f,c)=0;当频率值为f0时,一阶模态相速度中最大值与最小值之差δ等于某一预设小值ε0时,此时忽略冰层薄层对一阶模态相速度的影响,并且由于频率越大,δ越小,使得在大于频率值f0时,一阶模态相速度认为只受到冰层厚层的厚度值影响;当冰层厚层的厚度值连续取值时,绘制出一条曲线,其横坐标为上述频率值,纵坐标为冰层薄层对一阶模态相速度影响忽略时的一阶模态相速度;并且处于该曲线下方的频率与一阶模态相速度点只与冰层厚层的厚度对应,而与冰层薄层的厚度无关;
[0024] 根据f(h1,h2,f,c)=0,得出c=g(h1,h2,f),已知当f≥f0时,覆盖冰层薄层的厚度值对一阶模态相速度影响忽略,则c=g(h1,f),该函数即一阶模态相速度只与冰层厚层的厚度值有关的区间。
[0025] 优选地,上述步骤3)具体包括:根据检测时SH波的中心频率与所得到的一阶模态相速度值,若属于上述区间,则满足此时一阶模态相速度与冰层厚层的厚度的对应关系,即若不属于上述区间,则此时一阶模态相速度受到厚层厚度与薄层厚度共同影响,无法准确判断覆盖冰层厚度,应当增加中心频率,使得频率与一阶模态相速度值处于上述区间内;若所得到的一阶模态相速度等于冰层中体剪切波波速时,可知此时频率过高,应降低检测时SH波的中心频率。
[0026] 优选地,上述步骤4)具体包括:由控制方程f(h1,h2,f,c)=0得到的位移的表达式u=f1(h1,h2,f);利用MATLAB求出在每一个频率下两个厚度值取任意值时的位移情况;
[0027] 根据u=f1(h1,h2,f),通过对h1、h2在各频率下不同的取值,得到不同h1、h2取值时上下冰层表面同一x1坐标处点的一阶模态位移幅值的大小关系;该幅值大小关系对应于厚度关系,即该位移幅值大的那层覆盖冰层厚度大,而位移幅值小的那层覆盖冰层厚度小,若位移幅值相等,则上下表面覆盖冰层厚度也相等。
[0028] 本发明的有益效果:
[0029] 本发明利用SH波的频散特性,对旋翼上下表面覆盖不同厚度冰层时较厚冰层厚度值的检测,为旋翼结冰传感器的设计提供理论依据。从而实现利用SH波检测旋翼上下表面均覆盖冰层时较厚一层厚度的目的。附图说明
[0030] 图1为本发明实施例提供的利用SH波检测旋翼结冰冰层厚度方法的物理模型图。
[0031] 图2为本发明实施例提供的波长与旋翼厚度相等时厚度与一阶模态的相速度关系图。
[0032] 图3为本发明实施例提供的旋翼结冰不同冰层厚度下一阶模态频散关系图。
[0033] 图4为本发明实施例提供的旋翼结冰忽略薄层影响临界频率与一阶模态相速度关系图。
[0034] 图5为本发明实施例提供的频率为200KHz冰层厚层的厚度值与一阶模态相速度关系图。
[0035] 图6为本发明实施例提供的旋翼结冰任意冰层厚度时上下表面一阶模态位移幅值比值图。

具体实施方式

[0036] 为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
[0037] 参照图1所示,本发明的一种利用SH波检测旋翼结冰冰层厚度的方法,包括步骤如下:
[0038] 1)建立物理模型及控制方程;根据旋翼上下表面覆盖的冰层建立三层板物理模型;根据SH波的波动方程与边界条件得到关于上下冰层层厚、频率和波速的控制方程;
[0039] 2)利用控制方程计算得出一阶模态相速度与冰层厚层的厚度值相关的区间;
[0040] 3)利用上述得到的区间,根据SH波频率与上述得到的一阶模态相速度值确定冰层厚层的厚度值;
[0041] 4)根据上下冰层表面关于旋翼对称处点的一阶模态位移幅值大小关系确定冰层厚层与冰层薄层在旋翼上下表面的分布情况。
[0042] 具体实施例中:1、选取如图1所示物理模型,旋翼为金属基层,厚度2h=10mm,材料参数为:μ=54.675GPa,ρ=2700kg/m3;上下层为旋翼结冰时覆盖的冰层,厚度分别为h1和3
h2,材料参数为:μ=54.675GPa,ρ=2700kg/m,坐标系选取如图所示。
[0043] SH波的波动方程为:
[0044]
[0045] 其中,u3为x3方向位移,cTi(i=1,2,3)为层中的体剪切波波速: μi为拉梅常数,ρi为各材料密度;由SH波的波动方程可知,各层中x3方向的位移表示为:
[0046] u3i=[Ai sin(qix2)+Bi cos(qix2)]exp[i(kx1-ωt)]
[0047] 其中, k为波数:k=ω/c,c为沿x方向行波速度,ω为声波角频率;Ai、Bi为未知常数;
[0048] 利用位移和应力的关系 表示出应力,再利用边界条件得到关于Ai、Bi为未知数的齐次线性方程组:
[0049]
[0050] 若要求此线性方程组有非零解,则其系数行列式为零,即|amn|6×6=0,其中amn元素为线性方程组的系数项;由于旋翼和冰层的材料参数已知,旋翼厚度为已知,可知控制方程:f(h1,h2,f,c)=0只与上下表面冰层厚度h1、h2、频率f和波速c有关,即f(h1,h2,f,c)=0,将参数带入控制方程,通过固定任意两个变量,利用MATLAB分析剩余两个变量的关系。
[0051] 2、取SH波波长λ与基层2h厚度相等 即H=1,当基层只有单面覆盖冰层时,即h2=0时,可以得出冰层厚度与一阶模态相速度关系图,如图2中o点所示;当基层上下表面都覆盖冰层时,取h2=0.1,同样可以得出冰层厚度与一阶模态相速度关系图,如图2中·点所示。从图2中可以看出,当基层上下表面均有冰层覆盖时,在SH波波长与金属基层厚度相等时,一阶模态相速度表现出只受冰层厚层的厚度值的影响。
[0052] 3、将h1=2、4、6、8(mm),h2=0、2、4、6、8(mm)分别对应取值,同时保证h1>=h2,根据控制方程,绘制出不同厚度冰层时频散曲线后得到图3。从图3中我们可以得出,频率在0-50KHz时,一阶模态相速度受到上下冰层厚度的共同影响,因此检测用SH波频率应大于
50KHz。
[0053] 冰层厚层的厚度值可能出现的4个取值将曲线划分为4簇,即对一阶模态相速度值的影响起到了明显占优的地位,每1簇曲线由于冰层薄层的厚度值的变化使得一阶模态相速度值发生微小变化,同时表现出随着频率的增加,将逐渐趋向于薄层厚度为零的曲线,即在频率超过某值后,冰层薄层的厚度值对于相速度的影响可忽略。
[0054] 4、h1固定时,h2=0~h1,频率值的取值逐渐增大,计算出h2变化中得到的一阶模态相速度值的最大值与最小值之差等于ε0时(本例中ε0取值为10)的频率值,同时记录下此时一阶模态相速度最大值的数值c0;将h1在0~8内取值时得到的所有(f0,c0)点绘制于同一张图中得到图4。
[0055] 图4所表明的意义:在图中曲线上及其下方的频率与一阶模态相速度关系,只与冰层厚层的厚度值有关,而与冰层薄层的厚度值无关。
[0056] 5、利用中心频率为200KHz的SH波对厚度为10mm,上下表面都覆盖有冰层的旋翼进行检测,若得出一阶模态相速度为2036m/s。首先根据图4可知,点(200,2036)处于曲线下方,即满足相速度只与冰层厚层的厚度值有关。
[0057] 图5为f=200KHz时,2h=10mm,h2=0时,h1与一阶模态相速度关系曲线;从图5中可以知道一阶模态相速度为2036m/s所对应的h1为4mm。由此可以确定此时旋翼结冰覆盖冰层厚层的厚度为4mm。
[0058] 若经过信号处理后得到的相速度点处于图4曲线上方,则可知此时一阶模态相速度受到旋翼上下冰层厚度共同影响,无法做出准确判断,此时应当适当增加SH波中心频率,使得得到的相速度点处于图4曲线下方。
[0059] 若经过信号处理后得到的一阶模态相速度为1729m/s,此时相速度已经为冰层中体剪切波波速,此时应当降低频率。
[0060] 6、图6为不同频率值时,h1,h2随机取值,并且始终保持h1>=h2,同一x1坐标下,旋翼下方冰层下表面一阶模态位移幅值与旋翼上方冰层上表面一阶模态位移幅值的比值,即纵坐标0~1表示下表面幅值小于等于上表面幅值。由图6可知,该比值保持在0~1之间,说明与h1,h2的具体取值无关,即位移幅值大的冰层厚度也大;位移幅值小的冰层厚度也小;位移幅值相等则冰层厚度也相等。
[0061] 综合上述的步骤,可以确定旋翼结冰时上下冰层中较厚一层的厚度及分布,从而实现利用SH波检测旋翼结冰冰层厚度。
[0062] 本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈