Processing math: 0%
高级检索

激光云粒子成像探测仪研制

史晓丁, 冯力天, 刘英, 杨泽后, 伍波, 石磊

史晓丁, 冯力天, 刘英, 杨泽后, 伍波, 石磊. 激光云粒子成像探测仪研制[J]. 激光技术, 2023, 47(1): 67-72. DOI: 10.7510/jgjs.issn.1001-3806.2023.01.010
引用本文: 史晓丁, 冯力天, 刘英, 杨泽后, 伍波, 石磊. 激光云粒子成像探测仪研制[J]. 激光技术, 2023, 47(1): 67-72. DOI: 10.7510/jgjs.issn.1001-3806.2023.01.010
SHI Xiaoding, FENG Litian, LIU Ying, YANG Zehou, WU Bo, SHI Lei. Laser cloud particle imaging probe[J]. LASER TECHNOLOGY, 2023, 47(1): 67-72. DOI: 10.7510/jgjs.issn.1001-3806.2023.01.010
Citation: SHI Xiaoding, FENG Litian, LIU Ying, YANG Zehou, WU Bo, SHI Lei. Laser cloud particle imaging probe[J]. LASER TECHNOLOGY, 2023, 47(1): 67-72. DOI: 10.7510/jgjs.issn.1001-3806.2023.01.010

激光云粒子成像探测仪研制

基金项目: 

国家重大科学仪器设备开发专项基金资助项目 2011YQ110059

详细信息
    作者简介:

    史晓丁(1984-), 男, 高级工程师, 主要从事激光大气探测技术的研究。E-mail: 13882266021@139.com

  • 中图分类号: TN249;O436.2

Laser cloud particle imaging probe

  • 摘要: 为了研究云粒子微观特性, 采用激光成像技术、高斯光束匀化处理技术、高分辨率阵列探测器以及图像抽帧还原算法, 研制出激光云粒子成像仪。通过地面模拟测量及实际飞行试验, 对仪器性能进行了测试。结果表明, 研制仪器可实现对25 μm~1550 μm云粒子的成像, 分辨率为25 μm, 粒径测量误差小于9.5%;实测的冰晶粒子形态符合云粒子自然分布规律, 粒子尺度谱符合指数分布, 决定系数在0.58~0.97之间。该研制仪器可为云降水微物理研究提供基础观测数据。
    Abstract: For studying the micro features of cloud particles, a laser cloud particle imaging probe was developed by using laser imaging technology, homogenizing technology of Gaussian beam, and a high resolution array detector and a frame-extracting algorithm of image. The performance of the instrument was tested through simulation measure and actual flight experiment. The results show that for cloud particles of 25 μm~1550 μm with 25 μm resolution can be imaged with the developed instrument, and the error of particle-size measured is less than 9.5%. The particle shape of ice crystal watched conforms to the natural distribution of cloud particles, the particle-size distribution conforms to the exponential function with a determination coefficient between 0.58 and 0.97.The developed instrument can provide basic observation data for the study of cloud and precipitation microphysics.
  • 激光焊接具有焊接变形小、热影响区小、能量密度高、效率高等优点,因此日益广泛地应用于工业生产中[1-4]。当高强度激光作用于熔融金属时,工件表面发生强烈的蒸发,当蒸发产生的压力超过流体静压和熔体的表面张力,熔体就会被推开产生匙孔,同时金属蒸汽从匙孔喷出[5-6], 因此产生了液相和气相的剧烈振荡,导致焊接过程极易出现不稳定状态。传统方法提出存在两种焊接状态:稳态和非稳态。FANG对钇铝石榴石固体激光焊接过程的研究表明,在全熔透基础上,匙孔穿透及未穿透两种模式并非突变,而是存在一个两种模式交替出现的过渡区[7]。这是因为激光焊接过程在多相耦合作用下,存在一种准稳态的自平衡机制,当外界干扰导致参量产生波动时(如激光能量密度减小),由于干扰幅度较小,不会导致焊接状态直接从稳态突变为非稳态,而是首先体现在某个相态的波动(如金属蒸汽体量减小),由于该相态对其它相态的作用机制(金属蒸汽减小导致孔内压力降低,引起匙孔深度及开口面积减小),并最终反作用于参量(匙孔深度减小,导致光束聚焦面积减小,能量密度提高),使得焊接过程再次恢复到稳定状态。只有当外界干扰持续对焊接过程产生影响,并最终突破准稳态的调控阈值,焊接过程才转为非稳态。在准稳态下,即使没有出现焊接缺陷,各种信号也会出现明显的波动。准确识别焊接过程中的准稳态,对于提早调整工艺参量、避免出现焊接缺陷具有重大意义。

    液态熔池和金属蒸汽不是单一地影响焊接稳定性,二者的振荡互相耦合,具有相关性。WANG等人认为,匙孔振荡是诱发金属蒸汽振荡的主要因素[8]。PANG等人通过建模方法研究了金属蒸汽和匙孔的动态行为,模拟结果表明,金属蒸汽的振荡与匙孔内的气体动力学密切相关[9]。金属蒸汽蒸发时,从匙孔喷出的颗粒对匙孔壁的能量吸收有影响,这可能影响焊接过程中匙孔的稳定性[10]。对多相振荡相关性的研究需要获取焊接过程中多种相态的信号,涉及到传感信号的采集技术。光学视觉监测方法目前被广泛使用[11-12],通过这种方法可以获得大量可靠的数据,诸如熔池、飞溅和金属蒸汽等多种形态特征。采用深度学习的方法研究激光焊接成为当今研究的热点[13-14]。GAO等人建立了包括光学信号、光谱信号、X光信号等6路信号的同步传感系统,通过关联分析,从16个关键特征中提取42个原始特征,建立基于关键特征的贝叶斯网络缺陷检测模型,实现对激光焊接状态的实时检测[15]。MA等人使用卷积神经网络和主动视觉技术自动检测焊缝的典型咬边、炸孔和焊穿缺陷[16]。HUANG等人通过视觉传感提取熔池正面的形状特征参量和背面熔宽信息作为训练集,建立反向传播神经网络模型有效地预测焊缝的熔透状态[17]。ZHANG等人利用包括激光成像系统、视觉成像系统、光谱测量仪、光电传感器的多传感器系统分析焊接状态信息,建立深信度网络模型对焊接状态进行监测,该模型具有更高的精度和鲁棒性[18]

    本文中提出了一种基于多传感信号的相关性及神经网络的焊接稳态识别方法,该方法通过光电传感器获得可见光强度和反射光强度,通过高速摄像机获得匙孔面积和金属蒸汽面积,利用相关分析法计算4种信号特征两两之间的相关系数,研究了信号特征的相关性随焊缝的稳态类型变化的规律。根据稳态类型与信号特征的相关系数的对应关系,对每一组相关系数进行标签分类,相关系数与分类标签作为训练集,构建全连接神经网络进行模型训练,实现了对焊接稳态类型准确的预测。

    实验在配备了大功率盘形激光焊接装置通快16002的平台上进行,执行机构为安川六轴机器人,激光焊接装置结构如图 1所示。盘型激光波长为1030nm,激光束直径为300μm,激光功率的范围为2kW~16kW,聚焦位置范围为-4mm~4mm,焊接速率范围为2m/min~4m/min,保护气体为氩气。焊接材料为304不锈钢钢板,厚度为6mm。数据采集系统主要由光电传感器模块和工业成像模块组成。具体来说,扫描激光头通过单向反射镜和聚焦透镜收集来自焊接区的光辐射,并将其传输到光电传感器模块,两个光电传感器为了捕获可见光和反射激光的光强信号,在激光头中预先安装了一个分光镜,来自焊接区域的光信号通过光电二极管盒内的分光镜传输到可见光光电二极管(波长380nm~780nm)和反射激光光电二极管(波长1030nm), 最后经放大后由示波器接收,信号采样频率为500kHz。对于工业成像模块,将一个配备了窄带滤光片(波长为976nm)的高速黑白摄像机安置于与水平方向成60°角的位置,以捕捉焊件上方的匙孔和熔池的图像。为了得到清晰的金属飞溅与金属蒸汽动态图像,将配备了可见光感应滤光片的高速彩色摄像机安装在垂直于焊接方向的位置。高速摄像机的采集帧频为5000frame/s,图像分辨率为512pixel×512pixel。

    图 1 Experimental equipment for laser welding
    图  1  Experimental equipment for laser welding

    控制功率变化的焊接实验参量如表 1所示。焊接速率为3m/min,离焦量为-3mm, 保护气流量为30L/min,通过变化功率,获得稳态、准稳态、非稳态3种焊接状态。

    表  1  Welding experimental parameters for controlling power variation
    experimental number power/ kW speed/ (m·min-1) defocus/ mm gas/ (L·min-1)
    1 6 3 -3 30
    2 6 3 -3 30
    3 6 3 -3 30
    4 7 3 -3 30
    5 8 3 -3 30
    6 9 3 -3 30
    7 10 3 -3 30
    8 11 3 -3 30
    9 12 3 -3 30
    10 13 3 -3 30
    11 14 3 -3 30
    下载: 导出CSV 
    | 显示表格

    金属蒸汽面积提取的图像处理流程如图 2a所示。首先对彩色摄像机拍摄的图片进行感兴趣区域(region of interest, ROI)提取,得到金属蒸汽的图像,周围还有少量细小飞溅;然后对ROI区域进行中值滤波,中值滤波将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值[19],从而消除大部分飞溅造成的孤立噪声点;但是经过中值滤波后无法消除所有飞溅噪声,因此再对滤波后的图像进行二值化和开运算处理。开运算对图像进行先腐蚀后膨胀的操作。腐蚀可以消除边界点,使得边界向内部收缩。膨胀可以将与边界接触的点合并进来,使得边界向外部扩张。经过开运算后可以平滑边界,消除剩余的细小飞溅,但不会明显改变金属蒸汽的面积大小。金属蒸汽的面积通过处理后图像的像素个数和彩色摄像机的像素分辨率(0.012345mm/pixel)的乘积来估计。匙孔面积提取的图像处理流程与金属蒸汽面积提取的流程类似,如图 2b所示。因为匙孔图像中由反光引起的大噪声较多,无法通过中值滤波和开运算完全消除,所以匙孔面积提取比金属蒸汽面积提取多一个ROI提取的步骤。首先获得多幅匙孔区域的图像,然后计算每一幅图像中匙孔区域的外接圆的圆心和半径并得出圆心和半径的均值,即可构建出匙孔区域的外接圆掩模,以实现匙孔的分割和提取。匙孔面积通过处理后图像的像素个数和黑白摄像机的像素分辨率(0.000597mm/pixel)的乘积来估计。

    图 2 Feature extraction
    图  2  Feature extraction
    a—extraction of keyhole area b—extraction of plume area

    观察上述每一组实验的焊缝表面质量与可见光强、反射光强、匙孔面积以及金属蒸汽面积的幅值波动程度,将焊接状态划分为3种状态。图 3a所示为稳态下4种信号特征的波动图像。每一幅图中间的细窄型曲线显示了每1000个原始数据点幅值的平均值,当焊缝处于稳态时,焊缝表面质量良好,4种信号特征的幅值和均值的波动较为平稳,金属蒸汽面积的波动较其它3种信号特征的波动稍微大一些,但金属蒸汽面积的均值大部分比较平稳,均值的波峰和波谷横跨的时间范围较小。图 3b所示为稳态的状态机理图。匙孔内的压力维持平衡,匙孔与金属蒸汽振荡均匀。

    图 3 Laws of signals in steady status
    图  3  Laws of signals in steady status
    a—features laws b—status mechanism

    图 4a所示为准稳态下4种信号特征的波动图像。焊缝质量良好,没有明显的缺陷,但信号特征的波动较为剧烈,可见光光强的幅值和均值在0.1s~0.3s时波动剧烈,出现较多较大的峰值;匙孔面积在0.4s~0.65s时波动剧烈;对比稳态下的金属蒸汽面积,准稳态下金属蒸汽面积在整个焊接过程中波动都较为剧烈,均值离散程度大,均值的波峰和波谷横跨的时间范围较大。图 4b所示为准稳态的状态机理图。激光功率增加使匙孔内的压力增大,匙孔的振荡和金属蒸汽振荡更剧烈,金属蒸汽体量有所增加,可见光强也增加。

    图 4 Laws of signals in quasi-steady status
    图  4  Laws of signals in quasi-steady status
    a—features laws b—status mechanism

    图 5a所示为非稳态下4种信号特征的波动图像。焊缝表面质量明显不佳,4种信号特征的波动非常剧烈,方框区域1中4种信号特征的幅值先增后减,焊缝出现焊瘤缺陷,方框区域2中4种信号特征的幅值上下来回振荡,焊缝出现凹陷缺陷。图 5b所示为非稳态的状态机理图。匙孔内的压力增大,使得匙孔的振荡非常剧烈,匙孔的振荡导致喷出的金属蒸汽振荡也非常剧烈,并且熔池向上凸起和向下收缩的幅度更大,易产生焊瘤和凹陷缺陷。根据焊缝的表面质量与各种信号特征的幅值和均值波动的情况,可以基本判别焊缝所处的稳态类型。

    图 5 Laws of signals in unsteady status
    图  5  Laws of signals in unsteady status
    a—features laws b—status mechanism

    在焊接过程中,各种信号会相互耦合,信号的变化具有一定的相关性,而焊接信号的相关性间接反映了焊接质量。为了量化焊接过程的不同状态,利用信号相关分析法。相关性系数可以用来统计两个变量之间的相关性。信号X与信号Y的相关系数可以描述为:

    r = \frac{{\sum\limits_{i = 1}^n {\left( {{\mathit{\boldsymbol{X}}_i} - \bar X} \right)\left( {{\mathit{\boldsymbol{Y}}_i} - \bar Y} \right)} }}{{\sqrt {\sum\limits_{i = 1}^n {{{\left( {{\mathit{\boldsymbol{X}}_i} - \bar X} \right)}^2}} } \sqrt {\sum\limits_{i = 1}^n {{{\left( {{\mathit{\boldsymbol{Y}}_i} - \bar Y} \right)}^2}} } }} (1)

    式中,XY分别为矩阵Xi和矩阵Yi的均值,i为当前变量的数据点序号,n为当前变量的数据点总个数,相关系数r的取值范围是[-1, 1]。

    根据上述焊缝状态分析的方法,对11组实验的焊后状态进行划分,图 6所示为各信号特征之间的相关性随激光功率的变化情况。除了可见光强与反射光强的相关系数随功率变化没有明显规律外,其它信号特征之间的相关系数随功率的增加总体上呈上升趋势,但各信号特征之间的相关系数并非单调递增。激光功率为7kW和12kW时焊缝各信号特征之间的相关系数出现拐点,因为激光功率为7kW和12kW时焊缝处于准稳态,7kW前后的两组焊缝处于稳态,12kW前后的两组焊缝处于非稳态。关注激光功率为7kW处的前后图像,当焊缝从稳态变为准稳态,各信号之间的相关系数增加;当焊缝从准稳态变为稳态,各信号之间的相关系数减少。关注激光功率为12kW处的前后图像,当焊缝从非稳态变为准稳态,各信号之间的相关系数减少;当焊缝从准稳态变为非稳态,各信号之间的相关系数增加。值得注意的是,当激光功率为8kW, 9kW和10kW时,焊缝分别处于稳态、准稳态和非稳态,各信号特征之间的相关系数均单调递增,其变化趋势与前述7kW和12kW处相关系数随不同稳态的变化趋势一致。因此可以看出,焊缝各信号特征的相关系数与焊缝的稳态有着密切的联系。

    图 6 Correlation coefficients between the signals varies with the laser power
    图  6  Correlation coefficients between the signals varies with the laser power
    a—correlation between visible light and reflected light b—correlation between visible light and keyhole area c—correlation between visible light and plume area d—correlation between reflected light and keyhole area e—correlation between reflected light and plume area f—correlation between keyhole area and plume area

    焊缝各信号特征之间的相关系数随功率的变化规律不显著,随稳态的变化具有一定规律。为了研究焊缝稳态类型的变化对各信号特征相关性的影响,对图 6中激光功率从7kW~12kW之间的数据点进行平移,处于稳态的数据点向前移动,处于非稳态的数据点向后移动,处于准稳态的数据点置于稳态和非稳态之间,根据稳态的类型重新排布图像,使处于同一稳态类型的数据点相邻显示,图 7所示为各信号特征之间的相关系数随焊缝稳态类型的变化情况。

    图 7 Correlation coefficients between the signals varies with the statuses
    图  7  Correlation coefficients between the signals varies with the statuses
    a—correlation between visible light and reflected light b—correlation between visible light and keyhole area c—correlation between visible light and plume area d—correlation between reflected light and keyhole area e—correlation between reflected light and plume area f—correlation between keyhole area and plume area

    图 7a为可见光强与反射光强的相关系数随稳态的变化情况。可见光强与反射光强的相关系数随功率变化的范围很小,数值在0.3~0.4之间波动,属于弱相关,可见光强与反射光强的相关性无法区分焊缝处于哪一种稳态。图 7b为可见光强与匙孔面积的相关系数随稳态的变化情况。当焊缝处于稳态,相关系数的范围为(0,0.1);当焊缝处于准稳态,相关系数的范围为(0.1,0.2),焊缝处于稳态和准稳态时可见光强与匙孔面积均极弱相关;当焊缝处于非稳态,相关系数在0.4上下细微波动,可见光强与匙孔面积中等程度相关。图 7c为可见光强与金属蒸汽面积的相关系数随稳态的变化情况。当焊缝处于稳态,相关系数趋近于0,可见光强与金属蒸汽面积不相关;当焊缝处于准稳态,相关系数的范围为(0.1,0.25),可见光强与金属蒸汽面积极弱相关;当焊缝处于非稳态,相关系数的范围为(0.4,0.5),可见光强与金属蒸汽面积中等程度相关。图 7d为反射光强与匙孔面积的相关系数随稳态的变化情况。当焊缝处于稳态和稳准态,相关系数趋近于0,反射光强与匙孔面积不相关;当焊缝处于非稳态,相关系数的范围为(0.1,0.2),反射光强与匙孔面积极弱相关。图 7e为反射光强与金属蒸汽面积的相关系数随稳态的变化情况。当焊缝处于稳态,相关系数趋于0,反射光强与金属蒸汽面积不相关;当焊缝处于准稳态,相关系数的范围为(0,0.2),反射光强与金属蒸汽面积极弱相关;当焊缝处于非稳态,相关系数的范围为(0.1,0.2),反射光强与金属蒸汽面积极弱相关。图 7f为匙孔面积与金属蒸汽面积的相关系数随稳态的变化情况。当焊缝处于稳态,相关系数的范围为(0.2,0.3),匙孔面积与金属蒸汽弱相关;当焊缝处于准稳态,相关系数的范围为(0.4,0.5),匙孔面积与金属蒸汽面积中等程度相关;当焊缝处于非稳态,相关系数的范围为(0.6,0.7),匙孔面积与金属蒸汽面积强相关。

    根据以上分析发现,除了可见光强与反射光强的相关性无法区分不同稳态外,焊缝的其它信号特征之间的相关性随着稳态、准稳态、非稳态3个层级依次递增,在各自的稳态类型内,各信号特征之间的相关性在一定的范围内波动。匙孔面积和金属蒸汽面积的相关性能明显区分稳态和准稳态。可见光强与匙孔面积的相关性、可见光强与金属蒸汽面积的相关性、匙孔面积与金属蒸汽面积的相关性都能明显区分准稳态和非稳态。

    传统的激光焊接状态预测模型通常将传感器采集的熔池、飞溅、金属蒸汽等信号特征直接作为模型的输入来进行训练,根据第2节中的分析可知,当焊接状态从稳态向准稳态、非稳态转变时,焊接过程中的各种信号的相关性会不断增强。当输入的特征之间的相关性越高,则模型的泛化性能越差[20]。如果将传感器捕获的原始信号物理特征作为模型的输入,则模型对稳态的预测效果较好,但对非稳态的预测效果会较差。因此本文中将采用不同信号特征之间的相关系数作为模型的输入,而不选择用原始信号的物理特征,以避免不同稳态时原始信号物理特征之间相关性不同而导致的模型对不同稳态类型的预测性能不同。

    为了获得足够的训练集和测试集,每500个数据点计算一次不同信号特征之间的相关系数,由于在焊接开始和结束时焊接速度会发生变化,导致相关系数异常,因此这些相关系数被剔除,最终提取出5926组训练数据,643组测试数据。不同稳态类型的数据集数量如表 2所示。

    表  2  Number of data sets in different steady statuses
    status training sets number test sets number
    steady 3627 407
    quasi-steady 1095 108
    unsteady 1204 128
    total 5926 643
    下载: 导出CSV 
    | 显示表格

    神经网络模型的结构如图 8所示。输入层有6个神经元,分别为可见光强、反射光强、匙孔面积、金属蒸汽面积4种信号特征两两之间的相关系数,隐藏层有3层,神经元个数分别为512个、256个和128个,输出层为焊缝3种不同的稳态类型。模型的超参量分别是:训练轮数40轮,学习率0.001,每批次的数量(batch size)为24。此外,训练集划分出1/10用作验证集,验证集不参与训练,每一轮训练后,用验证集验证模型的预测性能。

    图 8 Fully connected neural network model structure
    图  8  Fully connected neural network model structure

    每一个输入数据在每一层神经元中传递时需经过如下表达式的变换:

    \mathit{\boldsymbol{Y}} = \varphi \left( {\sum\limits_{j = 1}^m {{\mathit{\boldsymbol{W}}_j}{\mathit{\boldsymbol{X}}_j} + b} } \right) (2)

    式中,矩阵Y为一层神经元的输出,矩阵Xj为输入的相关系数值,矩阵Wj为权重,b为偏置,j为当前的样本序号,m为训练集总样本数,φ为激活函数。3个隐藏层之间使用的激活函数是线性整流函数(rectified linear units,ReLU),其表达式为:

    \varphi \left( x \right) = {\rm{max}}(0, x) (3)

    式中,x是自变量。最后一个隐藏层与输出层之间使用的激活函数是归一化指数函数(别名softmax函数),其表达式为:

    {q_k} = \frac{{{\rm{exp}}{y_k}}}{{\sum\limits_{l = 1}^C {{\rm{exp}}{y_l}} }} (4)

    式中,C为输出类别的个数,yk为第k个输出类别的输出值,yl为第l个输出类别的输出值,qk为第k个输出类别的预测概率值。一组相关系数值输入后经过隐藏层计算,再通过softmax函数得出该组输入数据对应稳态、准稳态和非稳态3个输出类别的预测概率值,3个概率值加和为1,概率值最大的输出类别即为模型预测的类别。每经过一轮的训练,都要计算一次损失函数,用以调整权重矩阵Wj和偏置b。损失函数采用交叉熵损失函数,表达式为:

    L\left( {p, q} \right) = - \sum\limits_{k = 1}^C {{p_k}{\rm{ln}}{q_k}} (5)

    式中,L(p, q)为损失函数值,pk为第k个输出类别的实际概率值(1或者0)。从后往前逐层计算损失函数分别对权重和偏置的偏导数,原始权重减去损失函数对权重Wj的偏导数与学习率的乘积得到更新的权重,原始偏置减去损失函数对偏置b的偏导数与学习率的乘积得到更新的偏置。焊接稳态识别模型重复以上权重和偏置更新的过程,使损失函数值不断减小,以完成模型对稳态类型预测的学习过程。

    焊接稳态类型识别模型的训练过程如图 9所示。在训练28轮之后,训练准确率和验证准确率已经达到收敛,此时训练准确率为99.4%,验证准确率为98.96%。

    图 9 Training process of welding steady status recognition model
    图  9  Training process of welding steady status recognition model

    使用该焊接稳态识别模型对643组测试数据进行预测,预测结果如表 3所示。对稳态的预测准确率达到99.75%,对准稳态的预测准确率达到99.44%,对非稳态的预测准确率达到99.23%,结果表明,此模型对不同稳态类型的预测性能十分接近,并且整体的预测准确率达到98.76%。经检验,基于相关系数和神经网络构建的焊接稳态识别模型能够准确地根据焊缝各种信号特征之间的相关性,对焊缝所处的稳态类型进行识别。

    表  3  Prediction accuracy rate of welding steady status recognition model
    status test sets number correct prediction number accuracy/ %
    steady 407 406 99.75
    quasi-steady 108 102 99.44
    unsteady 128 127 99.23
    total 643 635 98.76
    下载: 导出CSV 
    | 显示表格

    为了解决在焊接过程中识别准稳态、避免出现非稳态的问题,利用相关分析的方法,对焊缝的光学和视觉特征信号的相关性进行研究,结果表明,随着稳态类型从稳态向准稳态、非稳态转变时,除了可见光强与反射光强的相关性没有明显变化外,焊缝其它信号特征之间的相关性依次递增,可见光强、匙孔面积、金属蒸汽面积3种信号特征两两之间的相关性可明显区分3种稳态。将信号特征的相关系数与稳态类别作为训练集,基于全连接神经网络构建了焊缝稳态识别模型,该模型对焊缝稳态类型的识别准确率在测试集上整体达到98.76%,具有较高的识别精度。

  • 图  1   线阵成像原理

    Figure  1.   Principle of linear array imaging

    图  2   不同帧率下的粒子图像

    Figure  2.   Particle images at different frame rate

    图  3   激光云粒子成像仪组成框图

    Figure  3.   Composition block diagram of laser cloud particle imaging probe

    图  4   匀化激光束

    Figure  4.   Homogenized beam

    图  5   高分辨率线阵探测器

    Figure  5.   Linear array detector of high resolution

    图  6   激光云粒子成像仪

    Figure  6.   Laser cloud particle imaging probe

    图  7   粒子模拟圆盘示意图

    Figure  7.   Diagram of simulated particle disk

    图  8   模拟成像测试

    Figure  8.   Simulated imaging test

    图  9   飞行轨迹(2021-05-31,内蒙古赤峰市)

    Figure  9.   Flight path(2021-05-31, Chifeng, Inner Mongolia)

    图  10   云粒子图像实测

    a—针状  b—片状  c—辐枝状  d—空间辐枝状  e—球形

    Figure  10.   Measured image of cloud particle

    a—acicular   b—schistic   c—dendritic   d—spatial dendritic   e—spherical

    图  11   飞行过程高度及气温变化曲线

    Figure  11.   Altitude and temperature variation curve during flight

    图  12   粒子数浓度曲线

    Figure  12.   Particle number concentration curve

    图  13   各粒径通道粒子总数

    Figure  13.   Total number of particles in each particle size channel

    图  14   粒子谱分布数据拟合(09:56:03~09:59:53)

    a—各组样本数据拟合曲线  b—决定系数

    Figure  14.   Particle spectrum distribution data fitting(09:56:03~09:59:53)

    a—fitting curve of each sampling data   b—coefficient of determination

    图  15   粒子谱分布数据拟合(10:38:03~11:46:46)

    a—各组样本数据拟合曲线  b—决定系数

    Figure  15.   Particle spectrum distribution data fitting(10:38:03~11:46:46)

    a—fitting curve of each sampling data   b—coefficient of determination

    表  1   阵列探测器参数

    Table  1   Parameters of array detector

    pixels of detector pixel size pitch responsivity
    64×1 120 μm×150 μm 150 μm 0.4 A/W
    下载: 导出CSV

    表  2   仪器主要技术指标

    Table  2   Specifications of the instrument

    particle size range resolution sampling area length wavelength imaging magnification data products
    25 μm~1550 μm 25 μm 80 mm 635 nm 6× particle image,size distribution
    下载: 导出CSV

    表  3   模拟成像粒径计算

    Table  3   Particle size calculation of simulated imaging

    nominal value/μm mean value/μm absolute error/μm absolute error/%
    100 90.5 -9.5 -9.5
    250 233.3 -16.7 -6.7
    500 494.1 -5.9 -1.2
    1000 1000.6 0.6 0.1
    1550 1550.0 0 0
    下载: 导出CSV

    表  4   不同条件下的云粒子分布

    Table  4   Cloud particles distribution under different conditions

    time altitude/m temperature/℃ phase particle morphology
    09:44~11:00 3900~4000 -6~-3 solid acicular,schistic,dendritic,spatial dendritic
    11:40~11:45 3700~1000 -1~5 liquid spherical
    下载: 导出CSV

    表  5   两次入云测量数据拟合结果

    Table  5   Data fitting results of two measurements through clouds

    No. time quantum number of samples σ λ coefficient of determination
    range average range average
    1 09:56:03~09:59:53 171 4.6~72 23 0.0036~0.0097 0.0059 0.67~0.97
    2 10:38:03~11:46:46 178 2.7~96 28 0.0032~0.019 0.0053 0.58~0.97
    下载: 导出CSV
  • [1] 郭学良, 方春刚, 卢广献, 等. 2008-2018年我国人工影响天气技术及应用进展[J]. 应用气象学报, 2019, 30(6): 641-650. https://www.cnki.com.cn/Article/CJFDTOTAL-YYQX201906001.htm

    GUO X L, FANG Ch G, LU G X, et al. Progresses of weather modification technologies and applications in China from 2008 to 2018[J]. Journal of Applied Meteorological Science, 2019, 30(6): 641-650(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-YYQX201906001.htm

    [2] 黄美元, 何珍珍, 沈志来. 暖性层积云中大云滴分布特征[J]. 气象学报, 1983, 41(3): 358-364. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB198303013.htm

    HUANG M Y, HE Zh Zh, SHEN Zh L. Distribution characteristics of large-size cloud droplets in warm straticumulus clouds[J]. Acta Meteor Sinica, 1983, 41(3): 358-364(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB198303013.htm

    [3] 郭学良, 于子平, 杨泽后, 等. 高性能机载云粒子成像仪研制及应用[J]. 气象学报, 2020, 78(6): 1050-1064. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB202006013.htm

    GUO X L, YU Z P, YANG Z H, et al. Development and application of the high-performance airborne cloud particle imager[J]. Acta Meteorologica Sinica, 2020, 78(6): 1050-1064(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB202006013.htm

    [4] 郑国光, 郭学良. 人工影响天气科学技术现状及发展趋势[J]. 中国工程科学, 2012, 14(9): 20-27. https://www.cnki.com.cn/Article/CJFDTOTAL-GCKX201209005.htm

    ZHENG G G, GUO X L. Status and development of sciences and technology for weather modification[J]. Strategic Study of CAE, 2012, 14(9): 20-27(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GCKX201209005.htm

    [5] 史晓丁, 伍波, 樊冬. 基于前向Mie散射的激光雾滴谱仪[J]. 激光杂志, 2020, 41(4): 116-119. https://www.cnki.com.cn/Article/CJFDTOTAL-JGZZ202004023.htm

    SHI X D, WU B, FAN D. Laser fog droplet spectrometer based on forward Mie-scattering[J]. Laser Journal, 2020, 41(4): 116-119(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JGZZ202004023.htm

    [6]

    KOROLEV A V, STRAPP J W, ISAAC G A. Evaluation of the accuracy of PMS optical array probes[J]. Journal of Atmospheric and Oceanic Technology, 1998, 15(3): 708-720.

    [7]

    NAGEL D, MAIXNER U, STRAPP W, et al. Advancements in techniques for calibration and characterization of in situ optical particle measuring probes and application to the FSSP-100 probe[J]. Journal of Atmospheric and Oceanic Technology, 2007, 24(5): 745-760.

    [8]

    FUGAL J P, SHAW R A. Cloud particle size distributions measured with an aiborne digital in-line holographic instrument[J]. Atmospheric Measurement Techniques, 2009, 2(1): 259-271.

    [9]

    BRENGUIER J L, BOURRIANNE T, COELHO A A, et al. Improvements of droplet size distribution measurements with the Fast-FSSP (Forward Scattering Spectrometer Probe)[J]. Journal of Atmospheric and Oceanic Technology, 1998, 15(5): 1077-1090.

    [10]

    LAWSON R P, O'CONNOR D, ZMARZLY P, et al. The 2D-S(stereo) probe: Design and preliminary tests of a new airborne, high-speed, high-resolution particle imaging probe[J]. Journal of Atmospheric and Oceanic Technology, 2006, 23(11): 1462-1477.

    [11]

    KNOLLENBERG R G. The optical array: An alternative to scattering or extinction for airborne particle size determination[J]. Journal of Applied Meteorology and Oceanic Technology, 1970, 9(1): 86-103.

    [12] 曾星, 伍波, 史晓丁, 等. 机载激光云粒子成像仪研制与校准研究[J]. 激光技术, 2015, 39(6): 798-801. DOI: 10.7510/jgjs.issn.1001-3806.2015.06.014

    ZENG X, WU B, SHI X D, et al. Study on airborne laser cloud particle imaging probe and its calibration[J]. Laser Technology, 2015, 39(6): 798-801(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2015.06.014

    [13] 杨军, 陈宝君, 银燕. 云降水物理学[M]. 北京: 气象出版社, 2011: 87-88.

    YANG J, CHEN B J, YIN Y. Physics of clouds and precipitation[M]. Beijing: China Meteorological Press, 2011: 87-88(in Chin-ese).

    [14]

    ZHU S C, GUO X L, LU G X, et al. Ice crystal habits and growth processes in stratiform clouds with embedded convection examined through aircraft observation in northern China[J]. Journal of the Atmospheric Science, 2015, 72(5): 2011-2032.

    [15] 牛生杰, 马铁汉, 管月娥, 等. 宁夏夏季降水性层状云微结构观测分析[J]. 高原气象, 1992, 11(3): 241-248. https://www.cnki.com.cn/Article/CJFDTOTAL-GYQX199203002.htm

    NIU Sh J, MA T H, GUAN Y E, et al. Observation analysis of the microstructures of precipitation sheet clouds in summer over Ningxia[J]. Plateau Meteor, 1992, 11(3): 241-248 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GYQX199203002.htm

    [16] 秦彦硕, 蔡淼, 刘世玺, 等. 华北秋季一次低槽冷锋积层混合云宏微物理特征与催化响应分析[J]. 气象学报, 2017, 75(5): 835-849. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB201705013.htm

    QIN Y Sh, CAI M, LIU Sh X, et al. A study on macro and micro physical structures of convective-stratiform mixed clouds associated with a cold front in autumn and their catalytic responses in North China[J]. Acta Meteor Sinica, 2017, 75(5): 835-849(in Chin-ese). https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB201705013.htm

    [17] 张佃国, 郭学良, 付丹红, 等. 2003年8~9月北京及周边地区云系微物理飞机探测研究[J]. 大气科学, 2007, 31(4): 596-610. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK200704004.htm

    ZHANG D G, GUO X L, FU D H, et al. Aircraft observation on cloud microphysics in Beijing and its surrounding regions during August-September 2003[J]. Chinese Journal of Atmospheric Sciences, 2007, 31(4): 596-610 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK200704004.htm

    [18]

    KOROLEV A V, ISAAC G A. Shattering during sampling by OAPs and HVPS. PartⅠ: Snow particles[J]. Journal of Atmospheric and Oceanic Technology, 2005, 22(5): 528-542.

    [19] 黄敏松, 雷恒池, 金玲. 机载云降水粒子成像仪所测数据中伪粒子的识别[J]. 大气科学, 2017, 41(5): 1113-1124. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK201705016.htm

    HUANG M S, LEI H Ch, JIN L. Pseudo particle identification in the image data from the airborne cloud and precipitation particle image probe[J]. Chinese Journal of Atmospheric Sciences, 2017, 41(5): 1113-1124(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK201705016.htm

    [20] 陈万奎, 严采蘩. 云中冰晶尺度谱特征与分布函数[J]. 气象, 1987, 13(11): 13-17. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXX198711002.htm

    CHEN W K, YAN C F. The size distribution of ice crystal within clouds[J]. Meteorological Monthly, 1987, 13(11): 13-17(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-QXXX198711002.htm

  • 期刊类型引用(10)

    1. 贾娜,余本军,张纯朴,王春昕,刘九庆. 选区激光熔化WC-12Co单道成型工艺参数优化. 激光技术. 2025(01): 113-120 . 本站查看
    2. 卞宏友,王美男,刘伟军,邢飞,王慧儒,徐效文,霍庆生. DZ125合金激光沉积CoCrW涂层的组织与性能. 热加工工艺. 2024(19): 121-127+131 . 百度学术
    3. 李镭昌,魏昕. 激光熔覆复合涂层WC对裂纹产生机理影响研究. 激光技术. 2023(01): 52-58 . 本站查看
    4. 石圆圆,罗玉凤. 轻轨建筑钢结构的表面防护与性能研究. 电镀与精饰. 2023(03): 60-67 . 百度学术
    5. 杨文斌,李仕宇,肖乾,陈道云,王溯,张博. 减摩耐磨激光熔覆涂层的研究现状及发展趋势. 润滑与密封. 2023(04): 171-182 . 百度学术
    6. 蒋瑞鑫,牛宗伟,史程程,任智强,韩国峰,杨保伟,王文宇,杨善林,陈贺连. 镍基高温合金载能束增材修复技术研究现状. 材料导报. 2023(15): 188-199 . 百度学术
    7. 晁祥瑞,黄勇,陈子鹏,许学虎,李文建,王宁,张志虎. 激光重熔对In718熔覆层组织与性能的影响. 激光技术. 2023(04): 506-512 . 本站查看
    8. 陆靖,孙文磊,陈子豪,邢学峰,杨凯欣,周浩南,刘德明. 热作模具表面激光熔覆H13的数值模拟及实验研究. 激光技术. 2023(04): 558-564 . 本站查看
    9. 胡桂领,师鹏,张磊. 数控机床高速钢刀具激光熔覆Co-WC的组织与切削加工性能. 激光与光电子学进展. 2022(11): 350-357 . 百度学术
    10. 刘琛,穆星宇,李金华,刘斌. 基于灰色理论激光熔覆对形貌影响与优化. 辽宁工业大学学报(自然科学版). 2022(06): 367-372 . 百度学术

    其他类型引用(5)

图(15)  /  表(5)
计量
  • 文章访问数:  5
  • HTML全文浏览量:  1
  • PDF下载量:  5
  • 被引次数: 15
出版历程
  • 收稿日期:  2021-11-29
  • 修回日期:  2022-01-03
  • 发布日期:  2023-01-24

目录

/

返回文章
返回