双通道旋转输流管临界流速和振动模态分析1)
时间：2023-04-08 20:10:05 来源：柠檬阅读网 本文已影响 人 
张 博 郑昊楷 孙东生 丁 虎 陈立群,3)
* (长安大学理学院,西安 710064)
† (上海大学力学与工程科学学院,上海 200444)
航空发动机是飞机的心脏,是国家工业基础、科技水平和国防实力的重要标志.涡轮叶片是航空发动机中最重要的关键部件之一,其可靠性直接影响发动机能否正常.已有研究表明[1-4],航空发动机故障原因中62%以上是由于涡轮叶片振动过量,且由于高温、高压工况下通过蛇形冷却通道进行降温的设计,叶片通道内部受到科氏力、离心力等因素的相互作用,具有复杂的流固耦合效应.因此,研究含蛇形冷却通道的旋转叶片动力学特性,对发动机叶片的合理化设计与减振降噪具有重要意义.
由于旋转叶片结构在工程中的广泛应用,尤其是在航空发动机领域的核心地位,从上世纪四50 年代以来,旋转叶片的动力学行为研究吸引了国内外学者的广泛关注.如Carnegie[5]在1959 年首次研究了旋转叶片的自振频率.张伟等[1]将叶片简化为功能梯度材料薄壁梁,通过伽辽金法研究了系统在空气热弹性载荷下的动力学行为.此外,许多学者[6-11]将叶片简化为大变形柔性薄板.如赵飞云等[12]基于连续介质力学理论和Jourdain 速度变分原理,研究高速旋转下的柔性矩形薄板的耦合动力学模型.郑彤等[13]、蒋建平等[14],利用Langrage 方程法得到系统的一次近似耦合动力学方程.
以上研究,将旋转叶片简化为旋转的悬臂梁或旋转的大变形柔性薄板,但均未考虑内部流体对叶片动力学特性的影响.由前人的研究[15-21]可知内部流体会对叶片的振动产生复杂的影响,将叶片简化为旋转悬臂输液管模型可以从宏观上更深入地研究流体的某一特性对旋转叶片的动力学影响.Chen[22]在1971 年对输送脉动流体的管道振动进行研究,并通过数值方法确定稳定性边界.Paidoussis[23]在1987 年对输液管的振动作了详细的论述,并指出发散失稳和颤振失稳两种值得深入研究的现象.后续有Rousselet 等[24]、Lundgren 等[25-27]和徐鉴等[28-29]先后导出了悬臂输液管的非线性运动方程.
Paidoussis 等[30]、黄玉盈等[31]和Ibrahim[32-33]分别在1993 年、1998 年和2010 年对输液管模型动力学模型的研究现状进行了详细的综述.王乙坤等[34]研究了脉动内流作用下的输液管振动问题.Wang等[35]用牛顿法和龙格-库塔法对波纹输液管的动力学行为进行理论和数值分析,发现波纹总数和波纹幅度对输液管的动力学行为有较大影响.宫亚飞等[36]利用哈密顿原理和伽辽金法推导了含有初始弯曲的功能梯度输液管的非线性动力学方程,分析超临界功能梯度输液管的自由振动特性.Oh 等[37]利用瑞丽利兹法推导含冷却通道的旋转叶片的热弹动力耦合模型,研究冷却效应下系统的固有频率与拉伸特性.张博等[38]基于Langrage 原理和假设模态法研究不同端部集中质量、转速对旋转输液管临界流速的影响及特定参数组合下的内共振现象.
通过文献调研不难发现,在内冷叶片研究领域,考虑内流作用对系统振动特性影响的报道较少,文献[38]中虽开展了相关研究,但只局限于单冷却通道模型.为了更贴合包含多条冷却通道且通有流体介质的先进内冷叶片工程实例[39],本文在文献[38]基础上,引入叶片换热特性研究[40]中对内部通道的简化方式,将叶片简化为含双输液通道的旋转悬臂输液管,采用能量法建立其动力学模型,揭示包含多条冷却通道的旋转叶片结构中,冷却通道内流体流动对叶片动力学特性的影响规律.
本文将含蛇形冷却通道的涡轮叶片简化为如图1 所示长为L的两通道旋转悬臂输流管,固接在半径r的中空圆柱上以常数Ω旋转,截面高为h,宽为b,单位长度质量为m,两通道半径均为r0,两通道轴线与输流管轴线距离均为d,两通道内流体的单位长度质量分别为M1,M2,相对通道的流速分别为常数U1与U2.以中空圆柱中心O点为原点,建立全局坐标系XYZ,单位方向向量分别为i,j,k.以输流管中轴线与圆柱侧面的交点为o点,沿轴线方向为x轴,建立随转坐标系xyz,单位方向向量分别为i′,j′,k′,轴线上任意一点在t时刻的x,y方向位移分量分别用w1(x,t)和w2(x,t)表示,由变形引起的转角为θ.为了简化分析,引入以下假设: (1)忽略输流管剪切变形和转动惯量的影响;(2)忽略输流管在z方向的变形,即振动发生在XY平面内;(3)两通道内的流体均为为定常不可压缩的无黏流体;(4)输流管是均匀、各向同性的线弹性材料.
图1 含两通道的旋转输流管Fig.1 Sketch of a rotating cantilever double channel pipe
在全局坐标系下,输流管上任意一点的速度矢量为
令y=0 可表示输流管中轴线上任意一点的速度矢量为
则输流管的动能为
输流管因旋转运动而产生的轴向收缩势能[41]为
输流管的变形势能为
式中,EI为输流管的抗弯刚度,Ap为截面面积.在全局坐标系下,两通道内流体的速度分别为
两通道内流体的动能为
根据Euler-Bernoulli 梁理论,在随转坐标系下,通道1(y=d)和通道2(y=−d)中轴线上任意一点的位置矢量为
在输流管开放系统中,输流管做无穷小变形运动下,流体对输流管做的功可由文献[42]得
则旋转输流管开放系统的Lagrange 函数为
本文两通道旋转悬臂输流管的边界条件为
应用假设模态法离散系统的两个变形分量
式中,qi(t) 代表系统在相应方向的广义位移,N1,N2为选取的模态函数个数,ϕi(x)为满足系统位移边界条件的模态函数
其中参数βl为频率方程 c os(βlL)ch(βlL)=−1的根.
根据Lagrange 原理,得到系统矩阵形式的动力学方程
为了使研究结果具有普遍性,定义以下无量纲参数
其中,α0,α1,α2分别表示输流管及两通道内流体的单位长度质量权重;τ,κ,Ω*分别表示无量纲时间、长细比与无量纲转速;u1和u2分别表示两通道中流体的无量纲流速.
得到无量纲化后的系统动力学方程为
其中质量子阵、阻尼子阵、陀螺子阵、刚度子阵和广义力子阵见附录.
采用降阶扩维的方法求解系统的特征值,为简化方程,引入以下向量与矩阵
由于本文研究的是自由振动问题,方程可简化为
设其通解为
将其代入式(22),得到关于实值矩阵的一般特征值问题
求解得到系统呈复共轭出现的特征根和特征向量
式中,σr为阻尼因子,ωr为固有频率.因此,特征矢量矩阵为
由式(16)、式(29)与式(30)可得输液管横向位移的第r阶模态响应为
式中单划线与双划线部分分别代表了阻尼因子和相位变化对输流管横向位移模态响应的影响,为表述方便,下文使用符号和分别表示上述两因素归一化模态.
本文的系统参数选取参考文献[38],在后文中,若无特别说明,系统的无量纲参数设置为:α0=0.675,α1=0.1625,α2=0.1625,γ=0.5,κ=61.94,Ω*=5,δ=0.012.
图2 研究了假设模态法的收敛性,在不同试探函数个数下,系统的特征根轨迹图相似,临界流速数值存在微小差别.当试探函数个数取15 时,临界流速收敛.在后文的数值算例中均取N1=N2=15.由图可知,当流体流速较小时,流体引起了系统各阶模态的阻尼效应.随着流速的增加,第1 和3 阶模态轨迹相继穿越横轴,预示着输流管1 阶和3 阶模态将发生颤振失稳,但第2 阶模态始终未发生失稳.
图2 不同试探函数个数下前3 阶特征根轨迹曲线 (Ω*=0)Fig.2 The trajectories of the first three eigenvalues for different trail function numbers (Ω*=0)
以文献[43]中两端简支的单通道输流管系统为研究对象,得到该系统的前两阶特征根轨迹曲线,如图3 所示.将所得的前两阶临界流速与文献[43]的结果相比较,如表1 所示,两者结果吻合,验证了本文建模方法的正确性.
表1 两端简支单通道输流管系统的前两阶临界流速对比Table 1 Comparison of the first two critical flow velocities for simply supported flow pipe
图3 两端简支的单通道输流管系统的前两阶特征根轨迹Fig.3 The trajectories of the first two eigenvalues for simple supported flow pipe
采用双通道模型更贴合含多条冷却通道的先进内冷叶片工程实际,在此基础上开展的动力学特性研究更具有参考价值.图4 对比了单位时间流量相同时,中心单通道模型与对称双通道模型的系统临界流速值.其中,单通道的截面半径为(r0为双通道模型的截面半径),两模型的长细比分别为κsingle=59.75,κtwin=74.05,其余参数与前文一致.由图可见,在不同转速下,双通道模型的临界流速值均大于单通道模型,且随着转速的增加,双通道模型的临界流速值也比单通道模型增长得更快.
图4 单、双通道模型的临界流速随转速的变化曲线Fig.4 The curves of critical flow velocity with rotating speed for the single/twin-channel models
图5 研究了无量纲流速u1和u2的比例分别为1:1,1:2,1:3,1:4 时系统的特征根轨迹图,其中图中所标识的流速代表通道2 中流体流速.可以发现在不同流速比值下,得到的系统特征根轨迹图曲线形状大致相同,随着比值减小,第1 阶模态和第3 阶模态在发生失稳时的临界流速值增大,第2 阶模态始终未发生失稳.
图5 不同流速比下前3 阶特征根轨迹(ucr=u2,Ω*=5)Fig.5 The trajectories of the first three eigenvalues for different velocity ratios (ucr=u2,Ω*=5)
从系统的动力学方程不难发现,通道内流体的流动与管道整体的旋转运动均会给系统引入陀螺效应,即本文建立的动力学模型为典型的双陀螺系统.图6 研究了转速对系统特征根轨迹的影响.由图6 可见,转速第1 阶模态的特征根轨迹影响较小,并始终只存在一个临界流速.而随着转速的提高,系统的第2,3 阶特征根轨迹绕圈现象越来越显著,并会多次穿越虚轴,预示流速变化过程中系统的第2,3 阶模态将经历失稳→稳定→失稳的复杂动力学历程,产生该现象是由于旋转运动引入的陀螺效应对系统响应的影响.为了进一步研究这一现象,控制转速不变,在图7 中绘制了系统特征根轨迹随长细比的演变规律.可以发现,随着κ增大,系统第2,3 阶模态特征根轨迹的绕圈现象逐渐消失,当κ足够大时,第2 阶模态不再发生失稳,第3 阶特征根轨迹只穿越一次虚轴,只存在一个临界流速,类似于经典的悬臂输流管模型[15].实际上,随着κ增大,旋转运动引起的管道轴向与横向间的陀螺耦合效应逐渐减弱[44],系统双陀螺效应中通道内流体的流动引起的陀螺效应成为主导因素,因而系统表现出类似非旋转的悬臂输流管的动力学行为.
图6 不同转速下前3 阶特征根轨迹曲线(κ=35)Fig.6 The trajectories of the first three eigenvalues for different rotating speeds (κ=35)
图7 不同长细比下前3 阶特征根轨迹曲线(Ω*=5)Fig.7 The trajectories of the first three eigenvalues for different slenderness ratios (Ω*=5)
图8 揭示了相位变化对输流管横向位移前3 阶模态响应的影响.由图可见,陀螺效应使得系统不同位置的响应出现相位差.与实模态系统振动呈现驻波特性不同,复模态系统振动呈现出行波特性,即出现“节点”移动现象.
图8 相位变化对前3 阶模态响应的影响(u1=u2=8,Ω*=10)Fig.8 The curves of the first three modal responses for different phase changes (u1=u2=8,Ω*=10)
图9(a)～ 图9(c)给出了3 组不同流速和转速组合下,阻尼因子对输流管横向位移前3 阶模态响应的影响.当u1=u2=6,Ω*=15 时,前3 阶模态响应的幅值均随时间减弱;当u1=u2=10,Ω*=12 时,前3 阶模态响应的幅值随时间分别增强、减弱、减弱;当u1=u2=15,Ω*=15 时,前3 阶模态响应的幅值随时间分别增强、增强、减弱.这表明在不同的流速和转速组合下,系统不同阶模态响应会出现不同的增强或减弱现象.
图9 不同流速和转速下阻尼因子对前3 阶模态响应的影响Fig.9 The effects of damping factor on the first three mode responses for different flow rates and rotating speeds
本文通过能量法建立了含两通道的旋转输流管动力学控制方程,通过与文献结果对比验证了本文方法的正确性,研究了流体单位长度质量、不同流速比和输流管长细比对系统临界流速的影响,并分别分析了相位变化与阻尼因子对前3 阶模态响应的影响,得出了以下结论.
(1)在单位时间流量相同时,双通道模型的临界流速值大于单通道模型的,且随着转速的增大,双通道模型的临界流速值增长得更快;
(2)随着两通道内流速比值减小,第1 阶模态和第3 阶模态在发生失稳时的临界流速值增大,且第2 阶模态始终未发生失稳;
(3)旋转运动引入的陀螺效应对系统第1 阶特征根轨迹影响较小,但会显著影响第2,3 阶特征根轨迹,使其发生绕圈现象,并多次穿越虚轴,从而导致第2,3 阶模态振动存在多个临界流速.
(4)输流管横向位移模态响应出现相位差,呈现出行波特性,即出现"节点"移动现象,在不同的流速和转速组合下,系统不同阶模态响应会出现不同的增强或减弱现象.
本文首次研究了含两通道的旋转输流管动力学特性,为含复杂冷却通道布局的旋转叶片的设计提供一定理论参考.但本文中只考虑两通道为简单的对称式布局,输流管内部复杂不对称的通道布局等情况有待进一步的研究与讨论.
附录
猜你喜欢流管输液管流速直肠癌行可降解支架法肠道完全转流术患者的护理护理与康复(2022年10期)2022-10-22基础激励作用下悬臂输流管的振动实验研究*动力学与控制学报(2022年3期)2022-08-25液体压强与流速的关系初中生学习指导·提升版(2022年4期)2022-05-11『流体压强与流速的关系』知识巩固中学生数理化·八年级物理人教版(2022年4期)2022-04-26一次性防脱落输液管针头的设计临床护理杂志(2021年6期)2022-01-05山雨欲来风满楼之流体压强与流速大众科学(2020年7期)2020-10-26颈动脉转流管在复杂颈动脉体瘤手术切除中的应用中国普通外科杂志(2018年12期)2019-01-09爱虚张声势的水小天使·六年级语数英综合(2018年1期)2018-10-08输液管插入法治疗牛非泡沫性瘤胃胀气畜牧兽医科技信息(2016年12期)2016-02-23能自动过滤气泡的医用输液管发明与创新·中学生(2015年11期)2015-11-30
