(李双平)AutoPFA解决管路系统恒定流瞬变流的特征线方法应用.doc

AutoPFA解决管路系统恒定流、瞬变流的特征线方法应用 摘要：本文以特征线方法为基础，研究优易管网流体AutoPFA如何使用扬开林提出的一个简单有效的方法将恒定流和瞬变流的求解方法统一起来。这一方法尤其适用于工程上对具有复杂边界条件的系统仿真。为了适用给排水、石油、化工、AutoPFA将不同行业使用管道损失系数模型统一转化成达西损失系数，使AutoPFA能够使用一个解决方案适用各种管路系统模型。 引言 管路系统是工业和民用设备及船舶、航天工程的重要组成部分，管路系统的安全运行是整个系统安全性的重要环节。水锤冲击在管路系统经常发生，通过对管路系统各种工况水锤产生的压力波的传播和衰减的分析，确保管路系统在各种运行条件安全。 选择计算方法 目前求解非恒定流的数值方法可以分为三大类：有限差分法、有限元法和有限分析法。有限差分法是数值计算方法中最古老的方法，直至今日，在求解流体运动的数值分析方法中仍占主导地位。基于有限差分法的特征线方法，首先将偏微分方程转化成常微分方程然后使用显示差分法解常微分方程。水力过渡的特征线离散方程为： (1.1) (1.2) 其中： 式中：和为计算节点上的瞬态压力和瞬态质量流量，为前一时间段计算的常量数。由于每个边界条件和管段内部节点都是在一个时间段中独立的分析，所以这种方法尤其适用计算机仿真具有复杂边界条件的管路系统。 在一般情况下，管路水力瞬变的求解是由恒定流动开始的，所以知道恒定流动条件下每一根管道各个计算断面的压力和流量是瞬变计算的基础。在以往的分析计算中，通常采用不同的方法求解恒定流和瞬变流，这使得计算机程序过分复杂。采用特征线方法计算管路系统的恒定流动。这样做的好处是将恒定流和水力瞬变流的求解方法统一起来，使计算机程序简化。 求解恒定流我们关心的是瞬变末了恒定流的解，不是瞬变的过程，所以我们可以采用一些特殊方法以提高特征线方法求解恒定流的计算效率。为了减少计算量，可取管系每一根管道的分段数为1。我们可以取任何一个假想的状态作为特征线方法计算恒定流的初始条件。根据计算精度要求计算出各节点的恒定压力和流量。此外，由于在恒定流条件下调压室、空气罐和空气阀对管道系统的流动无影响，因此在计算中可以不考虑它们，即在计算中可以认为系统中没有它们，以加快计算的收敛速度。 数学模型 管路系统水力瞬变数学模型的建立主要是边界条件的处理。边界条件都是由相邻管段的正，负特征线方程和边界点的水利特性方程（连续方程，能量方程）联立求解。下面介绍优易管网流体AutoPFA管段和各种节点边界条件的数学模型： 管道摩擦模型 由公式其中为达西损失系数，K为管道附加损失K系数。特征线使用达西损失系数计算管道沿程了损失，为了适用不同行业标准需要优易管网流体AutoPFA提供8中计算管道摩擦损失模型，包括低速给水使用的海曾-威廉公式，原油输送MIT Equation及化工行业使用的MIT Equation等计算管道损失的模型。AutoPFA将这些模型计算损失使用的损失系数转化成达西损失系数，使用一个数学模型满足不同行业的需要。 下面讨论海曾-威廉系数转化达西损失系数。 其中： 为海曾-威廉系数，达西损失系数，流速V的单位为feet/sec直径D的单位为feet 管道波速计算 AutoPFA考虑了管道结构对波速的影响，提供7种波速修正模型计算管道波速，以下是管道波速计算公式： 其中 波速 K流体体积模量 P流体密度 E管道杨式模量 D管道内直径 e管道壁厚 u泊松比 c1修正系数 薄壁管只有上游端固定 薄壁管两端固定 薄壁管膨胀接头连接 厚壁管只有上游端固定 厚壁管两端固定 厚壁管膨胀接头连接 圆隧道 已知压力边界条件 水库、指定滞止压力、指定静止压力为已知压力边界条件，根据节点特性得到相连管道端点的滞止压力。管道端点的流量未知，我们可以根据特征线相容方程计算管道端点的流量。 如果已知管道上游端压力 如果已知管道下游端压力 如果已知压力随时间改变，可以通过当前系统时间步长计算得到当前压力再计算当前流量。 已知流量边界条件 已知流量边界条件知道管道入口或出口的流量。管道端点的压力未知，我们可以根据特征线相容方程计算管道端点的压力。 如果已知管道上游端流量 如果已知管道下游端流量 如果已知流量随时间改变，可以通过当前系统时间步长计算得到当前流量再计算当前流量。 分支元件边界条件 分支元件可以连接多根管道，流入分支元件的流量等于流出分支元件的流量。各连接管道的端点压力在分支元件处相等。 根据特征线相容方程得到管道流入分支元件的流量方程 同理得管道流出分支元件的流量方程 将所有管道的流量相加得： 如果存在一个已知的流量流入或流出分支元件 现在可以得到分支元件的压力