低温阀门闸板应场的数值计算及分析.docVIP

低温阀门闸板应力场的数值计算及分析 低温阀门闸板应力场的数值计算及分析李秀峰陈宗华（天华化工机械及自动化研究设计院）值模拟，并绘制了应力等值线图、位移图和应力变化曲线图，探讨了应力与变形的规律，为进一步完善低温闸板的设计提供了理论依据。 　　低温闸阀作为特殊用途的阀门，工作温度一般为 196　40 ，最高使用压力达10M Pa，且通常用于易燃、易爆及渗透性强等介质场合，工作条件要比普通阀门恶劣得多，闸板作为阀门密封件，其受力引起的变形将直接影响阀门的工作性能及可靠性。虽然其温度变化范围不大（完全处于低温环境中），但低温时流体分子对密封面楔入能力的增强，使阀门的密封十分困难，因而必须要提高阀头和阀座之间的密封压力，这样整个闸板的受力状况也变得更为复杂，因此对其进行受力分析是十分必要的。但是，低温阀门闸板在工作过程中的受力状况一般是没有资料可查，并难以进行实际测量的。笔者认为利用数值模拟方法研究应力变化规律是行之有效的手段。本文建立了应力场的数学模型，采用有限元法进行求解并编制了相应的计算程序。通过对闸板应力场的分析研究，揭示了闸板内部各向应力变化过程的基本规律，指出了结构的薄弱环节和改进方向，对设计计算具有一定的指导意义。 　　1模型的建立1. 1物理模型的简化低温阀门闸板与普通闸板在结构上没有本质的区别。闸板T形槽的形状不是很规则，其主要作用是对闸板施加闸紧力，在一般情况下它的结构形式不会对闸板密封面受力产生大的影响。所以，为了简化实物模型，在进行分析时可以把T形槽部位略去，而代之以一定的边界应力（图1a）。这样，楔式闸板就简化为一个大致的楔形圆台的形状。此时，除了密封面凸台有大约5的倾角外，其平板部分为均匀等厚，考虑到此倾角很小，故可以略去其影响，而将密封面凸台的高度取其平均值。这样，闸板从几何形状上就简化为轴对称图形（图1b）1. 2力学模型的建立闸板受力可以分解为水平和垂直两个方向（图2）。计算结果显示，在水平方向上，由于低温条件下难以密封的特点， N将比普通阀门大得多，实际上它们也是影响密封效果的主要因素，q为面均匀分布载荷垂直力（主要是摩擦力）引起的应力值约占总应力值的10，通过电测试验验证得知，摩擦力的影响不大于0. 1.为了简化计算，在系列计算时忽略了摩擦力。 　　化工机械李秀峰，男， 1978年8月生，硕士研究生。甘肃省兰州市， 730060. 　　2有限元法建立数学模型由于低温条件下闸板密封比压过大，在凹槽附近会产生明显的应力集中，为了正确反映此处应力，必须把该处的网格划分进行加密处理。但是网格加密以后会加大计算量，耗费过多的计算时间，而且如果相邻单元的尺寸相差过于悬殊，可能还会引起很大的计算误差。为了解决上述问题，本文采用了子模型技术，将计算分成两次进行：第1次计算时，把所有网格划分得大致相同疏密，主要的目的在于计算出别处的应力，并计算出此部分边缘一线上各结点的位移第2次计算时，只以凹槽部分为计算对象，把附近的网格划分得充分细密，而将第1次计算所得的边缘一线上各结点的位移作为已知量输入，即可将凹槽附近的应力集中区域计算得充分精密。 　　本文采用三结点三角形环状单元，这种单元适应性好、便于分析且具有很强的边界拟合能力，是分析轴对称体最常用的单元。网格基本大小为1mm，总数为1 676个，阀体内腔施加2. 5M Pa的压力。以对称轴作为z轴，所有应力、应变和位移都与方向无关，只是r和z的函数，因此原来的空间三维问题就可简化为以r、z为自变量的二维问题。 　　2. 1构造插值函数选择结点位移作为基本未知量，利用结点位移分量进行函数插值，在局部坐标系中构造单元的位移函数，单元中任意一点位移可以用结点的位移表示：式（ 1）可简写为：式中A三角形面积插值函数单元的结点位移值。 　　2. 2求单元应变利用平面问题的几何方程可以求得以结点位移表示的单元应变：式（ 4）简写为：几何矩阵。 　　2. 3求单元应力利用平面问题的物理方程，可以得出以结点位移表示的单元应力：化工机械弹性矩阵应力矩阵弹性模量泊松比。 　　2. 4建立结点平衡方程由线弹性力学的泛函变分原理，可得到系统总位能的离散形式，它是弹性体变形位能和外力位能之和由总位能的变分= 0，整理上式得出：组合各单元刚度矩阵形成总体刚度矩阵，得出如下公式：式中K总刚度矩阵结点位移列阵载荷列阵。 　　2. 5应力光滑处理2. 5. 1数值计算误差产生的原因只有在位移变分完全任意的情况下（单元尺寸趋于零时），平衡方程、力的边界条件以及单元交界面上内力的连续条件形成泛函的欧拉方程才能精确地满足。且由于以位移作为基本未知量，应变矩阵B是插值函数N对坐标进行求导后得到的矩阵，求导一次插值，多项式的次数就降低一次。所以通过导数运算得到的应变和应力可能具有较大的误差。 　　根据最小位能原