CN110008505B - 一种解决斜浪中船舶水平运动响应发散的数值模型方法 - Google Patents
一种解决斜浪中船舶水平运动响应发散的数值模型方法 Download PDFInfo
- Publication number
- CN110008505B CN110008505B CN201910110673.5A CN201910110673A CN110008505B CN 110008505 B CN110008505 B CN 110008505B CN 201910110673 A CN201910110673 A CN 201910110673A CN 110008505 B CN110008505 B CN 110008505B
- Authority
- CN
- China
- Prior art keywords
- rudder
- force
- motion
- ship
- generalized
- 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
- 230000033001 locomotion Effects 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000004044 response Effects 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 239000012530 fluid Substances 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 11
- 238000005452 bending Methods 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 4
- 238000010008 shearing Methods 0.000 claims description 4
- 230000007613 environmental effect Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000005012 migration Effects 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 2
- 229920001971 elastomer Polymers 0.000 abstract 1
- 239000000806 elastomer Substances 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 230000002706 hydrostatic effect Effects 0.000 description 4
- 230000001788 irregular Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000005489 elastic deformation Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004886 head movement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- 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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种解决斜浪中船舶水平运动响应发散的数值模型方法,具体包含以下步骤,步骤一、在PID(proportional,integral and derivative)自动舵基本控制方程的基础上,建立一种以矩阵的形式将横荡和首摇运动上的舵力分量表达式步骤二、将该项力(矩)计入到刚性/弹性船体运动方程中;步骤三、将方程内各项系数、作用力分别求解;步骤四、最后基于龙格‑库塔(Runge‑Kutta)法求解方程得到船体运动响应,再借助模态叠加原理,可得到弹性体载荷响应。本发明既可以从根本上解决斜浪中船体水平时域运动方程求解带来的响应发散的缺陷,其适用于船舶和海洋平台的载荷计算中,也可为随机海浪的载荷响应预报提供基础。
Description
技术领域
本发明涉及船舶与海洋工程波浪载荷技术水动力领域,可用于刚体船型及弹性体船型的斜浪中船舶水平运动响应发散的数值模型方法。
背景技术
众所周知,真实海洋环境是斜浪,海浪可能来自任意方向。在实际航行中,波浪中的船舶都是通过不断地调整舵来保证航向。专利申请人曾在时域内对船舶水平运动进行数值模拟时发现,由于船舶水平运动不具备恢复力,所以在进行船舶水平运动的时历模拟时,一旦发生漂移,没有机制能使船舶回到原来的航向上,即使在规则波中也很难保证稳态解。水平运动的解会在时间步内发生偏移,这种偏移也将严重影响其他运动模态。
随着近年来船舶主尺度不断加大,其结构的弹性变形与流体的耦合作用效果愈加显著,本发明也同样能够应用到弹性船体结构,即水弹性方法当中。
但是现有技术方法(如人工弹簧或采用数字滤波技术)仍然存在如下问题:1.通过施加人工恢复力,即安装人工弹簧来抑制船舶横荡和首摇运动的发散。这种方法可以得到运动的稳态解,但是弹簧刚度的选取带有主观性,会影响运动的幅值和相位,因此这种方法对于实船运动预报并不可取;2.采用数字滤波技术消除发散成分。首先进行滤波所必须的最小步长数的计算,把滤波加到得到的运动位移和速度分量上,在载荷平衡中决定加速度,再把这些值作为出发值再次返回操作,采用时间积分的计算法。这种方法可以消除横荡和首摇的发散,但是该方法耗费机时过多,实际操作不方便。
发明内容
本发明的目的就是克服上述现有方法的不足,从船舶操纵性出发,借助自动舵考虑舵力效应,将其计入到运动方程中进行求解,以消除漂移现象,从而解决在时域内求解斜浪中船舶水平运动响应解发散的问题。
为了解决上述问题,本发明提供了种解决斜浪中船舶水平运动响应发散的数值模型方法,具体过程包含以下步骤:
步骤一、建立PID自动舵航向偏差修正信号模型:
;
式中δ为舵角信号,为航向偏差信号;k为比例常数;
步骤二、建立下风或上风力矩干扰情况下为使航向保持不变的PID自动舵模型:
式中,δR为偏转舵角;η2为船舶的横荡运动幅值;η6为船舶的摇首运动角度;为船舶的摇首角速度;U——船体航速;k1~k3为比例系数,根据使自动舵能适应不同的装载、航速和海况来调节;舵的有效攻角为
式中,UR——舵后有效流速;(xR,yR,zR)——舵的位置坐标;
——船舶的横荡运动速度;/>——船舶的纵荡运动速度。
步骤三、由于横荡和首摇运动自身不具备恢复力,因此本发明只考虑舵力对于横荡和首摇运动方面的影响,引入比例系数K,则横荡和首摇运动上的舵力分量以PID自动舵数值模型的形式表示为:
其中,K为比例常数,它应被整定以根据适应载重和环境要求以及避免振荡来整定。
由式(5-6),可将横荡和首摇运动上的舵力分量进一步表示为:
步骤四、最终,基于PID自动舵的舵力/矩表达式,将式(8)中横荡和首摇运动上的舵力分量以矩阵的形式表示成如下形式:
其中,t指时间。则摇首和横荡运动的阻尼系数BRudder、恢复力系数CRudder,积分项系数ERudder,以矩阵的形式表示如下:
步骤五、在时域内,船体运动的非线性水弹性力学方程可以写作如下形式:
其中,[a]、[b]、[c]分别为结构广义质量矩阵、广义阻尼矩阵、广义刚度矩阵;[A]、[B]、[C]分别为广义流体附加质量、广义流体附加阻尼和广义流体恢复力系数矩阵;
Fr I(t)——第r阶入射波力;Fr D(t)——第r阶绕射力;Fr slammin g(t)——第r阶砰击力;
pr(t)——主坐标,泛指前面提到的运动,和/>分别指运动速度和加速度;r是模态数量,r=1~6是刚体的运动模态,当r=7~m时,则为船体的弹性振动模态。
将舵力即步骤四中的于PID自动舵的舵力/矩表达式(10-12)加入到水弹性时域运动方程(13)中,最终时域方程形式如下
方程采用四阶龙格-库塔法求解,在得到主坐标后,利用模态叠加原理,得到船体结构任意剖面x处的位移w(x,t)、弯矩M(x,t)动剪切力V(x,t)。
其中,wr(x)、Mr(x)、Vr(x)分别为船体梁位移、弯矩、剪力的第r阶固有振型,可由有限元法或迁移矩阵法得到。
附图说明
图1是自动舵验证:横荡响应时历曲线;自动舵验证:横荡响应时历(航速U=14knots,航向角β=45°,λ/L=1.0;
图2是自动舵验证:艏摇响应时历曲线;艏摇响应时历(航速U=14knots,航向角β=45°,λ/L=1.0;
图3是借助自动舵得到的规则波中船舯弯矩与试验值的比较结果,U=18knots,H=7m,β=45°,λ/L=1.0;
图4是借助自动舵得到的不规则波中船舯弯矩与试验值的比较结果,不规则波中垂向弯矩结果展示U=18knots,H1/3=14m,β=45°,TZ=9.88s。
具体实施方式
下面具体说明本发明的实施方式:
以弹性船体为例,一般而言船舶运动可看做是一种波浪激励力下的强迫振动,所以有,
其中,[a],[b],[c]——结构广义质量矩阵、广义阻尼矩阵、广义刚度矩阵.
{Fr(t)}——每一瞬时第r阶总的流体力,可看做是作用在船体上一系列流体力的组合,可表示为:
r=1,2...6——表示刚体的六个自由度的运动模态;r≥7——弹性模态。
{Fr(t)}={Fr S(t)}+{Fr I(t)}+{Fr D(t)}+{Fr R(t)}+{Fr slamming(t)}+{Fr δ(t)} (12)
其中,{Fr s(t)}——第r阶静水回复力;{Fr δ(t)}——第r阶舵力(矩)。
Fr I(t)、Fr D(t)——第r阶船体所受的入射波力、绕射波力;
{Fr slamming(t)}——第r阶船体所受的砰击力;{Fr R(t)}——第r阶辐射力;
·静水恢复力
根据三维水弹性理论,船体广义水动力在第r阶的分量可表示为:
其中,P——压力;——法向向量;/>——弹性结构第r阶位移矢量。
船体所受静水压力为:
其中,ρ——水密度;g——重力加速度;
pka(t)——第k阶主坐标幅值;wk——第k垂向位移。
借助三维船体瞬时网格截取,第r阶船体瞬时平均湿表面上的非线性静水恢复力载荷可表示为瞬时湿表面S(t)上的流体静压力的积分:.
其中,[Crk]——广义流体恢复力系数矩阵。
·入射力和绕射力
同理,对于入射波力和绕射波力的计算,同样借助船体瞬时网格的截取,每时每刻将入射波、绕射波压力在瞬时湿表面网格上进行积分:
·辐射力
第r阶辐射力表示为:
式中,广义流体附加质量、广义流体附加阻尼可表示为
Re和Im分别表示取结果的实部和虚部。
·砰击力
基于动量冲击理论,砰击载荷可表示为:
其中,m(x,t)——瞬时附加质量;s(x,t)——瞬时浸湿面积;
ZR(x,t)——船体与波浪的垂向相对位移。船体与波浪的垂向相对速度可表示为
其中,——入射波面垂向速度。
然后将砰击载荷在瞬时湿表面S(t)进行积分,则第r阶砰击力为:
·舵力
PID自动舵的舵力/矩表达式,可以将横荡和首摇运动上的舵力分量以矩阵的形式表达出来,即有:
对于摇首和横荡运动的阻尼系数BRudder及恢复力系数CRudder,积分项系数ERudder,以矩阵的形式表示如下:
式中:K为比例常数,它应被整定以根据适应载重和环境要求以及避免振荡来整定;
UR——计入船体、螺旋桨影响的舵后有效流速。舵后水流的速度可采用下式确定:
UR=1.2U (26)
最终,斜浪中有航速船体运动和波浪载荷的水弹性非线性时域运动方程形式如下:
为解决这一问题,采用龙格-库塔(Runge-Kutta)法。该法为显式单步法,具有4阶精度,其步进求解的格式为:
在得到主坐标后,利用模态叠加原理,得到时域内船体上任一横剖面的位移w(x,t)、弯矩M(x,t)、剪切力V(x,t)。
其中,wr(x)、Mr(x)、Vr(x)分别为船体梁位移、弯矩、剪力的第r阶固有振型。r=1~6是刚体的运动模态,当r≥7时,则为船体的弹性振动模态。
本发明的技术方案过程简单,易于实现,能够计及船体弹性变形的影响,已有试验证明,相比现有方法结果更为精确且计算效率明显提高,可用于船体和海洋工程结构的载荷响应预报,因此具有较大的应用前景。
附图1~2是对某船自动舵解决运动响应发散的验证。
为了验证本发明方法的优越性,采用本发明的方法现对试验中的规则波和不规则波中目标船的船舯弯矩进行计算,并与试验结果进行比较,参见附图3-4。
通过以上比较可以发现,对于不规则波,由于主峰频率法中对于各项作用力、水动力系数等都是基于主峰频率进行计算的,无法体现出不规则波响应特性。而本发明得到的载荷响应却能够与试验结果较好的吻合,同时也证明了本发明能够满足工程上的要求。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明的保护范围,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种解决斜浪中船舶水平运动响应发散的数值模型方法,包含以下步骤,
步骤一、建立PID自动舵航向偏差修正信号模型:
;
式中δ为舵角信号,为航向偏差信号;k为比例常数;
步骤二、建立下风或上风力矩干扰情况下为使航向保持不变的PID自动舵模型:
式中,δR为偏转舵角;η2为船舶的横荡运动幅值;η6为船舶的首摇运动角度;为船舶的首摇角速度;U——船体航速;k1~k3为比例系数,根据使自动舵能适应不同的装载、航速和海况来调节;舵的有效攻角为
式中,UR——舵后有效流速;(xR,yR,zR)——舵的位置坐标;
——船舶的横荡运动速度;/>——船舶的纵荡运动速度;
步骤三、由于横荡和首摇运动自身不具备恢复力,只考虑舵力对于横荡和首摇运动方面的影响,引入比例系数K,则横荡和首摇运动上的舵力分量以PID自动舵数值模型的形式表示为:
其中,K为比例常数,它应被整定以根据适应载重和环境要求以及避免振荡来整定;
将横荡和首摇运动上的舵力分量进一步表示为:
步骤四、最终,基于PID自动舵的舵力/矩表达式,将式(8)中横荡和首摇运动上的舵力分量以矩阵的形式表示成如下形式:
其中,t指时间;则首摇和横荡运动的阻尼系数BRudder、恢复力系数CRudder,积分项系数ERudder,以矩阵的形式表示如下:
步骤五、在时域内,船体运动的非线性水弹性力学方程写作如下形式:
其中,[a]、[b]、[c]分别为结构广义质量矩阵、广义阻尼矩阵、广义刚度矩阵;[A]、[B]、[C]分别为广义流体附加质量、广义流体附加阻尼和广义流体恢复力系数矩阵;Fr I(t)——第r阶入射波力;Fr D(t)——第r阶绕射力;Fr slamming(t)——第r阶砰击力;pr(t)——主坐标,泛指前面提到的运动,和/>分别指运动速度和加速度;r是模态数量,r=1~6是刚体的运动模态,当r=7~m时,则为船体的弹性振动模态;
将舵力即步骤四中的于PID自动舵的舵力/矩表达式加入到水弹性时域运动方程(13)中,最终时域方程形式如下
其中,[a]、[b]、[c]分别为结构广义质量矩阵、广义阻尼矩阵、广义刚度矩阵;[A]、[B]、[C]分别为广义流体附加质量、广义流体附加阻尼和广义流体恢复力系数矩阵;
Fr I(t)——第r阶入射波力;Fr D(t)——第r阶绕射力;Fr slamming(t)——第r阶砰击力;pr(t)——主坐标,泛指前面提到的运动,和/>分别指运动速度和加速度;r是模态数量,r=1~6是刚体的运动模态,当r=7~m时,则为船体的弹性振动模态;
方程采用四阶龙格-库塔法求解,在得到主坐标后,利用模态叠加原理,得到船体结构任意剖面x处的位移w(x,t)、弯矩M(x,t)动剪切力V(x,t);
其中,wr(x)、Mr(x)、Vr(x)分别为船体梁位移、弯矩、剪力的第r阶固有振型,由有限元法或迁移矩阵法得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910110673.5A CN110008505B (zh) | 2019-02-12 | 2019-02-12 | 一种解决斜浪中船舶水平运动响应发散的数值模型方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910110673.5A CN110008505B (zh) | 2019-02-12 | 2019-02-12 | 一种解决斜浪中船舶水平运动响应发散的数值模型方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110008505A CN110008505A (zh) | 2019-07-12 |
CN110008505B true CN110008505B (zh) | 2023-07-25 |
Family
ID=67165726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910110673.5A Active CN110008505B (zh) | 2019-02-12 | 2019-02-12 | 一种解决斜浪中船舶水平运动响应发散的数值模型方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110008505B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112115652B (zh) * | 2020-09-23 | 2022-04-05 | 中国海洋大学 | 一种层合型浮体结构物的水动力特性分析方法 |
CN113704965B (zh) * | 2021-07-22 | 2024-06-11 | 中国海洋大学 | 浮托安装荷载转移过程动态模拟分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246515A (zh) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | 基于cfd软件的数字船模平面运动机构实验方法 |
US7756689B1 (en) * | 2006-02-28 | 2010-07-13 | The United States Of America As Represented By The Secretary Of The Navy | Numerical modeling of six-degree-freedom ship motion |
-
2019
- 2019-02-12 CN CN201910110673.5A patent/CN110008505B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7756689B1 (en) * | 2006-02-28 | 2010-07-13 | The United States Of America As Represented By The Secretary Of The Navy | Numerical modeling of six-degree-freedom ship motion |
CN101246515A (zh) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | 基于cfd软件的数字船模平面运动机构实验方法 |
Non-Patent Citations (1)
Title |
---|
不同浪向下超大型船舶载荷响应特征的模型试验研究;陈占阳;李志鹏;;振动与冲击;第36卷(第19期);112-118 * |
Also Published As
Publication number | Publication date |
---|---|
CN110008505A (zh) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Perez | Ship motion control: course keeping and roll stabilisation using rudder and fins | |
CN104527943B (zh) | 一种基于能量优化的综合减摇装置双重神经网络自整定pid控制方法 | |
CN109917657B (zh) | 动力定位船舶的抗干扰控制方法、装置与电子设备 | |
CN110008505B (zh) | 一种解决斜浪中船舶水平运动响应发散的数值模型方法 | |
CN116911135B (zh) | 计及波浪下砰击载荷的非线性水弹性时域计算评估方法 | |
CN109739248B (zh) | 基于adrc的船载三自由度并联稳定平台稳定控制方法 | |
CN110794843B (zh) | 基于观测器的非线性船舶时滞动力定位船鲁棒镇定系统 | |
CN111552182B (zh) | 一种基于干扰观测器的船舶动力定位系统的抗干扰饱和控制方法 | |
Zhu et al. | Numerical modeling of a spar platform tethered by a mooring cable | |
Awad et al. | Ship roll damping via direct inverse neural network control system | |
CN107037821A (zh) | 连续冲击载荷下的水下潜器运动姿态估计及控制方法 | |
Zhao et al. | Simulation analysis of rudder roll stabilization during ship turning motion | |
Song et al. | Development of a control algorithm for active control of rolling motion of a ship using a gyrostabilizer | |
Munif et al. | Modeling extreme roll motions and capsizing of a moderate-speed ship in astern waves | |
Wang et al. | Motion response of a moored semi-submersible-type single module of a vlfs in multi-slope shallow water | |
Ghassemi et al. | Neural network-PID controller for roll fin stabilizer | |
Patil et al. | A time-domain method for analyzing the ship roll stabilization based on active fin control | |
Guan et al. | Concise robust fin roll stabilizer design based on integrator backstepping and CGSA | |
Auestad et al. | Motion compensation system for a free floating surface effect ship | |
Demirel et al. | LMI-BASED H2 AND H STATE-FEEDBACK CONTROLLER DESIGN FOR FIN STABILIZER OF NONLINEAR ROLL MOTION OF A FISHING BOAT | |
Wang et al. | Nonsingular fast terminal sliding mode-based robust adaptive structural reliable position-mooring control with uncertainty estimation | |
JP5296737B2 (ja) | 入射波の波高及び波向き推定方法、自動位置保持制御方法、自動位置保持システム及び船舶と洋上構造物 | |
Uyar et al. | Dynamic modelling, investigation of manoeuvring capability and navigation control of a cargo ship by using matlab simulation | |
Zhang et al. | Multi-dimensional Taylor network optimal control for ship roll stabilization | |
Sclavounos | Intersections between marine hydrodynamics and optimal control theory |
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 |