CN113255185B - 一种阀门调节特性自动仿真方法 - Google Patents
一种阀门调节特性自动仿真方法 Download PDFInfo
- Publication number
- CN113255185B CN113255185B CN202110560101.4A CN202110560101A CN113255185B CN 113255185 B CN113255185 B CN 113255185B CN 202110560101 A CN202110560101 A CN 202110560101A CN 113255185 B CN113255185 B CN 113255185B
- Authority
- CN
- China
- Prior art keywords
- valve
- module
- opening
- solving
- model
- 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
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/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
- 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)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及阀门设计领域,具体涉及一种阀门调节特性自动仿真方法,包括以下步骤:S1、建立阀门三维模型;S2、将阀门三维模型导入模型模块中,并进行阀门开度参数创建和流体域三维模型创建;S3、网格划分,并生成网格文件;S4、进入求解模块,对流体介质、湍流模型、求解器和求解控制参数进行设置;并创建输出参数;S5、进入参数集模块,设置阀门调节特性所需仿真的多个开度参数,并求解计算出输出参数;S6、求解完成后,读取参数集模块中的阀门各开度下的参数值,包括阀门开度参数值和输出参数;S7、按照流量系数计算公式计算出各开度流量系数,绘制出阀门调节特性曲线图。本方案降低了仿真人员的仿真操作时间及劳动强度,提高了仿真的效率。
Description
技术领域
本发明涉及阀门设计领域,具体涉及一种阀门调节特性自动仿真方法。
背景技术
阀门在现代化工工艺***中,起着调节介质压力、流量等参数的作用,是工艺***中重要的控制元件;而阀门调节特性则是工艺***参数设计的主要依据。在有限元技术成熟前,阀门调节特性仅能通过流量试验获取;有限元技术成熟后,可通过有限元仿真软件获取阀门调节特性曲线,缩短了阀门设计周期和节省了阀门试验成本。
中国专利申请号为CN201810742455.9的发明专利公开了一种调节阀内部尺寸设计方法,通过该专利可知,此专利提供了基于ANSYS平台的阀门Cv值模拟测试方法,解决了现有技术中调节阀内部尺寸与理论值误差较大且产品周期较长的问题,提高了内部尺寸的设计精度,缩短了产品研发周期。但是,此方法每次只能测试某个特定阀门开度的Cv值,若需获得阀门调节特性,则需从头开始对相关模型、模块进行更改、设置,对每个开度进行重复性测试,这样增加了仿真人员的仿真操作时间及劳动强度,仿真效率低下。
中国专利申请号为CN202010549961.3发明专利公开了一种阀门流量特性模拟实验方法,此专利通过建立三维流道模型、录制脚本、程序编写、动网格重构技术,实现了阀门流量特性的获取。但此方法存在以下缺点:(1)三维流道模型即阀门内部流体模型,这属于逆向建模,需要极强的三维建模能力和空间想象力,建模难度较大;(2)录制脚本及程序编写,需仿真人员具备一定编程能力及熟悉仿真软件程序编写特点,因此对仿真人员的能力具有极高的要求;(3)动网格重构技术是仿真软件对流道模型变动后自动网格划分技术,是动态过程,对计算机配置要求高且求解效率低下。
发明内容
本发明意在提供一种阀门调节特性自动仿真方法,以降低仿真人员的仿真操作时间及劳动强度,提高仿真的效率。
为达到上述目的,本发明采用如下技术方案:一种阀门调节特性自动仿真方法,包括以下步骤:
S1、使用三维软件建立阀门三维模型;
S2、在建好阀门三维模型之后,将阀门三维模型导入模型模块中,并进行阀门开度参数创建和流体域三维模型创建;
S3、将流体域三维模型进行网格划分,并生成网格文件;
S4、进入求解模块,对流体介质、湍流模型、求解器和求解控制参数进行设置;并创建输出参数:阀门进口压力、阀门出口压力及阀门进口流量;
S5、进入参数集模块,设置阀门调节特性所需仿真的多个开度参数,并求解计算出阀门进口压力、阀门出口压力及阀门进口流量等参数;
S6、求解完成后,读取参数集模块中的阀门各开度下的参数值,包括阀门开度参数值、阀门进口压力参数值、阀门出口压力参数值及阀门进口流量参数值;
S7、将阀门各开度下的参数值输入到数据处理工具中,按照流量系数计算公式计算出各开度流量系数,最后绘制阀门调节特性曲线图。
优选的,作为一种改进,S2中,还进行流场仿真几何清理,流场仿真几何清理包括倒角。通过仿真几何清理降低了网格数量,提高了计算效率。
优选的,作为一种改进,S3中,采用ANSYS mesh将流体域三维模型进行网格划分,通过设置全局网格尺寸以及细化节流处网格,完成网格划分并生成网格文件。
优选的,作为一种改进,S4中,流体介质为水、湍流模型为RNS K-e,求解器为Pressure-based,求解控制参数保持默认。
优选的,作为一种改进,S7中,流量系数计算公式为:
优选的,作为一种改进,S7中,数据处理工具为Excel。Excel在数据处理比较常用,无需借助其他高端的数据处理工具,操作简单方便易学。
优选的,作为一种改进,S2中,模型模块为DesignModeler模块。
优选的,作为一种改进,阀门的初始开度大于阀门的死区开度。因阀门结构特性,存在一定开启死区,当阀门开度大于阀门的死区开度时,阀门才会开启。
优选的,作为一种改进,模型模块、求解模块和参数集模块均位于Fluid Flow模块中。
优选的,作为一种改进,S2中,阀门三维模型转化为中间格式*.x_t,然后再导入模型模块中。
采用上述技术方案时,具有以下优点:
1、本方案基于有限元法,提出了一种阀门调节特性自动仿真方法,简单易行,仿真人员无需具备极强的三维建模能力和空间想象力,也无需具备一定编程能力及熟悉仿真软件程序编写特点,适用于具有不同仿真基础的仿真人员,对仿真人员的所具备的能力要求不高。
2、通过本方案,仿真人员只需完成阀门一个开度仿真分析设置及阀门开度参数、输出参数设置,其余开度流场仿真,软件依定义自动根据不同的开度而求解计算出相应的输出参数的值,无需对每个开度进行重复性测试,减少了仿真人员的仿真操作时间及劳动强度,仿真效率提高。
3、此外,每个开度流场仿真均采用稳态仿真,极大地降低了计算机配置,对计算机的配置没有过高的要求,从而提高了仿真效率。
附图说明
图1为本方法操作流程图。
图2为球阀三维模型图。
图3为阀门调节特性曲线图。
图4为开度18°下流体域三维模型图。
图5为开度18°下流体域网格图。
具体实施方式
下面通过具体实施方式进一步详细说明:
说明书附图中的附图标记包括:阀座1、阀杆2、球芯3、左阀体4、右阀体5、阀前管道6、阀后管道7。
一种阀门调节特性自动仿真方法,结合图1所示,包括以下步骤:
S1、建立阀门三维模型:根据现有的ZJHR-V DN100近似对数特性的二维图纸,运用三维软件建立阀门整机三维模型。请参阅图2,阀前阀后管道按照Cv(流量系数)测试标准建立,阀前2d、阀后6d(d为管道公称通径),只需要建立阀座1、阀杆2、球芯3、左阀体4、右阀体5、阀前管道6、阀后管道7。
S2、建立阀门开度参数及流体域模型:将S1中的阀门三维模型转化为中间格式*.x_t,导入Fluid Flow(Fluent)中的模型模块中,本实施例中模型模块为DesignModeler模块,建立阀门开度参数angle及创建流体域三维模型,由于此阀门为对称模型,采用一半模型,如图4所示。
S3、网格模型建立:进入Fluid Flow(Fluent)中的网格划分模块,对步骤S2中创建的流体域三维模型进行网格划分,通过设置全局网格尺寸、细化节流处网格,完成网格划分。对阀前、后管道采用六面体网格,阀体处采用四面体网格,六面体网格和四面体网格交界处网格节点自动对齐;对阀体处细小流体域采用网格细化,网格数量在100万左右,如图5所示。
S4、数值仿真设置及建立输出参数:进入Fluid Flow(Fluent)中的求解设置模块,设置流体介质为水、湍流模型RNS K-e、压力进口500000Pa和压力出口0pa,其求解器采用默认设置Pressure-based,求解控制参数保持默认。此外,建立求解输出参数:进口压力P_in、出口压力P_out及进口质量流量flow_in。
S5、设置阀门开度参数:进入Fluid Flow(Fluent)中的参数集模块,设置阀门调节特性所需仿真的其余开度参数值,并求解计算出各个开度参数下的输出参数:进口压力P_in、出口压力P_out及进口质量流量flow_in。设置阀门调节特性所需仿真的其余开度参数值时,因阀门结构特性,存在一定开启死区,阀门开度为18°时,阀门才开启,因此这里设置调节特性的阀门开度为18°-90°,且开度每隔10%(即9°)设置仿真求解点。
S6、求解完成后,读取参数集模块中的阀门各开度下的参数值,包括阀门开度参数angle、进口压力参数P_in、出口压力参数P_out及进口流量参数flow_in,具体计算的参数值如下表所示。
计算结果参数表:
S7、将S6中读取的数据,输入excel中,并按照公式①在Excel中编写流量系数计算公式,计算各开度Cv值,绘制调节特性曲线图,如图3所示。
以上所述的仅是本发明的实施例,方案中公知的具体技术方案和/或特性等常识在此未作过多描述。应当指出,对于本领域的技术人员来说,在不脱离本发明技术方案的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (9)
1.一种阀门调节特性自动仿真方法,其特征在于:包括以下步骤:
S1、使用三维软件建立阀门三维模型;
S2、在建好阀门三维模型之后,将阀门三维模型导入模型模块中,并进行阀门开度参数创建和流体域三维模型创建;
S3、将流体域三维模型进行网格划分,并生成网格文件;
S4、进入求解模块,对流体介质、湍流模型、求解器和求解控制参数进行设置;并创建输出参数:阀门进口压力、阀门出口压力及阀门进口流量;
S5、进入参数集模块,设置阀门调节特性所需仿真的多个开度参数,并求解计算出阀门进口压力、阀门出口压力及阀门进口流量参数,软件依定义自动根据不同的开度而求解计算出相应的输出参数的值,无需对每个开度进行重复性测试;
S6、求解完成后,读取参数集模块中的阀门各开度下的参数值,包括阀门开度参数值、阀门进口压力参数值、阀门出口压力参数值及阀门进口流量参数值;
S7、将阀门各开度下的参数值输入到数据处理工具中,按照流量系数计算公式计算出各开度流量系数,最后绘制阀门调节特性曲线图,流量系数计算公式为:
;式中,Cv为仿真流量系数,P1为阀门进口压力,单位为Pa;P2为阀门出口压力,单位为Pa;Q为质量流量,单位为kg/s。
2.根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:S2中,还进行流场仿真几何清理,流场仿真几何清理包括倒角。
3. 根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:S3中,采用ANSYS mesh将流体域三维模型进行网格划分,通过设置全局网格尺寸以及细化节流处网格,完成网格划分并生成网格文件。
4. 根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:S4中,流体介质为水、湍流模型为RNS K-e,求解器为Pressure-based,求解控制参数保持默认。
5.根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:S7中,数据处理工具为Excel。
6.根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:S2中,模型模块为DesignModeler模块。
7.根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:阀门的初始开度大于阀门的死区开度。
8. 根据权利要求1所述的一种阀门调节特性自动仿真方法,其特征在于:模型模块、求解模块和参数集模块均位于Fluid Flow模块中。
9.根据权利要求6所述的一种阀门调节特性自动仿真方法,其特征在于:S2中,阀门三维模型转化为中间格式*.x_t,然后再导入模型模块中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110560101.4A CN113255185B (zh) | 2021-05-21 | 2021-05-21 | 一种阀门调节特性自动仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110560101.4A CN113255185B (zh) | 2021-05-21 | 2021-05-21 | 一种阀门调节特性自动仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113255185A CN113255185A (zh) | 2021-08-13 |
CN113255185B true CN113255185B (zh) | 2023-05-02 |
Family
ID=77183764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110560101.4A Active CN113255185B (zh) | 2021-05-21 | 2021-05-21 | 一种阀门调节特性自动仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113255185B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114963014A (zh) * | 2022-06-10 | 2022-08-30 | 国家石油天然气管网集团有限公司 | 一种优化阀门开度降低波动压力的天然气输送减耗方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10114911B2 (en) * | 2010-05-24 | 2018-10-30 | Fujitsu Limited | Fluid structure interaction simulation method and apparatus, and computer-readable storage medium |
CN103729505B (zh) * | 2013-12-23 | 2017-01-18 | 苏州纽威阀门股份有限公司 | 一种基于cfd的阀门当量长度计算方法 |
US20160071004A1 (en) * | 2015-10-23 | 2016-03-10 | Sarkhoon and Qeshm LLC | Method and system for predictive maintenance of control valves |
CN110762268A (zh) * | 2018-07-09 | 2020-02-07 | 重庆川仪调节阀有限公司 | 一种调节阀内部尺寸设计方法 |
CN111767663A (zh) * | 2020-05-29 | 2020-10-13 | 江苏神通阀门股份有限公司 | 一种基于cfd仿真的阀门流量系数便捷计算方法 |
-
2021
- 2021-05-21 CN CN202110560101.4A patent/CN113255185B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113255185A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826130B (zh) | Bim建模在异形曲面结构形体设计建造中的使用方法 | |
CN105677964A (zh) | 一种基于cfd仿真和网格自适应的阀门流量系数计算方法 | |
CN108875213B (zh) | 一种反应堆堆芯热工水力多尺度分析方法 | |
CN111881524B (zh) | 一种阀门流量特性模拟实验方法 | |
CN109088776B (zh) | 一种基于超级计算机的并行cfd计算与通信重叠优化方法 | |
CN110762268A (zh) | 一种调节阀内部尺寸设计方法 | |
CN103135446B (zh) | 一种多轴数控机床刀具运动轨迹验证装置 | |
CN113255185B (zh) | 一种阀门调节特性自动仿真方法 | |
CN111062162B (zh) | 一种岩土材料精确本构模型的数值建模与应用方法 | |
CN105320806A (zh) | PML结合Excel-VBA处理数据汇总核电工艺管道材料的方法 | |
CN106650086A (zh) | 基于fluent软件的间接空冷***数值模拟平台 | |
CN115906718A (zh) | 一种旋转机械cfd*** | |
CN107506297A (zh) | 基于仿真平台的反应堆保护***功能测试方法和装置 | |
CN108763741A (zh) | 一种液压软管流固耦合数值预测方法 | |
CN112380795A (zh) | 一种基于预置模型库的试车台流场数值仿真*** | |
CN108197381B (zh) | 基于寻优空间形态分析的参数辨识方法 | |
CN112554864A (zh) | 一种计算产水气井单井控制储量的方法 | |
CN107808021B (zh) | 基于cfd的流体设备的阻力计算方法 | |
CN111680370B (zh) | 液压阀的设计方法及设计装置 | |
CN101561833B (zh) | 专用指令集处理器的设计方法 | |
CN105631135A (zh) | 基于abaqus基坑开挖快速建模分析***及其应用方法 | |
CN109684723A (zh) | 一种二维结构内部声学性能分析方法 | |
CN110728072B (zh) | 一种确定数字反应堆计算流体力学分析网格尺寸的方法 | |
CN102360400B (zh) | 自动模拟计算多方案城市大型缓流景观水体流态的设计方法 | |
CN113536640A (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 |