一种管路系统的流固耦合动力学特性预测方法
【技术领域】
[0001] 本发明涉及一种管路系统的流固耦合动力学特性预测方法，更具体的，本发明涉 及一种用于含有较大场传递管路元件的大型复杂管路系统的流固耦合动力学问题预测。
【背景技术】
[0002] 大型复杂管路系统广泛应用于船舶与海洋工程、石油化工、能源与电力工业等诸 多领域中，由于管路系统中存在栗、阀门、分支管路等内部激励以及基础或铺板振动等外部 激励，因而会诱发管路结构及管内流体的振动，这些管路结构及流体的振动还会相互耦合 而产生流固耦合（FSI)振动，因此这些看似安静的管路系统在输送流体的同时，还将振动 能量沿管壁和管内流体传播到管路系统的各个位置，造成管路及与之相连的其他元件和精 密仪器的破坏，影响管路系统安全和管路动力系统的正常运行，严重时会造成巨大的经济 损失，如何实现管路系统流固耦合动力学特性的快速高精度预测一直是国内外研究的热点 和难点。
[0003] 目前管路系统的流固耦合动力学特性的预测方法较多，常见的方法有四种：（1) 特征线方法（MOC)，MOC是一种时域数值求解方法，适用于简单管路的时域响应计算，特别 是流体压力波的瞬态响应分析，但MOC需要在时间和空间中离散求解，较为复杂，且很难考 虑管道的弹性支撑条件，多段管路计算时还存在插值误差和不同特征线的相交等问题；（2) 有限元法（FEM)，FEM是通过节点、单元构建管道及流体的运动形式，计算效率取决于模型 尺寸的大小或分析频率的高；（3)特征线-有限元法方法（MOC-FEM) ,MOC-FEM是通过MOC求 解管内流体运动，通过FEM求解管道结构运动，然后根据接触面的平衡关系实现FSI问题的 求解，MOC-FEM为输流管路动力学问题的求解提供了一个新的途径，但该方法的求解效率较 MOC低得多；(4)传递矩阵方法（TMM)，TMM是通过波动方程直接描述管路系统的运动。与 其他方法相比，TMM求解过程简单，易于编程且计算效率高，因而在管路系统结构振动预测， 管内流体压力脉动的传播特性分析及流固耦合动力学预测等问题中具有较多应用，但是目 前此方法尚难以应用于具有较长管路元件的大型管路系统的流固耦合计算中，原因之一是 TMM在计算长管路时会出现计算数值溢出和预测结果不稳定的现象。
[0004] 针对TMM的计算不稳定、积累数值误差大及数值溢出等"病态问题"，芮筱亭等人 在《多体系统传递矩阵法及其应用》一书中集中介绍和讨论了传递矩阵的无量纲化，特征值 问题的改进算法以及Riccati传递矩阵方法等几种常见的处理方法（芮筱亭，员来峰，陆毓 琪，何斌，王国平.多体系统传递矩阵法及其应用[M].北京：科学出版社，2008)。胡培民 通过一个一维数组和两个变量来存储一个数据的方法解决了超精度数据的存储问题（胡 培民.传递矩阵法在高频振动分析中的应用[J].振动与冲击，1996, 15(4):50-52.)。Liu 和 Li (Liu Gongmin, Li Yanhua. Vibration analysis of liquid-filled pipelines with elastic constraints[J]. Journal of Sound and Vibration, 2011, 330(13):3166-3181.) 根据管路系统的流固耦合动力学特性，推导得到了 12个一元四阶常微分方程和2个一元二 阶常微分方程，进而求解得到了管路系统的频域解析解，该方法在一定程度上能够减少由 于矩阵变化带来的计算误差，上述这些的改进措施，在一定程度上均提高了传递矩阵在大 型复杂管路系统流固耦合动力学特性预测精度和预测结果的稳定性，但利用这些方法在计 算管路系统的流固耦合响应时并无优势，并且难以实现模块化编程预测。
[0005] 国内外目前利用传递矩阵方法进行管路系统动力学预测的流程是：首先分析管路 系统的构成形式，然后将各管路元件的数理模型带入管路系统中，并根据各管路元件连接 处平衡及传递关系，将其组合为整体传递矩阵形式，进而利用定解条件，预测管路系统的动 力学特性，但该方法在预测管路系统的流固耦合振动响应时会出现计算数值溢出和预测结 果不稳定的现象。
【发明内容】
[0006] 针对现有技术中存在的缺陷，本发明的目的是提供一种精确、稳定、尚效的管路系 统的流固耦合动力学特性预测方法，其包括以下步骤：
[0007] 步骤一、根据管路系统中每个管路元件的属性，将管路系统离散为直管元件、弯管 元件和管路附件元件三种类型；
[0008] 步骤二、根据离散后的管路元件的类型，建立各管路元件相应的场传递流固耦合 矩阵或点传递流固耦合矩阵；
[0009] 步骤三、将离散后的管路系统重新组合为若干传递矩阵单元体，且根据所述管路 系统预测的管壁弯曲波最小波长或者预测的上限频率计算传递矩阵单元体的特征长度范 围；
[0010] 步骤四、结合管路系统边界处流体和结构已知条件组成边界条件矩阵和外部激励 列向量，利用传递矩阵方法，预测管路系统任意位置的流固耦合动力学特性。
[0011] 在上述技术方案的基础上，离散管路系统时，管路系统包括场传递元件和点传递 元件，所述场传递元件包括直管元件和弯管元件，所述点传递元件为管路附件元件。
[0012] 在上述技术方案的基础上，所述管路附件元件包括管路支撑、分支接头和阀门。
[0013] 在上述技术方案的基础上，根据所述管路系统预测的管壁弯曲波最小波长计算传 递矩阵单元体的最大特征长度的计算方法为：
[0015] 其中，Lniax表示重新组合后的单元体的最大特征长度，C表示比例系数，
:表示与管壁杨氏模量（E)，管路内径（R1),管路外径（R。)，管内流体密 度（Pf)，管壁密度（Pp)和泊松比（u)相关的计算式，λ_为预测的管壁的弯曲波最小波 长。
[0016] 在上述技术方案的基础上，根据所述管路系统预测的上限频率计算传递矩阵单元 体的最大特征长度的计算方法为：
[0018] 其中，表示组合后的单元体的最大特征长度，C表示比例系数，
表示一个与管壁杨氏模量（E)，管路内径（R1),管路外径（R。)，管内流体 密度（Pf)和管壁密度（Pp)相关的计算式，s_表示可以预测的上限频率。
[0019] 在上述技术方案的基础上，对于任意的管路系统，其任意考察位置的状态向量的 预测方法为：
[0021] 其中，表示在管路系统中选取的任意考察位置的状态向量，U aA表示从任意考 察位置所在的单元体的始端到考察位置的传递矩阵，Utq表示提取矩阵，其作用是将考察位 置所在单元体的始端的状态矩阵从管路系统的整体传递矩阵中提取出来，D表示由流体已 知条件和结构已知条件组成的边界条件矩阵，U tt表示由各单元体共同组成的整体传递矩 阵，[]1表示矩阵的逆阵，Fm1表示管路系统的外部激励列向量。
[0022] 特别地，当管路系统由单个单元体组成时，计算管路系统任意考察位置的状态向 量的方法为：
[0024] 其中，N = 2,3,4,5,6,7……，式中，表示在管路系统中选取的任意考察位置 的状态向量，Onxn表示NXN维的零矩阵，I ΝΧΝ表示NXN维的单位矩阵，D ΝΧ3Ν表示NX 3Ν维 的边界条件矩阵，UInxn表示从管路始端到考察位置的NXN维的传递矩阵，U2 ΝΧΝ表示从考 察位置到管路末端的ΝΧ3Ν维传递矩阵，[]1表示矩阵的逆阵，F 3nxn表示3ΝΧΝ维的激励 列向量，N表示用于描述管路系统流固耦合振动的变量数，其值为大于等于2的自然数。
[0025] 与现有技术相比，本发明以重组优化后的传递矩阵单元体为计算对象，可快速地 建立管路系统的流固耦合动力学问题的预测模型，通过划分单元体分段计算的技术，能够 消除管路系统流固耦合动力学问题数值溢出和预测结果不稳定问题，预测过程简便，非常 有利于编程计算，并且预测精度高。
[0026] 本发明所涉及的预测技术通过管路系统预测的管壁弯曲波最小波长或者预测的 上限频率计算传递矩阵单元体的特征长度范围，其在保证了预测精度的同时，降低计算量， 使得本发明的管路系统的流固耦合动力学特性预测方法更为简单、可以实现程式化建模、 具有预测精度高和速度快等显著特点。
【附图说明】
[0027] 图1为大型复杂管路系统流固耦合动力学预测步骤框图；
[0028] 图2为管路模型示意图；
[0029] 图3为传统传递矩阵方法计算的管路系统横向流固耦合振动响应；
[0030] 图4为本发明预测的管路系统横向流固耦合振动响应；
[0031] 图5为复杂管路系统示意图；
[0032] 图6为管路系统A点X方向振动加速度响应曲线；
[0033] 图7为管路系统B点X方向振动加速度响应曲线；
[0034] 图8为管路系统C点x方向振动加速度响应曲线。
【具体实施方式】
[0035] 以下结合附图及实施例对发明作进一步详细说明。
[0036] 参见图1所示，本发明管路系统的流固耦合动力学特性预测方法具体包括以下步 骤：
[0037] 步骤一、根据管路系统中每个管路元件的属性，将管路系统离散为直管元件、弯管 元件和管路附件元件三种类型；离散管路系统时，管路系统包括场传递元件和点传递元件， 所述场传递元件包括直管元件和弯管元件，所述点传递元件为管路附件元件。
[0038] 步骤二、根据离散后的管路元件的类型，建立各管路元件相应的场传递流固耦合 矩阵或点传递流固耦合矩阵；其具体为，对场传递元件建立场传递流固耦合矩阵，对点传递 元件建立点传递流固親合矩阵。
[0039] 步骤三、将离散后的管路系统重新组合为若干传递矩阵单元体，且根据所述管路 系统预测的管壁弯曲波最小波长或者预测的上限频率计算传递矩阵单元体的特征长度范 围；
[0040] 当已知预测的管壁弯曲波最小波长时，则采用最小波长方法计算传递矩阵单元体 的特征长度范围；根据所述管路系统预测的管壁弯曲波最小波长计算传递矩阵单元体的最 大特征长度的计算方法为：
[0042] 其中，Lniax表示重新组合后的单元体的最大特征长度，C表示比例系数，
表示与管壁杨氏模量（E)，管路内径（R1),管路外径（R。)，管 内流体密度（Pf)，管壁密度（Pp)和泊松比（u)相关的计算