CN112560361B - 一种基于仿真结果修正的离心泵扬程值预测方法 - Google Patents

一种基于仿真结果修正的离心泵扬程值预测方法 Download PDF

Info

Publication number
CN112560361B
CN112560361B CN202011464956.9A CN202011464956A CN112560361B CN 112560361 B CN112560361 B CN 112560361B CN 202011464956 A CN202011464956 A CN 202011464956A CN 112560361 B CN112560361 B CN 112560361B
Authority
CN
China
Prior art keywords
pump
value
new
flow
centrifugal pump
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011464956.9A
Other languages
English (en)
Other versions
CN112560361A (zh
Inventor
张翮辉
常春平
游浩林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dragon Totem Technology Hefei Co ltd
Hubei Tianmen Yongqiang Pump Industry Co ltd
Zhongou Hubei Intellectual Property Service Co ltd
Original Assignee
Xiangtan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xiangtan University filed Critical Xiangtan University
Priority to CN202011464956.9A priority Critical patent/CN112560361B/zh
Publication of CN112560361A publication Critical patent/CN112560361A/zh
Application granted granted Critical
Publication of CN112560361B publication Critical patent/CN112560361B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Control Of Non-Positive-Displacement Pumps (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明提供了一种基于仿真结果修正的离心泵扬程值预测方法，包括离心泵原泵数据矩阵的获得，原泵流量‑效率计算机仿真值一元二次函数关系式的拟合，原泵第一、第二扬程修正矩阵的获得，新泵数据矩阵的获得，新泵流量～效率计算机仿真值一元二次函数关系式的拟合，新泵第一、第二数据矩阵的获得，新泵扬程预测矩阵的获得和新泵流量～扬程预测值函数关系式的拟合等步骤。本发明从真实实验测试数据出发，只依靠一次离心泵实物制作，基于有限的测试数据获得离心泵扬程计算机仿真结果的修正方法，而不需进行成本较高的效率实验测试，具有客观科学、逻辑严密、测试简便等优点，有助于提高离心泵扬程仿真结果的可靠性进而节约开发时间与资金成本。

Description

一种基于仿真结果修正的离心泵扬程值预测方法
技术领域
本发明涉及离心泵领域，特别涉及一种基于仿真结果修正的离心泵扬程值预测方法。
背景技术
离心泵是一种应用最为广泛的泵类设备，主要起到液体输送和增压的作用。离心泵大量用于国民经济的各个部门，在工农业生产和人们生活中具有不可替代的重要价值。由于离心泵的各种应用场合千差万别，为了满足应用需要并尽可能地提高运行能效水平，设计单位时常需要对离心泵进行定制化设计：即以原有的离心泵设计方案为基础进行修改优化，并获得新的离心泵设计方案。
在离心泵的定制化设计过程中，为了加快工作进度，减少试错次数，降低实验和测试成本，设计单位越来越广泛地借助计算机仿真的手段来进行离心泵性能的预测，主要目的是得到离心泵设计方案下不同流量点所对应的效率值和扬程值。其中，准确获得离心泵不同流量点所对应的扬程值是评价离心泵设计方案优劣的必备条件——因为保证某个流量运行条件下的扬程是离心泵正常工作的前提，如果实际扬程与计算机仿真得到的扬程值相差较大，就意味着离心泵偏离所需工况，造成增压不足或增压过度，带来输送困难或管道破坏等危险。然而，由于诸多复杂因素的影响，离心泵扬程的计算机仿真结果通常会与真实的实验测试结果存在或多或少的差异，如何对离心泵扬程的计算机仿真结果进行修正，使之更为接近实验测试中的真实情况，是离心泵优化设计领域十分重要的课题。
目前公知的各种技术方案中，虽然广泛应用计算机仿真手段进行离心泵不同流量工况下扬程值的预测，也有相关文献报道对离心泵扬程的计算机仿真值和实验测试值进行比较和分析，但对于如何根据离心泵扬程的实验测试值对计算机仿真结果进行精确修正，则缺乏***深入的研究。故当前亟待发展离心泵扬程仿真结果的修正方法，以更为精确地预测不同流量工况下离心泵的扬程值，从而提高离心泵计算机仿真技术的有效性和精确性。
因此，在现有的离心泵计算机仿真技术的基础上，针对离心泵计算机仿真获得的扬程值，如何通过有限的真实实验测试数据，并充分利用离心泵计算机过程中可以获得的其他信息，建立通用而***的修正方法，基于计算机仿真和修真方法实现各种离心泵设计方案下更为精确的离心泵扬程值预测，无疑是一个亟待解决的技术问题。
发明内容
为了解决上述技术问题，本发明提供一种客观科学、操作方便的基于仿真结果修正的离心泵扬程值预测方法。
本发明解决上述问题的技术方案是：一种基于仿真结果修正的离心泵扬程值预测方法，包括以下步骤：
步骤1、获得离心泵原泵不同流量点对应的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式。
根据离心泵原泵设计方案，确定考察的流量范围，分别对原泵方案进行实验测试和计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure BDA0002833755150000021
式(1)中d0为原泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵原泵流量值，n为不少于7的正整数，其中Q1和Qn分别为考察流量范围内的最小流量值和最大流量值；第二行至第四行分别是与同列流量值相对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值。
步骤2、拟合得到离心泵原泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵原泵效率最高值及其对应的设计流量。
针对原泵的数据矩阵d0中离心泵原泵不同流量点对应的效率计算机仿真值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵原泵流量～效率计算机仿真值函数关系式：
E(Q)＝aQ2+bQ+c (2)
式(2)中Q为离心泵原泵流量值，E(Q)为流量为Q时的离心泵原泵效率计算机仿真值，a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
针对式(2)，计算出最高效率点对应的流量值即离心泵原泵流量设计Qd，并计算出离心泵原泵最高效率值Ed：
Figure BDA0002833755150000031
Figure BDA0002833755150000032
式(3)中Qd为离心泵原泵的设计流量值，式(4)中Ed为离心泵原泵的最高效率值，式(3)和式(4)中的a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
步骤3、以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，并在此基础上获得相对应的原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII。
在原泵数据矩阵d0第一行的流量值中确定离心泵原泵临界流量值Qm，并以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000033
Figure BDA0002833755150000034
分别以原泵第一数据矩阵dI和原泵第二数据矩阵dII为基础，计算得到原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII：
Figure BDA0002833755150000041
Figure BDA0002833755150000042
式(7)和式(8)中的REi和ΔHi分别为原泵效率相对值和原泵扬程修正值，i＝1，2，…，m，…，n，原泵效率相对值REi和原泵扬程修正值ΔHi的计算式分别为：
Figure BDA0002833755150000043
ΔHi＝HTi-HSi (10)
式(7)中Ed为步骤2得到的原泵最高效率值，Ei为原泵效率计算机仿真值；式(8)中的HTi和HSi分别为原泵扬程实验测试值和原泵扬程计算机仿真值。
步骤4、确定离心泵新泵设计方案，获得离心泵新泵不同流量点对应的扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式。
在离心泵原泵设计方案基础上修改得到离心泵新泵设计方案，新泵设计方案考察的流量范围与原泵保持一致，对离心泵新泵设计方案进行计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure BDA0002833755150000044
式(11)中d′0为新泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵新泵流量值，q为不少于7的正整数，其中Q′1和Q′q分别为考察流量范围内的最小流量值和最大流量值；第二行至第三行分别是与同列流量值相对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值。
步骤5、拟合得到离心泵新泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵新泵效率最高值及其对应的设计流量。
针对离心泵新泵数据矩阵d′0中不同流量点对应的效率计算机仿真值，基于最小二乘法进行一元二次函数拟合，并获得离心泵新泵流量～效率计算机仿真值函数关系式：
E′(Q′)＝a′Q′2+b′Q′+c′ (12)
式(12)中Q′为离心泵新泵流量值，E′(Q′)为流量为Q′时的离心泵新泵效率仿真计算值，a′、b′、和c′分别为拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
针对式(12)，计算出最高效率点对应的流量值即离心泵新泵设计流量Q′d，并计算出离心泵新泵最高效率值E′d：
Figure BDA0002833755150000051
Figure BDA0002833755150000052
式(13)中Q′d为离心泵新泵的设计流量值，式(14)中E′d为离心泵新泵的最高效率值，式(13)和式(14)中的a′、b′、和c′分别为本步骤拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
步骤6、以离心泵新泵临界流量值Q′p为分隔界限将新泵数据矩阵d′0拆分，并补充效率相对值数据形成新泵第一数据矩阵d′Ι和新泵第二数据矩阵d′ΙΙ。
在新泵数据矩阵d′0第一行的流量值中确定临界流量值Q′p，以临界流量值Q′p为分隔界限将新泵数据矩阵d′0拆分，并对拆分后的两个数据矩阵补充效率相对值数据后形成新泵第一数据矩阵d′Ι和新泵第二数据矩阵d′ΙΙ，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000061
Figure BDA0002833755150000062
式(15)和式(16)中最后一行的RE′i为新泵效率相对值，i＝1，2，…，p，…，q，RE′i是在新泵数据矩阵d′0拆分后的基础上进一步补充的数据，新泵效率相对值RE′i的计算式分别为：
Figure BDA0002833755150000063
式(17)中E′d为步骤5得到的新泵最高效率值，E′i为新泵效率计算机仿真值。
步骤7、计算得到新泵第一扬程修正矩阵和新泵第二扬程修正矩阵，并最终获得新泵扬程预测矩阵，拟合得到新泵流量～扬程预测值函数关系式。
根据步骤3获得的原泵第一扬程修正矩阵hI和步骤6获得的新泵第一数据矩阵d′Ι，计算得到新泵第一扬程修正矩阵h′Ι：
Figure BDA0002833755150000064
式(18)中，新泵第一扬程修正矩阵h′Ι第二行的ΔH′i为新泵扬程修正值，i＝1，2，…，p。
根据步骤3获得的原泵第二扬程修正矩阵hII和步骤6获得的新泵第二扬程预测矩阵d′ΙΙ，计算得到新泵第二扬程修正矩阵h′ΙΙ：
Figure BDA0002833755150000071
式(19)中，新泵第二扬程修正矩阵h′ΙΙ第二行的ΔH′i为新泵扬程修正值，i＝p+1，p+2，…，q。
由新泵第一扬程修正矩阵h′Ι和步骤6获得的新泵第一数据矩阵d′Ι计算得到新泵第一扬程预测矩阵yh′Ι：
Figure BDA0002833755150000072
式(20)中最后一行的YH′i为新泵扬程预测值，i＝1，2，…，p，由同流量下对应的新泵第一扬程修正矩阵h′Ι中的新泵扬程修正值ΔH′i和新泵第一数据矩阵d′Ι中的新泵扬程计算机仿真值H′Si相加得到：
YH′i＝ΔH′i+H′Si (21)
由新泵第二扬程修正矩阵h′ΙΙ和步骤6获得的新泵第二数据矩阵d′ΙΙ计算得到新泵第二扬程预测矩阵yh′ΙΙ：
Figure BDA0002833755150000073
式(22)中最后一行的YH′i为新泵扬程预测值，i＝p+1，p+2，…，q，由同流量下对应的新泵第二扬程修正矩阵h′ΙΙ中的新泵扬程修正值ΔH′i和新泵第二数据矩阵d′ΙΙ中的新泵扬程计算机仿真值H′Si相加得到。
将新泵第一扬程预测矩阵yh′Ι和新泵第二扬程预测矩阵yh′ΙΙ横向合并为新泵扬程预测矩阵yh′：
Figure BDA0002833755150000081
式(23)中第一行为新泵的流量值，第二行是与第一行的流量值相对应的新泵的扬程预测值。
针对式(23)新泵扬程预测矩阵yh′中离心泵新泵不同流量点对应的扬程预测值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵新泵流量～扬程预测值函数关系式：
YH′(Q′)＝a″Q′2+b″Q′+c″ (24)
式(24)中Q′为离心泵新泵流量值，H′(Q′)为流量为Q′时的离心泵新泵扬程预测值，a″、b″、和c″分别为拟合得到的离心泵新泵流量～扬程预测值一元二次函数关系式的二次项系数、一次项系数和常数项。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤1或步骤4中离心泵原泵或新泵设计方案，包括但不限于离心泵的几何形状和尺寸、额定转速、叶轮旋转方向、运行流量范围等信息。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤1中离心泵原泵扬程实验测试值的获取包括如下步骤：
步骤S1：根据离心泵原泵设计方案，加工制作离心泵实物，并将离心泵实物安装至离心泵水力性能测试台；
步骤S2：在离心泵水力性能测试台上，调节离心泵出口管道阀门设置离心泵运行流量为其考察的最小流量值Q1，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S3：增大离心泵出口管道阀门开度，适当增大离心泵运行流量值，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S4：重复步骤S3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程实验测试值。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤1中离心泵原泵的扬程计算机仿真值和效率计算机仿真值的获取包括如下步骤：
步骤O1：根据离心泵原泵设计方案，绘制离心泵水体部分三维几何图，并划分数值仿真用网格；
步骤O2：将离心泵数值仿真用网格导入计算流体力学软件，设置通过离心泵的流量为其考察的最小流量值Q1，添加离心泵额定转速和叶轮旋转方向等必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O3：调整计算流体力学软件设置，适当增大通过离心泵的流量值，添加离心泵额定转速和叶轮旋转方向等必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O4：重复步骤O3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程值和效率值。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤3中离心泵原泵临界流量值Qm是原泵数据矩阵d0第一行中满足如下条件的某一个具有唯一性的流量值：Qm≤Qd且Qm+1>Qd，其中Qd为步骤3确定的离心泵原泵的设计流量值，Qm+1是原泵数据矩阵d0第一行中Qm右侧相邻的流量值。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤6中离心泵新泵临界流量值Q′p是新泵数据矩阵d′0第一行中满足如下条件的某一个具有唯一性的流量值：Q′p≤Q′d且Q′p+1＞Q′d，其中Q′d为步骤5确定的离心泵新泵的设计流量值，Q′p+1是新泵数据矩阵d′0第一行中Q′p右侧相邻的流量值。
上述基于仿真结果修正的离心泵扬程值预测方法，所述步骤7中新泵第一扬程修正矩阵h′Ι(或新泵第二扬程修正矩阵h′ΙΙ)中的新泵扬程修正值ΔH′i，i＝1，2，…，p(或i＝p+1，p+2，…，q)，由新泵第一数据矩阵d′Ι(或新泵第二数据矩阵d′ΙΙ)与原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中的相关数据经过插值运算得到：
Figure BDA0002833755150000101
式(25)中ΔH′i为新泵第一扬程修正矩阵h′Ι(或新泵第二扬程修正矩阵h′ΙΙ)中与流量Q′i对应的新泵扬程修正值，i＝1，2，…，p(或i＝p+1，p+2，…，q)，RE′i为新泵第一数据矩阵d′Ι(或新泵第二数据矩阵d′ΙΙ)中流量Q′i对应的新泵效率相对值，REf和REg分别为原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中与RE′i相差最小和相差第二小的两个原泵效率相对值，ΔHf和ΔHg分别为原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中与REf和REg相对应的原泵扬程修正值。
本发明的有益效果在于：
1、本发明针对计算机仿真与实验测试这两种途径得到的离心泵扬程值存在偏差这一事实，通过各种矩阵的构造和计算进行修正，建立起离心泵仿真结果的修正方法，以便更为精确地预测各种新的设计方案下离心泵扬程值随流量的变化情况，这种方法客观科学，通用性强。
2、本发明充分利用了计算机仿真过程得到的离心泵效率数据，基于最小二乘法分别对离心泵原泵和新泵的计算机仿真结果进行拟合获得各自的流量～效率计算机仿真值一元二次函数式，在此基础上获得原泵和新泵的设计流量值，并对离心泵原泵和新泵的相关数据矩阵进行拆分以便分别处理，最终预测获得新泵不同流量下的扬程修正值，整个过程步骤清晰，逻辑严密，操作方便。
3、本发明提供的基于仿真结果修正的离心泵扬程值预测方法，只需要在起始阶段进行一次离心泵实物样机的制作和有限次的扬程测试，且避免了成本较高而精度通常难以保证的离心泵效率的实验测试，故消耗成本很低，数据处理量较小，数据字段清晰且不易出错，整个计算过程易于编程实现，可以根据扬程修正矩阵等信息进行多个不同新的设计方案的离心泵的扬程计算机仿真值的预测修正。
附图说明
图1为本发明的流程图。
图2为其中本发明实施例中离心泵新泵不同流量下扬程值的比较，曲线1为离心泵新泵不同流量点对应的扬程计算机仿真值的平滑连线，曲线2为离心泵新泵不同流量点对应的扬程预测值的平滑连线，曲线3为离心泵新泵不同流量点对应的扬程实验测试值的平滑连线。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示，一种基于仿真结果修正的离心泵扬程值预测方法，包括以下步骤：
步骤1、获得离心泵原泵不同流量点对应的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式。
根据离心泵原泵设计方案，确定考察的流量范围，分别对原泵方案进行实验测试和计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure BDA0002833755150000111
式(1)中d0为原泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵原泵流量值，n为不少于7的正整数，其中Q1和Qn分别为考察流量范围内的最小流量值和最大流量值；第二行至第四行分别是与同列流量值相对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值。
其中离心泵原泵设计方案，包括但不限于离心泵的几何形状和尺寸、额定转速、叶轮旋转方向、运行流量范围等信息。
所述离心泵原泵扬程实验测试值的获取包括如下步骤：
步骤S1：根据离心泵原泵设计方案，加工制作离心泵实物，并将离心泵实物安装至离心泵水力性能测试台；
步骤S2：在离心泵水力性能测试台上，调节离心泵出口管道阀门设置离心泵运行流量为其考察的最小流量值Q1，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S3：增大离心泵出口管道阀门开度，适当增大离心泵运行流量值，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S4：重复步骤S3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程实验测试值。
所述离心泵原泵的扬程计算机仿真值和效率计算机仿真值的获取包括如下步骤：
步骤O1：根据离心泵原泵设计方案，绘制离心泵水体部分三维几何图，并划分数值仿真用网格；
步骤O2：将离心泵数值仿真用网格导入计算流体力学软件，设置通过离心泵的流量为其考察的最小流量值Q1，添加离心泵额定转速和叶轮旋转方向等必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O3：调整计算流体力学软件设置，适当增大通过离心泵的流量值，添加离心泵额定转速和叶轮旋转方向等必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O4：重复步骤O3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程值和效率值。
步骤2、拟合得到离心泵原泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵原泵效率最高值及其对应的设计流量。
针对原泵的数据矩阵d0中离心泵原泵不同流量点对应的效率计算机仿真值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵原泵流量～效率计算机仿真值函数关系式：
E(Q)＝aQ2+bQ+c (2)
式(2)中Q为离心泵原泵流量值，E(Q)为流量为Q时的离心泵原泵效率计算机仿真值，a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
针对式(2)，计算出最高效率点对应的流量值即离心泵原泵流量设计Qd，并计算出离心泵原泵最高效率值Ed：
Figure BDA0002833755150000131
Figure BDA0002833755150000132
式(3)中Qd为离心泵原泵的设计流量值，式(4)中Ed为离心泵原泵的最高效率值，式(3)和式(4)中的a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
步骤3、以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，并在此基础上获得相对应的原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII。
在原泵数据矩阵d0第一行的流量值中确定离心泵原泵临界流量值Qm，并以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000133
Figure BDA0002833755150000134
分别以原泵第一数据矩阵dI和原泵第二数据矩阵dII为基础，计算得到原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII：
Figure BDA0002833755150000141
Figure BDA0002833755150000142
式(7)和式(8)中的REi和ΔHi分别为原泵效率相对值和原泵扬程修正值，i＝1，2，…，m，…，n，原泵效率相对值REi和原泵扬程修正值ΔHi的计算式分别为：
Figure BDA0002833755150000143
ΔHi＝HTi-HSi (10)
式(7)中Ed为步骤2得到的原泵最高效率值，Ei为原泵效率计算机仿真值；式(8)中的HTi和HSi分别为原泵扬程实验测试值和原泵扬程计算机仿真值。
所述离心泵原泵临界流量值Qm是原泵数据矩阵d0第一行中满足如下条件的某一个具有唯一性的流量值：Qm≤Qd且Qm+1>Qd，其中Qd为步骤3确定的离心泵原泵的设计流量值，Qm+1是原泵数据矩阵d0第一行中Qm右侧相邻的流量值。
步骤4、确定离心泵新泵设计方案，获得离心泵新泵不同流量点对应的扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式。
在离心泵原泵设计方案基础上修改得到离心泵新泵设计方案，新泵设计方案考察的流量范围与原泵保持一致，对离心泵新泵设计方案进行计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure BDA0002833755150000144
式(11)中d′0为新泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵新泵流量值，q为不少于7的正整数，其中Q′1和Q′q分别为考察流量范围内的最小流量值和最大流量值；第二行至第三行分别是与同列流量值相对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值。
所述离心泵新泵设计方案，包括但不限于离心泵的几何形状和尺寸、额定转速、叶轮旋转方向、运行流量范围等信息。
步骤5、拟合得到离心泵新泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵新泵效率最高值及其对应的设计流量。
针对离心泵新泵数据矩阵d′0中不同流量点对应的效率计算机仿真值，基于最小二乘法进行一元二次函数拟合，并获得离心泵新泵流量～效率计算机仿真值函数关系式：
E′(Q′)＝a′Q′2+b′Q′+c′ (12)
式(12)中Q′为离心泵新泵流量值，E′(Q′)为流量为Q′时的离心泵新泵效率仿真计算值，a′、b′、和c′分别为拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
针对式(12)，计算出最高效率点对应的流量值即离心泵新泵设计流量Qd′，并计算出离心泵新泵最高效率值Ed′：
Figure BDA0002833755150000151
Figure BDA0002833755150000152
式(13)中Qd′为离心泵新泵的设计流量值，式(14)中Ed′为离心泵新泵的最高效率值，式(13)和式(14)中的a′、b′、和c′分别为本步骤拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项。
步骤6、以离心泵新泵临界流量值Qp′为分隔界限将新泵数据矩阵d0′拆分，并补充效率相对值数据形成新泵第一数据矩阵dΙ′和新泵第二数据矩阵dΙΙ′。
在新泵数据矩阵d0′第一行的流量值中确定临界流量值Qp′，以临界流量值Qp′为分隔界限将新泵数据矩阵d0′拆分，并对拆分后的两个数据矩阵补充效率相对值数据后形成新泵第一数据矩阵dΙ′和新泵第二数据矩阵dΙΙ′，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000161
Figure BDA0002833755150000162
式(15)和式(16)中最后一行的RE′i为新泵效率相对值，i＝1，2，…，p，…，q，RE′i是在新泵数据矩阵d0′拆分后的基础上进一步补充的数据，新泵效率相对值RE′i的计算式分别为：
Figure BDA0002833755150000163
式(17)中E′d为步骤5得到的新泵最高效率值，E′i为新泵效率计算机仿真值。
所述离心泵新泵临界流量值Q′p是新泵数据矩阵d0′第一行中满足如下条件的某一个具有唯一性的流量值：Q′p≤Q′d且Q′p+1＞Q′d，其中Qd′为步骤5确定的离心泵新泵的设计流量值，Q′p+1是新泵数据矩阵d0′第一行中Qd′右侧相邻的流量值。
步骤7、计算得到新泵第一扬程修正矩阵和新泵第二扬程修正矩阵，并最终获得新泵扬程预测矩阵，拟合得到新泵流量～扬程预测值函数关系式。
根据步骤3获得的原泵第一扬程修正矩阵hI和步骤6获得的新泵第一数据矩阵dΙ′，计算得到新泵第一扬程修正矩阵h′Ι：
Figure BDA0002833755150000171
式(18)中，新泵第一扬程修正矩阵h′Ι第二行的ΔH′i为新泵扬程修正值，i＝1，2，…，p。
根据步骤3获得的原泵第二扬程修正矩阵hII和步骤6获得的新泵第二扬程预测矩阵dΙΙ′，计算得到新泵第二扬程修正矩阵h′ΙΙ：
Figure BDA0002833755150000172
式(19)中，新泵第二扬程修正矩阵h′ΙΙ第二行的ΔHi′为新泵扬程修正值，i＝p+1，p+2，…，q。
由新泵第一扬程修正矩阵h′Ι和步骤6获得的新泵第一数据矩阵dΙ′计算得到新泵第一扬程预测矩阵yh′Ι：
Figure BDA0002833755150000173
式(20)中最后一行的YH′i为新泵扬程预测值，i＝1，2，…，p，由同流量下对应的新泵第一扬程修正矩阵hΙ′中的新泵扬程修正值ΔH′i和新泵第一数据矩阵d′Ι中的新泵扬程计算机仿真值H′Si相加得到：
YH′i＝ΔH′i+H′Si (21)
由新泵第二扬程修正矩阵h′ΙΙ和步骤6获得的新泵第二数据矩阵dΙΙ′计算得到新泵第二扬程预测矩阵yh′ΙΙ：
Figure BDA0002833755150000181
式(22)中最后一行的YH′i为新泵扬程预测值，i＝p+1，p+2，…，q，由同流量下对应的新泵第二扬程修正矩阵h′ΙΙ中的新泵扬程修正值ΔHi′和新泵第二数据矩阵dΙΙ′中的新泵扬程计算机仿真值H′Si相加得到。
将新泵第一扬程预测矩阵yh′Ι和新泵第二扬程预测矩阵yh′ΙΙ横向合并为新泵扬程预测矩阵yh′：
Figure BDA0002833755150000182
式(23)中第一行为新泵的流量值，第二行是与第一行的流量值相对应的新泵的扬程预测值。
针对式(23)新泵扬程预测矩阵yh′中离心泵新泵不同流量点对应的扬程预测值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵新泵流量～扬程预测值函数关系式：
YH′(Q′)＝a″Q′2+b″Q′+c″ (24)
式(24)中Q′为离心泵新泵流量值，H′(Q′)为流量为Q′时的离心泵新泵扬程预测值，a″、b″、和c″分别为拟合得到的离心泵新泵流量～扬程预测值一元二次函数关系式的二次项系数、一次项系数和常数项。
所述新泵第一扬程修正矩阵h′Ι(或新泵第二扬程修正矩阵h′ΙΙ)中的新泵扬程修正值ΔH′i，i＝1，2，…，p(或i＝p+1，p+2，…，q)，由新泵第一数据矩阵d′Ι(或新泵第二数据矩阵d′ΙΙ)与原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中的相关数据经过插值运算得到：
Figure BDA0002833755150000191
式(25)中ΔH′i为新泵第一扬程修正矩阵h′Ι(或新泵第二扬程修正矩阵h′ΙΙ)中与流量Qi′对应的新泵扬程修正值，i＝1，2，…，p(或i＝p+1，p+2，…，q)，RE′i为新泵第一数据矩阵d′Ι(或新泵第二数据矩阵d′ΙΙ)中流量Qi′对应的新泵效率相对值，REf和REg分别为原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中与RE′i相差最小和相差第二小的两个原泵效率相对值，ΔHf和ΔHg分别为原泵第一扬程修正矩阵hI(或原泵第二扬程修正矩阵hII)中与REf和REg相对应的原泵扬程修正值。
实施例
某离心泵，原泵设计方案为：进口直径0.15m，出口直径0.12m，额定转速1480r/min，叶轮旋转方向为逆时针方向，额定流量350m3/h，运行流量范围为200～500m3/h。亦即考察的最小流量值Qmin＝700m3/h、最大流量值Qmax＝1700m3/h。
根据该离心泵原泵设计方案，加工制作离心泵实物，并通过水力性能测试台测试，获得的不同流量工况下的扬程实验测试值见表1。
表1离心泵原泵设计方案不同流量工况下的扬程实验测试值
流量(m<sup>3</sup>/h) 200 250 300 350 400 450 500
扬程实验测试值(m) 55.2 54.7 53.9 52.3 49.2 45.8 38.6
根据离心泵原泵设计方案，绘制离心泵水体部分三维几何图，并划分数值仿真用网格，导入计算流体力学软件CFX设置边界条件如下：进口为某个流量值，出口静压力为0，叶轮区域旋转速度1480r/min，叶轮旋转方向为逆时针方向，采用有限体积法进行迭代计算，基于计算结果获得离心泵某个流量值下的进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值。
计算机仿真得到的原泵不同流量工况下的扬程计算机仿真值和效率计算机仿真值见表2。
表2离心泵原泵设计方案不同流量工况下的扬程计算机仿真值和效率计算机仿真值
流量(m3/h) 200 250 300 350 400 450 500
扬程计算机仿真值(m) 56.3 55.6 54.6 52.8 49.3 44.7 37.5
效率计算机仿真值 0.375 0.406 0.423 0.438 0.425 0.408 0.377
根据表1和表2中的数据，得到离心泵原泵数据矩阵d0为：
Figure BDA0002833755150000201
上述离心泵原泵数据矩阵d0中第一行为从小到大排序的离心泵原泵流量值，考察流量范围内的最小流量值为Q1＝200m3/h、最大流量值为Q7＝500m3/h；第二行至第四行分别是与同列流量值相对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值。
针对原泵的数据矩阵d0中离心泵原泵不同流量点对应的效率计算机仿真值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵原泵流量～效率计算机仿真值函数关系式为：
E(Q)＝aQ2+bQ+c＝-2.55E-6Q2+1.80E-3Q+1.17E-1
式中Q为离心泵原泵流量值，E(Q)为流量为Q时的离心泵原泵效率计算机仿真值，a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项，即a＝-2.55E-6，b＝1.80E-3，c＝1.17E-1。
由此计算出最高效率点对应的流量值即离心泵原泵流量设计Qd，并计算出离心泵原泵最高效率值Ed：
Figure BDA0002833755150000211
Figure BDA0002833755150000212
因为Qd＝353m3/h，故在原泵数据矩阵d0第一行的流量值中确定离心泵原泵临界流量值Qm为Qm＝Q4＝350m3/h，Q4右侧相邻的流量值Q5＝400m3/h，这时有Q4≤Qd且Q5>Qd。
以离心泵原泵临界流量值Qm＝Q4为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000213
Figure BDA0002833755150000214
分别以原泵第一数据矩阵dI和原泵第一数据矩阵dII为基础，计算得到原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII：
Figure BDA0002833755150000215
Figure BDA0002833755150000216
对于原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII，二者的第一行为原泵效率相对值，第二行为原泵扬程修正值。
在离心泵原泵设计方案基础上修改得到离心泵新泵设计方案，新泵设计方案考察的流量范围与原泵保持一致，对离心泵新泵设计方案进行计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值，得到新泵的数据矩阵d′0为：
Figure BDA0002833755150000221
新泵的数据矩阵d′0中第一行为从小到大排序的离心泵新泵流量值，其中Q′1＝200m3/h和Q′7＝500m3/h分别为考察流量范围内的最小流量值和最大流量值；第二行至第三行分别是与同列流量值相对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值。
针对离心泵新泵数据矩阵d0′中不同流量点对应的效率计算机仿真值，基于最小二乘法进行一元二次函数拟合，并获得离心泵新泵流量～效率计算机仿真值函数关系式：
E′(Q′)＝a′Q′2+b′Q′+c′＝-2.13E-6Q′2+1.65E-3Q′+1.17E-1
拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数为a′＝-2.13E-6，一次项系数为b′＝1.65E-3，常数项为c′＝1.17E-1。
计算出最高效率点对应的流量值即离心泵新泵设计流量Qd′，并计算出离心泵新泵最高效率值Ed′：
Figure BDA0002833755150000222
Figure BDA0002833755150000223
因为Qd′＝387m3/h，故在新泵数据矩阵d0′第一行的流量值中确定离心泵新泵临界流量值Qp′为Qp′＝Q4′＝350m3/h，Q4′右侧相邻的流量值Q5′＝400m3/h，这时有Q4′≤Qd′且Q5′>Qd′。以临界流量值Qp′＝Q4′＝350m3/h为分隔界限将新泵数据矩阵d′0拆分，并对拆分后的两个数据矩阵补充效率相对值数据后形成新泵第一数据矩阵dΙ′和新泵第二数据矩阵dΙΙ′，而数据矩阵中原有的数据排序和数值均保持不变：
Figure BDA0002833755150000231
Figure BDA0002833755150000232
新泵第一数据矩阵dΙ′和新泵第二数据矩阵dΙΙ′中最后一行为新泵效率相对值。
根据原泵第一扬程修正矩阵hI和新泵第一数据矩阵dΙ′，计算得到新泵第一扬程修正矩阵h′Ι：
Figure BDA0002833755150000233
以新泵流量Q′1＝200m3/h为例，介绍该流量点下其新泵扬程修正值的计算过程。
新泵流量Q′1＝200m3/h时，由新泵第一数据矩阵dΙ′查得其对应的新泵效率相对值RE1′＝0.838，又由原泵第一扬程修正矩阵hI查得与RE1′＝0.838相差最小和相差第二小的两个原泵效率相对值分别为REf＝RE1＝0.862、REg＝RE2＝0.933，二者对应的原泵扬程修正值分别为ΔHf＝ΔH1＝-1.1m、ΔHg＝ΔH2＝-0.9m，由此计算出该流量点下新泵的扬程修正值为：
Figure BDA0002833755150000241
类似地，根据原泵第二扬程修正矩阵hII和新泵第二扬程预测矩阵d′ΙΙ，计算得到新泵第二扬程修正矩阵h′ΙΙ：
Figure BDA0002833755150000242
由新泵第一扬程修正矩阵h′Ι和新泵第一数据矩阵d′Ι计算得到新泵第一扬程预测矩阵yh′Ι：
Figure BDA0002833755150000243
该矩阵第一行中的数值为流量值，最后一行中的数值为新泵扬程预测值，由同流量下对应的新泵第一扬程修正矩阵h′Ι中的新泵扬程修正值ΔH′i和新泵第一数据矩阵d′Ι中的新泵扬程计算机仿真值H′Si相加得到。
类似地，由新泵第二扬程修正矩阵h′ΙΙ和步骤6获得的新泵第二数据矩阵d′ΙΙ计算得到新泵第二扬程预测矩阵yh′ΙΙ：
Figure BDA0002833755150000244
将新泵第一扬程预测矩阵yh′Ι和新泵第二扬程预测矩阵yh′ΙΙ横向合并为新泵扬程预测矩阵yh′：
Figure BDA0002833755150000251
该新泵扬程预测矩阵yh′中第一行为新泵的流量值，第二行是与第一行的流量值相对应的新泵的扬程预测值。
由新泵扬程预测矩阵yh′中离心泵新泵不同流量点对应的扬程预测值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵新泵流量～扬程预测值函数关系式为：
YH′(Q′)＝a″Q′2+b″Q′+c″＝-2.07E-4Q′2+8.82E-2Q′+4.63E1
为了验证本发明提供的离心泵扬程值预测方法的可靠性，根据该离心泵新泵设计方案，加工制作离心泵新泵实物，并通过水力性能测试台测试，获得的不同流量工况下的扬程实验测试值见表3。
表3离心泵新泵设计方案不同流量工况下的扬程实验测试值
流量(m3/h) 200 250 300 350 400 450 500
扬程实验测试值(m) 56.3 55.6 54.6 52.8 49.3 44.7 37.5
为了更方便观察以上离心泵新泵不同流量下的扬程计算机仿真值、扬程预测值和扬程实验测试值之间的差别情况，在考察的流量范围200～500m3/h内，对离心泵新泵相关数据点予以展示，如图2所示。
图2为本发明实施例中离心泵新泵不同流量下扬程值的比较，曲线1为离心泵新泵不同流量点对应的扬程计算机仿真值的平滑连线，曲线2为离心泵新泵不同流量点对应的扬程预测值的平滑连线，曲线3为离心泵新泵不同流量点对应的扬程实验测试值的平滑连线。
图2中容易发现，曲线2比曲线1更为接近曲线3即离心泵新泵不同流量点对应的扬程实验测试值的平滑连线。这充分说明使用本发明提供的离心泵扬程值预测方法，在曲线1即离心泵新泵不同流量点对应的扬程计算机仿真值的平滑连线的基础上，经过相关修正预测得到曲线2即离心泵新泵不同流量点对应的扬程预测值的平滑连线，与真实的实验测试结果即曲线3明显更为接近。
1、本发明针对计算机仿真与实验测试这两种途径得到的离心泵扬程值存在偏差这一事实，通过各种矩阵的构造和计算进行修正，建立起离心泵仿真结果的修正方法，以便更为精确地预测各种新的设计方案下离心泵扬程值随流量的变化情况，这种方法客观科学，通用性强。
2、本发明充分利用了计算机仿真过程得到的离心泵效率数据，基于最小二乘法分别对离心泵原泵和新泵的计算机仿真结果进行拟合获得各自的流量～效率计算机仿真值一元二次函数式，在此基础上获得原泵和新泵的设计流量值，并对离心泵原泵和新泵的相关数据矩阵进行拆分以便分别处理，最终预测获得新泵不同流量下的扬程修正值，整个过程步骤清晰，逻辑严密，操作方便。
3、本发明提供的基于仿真结果修正的离心泵扬程值预测方法，只需要在起始阶段进行一次离心泵实物样机的制作和有限次的扬程测试，且避免了成本较高而精度通常难以保证的离心泵效率的实验测试，故消耗成本很低，数据处理量较小，数据字段清晰且不易出错，整个计算过程易于编程实现，可以根据扬程修正矩阵等信息进行多个不同新的设计方案的离心泵的扬程计算机仿真值的预测修正。
本实施例提供的基于仿真结果修正的离心泵扬程值预测方法，针对计算机仿真与实验测试这两种途径得到的离心泵扬程值存在偏差这一事实，只需要在起始阶段进行一次离心泵实物样机的制作和有限次的扬程测试，避免成本较高而精度通常难以保证的离心泵效率的实验测试，同时充分利用计算机仿真过程得到的离心泵效率数据，通过各种矩阵的构造和计算进行修正，并对离心泵原泵和新泵的相关数据矩阵进行拆分以便分别处理，最终建立起离心泵仿真结果的修正方法，以便更为精确地预测各种新的设计方案下离心泵扬程值随流量的变化情况，这种方法客观科学，通用性强，成本低廉，数学处理过程步骤清晰，逻辑严密，操作方便，便于编程实现。

Claims (7)

1.一种基于仿真结果修正的离心泵扬程值预测方法，包括以下步骤：
步骤1、获得离心泵原泵不同流量点对应的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式；
根据离心泵原泵设计方案，确定考察的流量范围，分别对原泵方案进行实验测试和计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure FDA0003556094750000011
式(1)中d0为原泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵原泵流量值，n为不少于7的正整数，其中Q1和Qn分别为考察流量范围内的最小流量值和最大流量值；第二行至第四行分别是与同列流量值相对应的离心泵原泵的扬程实验测试值、扬程计算机仿真值和效率计算机仿真值；
步骤2、拟合得到离心泵原泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵原泵效率最高值及其对应的设计流量；
针对原泵的数据矩阵d0中离心泵原泵不同流量点对应的效率计算机仿真值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵原泵流量～效率计算机仿真值函数关系式：
E(Q)＝aQ2+bQ+c (2)
式(2)中Q为离心泵原泵流量值，E(Q)为流量为Q时的离心泵原泵效率计算机仿真值，a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
针对式(2)，计算出最高效率点对应的流量值即离心泵原泵流量设计Qd，并计算出离心泵原泵最高效率值Ed：
Figure FDA0003556094750000012
Figure FDA0003556094750000013
式(3)中Qd为离心泵原泵的设计流量值，式(4)中Ed为离心泵原泵的最高效率值，式(3)和式(4)中的a、b、和c分别为拟合得到的离心泵原泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
步骤3、以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，并在此基础上获得相对应的原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII；
在原泵数据矩阵d0第一行的流量值中确定离心泵原泵临界流量值Qm，并以离心泵原泵临界流量值Qm为分隔界限将原泵数据矩阵d0拆分为原泵第一数据矩阵dI和原泵第二数据矩阵dII，而数据矩阵中原有的数据排序和数值均保持不变：
Figure FDA0003556094750000021
Figure FDA0003556094750000022
分别以原泵第一数据矩阵dI和原泵第二数据矩阵dII为基础，计算得到原泵第一扬程修正矩阵hI和原泵第二扬程修正矩阵hII：
Figure FDA0003556094750000023
Figure FDA0003556094750000024
式(7)和式(8)中的REi和ΔHi分别为原泵效率相对值和原泵扬程修正值，i＝1，2，…，m，…，n，原泵效率相对值REi和原泵扬程修正值ΔHi的计算式分别为：
Figure FDA0003556094750000025
ΔHi＝HTi-HSi (10)
式(7)中Ed为步骤2得到的原泵最高效率值，Ei为原泵效率计算机仿真值；式(8)中的HTi和HSi分别为原泵扬程实验测试值和原泵扬程计算机仿真值；
步骤4、确定离心泵新泵设计方案，获得离心泵新泵不同流量点对应的扬程计算机仿真值和效率计算机仿真值，并写成数据矩阵形式；
在离心泵原泵设计方案基础上修改得到离心泵新泵设计方案，新泵设计方案考察的流量范围与原泵保持一致，对离心泵新泵设计方案进行计算机仿真，获得考察流量范围内不少于7个不同流量点对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值，并将以上获得的所有值写成矩阵形式：
Figure FDA0003556094750000031
式(11)中d′0为新泵的数据矩阵，该矩阵中第一行为从小到大排序的离心泵新泵流量值，q为不少于7的正整数，其中Q′1和Q′q分别为考察流量范围内的最小流量值和最大流量值；第二行至第三行分别是与同列流量值相对应的离心泵新泵的扬程计算机仿真值和效率计算机仿真值；
步骤5、拟合得到离心泵新泵流量～效率计算机仿真值一元二次函数关系式，并计算得到离心泵新泵效率最高值及其对应的设计流量；
针对离心泵新泵数据矩阵d′0中不同流量点对应的效率计算机仿真值，基于最小二乘法进行一元二次函数拟合，并获得离心泵新泵流量～效率计算机仿真值函数关系式：
E′(Q′)＝a′Q′2+b′Q′+c′ (12)
式(12)中Q′为离心泵新泵流量值，E′(Q′)为流量为Q′时的离心泵新泵效率仿真计算值，a′、b′、和c′分别为拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
针对式(12)，计算出最高效率点对应的流量值即离心泵新泵设计流量Q′d，并计算出离心泵新泵最高效率值E′d：
Figure FDA0003556094750000041
Figure FDA0003556094750000042
式(13)中Q′d为离心泵新泵的设计流量值，式(14)中E′d为离心泵新泵的最高效率值，式(13)和式(14)中的a′、b′、和c′分别为本步骤拟合得到的离心泵新泵流量～效率计算机仿真值一元二次函数关系式的二次项系数、一次项系数和常数项；
步骤6、以离心泵新泵临界流量值Q′p为分隔界限将新泵数据矩阵d′0拆分，并补充效率相对值数据形成新泵第一数据矩阵d′Ι和新泵第二数据矩阵d′ΙΙ；
在新泵数据矩阵d′0第一行的流量值中确定临界流量值Q′p，以临界流量值Q′p为分隔界限将新泵数据矩阵d′0拆分，并对拆分后的两个数据矩阵补充效率相对值数据后形成新泵第一数据矩阵d′Ι和新泵第二数据矩阵d′ΙΙ，而数据矩阵中原有的数据排序和数值均保持不变：
Figure FDA0003556094750000043
Figure FDA0003556094750000044
式(15)和式(16)中最后一行的RE′i为新泵效率相对值，i＝1，2，…，p，…，q，RE′i是在新泵数据矩阵d′0拆分后的基础上进一步补充的数据，新泵效率相对值RE′i的计算式分别为：
Figure FDA0003556094750000051
式(17)中E′d为步骤5得到的新泵最高效率值，E′i为新泵效率计算机仿真值；
步骤7、计算得到新泵第一扬程修正矩阵和新泵第二扬程修正矩阵，并最终获得新泵扬程预测矩阵，拟合得到新泵流量～扬程预测值函数关系式；
根据步骤3获得的原泵第一扬程修正矩阵hI和步骤6获得的新泵第一数据矩阵d′Ι，计算得到新泵第一扬程修正矩阵h′Ι：
Figure FDA0003556094750000052
式(18)中，新泵第一扬程修正矩阵h′Ι第二行的ΔH′i为新泵扬程修正值，i＝1，2，…，p；
根据步骤3获得的原泵第二扬程修正矩阵hII和步骤6获得的新泵第二扬程预测矩阵d′ΙΙ，计算得到新泵第二扬程修正矩阵h′ΙΙ：
Figure FDA0003556094750000053
式(19)中，新泵第二扬程修正矩阵h′ΙΙ第二行的ΔH′i为新泵扬程修正值，i＝p+1，p+2，…，q；
由新泵第一扬程修正矩阵h′Ι和步骤6获得的新泵第一数据矩阵d′Ι计算得到新泵第一扬程预测矩阵yh′Ι：
Figure FDA0003556094750000054
式(20)中最后一行的YH′i为新泵扬程预测值，i＝1，2，…，p，由同流量下对应的新泵第一扬程修正矩阵h′Ι中的新泵扬程修正值ΔH′i和新泵第一数据矩阵d′Ι中的新泵扬程计算机仿真值H′Si相加得到：
YH′i＝ΔH′i+H′Si (21)
由新泵第二扬程修正矩阵h′ΙΙ和步骤6获得的新泵第二数据矩阵d′ΙΙ计算得到新泵第二扬程预测矩阵yh′ΙΙ：
Figure FDA0003556094750000061
式(22)中最后一行的YH′i为新泵扬程预测值，i＝p+1，p+2，…，q，由同流量下对应的新泵第二扬程修正矩阵h′ΙΙ中的新泵扬程修正值ΔH′i和新泵第二数据矩阵d′ΙΙ中的新泵扬程计算机仿真值H′Si相加得到；
将新泵第一扬程预测矩阵yh′Ι和新泵第二扬程预测矩阵yh′ΙΙ横向合并为新泵扬程预测矩阵yh′：
Figure FDA0003556094750000062
式(23)中第一行为新泵的流量值，第二行是与第一行的流量值相对应的新泵的扬程预测值；
针对式(23)新泵扬程预测矩阵yh′中离心泵新泵不同流量点对应的扬程预测值，基于最小二乘法分别进行一元二次函数拟合，并获得离心泵新泵流量～扬程预测值函数关系式：
YH′(Q′)＝a″Q′2+b″Q′+c″ (24)
式(24)中Q′为离心泵新泵流量值，H′(Q′)为流量为Q′时的离心泵新泵扬程预测值，a″、b″、和c″分别为拟合得到的离心泵新泵流量～扬程预测值一元二次函数关系式的二次项系数、一次项系数和常数项。
2.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤1或步骤4中离心泵原泵或新泵设计方案，包括离心泵的几何形状和尺寸、额定转速、叶轮旋转方向、运行流量范围这些信息。
3.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤1中离心泵原泵扬程实验测试值的获取包括如下步骤：
步骤S1：根据离心泵原泵设计方案，加工制作离心泵实物，并将离心泵实物安装至离心泵水力性能测试台；
步骤S2：在离心泵水力性能测试台上，调节离心泵出口管道阀门设置离心泵运行流量为其考察的最小流量值Q1，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S3：增大离心泵出口管道阀门开度，适当增大离心泵运行流量值，记录该流量值下离心泵的进口压力和出口压力，并由进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程实验测试值；
步骤S4：重复步骤S3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程实验测试值。
4.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤1中离心泵原泵的扬程计算机仿真值和效率计算机仿真值的获取包括如下步骤：
步骤O1：根据离心泵原泵设计方案，绘制离心泵水体部分三维几何图，并划分数值仿真用网格；
步骤O2：将离心泵数值仿真用网格导入计算流体力学软件，设置通过离心泵的流量为其考察的最小流量值Q1，添加离心泵额定转速和叶轮旋转方向这两方面必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O3：调整计算流体力学软件设置，适当增大通过离心泵的流量值，添加离心泵额定转速和叶轮旋转方向这两方面必要信息，采用有限体积法进行迭代计算，基于计算结果读出离心泵进口压力、出口压力和主轴扭矩，由离心泵进出口压差再加上离心泵进出口面的中心高度差得到离心泵的扬程计算机仿真值，由主轴扭矩、额定转速、流量和扬程计算机仿真值计算该运行流量点下离心泵的效率值；
步骤O4：重复步骤O3至少5次，直至获得考察的最大流量值Qn所对应的离心泵的扬程值和效率值。
5.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤3中离心泵原泵临界流量值Qm是原泵数据矩阵d0第一行中满足如下条件的某一个具有唯一性的流量值：Qm≤Qd且Qm+1>Qd，其中Qd为步骤3确定的离心泵原泵的设计流量值，Qm+1是原泵数据矩阵d0第一行中Qm右侧相邻的流量值。
6.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤6中离心泵新泵临界流量值Q′p是新泵数据矩阵d′0第一行中满足如下条件的某一个具有唯一性的流量值：Q′p≤Q′d且Q′p+1＞Q′d，其中Q′d为步骤5确定的离心泵新泵的设计流量值，Q′p+1是新泵数据矩阵d′0第一行中Q′p右侧相邻的流量值。
7.权利要求1所述的基于仿真结果修正的离心泵扬程值预测方法，其特征在于，所述步骤7中新泵第一扬程修正矩阵h′Ι或新泵第二扬程修正矩阵h′ΙΙ中的新泵扬程修正值ΔH′i，i＝1，2，…，p或i＝p+1，p+2，…，q，由新泵第一数据矩阵d′Ι或新泵第二数据矩阵d′ΙΙ与原泵第一扬程修正矩阵hI或原泵第二扬程修正矩阵hII中的相关数据经过插值运算得到：
Figure FDA0003556094750000081
式(25)中ΔH′i为新泵第一扬程修正矩阵h′Ι或新泵第二扬程修正矩阵h′ΙΙ中与流量Q′i对应的新泵扬程修正值，i＝1，2，…，p或i＝p+1，p+2，…，q，RE′i为新泵第一数据矩阵d′Ι或新泵第二数据矩阵d′ΙΙ中流量Q′i对应的新泵效率相对值，REf和REg分别为原泵第一扬程修正矩阵hI或原泵第二扬程修正矩阵hII中与RE′i相差最小和相差第二小的两个原泵效率相对值，ΔHf和ΔHg分别为原泵第一扬程修正矩阵hI或原泵第二扬程修正矩阵hII中与REf和REg相对应的原泵扬程修正值。
CN202011464956.9A 2020-12-14 2020-12-14 一种基于仿真结果修正的离心泵扬程值预测方法 Active CN112560361B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011464956.9A CN112560361B (zh) 2020-12-14 2020-12-14 一种基于仿真结果修正的离心泵扬程值预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011464956.9A CN112560361B (zh) 2020-12-14 2020-12-14 一种基于仿真结果修正的离心泵扬程值预测方法

Publications (2)

Publication Number Publication Date
CN112560361A CN112560361A (zh) 2021-03-26
CN112560361B true CN112560361B (zh) 2022-04-26

Family

ID=75064211

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011464956.9A Active CN112560361B (zh) 2020-12-14 2020-12-14 一种基于仿真结果修正的离心泵扬程值预测方法

Country Status (1)

Country Link
CN (1) CN112560361B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113297809B (zh) * 2021-05-06 2022-04-26 湘潭大学 一种离心泵叶轮口环间隙泄漏量的仿真估算方法
CN116259428A (zh) * 2022-12-12 2023-06-13 中广核研究院有限公司 反应堆一回路中流量确定方法、装置和计算机设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189472A (ja) * 1994-10-31 1996-07-23 Meidensha Corp ポンプ吐出量演算方法および装置
CN107315858A (zh) * 2017-06-01 2017-11-03 哈尔滨工程大学 一种面向流网仿真应用基于泵内特性的离心泵参数化仿真方法
CN108506223A (zh) * 2018-03-05 2018-09-07 浙江理工大学 基于弱可压流动分析的大功率离心泵水力性能预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189472A (ja) * 1994-10-31 1996-07-23 Meidensha Corp ポンプ吐出量演算方法および装置
CN107315858A (zh) * 2017-06-01 2017-11-03 哈尔滨工程大学 一种面向流网仿真应用基于泵内特性的离心泵参数化仿真方法
CN108506223A (zh) * 2018-03-05 2018-09-07 浙江理工大学 基于弱可压流动分析的大功率离心泵水力性能预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《液货泵输送粘液的性能仿真与经验修正对比分析》;沈飞等;《船海工程》;20170228;第83-86页 *
《离心泵理论扬程修正系数的再研究》;何希杰等;《水泵技术》;19970430;第13-15，33页 *

Also Published As

Publication number Publication date
CN112560361A (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN112560361B (zh) 一种基于仿真结果修正的离心泵扬程值预测方法
CN112268011B (zh) 多级轴流压气机气动性能预估方法和装置
CN103226635B (zh) 基于三维动网格的旋转叶轮机械非定常流场的计算方法
CN112032032B (zh) 湿冷机组开式循环水泵变频运行方式的寻优方法
CN107203364B (zh) 一种用于压气机全工况特性的预测和辨识方法
CN111400941A (zh) 一种叶片泵内部回流及回流旋涡空化的数值预测方法
CN115552150A (zh) 确定静液压机当前磨损状况的方法
CN110276115B (zh) 基于燃机叶片型线参数的气路故障诊断方法
CN112412829A (zh) 基于物质平均大涡模拟的离心泵水力性能测量方法
CN108009382A (zh) 离心压缩机设计***
CN116595874A (zh) 叶轮机械性能预测模型参数优化方法及装置、存储介质
CN104895777B (zh) 一种确定压缩机动态机械效率的方法及装置
CN105604953A (zh) 一种基于振动优化的多级泵叶轮错位布置方法
CN115935833A (zh) 一种获取压气机特性的方法和装置
CN113359450B (zh) 一种阀门流量特性曲线拟合方法及***
JP2023024874A (ja) 管理システム、管理方法、および管理プログラム
CN104361258A (zh) 流体机械***变速调节相对阀门调节的节能效果检测方法
CN115114865A (zh) 流场数据计算方法、装置、计算设备及存储介质
Gileva et al. The study of Y+ influence on the results of ANSYS CFX flow simulation in the centrifugal compressor radial inlet
CN112580238B (zh) 一种基于仿真结果修正的离心泵效率值预测方法
CN114810728A (zh) 一种多路阀数字孪生模型及其建模方法
CN107299904A (zh) 一种用于描述离心泵正转逆流制动工况外特性的数学模型
Neelambika et al. CFD analysis of mixed flow impeller
CN114297793B (zh) 海水淡化泵叶轮结构多学科优化设计方法
CN114492247B (zh) 基于叶片欧拉扬程标准函数的叶轮高效低振动优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231101

Address after: Baimao Lake Industrial Concentration Zone, Tianmen City, Hubei Province 431700

Patentee after: HUBEI TIANMEN YONGQIANG PUMP INDUSTRY CO.,LTD.

Address before: 430070 room 01, R & D No. 3, 4 / F, building C5, phase III, Rongke Zhigu industrial project, Liqiao village, Hongshan District, Wuhan City, Hubei Province

Patentee before: Zhongou (Hubei) Intellectual Property Service Co.,Ltd.

Effective date of registration: 20231101

Address after: 430070 room 01, R & D No. 3, 4 / F, building C5, phase III, Rongke Zhigu industrial project, Liqiao village, Hongshan District, Wuhan City, Hubei Province

Patentee after: Zhongou (Hubei) Intellectual Property Service Co.,Ltd.

Address before: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee before: Dragon totem Technology (Hefei) Co.,Ltd.

Effective date of registration: 20231101

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Address before: 411105 College of mechanical engineering, Xiangtan University, yanggutang street, Yuhu District, Xiangtan City, Hunan Province

Patentee before: XIANGTAN University