基于传递矩阵法的裂纹叶片振动特性分析.docxVIP

基于传递矩阵法的裂纹叶片振动特性分析 1 裂纹叶片的振动特性 裂缝打开机制的描述。以往的研究, 基本上都将裂纹认为开口裂纹。振动过程中裂纹的开合现象, 使裂纹截面刚度呈现出时变的特征。任一瞬时的裂纹截面弯曲刚度可表示为k(ωt) =k-λ(ωt) [k-k1]。其中k为裂纹全闭状态时 (即无裂纹) 的截面弯曲刚度,k1为裂纹全开状态时的截面弯曲刚度,λ(ωt) 为描述裂纹开合机制的开合函数。本文选用方波函数作为开合函数, 即 λ(ωt)={00＜ωt＜π1π＜ωt＜2π 方波函数表明, 裂纹在振动的前半周处于全闭状态, 后半周则全开。 裂纹叶片的振动微分方程。由于裂纹的开合作用, 致使裂纹处的横截面惯性矩J和面积A不再是常数。据此, 可导得无扭曲裂纹悬臂梁的振动微分方程 ?[EJ(x,y,t)?2y/?x2]/?x2+ρA(x,y,t)?2y/?x2=0(1) 裂纹截面的等效刚度。根据平面应力状态的应变能理论, 可导出矩形裂纹截面的等效刚度为 kz=BH2E/72πr2(0.6384?1.035r+3.7201r2?5.1773r3+7.553r4?7.332r5+2.4909r6)(2) 式中:B为截面宽度,H为截面高度,a为裂纹深度,r=a/H为相对裂纹深度,E为弹性模量。裂纹截面的相对等效柔度 (有无裂纹的截面柔度之比)Cˉˉˉz=EJ/kzL。式中,L为叶片长度。 应用传递矩阵法求解裂纹叶片的振动特性。应用传递矩阵法求解裂纹叶片振动特性时, 要把叶片分成许多等截面段 (叶片段) , 凡是裂纹所在位置必须取为叶片段截面。 裂纹叶片的固有频率和强迫振动响应。固有频率以及均布 (沿叶高) 正弦型扰力f(x,t) =f(x) sinωt作用下的裂纹叶片强迫振动响应的求解过程与无裂纹叶片算法基本相同, 不同之处在于总传递矩阵中应包含裂纹截面处的段间传递矩阵, 此时, 裂纹截面可看作弹性铰链形式的结构单元, 弹性铰链的力矩刚性系数取为Ch=1/Cˉˉˉz。开合函数为方波函数的裂纹叶片的振动历程和频谱分析。由于裂纹的开合作用, 正弦型扰力激励的裂纹叶片振动历程 (时域信号) 将呈现前、后半周的非对称性。非对称性的程度取决于裂纹的开合过程 (开合函数) 和裂纹深度等因素。若开合函数为方波函数, 则时域信号的前、后半周仍各为半个正弦波, 但其振幅不同, 从而形成整个周期内的非正弦性。应用前述强迫振动响应的求解方法, 可计算出前、后半周的振幅比。将这种前、后半周不同振幅半个正弦波的周期信号展成傅氏级数, 得: y(t)=(1/π)(A1?A2)+(1/2)(A1+A2)sinωt+(2/π)(A2?A1)[(1/3)cos2ωt+(1/15)cos4ωt+(1/35)cos6ωt+(1/63)cos8ωt+??] 式中:A1和A2为前、后半周的振幅。由上式看出, 在正弦型扰力作用下, 开合函数为方波函数的裂纹叶片的振动信号中, 除基波成份外, 还含有偶数项谐波的高阶谐波分量, 而且, 随着前、后半周振幅差的增加, 高阶振动分量加大。 旋转裂纹叶片的固有频率。计算旋转裂纹叶片固有频率时, 需将各叶片段中央以上叶身部分产生的离心轴向力, 附加到各叶片段上, 从而导出各叶片段的传递矩阵。 计算实例。应用文中建立的方法, 计算了一钢制等截面矩形裂纹悬臂梁的有关振动特性。梁的几何尺寸:长度L=0.1 m, 宽度B=0.012 m, 高度H=0.0025 m。 图1给出一阶弯曲固有频率相对值ω/ωno(ωno为无裂纹叶片的一阶弯曲固有频率) 随裂纹的相对深度r、裂纹的相对位置S(裂纹至叶根距离/叶片长度) 的变化关系曲线。曲线表明:固有频率随裂纹深度的增加而下降 (在小裂纹深度范围内, 呈现线性变化) ;裂纹深度相同时, 固有频率将随裂纹位置由叶根向叶尖的移动而上升。 后、前半周振幅比A2/A1与裂纹相对深度间的关系示于图2。由图得知, 在扰力 (幅值和频率) 一定的情况下,A2/A1随r的增加而加大, 即裂纹深度增加时, 将加大前、后半周的非对称性, 这表明叶片振动信号中的高阶谐波分量, 随裂纹深度的增大而加大。 图3给出不同裂纹深度下一阶弯曲固有频率 (动频) 随转速的变化曲线。其特征是, 在高转速下, 不同裂纹深度下的动频趋于无裂纹叶片的动频。造成这种趋势的原因在于:随着裂纹深度的增加, 离心力对固有频率的影响将加剧。但是, 随着离心力 (转速) 的增加, 裂纹的影响将减弱 (相对于离心力的影响) 。 2 疲劳裂纹测试 对计算实例的等截面裂纹悬臂梁进行了振动特性实验。实验前, 分别在4个试件的不同位置 (沿叶高) 上, 预制出一条横向疲劳裂纹 (非缝隙裂纹) 。预制过程是, 事先用一尖锐刀具轻微刻划一尖劈缺口, 然后将试件安装在振动台上振动,