CN114850658B - 基于反射光谱的激光器参数自动调控方法 - Google Patents

基于反射光谱的激光器参数自动调控方法 Download PDF

Info

Publication number
CN114850658B
CN114850658B CN202210531991.0A CN202210531991A CN114850658B CN 114850658 B CN114850658 B CN 114850658B CN 202210531991 A CN202210531991 A CN 202210531991A CN 114850658 B CN114850658 B CN 114850658B
Authority
CN
China
Prior art keywords
laser
target
parameters
spectrum
reflection spectrum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210531991.0A
Other languages
English (en)
Other versions
CN114850658A (zh
Inventor
季凌飞
曹丽杰
张犁天
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202210531991.0A priority Critical patent/CN114850658B/zh
Publication of CN114850658A publication Critical patent/CN114850658A/zh
Application granted granted Critical
Publication of CN114850658B publication Critical patent/CN114850658B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23KSOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00Working by laser beam, e.g. welding, cutting or boring
    • B23K26/36Removing material
    • B23K26/38Removing material by boring or cutting
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23KSOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00Working by laser beam, e.g. welding, cutting or boring
    • B23K26/02Positioning or observing the workpiece, e.g. with respect to the point of impact; Aligning, aiming or focusing the laser beam
    • B23K26/06Shaping the laser beam, e.g. by masks or multi-focusing
    • B23K26/062Shaping the laser beam, e.g. by masks or multi-focusing by direct control of the laser beam
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23KSOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00Working by laser beam, e.g. welding, cutting or boring
    • B23K26/70Auxiliary operations or equipment
    • B23K26/702Auxiliary equipment

Landscapes

  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Engineering & Computer Science (AREA)
  • Plasma & Fusion (AREA)
  • Mechanical Engineering (AREA)
  • Laser Beam Processing (AREA)

Abstract

基于反射光谱的激光器参数自动调控方法,涉及激光加工和自动化控制领域。该方法可在靶材被激光消融之前,精准调试出适合该靶材加工的激光参数。根据激光频率及通讯频率自定义光谱积分参数,采集激光消融多组分靶材过程中的反射光谱数据;提取光谱特征参数,以95%置信区间为标准建立光谱特征参数与此时待消融靶材所需激光参数的主成分回归模型;当激光消融过程中暴露出待去除靶材时,重复前述步骤,再次调试出适用于新靶材消融的激光参数;在激光消融过程中根据反射光谱数据对激光参数进行实时闭环调控,实现靶材的高度选择性、无损去除。在高效、高质、高精激光材料加工等领域具有重大的应用前景。

Description

基于反射光谱的激光器参数自动调控方法
技术领域
本发明涉及激光加工和自动化控制领域,尤其涉及基于反射光谱的激光器参数自动调控方法。
背景技术
激光加工具有精度高,效率快等优势。然而在多组分材料的消融过程中,由于加工靶材往往为非连续介质,其结构复杂多变、成分种类繁多,如果可以在待去除靶材被消融之前及时调整激光参数,有望大幅提升激光去除效果。传统加工过程中,操作人员需要关闭激光,根据经验区分待去除的靶材性质之后,再调整激光参数,虽然操作人员具备丰富的专业知识和很强的动手能力,但整个过程仍然耗时耗力,尤其是,无法保证激光参数调控过程的及时性和准确性。因此,如何实现靶材去除过程中激光参数的自动化调控,如何确保激光参数调控的及时性和靶材去除的无损性,成为当前多组分材料激光自动化加工面临的技术瓶颈。
已报道的方法通过监测激光与靶材之间相互作用造成的光学和声学过程发射,进而调控激光参数。文章《Real-time spectral response guided smart femtosecondlaser bone drilling》使用光谱仪对激光骨加工过程中喷射出的等离子体羽流进行在线监测,实时优化工艺参数;专利CN111136382A(公布日期:2020.05.12)中通过监测激光加工过程中的光声信号,调控优化激光加工参数。但是此类方法中监测信号的获取晚于待去除靶材的消融,即靶材被去除后才能对激光加工参数进行调整,无法实现监测的无损化和调控的及时性,如用在激光清洗文物表面的脏污和锈迹时及易伤及文物本身,造成极大的经济损失;并且少量光谱检测数据无法代表靶材的广域特征,现有的检测方法分辨率低,难以实现靶材的准确识别。非损伤性的探测方式,除了需要准确识别激光加工过程中靶材的变化,还需要在新靶材被消融前对激光参数做出调控,由于需要精确匹配监测反馈响应时间与激光调控响应时间的对应关系,难度极大,目前研究尚未报道。
为了克服上述问题,本发明提出一种基于反射光谱的激光器参数自动调控方法。本方法基于反射光谱探测原理,通过光谱仪监测宽带探测光源在待去除靶材上形成的反射光谱,建立光谱数据与激光参数的映射关系,实时调控激光器参数。由于所使用的宽带探测光源能量远低于靶材的损伤阈值,并且根据激光脉冲频率和通讯频率设置特定的光谱积分时间,可在靶材被消融之前及时调整工艺参数,保证了监测的无损性和调控的前置性,突破现有技术滞后性调控问题。
发明内容
本发明要解决的技术问题是克服激光加工过程中工艺参数无法在靶材被消融之前做出调控的问题,发明了一种基于反射光谱的激光器参数自动调控方法。本方法基于反射光谱探测原理,可根据靶材的变化及时调整激光参数,具有监测无损性、参数调控前置性、材料去除高效性等优点。
本发明专利为解决其技术问题所采用的方案是:
1.一种基于反射光谱的激光器参数自动调控方法,该方法包括如下步骤:
S1、根据激光频率与通讯频率自定义光谱采集积分参数,利用脉冲激光器对多组分靶材进行消融,采集反射光谱数据;
S2、对采集到的反射光谱进行预处理,剔除反射率<0的光谱波段,对筛选出的光谱进行非线性拟合,以95%置信区间为标准,筛选出建模谱带,再建立主成分回归模型,调试出适用于待消融靶材的激光参数;
S3、当脉冲激光消融过程中再次暴露出新的待去除靶材时,重复前述步骤,获取该待消融新靶材光谱数据,再次调试出适用于新靶材消融的激光参数;
S4、根据反射光谱数据调试出适用于靶材消融的激光参数后,即可根据所采集反射光谱数据对消融过程中的激光参数进行实时闭环调控,实现靶材的高度选择性、无损消融。
2.优选的,光谱采集系统的积分时间上限可由如下公式确定:
t≤(1/fL-1/fT)/2
其中,t为光谱采集的积分时间,fL为激光器输出的脉冲频率,fT为通讯频率。
3.优选的,所述对采集到的光谱进行预处理的流程为:采用平滑滤波方法对采集到的N组靶材反射光谱数据分别进行平滑处理,以降低各种原因导致的随机噪声,提高信噪比。
4.优选的,对光谱数据进行归一化处理,获取光谱特征信息中的波峰、波谷、斜率、包络面积、反射率。再利用收集到的光谱特征信息描述靶材信息,每个光谱样本信息用p个变量描述,N组靶材数据包含q个类别,该类有nq个样本,整个训练集样本数据矩阵可记为Mp×N,训练集光谱样本采用主成分模型进行分类,模型表示如下:
其中,为该类中对应的元素(i=1,2,…,p;k=1,2,…,nq),/>为q类中nq个样本中第i个变量的均值,m为该类样本数据主成分数,/>为训练集的光谱残差,b和cαk是为了使光谱残差达到极小值的参数。模型建立后,计算该q类样本的总体标准偏差,采用多元线性回归法拟合,再判别未知样本u是否属于q类,区分此模型与原训练集构建模型的不同,以/>表示拟合后的光谱残差,/>为该类中对应的元素(i=1,2,…,p;u=1,2,…,nq),/>为q类中nq个样本中第i个变量的均值,b和Cαk是为了使光谱残差达到极小值的参数。
随后,以主成分模型和拟合后的模型的标准差比值为检测计量M,当M≥0.8时样本u属于q类。
5.优选的,在建立光谱特征信息与待消融靶材的消融检测模型之前,需要根据靶材种类、切割尺寸要求、切割精度要求、切割瑕疵要求、切割效率要求确定适用于待消融靶材的最佳激光参数,这属于现有技术范畴;然后基于确立的主成分数,通过相似分析法建立激光消融检测模型,同时实时采集光谱信息,直至采集到光谱数据的主成分模型和拟合后的模型的标准差比值M≥0.8时为止,即可调制出适用于该靶材的最佳激光消融参数。
6.优选的,所述光谱采集系统所使用的非损伤性宽带探测光源的光谱范围为350-3000nm。
与现有技术相比,本发明有如下优点:
本方法基于反射光谱探测原理,通过光谱仪监测宽带探测光源在待去除靶材上形成的反射光谱,进而调控激光工艺参数。由于所使用的宽带探测光源能量远低于靶材的损伤阈值,并且根据激光脉冲频率和通讯频率设置特定的光谱积分时间,可在靶材被消融之前及时调整工艺参数,保证了监测的无损性和调控的前置性。
附图说明
图1为本发明提供的基于反射光谱的激光器参数自动调控方法流程图。
具体实施方式
以下结合实例对本发明进行详细说明。应当理解,此处所描述的具体实施列仅用于解释本发明,并不用于限定本发明。
实施例1:
本发明用于不锈钢-碳钢层合板的切割,不锈钢板的最佳切割激光工艺参数为:激光波长1064nm,脉宽600μs,平均功率2000W,碳钢板的最佳切割激光工艺参数为:激光波长10.6μm,脉宽300μs,平均功率1500W,可通过反射光谱在靶材被消融前及时调控激光工艺参数,以达到最佳切割效率及切割质量。具体实施方式如下:
加工激光的脉冲频率fL为100Hz,系统的通讯频率fT为125Hz,根据公式t≤(1/fL-1/fT)/2,光谱采集系统的积分时间t≤1ms,为提升光谱收集的信号反射率、提高光谱峰前切换速率,光谱采集的积分时间自定义为1ms;
利用脉冲激光器对不锈钢-碳钢层合板进行切割,通过光谱采集系统收集切割过程中待消融不锈钢的反射光谱数据;采用Savitzky-Golay卷积平滑滤波方法对采集到的1600组光谱数据进行预处理,剔除所选靶材反射率<0的光谱波段(350-450nm),再对筛选出的光谱带曲线进行非线性拟合、归一化处理,采用一阶微分方法和多元散射校正方法对数据进行裁剪,以降低环境因素导致的随机噪声,提高信噪比;采用原始光谱和经过一阶微分、二阶微分等方法预处理后的数据进行建模,以正确率位置表进行对比,结果如表1所示,原始光谱经过一阶微分方法和多元散射校正方法处理后,建立的模型光谱预测结果正确率较高。
表1不同预处理方法的正确率
然后,选取不同主成分数量后用交互验证方法确立主成分个数,发现主成分为m=3时,光谱残差最小。随后以95%置信区间为标准,筛选出建模谱带,采用偏差权重法选取光谱段位于1150~1550nm范围内,直至采集到光谱数据的主成分模型和拟合后的模型的标准差比值M=0.9时为止,调制出适用于该不锈钢靶材的最佳激光消融参数(激光波长1064nm,脉宽600μs,平均功率2000W)。模型的标准差比值M为:
其中,每个样本信息用p个变量描述,N组靶材数据包含q个类别,每类有nq个样本,m为该类样本数据主成分数,为训练集的光谱残差,/>是拟合后的光谱残差。
当激光消融过程中暴露出碳钢时,重复前述步骤,实时调控出适用于待消融碳钢的最佳激光参数(激光波长10.6μm,脉宽300μs,平均功率1500W),实现不锈钢-碳钢层合板的高度选择性、无损消融。
实施例2:
本发明用于离体羊的口腔颌面骨激光切割,可通过反射光谱调控激光参数:激光波长2.1μm,脉宽500μs,平均功率30W,功率密度150W/mm2,在激光消融神经等软组织之前及时关闭激光,达到只去除离体骨组织,不会损伤离体组织中的神经等软组织的目的。具体实施方式如下:
切割激光的脉冲频率fL为20Hz,系统的通讯频率fT为125Hz,根据公式t≤(1/fL-1/fT)/2,光谱采集系统的积分时间t≤21ms,为提升光谱收集的信号反射率、提高光谱峰前切换速率,光谱采集的积分时间自定义为20ms;
将离体的口腔颌面置于生理盐水溶液中,利用脉冲激光器对口腔颌面的离体组织进行去除,通过光谱采集系统收集去除过程中待消融骨组织的反射光谱数据;对采集到的1500组光谱数据进行预处理,剔除所选靶材反射率<0的光谱波段(550-650nm),再对筛选出的光谱带曲线进行非线性拟合、归一化处理,采用一阶微分方法和Whittaker平滑滤波方法对数据进行裁剪,以降低环境因素导致的随机噪声,提高信噪比;
然后,选取不同主成分数量后用交互验证方法确立主成分个数,发现主成分为m=5时,光谱残差最小。随后以95%置信区间为标准,筛选出建模谱带,采用偏差权重法选取光谱段位于1350~1800nm范围内直至采集到光谱数据的主成分模型和拟合后的模型的标准差比值M=0.86时为止,调制出适用于该靶材的最佳激光消融参数(激光波长2.1μm,脉宽500μs,平均功率30W,功率密度150W/mm2)。模型的标准差比值M为:
其中,每个样本信息用p个变量描述,N组靶材数据包含q个类别,每类有nq个样本,m为该类样本数据主成分数,为训练集的光谱残差,/>是拟合后的光谱残差。
当激光消融过程中暴露出离体组织的神经时,重复前述步骤,构建光谱特性与激光参数的动态模型,将激光参数调节至激光波长2.1μm,脉宽500μs,平均功率0W,功率密度0W/mm2
通过此方法,在离体组织被消融之前,即可实时调控出适用于待消融组织的最佳激光参数,为实现离体口腔颌面的高度选择性、无损消融提供了新方法。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (1)

1.基于反射光谱的激光器参数自动调控方法,其特征在于,包括如下步骤:
S1、根据激光频率与通讯频率自定义光谱采集积分参数,利用脉冲激光器对多组分靶材进行消融,采集反射光谱数据;
S2、对采集到的反射光谱进行预处理,剔除反射率<0的光谱波段,对筛选出的光谱进行非线性拟合,以95%置信区间为标准,筛选出建模谱带,再建立主成分回归模型,调试出适用于待消融靶材的激光参数;
S3、当脉冲激光消融过程中再次暴露出新的待去除靶材时,重复前述步骤,获取该待消融新靶材光谱数据,再次调试出适用于新靶材消融的激光参数;
S4、根据反射光谱数据调试出适用于靶材消融的激光参数后,根据所采集反射光谱数据对消融过程中的激光参数进行实时闭环调控;
对光谱数据进行归一化处理,获取光谱特征信息中的波峰、波谷、斜率、包络面积、反射率;再利用收集到的光谱特征信息描述靶材信息,每个光谱样本信息用p个变量描述,N组靶材数据包含q个类别,该类有nq个样本,整个训练集样本数据矩阵记为Mp×N,训练集光谱样本采用主成分模型进行分类,模型表示如下:
其中,为该类中对应的元素,其中i=1,2,…,p;k=1,2,…,nq;/>为q类中nq个样本中第i个变量的均值,m为该类样本数据主成分数,/>为训练集的光谱残差,b和cαk是为了使光谱残差达到极小值的参数;模型建立后,计算该q类样本的总体标准偏差,采用多元线性回归法拟合,再判别未知样本u是否属于q类,区分此模型与原训练集构建模型的不同,以表示拟合后的光谱残差,/>为该类中对应的元素,其中i=1,2,…,p;u=1,2,…,nq为q类中nq个样本中第i个变量的均值,b和Cαk是为了使光谱残差达到极小值的参数;
随后,以主成分模型和拟合后的模型的标准差比值为检测计量M,当M≥0.8时样本u属于q类;
CN202210531991.0A 2022-05-17 2022-05-17 基于反射光谱的激光器参数自动调控方法 Active CN114850658B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210531991.0A CN114850658B (zh) 2022-05-17 2022-05-17 基于反射光谱的激光器参数自动调控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210531991.0A CN114850658B (zh) 2022-05-17 2022-05-17 基于反射光谱的激光器参数自动调控方法

Publications (2)

Publication Number Publication Date
CN114850658A CN114850658A (zh) 2022-08-05
CN114850658B true CN114850658B (zh) 2024-05-24

Family

ID=82638083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210531991.0A Active CN114850658B (zh) 2022-05-17 2022-05-17 基于反射光谱的激光器参数自动调控方法

Country Status (1)

Country Link
CN (1) CN114850658B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6700661B1 (en) * 1999-10-14 2004-03-02 Cme Telemetrix, Inc. Method of optimizing wavelength calibration
CN106248652A (zh) * 2016-07-06 2016-12-21 湖南大学 基于工艺参数修正的svr对激光金属沉积中元素浓度实时监测方法
CN108051410A (zh) * 2017-10-31 2018-05-18 浙江农林大学 一种基于激光诱导击穿光谱技术的烟草根部镉含量快速检测方法
CN113155809A (zh) * 2021-03-12 2021-07-23 北京理工大学 一种矿石分类与实时定量分析的光谱检测新方法
CN113295674A (zh) * 2021-04-29 2021-08-24 中国科学院沈阳自动化研究所 一种基于s变换的激光诱导击穿光谱特征非线性处理方法
CN113751429A (zh) * 2021-08-11 2021-12-07 北京工业大学 一种机械臂辅助激光加工的高精度定位及运动补偿方法
CN114252430A (zh) * 2020-09-21 2022-03-29 中国科学院沈阳自动化研究所 一种基于激光诱导击穿光谱技术的矿浆品位在线检测方法
CN114393313A (zh) * 2021-12-01 2022-04-26 华中科技大学 一种基于同轴检测的激光扫描智能加工装置及方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8655807B2 (en) * 2010-04-05 2014-02-18 Applied Research Associates, Inc. Methods for forming recognition algorithms for laser-induced breakdown spectroscopy
BR112017003169A2 (pt) * 2014-08-19 2018-07-03 Ingrain Inc método e sistema para obter informação geoquímica a partir de pirólise induzida por espectroscopia de degradação induzida por laser
US10232473B2 (en) * 2016-02-26 2019-03-19 General Electric Company System and method for performing laser induced breakdown spectroscopy during laser ablation coating removal

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6700661B1 (en) * 1999-10-14 2004-03-02 Cme Telemetrix, Inc. Method of optimizing wavelength calibration
CN106248652A (zh) * 2016-07-06 2016-12-21 湖南大学 基于工艺参数修正的svr对激光金属沉积中元素浓度实时监测方法
CN108051410A (zh) * 2017-10-31 2018-05-18 浙江农林大学 一种基于激光诱导击穿光谱技术的烟草根部镉含量快速检测方法
CN114252430A (zh) * 2020-09-21 2022-03-29 中国科学院沈阳自动化研究所 一种基于激光诱导击穿光谱技术的矿浆品位在线检测方法
CN113155809A (zh) * 2021-03-12 2021-07-23 北京理工大学 一种矿石分类与实时定量分析的光谱检测新方法
CN113295674A (zh) * 2021-04-29 2021-08-24 中国科学院沈阳自动化研究所 一种基于s变换的激光诱导击穿光谱特征非线性处理方法
CN113751429A (zh) * 2021-08-11 2021-12-07 北京工业大学 一种机械臂辅助激光加工的高精度定位及运动补偿方法
CN114393313A (zh) * 2021-12-01 2022-04-26 华中科技大学 一种基于同轴检测的激光扫描智能加工装置及方法

Also Published As

Publication number Publication date
CN114850658A (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
RU2368931C2 (ru) Способ для управления качеством промышленных процессов, в частности лазерных сварочных процессов
CN102147602B (zh) 用于监控激光加工工艺的质量的方法及相应的系统
CN110653223B (zh) 一种激光清洗监测装置及其监测方法和激光清洗机
CN107025651B (zh) 激光清洗能量的确定方法及装置
WO2009062887A1 (de) Verfahren und vorrichtung zum ultraschallbonden
CN104374752A (zh) 基于共线激光诱导击穿光谱的作物营养元素快速检测方法
CN110349153A (zh) 基于多特征融合的图像法烟草含梗率检测方法
CN107597737A (zh) 一种激光清洗过程实时监测方法
CN109142251B (zh) 随机森林辅助人工神经网络的libs定量分析方法
CN113049507A (zh) 多模型融合的光谱波长选择方法
DE102019121446A1 (de) Quantitatives Charakterisierungsverfahren für die Fläche und für den Inhalt unterschiedlicher Typen von Einschlüssen in Stahl
CN114850658B (zh) 基于反射光谱的激光器参数自动调控方法
CN112200000B (zh) 一种焊接稳定性识别模型训练方法、焊接稳定性识别方法
CN114486687B (zh) 飞秒激光加工细胞的多尺度连续观测反馈方法及装置
CN112115856A (zh) 脑电质量评估方法、存储介质及系统
DE102005046129A1 (de) Vorrichtung, System und Verfahren zur Steuerung und Überwachung der Energie eines Lasers
CN109164069B (zh) 一种果树叶面病害级别的鉴定方法
CN117102630B (zh) 电弧焊焊接质量监测方法、装置、电子设备和存储介质
CN117852949A (zh) 一种烟叶烘烤自动评价考核方法
CN114504116A (zh) 基于激光诱导击穿光谱的烟叶分类设备及方法
CN100464270C (zh) 控制工业生产过程质量的方法及其系统
CN110517226A (zh) 基于双边滤波的多特征纹理图像融合的烟梗区域提取方法
CN106770152A (zh) 一种基于特征峰及算法参数选取的奶粉品牌快速鉴别方法
CN110000162A (zh) 利用振动信号判定激光除漆效果的方法
CN112611745A (zh) 一种水果可溶性固形物间接无损检测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant