[发明专利]一种基于CFD的可调导叶对泵能量性能预测的方法

权利要求书

1.一种基于CFD的可调导叶对泵能量性能预测的方法，其特征是，基于泵的CFD数值计算构建泵的能量性能数据库；采用泵综合特性绘制Visual Fortran程序对泵能量性能数据进行处理绘制泵的综合特性曲线，并输出CFD数值计算的可调导叶安放角时任意数量的泵性能数据；以泵的效率和泵的扬程为因变量，以可调导叶的安放角和泵的流量为自变量采用多元非线性回归分析方法分别构建通过2个自变量预测因变量的方法；具体包括以下步骤：

1)基于CFD数值计算构建不同可调导叶安放角时泵的能量性能数据库；

采用CFD数值计算程序，对配不同可调导叶安放角的泵进行8～10个工况点的能量性能的三维定常数值计算，基于CFD数值计算结果，分别求出泵的扬程、效率；

泵的扬程H的计算式：

H=[Σi=1Noutlet(Psoutletρg)iNoutlet+Σi=1Noutlet(voutlet22g)iNoutlet]out-[Σi=1Nin(Psinρg)iNin+Σi=1Nin(vin22g)iNin]in+Zout-Zin---(1)]]>

式(1)中：P_s为断面各节点的静压值；Z为断面几何中心的位能；v为断面各网格节点的绝对速度；N为断面网格节点数总和；下表outlet表示出口断面；下表in表示进口断面；

泵的效率η的计算式：

η=30ρgQHZZπnTp;TP=Σi=1NΔAi{[r→×(τ→·n→)]·ez}-Σi=1NPΔAi[(r→×n→)·e→z]---(2)]]>

式(2)中：为转轴方向的单位向量；△A_i为压力面或吸力面上第i单元的面积；为△A_i上的单位向量；为向径；为不含静压力P的应力张量；Q为泵装置流量；n为叶轮转速。

2)基于Visual Fortran程序绘制不同可调导叶时泵的综合特性曲线，并输出CFD数值计算的可调导叶安放角时任意数量的泵能量性能数据；

为绘制不同可调导叶安放角时泵的能量性能综合特性曲线，规定：可调导叶出口水流方向与转轮旋转方向相同时，可调导叶安放角为正，反之为负；采用Visual Fortran程序对步骤1建立的泵能量性能数据库进行自动处理绘制出泵的能量性能综合特性曲线，即泵的扬程、流量和效率3者的复杂关系曲线，并采用内插值的方法给出CFD数值计算的不同可调导叶安放角时任意数量的泵能量性能数据，为最终提出可调导叶对泵能量性能预测的数学模型奠定足够的泵能量性能样本数据；

3)构建可调导叶对泵能量性能预测的数学模型；

采用多元非线性回归分析方法(Multiple Non-liner Regression)构造可调导叶对泵能量性能预测的数学模型，该数学模型一种快捷方法；

采用Gauss-Newton算法对非线性回归方程系数进行求解；将可调导叶安放角θ记为x₁，流量Q记为x₂，扬程H记为y₁，可得非线性回归方程1；将进口导叶片安放角θ记为x₁，流量H记为x₂，泵效率η记为y₂，可得非线性回归方程2，回归预测模型的准确度可通过数学模型的判定系数R₂进行评判(判定系数R₂需大于0.98)，若判定系数R₂不满足要求，则需回到步骤2，重新输出更多的泵性能数据样本；

建立的2元3次非线性回归预测数学模型的通用表达式如下：

y=a+bx1+cx2+dx12+ex22+fx1x2+gx13+hx23+ix1x22+jx12x2]]>

式中：a₁、b₁、c₁、d₁、e₁、f₁、g₁、h₁、i₁、j₁为回归方程的常数及系数。

2.根据权利要求1所述的一种基于CFD的可调导叶对泵能量性能预测的方法，其特征是，所述Visual Fortran程序包括如下流程：

步骤1、输入泵的流量、扬程和效率的组数；

步骤2、分别依次输入流量、扬程、效率坐标的间隔值；

步骤3、分别依次输入流量起点和终点坐标、扬程的起点和终点坐标、效率的起点和终点坐标(提示：插值最好为公约数的倍数)；

步骤4、输入划分区间数(提示：区间数越大曲线拟合越接近)；

步骤5、采用中点分割和拐点分割相混合的Bezier曲线降阶方法分别绘制流量-扬程曲线，流量-效率曲线；

步骤6、在拟合的流量-扬程线上采用插值法寻找效率值相等的点，并采用Bezier曲线降阶方法对效率值相等的点进行连接拟合，最终获得泵的综合特性曲线。