埋地管道关阀水锤数值模拟及结构响应流固耦合分析.docxVIP

? ? 埋地管道关阀水锤数值模拟及结构响应流固耦合分析 ? ? 虞学军 郑舟斌 熊伟东 金建波 张文奇 （1.舟山市特种设备检测研究院；2.杭州市特种设备检测研究院；3.浙江大学能源工程学院） 管道广泛应用于城市供水/暖、化工及核电等行业中，对于长距离管道,为了布设需求，常常以埋地的方式进行布管。 但有时管道部件运行状态的突然改变（如下游阀门的突然关闭）会引起管道内流体流动状态的瞬间变化， 产生流体冲击力，并在介质弹性作用下持续较长时间。 这种管道流体动态过渡的过程， 会导致系统部件振动，严重时会引起管道的破坏和断裂，这种现象称为水锤［1～4］。 水锤引起的埋地管道损坏，使埋地管道的检测频率和成本大幅提高。 因此，研究水锤造成的冲击对管道的可靠性分析和事故预测至关重要。 水锤导致埋地管道振动甚至损坏，是一种典型的流体流动现象引起的结构振动。 需要在研究水锤的水动力学特性的基础上，再进行管道结构力学的研究。 对于水锤的水动力学特性的试验研究，Ruus E和Karney B B在考虑管道常数和摩擦参数的基础上， 探讨止回阀关闭后引起的水锤波特性，确定了水锤影响下泵不同位置压头的升高和降低的规律［5］。在理论研究的方法方面，早期主要基于特征线法进行。 如Tian W X等通过特征线法研究并联泵交替启动过程中的阀门水锤现象，并提出采用阻尼矩阵降低止回阀的关闭速度的方法，以减轻阀门水锤的潜在危害［6］。 由于特征线法主要针对管道系统研究，结果只可得到关键节点的动态数据，无法从内流场角度进行分析，也就无法进一步将内流场的压力等数据作为载荷进一步研究结构力学的流固耦合。Al-Khomairi A M通过实验证明基于特征线法计算瞬态管道流动存在误差［7］。 基于管道水锤的水动力学特性的研究，Keramat A等率先进行管道水锤的流固耦合特性研究，但采用的方法是特征线法-有限元法，只能计算管道的一维变形［8］。 随着计算流体动力学（CFD）的发展和计算机技术的进步。 CFD方法相比传统的一维方法，其结果具有更高的准确性且更加详细。 Wu D Z等采用Fluent软件对管道系统的阀门快速开启过程进行模拟， 模拟结果与PIV实验结果有良好的吻合度，得到过渡过程中内流畅的瞬态演化过程［9］。 将CFD应用在水锤的研究中，尤其是将CFD得到的压力等流场数据作为载荷，并通过计算结构动力学（CSD）方法研究管道结构在水锤影响下的响应特性，可以得到传统一维方法难以计算的管道水锤流固耦合特性。 综上所述， 探索一种将CFD和CSD耦合的计算方法，以准确预测埋地管道的关阀水锤水动力学特性，并得到管道结构的响应特性，在得到管道水锤的流固耦合机理的同时，还能对水锤引起管道破坏的位置进行预测，减少埋地管道检测的成本。 因此，以Bergant经典管道关阀水锤试验［2，3］为基础，通过多种CFD方法探究最优数值解法，再采用流固耦合的方法， 借由CSD软件分析埋地管道在水锤压力波载荷下的响应特性，揭示管道结构的振动机理。 1 数值计算模型 通过流固耦合方法分析管道关阀水锤的水动力学特性和结构响应特性， 需要分别建立CFD计算模型和CSD计算模型。 1.1 CFD计算模型 澳大利亚研究委员会（ARC）曾经针对管道的水锤进行了试验研究。 该水锤试验主要由Bergant A完成。 试验台由一根长37.2 m，内径22 mm，壁厚1.6 mm的管道作为试验段（图1），两端连接有球阀和水箱，球阀的下游设有段管道和水箱，两水箱间连有压力调节器，并由计算机控制水箱的压力。 此外，试验台还包括安装在管道试验段的压力传感器和数据采集仪。 Bergant A完成了多组试验，通过关阀使管道试验段中产生水锤波，并测量水锤波的特性。 其中，管道入口压力为0.32 MPa， 流场初始速度为0.1 m/s被武汉理工大学的杨成选取，并通过ANSYS CFX计算软件进行数值模拟，并分析了多种初始流速和关阀速度对水锤特性的影响，但其计算值与试验值吻合性较差［10］，而且没有得到压力波的衰减规律。 为此，笔者采用相同的工况，在上述研究的基础上继续进行数值方法的研究，以探索水锤数值模拟的最优方法。 图1 水锤试验台示意图 1.1.1 网格的收敛性分析 为了实现对发射系统的CFD模拟， 需要对管道试验段的求解区域进行网格剖分。 相比于非结构化网格， 结构化网格划分方法可在减少网格数的同时，大幅提高网格的质量，尤其可以有针对性地控制特定区域的节点分布， 提高计算的收敛性和光滑性。此外，为了避免网格划分对计算结果产生影响， 采用Roache提出的网格收敛系数法（GCI方法） 对网格收敛性进行检验。 GCI方法的指标——GCI值是模拟值偏离渐进数值的度量。选择两种网格划分方案，即粗糙网格（网格数为96 065）和