CN106874547B - 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 - Google Patents
一种预测纤维全缠绕复合材料气瓶爆破压力的方法 Download PDFInfo
- Publication number
- CN106874547B CN106874547B CN201710012011.5A CN201710012011A CN106874547B CN 106874547 B CN106874547 B CN 106874547B CN 201710012011 A CN201710012011 A CN 201710012011A CN 106874547 B CN106874547 B CN 106874547B
- Authority
- CN
- China
- Prior art keywords
- fiber
- stress
- gas cylinder
- matrix
- damage
- 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
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]
Abstract
本发明涉及高压储氢气瓶制造技术,旨在提供一种预测纤维全缠绕复合材料气瓶爆破压力的方法。包括下述过程:建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型;基于微观力学对复合材料建立多尺度模型;利用FORTRAN语言编写的ABAQUS‑UMAT即ABAQUS用户静态材料子程序模块,基于多尺度模型对复合材料气瓶进行渐进失效分析,求取爆破压力;本发明从微观角度对纤维和基体组分的损伤进行判别以及对损伤变量进行计算,相比于传统的唯象宏观失效准则,该方法能更准确的确定复合材料的失效模式,需要确定的材料参数则更少,因此能高效准确的对纤维全缠绕复合材料气瓶进行渐进失效分析,获取爆破压力。
Description
技术领域
本发明是关于高压储氢气瓶制造技术领域,特别涉及一种预测纤维全缠绕复合材料气瓶爆破压力的方法。
背景技术
当前,纤维全缠绕复合材料气瓶正广泛应用于航空航天、压力容器、新能源汽车等高新技术领域。复合材料气瓶爆破压力的确定是气瓶最为重要的方面,直接与复合材料气瓶的优化设计密切相关,但是现有复合材料气瓶的爆破压力主要是通过实验方式确定,这种方法代价昂贵且效率低下。
随着有限元技术的发展,采用有限元方法预测复合材料气瓶爆破压力逐渐兴起,但是复合材料气瓶的有着复杂的失效模式,比如内胆塑性变形、纤维断裂、基体开裂以及纤维/基体界面分离,因此需要可靠的失效理论和损伤演化方法预测纤维全缠绕复合材料气瓶的爆破压力。
复合材料气瓶的复合材料缠绕层具有明显的各向异性,目前预测纤维全缠绕复合材料气瓶的爆破压力主要是基于Tsai-Wu、Hashin、Chang-Chang等唯象宏观的失效准则和相应的损伤演化方式,这些失效准则需要经验确定气瓶的失效模式和相应的材料参数,因此难以准确可靠的预测复合材料气瓶的爆破压力。相反,从组分角度出发的微观力学失效理论可以准确的判别失效模式,需要确定的材料参数更少,因此提出一种基于宏-微观多尺度模型预测纤维全缠绕复合材料气瓶爆破压力的方法十分重要。
复合材料及复合材料气瓶的计算分析一般可采用通用的有限元软件ABAQUS进行,但是通过ABAQUS软件自带的模块仅可以基于二维Hashin唯象宏观失效准则对复合材料气瓶渐进失效分析,但无法直接采用宏-微观多尺度模型预测纤维全缠绕复合材料气瓶的爆破压力。
发明内容
本发明的主要目的在于克服现有技术中的不足,利用ABAQUS用户子程序界面提供一种基于宏-微观多尺度模型预测纤维全缠绕复合材料气瓶爆破压力的方法。该方法从纤维和基体组分角度判别纤维、基体的损伤以及更新应力应变,对气瓶进行渐进失效的分析,求取爆破压力。
为解决上述技术问题,本发明的解决方案是:提供一种预测纤维全缠绕复合材料气瓶爆破压力的方法,包括下述过程:
一、建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型;
二、基于微观力学对复合材料建立多尺度模型;
三、利用(使用FORTRAN语言编写)ABAQUS-UMAT即ABAQUS用户静态材料子程序模块,基于多尺度模型对复合材料气瓶进行渐进失效分析,求取爆破压力;
所述过程一中,建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型:基于ABAQUS建立复合材料气瓶部件,分别对内胆和复合材料纤维层设置材料属性和划分网格,再用ABAQUS-ASSEMBLY模块对其组装之后设置分析步、变量输出、加载方式以及边界条件;
所述过程二具体包括下述步骤:
步骤(1):建立微观应力和宏观应力的对应关系以及求取应力放大系数分量:
(a)微观应力和宏观应力的对应关系
其中,Mσ为应力放大系数,Mij(i,j=1,2,3,4,5,6)为应力放大系数的分量,σ为微观应力,为宏观应力。
(b)获取应力放大系数分量;
为了获得应力放大系数,本发明建立正六面体的复合材料微观胞元,其尺寸为针对三个正方向(1,2,3)和三个剪切方向(4,5,6)六种加载情况,分别对胞元施加单位应力载荷,得到微观胞元的应力分布。其中微观胞元每点的应力分布即为该点的应力放大系数。为了简化,分别从微观胞元的纤维和基体上选取若干点,计算出纤维和基体的应力放大系数并且存储在参数文件中。
步骤(2):建立宏观应力和应变的本构关系。
其中,和为宏观应力和宏观应变。Cd为含损伤刚度矩阵,df和dm为纤维损伤变量和基体损伤变量。为弹性刚度矩阵分量,由材料弹性模量和泊松比等材料参数计算得到。
步骤(3):建立基于微观力学的纤维和基体失效判据:
(c)针对纤维组分,损伤初始判据为:
-Cf<σf<Tf;
其中,Tf和Cf为纤维拉伸强度和纤维压缩强度,σf为纵向方向的纤维微观应力。
(d)针对基体组分,损伤初始判据为:
其中,Tm和Cm分别为基体拉伸和压缩强度,σi(i=1,2,3,4,5,6)为六个微观应力分量,I1和σvm分别为第一微观应力不变量和微观Mises应力。
步骤(4):采用参数化折减刚度矩阵的方法,对满足失效准则的单元应力进行应力更新。
当纤维达到初始失效判据时,将df设为1,当基体达到初始失效判据时,将dm设为1。
所述过程三具体包括下述步骤:
步骤(5):通过ABAQUS用户子程序UMAT实现上述过程,ABAQUS-UMAT是由应变驱动,把纤维损伤和基体损伤变量定为UMAT的状态变量。随着载荷增加,首先读取材料信息、应力放大系数和状态变量值SDV,根据步骤(2)计算宏观应力,再由步骤(1)计算微观应力,将计算的微观应力代入步骤(3)判别纤维和基体是否进入损伤,如果没有进入损伤,则直接进行应力应变的更新;如果进入损伤,则根据步骤(4)对单元刚度矩阵进行折减,使用折减的单元刚度矩阵进行应力应变的更新,同时将损伤数值存储为状态变量SDV;
步骤(6):将过程一建立的模型主文件和步骤(5)建立的ABAQUS-UMAT用户子程序联合,使用ABAQUS/Standard模块对纤维全缠绕复合材料气瓶进行渐进失效分析,得到爆破压力。
与现有技术相比,本发明的有益效果是:
本发明利用ABAQUS用户子程序界面来数值实现纤维全缠绕复合材料气瓶爆破压力的预测方法,该方法从微观角度对纤维和基体组分的损伤进行判别以及对损伤变量进行计算,相比于传统的唯象宏观失效准则,该方法能更准确的确定复合材料的失效模式,需要确定的材料参数则更少,因此能高效准确的对纤维全缠绕复合材料气瓶进行渐进失效分析,获取爆破压力。
附图说明
图1为本发明实施例纤维全缠绕复合材料气瓶有限元模型图。
图2为实施例中建立的微观胞元及用于应力放大系数计算的若干点分布图。
图3为本发明对所提出的基于宏-微观多尺度模型的UMAT数值实现流程图。
图4为实施例中纤维全缠绕复合材料气瓶的内压-位移示意图。
具体实施方式
首先需要说明的是,本发明是计算机技术在高压储氢气瓶领域的一种应用。在本发明的实现过程中,会涉及到多个软件功能模块的应用。申请人认为,如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后,在结合现有公知技术的情况下,本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。凡本发明申请文件提及的均属此范畴,申请人不再一一列举。
一种预测纤维全缠绕复合材料气瓶爆破压力的方法,其特征在于,包括下述过程:
一、建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型;
二、基于微观力学对复合材料建立多尺度模型;
三、利用FORTRAN语言编写的ABAQUS-UMAT即用户静态材料子程序模块,基于多尺度模型对复合材料气瓶进行渐进失效分析,求取爆破压力;
所述过程一具体包括下述步骤:基于ABAQUS建立复合材料气瓶部件,分别对内胆和复合材料纤维层设置材料属性和划分网格,再用ABAQUS-ASSEMBLY模块对其组装之后设置分析步、变量输出、加载方式以及边界条件;
所述过程二具体包括下述步骤:
步骤(1):建立微观应力和宏观应力的对应关系以及求取应力放大系数分量:
(a)微观应力和宏观应力的对应关系
其中,Mσ为应力放大系数,Mij(i,j=1,2,3,4,5,6)为应力放大系数的分量,σ为微观应力,为宏观应力;
(b)获取应力放大系数分量;
为获得应力放大系数,建立正六面体的复合材料微观胞元,其尺寸为针对三个正方向(1,2,3)和三个剪切方向(4,5,6)六种加载情况,分别对胞元施加单位应力载荷,得到微观胞元的应力分布;其中微观胞元每点的应力分布即为该点的应力放大系数;为了简化,分别从微观胞元的纤维和基体上选取若干点,计算出纤维和基体的应力放大系数并且存储在参数文件中;
步骤(2):建立宏观应力和应变的本构关系:
其中,和为宏观应力和宏观应变;Cd为含损伤刚度矩阵,df和dm为纤维损伤变量和基体损伤变量;为弹性刚度矩阵分量,由材料弹性模量和泊松比等材料参数计算得到;
步骤(3):建立基于微观力学的纤维和基体失效判据:
(c)针对纤维组分,损伤初始判据为:
-Cf<σf<Tf;
其中,Tf和Cf为纤维拉伸强度和纤维压缩强度,σf为纵向方向的纤维微观应力;
(d)针对基体组分,损伤初始判据为:
其中,Tm和Cm分别为基体拉伸和压缩强度,σi(i=1,2,3,4,5,6)为六个微观应力分量,I1和σvm分别为第一微观应力不变量和微观Mises应力;
步骤(4):采用参数化折减刚度矩阵的方法,对满足失效准则的单元应力进行应力更新:
当纤维达到初始失效判据时,将df设为1,当基体达到初始失效判据时,将dm设为1;
所述过程三具体包括下述步骤:
步骤(5):通过ABAQUS用户子程序UMAT实现上述过程,ABAQUS-UMAT是由应变驱动,把纤维损伤和基体损伤变量定为UMAT的状态变量;随着载荷增加,首先读取材料信息、应力放大系数和状态变量值SDV,根据步骤(2)计算宏观应力,再由步骤(1)计算微观应力,将计算的微观应力代入步骤(3)判别纤维和基体是否进入损伤,如果没有进入损伤,则直接进行应力应变的更新;如果进入损伤,则根据步骤(4)对单元刚度矩阵进行折减,使用折减的单元刚度矩阵进行应力应变的更新,同时将损伤数值存储为状态变量SDV;
步骤(6):将过程一建立的模型主文件和步骤(5)建立的ABAQUS-UMAT用户子程序联合,使用ABAQUS/Standard模块对纤维全缠绕复合材料气瓶进行渐进失效分析,得到爆破压力。
以下用基于上述过程应用的具体实施例子来说明本发明的实现效果:
在ABAQUS/CAE中建立纤维全缠绕复合材料气瓶的有限元模型,如图1所示。其中复合材料层的铺层方式为:[90°/±14°/90°/±14°/90°/±14°/90°],环向层的复合材料厚度为0.32mm,螺旋层的复合材料厚度为0.18mm,总共分为10个单层,均用减缩积分三维实体单元C3D8R来仿真,为了充分利用对称性,复合材料气瓶建立四分之一模型并且施加合适的边界条件。
利用ABAQUS/Standard计算复合材料气瓶的渐进失效,如图2采用微观胞元计算各点的应力放大系数(纤维用F表示,基体用M表示)并存储到参数文件当中,如图3利用用户子程序UMAT首先读取材料信息、应力放大系数和状态变量值SDV计算宏观应力和微观应力,再判别纤维和基体损伤状态和计算损伤变量,进行应力应变更新。将模型主文件和ABAQUS-UMAT用户子程序联合,使用ABAQUS/Standard模块对纤维全缠绕复合材料气瓶进行渐进失效分析,得到爆破压力。
图4分别为纤维全缠绕复合材料气瓶的内压-位移示意图,数值计算的爆破压力为71MPa,而实验爆破压力为69MPa,可见与试验值比较准确的吻合。所以本发明提出的宏-微观多尺度模型可以较为准确的预测纤维全缠绕复合材料气瓶的爆破压力。
本发明在ABAQUS软件的基础上进行用户子程序的开发,提出的宏-微观多尺度模型相比于传统的唯象宏观失效准则,该方法能更准确的确定复合材料的失效模式,需要确定的材料参数更少,能高效且较为准确预测复合材料气瓶的爆破压力,为提升轻量化强度设计水平提供了技术支撑。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (1)
1.一种预测铝内胆纤维全缠绕复合材料气瓶爆破压力的方法,其特征在于,包括下述过程:
一、建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型;
二、基于微观力学对复合材料建立多尺度模型;
三、利用FORTRAN语言编写的ABAQUS-UMAT即用户静态材料子程序模块,基于多尺度模型对复合材料气瓶进行渐进失效分析,求取爆破压力;
所述过程一中,建立含变厚度变角度封头的纤维全缠绕复合材料气瓶有限元模型:基于ABAQUS建立复合材料气瓶部件,分别对内胆和复合材料纤维层设置材料属性和划分网格,再用ABAQUS-ASSEMBLY模块对其组装之后设置分析步、变量输出、加载方式以及边界条件;
所述过程二具体包括下述步骤:
步骤(1):建立微观应力和宏观应力的对应关系以及求取应力放大系数分量:
(a)微观应力和宏观应力的对应关系
其中,Mσ为应力放大系数,Mij(i,j=1,2,3,4,5,6)为应力放大系数的分量,σ为微观应力,为宏观应力;
(b)获取应力放大系数分量;
为了获得应力放大系数,本发明建立正六面体的复合材料微观胞元,其尺寸为针对三个正方向(1,2,3)和三个剪切方向(4,5,6)六种加载情况,分别对胞元施加单位应力载荷,得到微观胞元的应力分布;其中微观胞元每点的应力分布即为该点的应力放大系数;为了简化,分别从微观胞元的纤维和基体上选取若干点,计算出纤维和基体的应力放大系数并且存储在参数文件中;
步骤(2):建立宏观应力和应变的本构关系;
其中,和为宏观应力和宏观应变;Cd为含损伤刚度矩阵,df和dm为纤维损伤变量和基体损伤变量;为弹性刚度矩阵分量,由材料弹性模量和泊松比材料参数计算得到;
步骤(3):建立基于微观力学的纤维和基体失效判据:
(c)针对纤维组分,损伤初始判据为:
-Cf<σf<Tf;
其中,Tf和Cf为纤维拉伸强度和纤维压缩强度,σf为纵向方向的纤维微观应力;
(d)针对基体组分,损伤初始判据为:
其中,Tm和Cm分别为基体拉伸和压缩强度,σi(i=1,2,3,4,5,6)为六个微观应力分量,I1和σvm分别为第一微观应力不变量和微观Mises应力;
步骤(4):采用参数化折减刚度矩阵的方法,对满足失效准则的单元应力进行应力更新;
当纤维达到初始失效判据时,将df设为1,当基体达到初始失效判据时,将dm设为1;
所述过程三具体包括下述步骤:
步骤(5):通过ABAQUS用户子程序UMAT实现上述过程,ABAQUS-UMAT是由应变驱动,把纤维损伤和基体损伤变量定为UMAT的状态变量;随着载荷增加,首先读取材料信息、应力放大系数和状态变量值SDV,根据步骤(2)计算宏观应力,再由步骤(1)计算微观应力,将计算的微观应力代入步骤(3)判别纤维和基体是否进入损伤,如果没有进入损伤,则直接进行应力应变的更新;如果进入损伤,则根据步骤(4)对单元刚度矩阵进行折减,使用折减的单元刚度矩阵进行应力应变的更新,同时将损伤数值存储为状态变量SDV;
步骤(6):将过程一建立的模型主文件和步骤(5)建立的ABAQUS-UMAT用户子程序联合,使用ABAQUS/Standard模块对纤维全缠绕复合材料气瓶进行渐进失效分析,得到爆破压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710012011.5A CN106874547B (zh) | 2017-01-08 | 2017-01-08 | 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710012011.5A CN106874547B (zh) | 2017-01-08 | 2017-01-08 | 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106874547A CN106874547A (zh) | 2017-06-20 |
CN106874547B true CN106874547B (zh) | 2019-07-16 |
Family
ID=59164503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710012011.5A Active CN106874547B (zh) | 2017-01-08 | 2017-01-08 | 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106874547B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742001B (zh) * | 2017-09-05 | 2021-11-02 | 兰州空间技术物理研究所 | 一种超高压多层缠绕复合材料气瓶的ansys仿真方法 |
CN110686816B (zh) * | 2019-10-30 | 2021-06-15 | 江苏理工学院 | 一种高压复合材料气瓶应力过载的预警检测方法 |
CN111209693B (zh) * | 2019-12-10 | 2023-09-15 | 合肥工业大学 | 一种气瓶受外物撞击后的爆破强度评估方法 |
CN112613220B (zh) * | 2021-01-05 | 2022-09-30 | 大连理工大学 | 一种关于金属内衬纤维缠绕气瓶的承载能力的预测方法 |
CN113158489B (zh) * | 2021-05-13 | 2022-08-26 | 中国石油大学(华东) | 一种基于等效载荷的抗爆管道壁厚的计算方法 |
CN113484102B (zh) * | 2021-07-12 | 2022-12-06 | 吉林化工学院 | 一种碳纤维复合材料的石油炼化气体取样器及取样方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101586951A (zh) * | 2009-05-31 | 2009-11-25 | 哈尔滨工业大学 | 一种纤维缠绕复合材料压力容器封头厚度预测方法 |
CN104133930A (zh) * | 2014-04-27 | 2014-11-05 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料层合板损伤过程模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2805097A4 (en) * | 2012-01-20 | 2015-12-02 | Lightsail Energy Inc | COMPRESSED GAS STORAGE UNIT |
-
2017
- 2017-01-08 CN CN201710012011.5A patent/CN106874547B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101586951A (zh) * | 2009-05-31 | 2009-11-25 | 哈尔滨工业大学 | 一种纤维缠绕复合材料压力容器封头厚度预测方法 |
CN104133930A (zh) * | 2014-04-27 | 2014-11-05 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料层合板损伤过程模拟方法 |
Non-Patent Citations (3)
Title |
---|
《UHMWPE纤维缠绕铝内衬复合材料气瓶爆破压力预测》;徐明林等;《纤维复合材料》;20120613;第19-22页 |
《低温复合材料气瓶爆破压力的有限元分析与试验研究》;王洪锐等;《压力容器》;20160601;第39-44页 |
《复合材料气瓶的多轴疲劳寿命预测研究》;黄其忠等;《玻璃钢/复合材料》;20161130;第39-45页 |
Also Published As
Publication number | Publication date |
---|---|
CN106874547A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106874547B (zh) | 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 | |
Rafiee et al. | Stochastic prediction of burst pressure in composite pressure vessels | |
Soto et al. | Low velocity impact and compression after impact simulation of thin ply laminates | |
Yang et al. | Structural investigation of composite wind turbine blade considering structural collapse in full-scale static tests | |
Liu et al. | Numerical simulation and optimal design for composite high-pressure hydrogen storage vessel: A review | |
Dahmani et al. | Crack identification in reinforced concrete beams using ANSYS software | |
CN108804735B (zh) | 低速冲击下复合材料层合板渐进失效的多尺度预测方法 | |
CN112362520B (zh) | 一种有限元-离散元耦合数值模拟程序(fdem)输入参数快速标定方法 | |
EP3051447A1 (en) | Strength prediction system and method for composite laminates | |
CN105740566B (zh) | 一种预测层状复合材料层内损伤和层间分层的有限元方法 | |
CN106909708B (zh) | 确定铝内胆纤维全缠绕复合材料气瓶最佳自紧压力的方法 | |
CN111368389A (zh) | 一种预测复合材料层合板失效强度的方法 | |
Mishnaevsky Jr et al. | Statistical modelling of compression and fatigue damage of unidirectional fiber reinforced composites | |
Liao et al. | Finite element analysis of dynamic responses of composite pressure vessels under low velocity impact by using a three-dimensional laminated media model | |
CN111209693A (zh) | 一种气瓶受外物撞击后的爆破强度评估方法 | |
Lam et al. | Study of stress intensity factor of a cracked steel plate with a single-side CFRP composite patching | |
Kroupa et al. | Improved nonlinear stress-strain relation for carbon-epoxy composites and identification of material parameters | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
Xie et al. | GWFMM model for bi-modulus orthotropic materials: application to mechanical analysis of 4D-C/C composites | |
Thionnet et al. | Fibre break failure processes in unidirectional composites. Part 1: Failure and critical damage state induced by increasing tensile loading | |
Orifici et al. | An analysis tool for design and certification of postbuckling composite aerospace structures | |
Foster et al. | Tensile strength of titanium matrix composites: direct numerical simulations and analytic models | |
CN112906263B (zh) | 含制孔分层损伤的复合材料层合板强度预测方法 | |
Borrelli et al. | A repair criterion for impacted composite structures based on the prediction of the residual compressive strength | |
Ratcliffe et al. | A model for simulating the response of aluminum honeycomb structure to transverse loading |
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 |