CN111596549B - 一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 - Google Patents
一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 Download PDFInfo
- Publication number
- CN111596549B CN111596549B CN202010557741.5A CN202010557741A CN111596549B CN 111596549 B CN111596549 B CN 111596549B CN 202010557741 A CN202010557741 A CN 202010557741A CN 111596549 B CN111596549 B CN 111596549B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- control
- adaptive
- control voltage
- processing
- 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
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P70/00—Climate change mitigation technologies in the production process for final industrial or consumer products
- Y02P70/10—Greenhouse gas [GHG] capture, material saving, heat recovery or other energy efficient measures, e.g. motor control, characterised by manufacturing processes, e.g. for rolling metal or metal working
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法。本发明首先,设计前馈补偿的方式,解决换能器温度、匹配电感温度的阻抗变化导致振幅变动问题。然后,由于超声振幅与谐振频率存在强耦合关联,为了保证系统的振幅控制精度,消除系统参数的不确定性的影响,提高系统的抗干扰性,设计了一种鲁棒自适应控制器在线调整自适应控制率,确保系统参数不确定时的振幅控制精度。本发明可解决在复杂加工环境下,由于加工过程中由于负载、工艺及环境的变化以及换能器温度、匹配电感温度等原因导致振幅波动的问题,设计的强耦合鲁棒自适应振幅控制实现振幅的高精度稳定输出,保证硬脆材料在实际加工过程中的加工精度、加工工艺要求以及加工成品率。
Description
技术领域
本发明涉及硬脆材料超声加工领域,涉及一种应用于硬脆材料超声加工过程中强耦合鲁棒自适应振幅控制方法。
背景技术
随着5G时代的来临,手机制造商正紧随运营商开始切入5G市场,2019年,华为、苹果、三星等手机大厂相继宣布即将发布5G手机,因为5G的通信频段覆盖了3.3GHz-4.2GHz、4.4GHz-5.0GHz的低频段以及26GHz/28GHz/39GHz的毫米波高频段。5G信号对金属屏蔽非常敏感,尤其是对于新一代5G手机,金属外壳已不适用于高频毫米波段,对手机外壳的信号传输能力提出了更高的要求。
陶瓷、蓝宝石、玻璃和复合材料等新型材料手机外壳材料具有良好的质感,优良的耐磨性能和散热性能,并且能很好地满足5G通信和无线充电技术对机身材料的要求。然而蓝宝石、陶瓷、玻璃和复合材料等具有硬度高、强度大、耐磨性好、脆性大、易碎等特点,导致加工困难、加工质量差、加工效率低下、刀具磨损严重等问题。超声加工利用高频振动的碰撞作用实现材料微去除,由于作用力较小,因此加工质量好,表面损伤小,是目前加工硬脆材料的最有效手段之一。
在硬脆材料超声加工过程中的振幅控制是决定硬脆材料加工效果的最关键因素,但是在实际加工过程中由于负载、工艺及环境的变化,会导致振动系统阻抗特性和谐振频率的变化,从而影响超声振幅。同时由于换能器温度、匹配电感温度等均会引起阻抗发生变化,也会导致振幅波动。另外,在硬脆材料的加工过程中,超声振幅的控制精度往往要达到0.1um级别。然而在现有的超声振幅控制方法中往往难以应对实际加工过程中的复杂环境以及多影响因素的变化,同时对于高精度的振幅控制要求难以达到稳定振幅控制。
发明内容
为解决在实际加工过程中由于负载、工艺及环境的变化,会导致振动系统阻抗特性和谐振频率的变化的问题,本发明提供了一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法。
本发明首先设计前馈补偿的方式,解决换能器温度、匹配电感温度的阻抗变化导致振幅变动问题;由于超声振幅与谐振频率存在强耦合关联,为了保证系统的振幅控制精度,消除系统参数的不确定性的影响,提高系统的抗干扰性,然后设计了一种鲁棒自适应控制器在线调整自适应控制率,确保系统参数不确定时的振幅控制精度;最终实现鲁棒自适应控制与前馈控制相结合的方法进行强耦合鲁棒自适应振幅控制。
本发明解决技术问题所采取的技术方案为:
本发明包括以下步骤:
步骤一:建立振幅控制系统模型
所建立的振幅控制系统模型如下:
x2=f(x1,u)+d(x1,t)=f(x1,t)+g(x1,t)u+d(x1,t)
y=T(x2,t)
其中:x1=[x11,x12,...,x1n]T为振动系统内部可测状态变量矩阵,包括电源驱动电压、移相角,补偿电容,x2=[x21,x22,...,x2k]为声音传感器及压电加速度传感器的可测状态变量矩阵;u为系统控制输入的电压参数;y为系统振幅输出,f(x1,t),g(x1,t)为非线性连续函数;d(x1,t)为系统的外部扰动,T(x2,t)为根据声音传感器及压电加速度传感器的数据得到振幅的计算函数。最终的控制目的是使系统振幅输出y(t)跟踪参考振幅r(t)。
步骤二:振幅软测量
在刀柄处安装压电加速度传感器和声音传感器,分别测量刀柄振动的加速度和超声磨削加工引起的声音信号,从而实时得到超声波刀柄振幅。
步骤三:设计前馈补偿器
采用RBF神经网络离线建立谐振频率、换能器温度与影响振幅的控制电压变化量的关系模型,通过对大量离线数据采用神经网络进行训练,逼近得到影响振幅的控制电压变量与谐振频率、换能器温度的非线性关系函数。
up(t)=h(x3,t)
其中x3=[x31,x32,...,x3m]为可测向量,包括换能器温度、谐振频率,up(t)为前馈补偿器得到的控制电压变量,h(x3,t)为RBF神经网络逼近得到的非线性函数。
步骤四:建立鲁棒自适应控制器模型
由步骤一中振幅控制系统模型进一步得到如下的非线性系统:
x2 (n)=[fn(x1)+Δf(x1,t)]+[gn(x1)+Δg(x1,t)]u+d(x1,t)
=fn(x1)+gn(x1)u+z(x1,t)+d(x1,t)
=fn(x1)+gn(x1)u+zn(x1)+zd(x1,t)+d(x1,t)
其中fn(x1)和gn(x1)分别为f(x1,t)和g(x1,t)的标称部分;Δf(x1,t)和Δg(x1,t)分别为f(x1,t)和g(x1,t)的不确定变动部分;
令z(x1,t)=Δf(x1,t)+Δg(x1,t)u=zn(x1)+zd(x1,t)为系统未知不确定部分,zn(x1)为z(x1)的标称部分,zd(x1,t)为z(x1)的不确定变动部分,|zd(x1,t)|≤εzd,εzd为给定常数。此时振幅控制系统问题是求得控制电压变量u使得系统的输出振幅y(t)跟踪参考振幅r(t)。
步骤五:构建滑模面:
定义系统的跟踪误差为:
则滑模面为
其中c=[c1,c2,…,cn-1,1]T为选择的相应的特征多项式。
步骤六:确定自适应学习率
采用RBF神经网络来逼近非线性函数fn(x),gn(x),zn(x),即
步骤七:确定RBF神经网络各项控制率
采用自适应PID控制,则采用自适应滑模PID控制得到的控制电压变量为:
其中ks为滑模控制系数且ks≥εzd;r(n)为参考振幅的n状态;若选取控制电压变量为u=uf+us+un+uPID,则uf,us,un,uPID分别为:
为保证系统跟踪误差渐进收敛,则自适应学习率为:
步骤八:调节PWM实现振幅稳定控制
将鲁棒自适应控制得到的控制电压变量与前馈控制器得到的控制电压变量进行叠加得到最终的控制电压变量u(t)+up(t),并通过PWM调节电压占空比进而实现振动系统的振幅稳定控制。
本发明的有益效果:
超声发生器输出占空比、输出频率以及输入电压幅值是影响超声振动系统振幅的影响因素,由于占空比调节实现较容易且具有一定的线性调节能力,通过调节超声发生器输出占空比的方法来实现对超声振幅的控制。
面对由于在新型材料加工过程中,振动系统振幅控制往往存在着不确定性,如参数不确定性、结构不确定性以及外界干扰的不确定性。采用鲁棒自适应控制与前馈控制相结合的方法进行强耦合鲁棒自适应振幅控制。其中前馈控制用于补偿谐振频率及换能器温度变化引起的振幅变化,鲁棒自适应控制用于解决参数不确定、干扰不确定等情况下的振幅稳定控制。
最终解决在复杂加工环境下,由于加工过程中由于负载、工艺及环境的变化以及换能器温度、匹配电感温度等原因导致振幅波动的问题,设计的强耦合鲁棒自适应振幅控制实现振幅的高精度稳定输出,保证硬脆材料在实际加工过程中的加工精度、加工工艺要求以及加工成品率。
附图说明
图1超声振幅软测量模型;
图2RBF神经网络示意图;
图3强耦合鲁棒自适应控制器系统框图。
具体实施方式
以下结合附图对本发明作进一步描述。
本发明包括以下步骤:
步骤一:建立振幅控制系统模型:
硬脆材料超声加工的振幅控制系统为多输入单输出系统,建立反馈线性化控制模型抵消系统的非线性,建立如下系统模型:
其中:x1=[x11,x12,...,x1n]T为振动系统内部电源驱动电压、移相角,补偿电容等可测状态变量矩阵,x2=[x21,x22,...,x2k]为声音传感器及压电加速度传感器可测的电压、频率等可测状态变量矩阵;u为系统控制输入的电压参数;y为系统振幅输出,f(x1,t),g(x1,t)为非线性连续函数;d(x1,t)为系统的外部扰动,T(x2,t)为根据声音传感器及压电加速度传感器的数据得到振幅的计算函数。最终的控制目的是使系统的输出振幅y(t)跟踪参考振幅r(t)。
步骤二:振幅软测量模型设计
由于在加工现场环境复杂,现有的一些超声振幅测量方法如激光测量、显微镜测量等一些方法很难应用到加工现场,且大多数方法无法做到实时测量。并且由于超声加工过程中存在多种噪声信号,为此设计采用多源传感器数据融合的方法进行超声振幅软测量。如附图1所示,在刀柄处安装压电加速度传感器和声音传感器,分别测量刀柄振动的加速度和超声磨削加工引起的声音信号并与换能器端电信号相结合,进行超声波刀柄振幅实时测量。并将该振幅信号实时反馈给振幅控制系统。
步骤三:前馈补偿器设计:
采用RBF神经网络离线建立建立谐振频率、换能器温度与影响振幅的控制电压变化量的关系模型,通过对大量离线数据采用神经网络进行训练,逼近得到影响振幅的控制电压量与谐振频率、换能器温度的非线性关系函数。RBF神经网络示意图如附图2所示。
up(t)=h(x3,t) (2)
其中x3=[x31,x32,...,x3m]为换能器温度,谐振频率等外部可测向量,up(t)为前馈补偿器得到的控制电压变量,h(x3,t)为RBF神经网络逼近得到的非线性函数。
步骤四:建立鲁棒自适应控制器模型
由公式(1)确定的振幅控制系统模型可以进一步得到如下的非线性系统:
其中fn(x1)和gn(x1)分别为f(x1,t)和g(x1,t)的标称部分;Δf(x1,t)和Δg(x1,t)分别为f(x1,t)和g(x1,t)的不确定变动部分。
令z(x1,t)=Δf(x1,t)+Δg(x1,t)u=zn(x1)+zd(x1,t)为系统未知不确定部分,zn(x1)为z(x1)的标称部分,zd(x1,t)为z(x1)的不确定变动部分,|zd(x1,t)|≤εzd,εzd为给定常数。此时振幅控制系统问题是求得控制输入电压参数u使得系统的输出振幅y(t)跟踪参考振幅r(t)。
步骤五:构建滑模面:
定义系统的跟踪误差为:
则滑模面为
其中c=[c1,c2,…,cn-1,1]T为选择的相应的特征多项式。
步骤六:确定自适应学习率:
采用RBF神经网络来逼近非线性函数fn(x),gn(x),zn(x),即
步骤七:确定RBF神经网络各项控制率:
为提高系统系统,采用自适应PID控制,则采用自适应滑模PID控制得到的控制电压变量为:
若选取控制电压变量为u=uf+us+un+uPID则uf,us,un,uPID分别为:
为保证系统跟踪误差渐进收敛,则自适应学习率为:
步骤八:调节PWM实现振幅稳定控制:
将鲁棒自适应控制得到的控制电压变量与前馈控制器得到的控制电压变量进行叠加得到最终的控制电压变量u(t)+up(t),并通过PWM调节电压占空比进而实现振动系统的振幅稳定控制,强耦合鲁棒自适应控制器系统框图如附图3所示。
Claims (1)
1.一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法,包括以下步骤:
步骤一:建立振幅控制系统模型
所建立的振幅控制系统模型如下:
x2=f(x1,u)+d(x1,t)=f(x1,t)+g(x1,t)u+d(x1,t)
y=T(x2,t)
其中:x1=[x11,x12,...,x1n]T为振动系统内部可测状态变量矩阵,其中的变量包括电源驱动电压、移相角,补偿电容,x2=[x21,x22,...,x2k]为声音传感器及压电加速度传感器的可测状态变量矩阵;u为系统控制输入的电压参数;y为系统振幅输出,f(x1,t),g(x1,t)为非线性连续函数;d(x1,t)为系统的外部扰动,T(x2,t)为根据声音传感器及压电加速度传感器的数据得到的振幅计算函数;
最终的控制目的是使系统振幅输出y(t)跟踪参考振幅r(t);
步骤二:振幅软测量
在刀柄处安装压电加速度传感器和声音传感器,分别测量刀柄振动的加速度和超声磨削加工引起的声音信号,从而实时得到超声波刀柄振幅;
步骤三:设计前馈补偿器
采用RBF神经网络离线建立谐振频率、换能器温度与影响振幅的控制电压变化量的关系模型,通过对大量离线数据采用神经网络进行训练,逼近得到影响振幅的控制电压变量与谐振频率、换能器温度的非线性关系函数;
up(t)=h(x3,t)
其中x3=[x31,x32,...,x3m]为可测向量,包括换能器温度可测向量、谐振频率可测向量,up(t)为前馈补偿器得到的控制电压变量,h(x3,t)为RBF神经网络逼近得到的非线性函数;
步骤四:建立鲁棒自适应控制器模型
由步骤一中振幅控制系统模型进一步得到如下的非线性系统:
x2 (n)=[fn(x1)+Δf(x1,t)]+[gn(x1)+Δg(x1,t)]u+d(x1,t)
=fn(x1)+gn(x1)u+z(x1,t)+d(x1,t)
=fn(x1)+gn(x1)u+zn(x1)+zd(x1,t)+d(x1,t)
其中fn(x1)和gn(x1)分别为f(x1,t)和g(x1,t)的标称部分;Δf(x1,t)和Δg(x1,t)分别为f(x1,t)和g(x1,t)的不确定变动部分;
令z(x1,t)=Δf(x1,t)+Δg(x1,t)u=zn(x1)+zd(x1,t)为系统未知不确定部分,zn(x1)为z(x1)的标称部分,zd(x1,t)为z(x1)的不确定变动部分,|zd(x1,t)|≤εzd,εzd为给定常数;
此时振幅控制系统问题是求得控制电压变量u使得系统的输出振幅y(t)跟踪参考振幅r(t);
步骤五:构建滑模面:
定义系统的跟踪误差为:
则滑模面为
其中c=[c1,c2,…,cn-1,1]T为选择的相应的特征多项式;
步骤六:确定自适应学习率
采用RBF神经网络来逼近非线性函数fn(x),gn(x),zn(x),即
步骤七:确定RBF神经网络各项控制率
采用自适应PID控制,则采用自适应滑模PID控制得到的控制电压变量为:
其中ks为滑模控制系数且ks≥εzd;r(n)为参考振幅的n状态;
若选取控制电压变量为u=uf+us+un+uPID,则uf,us,un,uPID分别为:
为保证系统跟踪误差渐进收敛,则自适应学习率为:
步骤八:调节PWM实现振幅稳定控制
将鲁棒自适应控制得到的控制电压变量与前馈控制器得到的控制电压变量进行叠加得到最终的控制电压变量u(t)+up(t),并通过PWM调节电压占空比进而实现振动系统的振幅稳定控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010557741.5A CN111596549B (zh) | 2020-06-18 | 2020-06-18 | 一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010557741.5A CN111596549B (zh) | 2020-06-18 | 2020-06-18 | 一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111596549A CN111596549A (zh) | 2020-08-28 |
CN111596549B true CN111596549B (zh) | 2022-09-02 |
Family
ID=72186218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010557741.5A Active CN111596549B (zh) | 2020-06-18 | 2020-06-18 | 一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111596549B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112507277A (zh) * | 2020-10-28 | 2021-03-16 | 绍兴精宸智能科技有限公司 | 一种基于电信号的快速振幅预测算法 |
CN113029324B (zh) * | 2021-02-20 | 2022-06-10 | 山东骏腾医疗科技有限公司 | 一种基于超声波的快速病理组织处理方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105458427A (zh) * | 2015-12-25 | 2016-04-06 | 华南理工大学 | 一种超声振动辅助放电磨削加工装置及方法 |
CN110045603A (zh) * | 2018-12-03 | 2019-07-23 | 华南理工大学 | 一种针对压电陶瓷驱动部件变载荷环境下的鲁棒自适应控制方法 |
CN110673472A (zh) * | 2019-09-09 | 2020-01-10 | 南京理工大学 | 基于神经网络补偿死区反演误差的自适应鲁棒控制方法 |
-
2020
- 2020-06-18 CN CN202010557741.5A patent/CN111596549B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105458427A (zh) * | 2015-12-25 | 2016-04-06 | 华南理工大学 | 一种超声振动辅助放电磨削加工装置及方法 |
CN110045603A (zh) * | 2018-12-03 | 2019-07-23 | 华南理工大学 | 一种针对压电陶瓷驱动部件变载荷环境下的鲁棒自适应控制方法 |
CN110673472A (zh) * | 2019-09-09 | 2020-01-10 | 南京理工大学 | 基于神经网络补偿死区反演误差的自适应鲁棒控制方法 |
Non-Patent Citations (2)
Title |
---|
三维椭圆振动切削系统建模及其控制的研究;张毫杰;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20161215(第12期);B022-145 * |
二维椭圆振动辅助加工新方法及其装置的研究;王荣奇;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20171115(第11期);C029-22 * |
Also Published As
Publication number | Publication date |
---|---|
CN111596549A (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111596549B (zh) | 一种用于硬脆材料超声加工的鲁棒自适应振幅控制方法 | |
WO1989005002A1 (en) | Industrial control systems | |
CN109828468B (zh) | 一种针对磁滞非线性机器人系统的控制方法 | |
CN111142404A (zh) | 一种基于压电陶瓷驱动的微定位平台及其建模与控制方法 | |
CN105652662A (zh) | 一种窄带自适应滤波的压电结构振动主动控制方法 | |
CN105159069A (zh) | 一种压电陶瓷作动器的位移控制方法 | |
CN112959321B (zh) | 基于改进pi结构的机器人柔性关节转换误差补偿方法 | |
CN101795110A (zh) | 直线音圈电机非线性特性的智能补偿方法及控制系统 | |
CN116127838B (zh) | 基于负载不确定性的光学精密设备主动隔振方法 | |
CN105227035A (zh) | 一种永磁直线电机控制方法 | |
CN110045603B (zh) | 压电陶瓷驱动部件变载荷环境下的鲁棒自适应控制方法 | |
CN113114128B (zh) | 基于广义Bouc-Wen逆模型的压电前馈补偿方法 | |
CN106505909A (zh) | 基于广义回归神经网络的超声波电机非线性死区补偿控制系统及方法 | |
Jingzhuo et al. | T–S fuzzy control of travelling-wave ultrasonic motor | |
CN113110064A (zh) | 基于Duhem逆模型的压电动态迟滞前馈补偿方法 | |
Lin et al. | Tuning PID control parameters for micro-piezo-stage by using grey relational analysis | |
CN116627028B (zh) | 交联电缆生产线控制方法 | |
Tzes et al. | Adaptive precompensators for flexible-link manipulator control | |
CN107807532B (zh) | 一种用于超磁致伸缩隔振平台的自适应逆隔振控制方法 | |
CN106940524A (zh) | 一种压电定位装置的振动和非线性抑制方法 | |
CN114690642B (zh) | 一种压电驱动系统的耦合控制方法 | |
CN113110105B (zh) | 一种基于逆补偿和扰动触发的压电执行器控制方法 | |
CN106026832A (zh) | 基于改进adrc控制算法的永磁同步直线电机控制方法 | |
CN113885314B (zh) | 一种具有未知增益和干扰的非线性系统跟踪控制方法 | |
Dong et al. | Adaptive Inverse Control of Piezoelectric Actuator with non-Smooth Hysteresis Model |
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 |