高压计量泵单向阀密封结构的有限元分析.docxVIP

高压计量泵单向阀密封结构的有限元分析 随着生产的发展和生活质量的提高，计量泵的应用范围也越来越大。它延伸到经济生产的各个领域。同时, 科学技术的发展和生产工艺控制过程要求的提高, 对计量泵的技术含量和技术要求日渐苛刻, 迫使国内计量泵的生产开发速度不断加快。计量泵的单向阀密封技术关系到计量泵的运行效率, 因此计量泵单向阀密封技术及失效机理的研究越来越引起人们的重视。可是, 目前对它的力学性能分析和计算涉及尚少。为此, 笔者采用接触有限元法, 力求在这方面做一些基础性工作。 问题的提出 图1所示为某计量泵 (工作压力250 MPa) 的单向阀结构示意图。 由其结构可知, 通过阀座之间的密封面及阀座与三通之间的密封面把流体限制在单向阀内流动, 当它们之间密封不严时, 就会导致单向阀的泄漏。为了保证密封, 需将接头螺纹上紧, 即加大螺纹的预紧力。但如果预紧力过大, 超过材料的屈服极限, 阀座发生塑性变形, 也会造成单向阀破坏, 从而使密封失效。因此, 需对高压计量泵单向阀密封结构的破坏和泄漏进行接触分析, 并确定高压接头的最大和最小上扣扭矩。 有限元计算与分析 1. 接触单元的生成 在实际应用过程中, 上阀座底部与下阀座接触的地方最易破坏, 所以力学模型可这样简化:取整个上阀座与部分下阀座, 这样既维持两者之间的接触不变, 又简化结构, 在计算结果不失真的情况下可减少计算时间。 考虑到阀座的结构和所受的载荷都是轴对称的, 笔者采用轴对称方法进行分析。由于采用平面模型模拟三维实体, 所以选用了plane42单元, 单元属性定义为轴对称。阀座的材料为0Cr17Ni12Mo2, 根据手册查到:弹性模量为2.05×105MPa, 屈服极限取为205 MPa, 强度极限取为520 MPa, 延伸率取为60%, 泊松比0.25。 考虑到三角形自由网格易于细化, 所以采用三角形网格对模型进行自由划分, 然后细化接触部分的网格, 最后共生成10 590个plane42单元。取下阀座的顶面作为目标面, 上阀座的底面作为接触面, 利用接触向导建立面-面接触, 最后生成targe169单元136个, conta172单元140个。接触分析的计算模型如图2所示。从图中可以看出, 接触部分的网格非常细密, 这样有利于提高计算精度。 2. 破坏元分析模式的计算和分析 (1) 大变形瞬态计算 上阀座顶部和螺纹接头接触, 考虑到阀座比螺纹接头材料硬度大, 所以可假设上阀座顶部是刚性位移, 在ANSYS中对上阀座顶部节点的轴向位移进行耦合。 由于下阀座取得足够长, 可以近似认为模型的底边受接触的影响很小, 已经没有轴向位移, 所以对底边的节点施加轴向约束。 对上阀座顶部的耦合集施加轴向载荷来模拟高压接头带来的压力, 预先不知道阀座被压坏时到底对应多大的轴向力, 所以此时的轴向力应该进行试算, 边界条件和载荷施加完毕后模型如图2。 由于是试算压力, 故可进行瞬态计算, 将载荷取作渐变载荷, 多设时间步, 将时间步长取小, 这样在后处理时可以观察不同载荷时的应力状况。考虑这种接触牵扯到塑性变形, 所以选择大变形瞬态计算。对于接触模型, 选用改进的拉格朗日方法, 摩擦因数取作0.2。 (2) 上扣扭矩的应力分析 笔者采用最大值为12 kN的瞬态渐变载荷进行计算, 计算结果分析如下。 图3为轴向载荷10 kN时的等效应力云图, 所用理论为第四强度理论。从图中可以看出, 施加载荷处及上、下阀座接触处应力集中, 其中下阀座上平面的内缘受力最大, 为203.582 MPa, 并未屈服。 图4为轴向载荷10 kN时的接触应力云图。从图中可以看出, 接触面的接触应力比较均匀。 图5为轴向载荷10 kN时接触面上各个节点的接触应力曲线。从图中可知, 第2个节点处的接触应力值最大, 为312.555 MPa。 图6为下阀座上平面的内缘处第1节点的等效应力随轴向载荷变化的曲线图。图中横坐标为载荷占所加总载荷12 kN的比率, 从图中可以看出, 该曲线与材料的特性曲线基本一致。该节点大约在载荷比率0.865时进入屈服, 此时该节点的等效应力为205.188 MPa。此时所对应的轴向载荷:F=10.38 kN。 图7为不同载荷时接触面上各个节点的接触应力的曲线图。从图中可以看出, 随着载荷的增大, 各个节点的接触应力比较均匀地增大, 没有接触应力集中。 通过以上分析, 可得最大上扣扭矩时所对应的轴向载荷为10.38 kN。确定出阀座屈服时的轴向载荷后, 便可通过螺纹的扭矩公式算出高压接头的最大上扣扭矩。由文献可得, 扭紧螺纹所需力矩T为螺纹摩擦力矩T1和支撑面摩擦力矩T2(T′2) 之和, 计算螺纹扭紧力矩的计算公式为 T=T1+T2′=F2d2tan(λ+ρυ)+12Ffldm=F2[d2