高级检索

焊缝表面缺陷激光视觉传感HOG-SVM的检测方法

胡丹, 吕波, 王静静, 高向东

胡丹, 吕波, 王静静, 高向东. 焊缝表面缺陷激光视觉传感HOG-SVM的检测方法[J]. 焊接学报, 2023, 44(1): 57-62, 70. DOI: 10.12073/j.hjxb.20211231001
引用本文: 胡丹, 吕波, 王静静, 高向东. 焊缝表面缺陷激光视觉传感HOG-SVM的检测方法[J]. 焊接学报, 2023, 44(1): 57-62, 70. DOI: 10.12073/j.hjxb.20211231001
HU Dan, LYU Bo, WANG Jingjing, GAO Xiangdong. Study on HOG-SVM detection method of weld surface defects using laser visual sensing[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 57-62, 70. DOI: 10.12073/j.hjxb.20211231001
Citation: HU Dan, LYU Bo, WANG Jingjing, GAO Xiangdong. Study on HOG-SVM detection method of weld surface defects using laser visual sensing[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 57-62, 70. DOI: 10.12073/j.hjxb.20211231001

焊缝表面缺陷激光视觉传感HOG-SVM的检测方法

基金项目: 国家自然科学基金(52275317);广东省自然科学基金项目(2023A1515012172);广州市科技计划项目(2023B03J1326).
详细信息
    作者简介:

    胡 丹,女,1995年出生,研究生,主要研究方向为无损检测. Email: 490048642@qq.com

    通讯作者:

    高向东,男,1963出生,教授,博士生导师;Email: gaoxd@gdut.edu.cn.

  • 中图分类号: TG 409

Study on HOG-SVM detection method of weld surface defects using laser visual sensing

  • 摘要: 为了实现对焊缝表面缺陷的自动检测与分类,研究一种有效识别焊缝表面缺陷的激光视觉检测方法. 通过激光视觉传感器采集焊缝图像并进行预处理,包括图像分割,灰度化,平滑去噪以及焊缝轮廓提取. 采用方向梯度直方图(Histogram of Oriented Gradient, HOG)提取焊缝激光条纹轮廓图像的特征向量. 其次,基于5折-交叉验证网格搜索方法进行模型参数寻优,最终建立了支持向量机(Support Vector Machine, SVM)智能模型识别与分类焊缝表面缺陷. 通过调整焊缝轮廓提取算法、HOG特征维度得到不同特征数据并进行对比、分析焊缝缺陷的识别效果. 在相同试验条件下,发现支持向量机比随机森林分类器、K最近邻分类器以及朴素贝叶斯分类器的识别率更高,达到97.86%. 基于HOG-SVM的焊缝表面缺陷智能识别方法可有效提高焊缝缺陷(气孔、凹陷、咬边)及无缺陷的分类精度.
    Abstract: In order to realize automatic detection and classification of weld surface defects, an effective laser vision detection method for weld surface defects was studied. First, the weld image is collected by a laser vision sensor and preprocessed, including image segmentation, grayscale, smooth denoising and weld contour extraction. Then the feature vectors of laser stripe contour image of weld seam are extracted by means of Histogram of Oriented Gradient (HOG), the model parameters were optimized based on the five-fold cross-validation grid search method. Finally, an intelligent model of Support Vector Machine (SVM) was established to identify and classify weld surface defects. Different feature data were obtained by adjusting the weld contour extraction algorithm and HOG feature dimension, and then the identification effect of weld defects was analyzed by comparing each other. Under same experimental conditions, it is found that the recognition rate of SVM is higher than that of random forest classifier, K-nearest neighbor classifier and naive BAYES classifier, reaching 97.86%. The proposed intelligent identification method of weld surface defects based on HOG-SVM can effectively improve the classification accuracy of weld defects (porosity, sag, undercut) and non-defects.
  • 电弧增材制造技术是根据零件的数字三维模型,通过电弧为热源将金属焊丝熔化,并逐层堆积出零件的新型先进数字化制造技术.电弧增材制造的热输入和熔池体积大于激光增材,其成形精度相对更低,但成形效率高,制造成本低,适用于铝合金、钛合金、高强钢等材料的中大型尺寸复杂金属零构件的直接制造[1-2].6061铝合金属于Al-Mg-Si系,具有中等强度、抗腐蚀性强、焊接性能好等特点;在固溶加人工时效处理的条件下,具有较高的硬度、较好的塑性以及优良的冷加工性,可用作结构材料[3],因此6061铝合金被广泛应用在航天航空、船舶、农业、汽车等领域中作为中大型、强度和耐蚀性能要求高的工业型材[4].

    目前铝合金的电弧增材制造的研究主要有两个方面:成形件的形貌控制以及其力学性能控制.在形貌控制上主要研究电流方式、焊枪角度等工艺条件对成形件形貌及尺寸的影响规律[5-7].在力学性能控制上研究通过辅助措施或热处理来优化成形件微观组织和力学性能[8-9].而关于6061铝合金电弧增材制造的研究至今未见文献报道,仅见一些关于6061铝合金焊接接头性能的研究[10-13],其焊接接头的拉伸性能一般小于基材,并且在焊道与基材交界处硬度高于焊道区域.

    文中开展对MIG电弧增材制造6061铝合金的工艺、微观组织和力学性能的研究,设计试验探索不同工艺参数下的焊道形貌和堆积成形可行性,采用优化的工艺参数堆积成形薄壁件,并研究薄壁件不同区域的微观组织与力学性能,为电弧增材制造6061铝合金中大型零构件的工业应用提供工艺研究基础.

    试验所用原材料包括:6061铝合金焊丝(直径是1.2 mm)和6061铝合金焊接基板(450 mm × 150 mm × 8 mm),焊丝化学成分如表1所示.电弧增材制造系统主要由ABB工业机器人(IRB1600)、沙福全数字控制MIG焊机(DIGIPULS III 420)等组成,焊接模式选用2T + 一元化焊接成形.堆积结束后,分别在上部、中部和下部3个区域用线切割机取立方体金相试样(15 mm × 10 mm × 10 mm)和拉伸试样(依据国家标准GB2651—89,试样总长50 mm,标距部分长38 mm,宽度3 mm,厚度2 mm).金相试样经过砂纸研磨和高氯酸电解抛光后,用凯乐试剂腐蚀处理并在光学显微镜下观察微观组织.采用维氏硬度计测量试样截面沿堆积方向的硬度变化.在万能材料试验机上进行拉伸试验.拉伸测试结束后使用扫描电镜分别观察3个区域拉伸试样的断口形貌.所有试样的取样位置如图1所示,其中上部、中部、下部试样沿着水平方向,与焊接方向平行;竖直试样沿着堆积方向,与焊接方向垂直.

    表  1  焊丝化学成分(质量分数,%)
    Table  1.  Chemical composition of welding wire
    SiFeCuMnMgCrZnTiAl
    0.4 ~ 0.80.70.15 ~ 0.40.150.8 ~ 1.200.04 ~ 0.350.250.15余量
    下载: 导出CSV 
    | 显示表格
    图  1  拉伸及金相试样位置
    Figure  1.  Positions of tensile and metallographic samples

    电弧增材制造的焊道形貌受到多个工艺参数影响,其中受送丝速度和焊接速度的影响最大.为了探究电弧增材6061铝合金堆积可行的工艺参数窗口,设计一系列的探究试验(表2).其中焊接电压固定为20 V,焊丝伸出长度为15 mm,保护气体(纯氩气)流量为15 L/min,每层堆积结束后调整焊丝长度,确保焊丝底端到工件的距离为3 mm.利用机器人控制焊枪垂直于基板进行堆积,电弧增材制造工艺试验结果如图3所示.

    表  2  探究试验工艺参数
    Table  2.  Used process parameters in experiments
    编号送丝速度
    vs/(m·min−1)
    焊接速度
    v/(mm·s−1)
    1 ~ 421,2,3,4
    5 ~ 831.5,3,4.5,6
    9 ~ 1242,4,6,8
    13 ~ 1652.5,5,7.5,10
    17 ~ 2063,6,9,12
    21 ~ 2473.5,7,10.5,14
    25 ~ 2884,8,12,16
    下载: 导出CSV 
    | 显示表格

    图2所示结果表明,当MIG送丝速度小于3 m/min时,改变焊接速度,都不能完整成形,并且极易烧坏导电嘴;当送丝速度为3 ~ 4 m/min时,由于电流太小导致焊接过程不稳定,焊道会呈不规则的间歇环状,增大焊接速度则不规则的程度减弱.

    图  2  焊接工艺参数匹配下的单层单道成形
    Figure  2.  Single-layer single-pass deposition with matched welding parameters

    当送丝速度介于5 ~ 7 m/min,且送丝速度和焊接速度的比值在0.5 ~ 1倍之间时可以获得宽度适中、成形良好的焊道. 当比值在0.5以下时,送丝速度相对焊接速度偏小,导致热输入量小,且单位时间内金属填充量不足,故焊道形状细窄,不利于多层堆积成形. 若比值增大至1以上,送丝速度相对焊接速度过高,导致单位时间内金属填充量偏多,这将造成焊道变粗和容易发生熔池塌陷,亦不利于堆积成形.综上所述,电弧增材6061铝合金MIG焊接参数在送丝速度和焊接速度的比值在0.5 ~ 1之间,且送丝速度选择5 ~ 7 m/min之间时可获得良好焊道形貌.

    堆积微观组织和力学性能成形试样时,选择堆积形貌最佳的焊接参数.其中送丝速度为6 m/min,焊接速度为6 mm/s,采用逐层往复堆积的方式,层数为60层.图3为堆积成形薄壁件宏观形貌,图中可以看出薄壁件起弧端和熄弧端略微坍塌,薄壁件侧壁有明显的层间条纹特征.薄壁件两端坍塌是因为焊接启停时焊接参数的非线性变化,焊道首尾部支撑不良而在电弧冲击下引起的熔池失稳.

    图  3  堆积成形件的宏观形貌
    Figure  3.  Macro morphology of deposition

    堆积前先对基板进行预热,有利于焊丝的充分熔化减少误差;起弧之前保护气提前打开,使堆积起弧端焊前被纯氩气充分保护;电弧熄灭阶段,停止送丝并且保护气体延迟关闭,能够避免多余的焊丝堆积,有利于形成均匀的堆积层;每层堆积结束后用钢刷清理焊道表面的黑黄色氧化物,当成形件整体温度降至75 ℃左右时再进行下一道堆积.

    图4为腐蚀过后的金相试样表面,焊道层与层之间的交界处为结合层,结合层大致位于下一层堆积对已堆积焊道的重熔区域,对应的主要是堆积熔池的底部区域.截面上,除结合层外的其余区域为沉积层,沉积层主要对应的是堆积熔池的中上部区域.如图4所示,在亮光照射下,结合层和沉积层呈现出沿堆积高度方向灰白色带依次交替的形貌.

    图  4  腐蚀后试样表面
    Figure  4.  Surface of the sample after corrosion

    图5表示不同区域金相试样的微观组织.从图中可以发现,无论是沉积层还是结合层都呈现出各种尺寸大小的气孔多发的状态,但二者的气孔分布具有不同的特点,这与堆积成形时熔池不同区域具有不同的冷却速率有关.结合层对应的是熔池底部,靠近已有的冷却焊道和基板,热传导面积大,冷却速率快,凝固界面上析出的气体不易逃逸出来,从而在结合层形成了数量多、以小孔为主的密集气孔分布.而沉积层对应的是熔池中上部,冷却速率相对较慢,熔池对流更充分,更有利于气体逸出到熔池外,因此沉积层气孔数量比结合层少.部分从结合层逃逸出来的小气泡以及沉积层中析出的气体,不断上浮、聚合、长大,若最终未能从熔池表面逃逸出来,则被俘获在沉积层中,形成体积更大、数量相对更少的气泡分布[14].

    图  5  试样上部、中部和下部区域的微观组织
    Figure  5.  Microstructure of the upper, middle and lower regions

    在试样下部区域还出现了“液滴”状物质,这是由于在堆积成形过程中,当焊道的温度处在共晶温度与液相线温度之间时,材料会发生共晶反应,造成一种局部“液化”现象[15].相比于其它金属,铝合金有着更快的冷却速率,且电弧增材制造过程中有着复杂的热环境,因此有更能满足“液化”现象所需的条件:堆积成形过程中,若已堆积焊道对下一层堆积焊道的散热速率较慢,当温度达到液化线温度时,固态扩散会使组织中的β相溶解到α-Al基体中.在6061铝合金的电弧增材制造过程中,已堆积焊道对新一层焊道的冷却速率很快,导致这一过程会很快完成,从而使β相不能完全溶解到α-Al基体中,当温度达到共晶反应温度时,未溶解的β相会和α-Al基体发生共晶反应,形成共晶液相,焊道冷却后就会在原β相位置处形成“液化”现象[16].这种现象一般发生在堆积成形件的下部区域.

    硬度试样的显微硬度测试结果如图6所示,由图可知,随着堆积高度的增加,成形件整体显微硬度存在一定的波动但变化不大,上、中、下部区域沉积层的平均显微硬度值分别为78.3,77.4和75.4 HV,而结合层硬度值波动更明显,最低硬度值为67 HV,这是由于结合层存在气孔缺陷,导致硬度值下降.对比全部测试点发现,硬度最高出现在最顶端和最底端的沉积层,这是因为堆积开始和结束时接触基板和空气,散热条件良好,焊道冷却速度快,显微组织稳定,晶粒细化导致其硬度略有提高[17].

    图  6  显微硬度
    Figure  6.  Microhardness chart

    图1所示,分别在成形件上中下3区域各取5个拉伸试样测试,测试结果计算平均值如表3所示.表4则列出了其它焊接工艺得到6061铝合金焊接接头强度作为对比[10-13].可以看出,MIG电弧增材堆积的薄壁件屈服强度会略低,但抗拉强度高于其它工艺方式.这是因为,与其它焊接工艺相比,MIG电弧增材制造过程的稳定性更差,气孔缺陷产生的机率更高,降低了疲劳性能;而电弧增材多层堆积的热循环过程类似于对材料的一个反复热处理,提高了强韧性.由表3可知,在水平方向上不同区域的强度没有明显差异,其平均抗拉强度约为244.3 MPa,屈服强度约为123.5 MPa.与水平方向的均值相比,竖直方向上的抗拉强度和屈服强度均略低,而断后伸长率更高.这说明MIG电弧增材制造的6061铝合金存在各向异性.

    表  3  不同位置拉伸性能对比
    Table  3.  Comparison of tensile properties in different position
    位置抗拉强度Rm/MPa屈服强度ReL/MPa断后伸长率A(%)
    竖直225.6107.622.6
    上部240122.719.2
    中部251131.318.1
    下部239116.716.8
    下载: 导出CSV 
    | 显示表格
    表  4  不同工艺焊接接头性能对比
    Table  4.  Properties comparison of welding joints by different processes
    工艺条件抗拉强度Rm/MPa屈服强度ReL/MPa断后伸长率
    A(%)
    文献[10]-TIG1751209.2
    文献[11]-CMT185.415510.1
    文献[12]-EBW1911565.2
    文献[13]-FRW195121.9610
    下载: 导出CSV 
    | 显示表格

    图7为堆积成形薄壁件不同区域的拉伸断口形貌图.图7a,7b为竖直拉伸的断裂位置,竖直试样的断口一般在结合层,此处气孔密集数量众多.图中可以看出竖直断口处不仅有大量的中小型气孔,放大后还能看到大量的直径3 μm以下的微细气孔.相对地,在图7c,7d所示的水平试样位于沉积层的断口截面上,表现出数量相对更少、但体积却更大的气孔分布.此外竖直和水平试样的断口表面都均匀分布着细密的韧窝,这是典型的韧性断裂特征.电弧增材制造过程中,由于氢在固液态铝中的溶解度差别大,而且铝材料的导热率高,熔池冷却快,熔池中的氢不能及时逸出,以气体的形式存在于熔池中形成气孔,故气孔缺陷是铝合金电弧增材制造常见的现象,拉伸试样的气孔聚集位置会引起应力集中,更容易造成断裂[18].

    图  7  拉伸试样断口形貌
    Figure  7.  Fracture morphology of tensile samples. (a) vertical sample(low); (b) vertical sample(high); (c) horizontal specimen(low); (d) horizontal specimen (high)

    通过变形加工 + 热处理的方式可以使多数气孔闭合,大幅提高增材制造铝合金构件的疲劳性能[19].MIG电弧增材制造6061铝合金的后处理工艺将在后续研究中进行.

    (1) 6061铝合金的MIG电弧增材制造中,当送丝速度和焊接速度的比值介于0.5 ~ 1之间,且送丝速度在5 ~ 7 m/min之间时,可以获得良好的多层堆积形貌.

    (2) 焊道层与层之间交界处为结合层,其余区域为沉积层,结合层和沉积层呈现出沿堆积高度方向灰白色带依次交替的形貌,并都呈现出各种尺寸大小的气孔多发的状态.下部区域成形过程中温度更快达到共晶温度,会使β相来不及溶解与α-Al基体发生共晶反应,形成“液化”现象.

    (3) 沿着堆积方向硬度变化不大,结合层硬度低于沉积层,且硬度波动性更大.不同区域水平方向强度差异不大,竖直和水平方向断后伸长率分别为22.6%和18%,成形件的力学性能存在各向异性.

  • 图  1   焊缝质量检测平台示意图

    Figure  1.   Schematic diagram of weld quality inspection platform

    图  2   扫描采集焊缝示意图

    Figure  2.   Diagram of weld scanning collection

    图  3   四种缺陷样本的实物图及原始激光条纹图

    Figure  3.   Physical image, fringe image of four kinds of defects. (a) porosity; (b) sag; (c) undercut; (d) no defect

    图  4   细胞单元cell的梯度方向&梯度方向直方图

    Figure  4.   Gradient direction & gradient direction histogram of cell(a)Gradient direction(b)gradient direction histogram

    图  5   焊接缺陷条纹图像的HOG特征可视化效果

    Figure  5.   HOG feature visualization effect of weld defect streak image. (a) Porosity; (b) Sag; (c) Undercut; (d) No defect

    图  6   模型参数的交叉验证结果

    Figure  6.   Cross-validation results for the model parameters

    图  7   图像预处理及特征提取过程示意图

    Figure  7.   Schematic diagram of process of image preprocessing and feature extraction

    图  8   Roberts算子处理后的焊接缺陷识别结果

    Figure  8.   Results of weld defect identification using Roberts

    图  9   Canny算子处理后的焊接缺陷识别结果

    Figure  9.   Results of weld defect identification using Canny

    图  10   平均识别率对比图

    Figure  10.   Comparison chart of average recognition rate

    表  1   四种缺陷样本的实物图及原始激光条纹图

    Table  1   Physical image, fringe image of four kinds of defects

    步骤气孔凹陷咬边无缺陷
    ROI
    灰度化
    去噪图
    二值化
    轮廓图
    下载: 导出CSV

    表  2   不同边缘检测算子的试验结果

    Table  2   Experimental results of different edge detection operators

    边缘检测算子平均识别率 B (%)
    Roberts97.73
    Sobel97.17
    Prewitt97.22
    Laplacian97.28
    Canny96.50
    下载: 导出CSV

    表  3   不同cell大小的试验结果

    Table  3   Experimental results of different cell size

    Cell大小特征维数平均识别率 B (%)训练时间t/s
    4 × 41674097.284.6256
    8 × 8378097.391.1064
    16 × 1675696.650.7897
    32 × 3210894.720.3178
    下载: 导出CSV
  • [1]

    Gao X D, Ma N J, Du L L. Magneto-optical imaging characteristics of weld defects under alternating magnetic field excitation[J]. Optics express, 2018, 26(8): 9972 − 9983. doi: 10.1364/OE.26.009972

    [2]

    Gantala T, Balasubramaniam K. Automated defect recognition for welds using simulation assisted TFM imaging with artificial intelligence[J]. Journal of Nondestructive Evaluation, 2021, 40(1): 1 − 24. doi: 10.1007/s10921-020-00734-w

    [3]

    Lu Y, Jiang H Q. Weld defect classification in radiographic images using unified deep neural network with multi-level features[J]. Journal of Intelligent Manufacturing, 2020, 32(3): 459 − 469.

    [4]

    Kumar D, Verma D, Suryanarayana B, et al. Analysis of welding characteristics on stainless steel for the process of TIG and MIG with dye penetrate testing[J]. International Journal of Engineering and Innovative Technology, 2012, 2(1): 283 − 290.

    [5]

    Gao X D, Ding D, Bai T, Katayama S. Weld-pool image centroid algorithm for seam-tracking vision model in arc-welding process[J]. IET image processing, 2011, 5(5): 410 − 419. doi: 10.1049/iet-ipr.2009.0231

    [6] 谢志孟, 高向东. 基于Canny算子的焊缝图像边缘提取技术[J]. 焊接学报, 2006, 27(1): 29 − 32. doi: 10.3321/j.issn:0253-360X.2006.01.008

    Xie Zhimeng, Gao Xiangdong. Edge detection of weld image based on Canny operator[J]. Transactions of the China Welding Institution, 2006, 27(1): 29 − 32. doi: 10.3321/j.issn:0253-360X.2006.01.008

    [7]

    Gao X D, Mo L, Xiao Z L, et al. Seam tracking based on Kalman filtering of micro-gap weld using magneto-optical image[J]. The International Journal of Advanced Manufacturing Technology, 2016, 83(1-4): 21 − 32. doi: 10.1007/s00170-015-7560-x

    [8]

    Chi D Z, Gang T. Defect detection method based on 2D entropy image segmentation[J]. China Welding, 2020, 29(1): 45 − 49.

    [9]

    Sun J, Wang P, Luo Y K, et al. Surface Defects Detection Based on Adaptive Multiscale Image Collection and Convolutional Neural Networks[J]. IEEE Transactions on Instrumentation and Measurement, 2019, 68(12): 4787 − 4797. doi: 10.1109/TIM.2019.2899478

    [10]

    Gao X D, Sun Y, Katayama S. Neural network of plume and spatter for monitoring high-power disk laser welding[J]. International Journal of Precision Engineering and Manufacturing-Green Technology, 2014, 1(4): 293 − 298. doi: 10.1007/s40684-014-0035-y

    [11]

    Wang H F, Wu Z J, He Z C, et al. Detection of HF-ERW Process by 3D Bead Shape Measurement with Line-Structured Laser Vision[J]. IEEE Sensors Journal, 2021, 21(6): 7681 − 7690. doi: 10.1109/JSEN.2021.3049396

    [12]

    Han Y Q, Fan J F, Yang X Z. A structured light vision sensor for on-line weld bead measurement and weld quality inspection[J]. The International Journal of Advanced Manufacturing Technology, 2020, 106(1): 2065 − 2078.

    [13] 杨国威, 闫树明, 王以忠. 基于方向梯度直方图粒子滤波的V型焊缝跟踪[J]. 中国激光, 2020, 47(7): 330 − 338.

    Yang Guowei, Yan Shuming, Wang Yizhong. V-Shaped Seam Tracking Based on Particle Filter with Histogram of Oriented Gradient[J]. Chinese Journal of Lasers, 2020, 47(7): 330 − 338.

    [14] 周晓晓, 王克鸿, 杨嘉佳, 等. 电压近似熵-SVM铝合金双丝PMIG焊过程稳定性评价[J]. 焊接学报, 2017, 38(3): 107 − 111.

    Zhou Xiaoxiao, Wang Kehong, Yang Jiajia, et al. Process stability evaluation on aluminum alloy twin-wire PMIG welding by approximate entropy based SVM of voltage signa[J]. Transactions of the China Welding Institution, 2017, 38(3): 107 − 111.

  • 期刊类型引用(5)

    1. 金东升,曹文田,张英杰. 双目极线约束下的储罐环焊缝定位算法. 机械科学与技术. 2025(01): 1-8 . 百度学术
    2. 于强,肖澜,郑大伟,彭昭鸿,宋凯宏. 基于计算机视觉的锂电池连接片焊接检测系统. 焊接学报. 2024(10): 38-49 . 本站查看
    3. 李春凯,潘宇,石玗,王文楷,赵中博. 基于视觉传感的热丝激光金属沉积熔滴—熔池多特征信息同步监测. 焊接学报. 2024(11): 115-120 . 本站查看
    4. 乐健,李发园,束志恒,曾明如,张华. 基于视觉传感的焊接电流及电压检测控制方法. 焊接学报. 2024(11): 85-89 . 本站查看
    5. 贾瑞燕,李海超,魏方锴,徐勇,周宇飞. 基于三维点云的磁极焊缝识别及机器人轨迹生成技术. 焊接学报. 2024(11): 50-54 . 本站查看

    其他类型引用(2)

图(10)  /  表(3)
计量
  • 文章访问数:  444
  • HTML全文浏览量:  152
  • PDF下载量:  72
  • 被引次数: 7
出版历程
  • 收稿日期:  2021-12-30
  • 网络出版日期:  2022-12-14
  • 刊出日期:  2023-01-24

目录

/

返回文章
返回