CN112613134B - 一种基于涡流分布的阀体结构优化方法 - Google Patents

一种基于涡流分布的阀体结构优化方法 Download PDF

Info

Publication number
CN112613134B
CN112613134B CN202011448557.3A CN202011448557A CN112613134B CN 112613134 B CN112613134 B CN 112613134B CN 202011448557 A CN202011448557 A CN 202011448557A CN 112613134 B CN112613134 B CN 112613134B
Authority
CN
China
Prior art keywords
valve
model
valve body
vortex
flow channel
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
CN202011448557.3A
Other languages
English (en)
Other versions
CN112613134A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202011448557.3A priority Critical patent/CN112613134B/zh
Publication of CN112613134A publication Critical patent/CN112613134A/zh
Application granted granted Critical
Publication of CN112613134B publication Critical patent/CN112613134B/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/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16K—VALVES; TAPS; COCKS; ACTUATING-FLOATS; DEVICES FOR VENTING OR AERATING
    • F16K27/00—Construction of housing; Use of materials therefor
    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Mechanical Engineering (AREA)
  • Lift Valve (AREA)

Abstract

本发明公开了一种基于涡流分布的阀体结构优化方法，包括以下步骤：首先根据阀门的设计文件进行三维建模，建立完整的阀门实体模型，通过进行阀门的流道抽取，建立阀门内部流道模型；通过网格划分，通过数值模拟方法，利用计算流体力学软件进行阀内的流场计算；通过对计算结果的后处理，实现阀体流场的可视化，确定阀内的涡流分布；通过对阀体流道特征参数的量化定义和结构参数优化，达到消除涡流和提高流通能力的效果。本发明能够优化消除阀内涡流，增大阀门的流通能力，避免了繁琐复杂的理论计算与实验校正，大大简化了阀体的设计过程，为消除阀内涡流和提高阀门流通能力提供了新的思路。

Description

一种基于涡流分布的阀体结构优化方法
技术领域
本发明属于工业阀门设计领域，具体涉及一种基于涡流分布的阀体结构优化方法。
背景技术
阀门是管路流体输送系统中的控制部件，用来改变通路截断和介质流动方向，具有截止、导流、防止逆流、调节、稳压。分流或溢流泄压等功能。阀门的流量系数是衡量阀门流通能力的指标，流量系数值越大说明流体流过阀门时的压力损失越小。国外工业发达国家的阀门生产厂家大多把不同压力等级、不同类型和不同公称通径的流量系数值列入产品样本，供设计部门和使用单位选用。流量系数值随阀门的尺寸、形式、结构而变化，不同类型和不同规格的阀门都要分别进行试验，才能确定该种阀门的流量系数值。涡流是影响阀门流通能力和调节精度的关键参数，在阀门性能改善进程中一般通过消除涡流的存在来提高阀门的流通能力。在阀门的结构基本相似的前提下，通过消除阀内存在的涡流数量，提高阀门的流通能力，从而达到改善阀门性能的目的，这为阀体的结构改善设计提供了新的思路。
数值模拟，又称为数值分析方法，是一种依靠计算机程序获得符合精度需要的数学模型近似解的方法。将数值模拟运用于阀门设计，一方面可以在理论计算中尽量减少简化和假设，从而使计算结果尽量贴近真实解，进行避免了反复试算与校正；另一方面，几乎不需要进行实验校正，从而节省了大量的经济成本。当前，虽然将数值模拟运用于阀门设计已有一定的基础，但是未有针对涡流分布的阀体结构优化方法。
发明内容
本发明的目的在于解决现有技术中的缺陷，并提供了一种基于涡流分布的阀体结构优化方法，通过改变流道特征系数，来改善阀内的流体流动，减小阀内涡流的数量，提高阀门的流通能力。
本发明所采用的具体技术方案如下：
本发明提供了一种基于涡流分布的阀体结构优化方法，具体如下：
S1：根据阀门的设计文件，通过三维建模软件构建待优化结构阀门的三维结构模型，再根据所述三维结构模型建立阀门实体模型；根据所述阀门实体模型的内部流道，抽取出阀门内部流道模型；通过网格划分软件，将所述阀门内部流道模型离散为内部流道网格模型；
根据所述设计文件，通过计算流体力学软件设置所述内部流道网格模型的进出口边界条件；同时，根据设计文件对所述内部流道网格模型占据的全部空间赋予流体属性；
S2：在S1的基础上，通过计算流体力学软件，选择相应的湍流模型并设置参数后进行计算；当计算收敛后，选择计算结果中的速度结果，以阀门的纵剖面为对象，绘制该纵剖面所在处的速度云图；根据所述速度云图绘制流线，并根据流线分布，确定速度云图上的涡流位置和数量；
S3：根据S2中速度云图上的涡流位置和数量，择一执行S31或S32；
S31：若S2中速度云图上没有涡流存在，则该阀体结构无需优化，结束优化过程；
S32：若S2中速度云图上有涡流存在，则过每个所述涡流的中心点作与其距离最近阀体内壁面的垂直线；将所述垂直线与阀体内壁面的交汇点记为点A，与阀体中轴线的交汇点记为点B；将点A和点B之间构成的线段AB记为与涡流相关的流道特征系数RAB；
S4：根据S32中与涡流中心点距离最近的阀体内壁面类型，择一执行S41～S45或S51～S54；
S41：若S32中与涡流中心点距离最近的阀体内壁面为曲面，且其曲率中心点在阀体内部时，则将该阀体内壁面初始的流道特征系数记为RABi，其中i的初始值为1；通过三维建模软件将该阀体内壁面变为平面，此时的流道特征系数记为RABj，其中j的初始值为2；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RABj时的涡流是否消失；
S42：若S41中的涡流消失，则结束优化过程；
S43：若S41中的涡流未消失，则定义ΔRij＝RABi-RABj，并将S41中原来的流道特征系数RABj变为RABj+1＝RABj-ΔRij；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RABj+1时的涡流是否消失；
S44：若S43中的涡流消失，则结束优化过程；
S45：若S43中的涡流未消失，则对i和j均执行加1操作，重新根据S43的操作继续优化流道特征系数，直至涡流消失；
S51：若S32中与涡流中心点距离最近的阀体内壁面为平面时，将该阀体内壁面初始的流道特征系数记为R’ABi，其中i的初始值为1；
S52：将该流道特征系数变为R’ABi+1＝R’ABi-0.05R’ABi；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为R’ABi+1时的涡流是否消失；
S53：若S52中的涡流消失，则结束优化过程；
S54：若S52中的涡流未消失，则对i执行加1操作，根据S52的操作继续优化流道特征系数，直至涡流消失。
作为优选，所述S1中，当阀门的结构相对其纵剖面几何对称时，所述三维结构模型为完整的阀门模型或相对于纵剖面的1/2阀门模型；当阀门的结构相对其纵剖面几何不对称时，所述三维结构模型为完整的阀门模型；
所述阀门的纵剖面为阀门进口面中心点和出口面中心点的连线与阀体中轴线所在的平面。
进一步的，所述S1中，当三维结构模型为完整的阀门模型时，阀门内部流道模型为阀门内壁面、流体入口面与流体出口面包围成的三维几何空间结构；当三维结构模型为相对于纵剖面的1/2阀门模型时，阀门内部流道模型为阀门内壁面、流体入口面、流体出口面与阀门纵剖面包围成的三维几何空间结构；
进一步的，所述阀门内壁面为阀门内部与流体的接触面；所述流体入口面为阀门内部流体流动空间与上游管道中流体流动空间的分隔面；所述流体出口面为套筒阀内部流体流动空间与下游管道中流体流动空间的分隔面。
作为优选，所述S1中，进出口边界条件包括压力进口边界条件、压力出口边界条件、速度进口边界条件和速度出口边界条件中的一种或多种。
进一步的，所述S1中，在赋予流体属性之前，若三维结构模型为完整的阀门模型时，将所述内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件；若三维结构模型为相对于纵剖面的1/2阀门模型时，将所述内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件，同时将阀门纵剖面所在处设置对称边界条件。
作为优选，所述三维建模软件为Solidworks、Creo、Inventor、UG/NX、CATIA、ANSYSWorkbench DesignModeler或ANSYS Workbench SpaceClaim中的一种。
作为优选，所述网格划分软件为ICEM CFD、HyperMesh、TGrid、PointWise、ANSA、GridPro或ANSYS Workbench Mesh中的一种。
作为优选，所述计算流体力学软件为ANSYS Fluent、ANSYS CFX、STAR-CD、STAR-CCM、NUMECA或OpenFOAM中的一种。
作为优选，所述湍流模型为Spalart-Allmaras模型、Standard k-ε模型、RNG k-ε模型、Realizable k-ε模型、Standard k-ω模型、BSL k-ω模型或SST k-ω模型中的一种。
本发明相对于现有技术而言，具有以下有益效果：
(1)本发明在无需反复理论计算与实验校正的前提下实现了对阀体结构的改善设计，大大推动了工业阀体的结构优化。
(2)本发明提出了针对阀门流道特征的流道特征系数，针对阀体的流道结构参数进行了量化的定义，将阀体的流道结构特征用流道特征系数来表征，通过改变流道特征系数，来改善阀内的流体流动，从而减少了阀内涡流的数量，提高了阀门的流通能力。
(3)本发明基于数值模拟方法，使阀内的流体流动实现了可视化，能够直观的确定阀内的涡流分布；同时，通过改变阀体的流道特征系数，消除了阀内的涡流，增大了阀门的流通能力，避免了繁琐复杂的理论计算与实验校正，大大简化了阀体的设计过程，为消除阀内涡流和提高阀门流通能力提供了新的思路。
附图说明
图1为实施例1中阀体的二维图纸；
图2为实施例1中根据阀体的二维模型建立的三维结构模型；
图3为实施例1中根据三维结构模型建立的内部流道网格模型；
图4为实施例1中初始结构下阀体结构的速度云图；
图5为图4中涡流位置的示意图；
图6为实施例1中流道特征系数的优化过程示意图；
图7为实施例1中第一次结构优化后的速度云图；
图8为实施例1中第二次结构优化后的速度云图。
具体实施方式
为使本发明实施例的目的、技术方法和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方法进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
针对阀体的流道结构，本发明提出了针对阀门流道特征的阀体流道特征系数，针对阀体的流道结构参数进行了量化的定义，将阀体的流道结构特征用阀体流道特征系数来表征，通过改变流道特征系数，来改善阀内的流体流动，减小阀内涡流的数量，提高阀门的流通能力。
该阀体结构优化方法的具体操作如下：
S1：根据步骤1)～7)依次进行，具体如下：
1)根据阀门的设计文件，通过三维建模软件构建待优化结构阀门的三维结构模型。在构建三维结构模型的过程中，当阀门的结构相对其纵剖面几何对称时，三维结构模型可以采用完整的阀门模型或者相对于纵剖面而言的1/2阀门模型。当阀门的结构相对其纵剖面几何不对称时，三维结构模型只能采用完整的阀门模型。
其中，阀门的纵剖面指的是阀门进口面中心点和出口面中心点的连线与阀体中轴线所在的平面。
2)根据三维结构模型建立阀门实体模型。若在构建三维结构模型时采用的是完整的阀门模型，则阀门实体模型也为完整的模型。若在构建三维结构模型时采用的是以纵剖面为边界的1/2阀门模型，则阀门实体模型也为一半的阀门模型。
3)根据阀门实体模型的内部流道，抽取出阀门内部流道模型。若在构建三维结构模型时采用的是完整的阀门模型，则此时构建的阀门内部流道模型为阀门内壁面、流体入口面与流体出口面包围成的三维几何空间结构。若在构建三维结构模型时采用的是以纵剖面为边界的1/2阀门模型，则此时构建的阀门内部流道模型为阀门内壁面、流体入口面、流体出口面与阀门纵剖面包围成的三维几何空间结构。
其中，阀门内壁面指的是阀门内部与流体的接触面；流体入口面指的是阀门内部流体流动空间与上游管道中流体流动空间的分隔面；流体出口面指的是套筒阀内部流体流动空间与下游管道中流体流动空间的分隔面。
4)将获得的阀门内部流道模型导入网格划分软件中，并将其离散为内部流道网格模型。
5)将获得的内部流道网格模型导入计算流体力学软件中，根据原始的设计文件，设置内部流道网格模型的进出口边界条件。
在该过程中，若文件中规定阀门入口面和出口面的压力值，则对流体入口面设置压力入口边界条件，压力值为设计文件中所述阀门的入口压力，对流体出口面设置压力出口边界条件，压力值为设计文件中所述阀门的出口压力；若文件中规定阀门入口面和出口面的速度值，则对流体入口面设置速度入口边界条件，速度值为设计文件中所述阀门的入口速度，对流体出口面设置速度出口边界条件，压力值为设计文件中所述阀门的出口速度；若文件中规定阀门入口面的压力值和出口面的速度值，则对流体入口面设置压力入口边界条件，压力值为设计文件中所述阀门的入口压力，对流体出口面设置速度出口边界条件，速度值为设计文件中所述阀门的出口速度；若文件中规定阀门入口面的速度和出口面的压力值，对流体入口面设置为速度入口边界条件，速度值为设计文件中所述阀门的速度值，对流体出口面设置为压力出口边界条件，压力值为设计文件中所述阀门的压力值。
6)若在构建三维结构模型时采用的是完整的阀门模型，则将内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件。若在构建三维结构模型时采用的是以纵剖面为边界的1/2阀门模型，则将内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件，同时将阀门纵剖面所在处设置对称边界条件。
7)对内部流道网格模型占据的全部空间赋予流体属性，其中，流体密度与流体粘性数值等物性参数均根据阀门的原始设计文件确定。
S2：在S1的基础上，通过计算流体力学软件，选择相应的湍流模型并设置其余的参数后开始进行计算。当计算收敛后，选择计算结果中的速度结果，以阀门的纵剖面为对象，绘制该纵剖面所在处的速度云图。在阀门纵剖面速度云图的基础上，绘制纵剖面上的流线，并根据流线分布，确定速度云图上的涡流位置和数量。
其中，涡流在速度云图上的特征为漩涡状，其与正常流动区域有明显区别，能够很容易被识别出。
S3：根据S2中速度云图上的涡流位置和数量，择一执行S31或S32，具体如下：
S31：若S2中的速度云图上没有涡流存在，则现有的该阀体结构无需优化，可以直接结束优化过程。
S32：若S2中的速度云图上有涡流存在，则过每个涡流的中心点作一条与其距离最近阀体内壁面的垂直线。该壁面在几何上为具有相同参数的一部分壁面，如具有相同锥度的圆锥面或者具有相同半径的圆柱面。
将该垂直线与阀体内壁面的交汇点记为点A，与阀体中轴线的交汇点记为点B。将点A和点B之间构成的线段AB记为与涡流相关的流道特征系数RAB。
S4：根据S32中与涡流中心点距离最近的阀体内壁面类型，择一执行S41～S45或S51～S54，具体如下：
S41：若S32中与涡流中心点距离最近的阀体内壁面为曲面，且其曲率中心点在阀体内部时，则将该阀体内壁面初始的流道特征系数记为RAB1。通过三维建模软件将该阀体内壁面变为平面，此时的流道特征系数记为RAB2。由于不同的流道特征系数对应的阀体结构也有所不同，在流道特征系数为RAB2的情况下，根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为RAB2时的涡流是否消失。
S42：若S41中的涡流消失，则结束该优化过程，此时流道特征系数对应下的阀体结构即为最优结构。
若S41中的涡流未消失，则定义ΔR12＝RAB1-RAB2，并将S41中原来的流道特征系数RAB2变为RAB3＝RAB2-ΔR12。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为RAB3时的涡流是否消失。
S43：若S42中的涡流消失，则结束该优化过程，此时流道特征系数对应下的阀体结构即为最优结构。
若S42中的涡流未消失，则定义ΔR23＝RAB2-RAB3，并将S42中原来的流道特征系数RAB3变为RAB4＝RAB3-ΔR23。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为RAB4时的涡流是否消失。
…
S4j：若S4i中的涡流消失，则结束该优化过程，此时流道特征系数对应下的阀体结构即为最优结构。
若S4i中的涡流未消失，则定义ΔRij＝RABi-RABj，并将S4i中原来的流道特征系数RABj变为RABj+1＝RABj-ΔRij。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为RABj+1时的涡流是否消失。其中，i的初始值为1，j的初始值为2。
S4j+1：若S4j中的涡流消失，则结束优化过程。
若S4j中的涡流未消失，则对i和j均执行加1操作，重新根据S4j的操作继续优化流道特征系数，直至涡流消失。
S51：若S32中与涡流中心点距离最近的阀体内壁面为平面时，则将该阀体内壁面初始的流道特征系数记为R’AB1。
S52：将该流道特征系数变为R’AB2＝R’AB1-0.05R’AB1。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为R’AB2时的涡流是否消失。
S53：若S52中的涡流消失，则结束优化过程。
若S52中的涡流未消失，则将流道特征系数变为R’AB3＝R’AB2-0.05R’AB2。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为R’AB3时的涡流是否消失。
…
S5i+1：若S5i中的涡流消失，则结束优化过程。
若S5i中的涡流未消失，则将流道特征系数变为R’ABi+1＝R’ABi-0.05R’ABi。根据此时的阀体结构重复S1～S2的操作，来判断当流道特征系数为R’ABi+1时的涡流是否消失。其中，i的初始值为1。
S5i+2：若S5 i+1中的涡流消失，则结束优化过程。
若S5i+1中的涡流未消失，则对i执行加1操作，根据S5i+1的操作继续优化流道特征系数，直至涡流消失。
在实际应用过程中，三维建模软件可以采用Solidworks、Creo、Inventor、UG/NX、CATIA、ANSYS Workbench DesignModeler或ANSYS Workbench SpaceClaim中的一种。网格划分软件可以采用ICEM CFD、HyperMesh、TGrid、PointWise、ANSA、GridPro或ANSYSWorkbench Mesh中的一种。计算流体力学软件可以采用ANSYS Fluent、ANSYS CFX、STAR-CD、STAR-CCM、NUMECA或OpenFOAM中的一种。湍流模型可以采用Spalart-Allmaras模型、Standard k-ε模型、RNG k-ε模型、Realizable k-ε模型、Standard k-ω模型、BSL k-ω模型或SST k-ω模型中的一种。
实施例1
S1：如图1所示，为待优化结构的减压阀阀体的二维图纸。根据阀门的设计文件，通过三维建模软件构建该阀门的三维结构模型，构建的减压阀三维结构模型如图2所示。由图2可以看出，该阀体结构是关于阀体的纵向剖切面对称的，故可建立关于剖切面对称的1/2阀体模型。在三维建模软件中，建立该减压阀的实体模型。根据该阀门实体模型的内部流道，抽取出阀门内部流道模型。由于本实施例中阀门三维结构模型为相对于纵剖面的1/2阀门模型，即阀门内部流道模型为阀门内壁面、流体入口面、流体出口面与阀门纵剖面包围成的三维几何空间。通过网格划分软件，将阀门内部流道模型离散为内部流道网格模型，如图3所示，为利用FLUENT自带的网格划分软件进行划分的内部流道网格模型。
将得到的内部流道网格模型导入计算流体力学软件FLUENT软件中，并根据设计文件，将内部流道网格模型的压力进口边界条件设置为入口压力10MPa，压力出口边界条件设置为入口压力1MPa。将内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件，同时将阀门纵剖面所在处设置对称边界条件。根据设计文件对内部流道网格模型占据的全部空间赋予流体属性。
S2：在计算流体力学FLUENT软件中，选择相应的湍流模型并设置其余参数后开始计算。计算收敛后，在所述计算流体力学FLUENT软件或者其配套后处理软件Tecplot中进行计算结构的后处理；所述的后处理步骤为：选择计算结果中的速度结果，选择所述阀门的纵剖面为对象，绘制在纵剖面上的速度云图中，如图4所示。在速度云图的基础上，绘制纵剖面上的流线，根据流线分布，确定涡流的位置，涡流在流线图上的特征为漩涡状，其与正常流动区域有明显区别，如图5所示。
S3：如图5所示，共有涡流3个，现选取涡流1举例进行优化过程。
如图6所示，过涡流1的中心点作与其距离最近阀体内壁面的垂直线。将垂直线与阀体内壁面的交汇点记为点A，与阀体中轴线的交汇点记为点B，将点A和点B之间构成的线段AB记为与涡流相关的流道特征系数RAB。
S4：由于与S3中涡流中心点距离最近的阀体内壁面为曲面，且其曲率中心点在阀体内部时，则将该阀体内壁面初始的流道特征系数记为RAB1。通过三维建模软件将该阀体内壁面变为平面，此时的流道特征系数记为RAB2。在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RAB2时的涡流是否消失，第一次结构优化后的速度云图如图7所示。
从图7可以看出，该涡流仍然存在。则定义ΔR12＝RAB1-RAB2，将流道特征系数变为RAB3，RAB3的值为RAB3＝RAB2+ΔR12。在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RAB3时的涡流是否消失，第二次结构优化后的速度云图如图8所示。
从图8中可以看出，该结构下的涡流已经全部消失，此时的阀体结构即为最佳结构，可以直接输出，无需继续优化。
本发明提供了一种基于涡流分布的阀体结构优化方法，能够在进行阀体结构优化的基础上，通过消除阀内存在的涡流，提高阀门的流通能力，从而达到改善阀门性能的目的，同时避免反复的理论计算与实验校正。
以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (10)

1.一种基于涡流分布的阀体结构优化方法，其特征在于，具体如下：
S1：根据阀门的设计文件，通过三维建模软件构建待优化结构阀门的三维结构模型，再根据所述三维结构模型建立阀门实体模型；根据所述阀门实体模型的内部流道，抽取出阀门内部流道模型；通过网格划分软件，将所述阀门内部流道模型离散为内部流道网格模型；
根据所述设计文件，通过计算流体力学软件设置所述内部流道网格模型的进出口边界条件；同时，根据设计文件对所述内部流道网格模型占据的全部空间赋予流体属性；
S2：在S1的基础上，通过计算流体力学软件，选择相应的湍流模型并设置参数后进行计算；当计算收敛后，选择计算结果中的速度结果，以阀门的纵剖面为对象，绘制该纵剖面所在处的速度云图；根据所述速度云图绘制流线，并根据流线分布，确定速度云图上的涡流位置和数量；
S3：根据S2中速度云图上的涡流位置和数量，择一执行S31或S32；
S31：若S2中速度云图上没有涡流存在，则该阀体结构无需优化，结束优化过程；
S32：若S2中速度云图上有涡流存在，则过每个所述涡流的中心点作与其距离最近阀体内壁面的垂直线；将所述垂直线与阀体内壁面的交汇点记为点A，与阀体中轴线的交汇点记为点B；将点A和点B之间构成的线段AB记为与涡流相关的流道特征系数RAB；
S4：根据S32中与涡流中心点距离最近的阀体内壁面类型，择一执行S41～S45或S51～S54；
S41：若S32中与涡流中心点距离最近的阀体内壁面为曲面，且其曲率中心点在阀体内部时，则将该阀体内壁面初始的流道特征系数记为RABi，其中i的初始值为1；通过三维建模软件将该阀体内壁面变为平面，此时的流道特征系数记为RABj，其中j的初始值为2；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RABj时的涡流是否消失；
S42：若S41中的涡流消失，则结束优化过程；
S43：若S41中的涡流未消失，则定义ΔRij＝RABi-RABj，并将S41中原来的流道特征系数RABj变为RABj+1＝RABj-ΔRij；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为RABj+1时的涡流是否消失；
S44：若S43中的涡流消失，则结束优化过程；
S45：若S43中的涡流未消失，则对i和j均执行加1操作，重新根据S43的操作继续优化流道特征系数，直至涡流消失；
S51：若S32中与涡流中心点距离最近的阀体内壁面为平面时，将该阀体内壁面初始的流道特征系数记为R’ABi，其中i的初始值为1；
S52：将该流道特征系数变为R’ABi+1＝R’ABi-0.05R’ABi；在该阀门结构的情况下重复S1～S2，判断当流道特征系数为R’ABi+1时的涡流是否消失；
S53：若S52中的涡流消失，则结束优化过程；
S54：若S52中的涡流未消失，则对i执行加1操作，根据S52的操作继续优化流道特征系数，直至涡流消失。
2.根据权利要求1所述的阀体结构优化方法，其特征在于，所述S1中，当阀门的结构相对其纵剖面几何对称时，所述三维结构模型为完整的阀门模型或相对于纵剖面的1/2阀门模型；当阀门的结构相对其纵剖面几何不对称时，所述三维结构模型为完整的阀门模型；
所述阀门的纵剖面为阀门进口面中心点和出口面中心点的连线与阀体中轴线所在的平面。
3.根据权利要求2所述的阀体结构优化方法，其特征在于，所述S1中，当三维结构模型为完整的阀门模型时，阀门内部流道模型为阀门内壁面、流体入口面与流体出口面包围成的三维几何空间结构；当三维结构模型为相对于纵剖面的1/2阀门模型时，阀门内部流道模型为阀门内壁面、流体入口面、流体出口面与阀门纵剖面包围成的三维几何空间结构。
4.根据权利要求3所述的阀体结构优化方法，其特征在于，所述阀门内壁面为阀门内部与流体的接触面；所述流体入口面为阀门内部流体流动空间与上游管道中流体流动空间的分隔面；所述流体出口面为套筒阀内部流体流动空间与下游管道中流体流动空间的分隔面。
5.根据权利要求1所述的阀体结构优化方法，其特征在于，所述S1中，进出口边界条件包括压力进口边界条件、压力出口边界条件、速度进口边界条件和速度出口边界条件中的一种或多种。
6.根据权利要求3所述的阀体结构优化方法，其特征在于，所述S1中，在赋予流体属性之前，若三维结构模型为完整的阀门模型时，将所述内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件；若三维结构模型为相对于纵剖面的1/2阀门模型时，将所述内部流道网格模型中阀门内壁面所在处设置无滑移壁面边界条件，同时将阀门纵剖面所在处设置对称边界条件。
7.根据权利要求1所述的阀体结构优化方法，其特征在于，所述三维建模软件为Solidworks、Creo、Inventor、UG/NX、CATIA、ANSYS Workbench DesignModeler或ANSYSWorkbench SpaceClaim中的一种。
8.根据权利要求1所述的阀体结构优化方法，其特征在于，所述网格划分软件为ICEMCFD、HyperMesh、TGrid、PointWise、ANSA、GridPro或ANSYS Workbench Mesh中的一种。
9.根据权利要求1所述的阀体结构优化方法，其特征在于，所述计算流体力学软件为ANSYS Fluent、ANSYS CFX、STAR-CD、STAR-CCM、NUMECA或OpenFOAM中的一种。
10.根据权利要求1所述的阀体结构优化方法，其特征在于，所述湍流模型为Spalart-Allmaras模型、Standard k-ε模型、RNG k-ε模型、Realizable k-ε模型、Standard k-ω模型、BSL k-ω模型或SST k-ω模型中的一种。
CN202011448557.3A 2020-12-09 2020-12-09 一种基于涡流分布的阀体结构优化方法 Active CN112613134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011448557.3A CN112613134B (zh) 2020-12-09 2020-12-09 一种基于涡流分布的阀体结构优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011448557.3A CN112613134B (zh) 2020-12-09 2020-12-09 一种基于涡流分布的阀体结构优化方法

Publications (2)

Publication Number Publication Date
CN112613134A CN112613134A (zh) 2021-04-06
CN112613134B true CN112613134B (zh) 2022-05-13

Family

ID=75233395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011448557.3A Active CN112613134B (zh) 2020-12-09 2020-12-09 一种基于涡流分布的阀体结构优化方法

Country Status (1)

Country Link
CN (1) CN112613134B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255246B (zh) * 2021-05-25 2023-10-27 核电运行研究(上海)有限公司 一种主蒸汽阀站中主蒸汽安全阀的开启压力设置方法
CN113626956A (zh) * 2021-08-20 2021-11-09 浙江理工大学 一种基于流固耦合分析的蝶阀振动预测方法
CN113722850B (zh) * 2021-08-25 2023-04-18 浙江大学 一种基于冲蚀磨损的阀体壁厚轻量化设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104662235A (zh) * 2012-08-02 2015-05-27 海德鲁国际公司 配置涡旋流控制装置的方法和涡旋流控制装置
CN105431641A (zh) * 2013-05-17 2016-03-23 卡梅伦国际有限公司 流调节器及优化方法
CN105624646A (zh) * 2014-11-26 2016-06-01 朗姆研究公司 借助可再入流路径的阀歧管盲管消除
CN108573089A (zh) * 2017-03-10 2018-09-25 通用电气公司 利用3d cad点云自动创建流体模型的系统和方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104662235A (zh) * 2012-08-02 2015-05-27 海德鲁国际公司 配置涡旋流控制装置的方法和涡旋流控制装置
CN105431641A (zh) * 2013-05-17 2016-03-23 卡梅伦国际有限公司 流调节器及优化方法
CN105624646A (zh) * 2014-11-26 2016-06-01 朗姆研究公司 借助可再入流路径的阀歧管盲管消除
CN108573089A (zh) * 2017-03-10 2018-09-25 通用电气公司 利用3d cad点云自动创建流体模型的系统和方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"活塞式调节阀流场特性分析与结构优化";邓君;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20160615(第06期);第C029-108页 *
"滑动式调节阀流场及调节性能研究";胡畏;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20190115(第01期);第C029-365页 *

Also Published As

Publication number Publication date
CN112613134A (zh) 2021-04-06

Similar Documents

Publication Publication Date Title
CN112613134B (zh) 一种基于涡流分布的阀体结构优化方法
Yao et al. Parametric design and optimization of high speed train nose
CN110534776B (zh) 一种基于fluent仿真的燃料电池多堆分配内腔结构优化的方法
CN106383955A (zh) 管道设计中应力分析与三维模型的数据相互转换方法
CN111767663A (zh) 一种基于cfd仿真的阀门流量系数便捷计算方法
CN110750855B (zh) 一种外形定尺寸限制下的蜗壳型线设计方法
CN112100761A (zh) 一种火箭发动机管路的动响应分析和减振优化设计方法
CN109063363A (zh) 一种基于数据库的喷水推进器的优化设计方法
CN108062444B (zh) 弯管参数确定方法
CN109214136A (zh) 一种轴流泵叶轮叶顶间隙空化数值预测方法
CN110196987B (zh) 基于代理模型的风道结构尺寸优化方法
CN108804791B (zh) 一种适用于埋入式进气道布局的飞行器参数化方法
CN109779962B (zh) 一种超低比转速离心泵叶轮的二元粘性水力设计方法
CN113503262B (zh) 一种两级增压系统高低压级涡轮的仿真方法
CN115358101A (zh) 一种基于声速解和特征线逆推的喷管设计方法
CN104537170A (zh) 用于泵站三维流场计算的网格划分方法及装置
CN112163288B (zh) 基于大涡模拟的流体机械叶片优化设计方法
CN113722850B (zh) 一种基于冲蚀磨损的阀体壁厚轻量化设计方法
CN107423481B (zh) 基于渗透边界反问题方法的全三维内转式进气道设计方法
CN110069808A (zh) 基于壁面压力梯度特征的s形进气道的优化设计方法
CN113033116B (zh) 一种能克服不平衡力矩的活塞式阀芯底面开孔优化方法
CN113536711B (zh) 一种基于田口法的橄榄形流量计结构优化设计方法
CN116305700B (zh) 一种用于长距离输水管网的水锤评估和消除模型
CN112199831B (zh) 一种高流通能力排气道设计方法
CN116796602A (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