CN109829252B - 一种基于影响线识别的桥梁状况快速评级方法 - Google Patents

一种基于影响线识别的桥梁状况快速评级方法 Download PDF

Info

Publication number
CN109829252B
CN109829252B CN201910198942.8A CN201910198942A CN109829252B CN 109829252 B CN109829252 B CN 109829252B CN 201910198942 A CN201910198942 A CN 201910198942A CN 109829252 B CN109829252 B CN 109829252B
Authority
CN
China
Prior art keywords
bridge
influence line
influence
theoretical
rating
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.)
Expired - Fee Related
Application number
CN201910198942.8A
Other languages
English (en)
Other versions
CN109829252A (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.)
FUZHOU HIGHWAY BUREAU
JSTI Group Co Ltd
Original Assignee
FUZHOU HIGHWAY BUREAU
JSTI Group 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 FUZHOU HIGHWAY BUREAU, JSTI Group Co Ltd filed Critical FUZHOU HIGHWAY BUREAU
Priority to CN201910198942.8A priority Critical patent/CN109829252B/zh
Publication of CN109829252A publication Critical patent/CN109829252A/zh
Application granted granted Critical
Publication of CN109829252B publication Critical patent/CN109829252B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种基于影响线识别的桥梁状况快速评级方法,属于桥梁检测领域,旨在提供一种有能够准确高效地识别桥梁位移影响的桥梁状况快速评级方法,其技术方案要点如下,它包括如下步骤:S1:计算理论影响线,通过有限元瞬态动力学求解得到;S2:获取实测影响线,实时获取跑车全过程挠度时程数据,测试截面为跨中位置,截面横向布置2‑3个测点;S3:提取静载影响线,将各频带上的噪声对应的小波系数去除,保留原始信号的小波分解系数,对处理后的系数进行小波系数重构,得到去燥后的信号;S4:构建校验系数指标δ,根据校验系数指标δ进行桥梁状况评级。本发明适用于桥梁状况的快速评级。

Description

一种基于影响线识别的桥梁状况快速评级方法
技术领域
本发明涉及一种桥梁状况快速评级方法,特别涉及一种基于影响线识别的桥梁状况快速评级方法。
背景技术
影响线是桥梁的固有特性,具有物理意义明确、对损伤敏感、对环境变化较不明显的特点,包含丰富的桥梁局部信息,新桥的成桥实验和旧桥的检测评估中常包含测量影响线。众所周知,一座桥梁服役了一段时间后,结构就有可能会发生劣化,随之出现局部刚度下降或其他固有特性的改变。基于影响线实测结果可以实现桥梁承载力的快速检测和对于桥梁的健康状况进行快速评级。
针对桥梁影响线的测定方法,目前,国内外学者开展了大量的研究工作和广泛的探索,但仍旧缺乏一种普遍认可的能够准确高效地识别桥梁位移影响线的方法。与此同时利用影响线进行桥梁健康状况快速评级方面的研究做的也还不是很多。
因此,针对上述现有技术的不足,桥梁影响线测定方法需要从准确性和工程实用性两个方面进一步深入研究,并且需要构建一个新的指标从而达到对桥梁的健康状况进行快速评级。
本设计人基于从事此类产品工程应用多年丰富的实务经验及专业知识,并配合学理的运用,积极加以研究创新,以期创设一种基于影响线识别的桥梁状况快速评级方法,使其更具有实用性。
发明内容
本发明的目的是提供一种基于影响线识别的桥梁状况快速评级方法,具有能够准确高效地识别桥梁位移影响的优点。
本发明的上述技术目的是通过以下技术方案得以实现的:
一种基于影响线识别的桥梁状况快速评级方法,包括如下步骤:
S1:计算理论影响线,通过有限元瞬态动力学求解得到,具体为,将结构沿跨径方向离散为若干节点,通过改变荷载作用的节点位置实现对移动荷载的模拟施加,每一荷载步内打开瞬态效应,循环加载求解的结果即为理论影响线;
S2:获取实测影响线,实时获取跑车全过程挠度时程数据,测试截面为跨中位置,截面横向布置2-3个测点;
S3:提取静载影响线,将各频带上的噪声对应的小波系数去除,保留原始信号的小波分解系数,对处理后的系数进行小波系数重构,得到去燥后的信号;
S4:构建校验系数指标δ,根据校验系数指标δ进行桥梁状况评级;
其中,
Figure BDA0001996735690000021
Figure BDA0001996735690000022
式中,w为静载影响线与理论影响线的差值指标,w为理论影响线指标。
进一步的,理论影响线的校验方法,具体为,通过车桥两子系统间的分离迭代实现,其中,车桥运动方程为:
Figure BDA0001996735690000031
Figure BDA0001996735690000032
式中:Mb为桥梁的质量(kg),ub为桥梁节点的位移(m),Cb为桥梁的阻尼[N/(m/s)],Kb为桥梁的刚度(N/m),Fgb作用在桥梁的与桥梁车辆运动无关的荷载(N),Fvb桥车间的相互作用力(N);
Figure BDA0001996735690000033
式中,Mv为车辆的质量(kg),uv为车辆节点的位移(m),Cv为车辆的阻尼[N/(m/s)],Kv为车辆的刚度(N/m),Fvg作用在车辆的与桥梁车辆运动无关的荷载(N),Fbv桥车间的相互作用力(N)。
进一步的,提取静载影响线,利用MATLAB进行计算,剔除动载效应,具体为,基于两尺度方程系数{hk}设计离散正交小波dbN;
其中,N为小波阶数,Ψ小波和尺度函数域中的支撑区为2N-1,
Figure BDA0001996735690000034
的消失矩为N,dbN不具对称性,且dbN无显式表达式;
其中,{hk}的传递函数的模的平方有显式表达式;
Figure BDA0001996735690000035
Figure BDA0001996735690000036
式中,
Figure BDA0001996735690000037
为二项式的系数;
Figure BDA0001996735690000038
Figure BDA0001996735690000041
Figure BDA0001996735690000042
进一步的,理论影响线指标w的计算具体为,
Figure BDA0001996735690000043
式中,n个点将桥梁划分为n-1等分,αi表示桥梁健康临界状态下i点处的跨中挠度影响线值,ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure BDA0001996735690000044
时,
Figure BDA0001996735690000045
Figure BDA0001996735690000046
Figure BDA0001996735690000047
时,
Figure BDA0001996735690000048
Figure BDA0001996735690000049
进一步的,静载影响线与理论影响线的差值指标w的计算具体为,
Figure BDA00019967356900000410
式中,n个点将桥梁划分为n-1等分,βi为实际情况下桥梁i点处的跨中挠度影响线值,
Figure BDA00019967356900000411
ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure BDA00019967356900000412
时,
Figure BDA00019967356900000413
Figure BDA0001996735690000051
时,
Figure BDA0001996735690000052
进一步的,当δ>0.6时,说明桥梁结构工作性能良好,为1类桥;
当0.5<δ≤0.6时,说明桥梁结构工作性能低于健康临界值,为2类桥;
当0.4<δ≤0.5时,说明桥梁结构工作性能低于健康临界值,为3类桥;
当0.3<δ≤0.4时,说明桥梁结构工作性能低于健康临界值,为4类桥;
当δ≤0.3时,说明桥梁结构工作性能已超过健康临界值,为5类桥。
本发明具有以下有益效果:
首先,计算理论影响线,然后对实测影响线进行去躁,剔除实测影响线中动载效应引起的局部波动,提取静载影响线,静载影响线为去躁后的实测影响线,最后构建校验系数指标并基于该指标对桥梁健康状况进行评级。本发明能够很快的提取影响线对桥梁状况进行评级,为准确高效地识别桥梁位移影响线提供了一种新的思路。
附图说明
图1是本实施例中用于体现基于影响线识别的桥梁状况快速评级方法的流程图;
图2是本实施例中用于体现理论影响线、实测影响线以及静载影响线识别的效果图;
图3是本实施例中用于体现校验系数指标与桥梁状况评级的对应关系图。
图中,1、实测影响线;2、静载影响线;3、理论影响线。
具体实施方式
以下结合附图对本发明作进一步详细说明。
其中相同的零部件用相同的附图标记表示。需要说明的是,下面描述中使用的词语“前”、“后”、“左”、“右”、“上”和“下”指的是附图中的方向,词语“底面”和“顶面”、“内”和“外”分别指的是朝向或远离特定部件几何中心的方向。
一种基于影响线识别的桥梁状况快速评级方法,如图1所示,包括如下步骤:
考虑动载效应的理论影响线计算可看作是桥梁结构在车辆荷载下的强迫振动问题,其特点为荷载是移动的,且车辆本身也是一个具有质量的振动体系,这使得车桥耦合系统的动力特征随荷载位置的移动而不断变化。
S1:计算理论影响线,通过有限元(ANSYS)瞬态动力学求解得到,具体为,将结构沿跨径方向离散为若干节点,通过改变荷载作用的节点位置实现对移动荷载的模拟施加,每一荷载步内若按静力求解所得的跨中节点位移时程曲线即为标准静力影响线;每一荷载步内打开瞬态效应,循环加载求解的结果即为理论影响线。最终理论影响线的获取通过上述有限元瞬态动力学求解方法得到,并通过下述理论方法校验。
理论影响线的校验方法,具体为,通过车桥两子系统间的分离迭代实现,其中,车桥运动方程为:
Figure BDA0001996735690000061
式中:Mb为桥梁的质量(kg),ub为桥梁节点的位移(m),Cb为桥梁的阻尼[N/(m/s)],Kb为桥梁的刚度(N/m),Fgb作用在桥梁的与桥梁车辆运动无关的荷载(N),Fvb桥车间的相互作用力(N);
Figure BDA0001996735690000071
式中,Mv为车辆的质量(kg),uv为车辆节点的位移(m),Cv为车辆的阻尼[N/(m/s)],Kv为车辆的刚度(N/m),Fvg作用在车辆的与桥梁车辆运动无关的荷载(N),Fbv桥车间的相互作用力(N)。
对于车桥系统,将桥梁与车辆运动方程分别独立求解,通过分离迭代来满足车桥两子系统间的几何、力学耦合关系。
S2:获取实测影响线,实时获取跑车全过程挠度时程数据,测试截面为跨中位置,截面横向布置2-3个测点;具体的,实测影响线通过现场定速跑车试验获取。车辆轴重与轴距与数值计算的加载车辆保持一致,匀速跑车。使用雷达挠度测量仪器获取跑车全过程挠度时程数据。
S3:提取静载影响线,根据噪声与信号在不同频带上的小波分解系数具有不同强度分布特点,将各频带上的噪声对应的小波系数去除,保留原始信号的小波分解系数,对处理后的系数进行小波系数重构,得到去燥后的信号。
基于dbN方法,利用MATLAB进行计算,剔除动载效应,即可得到提取静载影响线,具体为,基于两尺度方程系数{hk}设计离散正交小波dbN;
其中,N为小波阶数,Ψ小波和尺度函数域中的支撑区为2N-1,
Figure BDA0001996735690000072
的消失矩为N,dbN不具对称性,且dbN无显式表达式;
其中,{hk}的传递函数的模的平方有显式表达式;
Figure BDA0001996735690000073
式中,
Figure BDA0001996735690000081
为二项式的系数;
Figure BDA0001996735690000082
Figure BDA0001996735690000083
S4:构建校验系数指标δ,根据校验系数指标δ进行桥梁状况评级;其中,
Figure BDA0001996735690000084
式中,w为静载影响线与理论影响线的差值指标,w为理论影响线指标。
计算理论影响线指标,即在健康临界状态下桥梁的跨中挠度影响线,从桥端到桥尾依次取点1,2,…,n(桥端点为1号点,桥尾点为n号点),使这n个点将桥梁划分为n-1等分;
理论影响线指标w的计算具体为,
Figure BDA0001996735690000085
式中,n个点将桥梁划分为n-1等分,αi表示桥梁健康临界状态下i点处的跨中挠度影响线值,ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure BDA0001996735690000086
时,
Figure BDA0001996735690000087
Figure BDA0001996735690000088
时,
Figure BDA0001996735690000089
静载影响线与理论影响线的差值指标w的计算具体为,
Figure BDA00019967356900000810
式中,n个点将桥梁划分为n-1等分,βi为实际情况下桥梁i点处的跨中挠度影响线值,
Figure BDA0001996735690000091
ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure BDA0001996735690000092
时,
Figure BDA0001996735690000093
Figure BDA0001996735690000094
时,
Figure BDA0001996735690000095
根据校验系数指标δ进行桥梁状况评级,具体的,结合桥梁评定规范中对于桥梁状况的评估方法,校验系数指标δ是评定结构工作状况的重要指标。随着校验系数指标δ的增大,桥梁的损伤状况越小,反映桥梁健康状况越好。参照《公路桥梁技术状况评定标准》,以校验系数为指标对一般对于钢筋混凝土梁桥有如图3评定等级划分,评定等级分1~5类技术状况的桥梁,级别越高表示结构损伤越来越大。
本实施例中,当δ>0.6时,说明桥梁结构工作性能良好,为1类桥;
当0.5<δ≤0.6时,说明桥梁结构工作性能低于健康临界值,为2类桥;
当0.4<δ≤0.5时,说明桥梁结构工作性能低于健康临界值,为3类桥;
当0.3<δ≤0.4时,说明桥梁结构工作性能低于健康临界值,为4类桥;
当δ≤0.3时,说明桥梁结构工作性能已超过健康临界值,为5类桥。
首先,计算理论影响线,然后对实测影响线进行去躁,剔除实测影响线中动载效应引起的局部波动,提取静载影响线,静载影响线为去躁后的实测影响线,最后构建校验系数指标并基于该指标对桥梁健康状况进行评级。本发明能够很快的提取影响线对桥梁状况进行评级,为准确高效地识别桥梁位移影响线提供了一种新的思路,利用本发明得到的影响线识别效果图如图2所示。
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。

Claims (4)

1.一种基于影响线识别的桥梁状况快速评级方法,其特征在于,包括如下步骤:
S1:计算理论影响线,通过有限元瞬态动力学求解得到,具体为,将结构沿跨径方向离散为若干节点,通过改变荷载作用的节点位置实现对移动荷载的模拟施加,每一荷载步内打开瞬态效应,循环加载求解的结果即为理论影响线;
S2:获取实测影响线,实时获取跑车全过程挠度时程数据,测试截面为跨中位置,截面横向布置2-3个测点;
S3:提取静载影响线,将各频带上的噪声对应的小波系数去除,保留原始信号的小波分解系数,对处理后的系数进行小波系数重构,得到去躁后的信号;
S4:构建校验系数指标δ,根据校验系数指标δ进行桥梁状况评级;
其中,
Figure FDA0003658400890000011
式中,w为静载影响线与理论影响线的差值指标,w为理论影响线指标;
理论影响线指标w的计算具体为,
Figure FDA0003658400890000012
式中,n个点将桥梁划分为n-1等分,αi表示桥梁健康临界状态下i点处的跨中挠度影响线值,ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure FDA0003658400890000013
时,
Figure FDA0003658400890000014
Figure FDA0003658400890000015
时,
Figure FDA0003658400890000016
静载影响线与理论影响线的差值指标w的计算具体为,
Figure FDA0003658400890000017
式中,n个点将桥梁划分为n-1等分,βi为实际情况下桥梁i点处的跨中挠度影响线值,
Figure FDA0003658400890000021
ηi为i点处的位置加权系数;
其中,ηi求解方法:记跨中处加权系数为1,η1=0,ηn=0;
Figure FDA0003658400890000022
时,
Figure FDA0003658400890000027
Figure FDA0003658400890000023
时,
Figure FDA0003658400890000028
2.根据权利要求1所述的一种基于影响线识别的桥梁状况快速评级方法,其特征在于,理论影响线的校验方法,具体为,通过车桥两子系统间的分离迭代实现,其中,车桥运动方程为:
Figure FDA0003658400890000024
式中:Mb为桥梁的质量(kg),ub为桥梁节点的位移(m),Cb为桥梁的阻尼[N/(m/s)],Kb为桥梁的刚度(N/m),Fbg为作用在桥梁的与桥梁车辆运动无关的荷载(N),Fvb为桥车间的相互作用力(N);
Figure FDA0003658400890000025
式中,Mv为车辆的质量(kg),uv为车辆节点的位移(m),Cv为车辆的阻尼[N/(m/s)],Kv为车辆的刚度(N/m),Fvg为作用在车辆的与桥梁车辆运动无关的荷载(N),Fbv为桥车间的相互作用力(N)。
3.根据权利要求1所述的一种基于影响线识别的桥梁状况快速评级方法,其特征在于,提取静载影响线,利用MATLAB进行计算,剔除动载效应,具体为,基于两尺度方程系数{hk}设计离散正交小波dbN;
其中,N为小波阶数,Ψ小波和尺度函数域中的支撑区为2N-1,
Figure FDA0003658400890000026
的消失矩为N,dbN不具对称性,且dbN无显式表达式;
其中,{hk}的传递函数的模的平方有显式表达式;
Figure FDA0003658400890000031
式中,
Figure FDA0003658400890000032
为二项式的系数;
Figure FDA0003658400890000033
Figure FDA0003658400890000034
4.根据权利要求1所述的一种基于影响线识别的桥梁状况快速评级方法,其特征在于,当δ>0.6时,说明桥梁结构工作性能良好,为1类桥;
当0.5<δ≤0.6时,说明桥梁结构工作性能低于健康临界值,为2类桥;
当0.4<δ≤0.5时,说明桥梁结构工作性能低于健康临界值,为3类桥;
当0.3<δ≤0.4时,说明桥梁结构工作性能低于健康临界值,为4类桥;
当δ≤0.3时,说明桥梁结构工作性能已超过健康临界值,为5类桥。
CN201910198942.8A 2019-03-15 2019-03-15 一种基于影响线识别的桥梁状况快速评级方法 Expired - Fee Related CN109829252B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910198942.8A CN109829252B (zh) 2019-03-15 2019-03-15 一种基于影响线识别的桥梁状况快速评级方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910198942.8A CN109829252B (zh) 2019-03-15 2019-03-15 一种基于影响线识别的桥梁状况快速评级方法

Publications (2)

Publication Number Publication Date
CN109829252A CN109829252A (zh) 2019-05-31
CN109829252B true CN109829252B (zh) 2022-07-19

Family

ID=66870798

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910198942.8A Expired - Fee Related CN109829252B (zh) 2019-03-15 2019-03-15 一种基于影响线识别的桥梁状况快速评级方法

Country Status (1)

Country Link
CN (1) CN109829252B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110261305A (zh) * 2019-06-18 2019-09-20 南京东南建筑机电抗震研究院有限公司 基于影响线的多跨等跨径连续梁桥损伤识别方法
CN111400794A (zh) * 2020-03-10 2020-07-10 上海城建城市运营(集团)有限公司 一种城市高架桥梁总挠度快速智能化监测算法
CN111504245B (zh) * 2020-04-24 2022-01-28 山东高速工程检测有限公司 一种运营期桥梁线形测量方法
CN113627060A (zh) * 2021-08-09 2021-11-09 苏交科集团股份有限公司 一种基于实测挠度影响线的公路梁式桥承载能力评定方法
CN115791033B (zh) * 2022-11-28 2023-08-22 安徽省交通规划设计研究总院股份有限公司 基于通行车辆荷载及动挠度监测的桥梁承载能力评定方法
CN116127583B (zh) * 2023-03-03 2023-08-29 大连理工大学 一种用于桥梁结构抗弯刚度重构的逆单位荷载法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105973619A (zh) * 2016-04-27 2016-09-28 厦门大学 结构健康监测系统下基于影响线的桥梁局部损伤识别方法
CN107588915A (zh) * 2017-10-18 2018-01-16 厦门大学 一种桥梁影响线识别方法及系统
AU2017266832A1 (en) * 2016-05-18 2018-12-20 Heijmans N.V. Method for determining the structural integrity of an infrastructural element

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105973619A (zh) * 2016-04-27 2016-09-28 厦门大学 结构健康监测系统下基于影响线的桥梁局部损伤识别方法
AU2017266832A1 (en) * 2016-05-18 2018-12-20 Heijmans N.V. Method for determining the structural integrity of an infrastructural element
CN107588915A (zh) * 2017-10-18 2018-01-16 厦门大学 一种桥梁影响线识别方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
桥梁位移影响线提取方法及其精度分析;徐剑;《现代交通技术》;20181026(第05期);全文 *

Also Published As

Publication number Publication date
CN109829252A (zh) 2019-05-31

Similar Documents

Publication Publication Date Title
CN109829252B (zh) 一种基于影响线识别的桥梁状况快速评级方法
Zhu et al. Recent developments in inverse problems of vehicle–bridge interaction dynamics
Yu et al. Moving force identification based on the frequency–time domain method
CN105005694B (zh) 一种基于动态称重系统的桥梁疲劳寿命频域分析方法
Fauriat et al. Estimation of road profile variability from measured vehicle responses
CN109668610A (zh) 基于神经网络回归的车辆动态测重方法及其采用的系统
CN111581867B (zh) 一种桥梁损伤快速检测方法
CN108804814B (zh) 一种提取路面谱参数的方法
Chan et al. A theoretical study of force identification using prestressed concrete bridges
CN112857718B (zh) 一种基于移动车辆测试的桥梁承载能力快速评估方法
CN102721397A (zh) 基于垂直动载荷的现代时间序列提取路面特征参数的方法
Liu et al. A semi-convex function for both constant and time-varying moving force identification
CN108228994B (zh) 在越野道路随机激励下车辆及工作装置的应力计算方法
CN112347541A (zh) 车辆正常通行下沥青路面承载能力快速测试方法
CN102998133A (zh) 一种基于准分布式加速度数据的能量损伤识别方法
CN113420362B (zh) 一种适用于时变中小跨径桥梁的损伤识别方法及装置
CN110987499A (zh) 一种桥梁动载试验方法
CN114323512B (zh) 一种重载车辆识别方法及系统
CN105066959A (zh) 路面纵断面的高程信息获取方法
CN110717250A (zh) 钢筋混凝土梁桥累积损伤简化标准等效车模型制定方法
CN111523180B (zh) 一种车载设备加速试验谱的构建方法
CN114541222B (zh) 基于多车众筹振动数据的路网级路面平整度检测方法
Prasad et al. Study on correlation of commercial vehicle axle response with road profile for ISO road class categorization and durability analysis
WO2019232737A1 (zh) 一种基于迭代法的准静态桥梁影响线识别方法
CN104266849B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220719

CF01 Termination of patent right due to non-payment of annual fee