机械密封端面形貌的三维重建及其表征

机械密封端面形貌的三维重建及其表征
张浩;孙见君;马晨波;涂桥安
【摘要】In order to characterize the surface 3D topography and leakage path of mechanical seal end faces accurately and conveniently,a 3D reconstruction program for surface topography of mechan-ical seal end faces was developed and studied based on 2D sectional images.A new characterization pa-rameter for surface topography of mechanical seal end faces was proposed and based on voxel an algo-rithm was established to calculate the fractal dimensions and the porosity.The results show that it is effective to reflect the 3D box-counting dimension,porosity and other parameters of surface topogra-phy of mechanical seal end faces with 3D reconstruction based on sectional images.The surface porosi-ty of different layers has a law of change,so when the volume porosity is constant,the formation and the scale of leakage path of mechanical seal end faces will be dominated by the change of surface poros-ity distribution on different sectional images.The sampling number of layers have a great influence on the calculating results when it is fewer,however,the influence becomes smaller as the number rea-ches a certain value.%为了准确方便地描述机械密封端面的三维形貌及泄漏通道,基于密封环端面断层扫描二维图像,研究开发了机械密封端面形貌三维重建程序;提出了机械密封端面形貌的表征参数,并建立了基于体素的密封端面计盒维数和孔隙率算法.研究结果表明:基于断层图像的密封端面三维重建形貌能有效地反映密封端面的三维计盒维数和孔隙率等参数;层与层之间的面孔隙率具有一定的变化规律,在体孔隙率不变的情况下,改变各断层面孔隙率的分布,
可以控制动静环密封端面间泄漏通道的形成及泄漏通道的大小;当采样层数较少时,采样层数对计算结果影响较大,在采样层数达到一定数值后,采样层数对计算精度的影响越来越小.
【期刊名称】《中国机械工程》
【年(卷),期】2017(028)011
【总页数】6页(P1287-1291,1299)
【关键词】机械密封;表面形貌;三维重建;分形;孔隙率
【作者】张浩;孙见君;马晨波;涂桥安
【作者单位】南京林业大学机械电子工程学院,南京,210037;南京林业大学机械电子工程学院,南京,210037;南京林业大学机械电子工程学院,南京,210037;南京林业大学机械电子工程学院,南京,210037
【正文语种】中文
【中图分类】TH136
机械密封作为旋转设备中不可或缺的部件，在泵、压缩机及搅拌反应釜等工业装备上得到了广泛应用[1-3]。

机械密封对工业装备的可靠性、安全性有着重要影响，因其失效导致的泄漏不仅会引起资源浪费、环境污染，甚至还会引发火灾、人身伤亡事故。

为了揭示机械密封端面泄漏成因，人们提出了端面间隙锥角模型[4]、“波度”理论[5]和分形泄漏模型[6]等理论，将粗糙密封端面上的泄漏通道不同程度地简化为平均平行间隙、余弦间隙或锥度间隙。

事实上，具有一定粗糙度的机械密封动静环在配对时，除微凸体点接触外，接触点周围均为孔隙[7]；密封端面的泄漏可能性取决于各个孔隙组成的孔隙群能否形成贯穿密封端面的通道[8]。

如同
工程表面一样，机械密封端面形貌通常采用基于LONGUET[9]和NAYAK[10]的统计数据获得的2D粗糙度和3D粗糙度[11]来描述，但这些统计分析法存在分布的不稳定性、测量的多尺度性和仪器的依赖性等明显的缺点。

一些研究者探索采用轮廓线的二维分形维数[12]以及基于轮廓线二维分形维数D获得的三维形貌的分形
维数DB=D+1[13]来表征密封端面的形貌，克服了上述不足。

然而，机械密封规
律性的旋转运动会使动静环磨损表面周向和径向分形维数产生一定的差异，导致同一端面采用周向轮廓与采用径向轮廓的二维分形维数或基于此获得的三维分形维数不同，这大大降低了该表达方法的有效性。

为了准确和方便地描述密封界面的泄漏通道状态，本文利用密封端面的断层二维图像，通过体素绘制方法构建密封端面三维形貌；引入三维分形维数和孔隙率来表征密封端面几何形状复杂性和孔隙分布信息，研究建立基于体素的密封端面计盒维数和孔隙率计算方法，从而为机械密封的泄漏研究提供一种新的思路和方法。

三维重建包括面绘制和体绘制两类。

面绘制分为基于轮廓的表面重建和基于等值面的间接体素三维重建；体绘制包括体光线跟踪法和体单元投影法[14]。

LORENSEN等[15]在深入研究面绘制的算法后，提出了一种重建效果好、速度快、基于体素的等值面三维重建技术，并命名为“Marching Cubes”算法，即先将相邻断层间的8个网格点组成的单元定义为体素，然后确定一个表面阈值，计算每
一个体素内的梯度值，并与表面阈值进行比较判断，找出包含在等值面内的顶点，查找给定的索引表，确定体素的三角剖分形式，再利用差值的方法绘制构造体素的表面。

本文提出的密封端面三维重建流程如图1所示。

(1)测量仪器。

采用奥林巴斯OLS4100激光共聚焦显微镜对密封环端面进行断层
扫描，在扫描过程中需要保证激光共聚焦显微镜镜头的光轴垂直于样品表面，采样的高度为扫描范围内三维形貌最大高度，即从微凸体的最低谷到最高峰，图2为
扫描原理示意图。

测量范围：x、y向0～127 μm，z向0～1.6 μm；测量分辨率：
x、y向0.12 μm，z向0.01 μm。

(2)试样。

试样为浸渍酚醛树脂碳石墨环，端面尺寸为：内径53.0 mm，外径61.0 mm。

(3)测量结果。

在测量范围内对密封环三维形貌进行断层扫描，采样高度1.6 μm，共采集32幅图像，每层间距0.05 μm，图像大小为1024×1024像素，放大倍数为2041。

考虑到计算机的计算能力，将图像尺寸截取为512×512像素，图3为
部分断层图像。

图像预处理采用MATLAB完成。

首先将原始RGB图像转换为灰度图像；为了克
服采集过程中因扫描环境和扫描参数影响导致图像曝光过度或曝光不足产生的对比度不明显问题，对图像中每一个像素的灰度级进行标度变换，以扩大图像灰度范围，改善图像质量；然后进行直方图均衡化，将原始图像的灰度分布均匀化，在不影响整体对比度的情况下，增大目标与背景的反差，增强图像暗部信息；进而通过平滑滤波，消除图像在采集、传输和显示过程中受到干扰而产生的噪声，并保护图像的边缘信息；最后，将Otsu法计算的灰度值作为阈值对图像进行分割，分割后图像以只包含0和1的二值矩阵形式保存在计算机中，其中0代表密封端面的孔隙，1代表密封端面的固体骨架。

图4为预处理前和预处理后的图像。

将处理好的二维断层图像依照顺序读入MATLAB中，通过构造多维数组创建一个二值的三维体素矩阵M(体数据的长度l×宽度w×高度h为512×512×32)，计算
体数据集在显示平面的累计投影，并对碎片进行构造，定义图像颜色和光线，再设置图像的显示效果。

图5为浸渍酚醛树脂碳石墨环端面三维重建效果图。

图5所示的密封端面三维形貌不仅给出了研究表面的直观图像,还反映了表面几何
形状的复杂程度。

由于计盒数法数学表达简单，物理含义直观[16]，计算粗糙程度近似于磨削表面的三维分形维数最为准确[17]，所以本文采用计盒数法求解三维重建密封端面的分形维数；考虑到区别几何形状复杂程度相近的表面，以及为后续的
动静环接触端面间泄漏通道的研究提供便利，拟引入各断层面孔隙率和体孔隙率来描述密封端面的差异。

对于三维重建得到的密封端面二值化三维体素矩阵M，采用立方体覆盖法对体素
矩阵进行覆盖，从而得到三维分形维数。

具体步骤如下：①建立三维体素矩阵M；
②用边长为r的立方体依次去覆盖得到的三维体素矩阵，将包含目标体素(0表示
孔隙)的立方体计入盒子数Nr，遍历整个三维体素矩阵，得到某一盒子边长尺寸下的覆盖目标区域所需的盒子数；③根据计算公求出三维分形维数DB。

由于三维体素矩阵M是由二维断层图像叠加而成的，因此矩阵中的体素和二维图
像中的像素一样具有离散性，即当用边长为r的立方体覆盖时，其最小立方体为一个体素。

综合考虑计算机内存和计算精度，边长序列采用下式[19]构造：
利用MATLAB完成上述程序的设计，其流程如图6所示。

对密封端面上随机选取的4个采样区域A、B、C、D组成的三维体素矩阵进行计算，结果如表1所示，不同区域的拟合直线如图7所示。

对于A、B、C、D四组数据，拟合方程和相关数据如表2所示。

从表2中可以看出，4个区域的拟合方程都具有很强的线性相关度，4组分形维数之间的极差为0.0086，且与文献[17]中关于工程表面的分形维数计算结果相近。

上述数据表明，重建的密封端面具有较高的精度，能较准确地反映密封端面具有的分形特征，同时也表明一定精度范围内未磨损密封环的粗糙端面是各向同性的[20]。

对于图4b所示的二值图像，可以通过统计像素的方法计算目标像素(孔隙)占图像
总像素的比例，则密封端面二值化断层图像面孔隙率可利用下式计算：
推广到三维空间，对于图5所示的重建形貌，同样可利用统计体素的方法求取目
标体素(孔隙)占总体素数的比例，则三维密封端面体孔隙率可利用下式计算：
利用MATLAB开发的程序求取32幅二维断层扫描图像每一层图像的面孔隙率，
层数与面孔隙率的对应关系如图8所示。

可以看出，随着层数的增加，面孔隙率
逐渐增大，且在20层至30层内上升幅度明显增大，这说明密封端面形貌在这一层数区间内变化较大。

统计图5所示的三维体素矩阵中目标体素(孔隙)的个数V0和总的体素个数V，可以得到体孔隙率φV=V0/V=38.59%。

需要说明的是，体孔隙率与测量仪器的分辨率有关，不同分辨率仪器测量的密封环，其体孔隙率可能会存在较大的差异。

分形维数计算公式表明，在体孔隙率不变的情况下，改变各断层面孔隙率的分布，并不影响密封端面的三维分形维数，但可以通过改变各断层面孔隙率的分布来控制动静环密封端面间泄漏通道的形成以及泄漏通道的大小。

为了进一步探究采样层数与分形维数和孔隙率的关系，以相同的处理计算方法对区域D在不同采样层数下的分形维数和孔隙率进行计算，结果如图9所示。

可以看出，当采样层数较小时，无论是分形维数还是孔隙率的相对变化幅度都比较大；当层数超过16层之后，二者的变化幅度趋于平缓，即当采样层数大于16之后，增大层数并不能大幅提高计算结果的精度。

此外，较大的采样层数会成倍地增加计算机处理时间，降低计算效率。

因此，在综合考虑计算精度和计算效率的前提下，本文采样层数选取32是合理可行的。

(1)研究开发了基于二维断层扫描图像的密封端面形貌三维重建程序，通过MATLAB实现了密封端面形貌的三维重建；提出了基于体素的计盒维数和孔隙率表征方法，验证了密封环端面的各向同性特征。

(2)孔隙率随着层数的递增而变大，即密封端面底部形貌主要由固体骨架组成，顶部形貌主要由孔隙组成，同时靠近顶部区域的形貌变化幅度比底部区域的剧烈；在体孔隙率不变的情况下，通过改变各断层面孔隙率的分布，可以控制动静环密封端面间泄漏通道的形成及泄漏通道的大小。

(3)当采样层数小于16时，计算结果稳定性差，精度低；当采样层数大于16时，计算的分形维数逐步稳定在2.82左右，而体孔隙率逐步稳定在38.5%左右，此时
再过多地增加层数并不能大幅度提高计算精度。

(4)密封环端面的泄漏通道反映在本文的三维数据中即为相互连通的体素，因此本文提出的基于体素的描述方法能够为后续的动静环接触端面间泄漏通道的研究提供便利。