调节阀芯的静态不平衡力.docxVIP

调节阀芯的静态不平衡力 仪表误差是一种自动仪表执行器。随着过程控制的这一组件的自动化程度的提高，它已广泛应用于冶金、能源、化工、石油、纺织、建筑等行业。但是,由于阀门性能和质量问题造成重大事故,给工业生产的正常运行、人身安全和财产安全带来了不可估量的损失。据统计,调节阀在某些工况下产生振动是引起各种事故的主要原因,振动严重时甚至引起阀杆断裂,影响机组安全平稳地运行。在调节阀阀内流体压力高、压差大、流速快的情况下,由于不平衡力的作用,调节阀普遍存在着泄漏大、寿命短等重大问题,使用情况不理想。所谓调节阀的不平衡力就是指对直行程的阀芯所受到流体的轴向合力。不平衡力直接影响调节阀的行程位置与执行机构信号压力之间的关系。由于阀芯在中间位置时流体动压的影响,不平衡力难以用公式表示,因此,一般把调节阀全关时阀芯所受的静态不平衡力作为调节阀执行机构的设计依据。然而,阀芯匀速运动或在此基础上叠加振动情况下的动态不平衡力很少被考虑,这就有可能成为调节阀执行机构设计缺陷之一。 为了考察调节阀阀芯在调节过程中的动态不平衡力,本文针对直通单座调节阀,利用ANSYS软件CFD模块和任意拉格朗日-欧拉(Arbitrary Lagrangian-Eulerian, ALE)有限元法建立考虑调节阀阀芯移动边界流场模型,由此分别计算定流量和定压差下以匀速运动的阀芯在变开度时所受到的动态不平衡力。另外,根据阀芯系统运动方程,利用预测-校正流固耦合算法求得阀芯在变开度振动过程中的动态不平衡力。在此基础上,将变开度下动态不平衡力与定开度下静态不平衡力进行初步比较分析。 1 元模型的计算 建立调节阀流场有限元模型的目的是,通过流场有限元模型的计算获得作用在阀芯上的流体力,即不平衡力。本文采用ANSYS软件中的FLOTRAN CFD模块描述调节阀内单相粘性流体的三维流动。 1.1 anyss接口 调节阀流场的三维几何模型比较复杂,有许多复杂曲面,考虑到ANSYS建模功能的不足,几何模型没有用ANSYS本身的建模功能创建,而是用建模功能比较强大而且比较方便的SOLIDWORKS软件创建,如图1所示。将该几何模型另存为sat格式的文件,利用ANSYS的接口功能,把sat格式模型文件导入ANSYS。正确的流场模型应该是图1模型去掉内部阀芯后的结果。 1.2 阀芯底部流场变化梯度的网格 ANSYS软件提供了各种网格划分工具,如:自由网格划分、映射网格划分、扫略网格划分和过渡网格划分等。本文采用自由网格划分方法,其依据是:流体和阀芯接触的地方尤其是阀芯底部流场变化梯度较大,这些部位的网格要细化,以便更精确地计算和后处理;其它部位网格可以适当划分得粗一点,以便减小分析规模,提高计算效率。采用FLOTRAN CFD模块中的FLUID142单元对图1所示几何模型进行网格划分,流场网格划分结果如图2所示,即流场有限元模型,共有单元120 336个,节点24 016个,图3是流场从中间对称面剖开后的造型。 1.3 为减轻负荷和边界条件提供压力 对图2流场模型施加边界条件,具体的载荷及边界条件如下。 1 2进口定义的流量条件或压力条件 其值大小后面具体计算时给出,进口边界不能移动,即位移为0。 2 1出口仅定义压力条件，本工作采用0.1 MPa,出口边界不能动,即位移为0。 3 界接触面上的速度和位移 根据壁面无滑移的原理,假设所有与阀体边界接触面上的流动速度为0,位移为0,即除了进口端面、出口端面和阀芯周围的面以外,其余所有边界的流体速度为0,位移为0。 4 2阀芯和流场之间的界面移动壁面边界条件 流场在阀芯周围的固体接触面上的流体速度应该和阀芯的运动速度一致。 2 不同阀芯在不同工况下的流场分析 两种不同介质共同界面的运动,即结构与流体接触界面的移动是流固耦合问题中的一个特殊问题。调节阀变开度下阀芯上动态不平衡力计算涉及到具有移动边界的流固耦合场分析。针对不同情况下的调节阀流场计算,要用不同的方法。对于计算调节阀在一个稳态开度下的流体力相对比较容易,这个分析模型可以归纳到固定壁面类型分析。对于计算调节阀开度调整过程中的流体力相对就比较困难,而且计算量远远大于稳态开度。如果调节阀阀芯的运动规律已经给出,这个应该归纳到移动壁面类型分析;如果阀芯的运动规律没有给出,只是给出控制力的模型,或者只给出了阀芯受力情况,这个就应该归纳到流固耦合类型分析。本文讨论后一种情况下的流场分析和计算模型。 2.1 移动边界问题 当阀芯按给定速度作匀速运动时,阀芯与流体之间的接触边界也作匀速运动。这种移动边界问题是流固耦合问题中的一个特例。由于移动边界的运动规律已经给定,即可以直接给出移动边界任意时刻的位移和速度,因此,可以归结为移动壁面类型,用ALE有限元法进行分析计算不平衡力。 2.2 未赋予阀芯运动方程求解