CN112461358A - 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 - Google Patents
基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 Download PDFInfo
- Publication number
- CN112461358A CN112461358A CN202011320810.7A CN202011320810A CN112461358A CN 112461358 A CN112461358 A CN 112461358A CN 202011320810 A CN202011320810 A CN 202011320810A CN 112461358 A CN112461358 A CN 112461358A
- Authority
- CN
- China
- Prior art keywords
- bridge
- frequency
- time
- vehicle
- instantaneous frequency
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于车辆‑桥梁***瞬时频率的桥梁模态参数识别方法,是使两轴车辆以匀速在桥梁上行驶,并通过布置在桥梁上的少量传感器进行测试,获得两轴车辆‑桥梁耦合***的动力响应,再根据车辆‑桥梁耦合***的时变特性,通过同步提取变换进行时频分析,提取出与桥梁频率相关的瞬时频率,接着通过时频系数模极大值法确定时频脊线,选取时频脊线各阶梯段的中点作为有效数据点,然后通过瞬时频率与桥梁振型之间的物理关系,识别出桥梁振型,最后对车辆驶离桥梁后桥梁的动力响应进行傅里叶变换识别桥梁频率。本发明能将振型识别转化成瞬时频率识别,从而能够有效解决现有计算中处理非平稳信号精度和效率偏低等问题。
Description
技术领域
本发明涉及于桥梁安全检测领域,具体地说是一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法,识别结果可用于评估桥梁结构安全状态。
背景技术
随着经济的发展和时代变迁,越来越多的大型桥梁投入使用。与此同时人们也越来越来重视结构的安全性。模态分析在大型结构工程的安全诊断中起到了举足轻重的作用,模态参数可以用于监测和检测结构的健康状况,广泛应用于梁结构损伤识别和健康状态评估。桥梁的模态参数中至关重要时自然频率、模态振型,因而频率和振型识别成为了桥梁工程结构检测领域的重要任务之一
桥梁模态参数识别实际上是处理传感器采集到的各类动力响应的过程。传统直接测量方法是在桥梁结构上安装大量各种类型传感器,通过分析传感器采集的海量响应信号以获得其频率和振型,并利用模态参数变化规律来间接评估其健康状态。但是安装大量数目传感器耗资巨大且浪费人力;处理海量响应信号的工作量大且难度高。而且在实际应用中,在车辆桥梁质量相比并不小的情况下,车辆-桥梁***的响应信号是非平稳的,模态参数具有时变性。传统的方法忽略了车辆-桥梁***的时变特性,处理数据时精度偏低。模态参数识别精度直接影响其在桥梁安全诊断应用中的有效性。
发明内容
本发明是为避免上述传统桥梁模态识别方法所存在的缺陷,提供一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法,以期能通过少量(单个或者两个)传感器布置测试车辆-桥梁***瞬时频率的方法,将振型识别转化成瞬时频率识别,从而能够有效的解决传统的模态参数识别方法需要在桥梁上布置大量的传感器、耗资大、数据处理难度高、处理非平稳信号精度和效率偏低等问题。
本发明为达到上述发明目的,采用如下技术方案:
本发明一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法的特点是包含以下步骤:
步骤1:确定两轴车辆参数,包括:总重mv、两个轴距d1和d2和行驶速度v;
步骤2:在车辆以行驶速度v在桥梁上行驶时,利用加速度传感器在桥梁测点上实时获取t时刻的两轴车辆-桥梁耦合***的动力响应s1(t):
步骤3:通过时频处理方法对动力响应s1(t)进行时频分析,识别两轴车辆-桥梁耦合***的各阶瞬时频率;
步骤3.1:对所述动力响应s1(t)进行短时傅里叶变换,得到两轴车辆-桥梁耦合***的时频系数矩阵St(t,η),η为车辆-桥梁***的频率;
步骤3.2:利用式(1)初步估算两轴车辆-桥梁耦合***的二维瞬时频率η0(t,η):
式(1)中,i为复数单位;
步骤3.3:利用式(2)得到两轴车辆-桥梁耦合***的时频系数ω0(t,η):
步骤4:对时频系数ω0(t,η)进行模极大值脊线提取处理,获得每一列的模极大值并作为时频脊线,选取时频脊线图的每一段阶梯的中点,作为有效数据点,其中,第k个有效数据点记为(tk,ωn(tk)),其中,tk为第k个有效数据点所对应的车辆行驶时间,n代表桥梁的第n阶瞬时频率,ωn(tk)为车辆行驶时间为tk时桥梁第n阶瞬时频率;
步骤5:识别桥梁频率:
步骤5.1:当车辆离开桥梁后,利用加速度传感器采集桥梁的加速度响应s2(t);
步骤5.2:利用傅里叶变换对加速度响应s2(t)进行频谱分析,识别出桥梁的固有频率ωn;
步骤6:计算桥梁振型:
步骤6.1:利用式(3)计算时间为tk时车辆行驶到桥梁所对应的桥上位置与起始点的距离ak:
ak=v×tk-d1 (3)
步骤6.2:利用式(4)计算车辆行驶距起始点的距离为ak时的桥梁振型φn(ak):
以所述固有频率ωn和桥梁振型φn(ak)作为所识别的桥梁模态参数。
与现有技术相比,本发明的有益效果在于:
1、本发明通过使车辆匀速行驶过桥梁,仅需要布置少量(单个或者两个)数目的加速度传感器进行测试,通过采用先进的信号处理方法对采集到的加速度响应进行时频分析,通过得到的瞬时频率即可获得高精度的振型,从而有效解决了以往测试需要大量的传感器密集排列,安装设备昂贵耗时,需要多次重复操作,数据处理过程中工作量大的问题。
2、本发明很好的利用了车辆-桥梁***频率变化与该点处桥梁模态参数的物理关系,通化振型识别为瞬时频率识别,能显著提高识别的精度,整个过车辆过程采用先进的时频分析方法处理信号,操作简便且能识别出车辆-桥梁耦合***的时变特性,处理非平稳信号结果精度高。
3、本发明采用两轴车辆进行测试,在实际桥梁健康监测和检测中操作性更强。
附图说明
图1为本发明方法的识别过程示意图;
图2为本发明数值模拟等截面简支梁桥图;
图3为本发明车辆激励下等截面简支梁桥的典型动力响应图;
图4为本发明车辆激励下等截面简支梁桥的典型动力响应的时频分析图;
图5为本发明车辆激励下等截面简支梁桥的典型动力响应的时频脊线图
图6为本发明车辆驶离后等截面简支梁桥的典型动力响应的频谱图
图7为本发明等截面简支梁桥一阶振型识别结果图;
图8为本发明等截面简支梁桥二阶振型识别结果图;
图9为本发明数值模拟两跨等截面连续梁桥图;
图10为本发明两跨等截面连续梁桥一阶振型识别结果图;
图11为本发明两跨等截面连续梁桥二阶振型识别结果图。
具体实施方式
本实施例1中,如图2所示的等截面简支梁桥,桥梁跨长为30m,弹性模量27.5Gpa,每延米质量2000kg/m。采用有限元法模拟时,桥梁等距划分为30个平面欧拉梁单元。通过Newmark-β法计算桥梁的动力响应。采样频率为1000Hz。一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法,如图1所示,包含以下步骤:
步骤1:确定两轴车辆参数,包括:总重mv=1800kg、两个轴距d1=0.2m和d2=0.2m、行驶速度v=0.5m/s;
步骤2:桥梁测点位置选择距车辆起始位置10m处,在车辆以行驶速度v在桥梁上行驶时,利用加速度传感器在桥梁测点上实时获取t时刻的两轴车辆-桥梁耦合***的动力响应s1(t),两轴车辆-桥梁耦合***的动力响应s1(t)如图3所示;
步骤3:通过时频处理方法对动力响应s1(t)进行时频分析,识别两轴车辆-桥梁耦合***的各阶瞬时频率;
步骤3.1:对动力响应s1(t)进行短时傅里叶变换,得到两轴车辆-桥梁耦合***的时频系数矩阵St(t,η),η为车辆-桥梁***的频率;
步骤3.2:利用式(1)初步估算两轴车辆-桥梁耦合***的二维瞬时频率η0(t,η):
式(1)中,i为复数单位;
步骤3.3:利用式(2)得到两轴车辆-桥梁耦合***的时频系数ω0(t,η),时频图如图4所示:
步骤4:对时频系数ω0(t,η)进行模极大值脊线提取处理,获得每一列的模极大值并作为时频脊线,选取时频脊线的每一段阶梯的中点,作为有效数据点,其中,第k个有效数据点记为(tk,ωn(tk)),其中,tk为第k个有效数据点所对应的车辆行驶时间,n代表桥梁的第n阶瞬时频率,ωn(tk)为车辆行驶时间为tk时桥梁第n阶瞬时频率,时频脊线如图5所示;
步骤5:识别桥梁频率:
步骤5.1:当车辆离开桥梁后,利用加速度传感器采集桥梁的加速度响应s2(t);
步骤5.2:利用傅里叶变换对加速度响应s2(t)进行频谱分析,频谱如图6所示,识别出桥梁的固有频率ωn;
步骤6:计算桥梁振型:
步骤6.1:利用式(3)时间为tk时车辆行驶到桥梁所对应的桥上位置与起始点的距离ak,即可将ωn(tk)转换为ωn(ak):
ak=v×tk-d1 (3)
步骤6.2:利用车辆-桥梁耦合***瞬时频率与对应点处桥梁振型幅值的物理关系式(4)计算车辆行驶距起始点的距离为ak时的桥梁振型φn(ak),第一阶和第二阶的识别结果与参照值对比分别如表1、图7和图8所示。
以固有频率ωn和桥梁振型φn(ak)作为所识别的桥梁模态参数。
表1
本实施例2中,如图9所示的两跨等截面连续梁桥,每跨长度为15m,弹性模量27.5Gpa,每延米质量2000kg/m。采用有限元法模拟时,桥梁等距划分为30个平面欧拉梁单元。通过Newmark-β法计算桥梁的加速度响应。采样频率为1000Hz,车辆的参数与实施例1相同。求解方法和识别过程与实施例1过程一致,第一阶和第二阶频率和振型的识别结果与参照值对比分别如表2、图10和图11所示。
表2
从表1和表2中可以看出本发明提出的一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法能较精准地识别出不同类桥梁的模态参数,误差小,适用和实用性强。
Claims (1)
1.一种基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法,其特征是包含以下步骤:
步骤1:确定两轴车辆参数,包括:总重mv、两个轴距d1和d2和行驶速度v;
步骤2:在车辆以行驶速度v在桥梁上行驶时,利用加速度传感器在桥梁测点上实时获取t时刻的两轴车辆-桥梁耦合***的动力响应s1(t):
步骤3:通过时频处理方法对动力响应s1(t)进行时频分析,识别两轴车辆-桥梁耦合***的各阶瞬时频率;
步骤3.1:对所述动力响应s1(t)进行短时傅里叶变换,得到两轴车辆-桥梁耦合***的时频系数矩阵St(t,η),η为车辆-桥梁***的频率;
步骤3.2:利用式(1)初步估算两轴车辆-桥梁耦合***的二维瞬时频率η0(t,η):
式(1)中,i为复数单位;
步骤3.3:利用式(2)得到两轴车辆-桥梁耦合***的时频系数ω0(t,η):
步骤4:对时频系数ω0(t,η)进行模极大值脊线提取处理,获得每一列的模极大值并作为时频脊线,选取时频脊线图的每一段阶梯的中点,作为有效数据点,其中,第k个有效数据点记为(tk,ωn(tk)),其中,tk为第k个有效数据点所对应的车辆行驶时间,n代表桥梁的第n阶瞬时频率,ωn(tk)为车辆行驶时间为tk时桥梁第n阶瞬时频率;
步骤5:识别桥梁频率:
步骤5.1:当车辆离开桥梁后,利用加速度传感器采集桥梁的加速度响应s2(t);
步骤5.2:利用傅里叶变换对加速度响应s2(t)进行频谱分析,识别出桥梁的固有频率ωn;
步骤6:计算桥梁振型:
步骤6.1:利用式(3)计算时间为tk时车辆行驶到桥梁所对应的桥上位置与起始点的距离ak:
ak=v×tk-d1 (3)
步骤6.2:利用式(4)计算车辆行驶距起始点的距离为ak时的桥梁振型φn(ak):
以所述固有频率ωn和桥梁振型φn(ak)作为所识别的桥梁模态参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011320810.7A CN112461358B (zh) | 2020-11-23 | 2020-11-23 | 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011320810.7A CN112461358B (zh) | 2020-11-23 | 2020-11-23 | 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112461358A true CN112461358A (zh) | 2021-03-09 |
CN112461358B CN112461358B (zh) | 2022-04-26 |
Family
ID=74799555
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011320810.7A Active CN112461358B (zh) | 2020-11-23 | 2020-11-23 | 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112461358B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113639856A (zh) * | 2021-08-17 | 2021-11-12 | 合肥工业大学 | 一种考虑环境温度影响的梁桥模态频率识别方法 |
CN114323512A (zh) * | 2021-12-17 | 2022-04-12 | 中铁桥隧技术有限公司 | 一种重载车辆识别方法及*** |
CN114964456A (zh) * | 2022-04-20 | 2022-08-30 | 中国矿业大学 | 一种利用车辆速度监测的梁式桥钢结构振动安全监测*** |
CN115510724A (zh) * | 2022-10-27 | 2022-12-23 | 安徽省交通控股集团有限公司 | 一种基于移动车辆测试的桥梁损伤识别方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4901575A (en) * | 1988-11-30 | 1990-02-20 | Gp Taurio, Inc. | Methods and apparatus for monitoring structural members subject to transient loads |
JP2009204517A (ja) * | 2008-02-28 | 2009-09-10 | Tokyo Institute Of Technology | 振動解析方法及び振動解析装置 |
CN102565194A (zh) * | 2012-02-09 | 2012-07-11 | 东南大学 | 运营状态下大跨桥梁钢箱梁损伤预警方法 |
KR101284125B1 (ko) * | 2012-03-08 | 2013-07-10 | 인하대학교 산학협력단 | 교량의 동적응답 측정을 위한 이동식 측정시스템 |
CN103411743A (zh) * | 2013-07-04 | 2013-11-27 | 东南大学 | 利用可移动设备对桥梁进行移动式冲击振动的测试方法 |
CN104457681A (zh) * | 2014-09-18 | 2015-03-25 | 合肥工业大学 | 一种基于应变模态的梁结构动挠度监测方法 |
CN106197910A (zh) * | 2016-07-01 | 2016-12-07 | 东南大学 | 一种基于车桥耦合振动分析的桥梁检测方法与检测*** |
CN106802222A (zh) * | 2017-01-13 | 2017-06-06 | 重庆大学 | 一种基于车桥耦合***的桥梁损伤诊断方法 |
CN107009866A (zh) * | 2017-04-06 | 2017-08-04 | 北京航空航天大学 | 一种面向车辆发动机垂直方向振动的减振方法 |
JP2018031676A (ja) * | 2016-08-24 | 2018-03-01 | 公益財団法人鉄道総合技術研究所 | 鉄道橋の構造性能調査方法 |
US20180149554A1 (en) * | 2016-11-29 | 2018-05-31 | Korea Institute Of Civil Engineering And Building Technology | Safety inspection apparatus for bridge using expansion joint with load cell and method for the same |
CN109357822A (zh) * | 2018-08-13 | 2019-02-19 | 东南大学 | 一种基于车桥耦合***时变动力特征改变的桥梁快速测试与评估方法 |
CN109813511A (zh) * | 2019-01-31 | 2019-05-28 | 东南大学 | 基于移动车辆的桥梁快速测试与参数识别方法 |
CN109839441A (zh) * | 2019-03-20 | 2019-06-04 | 合肥工业大学 | 一种桥梁模态参数识别方法 |
-
2020
- 2020-11-23 CN CN202011320810.7A patent/CN112461358B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4901575A (en) * | 1988-11-30 | 1990-02-20 | Gp Taurio, Inc. | Methods and apparatus for monitoring structural members subject to transient loads |
JP2009204517A (ja) * | 2008-02-28 | 2009-09-10 | Tokyo Institute Of Technology | 振動解析方法及び振動解析装置 |
CN102565194A (zh) * | 2012-02-09 | 2012-07-11 | 东南大学 | 运营状态下大跨桥梁钢箱梁损伤预警方法 |
KR101284125B1 (ko) * | 2012-03-08 | 2013-07-10 | 인하대학교 산학협력단 | 교량의 동적응답 측정을 위한 이동식 측정시스템 |
CN103411743A (zh) * | 2013-07-04 | 2013-11-27 | 东南大学 | 利用可移动设备对桥梁进行移动式冲击振动的测试方法 |
CN104457681A (zh) * | 2014-09-18 | 2015-03-25 | 合肥工业大学 | 一种基于应变模态的梁结构动挠度监测方法 |
CN106197910A (zh) * | 2016-07-01 | 2016-12-07 | 东南大学 | 一种基于车桥耦合振动分析的桥梁检测方法与检测*** |
JP2018031676A (ja) * | 2016-08-24 | 2018-03-01 | 公益財団法人鉄道総合技術研究所 | 鉄道橋の構造性能調査方法 |
US20180149554A1 (en) * | 2016-11-29 | 2018-05-31 | Korea Institute Of Civil Engineering And Building Technology | Safety inspection apparatus for bridge using expansion joint with load cell and method for the same |
CN106802222A (zh) * | 2017-01-13 | 2017-06-06 | 重庆大学 | 一种基于车桥耦合***的桥梁损伤诊断方法 |
CN107009866A (zh) * | 2017-04-06 | 2017-08-04 | 北京航空航天大学 | 一种面向车辆发动机垂直方向振动的减振方法 |
CN109357822A (zh) * | 2018-08-13 | 2019-02-19 | 东南大学 | 一种基于车桥耦合***时变动力特征改变的桥梁快速测试与评估方法 |
CN109813511A (zh) * | 2019-01-31 | 2019-05-28 | 东南大学 | 基于移动车辆的桥梁快速测试与参数识别方法 |
CN109839441A (zh) * | 2019-03-20 | 2019-06-04 | 合肥工业大学 | 一种桥梁模态参数识别方法 |
Non-Patent Citations (3)
Title |
---|
SUN, JH: "Influence of Variation of Damping Ratio on Vehicle-bridge Coupling Vibration Responses of Bridge with Corrugated Steel Webs", 《CIVIL ENGINEERING IN CHINA - CURRENT PRACTICE AND RESEARCH REPORT》 * |
申启坤: "车—桥耦合振动作用下高速公路高陡边坡动力响应研究", 《中国优秀硕士论文电子期刊网 工程科技Ⅱ辑》 * |
贺文宇: "基于间接法识别的桥梁振型的损伤定位方法", 《振动与冲击》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113639856A (zh) * | 2021-08-17 | 2021-11-12 | 合肥工业大学 | 一种考虑环境温度影响的梁桥模态频率识别方法 |
CN113639856B (zh) * | 2021-08-17 | 2022-10-14 | 合肥工业大学 | 一种考虑环境温度影响的梁桥模态频率识别方法 |
CN114323512A (zh) * | 2021-12-17 | 2022-04-12 | 中铁桥隧技术有限公司 | 一种重载车辆识别方法及*** |
CN114323512B (zh) * | 2021-12-17 | 2023-06-16 | 中铁桥隧技术有限公司 | 一种重载车辆识别方法及*** |
CN114964456A (zh) * | 2022-04-20 | 2022-08-30 | 中国矿业大学 | 一种利用车辆速度监测的梁式桥钢结构振动安全监测*** |
CN115510724A (zh) * | 2022-10-27 | 2022-12-23 | 安徽省交通控股集团有限公司 | 一种基于移动车辆测试的桥梁损伤识别方法 |
CN115510724B (zh) * | 2022-10-27 | 2023-06-27 | 安徽省交通控股集团有限公司 | 一种基于移动车辆测试的桥梁损伤识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112461358B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112461358B (zh) | 基于车辆-桥梁***瞬时频率的桥梁模态参数识别方法 | |
CN109357822B (zh) | 一种基于车桥耦合***时变动力特征改变的桥梁快速测试与评估方法 | |
CN102353509B (zh) | 一种基于分块冲击振动测试的中小型桥梁快速检测方法 | |
CN104165742B (zh) | 一种基于互谱函数的运行模态分析实验方法及装置 | |
CN109839440B (zh) | 一种基于静置车辆测试的桥梁损伤定位方法 | |
CN109839441B (zh) | 一种桥梁模态参数识别方法 | |
CN101504324A (zh) | 拉索运营状态智能判别方法及*** | |
CN105067239A (zh) | 基于扫频激励振动的梁裂纹故障检测装置及方法 | |
CN102721462B (zh) | 旋转机械启停车过程波德图/奈奎斯特图的快速计算方法 | |
CN102520070A (zh) | 基于非线性输出频率响应函数的结构损伤检测方法 | |
CN113639856B (zh) | 一种考虑环境温度影响的梁桥模态频率识别方法 | |
CN101871974B (zh) | 一种阻抗谱的测量方法 | |
CN113343541B (zh) | 长大跨桥梁的涡激振动预警方法、装置及终端 | |
CN111487318B (zh) | 一种时变结构瞬时频率提取方法 | |
CN104408322A (zh) | 综合多源故障概率似然信度的旋转机械设备故障诊断方法 | |
CN113295310A (zh) | 一种基于应变刚度代表值的桥梁损伤判定方法 | |
CN106918389A (zh) | 一种基于多普勒光学位移法的振型分析方法及其应用 | |
CN109543357B (zh) | 一种多元回归模型优化的故障程度量化评估方法 | |
CN109781863B (zh) | 基于快速振动测试的结构腐蚀二阶段探测方法及其*** | |
CN107330264B (zh) | 一种桥梁监测数据可靠性的验证方法 | |
CN105651537B (zh) | 一种高损伤敏感性的桁架结构损伤实时监测*** | |
CN116007930B (zh) | 一种汽车传动性能的测试方法及*** | |
CN112326787A (zh) | 一种基于专属测试车多点快速静态采集的梁式桥识别方法 | |
CN103698071B (zh) | 基于监测加速度的拉索时变索力历程识别的数据驱动方法 | |
CN113705721B (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 |