CN114034331B - 一种压电陶瓷圆管换能器的一致性控制方法 - Google Patents

一种压电陶瓷圆管换能器的一致性控制方法 Download PDF

Info

Publication number
CN114034331B
CN114034331B CN202111304653.5A CN202111304653A CN114034331B CN 114034331 B CN114034331 B CN 114034331B CN 202111304653 A CN202111304653 A CN 202111304653A CN 114034331 B CN114034331 B CN 114034331B
Authority
CN
China
Prior art keywords
transducer
piezoelectric ceramic
torque
calculating
limit
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
Application number
CN202111304653.5A
Other languages
English (en)
Other versions
CN114034331A (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.)
University of Jinan
Original Assignee
University of Jinan
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 University of Jinan filed Critical University of Jinan
Priority to CN202111304653.5A priority Critical patent/CN114034331B/zh
Publication of CN114034331A publication Critical patent/CN114034331A/zh
Application granted granted Critical
Publication of CN114034331B publication Critical patent/CN114034331B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D18/00—Testing or calibrating apparatus or arrangements provided for in groups G01D1/00 - G01D15/00

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Compositions Of Oxide Ceramics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供一种压电陶瓷圆管换能器的一致性控制方法，包括根据工作电压和压电陶瓷抗张强度确定最大力矩，通过谐振频率计算机电耦合系数，通过力矩调整机电耦合系数，进而对换能器一致性进行判断、比较和调整的方法。本发明利用螺栓紧固力矩、静态电容和机电耦合系数对压电陶瓷圆管换能器的一致性进行调控，具备方法简单、成本低、调控准确、适用范围大等优势。

Description

一种压电陶瓷圆管换能器的一致性控制方法
技术领域
本发明涉及超声波传感器领域，具体涉及一种压电陶瓷圆管换能器的一致性控制方法。
背景技术
压电陶瓷圆管换能器已被广泛应用于声波测井、基桩完整性检测等领域，作为超声系统的关键器件，压电陶瓷圆管换能器决定了超声系统性能的好坏。若换能器一致性较差，不利于换能器的配对使用，测量波形的幅度较小，进而严重影响测量精度。而机电耦合系数是评价换能器机械能和电能互换能力的的关键，保证机电耦合系数的一致性，对保证换能器性能的一致性至关重要。
目前，保障压电陶瓷换能器一致性的方法主要是一致性筛选，一方面是利用相关系数法、最大值法等对压电晶体的谐振频率进行筛选，获得幅相一致的用电晶体；另一方面是利用导纳圆等参数对压电陶瓷换能器进行一致性筛选，导纳圆可以通过阻抗分析仪或导纳仪获得。一致性筛选不可避免地存在人力物力和财力的浪费。因此，通过换能器制备工艺中的过程控制，从源头上提高换能器的一致性，就具有更高的经济价值和工程意义。
发明内容
针对已有的压电陶瓷一致性调控技术的不足，本发明提供一种压电陶瓷圆管换能器一致性的控制方法，根据换能器工作电压和压电陶瓷抗张极限，确定扭矩最大值，通过谐振频率计算机电耦合系数，通过力矩调整机电耦合系数，进而对换能器一致性进行判断、比较和调整，具备成本低、方法简单、调控准确、适用范围大等优势。
本发明是通过以下技术方案实现的：
一种压电陶瓷圆管换能器的一致性控制方法，包括以下步骤：
(1)计算压电陶瓷圆管换能器的径向工作应力PV；
(2)计算定径向应力安全极限PT；
(3)计算定量施加轴向应力的螺栓拧紧力矩极限Tmax；
(4)准备换能器原材料，压电陶瓷圆管表面完整，电极无缺损，垫片和金属螺栓的尺寸误差小于0.01mm；
(5)用台钳固定组装好的换能器，并用扭力扳手施加初始扭矩Ti；
(6)用阻抗分析仪分别测量n组换能器在扭矩Ti作用下的并联谐振频率fpi、串联谐振频率fsi和静态电容C0i；
(7)计算n组换能器测试参数的平均值，得到平均并联谐振频率fp、平均串联谐振频率fs和平均静态电容C，平均机电耦合系数k31n，满足
Figure BDA0003339717130000021
(8)取机电耦合系数一致性区间为(mk31n，nk31n)，其中m和n为经验常数，满足m＜l，n≥l；
(9)取待测换能器，若待测换能器静态电容与平均值偏差满足lCi-Cl/C＞a，可能存在陶瓷损伤、虚焊、短路、断路、导线表皮损伤等，需进行换能器检修，重复步骤(4)-(8)；若满足lCi-Cl/C≤a，继续下一步；
(10)计算待测换能器的机电耦合系数k31i，满足
Figure BDA0003339717130000022
(11)比较k31i是否在(mk31n，nk31n)区间内：若k31i介于(mk31n，nk31n)区间，直接进行换能器灌封；否则，继续下一步；
(12)若k31i＞nk31n，扭矩Ti需降低为Ti+1，若k31i＜mk31n，扭矩Ti需要增大为Ti+1；
(13)比较Ti+1是否在(0，Tmax)区间内：若Ti+1介于有效扭矩区间(0，Tmax)，则在换能器上重新施加扭矩Ti+1，重复步骤(10)-(11)；否则，继续下一步；
(14)更换压电陶瓷圆管，重新组装换能器，重复步骤(5)-(11)。
优选地，所述步骤(1)中径向工作应力PV的确定方法为：
①通过压电陶瓷材料确定压电常数d31、杨氏模量
Figure BDA0003339717130000024
泊松比σ；
②测量压电陶瓷圆管的厚度t；
③确定换能器的工作电压V；
④计算圆管换能器径向工作应力PV：
Figure BDA0003339717130000023
优选地，所述步骤(2)中径向应力安全极限的确定方法为：
①通过使用类型压电陶瓷材料，确定抗张极限PC；
②计算圆管换能器径向工作应力限安全极限PT，满足：PT＝PC-PV。
优选地，所述步骤(3)中螺栓拧紧力矩极限Tmax的确定方法为：
Tmax＝α+βPT+γPT 2，其中α、β和γ为与压电陶瓷材料和换能器结构相关的经验常数，取值分别为α＝0.4±0.1，β＝0.04±0.01，γ＝0.0017±0.0001。
优选地，所述步骤(5)中扭力扳手施加初始扭矩Ti的确定方法为Ti＝κTmax，κ典型取值范围为0.4～0.8。
优选地，所述步骤(8)中静态电容与平均值偏差值由换能器设计需求允许频率误差许可范围决定，可取a≤10％。
有益效果：
本发明公开了一种压电陶瓷圆管换能器一致性的控制方法，具有以下优点：
(1)对设备要求低，扭矩调节和参数测试方式简单；
(2)对压电陶瓷圆管使用效率最大化，降低原料成本；
(3)可根据实时谐振频率计算机电耦合系数，准确性高；
(4)适用范围大，适用于不同尺寸和材料的压电陶瓷圆管换能器。
附图说明
图1为本发明所述的压电陶瓷圆管换能器结构示意图；
图2为本发明所述的在换能器上施加力矩及测试阻抗曲线的现场布置示意图；
图3为利用本发明实施案例1中6个抽样样本在初始扭矩下的阻抗曲线及参数变化；
图4为利用本发明实施案例1中2#换能器根据调控扭矩前后的阻抗曲线及参数变化；
图5为利用本发明实施案例1中利用无机非金属检测仪测量首波电压的现场布置图；
图6为利用本发明实施案例1中利用无机非金属检测仪测量的首波电压曲线及首波参数变化；
图7为本发明流程图。
其中1-压电陶瓷圆管、2-定位垫片、3-左端帽、4-右端帽、5-螺杆、6-固定平台、7-导线、8-扭力扳手、9-台钳、10-阻抗分析仪、11-无机非金属检测仪、12-50m线缆、13-标准换能器、14-自制换能器、15-水、16-声测管、17-混凝土。
具体实施方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
参照附图，本申请的使用的螺栓结构包括压电陶瓷圆管l、套在压电陶瓷圆管1两端的左端帽3和右端帽4、压电陶瓷圆管1与端帽之间的定位垫片2，压电陶瓷圆管1穿过螺杆5，当采用多个压电陶瓷圆管时，始终满足相邻压电陶瓷圆管之间有定位垫片进行电气隔离。
左端帽3、右端帽4和螺杆5为不锈钢材料，左端帽3和右端帽4外侧预留固定平台6，便于台钳和扳手的固定。定位垫片2为尼龙材料，在沿半径的方向固定压电陶瓷圆管。定位垫片结构简单，便于加工，可以有效实现不同尺寸压电圆管的固定夹持。
一种压电陶瓷圆管换能器一致性的控制方法，包括以下步骤：
(1)计算压电陶瓷圆管换能器的径向工作应力PV：
①选用淄博宇海电子陶瓷有限公司生产的PZT-52压电陶瓷，通过厂家提供的资料查询压电陶瓷材料的压电常数d31为204*10-12m/V，杨氏模量
Figure BDA0003339717130000041
为17*10-12m2/N，泊松比σ为0.36；
②测量压电陶瓷圆管的厚度t为1.5mm；
③换能器应用的最大工作电压V为1000V；
④根据
Figure BDA0003339717130000042
计算圆管换能器径向工作应力PV为22.2MPa。
(2)计算定径向应力安全极限PT：通过陶瓷生产厂家提供的资料或PZT-52压电陶瓷材料的抗张极限PC为60MPa，因此根据PT＝PC-PV计算圆管换能器径向工作应力限安全极限为37.78MPa；
(3)根据Tmax＝0.4+0.04PT+0.0017PT 2，计算定量施加轴向应力的螺栓拧紧力矩极限Tmax为4.3N·m；
(4)准备换能器原材料，压电陶瓷圆管尺寸为
Figure BDA0003339717130000043
表面清洁，电极无缺损；压电陶瓷圆管、垫片和金属螺栓的尺寸误差小于0.01mm；
(5)用台钳固定组装好的换能器，并根据Ti＝0.6Tmax用扭力扳手施加初始扭矩Ti＝2.6N·m；
(6)用阻抗分析仪(Key sight，4900A)，分别测量6组换能器在扭矩Ti作用下的并联谐振频率fpi、串联谐振频率fsi和静态电容C0i，如附图3所示；
(7)计算6组换能器测试参数的平均值，得到平均并联谐振频率fp为46.20kHz、平均串联谐振频率fs为48.69kHz和平均静态电容C为22.29nF，根据
Figure BDA0003339717130000051
计算6组换能器的平均机电耦合系数k31n为0.32；
(8)取经验系数m为0.9，n为1.1，则机电耦合系数一致性区间为(0.288，0.352)；
(9)取待测换能器，静态电容Ci与平均值进行比较，所有换能器偏差满足|Ci-C|/C＜2％，均无虚焊或陶瓷损伤；
(10)根据
Figure BDA0003339717130000052
计算待测换能器的机电耦合系数k31i，如附图3所示；
(11)取待测换能器进行k31i比较，1#、3#-6#样品的k31i均位于置信区间内，可以直接进行换能器灌封，现2#样品的k31i＜0.288，需进行下一步；扭矩需要增大为Ti+1；
(12)取Ti+1为3N·m，Ti+1介于有效扭矩区间(0，4.3)内，则在2号换能器上重新施加3N·m扭矩，重复步骤9，静态电容为22.5nF，在偏差允许范围内，计算机电耦合系数为0.34，位于机电耦合系数一致性区间(0.288，0.352)内，可以进行换能器灌封；
(13)待测换能器在混凝土模型中进行声波跨孔测试。检测设备采用无机非金属检测仪(武汉中岩，RSM-SY6)，设备配套的换能器作为标准发射声源。声测管平均直径为38mm，使用前注满清水。待测换能器和标准声源通过50m线缆连接到检测设备上，然后分别悬垂于声测管中，换能器水平间距为135mm。6个自制换能器测得的首波电压平均值为65.51mV，偏差满足|Vi-V|/V＜3.4％，一致性较好。
本发明未尽事宜为公知技术。
上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种压电陶瓷圆管换能器的一致性控制方法，其特征在于，包括以下步骤：
（1）计算压电陶瓷圆管换能器的径向工作应力
Figure 897447DEST_PATH_IMAGE001
；
（2）计算定径向应力安全极限
Figure 412742DEST_PATH_IMAGE002
；
（3）计算定量施加轴向应力的螺栓拧紧力矩极限
Figure 866857DEST_PATH_IMAGE003
；
（4）准备换能器原材料，压电陶瓷圆管表面完整，电极无缺损，垫片和金属螺栓的尺寸误差小于0.01mm；
（5）固定组装好的换能器，并用扭力扳手施加初始扭矩
Figure 770091DEST_PATH_IMAGE004
；
（6）用阻抗分析仪分别测量n组换能器在扭矩
Figure 44078DEST_PATH_IMAGE004
作用下的并联谐振频率
Figure 441692DEST_PATH_IMAGE006
、串联谐振频率
Figure 547051DEST_PATH_IMAGE008
和静态电容
Figure 824449DEST_PATH_IMAGE009
；
（7）计算n组换能器测试参数的平均值，得到平均并联谐振频率
Figure 382469DEST_PATH_IMAGE011
、平均串联谐振频率
Figure 708408DEST_PATH_IMAGE013
和平均静态电容
Figure 481323DEST_PATH_IMAGE014
，平均机电耦合系数k31n，满足
Figure 601726DEST_PATH_IMAGE015
；
（8）取机电耦合系数一致性区间为（
Figure 240518DEST_PATH_IMAGE016
，
Figure DEST_PATH_IMAGE017
）,其中m和n为经验常数，满足m＜1，n≥1；
（9）取待测换能器，若待测换能器静态电容与平均值偏差满足
Figure 370148DEST_PATH_IMAGE018
，存在陶瓷损伤、虚焊、短路、断路、导线表皮损伤，需进行换能器检修，重复步骤（4）-（8）；若满足
Figure 263149DEST_PATH_IMAGE019
，继续下一步；
（10）计算待测换能器的机电耦合系数
Figure 554453DEST_PATH_IMAGE020
，满足
Figure 414961DEST_PATH_IMAGE021
；
（11）比较
Figure 879441DEST_PATH_IMAGE020
是否在（
Figure 751582DEST_PATH_IMAGE016
，
Figure 354733DEST_PATH_IMAGE017
）区间内：若
Figure 577904DEST_PATH_IMAGE020
介于（
Figure 908391DEST_PATH_IMAGE016
，
Figure 635038DEST_PATH_IMAGE017
）区间，直接进行换能器灌封；否则，继续下一步；
（12）若
Figure 409090DEST_PATH_IMAGE022
，扭矩
Figure 853978DEST_PATH_IMAGE004
需降低为
Figure 660260DEST_PATH_IMAGE023
，若
Figure 897206DEST_PATH_IMAGE024
，扭矩
Figure 701214DEST_PATH_IMAGE004
需要增大为
Figure 774344DEST_PATH_IMAGE023
；
（13）比较
Figure 322000DEST_PATH_IMAGE023
是否在（0，
Figure 413453DEST_PATH_IMAGE003
）区间内：若
Figure 388362DEST_PATH_IMAGE023
介于有效扭矩区间（0，
Figure 942928DEST_PATH_IMAGE003
），则在换能器上重新施加扭矩
Figure 294275DEST_PATH_IMAGE023
，重复步骤（10）-（11）；否则，继续下一步；
（14）更换压电陶瓷圆管，重新组装换能器，重复步骤（5）-（11）。
2.根据权利要求1所述的方法，其特征在于，所述步骤（1）中径向工作应力
Figure 974655DEST_PATH_IMAGE001
的确定方法为：
①通过压电陶瓷材料确定压电常数d31、杨氏模量
Figure 386045DEST_PATH_IMAGE025
、泊松比
Figure 168187DEST_PATH_IMAGE026
；
②测量压电陶瓷圆管的厚度t；
③确定换能器的工作电压V；
④计算圆管换能器径向工作应力
Figure 854383DEST_PATH_IMAGE001
：
Figure 795794DEST_PATH_IMAGE027
。
3.根据权利要求1所述的方法，其特征在于，所述步骤（2）中径向应力安全极限的确定方法为：
①通过压电陶瓷材料确定抗张极限
Figure 237140DEST_PATH_IMAGE028
；
②计算圆管换能器径向工作应力限安全极限
Figure 631212DEST_PATH_IMAGE002
，满足：
Figure 730066DEST_PATH_IMAGE030
。
4.根据权利要求1所述的方法，其特征在于，所述步骤（3）中螺栓拧紧力矩极限
Figure 607893DEST_PATH_IMAGE031
的确定方法为：
Figure 285999DEST_PATH_IMAGE032
，其中
Figure 517260DEST_PATH_IMAGE033
、
Figure 43050DEST_PATH_IMAGE034
和
Figure 701565DEST_PATH_IMAGE035
为与压电陶瓷材料和换能器结构的经验常数，取值分别为
Figure 194863DEST_PATH_IMAGE036
=0.4±0.1，
Figure 229815DEST_PATH_IMAGE037
=0.04±0.01，
Figure 875691DEST_PATH_IMAGE038
=0.0017±0.0001。
5.根据权利要求1所述的方法，其特征在于，所述的步骤（5）中扭力扳手施加初始扭矩
Figure 970686DEST_PATH_IMAGE039
的确定方法为
Figure 685701DEST_PATH_IMAGE040
，
Figure DEST_PATH_IMAGE041
典型取值范围为0.4~0.8。
6.根据权利要求1所述的方法，其特征在于，所述步骤（8）中静态电容与平均值偏差值由换能器设计需求允许频率误差许可范围决定，a≤10%。
CN202111304653.5A 2021-11-05 2021-11-05 一种压电陶瓷圆管换能器的一致性控制方法 Active CN114034331B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111304653.5A CN114034331B (zh) 2021-11-05 2021-11-05 一种压电陶瓷圆管换能器的一致性控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111304653.5A CN114034331B (zh) 2021-11-05 2021-11-05 一种压电陶瓷圆管换能器的一致性控制方法

Publications (2)

Publication Number Publication Date
CN114034331A CN114034331A (zh) 2022-02-11
CN114034331B true CN114034331B (zh) 2023-04-11

Family

ID=80136315

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111304653.5A Active CN114034331B (zh) 2021-11-05 2021-11-05 一种压电陶瓷圆管换能器的一致性控制方法

Country Status (1)

Country Link
CN (1) CN114034331B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3954015A (en) * 1975-07-17 1976-05-04 The United States Of America As Represented By The Secretary Of The Navy Method of determining piezoelectric constants of ceramic rings
CN102495299A (zh) * 2011-11-16 2012-06-13 中国科学院上海硅酸盐研究所 一种压电材料性能自动测试系统及方法
CN104535863A (zh) * 2014-12-23 2015-04-22 上海电机学院 一种压电特性参数动态扫频测试装置及方法
CN106990170A (zh) * 2017-04-19 2017-07-28 南京邮电大学 一种压电阵列Lamb波损伤监测中的信号一致性补偿方法
CN208313487U (zh) * 2018-05-31 2019-01-01 河南理工大学 夹心式纵振换能器装配在线监测及施加预紧力装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3954015A (en) * 1975-07-17 1976-05-04 The United States Of America As Represented By The Secretary Of The Navy Method of determining piezoelectric constants of ceramic rings
CN102495299A (zh) * 2011-11-16 2012-06-13 中国科学院上海硅酸盐研究所 一种压电材料性能自动测试系统及方法
CN104535863A (zh) * 2014-12-23 2015-04-22 上海电机学院 一种压电特性参数动态扫频测试装置及方法
CN106990170A (zh) * 2017-04-19 2017-07-28 南京邮电大学 一种压电阵列Lamb波损伤监测中的信号一致性补偿方法
CN208313487U (zh) * 2018-05-31 2019-01-01 河南理工大学 夹心式纵振换能器装配在线监测及施加预紧力装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
聂建华等.压电换能器一致性筛选技术.石油仪器.2009,第23卷(第06期),第7-9,12页. *

Also Published As

Publication number Publication date
CN114034331A (zh) 2022-02-11

Similar Documents

Publication Publication Date Title
Jiang et al. Detection of debonding between fiber reinforced polymer bar and concrete structure using piezoceramic transducers and wavelet packet analysis
US20080022773A1 (en) Ultrasonic material monitor for determining a characteristic of the material
Zhang et al. Shear flexoelectric coefficient μ1211 in polyvinylidene fluoride
Li et al. A novel acoustic emission monitoring method of cross‐section precise localization of defects and wire breaking of parallel wire bundle
US20230228718A1 (en) Non-baseline On-line Stress Monitoring System and Monitoring Method Based on Multi-mode Lamb Wave Data Fusion
Volkov et al. Second-and third-order elastic coefficients in polycrystalline aluminum alloy AMg6
CN111077030A (zh) 一种高应变率下混凝土动态力学性能试验装置及方法
CN110243521A (zh) 一种薄板应力测量方法及薄板应力测量系统
Liu et al. Relation between the shear stress distribution and the resulting acoustic emission variation in concrete beams
CN114034331B (zh) 一种压电陶瓷圆管换能器的一致性控制方法
Hosten et al. Identification of composite materials elastic moduli from Lamb wave velocities measured with single sided, contactless ultrasonic method
CN109283079B (zh) 一种测定及计算岩石衰减系数及非线性系数的方法
CN110988136B (zh) 一种固体材料单裂纹非线性超声系数表征方法
US20230048646A1 (en) Method for measuring deviation angle of fatigue microcrack based on nonlinear ultrasound
Bragov et al. Determination of the mechanical properties of concrete using the split Hopkinson pressure bar method
CN107894443A (zh) 一种用于监测锚固松紧的压电陶瓷系统及监测方法
CN114216583A (zh) 一种基于sh导波的温度应力在线监测系统及其监测方法
Korobov et al. Influence of a Static Reversible Loading on Mechanical and Elastic Properties of Polycrystalline Aluminum Alloy AMg6
Lollioz et al. Application of piezo-magnetism for the measurement of stress during an impact
Carpinteri et al. Structural monitoring and diagnostics by the acoustic emission technique: scaling of dissipated energy in compression
Xu et al. Performance-based Piezoelectric Material Applied in Health Inspection of Civil Engineering Structures
Wang et al. Application of nonlinear ultrasonic technique to characterize the damage evolution in structural steel after tensile deformation
CN115468866B (zh) 压电材料的霍普金森一维动态压缩力电特性的试验方法
Xiang et al. Damage Detection of an Aluminum Plate by Using Nonlinear Ultrasound with a Frequency-swept Excitation
Lu et al. Guided-wave-based detections of weld and crack in steel plates

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