考虑流固耦合作用旋转波纹管固有频率计算的有限元法 - 农业工程学报.pdfVIP

第 15 卷　第 2 期 农 业 工 程 学 报 V o l. 15　N o. 2 　　　 　　　　　　　 　　　　　　　 50 1999 年 6 月 T ransactions of the CSA E Jun. 1999 考虑流固耦合作用旋转波纹管固有频率 计算的有限元法 ① 张进国 　史　岩　刁新华　蒋金琳　王树英 (莱阳农学院) 摘　要　该文提出用有限元法计算考虑流固耦合作用的旋转波纹管的固有频率。将波纹管用两节 点旋转锥壳单元离散化, 单元之间以节点环相连接, 将载荷和位移沿环向展开为傅立叶谐波级数, 降 低了计算维数, 减少了计算自由度和计算机时。推导出考虑波纹管与管内流体动力耦合作用的旋转 波纹管动力方程, 计算出旋转波纹管的固有频率, 可为工程中波纹管的动力设计提供理论依据。 关键词　流固耦合　波纹管　有限元　固有频率 在石油、化工、冶金、仪表等工业领域中, 旋转波纹管得到了广泛的应用。文献[1、2、3 ]分别 用摄动法、差分法和有限元法对旋转波纹管在轴对称载荷作用下的静态非线性进行了分析, 文 献[4 ]给出了任意载荷作用下旋转波纹管的几何非线性有限元静态分析方法。为提高结构系统 [ 5 ] 的动力性能和抗震性能, 计算波纹管的固有频率是必要的 , 文献[ 6 ]用有限元法计算了空旋 转波纹管的固有频率值。本文考虑波纹管与管内流体动力耦合作用, 提出考虑流固动力耦合作 用旋转波纹管固有频率计算的有限单元法。 1　基本方程 旋转波纹管是旋转壳体, 在动力特征值分析中, 仍采用两节点旋转锥壳单元[ 4～ 7 ] 将其离散 化, 单元与单元之间以节点环相连接。 两节点旋转锥壳单元如图 1 所示。单元节点编号为 、 i , 单元径向长度为 。法向坐标 、径向坐标 、环向坐标 j L N s 切线方向组成右手螺旋直角坐标系; 径向坐标 、轴向坐标 r 、环向坐标 组成右手螺旋圆柱坐标系; 坐标 与坐标 夹 z s z [ 4～ 7 ] 角为 。将载荷和位移沿环向展开为傅立叶级数 , 、、 u v 和 分别为节点的径向位移、环向位移、法向位移和径向 w 转角傅立叶级数系数。 则对于谐波阶次 , 单元 节点位移列阵为 n e e T ( ) 图 1　两节点旋转锥壳单元 { } = [u 　v 　w 　 　u 　v 　w 　 ] 1 n in in in in j n j n j n j n 由哈密尔顿变分原理可以得到, 时刻 阶谐波单元 的动 F ig. 1　T he conical frustum t n e