CN107145620A - 一种基于随机减量技术的结构动力特性识别方法 - Google Patents
一种基于随机减量技术的结构动力特性识别方法 Download PDFInfo
- Publication number
- CN107145620A CN107145620A CN201710148179.9A CN201710148179A CN107145620A CN 107145620 A CN107145620 A CN 107145620A CN 201710148179 A CN201710148179 A CN 201710148179A CN 107145620 A CN107145620 A CN 107145620A
- Authority
- CN
- China
- Prior art keywords
- random decrement
- random
- measuring point
- model
- fem model
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/025—Measuring arrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种基于随机减量技术的结构动力特性识别方法,具体实施步骤如下:A.实物模型实验:进行系统搭载、调试,布置测点,输出在一定激励下的相应方向周加速度响应;B.根据实物模型的模型参数进行有限元模型建模,在不同的时程工况下,对有限元模型进行加速度激励,得出不同阶的振型及频率;C.测点布置及相应输出:根据有限元模型结构布置测点,在给定的采样频率下输出所需坐标轴方向的加速度响应;D.对响应数据进行处理:在给定的触发条件下,提取时间片段长度,采用随机减量技术和Ibrahim时域法对各测点进行处理,并对所得到的数据进行对比分析,来检验随机减量技术在动力特性方面的适用范围和精度。
Description
技术领域
本发明涉及到使用随机减量技术,并与模态识别方法如Ibrahim时域法、协方差驱动随机子空间法配合使用进行动力特性识别的方法。
背景技术
在突发事故、地震、风载和运营载荷的作用下,建筑和桥梁会产生损伤累积并由此引发破坏性事故而带来巨大的人员伤亡与财产损失。因此,对于运行中的建筑和桥梁实行长期健康状态的评估非常必要。
传统的建筑和桥梁的损伤测试和评估主要是无损探测和人工视察为主,这种损伤识别方法只能在人力所及范围内的结构和构件上出现肉眼可见的缺陷时方可见效。显然,对于近来大量建成的大跨度桥梁这样的大型工程结构,这样的检测和评估方法已经远远地落后于形势的要求。因此,大型工程结构的模态识别已经成为现代动态测试和复杂结构在线监测的核心技术之一。
早在60年代,环境激励下结构工作模态的研究就已经开始,经过这几十年的研究,特别是近几年来,人们已经提出了多种环境激励下的模态参数识别方法。按识别信号域不同可分为:时域识别方法、频域识别方法和时频域识别方法;按激励信号分为:平稳随机激励和非平稳随机激励(有的方法假设环境激励为白噪声激励);按信号的测取方法分为:单输入多输出和多输入多输出;按识别方法特性分为:时间序列法、随机减量法、NExT法、随机子空间法、峰值拾取法、频域分解法等,以下对随机减量技术进行阐述:
随机减量技术:这种方法是基于在NASA关于空间结构对环境激励动态响应的工作,该工作的主要目的是在环境激励下通过测得的响应对空间结构进行动力特性识别和在役损伤监测。
对一单自由度结构在激励下得到的平稳随机响应x(t),取幅值x(ti)=a与响应x(t)交于N个点。以x(ti)为初始采样值,并对其后时间片段长度为τ的响应进行采样,将N个采样响应片段平均,有:
由振动理论可知,对一线性结构,其任意激励作用下的响应由三部分组成:
1)结构对初始位移的响应;
2)结构对初始速率的响应;
3)结构对激励的响应(受迫振动暂态响应及稳态响应)。
当激励为随机过程时,结构对激励的响应也是随机的。通过在同样的初始条件下将大量的响应片段平均,响应中的激励部分会在均值中趋向于零。而由于初始速度在各片段中正负交替,因此在大量平均后也趋于零。因此最终得到的便是结构对初始位移的响应,也就是说随机减量函数即结构由初始位移激励而引起的自由响应。
随机减量技术的初步设想只包含自谱随机减量函数,其中触发条件和时间片段在同一个响应中定义。后来,互谱随机减量函数的概念被提出,其触发条件在一个响应中定义而时间片段从另一个同时测得的响应中取得。我们因此能够求得一个完整的随机减量函数矩阵而无需对大量监测点进行检测采样。
考虑两个同时测得的响应x(t)和y(t),随机减量函数的自谱DXX(τ)和互谱DXY(τ)可以通过以下表达式定义:
在这个例子中,N表示平均时间片段的数量;Tx(ti)表示定义于时程x(t)的触发条件。
为了对随机减量函数求值,我们需要考虑不同的触发条件Tx(ti)。随机减量技术中最常见的触发条件有:
1)水平穿越;
Tx(ti)={x(ti)=a} (3)
2)正值点;
Tx(ti)={a≤x(ti)<b} (4)
3)斜率为正的零值穿越;
4)局部极值;
随机减量技术应用的一个重要问题在于公式(3)和(4)中触发条件a和b的定义。一般而言,使用大量的触发点是有优势的,但数值较小的触发点相比于数值较大的更容易被噪声污染。因此,我们需要在较多的触发点和较高的触发条件中达成一个平衡。一个常用且合适的选择是在正值点触发条件中采用a=σx,b=∞,其中a=σx是分析响应的标准差。同样对最佳触发条件进行定义,即使计算所得的随机减量函数方差最小。在考虑互谱触发条件时,该值为
随机减量技术的另一个重要方面是从响应时程中提取时间片段的长度。在与时域模态识别方法配合应用时,随机减量函数需要包含足够的采样点以满足系统测得矩阵的超定公式(考虑最高阶)。在与频域模态识别方法配合应用时,随机减量函数需要有足够的长度以包含一个完整的自由衰减振动(一般而言,频域模态识别方法需要时间片段的长度比时域方法更长)。
一般认为,自谱随机减量函数比互谱随机减量函数受更少的噪声影响。为解决这一问题,矢量触发随机减量技术(Vector Triggering Random Decrement Technique)被提出,其中触发条件在两个或更多同时测得的响应中定义。
简而言之,随机减量技术是一种将随机响应转化为结构自由振动的方法。即在触发条件下,对测得响应以一定长度的时间片段采样平均,求得的随机减量函数即结构在初始位移条件下的自由振动响应。而在之后随着研究的不断深入,可以证明在基于响应是均值为零的平稳高斯随机过程的假设下,随机减量函数和相关函数成正比。
随机减量技术的大部分应用与时域模态识别方法相结合,如Ibrahim时域法(Ibrahim Time Domain(ITD)Method)、随机子空间法(Stochastic Subspace Method(SSI))和特征系统实现算法(Eigensystem Realization Algorithm(ERA))。考虑到随机减量技术是一个时域过程,所以很显然与时域模态识别方法配合使用。
随机减量技术同样可以应用于谱密度计算和频率响应函数测量,因此它可以和频域模态识别方法配合使用。系统响应的谱密度通过对随机减量函数进行傅立叶快速变换计算得出,而不是更常见的通过对响应时程窗口采样进行傅立叶快速变换算法,该方法被称为集合平均法或Welch法。这个利用傅立叶快速变换算法从随机减量函数中计算谱密度并在频域输出模态识别方法中应用的计算过程在一个例子中被采用,与传统方法相比表现出巨大的优势。
随机减量技术相比于其他识别方法更具优势,但是仍然存在一些缺陷。为了保证模态参数的准确识别,需要基于原有的随机减量技术进行改进。
发明内容
本发明要克服原有的模态识别方法——随机减量技术的应用范围的局限性,提出一种基于随机减量技术的结构动力特性识别方法。测试方法采用多种有限元一种基于Midas Civil 2011进行有限元建模,在一定的设定工况下采用振型叠加法对瞬态时程进行分析,并采用振型阻尼计算方法计算阻尼的方法。测试装置为计算机、软件与设计结构模型。
本发明要解决以下几个方面的问题:
一是解决检验随机减量技术对随机过程以高斯分布的前提的敏感性问题。这里通过研究高斯白噪声对于理想线性结构产生的响应在随机减量技术处理后得到的函数以及从中提取的模态参数,并以此为基准,添加非高斯分布激励,改变高斯白噪声在激励中的比例,比较不同激励下的随机减量函数及模态参数的准确性及稳定性。
二是解决验证随机减量技术对非线性结构的应用问题。这里仍然通过添加高斯分布荷载作用将非线性结构转化为非高斯分布,因此这可以视作前一种情况下仅环境噪声是非高斯分布的延伸。
三是解决减少噪声的影响以及剔除虚假模态的问题。这里将采用将噪声视作随机过程的模态提取方法使提取的模态参数更精确、减少模态的数量,并能在一定程度上帮助解决利用自然激励测定超低频结构时分离频率相近模态时的困难
四是解决通过随机减量技术对结构进行实时在线健康监测的问题。这里通过计算,大量长期观察数据可以被转化为一段较短的随机减量函数。从中提取得到的模态参数(或随机减量函数本身)可以添加至结构动力特性监测序列,通过比对监测结构是否出现损伤。
本发明所述的一种基于随机减量技术的结构动力特性识别方法,具体实施步骤如下(包括有限元模型试验与设计模型试验两个部分):
A.实物模型实验;
A1.准备结构模型,进行系统搭载、调试;
A2.根据模型的结构特性在其上布置测点;
A3.利用振动台对结构模型施加荷载工况所要求的信号,在计算机中输出所需坐标轴方向加速度响应;
B.有限元模型建模;
B1.基于Midas Civil 2011根据实物模型进行有限元建模,确定结构尺寸,材料数据,边界条件,时程荷载工况,线性分析类型,计算方法以及相关参数;
B2.在不同时程工况下,对有限元模型施加不同的加速度激励;
B3.模态参数理论值采用Midas Civil内置子空间迭代法进行计算,得到不同阶的振型及频率。
C.测点布置及响应输出;
C1.根据有限元模型结构特点布置测点;
C2.在给定的采样频率下输出所需坐标轴方向的加速度响应。
D.对响应数据进行处理;
D1.给定触发条件与提取时间片段长度;
D2.采用随机减量技术对实物模型与有限元模型的各工况各测点响应数据进行处理,得出高斯白噪声激励下有限元模型各测点的随机减量函数;
D3.采用Ibrahim时域法对各测点的随机减量函数进行处理,根据模型计算模态参数;
D4.对比分析有限元模型在不同激励下的各阶模态参数与理论值以及实际各参数相同的实物模型实验值之间的差异。
本发明中通过多组不同的对比——不同的激励间的对比、不同的有限元模型结构之间的对比、有限元模型与实物模型之间的对比,对随机减量技术的适用性进行了充分的检验。
与现有的技术相比,本技术有几下几个优点:
1.解决了检验随机减量技术对随机过程以高斯分布的前提的敏感性问题。
2.解决了验证随机减量技术对非线性结构的应用问题。
3.解决了减少噪声的影响以及剔除虚假模态的问题。
4.解决了通过随机减量技术对结构进行实时在线健康监测的问题。
5.随机减量-Ibrahim时域法在高斯白噪声激励下,对于简单结构或者复杂结构的前两阶振型的频率有较高的识别精度。
6.实验时多方面的对比更能检验随机减量函数在模态识别时的精度与实用性。
附图说明
图1本发明的测量案例,
图2本发明的测量流程图。
图例说明:图1中的代号分别表示:
1—二层框架设计模型
2—测点
3—振动台
4—计算机
5—二层框架有限元模型
6—模态处理数据
具体实施方式
以下结合图1中所示的案例和图2中所示的工作流程进一步阐述本发明。
参见图1和图2,本发明中所举案例为利用一种基于随机减量技术的结构动力特性识别方法对二层框架模型进行模态识别,以进而完成对其的动力特性识别,具体的实施步骤如下:
A.实物模型实验;
A1.准备二层框架模型(1),系统搭载、调试;
A2.在二层框架上根据其结构特性选择布置6个测点(2);
A3.用振动台(3)对模型施加内置地震信号,在计算机(4)中输出Y轴方向加速度响应;
B.有限元模型建模;
B1.基于Midas Civil 2011根据实物模型进行有限元建立二层框架模型(5),确定结构尺寸,材料数据,边界条件,时程荷载工况,线性分析类型,计算方法以及相关参数;
B2.在不同时程工况下,对二层框架模型施加不同的加速度激励;
B3.模态参数理论值采用Midas Civil内置子空间迭代法进行计算,得到不同阶的振型及频率。
C.测点布置及响应输出;
C1.根据二层框架模型布置测点(2);
C2.在给定的采样频率下输出所需坐标轴方向的加速度响应。
D.对响应数据进行处理(6);
D1.给定触发条件与提取时间片段长度;
D2.采用随机减量技术对各工况各测点响应数据对二层框架实物模型和其有限元模型进行处理,得出高斯白噪声激励下两种模型各测点的随机减量函数;
D3.采用Ibrahim时域法对各测点的随机减量函数进行处理,计算二层框架结构前两阶模态参数;
D4.对比分析各激励下二层框架有限元模型的各阶模态参数与理论值以及其与实物模型所得实验值的差异,来检验随机减量技术的适用范围和精度。
本说明书实施案例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施案例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (1)
1.一种基于随机减量技术的桥梁动力特性识别方法,具体实施步骤如下:
A.实物模型实验;
A1.准备结构模型,进行系统搭载、调试;
A2.根据模型的结构特性在其上布置测点;
A3.利用振动台对结构模型施加荷载工况所要求的信号,在计算机中输出所需坐标轴方向加速度响应;
B.有限元模型建模;
B1.基于Midas Civil 2011根据实物模型进行有限元建模,确定结构尺寸,材料数据,边界条件,时程荷载工况,线性分析类型,计算方法以及相关参数;
B2.在不同时程工况下,对有限元模型施加不同的加速度激励;
B3.模态参数理论值采用Midas Civil内置子空间迭代法进行计算,得到不同阶的振型及频率;
C.测点布置及响应输出;
C1.根据有限元模型结构特点布置测点;
C2.在给定的采样频率下输出所需坐标轴方向的加速度响应;
D.对响应数据进行处理;
D1.给定触发条件与提取时间片段长度;
D2.采用随机减量技术对实物模型与有限元模型的各工况各测点响应数据进行处理,得出高斯白噪声激励下有限元模型各测点的随机减量函数;
D3.采用Ibrahim时域法对各测点的随机减量函数进行处理,根据模型计算模态参数;
D4.对比分析有限元模型在不同激励下的各阶模态参数与理论值以及实际各参数相同的实物模型实验值之间的差异。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710148179.9A CN107145620B (zh) | 2017-03-14 | 2017-03-14 | 一种基于随机减量技术的结构动力特性识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710148179.9A CN107145620B (zh) | 2017-03-14 | 2017-03-14 | 一种基于随机减量技术的结构动力特性识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145620A true CN107145620A (zh) | 2017-09-08 |
CN107145620B CN107145620B (zh) | 2019-08-06 |
Family
ID=59784150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710148179.9A Active CN107145620B (zh) | 2017-03-14 | 2017-03-14 | 一种基于随机减量技术的结构动力特性识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145620B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109520611A (zh) * | 2018-11-08 | 2019-03-26 | 温州大学 | 一种地震模拟振动台工况的监测方法 |
CN113375786A (zh) * | 2021-05-17 | 2021-09-10 | 国网浙江海盐县供电有限公司 | 一种基于振动信号的变压器监测方法及铁芯临时修复装置 |
CN114721294A (zh) * | 2022-04-19 | 2022-07-08 | 北京博科测试系统股份有限公司 | 一种利用传递函数实现混合试验的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004069598A (ja) * | 2002-08-08 | 2004-03-04 | Yamato Sekkei Kk | 構造物の損傷推定システムおよびプログラム |
KR20040067484A (ko) * | 2003-01-23 | 2004-07-30 | 한국항공우주산업 주식회사 | 항공기 시험 분석 시스템 |
EP1809999A1 (en) * | 2004-09-16 | 2007-07-25 | Sekos, Inc. | Method for detecting the onset of structural collapse of burning structures |
CN101226078A (zh) * | 2008-01-30 | 2008-07-23 | 广厦建设集团有限责任公司 | 一种基于分布式光纤传感器的长距离线性结构异常振动的检测方法 |
CN103198184A (zh) * | 2013-03-27 | 2013-07-10 | 深圳大学 | 一种电力系统低频振荡特征类噪声辨识方法 |
CN103983412A (zh) * | 2014-05-30 | 2014-08-13 | 北京航空航天大学 | 用于振动有限元模型修正航空电子设备工作模态测量方法 |
CN104133959A (zh) * | 2014-07-28 | 2014-11-05 | 东北大学 | 桥梁有限元模型的修正方法 |
CN104778514A (zh) * | 2015-04-17 | 2015-07-15 | 重庆交通大学 | 基于复杂系统理论的桥梁或构件安全状态预测方法 |
CN105865735A (zh) * | 2016-04-29 | 2016-08-17 | 浙江大学 | 一种基于视频监控的桥梁振动测试与动力特性识别方法 |
-
2017
- 2017-03-14 CN CN201710148179.9A patent/CN107145620B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004069598A (ja) * | 2002-08-08 | 2004-03-04 | Yamato Sekkei Kk | 構造物の損傷推定システムおよびプログラム |
KR20040067484A (ko) * | 2003-01-23 | 2004-07-30 | 한국항공우주산업 주식회사 | 항공기 시험 분석 시스템 |
EP1809999A1 (en) * | 2004-09-16 | 2007-07-25 | Sekos, Inc. | Method for detecting the onset of structural collapse of burning structures |
CN101226078A (zh) * | 2008-01-30 | 2008-07-23 | 广厦建设集团有限责任公司 | 一种基于分布式光纤传感器的长距离线性结构异常振动的检测方法 |
CN103198184A (zh) * | 2013-03-27 | 2013-07-10 | 深圳大学 | 一种电力系统低频振荡特征类噪声辨识方法 |
CN103983412A (zh) * | 2014-05-30 | 2014-08-13 | 北京航空航天大学 | 用于振动有限元模型修正航空电子设备工作模态测量方法 |
CN104133959A (zh) * | 2014-07-28 | 2014-11-05 | 东北大学 | 桥梁有限元模型的修正方法 |
CN104778514A (zh) * | 2015-04-17 | 2015-07-15 | 重庆交通大学 | 基于复杂系统理论的桥梁或构件安全状态预测方法 |
CN105865735A (zh) * | 2016-04-29 | 2016-08-17 | 浙江大学 | 一种基于视频监控的桥梁振动测试与动力特性识别方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109520611A (zh) * | 2018-11-08 | 2019-03-26 | 温州大学 | 一种地震模拟振动台工况的监测方法 |
CN113375786A (zh) * | 2021-05-17 | 2021-09-10 | 国网浙江海盐县供电有限公司 | 一种基于振动信号的变压器监测方法及铁芯临时修复装置 |
CN113375786B (zh) * | 2021-05-17 | 2022-01-25 | 国网浙江海盐县供电有限公司 | 一种基于振动信号的变压器监测方法及铁芯临时修复装置 |
CN114721294A (zh) * | 2022-04-19 | 2022-07-08 | 北京博科测试系统股份有限公司 | 一种利用传递函数实现混合试验的方法 |
CN114721294B (zh) * | 2022-04-19 | 2022-10-28 | 北京博科测试系统股份有限公司 | 一种利用传递函数实现混合试验的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107145620B (zh) | 2019-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107389285B (zh) | 一种基于温度改变的桥梁快速测试与评估方法 | |
Li et al. | Cluster analysis of winds and wind-induced vibrations on a long-span bridge based on long-term field monitoring data | |
Zhu et al. | A sensitivity-based structural damage identification method with unknown input excitation using transmissibility concept | |
CN106092479B (zh) | 板梁结构载荷识别和损伤识别的多功能试验台 | |
Mao et al. | The construction and comparison of damage detection index based on the nonlinear output frequency response function and experimental analysis | |
CN105260568B (zh) | 基于离散型卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN107145620B (zh) | 一种基于随机减量技术的结构动力特性识别方法 | |
CN109035730A (zh) | 一种考虑服役环境影响的混凝土坝损伤动力预警方法 | |
CN101706355A (zh) | 基于NExT/ARMA的结构响应分析方法 | |
CN101158623A (zh) | 获取系统特征函数和信号特征值的方法 | |
CN101901209A (zh) | 基于改进emd和arma模型的结构响应分析方法 | |
CN103513018A (zh) | 混凝土抗裂性能系统化检测方法 | |
CN104236893A (zh) | 一种液压阻尼器性能参数测试系统及测试方法 | |
Zhang et al. | Field measurement study on time-varying characteristics of modal parameters of super high-rise buildings during super typhoon | |
CN105466661A (zh) | 基于改进卡尔曼滤波的超高层建筑风荷载反分析方法 | |
Li et al. | Covariance of dynamic strain responses for structural damage detection | |
CN101539493A (zh) | 用于结构损伤诊断的对称信号方法 | |
Su et al. | Experimental identification of exponential damping for reinforced concrete cantilever beams | |
Li et al. | Structural damage identification based on short-time temporal coherence using free-vibration response signals | |
CN115659737A (zh) | 一种基于曲率模态和离散小波变换的梁构件损伤识别方法 | |
Wang et al. | Discrete analytical mode decomposition with automatic bisecting frequency selection for structural dynamic response analysis and modal identification | |
CN105651537B (zh) | 一种高损伤敏感性的桁架结构损伤实时监测系统 | |
Wang et al. | Influence of wind turbulence on aerodynamic admittances of a streamlined bridge deck at different angles of attack | |
CN113627048B (zh) | 基于局部传递率函数与模式匹配的结构损伤快速识别方法 | |
WO2019200953A1 (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 |