[发明专利]一种磁性液体密封环截面形状及密封耐压的计算方法

说明书

本发明公开了一种磁性液体密封环截面形状的计算方法，包括如下步骤：求得密封结构中的磁场分布；将密封间隙的空间沿径向等间距划分，拟合轴向坐标与磁场强度间的关系式；计算密封环截面自由界面曲线与极靴极齿表面的交点纵坐标；求解磁性液体的周向流速方程；计算关于周向速度的积分，并将其代入密封环截面边界的应力平衡方程中，将该方程的解代入磁性液体体积方程中，经循环迭代收敛后得到最终的截面形状和密封耐压。本发明还提供了一种磁性液体密封环密封耐压的计算方法。本发明通过迭代计算密封环截面形状，并在流速方程中引入离心力、磁化弛豫项，考虑周向速度引起的二次流，能够表征高转速工况下的磁性液体密封环截面形状和密封耐压。

技术领域

本发明涉及磁性液体密封技术领域，尤其涉及一种磁性液体密封环截面形状及密封耐压的计算方法。

背景技术

磁性液体是一种由包覆有表面活性剂的纳米量级铁磁性颗粒悬浮于液体介质中形成的胶体体系，密封是其应用最广泛的领域。在旋转密封中，磁性液体除受到热力学压力、粘性力、重力和磁体积力外，还受到离心力和表面张力等的影响，这些力的共同作用使得磁性液体密封环截面发生变形，当截面畸变达到一定程度后发生密封失效。

现有技术中，一般假设密封环截面边界曲线为磁力线，根据磁性液体静力学公式计算密封耐压。但在高速旋转密封中，由于离心力、二次流、磁化弛豫等效应的影响，磁性液体密封环的截面形状和耐压都将偏离以静力学为基础的密封理论的预测结果。因此，需要一种能够预测高转速工况下密封环截面形状和密封耐压的方法。

发明内容

有鉴于此，本发明实施例的目的在于提出一种磁性液体密封环截面形状及密封耐压的计算方法，通过迭代计算密封环截面形状，能够表征高转速工况下的密封环截面形状和密封耐压。

本发明一方面提出了一种磁性液体密封环截面形状的计算方法，包括以下步骤：

S1、基于密封结构计算密封结构中的磁场分布；

S2、沿密封结构的转轴径向为x方向，轴向为y方向，在密封间隙内建立二维坐标系，在x方向等距划分x_j，在每一个x_j处拟合y与磁场强度H的关系式y(H)；

S3、假定密封环截面的初始形状，分别用y₁和y₂表示第一和第二密封环截面自由界面的形状曲线，计算曲线y₂与极靴极齿表面的交点的横坐标ξ；

S4、在第一和第二密封环截面自由界面、转轴表面、极齿表面限定的计算区域内，计算磁性液体的周向流速方程以得到第一和第二密封环截面自由界面上的周向速度分布；

S5、基于周向速度分布，计算每一个x_j位置处的积分；

S6、将积分值代入密封环截面边界的应力平衡方程中，得到两自由界面上每一个x_j位置处的磁场表达式，磁场表达式包含未知量C_i；

S7、基于磁场表达式、y(H)以及磁性液体体积方程求解C_i，并基于C_i求解第一和第二密封环截面自由界面形状曲线；

S8、循环执行S3-7计算第一和第二密封环截面自由界面形状曲线，基于每一次的结果进行收敛判断，相应于当前形状曲线通过收敛判断，基于截面曲线输出截面形状结果。

在一些实施方式中，在S1中，使用有限元法求解密封结构中的磁场分布。

在一些实施方式中，在S2中，应用多项式拟合y坐标与磁场强度H间的关系式y(H)，并根据拟合标准差大于0.99的限值决定多项式的阶数，在转轴与极齿间的空间采用分段函数以获得高精度拟合关系。

在一些实施方式中，在S3中，假定密封环截面初始形状的自由界面与过极齿顶点的磁力线一致。