CN106168603B - 一种Lamb波结构健康监测中的温度补偿方法 - Google Patents
一种Lamb波结构健康监测中的温度补偿方法 Download PDFInfo
- Publication number
- CN106168603B CN106168603B CN201610529008.6A CN201610529008A CN106168603B CN 106168603 B CN106168603 B CN 106168603B CN 201610529008 A CN201610529008 A CN 201610529008A CN 106168603 B CN106168603 B CN 106168603B
- Authority
- CN
- China
- Prior art keywords
- signal
- temperature
- lamb wave
- reference signal
- compensation
- 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
Links
- 235000019687 Lamb Nutrition 0.000 title claims abstract description 68
- 238000012544 monitoring process Methods 0.000 title claims abstract description 34
- 230000036541 health Effects 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 3
- 239000002131 composite material Substances 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000006735 deficit Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000032798 delamination Effects 0.000 description 1
- 229920006335 epoxy glue Polymers 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/041—Analysing solids on the surface of the material, e.g. using Lamb, Rayleigh or shear waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/32—Arrangements for suppressing undesired influences, e.g. temperature or pressure variations, compensating for signal noise
- G01N29/326—Arrangements for suppressing undesired influences, e.g. temperature or pressure variations, compensating for signal noise compensating for temperature variations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开了一种Lamb波结构健康监测中的温度补偿方法。所述Lamb波结构健康监测中的温度补偿方法包括如下步骤:步骤1:采集一组Lamb波信号作为基准信号,并获取采集该基准信号时的温度以及在采集一组Lamb波信号作为参考信号,并获取采集参考信号时的温度;步骤2:采集当前Lamb波信号,并记录当前温度;步骤3:将基准信号的相位补偿到当前信号当前Lamb波信号,从而得到补偿后的当前Lamb波信号;步骤4:对进行正则化处理;步骤5:将基准信号补偿到当前信号的温度。本发明中的Lamb波结构健康监测中的温度补偿方法不需要采集大量基准信号,不需要确定任何温度补偿参数,处理过程简单,能对Lamb波信号中包含的全部波形进行温度补偿,实用性强。
Description
技术领域
本发明涉及结构健康监测技术领域,特别是涉及一种Lamb波结构健康监测中的温度补偿方法。
背景技术
近年来,基于Lamb波的结构健康监测技术受到了广泛的关注,成为了一个研究热点。这是由于Lamb波能够在结构中传播较长的距离,并且对裂纹、腐蚀和分层等损伤敏感。目前,Lamb波结构健康监测中的一种常用工作模式是激励—接收模式。这种工作模式是在被监测结构上布置一传感器阵列,然后利用一个传感器作激励器向结构中激励Lamb波,其它传感器作接收器接收激励的Lamb波,直到每一个传感器都作了一次激励器为止。当Lamb波结构健康监测技术应用于监测复杂结构时,其测量的信号一般非常复杂,难以进行分析解释。一种简单有效地方法是在监测前测量一组基准信号,然后将目前测量的当前信号与基准信号作相减,从而将此差信号作为损伤信号。但是外界环境,尤其是温度,会对测量信号产生很大影响。温度变化会造成基准信号相减产生误差。这些误差的强度往往和损伤信号基本相当。这样就会使损伤监测产生错误的结果。因此,在Lamb波结构健康监测中应该对温度影响进行补偿。
目前,公认为较有效地Lamb波结构健康监测的温度补偿方法为最优基准信号相减和基准信号压缩/扩展的结合。但基准信号压缩/扩展需要在温度补偿前确定信号相位和幅值的温度补偿参数。由于温度对直达Lamb波和反射 Lamb波造成的影响不同,因此对这些波的温度补偿参数也会不同。当被监测结构复杂时,温度补偿参数的确定也会变的非常复杂。且当被监测结构不同时,一般温度补偿参数需要重新确定。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本发明的目的在于提供一种Lamb波结构健康监测中的温度补偿方法来克服或至少减轻现有技术的至少一个上述缺陷。
为实现上述目的,本发明提供了一种Lamb波结构健康监测中的温度补偿方法,用于采用基准信号与测量信号之差作为损伤信号的Lamb波结构健康监测方法,所述Lamb波结构健康监测中的温度补偿方法包括如下步骤:步骤 1:在温度补偿前,采集一组Lamb波信号作为基准信号b(t),并获取采集该基准信号b(t)时的温度T0以及采集一组Lamb波信号作为参考信号r(t),并获取采集参考信号r(t)时的温度T1;步骤2:采集当前Lamb波信号s(t),并记录当前温度T2;步骤3:根据所述步骤1及所述步骤2的数据,通过公式将基准信号b(t)的相位补偿到当前Lamb波信号s(t),从而得到相位补偿后的基准信号步骤4:通过公式对相位补偿后的基准信号进行正则化处理;步骤 5:通过公式将相位补偿后的基准信号的幅值补偿到当前Lamb信号s(t),从而得到温度补偿后的基准信号。
优选地,所述温度T0与所述温度T1两温度之差小于6℃。
优选地,所述步骤3中的公式为:
其中,
为相位补偿后的基准信号;
—基准信号b(t)的希尔伯特变换;
Re{·}—取复数实部的运算;
j—虚数单位;
arg—信号的瞬时相位。
优选地,所述步骤4中的公式为:
其中,
—相位补偿后的基准信号;
—对 正则化后的信号。
优选地,所述步骤5中的公式为:
其中,
<·>—为内积运算;
btb(t)—为将基准信号b(t)补偿到当前信号s(t)温度的补偿信号。
本发明中的Lamb波结构健康监测中的温度补偿方法不需要采集大量基准信号,不需要确定任何温度补偿参数,处理过程简单,能对Lamb波信号中包含的全部波形进行温度补偿,实用性强,可应用于复杂结构的结构健康监测。
附图说明
图1是根据本发明一实施例的Lamb波结构健康监测中的温度补偿方法的流程示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明保护范围的限制。
图1是根据本发明一实施例的Lamb波结构健康监测中的温度补偿方法的流程示意图。
如图1所示的Lamb波结构健康监测中的温度补偿方法用于采用基准信号与测量信号之差作为损伤信号的Lamb波结构健康监测方法,该Lamb波结构健康监测中的温度补偿方法包括如下步骤:
步骤1:在温度补偿前,采集一组Lamb波信号作为基准信号b(t),并获取采集该基准信号b(t)时的温度T0以及采集一组Lamb波信号作为参考信号r(t),并获取采集参考信号r(t)时的温度T1;
步骤2:采集当前Lamb波信号s(t),并记录当前温度T2;
步骤3:根据所述步骤1及所述步骤2的数据,通过公式将基准信号b(t)的相位补偿到当前Lamb波信号s(t),从而得到相位补偿后的基准信号
步骤4:通过公式对相位补偿后的基准信号进行正则化处理;
步骤5:通过公式将相位补偿后的基准信号的幅值补偿到当前Lamb 信号s(t),从而得到温度补偿后的基准信号。
有利的是,在温度T0与温度T1两温度之差小于6℃的范围内,本申请的 Lamb波结构健康监测中的温度补偿方法最为准确。
在本实施例中,步骤3中的公式为:
其中,
为相位补偿后的基准信号;
—基准信号b(t)的希尔伯特变换;
Re{·}—取复数实部的运算;
j—虚数单位;
arg—信号的瞬时相位。
在本实施例中,步骤4中的公式为:
其中,
—相位补偿后的基准信号;
—对 正则化后的信号。
在本实施例中,所述步骤5中的公式为:
其中,
<·>—为内积运算;
btb(t)—为将基准信号b(t)补偿到当前信号s(t)温度的补偿信号。
在本实施例中,步骤2中的当前Lamb波信号s(t),及其他地方提及的当前Lamb波信号s(t)是指需要判定有无损伤时,采集的当下的信号。
在本实施例中,基准信号为结构没有破损的情况下的Lamb波信号。
在本实施例中,步骤3中的公式中的T1、T2、T0为步骤1中的各个信号的温度。
下面以举例的方式本申请进行进一步阐述,可以理解的是,该举例并不构成对本申请的任何限制。
对一碳纤维复合材料层板的Lamb波结构健康监测进行了温度补偿。复合材料层板的材料为T700/BA9916,尺寸为150mm×100mm×3mm,铺层为 [45/0/-45/90/0/45/0/-45/0/45/90/-45]S,每层铺层的厚度为0.125mm。利用环氧胶在板上粘贴了2个压电片,其中一个作为激励器,另一个作传感器。压电材料为P51,直径为8mm,厚度为0.45mm。复合材料层板的尺寸和压电片的布置。复合材料层板被置于温度环境箱中,通过改变温度环境箱中的温度来改变层板及周围的环境温度。
Lamb波的激励信号为调制的5波峰正弦信号。激励频率为70KHz,信号采样率为10MHz。
按照本申请的Lamb波结构健康监测中的温度补偿方法进行温度补偿具体为:在环境箱温度为20℃时,采集了基准信号b(t),并记录温度T0。接着以 1℃/min的速率升高环境箱温度至23℃,并保温1小时,然后采集参考信号r(t),并记录温度T1。
然后以1℃/min的速率升高环境箱温度至29℃,并保温1小时,然后采集当前Lamb波信号s(t),并记录温度T2。
通过公式将基准信号b(t)的相位补偿到当前信号当前Lamb波信号s(t),从而得到相位补偿后的基准信号
对相位补偿后的基准信号进行正则化处理。
将相位补偿后的基准信号的幅值补偿到当前Lamb信号s(t),从而得到温度补偿后的基准信号。
由于温度会对Lamb波信号的幅值和相位有影响,本申请将基准信号的相位和幅值,补偿为与当前信号具有相同温度时所具有的特征。从而形成一种 Lamb波结构健康监测中的温度补偿方法。
本发明的Lamb波结构健康监测中的温度补偿方法相比现有技术的优点在于:不需要采集大量基准信号,不需要确定任何温度补偿参数,处理过程简单,能对Lamb波信号中包含的全部波形进行温度补偿,实用性强,可应用于复杂结构的结构健康监测。
最后需要指出的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (2)
1.一种Lamb波结构健康监测中的温度补偿方法,用于采用基准信号与测量信号之差作为损伤信号的Lamb波结构健康监测方法,其特征在于,所述Lamb波结构健康监测中的温度补偿方法包括如下步骤:
步骤1:在温度补偿前,采集一组Lamb波信号作为基准信号b(t),并获取采集该基准信号b(t)时的温度T0以及采集一组Lamb波信号作为参考信号r(t),并获取采集参考信号r(t)时的温度T1;
步骤2:采集当前Lamb波信号s(t),并记录当前温度T2;
步骤3:根据所述步骤1及所述步骤2的数据,通过公式将基准信号b(t)的相位补偿到当前Lamb波信号s(t),从而得到相位补偿后的基准信号
步骤4:通过公式对相位补偿后的基准信号进行正则化处理;
步骤5:通过公式将相位补偿后的基准信号的幅值补偿到当前Lamb信号s(t),从而得到温度补偿后的基准信号;其中,
所述步骤3中的公式为:
其中,
为相位补偿后的基准信号;
—基准信号b(t)的希尔伯特变换;
Re{·}—取复数实部的运算;
j—虚数单位;
arg—信号的瞬时相位;
所述步骤4中的公式为:
其中,
—相位补偿后的基准信号;
—对 正则化后的信号;
所述步骤5中的公式为:
其中,
<·>—为内积运算;
btb(t)—为将基准信号b(t)补偿到当前信号s(t)温度的补偿信号。
2.如权利要求1所述的Lamb波结构健康监测中的温度补偿方法,其特征在于,所述温度T0与所述温度T1两温度之差小于6℃。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610529008.6A CN106168603B (zh) | 2016-07-05 | 2016-07-05 | 一种Lamb波结构健康监测中的温度补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610529008.6A CN106168603B (zh) | 2016-07-05 | 2016-07-05 | 一种Lamb波结构健康监测中的温度补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106168603A CN106168603A (zh) | 2016-11-30 |
CN106168603B true CN106168603B (zh) | 2019-04-23 |
Family
ID=58065911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610529008.6A Active CN106168603B (zh) | 2016-07-05 | 2016-07-05 | 一种Lamb波结构健康监测中的温度补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106168603B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107782787A (zh) * | 2017-10-20 | 2018-03-09 | 杭州电子科技大学 | 一种超声波缺陷检测方法 |
CN107748208B (zh) * | 2017-10-24 | 2019-07-02 | 厦门大学 | 一种基于基准导波信号匹配的温度补偿方法 |
CN107917957B (zh) * | 2017-10-27 | 2021-06-01 | 中车青岛四方机车车辆股份有限公司 | 一种板形结构的损伤检测方法 |
CN108663411A (zh) * | 2018-03-24 | 2018-10-16 | 北京工业大学 | 一种机电阻抗结构损伤监测的温度补偿方法 |
CN113567564B (zh) * | 2021-05-17 | 2022-07-26 | 南京航空航天大学 | 考虑幅值补偿的时域弯折超声导波大范围温度补偿方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0330077A (ja) * | 1989-06-27 | 1991-02-08 | Nec Corp | 温度補償型複素乗算器 |
CN1892255A (zh) * | 2005-07-04 | 2007-01-10 | 中国科学院声学研究所 | 一种对海底多金属结核矿进行测量的声学方法及系统 |
CN1977160A (zh) * | 2004-06-12 | 2007-06-06 | 阿库比奥有限公司 | 带有传感器阵列和校准元件的分析仪器 |
CN101009150A (zh) * | 2006-01-26 | 2007-08-01 | 徐炎堂 | 用于解析系统电路稳定化的频谱器 |
CN101135670A (zh) * | 2007-10-12 | 2008-03-05 | 电子科技大学 | 一种具有湿度补偿功能的双通道声表面波气体传感器 |
CN101135674A (zh) * | 2007-10-12 | 2008-03-05 | 电子科技大学 | 一种提高声表面波气体传感器温度稳定性的方法 |
CN102440762A (zh) * | 2011-09-13 | 2012-05-09 | 深圳沃夫特影像技术有限公司 | 干涉相位参数定标反馈调控稳频方法及其系统 |
CN104181237A (zh) * | 2014-04-04 | 2014-12-03 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种结构件探伤监测温度补偿方法及系统 |
-
2016
- 2016-07-05 CN CN201610529008.6A patent/CN106168603B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0330077A (ja) * | 1989-06-27 | 1991-02-08 | Nec Corp | 温度補償型複素乗算器 |
CN1977160A (zh) * | 2004-06-12 | 2007-06-06 | 阿库比奥有限公司 | 带有传感器阵列和校准元件的分析仪器 |
CN1892255A (zh) * | 2005-07-04 | 2007-01-10 | 中国科学院声学研究所 | 一种对海底多金属结核矿进行测量的声学方法及系统 |
CN101009150A (zh) * | 2006-01-26 | 2007-08-01 | 徐炎堂 | 用于解析系统电路稳定化的频谱器 |
CN101135670A (zh) * | 2007-10-12 | 2008-03-05 | 电子科技大学 | 一种具有湿度补偿功能的双通道声表面波气体传感器 |
CN101135674A (zh) * | 2007-10-12 | 2008-03-05 | 电子科技大学 | 一种提高声表面波气体传感器温度稳定性的方法 |
CN102440762A (zh) * | 2011-09-13 | 2012-05-09 | 深圳沃夫特影像技术有限公司 | 干涉相位参数定标反馈调控稳频方法及其系统 |
CN104181237A (zh) * | 2014-04-04 | 2014-12-03 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种结构件探伤监测温度补偿方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106168603A (zh) | 2016-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106168603B (zh) | 一种Lamb波结构健康监测中的温度补偿方法 | |
Kannusamy et al. | Accurate baseline-free damage localization in plates using refined Lamb wave time-reversal method | |
Bochud et al. | Towards real-time assessment of anisotropic plate properties using elastic guided waves | |
Zahedi et al. | Time–frequency analysis of electro-mechanical impedance (EMI) signature for physics-based damage detections using piezoelectric wafer active sensor (PWAS) | |
CN113030671B (zh) | 一种电气设备局部放电故障现场检测与定位装置及方法 | |
CN103134606B (zh) | 差分式声表面波温度传感器 | |
TW201009881A (en) | Passive capacitively-coupled electrostatic (CCE) probe arrangement for detecting plasma instabilities in a plasma processing chamber | |
CN104458834A (zh) | 一种材料缺陷检测设备以及材料缺陷检测方法 | |
CN103615995A (zh) | 一种基于超声表面波的薄熔覆层厚度无损评价方法 | |
CN107656015B (zh) | 一种碳纤维复合材料加固钢结构的端部脱胶检测装置及方法 | |
CN105388190A (zh) | 基于耦合机电阻抗的复合材料机翼损伤定位方法 | |
CN114235971A (zh) | 一种基于超声导波的编织复合材料板损伤形状识别方法 | |
Xu et al. | Sparse estimation of propagation distances in Lamb wave inspection | |
Jiang et al. | Capacitively-coupled dual ring antennas for bolt loosening detection | |
RU2015130842A (ru) | Способ выявления аномальных скачков порового давления на границах разделов в непробуренных геологических формациях и система для осуществления этого способа | |
Chacón et al. | A Review of Non-Destructive Techniques for Lithium-Ion Battery Performance Analysis | |
CN106483525A (zh) | 全向超声信号接收组件、全向超声测距系统及方法 | |
Fromme | Monitoring of plate structures using guided ultrasonic waves | |
CN110375989A (zh) | 一种柴油机燃烧噪声检测系统 | |
CN104280460B (zh) | 一种复合材料水浸超声波检测方法 | |
CN110895258A (zh) | 一种带有温度补偿功能的压电阻抗监测系统及方法 | |
Zhou et al. | Damage detection in plate structures based on baseline construction of ultrasonic Lamb wave using designed mobile transducer set | |
CN107271558B (zh) | 一种基于稀疏阵列层析成像的绝热层脱粘损伤监测方法 | |
Janarthan et al. | Damage detection in stiffened composite panels using Lamb wave | |
CN211426351U (zh) | 一种带有温度补偿功能的压电阻抗监测系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |