CN114235141A - 基于mems传感器的振动监测方法及装置 - Google Patents

基于mems传感器的振动监测方法及装置 Download PDF

Info

Publication number
CN114235141A
CN114235141A CN202111491956.2A CN202111491956A CN114235141A CN 114235141 A CN114235141 A CN 114235141A CN 202111491956 A CN202111491956 A CN 202111491956A CN 114235141 A CN114235141 A CN 114235141A
Authority
CN
China
Prior art keywords
vibration
data
equipment
speed
acceleration
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.)
Pending
Application number
CN202111491956.2A
Other languages
English (en)
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.)
Sichuan Cric Technology Co ltd
Original Assignee
Sichuan Cric Technology Co ltd
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 Sichuan Cric Technology Co ltd filed Critical Sichuan Cric Technology Co ltd
Priority to CN202111491956.2A priority Critical patent/CN114235141A/zh
Publication of CN114235141A publication Critical patent/CN114235141A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及振动监测技术领域,公开了一种基于MEMS传感器的振动监测方法及装置,以提高振动监测的响应速度、稳定性和准确性,方案主要包括:接收MEMS传感器采集并发送的待测设备的振动加速度数据;对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。本发明适用于机械设备。

Description

基于MEMS传感器的振动监测方法及装置
技术领域
本发明涉及振动监测技术领域,具体来说涉及一种基于MEMS传感器的振动监测方法及装置。
背景技术
各种包括旋转部件的机械设备因设计制造、安装调试、运行维修中的一些缺陷和环境影响,在运行时均会产生机械振动,而振动的存在必然引起结构损伤及材料疲劳,并在相当短的时间迅速发展扩大,进而引起设备疲劳破坏并可能导致重大事故。因此通过监测机械设备的振动状态,对机械设备健康进行预测并对潜在的故障进行早期预警,以避免机械设备突发故障而造成重大损失,是智慧工业的关键环节。
机械设备健康预测主要的技术方法是通过监测设备运行中的振动情况,进而分析设备是否存在失衡、失准或其他可能导致设备故障的机械异常。传统的振动监测技术是基于压电加速度传感器,但该技术存在振动响应慢、稳定性较低、成本高,并且由于监测装置体积大,不便于将其安装在机械设备的关键监测点位,造成检测准确度较低。
发明内容
本发明旨在提出一种基于MEMS传感器的振动监测方法及装置,以提高振动监测的响应速度、稳定性和准确性。
本发明解决上述技术问题所采用的技术方案是:
第一方面,提供一种基于MEMS传感器的振动监测方法,包括以下步骤:
步骤1、接收MEMS传感器采集并发送的待测设备的振动加速度数据;
步骤2、对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;
步骤3、计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。
作为进一步优化,该方法还包括:
步骤4、对所述振动加速度进行数字积分运算,得到对应的振动速度数据;
步骤5、对所述振动速度数据进行数字高通滤波运算,过滤所述振动速度数据的直流分量和速度线性漂移,得到振动速度;
步骤6、计算所述振动速度的第二均方根值,并根据所述第二均方根值获得待测设备的载荷循环速度和机械形变情况。
作为进一步优化,该方法还包括:
步骤7、对所述振动速度进行数字积分运算,得到对应的振动位移数据;
步骤8、对所述振动位移数据进行数字高通滤波运算,过滤所述振动位移数据的直流分量和速度线性漂移,得到振动位移;
步骤9、计算所述振动位移的第三均方根值,并根据所述第三均方根值获得待测设备的旋转部件的位置精度和形变引起的失准数据和失衡数据。
作为进一步优化,该方法还包括:
分别将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度以及形变引起的失准数据和失衡数据与正常状态下的对应预设阈值进行比较,并根据比较结果确定待测设备的后续运行状态。
作为进一步优化,所述待测设备的后续运行状态至少包括:正常运行状态、需监护运行的预警状态、需停机检修状态和需停机状态。
作为进一步优化,该方法还包括:
将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度、形变引起的失准数据和失衡数据以及确定出的待测设备的后续运行状态发送至后台服务器。
第二方面,提供一种基于MEMS传感器的振动监测装置,包括:
MEMS传感器,用于采集并发送待测设备的振动加速度数据;
微处理器,用于对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。
作为进一步优化,所述微处理器,还用于对所述振动加速度进行数字积分运算,得到对应的振动速度数据;对所述振动速度数据进行数字高通滤波运算,过滤所述振动速度数据的直流分量和速度线性漂移,得到振动速度;以及计算所述振动速度的第二均方根值,并根据所述第二均方根值获得待测设备的载荷循环速度和机械形变情况。
作为进一步优化,所述微处理器,还用于对所述振动速度进行数字积分运算,得到对应的振动位移数据;对所述振动位移数据进行数字高通滤波运算,过滤所述振动位移数据的直流分量和速度线性漂移,得到振动位移;以及计算所述振动位移的第三均方根值,并根据所述第三均方根值获得待测设备的旋转部件的位置精度和形变引起的失准数据和失衡数据。
作为进一步优化,所述微处理器,还用于分别将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度以及形变引起的失准数据和失衡数据与正常状态下的对应预设阈值进行比较,并根据比较结果确定待测设备的后续运行状态。
本发明的有益效果是:本发明所述的基于MEMS传感器的振动监测方法及装置,因MEMS传感器具有体积小、功耗低、响应快、稳定性高、多轴向检测的特点,同时本发明根据检测到的振动加速度数据,计算分析了能表征设备机械疲劳损伤的加速度、速度、位移以及频谱等多维度的特征参数,以作为设备运行健康状态的预测诊断标准,不仅易于嵌入设备的任何待测部位和长时间实时监测,而且能对设备因旋转部件位置偏移、机械形变、惯性力、超载荷等多种因素导致的健康异常进行早期预测预警,从而避免设备突发故障而造成重大损失。本发明具有振动恢复快、频率响应宽、稳定性高及低功耗、小体积、低成本等优点,易于实施,具有很好的实际应用价值。
附图说明
图1为本发明实施例所述的基于MEMS传感器的振动监测装置的结构示意图;
图2为本发明实施例所述的基于MEMS传感器的振动监测方法的流程示意图。
具体实施方式
下面将结合附图对本发明的实施方式进行详细描述。
本发明提供的基于MEMS传感器的振动监测方法及装置,主要方案包括:接收MEMS传感器采集并发送的待测设备的振动加速度数据;对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。
MEMS即微电子机械系统(Microelectro Mechanical Systems),MEMS传感器是采用微电子和微机械加工技术制造出来的新型传感器。与传统的传感器相比,它具有体积小、重量轻、成本低、功耗低、可靠性高、适于批量化生产、易于集成和实现智能化的特点。本发明基于MEMS传感器的振动监测装置对待测设备进行实时振动监测,具体流程如下:首先将MEMS传感器固定在待测设备的振动监测点位,在获取MEMS传感器采集的待测设备的振动加速度数据之后,对振动加速度数据进行数字高通滤波运算,以过滤所述振动加速度数据中的重力加速度和静态直流分量得到振动加速度,然后再根据振动速度对应的第一均方根值以及离散傅里叶变换结果获得待测设备的振动频率、振动强度以及振动惯性力对待测设备造成的影响。
实施例
本发明实施例所述的基于MEMS传感器的振动监测装置,如图1所示,包括:锂离子电池1、电源管理模块2、加速度传感器3、微处理器4和通讯模块5。锂离子电池1与电源管理模块2的输入端相连,为监测装置提供电源,同时电源管理模块2对锂离子电池1进行充放电管理;电源管理模块2的输出端与加速度传感器3、微处理器4、通讯模块5的电源输入端相连,锂离子电池的电压经电源管理模块转换为所需要的供电电压输出;加速度传感器3用于检测设备的机械振动,输出加速度数据至微处理器4;微处理器4对加速度数据进行运算分析并对设备的运行健康状态进行预测诊断,分析诊断的结果输出至通讯模块5;通讯模块接收上位机或远程监控中心的指令启动或停止振动监测装置,并将监测分析的数据发送给上位机或远程监控中心。
本实施例中,所述加速度传感器为MEMS加速度传感器,且为三轴传感器,可以同时监测XYZ三个方向的振动加速度数据,比如意法半导体ST公司的IIS2DH。
所述微处理器优选支持浮点运算的32位MCU芯片,比如基于ARM Cortex-M4内核的STM32微处理器。
所述通讯模块可以是有线或无线通讯技术中的任一种,比如有线的RS485,无线的蓝牙或WIFI或4G通讯模块等。
所述振动监测装置通过外壳上的螺钉固定在待测设备的振动监测点位,也可以在振动监测装置的外壳上嵌入永久磁铁,使振动监测装置吸附在待测设备的振动监测点位。
基于上述振动监测装置,本实施例所述的基于MEMS传感器的振动监测方法,由于XYZ三个方向的加速度数据的运算处理和预测诊断的方法相同,为简化阐述,以下不再区分XYZ方向。
本实施例将所述振动监测装置固定在待测设备振动监测点位后,可以实时对待测设备的振动情况进行监测,也可以在接收到后台服务器的指令启动监测时,MEMS加速度传感器采集振动加速度数据并输入微处理器,微处理器对加速度数据进行运算处理和预测诊断,如图2所示,具体包括以下步骤:
步骤1、接收MEMS传感器采集并发送的待测设备的振动加速度数据;
步骤2、对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;
步骤3、计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响;
步骤4、对所述振动加速度进行数字积分运算,得到对应的振动速度数据;
步骤5、对所述振动速度数据进行数字高通滤波运算,过滤所述振动速度数据的直流分量和速度线性漂移,得到振动速度;
步骤6、计算所述振动速度的第二均方根值,并根据所述第二均方根值获得待测设备的载荷循环速度和机械形变情况;
步骤7、对所述振动速度进行数字积分运算,得到对应的振动位移数据;
步骤8、对所述振动位移数据进行数字高通滤波运算,过滤所述振动位移数据的直流分量和速度线性漂移,得到振动位移;
步骤9、计算所述振动位移的第三均方根值,并根据所述第三均方根值获得待测设备的旋转部件的位置精度和形变引起的失准数据和失衡数据。
步骤10、分别将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度以及形变引起的失准数据和失衡数据与正常状态下的对应预设阈值进行比较,,根据监测计算值是否超出预设健康阈值及超出的比例判定待测设备的后续运行状态,以便于后续对待测设备的操作,例如正常状态、需监护运行的预警状态、需停机检修状态、需停机禁止运行的状态等。
步骤11、将振动监测参数及健康预测诊断结果上传至后台服务器,等待下一次监测指令。
特别地,作为振动预测诊断判定标准的预设健康阈值,可以来自设备厂商,也可以根据对设备长期使用、观察、维修与测试后的智能学习总结。
本领域的普通技术人员应理解,本实施例中所涉及的数字高通滤波运算、数字积分运算、均方根计算、离散傅里叶变换均是基础的数学计算方法,在此不再赘述。
综上所述,本实施例所述的基于MEMS加速度传感器的设备振动监测方法及装置,因MEMS传感器具有体积小、功耗低、响应快、稳定性高、多轴向检测的特点,同时本实施例根据检测到的振动加速度数据,计算分析了能表征设备机械疲劳损伤的加速度、速度、位移以及频谱等多维度的特征参数,以作为设备运行健康状态的预测诊断标准,因此不仅所述装置易于嵌入设备的任何待测部位和长时间实时监测,而且能对设备因旋转部件位置偏移、机械形变、惯性力、超载荷等多种因素导致的健康异常进行早期预测预警,从而避免设备突发故障而造成重大损失,本实施例的技术方案,在智慧工业的设备运维中具有很好的实际应用推广价值。

Claims (10)

1.基于MEMS传感器的振动监测方法,其特征在于,包括以下步骤:
步骤1、接收MEMS传感器采集并发送的待测设备的振动加速度数据;
步骤2、对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;
步骤3、计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。
2.如权利要求1所述的基于MEMS传感器的振动监测方法,其特征在于,该方法还包括:
步骤4、对所述振动加速度进行数字积分运算,得到对应的振动速度数据;
步骤5、对所述振动速度数据进行数字高通滤波运算,过滤所述振动速度数据的直流分量和速度线性漂移,得到振动速度;
步骤6、计算所述振动速度的第二均方根值,并根据所述第二均方根值获得待测设备的载荷循环速度和机械形变情况。
3.如权利要求2所述的基于MEMS传感器的振动监测方法,其特征在于,该方法还包括:
步骤7、对所述振动速度进行数字积分运算,得到对应的振动位移数据;
步骤8、对所述振动位移数据进行数字高通滤波运算,过滤所述振动位移数据的直流分量和速度线性漂移,得到振动位移;
步骤9、计算所述振动位移的第三均方根值,并根据所述第三均方根值获得待测设备的旋转部件的位置精度和形变引起的失准数据和失衡数据。
4.如权利要求3所述的基于MEMS传感器的振动监测方法,其特征在于,该方法还包括:
分别将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度以及形变引起的失准数据和失衡数据与正常状态下的对应预设阈值进行比较,并根据比较结果确定待测设备的后续运行状态。
5.如权利要求4所述的基于MEMS传感器的振动监测方法,其特征在于,所述待测设备的后续运行状态至少包括:正常运行状态、需监护运行的预警状态、需停机检修状态和需停机状态。
6.如权利要求4所述的基于MEMS传感器的振动监测方法,其特征在于,该方法还包括:
将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度、形变引起的失准数据和失衡数据以及确定出的待测设备的后续运行状态发送至后台服务器。
7.基于MEMS传感器的振动监测装置,其特征在于,包括:
MEMS传感器,用于采集并发送待测设备的振动加速度数据;
微处理器,用于对所述振动加速度数据进行数字高通滤波运算,过滤所述振动加速度数据中的重力加速度和静态直流分量,得到振动加速度;计算所述振动加速度的第一均方根值,并对所述振动加速度进行离散傅里叶变换,根据所述第一均方根值获得待测设备的振动频率和振动强度,根据所述离散傅里叶变换结果获得振动惯性力对待测设备造成的影响。
8.如权利要求7所述的基于MEMS传感器的振动监测装置,其特征在于,所述微处理器,还用于对所述振动加速度进行数字积分运算,得到对应的振动速度数据;对所述振动速度数据进行数字高通滤波运算,过滤所述振动速度数据的直流分量和速度线性漂移,得到振动速度;以及计算所述振动速度的第二均方根值,并根据所述第二均方根值获得待测设备的载荷循环速度和机械形变情况。
9.如权利要求8所述的基于MEMS传感器的振动监测装置,其特征在于,所述微处理器,还用于对所述振动速度进行数字积分运算,得到对应的振动位移数据;对所述振动位移数据进行数字高通滤波运算,过滤所述振动位移数据的直流分量和速度线性漂移,得到振动位移;以及计算所述振动位移的第三均方根值,并根据所述第三均方根值获得待测设备的旋转部件的位置精度和形变引起的失准数据和失衡数据。
10.如权利要求8所述的基于MEMS传感器的振动监测装置,其特征在于,所述微处理器,还用于分别将所述待测设备的振动频率、振动强度、振动惯性力对待测设备造成的影响、载荷循环速度、机械形变情况、旋转部件的位置精度以及形变引起的失准数据和失衡数据与正常状态下的对应预设阈值进行比较,并根据比较结果确定待测设备的后续运行状态。
CN202111491956.2A 2021-12-08 2021-12-08 基于mems传感器的振动监测方法及装置 Pending CN114235141A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111491956.2A CN114235141A (zh) 2021-12-08 2021-12-08 基于mems传感器的振动监测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111491956.2A CN114235141A (zh) 2021-12-08 2021-12-08 基于mems传感器的振动监测方法及装置

Publications (1)

Publication Number Publication Date
CN114235141A true CN114235141A (zh) 2022-03-25

Family

ID=80753938

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111491956.2A Pending CN114235141A (zh) 2021-12-08 2021-12-08 基于mems传感器的振动监测方法及装置

Country Status (1)

Country Link
CN (1) CN114235141A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114766116A (zh) * 2022-05-27 2022-07-22 四川启睿克科技有限公司 基于mems传感器的物联网土壤修复系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104236877A (zh) * 2014-08-13 2014-12-24 北京航空航天大学 一种航空液压滑阀磨损寿命计算方法
CN106441899A (zh) * 2016-11-18 2017-02-22 上海卫星工程研究所 基于振动信号检测卫星动量轮轴承损伤诊断及筛选方法
CN108680244A (zh) * 2018-04-26 2018-10-19 浙江大学 一种旋转机械振动无线监测装置及方法
JP2021076619A (ja) * 2021-02-19 2021-05-20 住友電気工業株式会社 振動センサシステム、振動計測方法および振動計測プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104236877A (zh) * 2014-08-13 2014-12-24 北京航空航天大学 一种航空液压滑阀磨损寿命计算方法
CN106441899A (zh) * 2016-11-18 2017-02-22 上海卫星工程研究所 基于振动信号检测卫星动量轮轴承损伤诊断及筛选方法
CN108680244A (zh) * 2018-04-26 2018-10-19 浙江大学 一种旋转机械振动无线监测装置及方法
JP2021076619A (ja) * 2021-02-19 2021-05-20 住友電気工業株式会社 振動センサシステム、振動計測方法および振動計測プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
汪永华 等: "《机电设备故障诊断与维修》", vol. 2, 机械工业出版社, pages: 55 - 56 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114766116A (zh) * 2022-05-27 2022-07-22 四川启睿克科技有限公司 基于mems传感器的物联网土壤修复系统

Similar Documents

Publication Publication Date Title
CN107247230B (zh) 一种基于支持向量机和数据驱动的旋转电机状态监测方法
CN105352588B (zh) 一种无刷直流电机振动检测系统的设计方法
CN108376184A (zh) 一种桥梁健康监控的方法及系统
TWI732324B (zh) 工具機健康狀態監測方法
CN107144303A (zh) 一种矿用机电设备的故障检测系统、方法及装置
CN109738181B (zh) 具有机械卡涩检测功能的隔离开关操作机构及其应用方法
CN106706241B (zh) 一种风力机叶片损伤主动自检装置及方法
CN106482828A (zh) 一种航空发动机振动故障的检测诊断装置及方法
US20160091394A1 (en) Techniques for monitoring gear condition
CN103884482A (zh) 一种基于压缩机的振动测试方法和系统
CN111351567A (zh) 一种基于加速度传感器的船舶设备振动状态监测装置及方法
CN108918118A (zh) 一种基于人工激励的输电铁塔螺栓松动监测系统及方法
CN114235141A (zh) 基于mems传感器的振动监测方法及装置
CN105372537B (zh) 基于Labview/PDA的变频器智能测试平台
CN203396914U (zh) 电机运行状态监测和故障检测系统
CN108825447B (zh) 一种风力机监测方法及系统
CN109142961A (zh) 一种电连接器间歇失效机理研究方法
KR20210006832A (ko) 기계고장 진단 방법 및 장치
CN112901546B (zh) 电控硅油风扇性能间接测量方法及装置
JP2018207650A (ja) 回転電機の特徴量評価システムおよび回転電機の特徴量評価方法
TW202004010A (zh) 風力發電系統
CN113485898A (zh) 振动测点显示方法、装置、设备及存储介质
CN112524076A (zh) 磁悬浮离心式压缩机故障测试装置及测试方法
CN202091172U (zh) 气体压缩机状态监测和可靠性分析装置
CN111722106A (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20220325