基于ANSYS二次开发的输液曲管振动特性分析.pdfVIP

2005年 10月 沈阳航空工业学院学报 O ct. 2005 第 2 2卷第 5期 Journal of Shenyang In stitute of A eronau tical Engineering Vo l. 22　No. 5 文章编号 : 1007 1385 2005 05 002 1 03 基于 AN SYS二次开发的输液曲管振动特性分析 于秀坤 　朱 　虹 　金基铎 　王 　维 沈阳航空工业学院航空宇航工程学院 , 辽宁 沈阳　110034 摘 　要 :针对悬臂的边界条件 ,建立了在流体作用下输液曲管的振动方程及有限元分析模型 ,并 研究了曲管的振动特性 。通过对 AN SYS进行 U PF 二次开发 , 自定义曲管单元进行特征值求解 , 得出了流体流速对曲管特征频率影响的规律 , 同时考察了输液曲管的稳定性 ,求出了输液曲管发 生颤振失稳时的临界流速 。计算实例表明自定义单元具有很好的精度 。 关键词 :输液曲管 ; U PF; 振动特性 ; 临界流速 中图分类号 : TB 123; O24 1. 82 文献标识码 : A 输液管在航空工程 、水利工程 、石油工业等领 2　输液曲管的有限元模型 域应用广泛 ,在一定工作条件下 ,这些管道会发生 强烈振动至不稳定状态的现象 。过去对输液管振 2. 1　单元分析 动问题所做的研究 ,大都是针对输液直管进行的 第 i 个单元的平衡方程为 : 研究 ,对输液曲管的振动及稳定性研究相对较少 。 [M ] [ q¨] | [D ] q | [D ] q 0 2 i i i i i i 输液曲管与输液直管相比,描述运动状态需要 式中, M , D , K 分别是第 i号单元的质量矩阵、陀 i i i 的变量多,而且静平衡位置还与流速有关。运用有 螺矩阵和刚度矩阵 。图 1为单元节点编号 。 限元法[ 1, 2 ] ,研究中心线不可伸缩的平面输液曲管在 特定边界条件和不同流速下的振动特性 ,通过计算 结果 ,分析得出了边界条件和液体流速对曲管振动 特性的影响规律 ,并且通过特性值的求解 ,进一步求 出曲管发生颤振失稳时的临界速度 。 1　输液曲管的运动方程 根据 A. K. M isra[ 3, 4 ] 等人的研究成果 ,在不考 虑内外阻尼的影响下 ,对于中心线不可伸缩的输 图 1　曲管模型中某一单元 2. 2　整体分析 液曲管的平面内振动方程 无量纲形式 为 : 6 4 2 4 2 整体平衡方程为 : η η 4 η η 2 η 4 +2θ + θ +u +2θ + θη + 6 4 2 4 2 [M ] q¨ + [D ] q + [ K ] q 0 3 ξ ξ ξ ξ ξ i i 3 η η 2 2 η 2 式中M 、D 、K 分别是输液曲管的总质量矩阵、总 1/ 2 2 2 η β θ θ 2 u 3 + + 2 2 - 2 0 τ ξ τ ξ ξ β τ τ 陀螺矩阵和总刚度矩阵 。D 、K 元素与质量比 、 流速 u有关 。 1 1/ 2 ω s B I η ξ τ ; ; t 4 ; 3　算例 L L M +M L t f M f L 采用 AN SYS二次开发 , 自定义曲管单元 ,对 β θ - ; M +M R t f 半圆形悬臂输液曲管进行分析计算 。 分别代表无量纲的切 向位移, 曲线坐标, 时 图 2 是悬臂曲管模型 ,为验证计算结果 , 曲管 间, 流速, 质量比和中心角 。 各参数同文献 〔1 〕: 曲率半径 R 200cm , 曲管和 流体单位长度的质量 M t M f 0. 1kg/ cm ,抗弯 收稿日期 : 2005 07 0 1 8 2 刚度 E I 2 ×10 kgfcm ,为了便于与文献 〔1 〕比 作者简介 :于秀坤 1963 , 女, 辽宁沈阳人, 副教授 22 沈阳航空工业学院学报 第 22 卷 2 曲管各阶模态运动的 A rgand 图 图 3 图 2　曲管模型 较 ,取如下频率作为一阶、二阶、三阶频率 ,见表 1 π 注 :表 1中的频率是除以 2 之后的结果 。 1 曲管在不同流速下的频率 表 1 表 1　曲管不同流速下的频率 流速 一阶频率 二阶频率 三阶频率 无量纲 0. 0 0. 000000 + 1. 60833i 0. 000000 + 2. 88054i 0. 000000 + 6. 38902i 0. 5 - 0. 0333 + 1. 58488i - 0. 03479 +2.