高压先导轴流式天然气减压阀性能优化研究

本资料包含pdf文件1个，下载需要1积分

21世纪是世界天然气工业高速发展的时期，也是我国天然气工业兴起并走向成熟的时期。天然气减压阀是燃气输配系统中的关键设备，是安全供气的重要保障 j。在实际运行中，管 网会产生各种各样 的扰动，这些扰动会影响减压阀的稳定性，造成其出口压力的波动2 J。目前减压阀设计都是以启闭临界点静态设计为基础，无法描述减压阀动态特性。加强天然气减压阀动态特性的研究，对于提高天然气输配系统的稳定性和节省投资，有着重要的现实意义。

我项目组与国内某阀门企业合作展开了系列天然气减压阀研究，采用先导轴流式结构，具有体积孝流量大、流阻小等优点。减压阀参数为 PN100、DN50，其结构如图 1所示，实际工况最高进 口压力为 9.65MPa，最大流量可达 123330 Nm /h。针对天然气大型输配管网及门站对减压阀稳定性的高要求，建立了机械系统运动方程、气体状态和流动方程以及各腔压力控制方程，运用 Simulink工具箱进行动态特性仿真，绘制不同进口压力、弹簧刚度以及泄流孔直径下天然气高压先导式减压阀的动态特性响应曲线。并根据不同设计变量，建立以输出压力为目标函数的优化数学模型，基于遗传算法理论得出最佳适应度和最佳均值轨迹，实现了减压阀的性能优化。

图 1 减压阀结构 图1 减压阀动态数学模型的建立1.1 减压阀结构及假设由减压阀结构画出其原理示意图如图2所示。根据气体在减压阀内的流动及状态变化规律，为简化分析对先导式减压阀模型作如下假设 ：(1)流通介质为理想气体，具有恒定比热容；(2)系统内高流速气体不与外界换热，各腔内气体发生等熵变化，且同-腔体内气体分布均匀，物性收稿日期：2012.10.31基金项 目：校优秀青年教师培养计划 (Q201002)；校企横向合作项目(H1205cc005)作者简介：李树勋(1971-)，男，陕西榆林人 ，副教授 ，工学硕士，主要从事特种阀门技术、过程测试技术的科研工作。

2013年第4期 液压与气动 33相同；(3)忽略气体流动时的沿程阻力。

导阀信号腔 导阀载荷腔号腔压腔图 2 减压阀原理示葸 图1.2 机械 系统运动方程主阀机械系统运动方程为：m- kl(z ) 面dx(pz-P1) (1)式中，m 为主阀机械系统的有效质量； 为主阀薄膜阀杆连接件的位移； 为时间；i 为主阀弹簧刚度；Z 为主阀阀芯全闭时弹簧压缩量；4 为主阀薄膜有效面积； 为主阀薄膜阀杆连接件运动时的阻力系数；P。

为主阀信号腔气体绝对压力；P：为主阀载荷腔气体绝对压力。

导阀机械系统运动方程为：m ： dt2 k2(f2 (p )以以(2)式中，m 为导阀机械系统的有效质量；y为导阀薄膜阀杆连接件的位移；g为重力加速度；j 为导阀弹簧刚度；f：为导阀阀芯全闭时弹簧压缩量； 为导阀薄膜有效面积；厂2为导阀薄膜阀杆连接件运动时的阻力系数；P 为导阀信号腔气体绝对压力；P 为导阀载荷腔气体绝对压力。

1.3 气体状 态方程由模型假设，忽略流体介质经过先导式减压阀时与周围环境的换热，则在各个腔体内进行气体绝热过程 。

1)主阀载荷腔气体状态方程为：dp， ( Ad dxAd ) (3) d l ， -式中，q 为主阀低压腔流入主阀载荷腔的流量；K为气体等熵指数；V1为主阀载荷腔体积。

信号腔气体状态方程为：dp2 ( Ad dxAd ) (4) d -， 2 f， 、。

式中，9 为主阀低压腔流入主阀载荷腔的流量；v2为主阀载荷腔体积。

低压腔气体状态方程为：： ， (5)d -t 93 L)式中，q 为主阀低压腔输人体积净流量；P 为主阀低压腔绝对压力；v3为主阀低压腔体积。

2)导阀信号腔气体状态方程为： (g4Ad： dtAd ) (6) d - ，、”。 ，式中，g 为系统流入导阀信号腔的流量；v4为导阀信号腔体积。

载荷腔气体状态方程为：dP5 ( Ad dtAd ) (7) d ，5 ， 、式中，g 为导阀载荷腔输入体积净流量； 为导阀载荷腔体积。

1.4 各腔流量方程整个系统流量为：g q1q2g3q4q5 (8)式中，g为经过系统流人减压阀的总流量。

1)主阀载荷腔气体流量方程为：91 1A1 3-P1 (9)式中， 。为载荷腔凶的流通面积； 。为载荷腔凶处流量计算系数。

信号腔气体流量方程为：q2(2A2/p4-P2 (10)式中，A：为信号腔凶的流通面积； 为信号腔凶处流量计算系数。

低压腔气体流量方程为：q3q- 1A1，/p3-P1- 2A2 4-P2- A √ -P - A /fp5-P。 (11)式中，A，为低压腔凶的流通面积； ，为低压腔小孑L处流量计算系数； 为导阀信号腔凶的流通面积；为导阀载荷腔凶处流量计算系数；A 为导阀信号腔凶的流通面积； 为导阀载荷腔凶处流量计算系数。

34 液压与气动 2013年第4期2)导阀信号腔气体流量方程为：q4g"4A4/p-P4 (12)式中，P为减压阀入口绝对压力。

载荷腔气体流量方程为：q5(sA5 5-Po (13)式中，P 为标准大气压力。

1.5 各腔连续性方程联立各个腔体气体状态方程和流量方程，可得减压阀系统连续性微分方程。

1)主阀低压腔连续性微分方程为：-dP3： fq- 4 - zAz 1(14) dt'/3- (4A -载荷腔连续性微分方程为：dpl ( A -Ad dx)(15)信号腔连续性微分方程为：ap2 ( Ad。 dx)(16)2)导阀信号腔连续性微分方程为：dp4 ( A Ad： ) (17)载荷腔连续性微分方程为：ap5 ( A -Ad： dt)(18)由式(1)、(2)、(14)～(18)构成的方程组为所求先导式减压阀的动态模型，其中 、Y，P 、P：、P。、P 、P 为连续状态变量，其余参数如减压阀内部结构参数，进口压力和 口压力等为已知量。

2 动态特性仿真利用 MATLAB仿真拈 Simulink建立系统方框图，选取不同的结构参数并进行仿真得到减压阀系统动态特性曲线 。]▲行仿真前设定的工况是，0 S之前先导式减压阀以入 口压力4 MPa稳定运行，流量为36100 Nm /h，0 s时先导式减压阀后管路系统负载突然增大，流量增至43710 Nm /h。

如图3所示，该图反映的是不同先导阀弹簧刚度k，取值下对主阀薄膜阀杆连接件的位移响应曲线。其中3条曲线分别代表了弹簧刚度 k 取值为 165 N/mm、175 N/mm、185 N/mm时薄膜阀杆连接件位移的响应曲线。从图中可以看出，当存在干扰时，随着弹簧刚度k 的增大，薄膜阀杆连接件位移波动程度减小，超调量逐渐减小，阀后压力在较短时间内达到稳定状态，有效地减小了阀后管路的振动程度，使减压阀获得了较好的稳定性。

3.02.52.01.5越 1 o时间∥图 3 不同先导阀弹簧刚度对阀芯位移的影响如图4所示，该图反映的是泄流孔直径 d与主阀薄膜阀杆连接件位移的响应曲线，其中3条曲线分别代表了泄流孔直径 d为4 mm、5 mm、6 mm时薄膜阀杆连接件位移的响应曲线，从图中可以看出当泄流孔直径取 5 mm时，既能保持减压阀对与系统压力波动的敏感性，又能保持较小的阀芯位移超调量，使减压阀获得稳定的动态性能。

9.o8.0宕7.o6.o簿时间f/0图4 不同泄流孔直径对阀芯位移的影响如图5所示，该图反映的是不同主阀信号腔体积与减压阀阀后输出压力的响应曲线。其中3条曲线分别代表了信号腔体积 取值为0.00014 m 、0.00015m 、0.00016 m 时阀后输出压力的响应曲线。从图中可以看出随着信号腔体积 的增大，阀后输出压力波动程度减小，超调量逐渐减小，阀后压力在较短时问内达到稳定状态，有效地减小了阀后管路的振动程度，使减压阀获得稳定动态性能。

善出丑箨O时间∥s图5 不同信号腔体积对阀后压力的影响2013年第4期 液压与气动 353 性能优化设计函数优化问题是遗传算法的经典应用领域，也是对遗传算法进行性能评价的常用算例 J。以阀后输出压力P 为目标函数， ，为阀后压力稳态值，这里选择弹簧刚度与主阀信号腔体积为优化对象，并为其选择-个适度的优化范围，建立天然气高压先导轴流式减压阀的优化数学模型为：min 二 d (19).to wp3 ，S.t. 1× 10- ≤ V ≤ 2 ×10-150 ≤ k2≤ 250在 MATLAB工具箱 gatool中输入编好的适应度函数名 myfun，以预设的规则对其 目标函数进行优化。

变量个数设置为2个，分别为弹簧刚度与信号腔体积，在约 束 栏 中输 入 边 界 条 件 [150 0.0002]，[2500.0005]。在图形栏中选择 Best fitness和Best indi。

vidual选项，以便以图形的方式判断目标函数是否收敛，右边的选项栏填人初始值并改变变异因子。遗传算法工具箱中的 Population type(种群类型)选择 doublevector选项，Population size(种群尺度)选择20，其他选项选择 default(默认设置)，所得结果如图6所示，随着代数的增加，最佳适应度值趋于平稳，到50代时几乎重合，说明找到了最佳个体，满足目标函数的最小值要求。

60 ·最佳适应度值籁30. ·平均适应度值 熠阔。0 -30-40-50 60 10 20 当前最佳个体10。00[ . 韫 1.图6 最佳适应度值与最佳个体X：[1.941139e-004 1.7882e002]fval：67.283exitmessage：[1 x86 char]output：[1xl struct]根据遗传优化的结果，将减压阀主阀信号腔体积由1.44×10 I4 m。增大为 1.94×10 I4 m ，将先导阀的弹簧刚度 k2由 1.67×10 N/ram提高为 1.79×10N/mm。改进的参数带入 Simulink模型当中得到优化后的减压阀动态特性曲线。由图7可以看出，在人 口压力 4 MPa的工况下，先导式减压阀内部结构参数经过优化后 ，阀后输 出压力最大超调量 由之前的 1.7MPa减小为 0.5 MPa，系统动态响应时间从 0.65 S缩短为 0.42 S，减压阀的动态特性得到了优化。

签鏖图 7 减压 阀结构参数 改进前后动态特性对 比4 结论根据天然气先导式减压阀工作原理，结合气体动力学基础知识，分析了不同的弹簧刚度、泄流孔直径以及信号腔体积等结构参数对减压阀的动态特性的影响。应用遗传算法对天然气先导式减压阀结构参数进行优化，结果表明：在-定的工况下，适当增大主阀信号腔体积和先导阀的弹簧刚度能够提高减压阀稳定性与灵敏性，实现了动态性能的优化，为今后进-步的研究工作提供了参考。