斯特林制冷机间隙密封泄漏量的影响因素权重分析

徐安波, 陈晓屏, 李昊岚, 孙皓

徐安波, 陈晓屏, 李昊岚, 孙皓. 斯特林制冷机间隙密封泄漏量的影响因素权重分析[J]. 红外技术, 2021, 43(1): 68-72.
引用本文: 徐安波, 陈晓屏, 李昊岚, 孙皓. 斯特林制冷机间隙密封泄漏量的影响因素权重分析[J]. 红外技术, 2021, 43(1): 68-72.
XU Anbo, CHEN Xiaoping, LI Haolan, SUN Hao. Weight Analysis of Influencing Factors of Clearance Seal Leakage in Stirling Cryocoolers[J]. Infrared Technology , 2021, 43(1): 68-72.
Citation: XU Anbo, CHEN Xiaoping, LI Haolan, SUN Hao. Weight Analysis of Influencing Factors of Clearance Seal Leakage in Stirling Cryocoolers[J]. Infrared Technology , 2021, 43(1): 68-72.

斯特林制冷机间隙密封泄漏量的影响因素权重分析

详细信息
    作者简介:

    徐安波(1992-),男,云南曲靖人,硕士研究生,研究方向:小型低温制冷机。E-mail:1006514396@qq.com

    通讯作者:

    孙皓(1979-),男,云南通海人,研高工,博士,主要从事小型低温制冷机研究。E-mail:652700@sina.com

  • 中图分类号: TH136

Weight Analysis of Influencing Factors of Clearance Seal Leakage in Stirling Cryocoolers

  • 摘要: 为分析斯特林制冷机间隙密封泄漏量影响因素权重,结合同心环形缝隙流动理论,确定影响因素包括间隙厚度、密封长度、间隙两端压差、转速、间隙内径以及活塞冲程。通过正交试验搭建32组数学模型,基于计算流体CFD数值模拟的方法,研究分析不同密封间隙的泄漏情况。利用方差分析对仿真结果进行处理,结果表明间隙厚度对斯特林制冷机的泄漏量影响最大;压差、密封长度、活塞冲程以及转速对泄漏量影响依次减小;间隙内径对泄漏量影响最小。这对斯特林制冷机间隙密封减小泄漏量的设计提供了参考。
    Abstract: In this study, the weight analysis of the influencing factors of clearance seal leakage of Stirling cryocoolers was conducted based on the concentric annular gap flow theory. Influencing factors, including sealing gap height, length, pressure difference, motor speed, clearance diameter, and piston stroke, were determined. Through an orthogonal test, 32 sets of mathematical models were established, and the leakages of different seal clearances were analyzed based on the numerical simulation method of computational fluid dynamics. The simulation results were processed using variance analysis, and indicated that the clearance thickness exhibited the most significant influence on the leakage of a Stirling refrigerator. Pressure difference, sealing gap length, piston stroke, and motor speed exhibited a lower influence, where as clearance diameter demonstrated a minimal impact on leakage. This study provides a reference for the design of a Stirling cryocooler clearance seal to reduce leakage.
  • 红外成像具有非接触直接测量输变电设备状态的特征,无需停电,在电力领域具有广泛的应用[1-3]。输变电设备在野外严峻环境中长期运行,由于材料老化、污秽闪络、机械受损等因素,常伴随有局部放电、温度增高等现象。绝缘子作为输电线路重要部件,起到电气隔离和机械支撑的作用,当红外图像显示的温度较高,表明其有异常缺陷,这是因为场强不均匀发生局部放电,严重时可能导致线路故障,甚至造成停电故障。文献[4-5]开展了变压器高压套管的红外诊断研究,通过热像特征谱图分析了故障原因。红外图像可应用于交流滤波器的故障分析中,提取其跳闸的典型故障特征[6-8]

    输变电设备红外图谱特征提取主要为图像处理方法,诸如纹理、色彩、边框等图像特征[9-10]。随着无人机航拍技术的发展,输变电设备红外图谱呈现爆发式增长,深度学习提供了一种良好的手段[11-13],采用卷积神经网络(Convolutional Neural Networks,CNN)的方法,对海量图片进行训练学习,提取特征进行测试和验证。CNN通过卷积层计算,输入较深层的特征图,对于小目标绝缘子、套管,权重值相对较少,无法实现小目标的有效提取[14-15]。针对这一缺点,本文对Faster R-CNN方法进行改进,提高绝缘子红外图谱诊断的精度。

    区别于普通的CNN,Faster R-CNN增加了一个区域提取网络(Region Proposal Network,RPN),即图 1中的候选区域,摒弃传统的滑动窗口,可在GPU直接运行计算,极大地加快了计算速度。RPN判断每个像素点对应的多个不同尺度和宽高比的锚框是否为前景目标的二分类,形成候选区域。

    图  1  Faster R-CNN的算法流程
    Figure  1.  Algorithm flow of Faster R-CNN

    Faster R-CNN一般采用随机梯度下降法(Stochastic Gradient Descent,SGD)训练神经网络,见式(1):

    $$ h(x) = \sum\limits_{i = 0}^n {{w_i}} {x_i} = {\mathit{\boldsymbol{W}}^{\rm{T}}}\mathit{\boldsymbol{X}} $$ (1)

    式中:X为输入;W为权重;wixi分别表示第i个权重和输入;h(x)为对应的输出。

    损失函数S(W)基于平方误差实现,见式(2):

    $$ S(\mathit{\boldsymbol{W}}) = \frac{1}{2}\sum\limits_{i = 0}^m {{{({h_W}({x_i}) - {y_i})}^2}} $$ (2)

    式中:yi为真实输出。

    W的更新函数见式(3):

    $$ {W_j} = {W_j} - \alpha \frac{\partial }{{\partial {W_j}}}S(\mathit{\boldsymbol{W}}) $$ (3)

    式中:α为学习率,可设置步长。W通过梯度下降法进行求解,首先正向计算样本输出值,接着根据反向传递的误差迭代计算,常用在CNN训练中。

    RPN的选择本质是通过平移或者尺度变化的方法将合适的区域提取网络R变成$\hat C$,从而接近实际的候选框C

    $$ f({R_x}, {R_y}, {R_w}, {R_h}) = \left( {{{\hat C}_x}, {{\hat C}_y}, {{\hat C}_w}, {{\hat C}_h}} \right) \approx \left( {{C_x}, {C_y}, {C_w}, {C_h}} \right) $$ (4)

    式中:(x, y)、(w, h)分别为矩形区域的中心坐标和宽、高。

    令${t_*}$为矩形区域的平移和缩放量,则有:

    $$ \left\{ {\begin{array}{*{20}{l}} {{t_x} = ({C_x} - {R_x})/{R_{\rm{w}}}}\\ {{t_y} = ({C_y} - {R_y})/{R_{\rm{h}}}}\\ {{t_{\rm{w}}} = \lg ({C_{\rm{w}}}/{R_{\rm{w}}})}\\ {{t_{\rm{h}}} = \lg ({C_{\rm{h}}}/{R_{\rm{h}}})} \end{array}} \right. $$ (5)

    式中:txty为矩形区域的中心坐标平移量;twth分别为矩形区域的宽、高的缩放量。

    预测值计算过程为:

    $$ {d_*}(R) = w_*^{\rm{T}}\phi (R) $$ (6)

    式中:ϕ是最后一次卷积计算。

    损失函数的目标值计算见式(7),通过调整平移和缩放的尺度,确定最终的候选区域[16]

    $$ {L_{{\rm{oss}}}} = \sum\limits_{i = 1}^N {{{(t_*^i - w_*^{\rm{T}}\phi ({C^i}))}^2}} $$ (7)

    为了增强Faster R-CNN的小目标特征提取能力,引入压缩激励结构,即压缩和激励两大操作,设图像的特征参数设置为(H, W, K),分别表示为长、宽和通道数。

    压缩操作Fsq(·)基于各个通道实现特征图空间信息的压缩,见式(8):

    $$ {h_c} = {F_{{\rm{sq}}}}\left( {{k_c}} \right) = \frac{1}{{H \times W}}\sum\limits_{i = 1}^H {\sum\limits_{j = 1}^W {{k_c}(i, j)} } $$ (8)

    式中:kc表示第c个通道;hc表示经过压缩后输出向量h的第c个元素。

    激励操作分为激励Fex(·)和校准Fscale(·)两个过程,分别见式(9)和式(10):

    $$ s = {F_{{\rm{ex}}}}(h, w) = \sigma (g(z, w)) = \sigma ({w_2}\delta ({w_1}z)) $$ (9)

    式中:σ是sigmoid激活函数;w1为$\frac{c}{r} \times \vec C$的实数矩阵,表示通道的缩减,r为缩减因子;δ是ReLU激活函数,w2为$\vec C \times \frac{c}{r}$的实数矩阵,表示通道的恢复。

    $$ {\tilde h_c} = {F_{{\rm{scale}}}}({h_c}, {s_c}) = {s_c}.{h_c} $$ (10)

    式中:sc表示激活向量s的第c个元素;${\tilde h_c}$表示校准后的对应元素。

    本文基于改进的Faster R-CNN方法,对平台的环境搭建要求较高,具体配置见表 1。操作系统为开源Linux,数据库为MySQL;硬件配置较高,CPU采用Intel高端系列,内存和硬盘容量均较大,保证大量数据的高效运算。框架采用2018年初公开的目标检测平台Detectron,包含最具代表性的目标检测、图像分割、关键点检测算法。

    表  1  软硬件配置
    Table  1.  Hardware and software configuration
    Name Model
    Operating system Ubuntu 16.04.1
    Database mysql 5.5.20
    CPU Intel Xeon Silver 4114T 12C
    GPU NVIDIA GTX1080Ti
    Memory 32 G
    Hard disk 1 T
    Frame Detectron
    下载: 导出CSV 
    | 显示表格

    图像数据来源于多条输电线路无人机拍摄的大量绝缘子照片。在神经网络的训练过程中,对正负样本的判定见图 2,主要基于锚框映射图与真实目标框的交并比(Intersection over Union,IoU)来进行计算分类。首先对RPN形成的锚框进行排序筛选形成锚框序列,接着利用边框回归参数向量修正锚框的位置形成候选区域集合,然后计算所有感兴趣区域(Region of Interest,RoI)与真实目标框的IoU,求最大值,并判断其是否大于0.5,若满足,则为正样本,否则为负样本。

    图  2  正负样本判定
    Figure  2.  Positive and negative sample decision

    完成正负样本的判定后,为了使样本的采样尽量均衡,保证双方的训练集和验证集数量一致,同时采用迁移学习的方法,经过相关修正与补偿,扩充样本总量至2375,样本信息配置见表 2

    表  2  样本配置信息
    Table  2.  Information of sample configuration
    sample type training set verification set test set total
    positive 500 250 750 1500
    negative 500 250 125 875
    total 1000 500 875 2375
    下载: 导出CSV 
    | 显示表格

    普通的CNN方法,原始图像经过卷积层和池化层后,全链接层输出结果,本文方法的结构如图 3所示,引入压缩激励的过程,压缩特征图的空间信息,并通过激励操作学习通道间的依赖关系,可自适应分配每个通道的权重值,提取有利于任务的重要特征通道,最终能进一步增强网络模型的特征提取能力,采用SE-DenseNet-169框架的Faster R-CNN模型。

    图  3  改进的Faster R-CNN结构
    Figure  3.  Structure of improved Faster R-CNN

    改进模型主要完成绝缘子异常状态的精准识别,首先对红外原始图像进行相关修正与补偿实现样本扩充,然后采用本文方法进行训练,收敛后,获得最终的改进Faster R-CNN模型。

    CNN学习中,精确度的衡量一般会采用准确率(Precision)和召回率(Recall),其计算过程分别见式(11)和式(12):

    $$ {P_{{\rm{re}}}} = \frac{{{\rm{TP}}}}{{{\rm{TP}} + {\rm{FP}}}} $$ (11)
    $$ {P_{{\rm{ca}}}} = \frac{{{\rm{TP}}}}{{{\rm{TP + FN}}}} $$ (12)

    式中:TP表示是实际值和预测值均是异常绝缘子的个数;FP表示预测值是异常绝缘子,实际值却不是的个数;FN表示是实际值是异常绝缘子,预测值却不是的个数。

    为进一步衡量改进模型的优劣,这里采用平均检测精度(mean Average Precision,mAP),其中AP等价于召回率和准确率形成曲线与横轴包围的几何图形的面积,对所有类别的AP求平均值即可得到mAP。

    基于样本数据,开展BP、Faster R-CNN以及本文方法的精确度和效率对比研究,不同方法的实验结果统计见表 3。Faster R-CNN和本文方法在Precision方面均明显优于BP方法,本文方法的Recall最高,mAP也最高,相对于BP提高了近10%,这说明经过改进的Faster R-CNN对于小目标的特征提取具有明显的优势。本文方法通过压缩激励结果,减少了数据量的计算,所以相对于其他方法,有更高的效率。

    表  3  不同方法的实验结果统计
    Table  3.  Statistics of experimental results by different methods
    Name Precision Recall mAP Time/s
    BP 93.5% 90.4% 80.3% 2.3
    Faster R-CNN 98.7% 95.3% 88.7% 1.2
    BFEM 99.2% 97.6% 90.2% 0.9
    下载: 导出CSV 
    | 显示表格

    绘制其准确率-召回率关系曲线,如图 4所示,可更加形象直观地反映出本文方法对绝缘子异常特征的提取优势,因为另外两种方法的曲线均被完整的覆盖,说明本文方法改进效果明显。

    图  4  准确率-召回率关系曲线
    Figure  4.  Relation curves of precision and recall

    常见绝缘子排列有单Ⅰ型、双Ⅰ型以及Ⅴ型。本文开展这3类绝缘子的红外图像研究,比较不同排列方式的诊断准确率,如图 5所示。根据电力标准DL/T 664-2008[19],图(b)和图(c)绝缘子端部明显发热,属于异常情况。

    图  5  不同类型绝缘子的红外图像
    Figure  5.  Infrared image of different types of insulators

    不同类型绝缘子的异常诊断准确率见表 4,准确率均较高,均在90%以上;Ⅰ型和Ⅴ型绝缘子的准确率明显优于双Ⅰ型绝缘子,这是因为双Ⅰ型绝缘子会出现两排绝缘子重叠的情况,对红外图像的研究造成一定的影响,为此无人机对于该种类型绝缘子的线路需开展多种角度的拍摄。

    表  4  绝缘子异常诊断的准确率
    Table  4.  Accuracy of insulator anomaly diagnosis
    Insulator type Abnormal total Detected number Accuracy
    Single Ⅰ 62 61 98.4%
    Double Ⅰ 47 44 93.6%
    31 31 100.0%
    下载: 导出CSV 
    | 显示表格

    本文提出一种改进的Faster R-CNN方法,引入激励压缩环节,搭建训练模型,完成绝缘子红外图像的异常诊断,并成功应用于电力现场运维。本文方法可高效并精准地识别出绝缘子的异常缺陷,mAP达到90.2%。研究结果可为输电线路绝缘子缺陷识别研究提供一定的参考。

  • 图  1   活塞间隙密封示意图

    Figure  1.   The model of seal clearance

    图  2   环形间隙二维轴对称模型

    Figure  2.   Two dimensional axisymmetric model of annular clearance

    图  3   流体域网格划分

    Figure  3.   Computational gride of fluid domain

    图  4   因素水平效应曲线

    Figure  4.   The effect curves at the factor level

    表  1   影响因素设置

    Table  1   Factors setting table

    Levels A/mm B/mm C/Pa D/(r/min) E/mm F/mm
    1 4 9.5 200000 1000 12 2
    2 8 13.5 500000 2000 14 2.6
    3 12 17.5 800000 3000 17 3.2
    4 16 21.5 1100000 4000 21 4
    下载: 导出CSV

    表  2   正交试验表

    Table  2   Orthogonal test table

    Test numbers Clearance thickness Sealing gap length Pressure difference Motor speed Clearance diameter Piston stroke Nullcolumn Leakage Sum of evaluation indexes yi
    1 1 1 1 1 1 1 1 1.27566E-05 0.99825
    2 1 2 2 2 2 2 2 4.06038E-05 0.98505
    3 1 3 3 3 3 3 3 1.78506E-04 0.91969
    4 1 4 4 4 4 4 4 3.01509E-04 0.86139
    5 2 1 1 2 2 3 3 1.25877E-04 0.94463
    6 2 2 2 1 1 4 4 1.39146E-04 0.93834
    7 2 3 3 4 4 1 1 2.89947E-04 0.86687
    8 2 4 4 3 3 2 2 1.01106E-03 0.52507
    9 3 1 2 3 4 1 2 9.60161E-04 0.54920
    10 3 2 1 4 3 2 1 5.26293E-04 0.75484
    11 3 3 4 1 2 3 4 9.09569E-04 0.57318
    12 3 4 3 2 1 4 3 4.62742E-04 0.78497
    13 4 1 2 4 3 3 4 8.76374E-04 0.58891
    14 4 2 1 3 4 4 3 8.34327E-04 0.60884
    15 4 3 4 2 1 1 2 1.31809E-03 0.37955
    16 4 4 3 1 2 2 1 8.93187E-04 0.58094
    17 1 1 4 1 4 2 3 8.39161E-04 0.60655
    18 1 2 3 2 3 1 4 6.59262E-05 0.97305
    19 1 3 2 3 2 4 1 2.42253E-05 0.99281
    20 1 4 1 4 1 3 2 9.06550E-06 1.00000
    21 2 1 4 2 3 4 1 1.86725E-03 0.11925
    22 2 2 3 1 4 3 2 6.24450E-04 0.70832
    23 2 3 2 4 1 2 3 3.42238E-04 0.84208
    24 2 4 1 3 2 1 4 6.29290E-05 0.97447
    25 3 1 3 3 1 2 4 7.04118E-04 0.67056
    26 3 2 4 4 2 1 3 1.02023E-03 0.52073
    27 3 3 1 1 3 4 2 2.47633E-04 0.88692
    28 3 4 2 2 4 3 1 6.38515E-04 0.70165
    29 4 1 3 4 2 4 2 2.11886E-03 0.00000
    30 4 2 4 3 1 3 1 1.82858E-03 0.13758
    31 4 3 1 2 4 2 4 9.14909E-04 0.57065
    32 4 4 2 1 3 1 3 6.08598E-04 0.71583
    下载: 导出CSV

    表  3   极差分析表

    Table  3   Range analysis table

    Clearance thickness/μm Sealing gap length/mm Pressure difference/Pa Motor speed/(r/min) Clearance diameter/mm Piston stroke/mm
    T1 7.3368 4.4774 6.7386 6.0083 5.7513 5.9779
    T2 5.919 5.6268 6.3139 5.4588 5.5718 5.5358
    T3 5.4421 6.0318 5.5044 5.3782 5.4836 5.574
    T4 3.5823 6.1443 3.7233 5.4348 5.4735 5.1925
    Range R 3.7545 1.667 3.0153 0.63011 0.27787 0.78542
    下载: 导出CSV

    表  4   方差分析表

    Table  4   Variance analysis table

    Sources of variation Sum of quadratic Degree of freedom Mean squared F value F Critical value Significant test
    A 0.9013 3 0.3004 18.222 (2.9912;4.6755) Highly significant
    B 0.2175 3 0.0725 4.3975 (2.9912;4.6755) Significant
    C 0.6667 3 0.2222 13.4785 (2.9912;4.6755) Highly significant
    D 0.0324 3 0.0108 0.6559 (2.9912;4.6755) Non-significant
    E 0.0062 3 0.0021 0.1255 (2.9912;4.6755) Non-significant
    F 0.0388 3 0.0129 0.7836 (2.9912;4.6755) Non-significant
    Null column 0.1173 3 0.0391 2.3707 (2.9912;4.6755) Non-significant
    Error 0.4122 25 0.0165 - - -
    Sum 1.9802 31 - - - -
    下载: 导出CSV
  • [1]

    Margaret M, Asantha C, David T, et al. Overview of the Origins Space telescope: science drivers to observatory requirements (Conference Presentation)[C]//Space Telescopes and Instrumentation, 2018: 1-8.

    [2]

    Timothy P M, Rishi N P, Jeremy D W, et al. Cryogenic packaging of an optomechanical crystal[J]. Optics Express, 2019, 35: 28782-2871.

    [3]

    ZHOUWJ, GAN Z H, ZHANG X, et al. Clearance loss analysis in linear compressor with CFD Method[C]//International Cryocooler Conference, 2010: 353-359.

    [4]

    Reed J, Bailey P B, Dadd M W, et al. Motor and thermodynamic losses in linear cryocooler compressors[J]. Aip Conference Proceedings, 2005, 823: 361-368.

    [5]

    Bailey P, Dadd M, Stone C, et al. Gas spring losses in linear clearance-seal compressors[C]//International Cryocooler Conference, 2007: 345-353.

    [6] 陈曦, 武卫东, 周志刚, 等.自由活塞式斯特林制冷机间隙密封技术研究[J].低温与超导, 2008, 36(5): 5-8. https://www.cnki.com.cn/Article/CJFDTOTAL-DWYC200805001.htm

    CHEN Xi, WU Weidong, ZHOU Zhigang, et al. Study on the clearance seal in free-piston Stirling cooler[J]. Cryogenic & Superconductivity, 2008, 36(5): 5-8. https://www.cnki.com.cn/Article/CJFDTOTAL-DWYC200805001.htm

    [7] 卢明.间隙密封式伺服液压缸密封特性研究及仿真[D].武汉: 武汉科技大学, 2010.

    LU Ming. Research and Simulation on The Gap Type Servo Hydraulic Cylinder Seals[D]. Wuhan: Wuhan University of Science and Technology, 2010.

    [8] 马诗曼, 陈曦, 李静.自由活塞式斯特林制冷机压缩活塞间隙密封泄漏的数值模拟[J].低温工程, 2011, 181(3): 24-28. https://www.cnki.com.cn/Article/CJFDTOTAL-DWGC201103008.htm

    MA Shiman, CHEN Xi, LI Jing. Numerical simulations of seal leakage of compression piston in free piston Stirling cooler[J]. Cryogenics, 2011, 181(3): 24-28. https://www.cnki.com.cn/Article/CJFDTOTAL-DWGC201103008.htm

    [9]

    QI Y X, MENG X Q, MU D F, et al. Study on mechanism and factors affecting the gas leakage through clearance seal at nano-level by molecular dynamics method[J]. Energy, 2016, 102: 252-259. DOI: 10.1016/j.energy.2016.02.087

  • 期刊类型引用(8)

    1. 李明超,闫宽,张聪,胡记伟,欧锴,陈绪兵. 用于激光软钎焊温度测量的高精度红外辐射测温装置. 红外技术. 2025(01): 108-114 . 本站查看
    2. 秦沛强,聂传杰,吝曼卿,卢永雄,张岸勤,何家冰. 磷矿巷道岩爆风险的可视化及特征增强研究. 矿业研究与开发. 2025(02): 123-131 . 百度学术
    3. 李贞,魏勇. 基于BP神经网络的红外测温补偿算法研究. 机械制造与自动化. 2023(01): 170-172+176 . 百度学术
    4. 曾飞,胡文祥,高彦鑫,宋杰杰. 基于激光扫描的输送带横向跑偏检测系统. 制造业自动化. 2023(05): 21-24 . 百度学术
    5. 杨帆,徐军,吴振生,孙明晓,金添. 基于Web端多节点红外热成像传感系统设计. 激光杂志. 2022(02): 154-157 . 百度学术
    6. 何翔. 非接触式检测装置综合实验设计. 电子技术与软件工程. 2022(07): 148-151 . 百度学术
    7. 谢彬棽. 基于双CAN总线的露天矿带式输送机速度自动化控制方法. 煤矿机械. 2022(12): 214-217 . 百度学术
    8. 武存喜. 回转窑焚烧设备退火温度模糊自适应控制技术. 工业加热. 2022(12): 27-31 . 百度学术

    其他类型引用(2)

图(4)  /  表(4)
计量
  • 文章访问数:  288
  • HTML全文浏览量:  63
  • PDF下载量:  70
  • 被引次数: 10
出版历程
  • 收稿日期:  2020-10-21
  • 修回日期:  2020-12-30
  • 刊出日期:  2021-01-19

目录

/

返回文章
返回
x 关闭 永久关闭

尊敬的专家、作者、读者:

端午节期间因系统维护,《红外技术》网站(hwjs.nvir.cn)将于2024年6月7日20:00-6月10日关闭。关闭期间,您将暂时无法访问《红外技术》网站和登录投审稿系统,给您带来不便敬请谅解!

预计6月11日正常恢复《红外技术》网站及投审稿系统的服务。您如有任何问题,可发送邮件至编辑部邮箱(irtek@china.com)与我们联系。

感谢您对本刊的支持!

《红外技术》编辑部

2024年6月6日