CN110852015B - 一种陶瓷基复合材料模态的非线性计算方法 - Google Patents
一种陶瓷基复合材料模态的非线性计算方法 Download PDFInfo
- Publication number
- CN110852015B CN110852015B CN201911138360.7A CN201911138360A CN110852015B CN 110852015 B CN110852015 B CN 110852015B CN 201911138360 A CN201911138360 A CN 201911138360A CN 110852015 B CN110852015 B CN 110852015B
- Authority
- CN
- China
- Prior art keywords
- matrix
- damage
- matrix composite
- ceramic matrix
- strain
- 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
- 239000011153 ceramic matrix composite Substances 0.000 title claims abstract description 46
- 239000000463 material Substances 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000013016 damping Methods 0.000 claims abstract description 28
- 230000004044 response Effects 0.000 claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims abstract description 16
- 230000010354 integration Effects 0.000 claims abstract description 7
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 51
- 230000008859 change Effects 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 230000001052 transient effect Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000003938 response to stress Effects 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000000844 transformation Methods 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 239000011825 aerospace material Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种陶瓷基复合材料模态的非线性计算方法,包括:1、建立三维有限元模型,获得单元、结点信息;2、基于损伤建立陶瓷基复合材料刚度模型;3、求解广义特征值问题,获得固有频率和固有振型;4、基于损伤建立陶瓷基复合材料阻尼模型;5、基于模态叠加法计算陶瓷基复合材料在振动载荷下的位移、应力、应变响应;6、基于陶瓷基复合材料的非线性应力‑应变曲线,判断材料的损伤情况,若损伤饱和,输出模态和应力、应变,否则重复执行步骤2‑6。本发明运用变刚度模态叠加法,并结合有限元,快速给出振动载荷下陶瓷基复合材料的模态响应。相比于现有计算方法,考虑了陶瓷基复合材料的非线性本构,取得比直接积分法高的计算效率,节省了大量时间。
Description
技术领域
本发明属于复合材料力学分析技术领域,具体涉及一种陶瓷基复合材料模态的非线性计算方法。
背景技术
陶瓷基复合材料高温下具有良好的性能,同时具有重量轻、模量高、抗拉强度高、不易疲劳破坏等优点。陶瓷基复合材料在热端部件的应用潜力巨大,是未来航空航天材料的必然选择。在发动机实际工况下,陶瓷基复合材料往往承受复杂振动载荷的作用,振动载荷产生损伤导致刚度衰退,并最终引起结构的破坏。因此,研究陶瓷基复合材料的动力学特性具有十分重要的工程意义。而模态分析是动力特性研究的一种方法,可为结构的振动特性分析、振动故障诊断以及动力特性的优化设计提供依据。
由于陶瓷基复合材料是一种新型非线性材料,国内外还没有高效的方法预测振动载荷下的模态响应,也未见公开的发明专利。Birman(Birman V,Byrd L W.Damping inceramic matrix composites with matrix cracks[J].International Journal ofSolids&Structures,2003,40(16):4239-4256.)将线弹性振动理论与单向CMCs细观损伤力学模型结合,但没有深入分析由变刚度引起的非线性振动现象。直接积分法直接进行逐步数值积分,以中心差分法的显式算法和Newmark法的隐式算法为代表,在非线性系统中得到广泛的应用。Gao(Gao X,Han D,Chen J,et al.Numerical and experimental study onthe nonlinear dynamic response of a ceramic matrix composites beam[J].Ceramics International,2018,44(6).)将中心差分法与有限元相结合,研究了单向陶瓷基复合材料梁的非线性动力响应,但计算非常耗时,效率不高。当前,如何快速预测陶瓷基复合材料的非线性模态响应是本技术领域重要的技术问题。
发明内容
本发明针对现有技术中的不足,提供一种可以快速获得陶瓷基复合材料振动载荷下非线性模态响应的方法。
为实现上述目的,本发明采用以下技术方案:
一种陶瓷基复合材料模态的非线性计算方法,其特征在于,包括如下步骤:
步骤1、建立陶瓷基复合材料的三维有限元模型,获得单元、结点信息;
步骤2、运用有限元法建立基于损伤的陶瓷基复合材料刚度模型;
步骤3、基于步骤1和步骤2求解广义特征值问题,获得固有频率和固有振型;
步骤4、由步骤2建立的刚度模型对Rayleigh阻尼进行同步更新变化,最终建立基于损伤的陶瓷基复合材料阻尼模型;
步骤5、将步骤3获得的模态结果代入步骤2建立的刚度模型和步骤4建立的阻尼模型,运用振型叠加法计算陶瓷基复合材料在振动载荷下的位移、应力、应变响应;
步骤6、基于陶瓷基复合材料的非线性应力应变曲线,判断材料的损伤情况,若损伤饱和,输出模态和应力、应变,否则更新弹性参数,重复执行步骤2至步骤6。
为优化上述技术方案,采取的具体措施还包括:
进一步地,步骤1中,建立陶瓷基复合材料的三维有限元模型并进行网格划分,提取有限元程序可识别的单元、结点信息;基于有限元原理,组装与损伤无关的形函数矩阵N、几何矩阵B、总体质量矩阵M、等效节点载荷向量F。
进一步地,步骤2中,有限元法中单元刚度矩阵Ke的计算公式为
式中,Ve表示单元体积,V表示体积,D表示弹性矩阵,由陶瓷基复合材料非线性本构关系结合损伤情况实时确定;初始状态下,假设材料未损伤;随着迭代的进行,对每个损伤状态进行刚度的折减,弹性矩阵D中的弹性参数由步骤6不断更新;
运用高斯积分法在单元体积Ve内积分获得单元刚度矩阵Ke,进一步组装成总体刚度矩阵K,作为基于损伤的陶瓷基复合材料刚度模型。
进一步地,步骤3中,求解广义特征值问题,即
Kφ-ω2Mφ=0 (2)
式中,φ表示固有振型,ω表示固有频率;采用C++语言编程并调用Eigen矩阵库,求解式(2)得到各阶固有频率ωi及对应的固有振型向量φi。
进一步地,步骤4中,在Rayleigh阻尼的假设下,将阻尼矩阵C视为刚度矩阵和质量矩阵的线性组合,即
C=αM+βK (3)
已知两阶的固有频率和阻尼比ωa、ζb和ωb、ζb,通过下式求得常数α、β
α和β确定后,组装Rayleigh阻尼矩阵,并运用正则振型使其对角化;则每阶固有频率ωi对应的阻尼比ζi也得以确定,即
其中,刚度矩阵随损伤情况实时变化,阻尼矩阵进行同步动态变化。
进一步地,步骤5中,已知振动系统运动控制方程为
式中,F为系统承受的外激励载荷,a(t)是系统的位移总响应;通过固有振型引入坐标变换
式中,Φ是由各阶固有振型向量φi组成的矩阵,x(t)是各阶广义位移值xi组成的向量;代入式(6)中解耦得到n个二阶非耦合微分方程
对结构施加正弦体积力载荷,即F=a cos ft,a和f分别是正弦载荷的幅值和频率,t表示时间;对于单自由度系统的振动,其解分为瞬态解和稳态解;瞬态解对应着齐次方程的解,会随着时间趋于无穷小;稳态解,系统做和外激励同频的振动;不考虑激励载荷开始时的瞬态响应,求解式(8)获得稳态响应为xi(t)=Xi cos(ft-θi) (9)
各阶位移响应的振幅Xi和相位角θi分别为
式中,ri=f/ωi是频率比;最终系统的位移总响应通过式(7)得到;基于位移响应,运用有限元同步获得应变、应力响应。
进一步地,步骤6中,基于陶瓷基复合材料的非线性应力-应变曲线,判断材料的损伤情况,计算割线刚度
E1=σx t/εx t (11)
式中,σx和εx分别是材料主方向的应力、应变;计算出的割线刚度E1用于更新下一迭代步中的弹性矩阵D;随着应变的变化,若当前应变大于历史最大应变值,损伤程度增大,采用当前应变计算割线刚度;若当前应变小于历史最大应变值,根据损伤不可逆性,采用历史最大应变值计算割线刚度;基于以上原则,实现当前损伤状态下的刚度折减;若损伤饱和,位移、应力、应变值趋于稳定,输出模态和应力、应变,否则重复执行步骤2-步骤6。迭代过程如图1所示。
本发明的有益效果是:本发明提供的一种陶瓷基复合材料模态的非线性计算方法,建立了考虑损伤的刚度和阻尼模型。运用变刚度模态叠加法,并结合有限元,能够快速给出振动载荷下陶瓷基复合材料的模态响应。相比于现有的计算方法,本发明考虑了陶瓷基复合材料的非线性本构,且取得比直接积分法高的计算效率,节省了大量时间。
附图说明
图1是本发明中割线刚度迭代示意图。
图2是本发明中陶瓷基复合材料悬臂梁模型。
图3是本发明中第一阶固有频率迭代过程。
图4是本发明中第一阶固有振型迭代过程。
图5是本发明中弹性参数迭代过程。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
实施例:单向陶瓷基复合材料悬臂梁,长、宽、高分别为0.16、0.015、0.00428m,密度为1771kg/m3。第一阶和第二阶阻尼比ζ1、ζ2分别为1.59×10-3、2.387×10-3。施加正弦体积力,大小为40m/s2,频率为1479Hz。
步骤1:建立陶瓷基复合材料悬臂梁,并划分网格,如图2所示。提取有限元程序可识别的单元、结点等信息。基于有限元原理,组装与损伤无关的形函数矩阵、几何矩阵、质量矩阵、载荷矩阵。
步骤2:复合材料的等效弹性参数如下,代表无损伤状态:
E1=137.37GPa,E2=E3=63.16GPa,G12=G13=20.62GPa,G23=24.61GPa,υ12=υ13=0.27,υ23=0.25。
单向陶瓷基复合材料可视为横观各向同性材料,其柔度矩阵为
对柔度矩阵求逆,可以得到弹性矩阵D=S-1。运用高斯积分法计算式(1)中的单元刚度矩阵,并组装结构刚度矩阵。其中拉伸方向的弹性参数E1由步骤6确定损伤状态实时更新。
步骤3:求解广义特征值问题,可获得初始未损伤状态的模态。随着损伤的增大,刚度、固有频率同步降低,最终迭代获得稳态的模态数据。迭代过程如图3~图4所示。
步骤4:由式(4)计算出初始阻尼系数α=2.93、β=8.153×10-7,组装Rayleigh阻尼,并运用正则振型使其对角化。则每阶固有频率ωi对应的阻尼比ζi也可确定。刚度矩阵随损伤情况实时变化,阻尼矩阵进行同步动态变化。
步骤5:运用振型叠加法,解耦振动微分方程。由式(8)-(10),获得当前状态下的稳态响应。
步骤6:基于陶瓷基复合材料的非线性细观本构曲线,判断材料的当前损伤状况。根据式(11)求解割线刚度,用于计算步骤2中的弹性矩阵D。随着应变的变化,若当前应变大于历史最大应变值,损伤程度增大,采用当前应变计算割线刚度;若当前应变小于历史最大应变值,根据损伤不可逆性,采用历史最大应变值计算割线刚度。基于以上原则,实现当前损伤状态下的刚度折减。若损伤饱和,位移、应力、应变值趋于稳定,输出模态和各响应,否则执行步骤2。其中,弹性参数E1的迭代过程如图5所示。
需要注意的是,发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (4)
1.一种陶瓷基复合材料模态的非线性计算方法,其特征在于,包括如下步骤:
步骤1、建立陶瓷基复合材料的三维有限元模型,获得单元、结点信息;步骤1中,建立陶瓷基复合材料的三维有限元模型并进行网格划分,提取有限元程序可识别的单元、结点信息;基于有限元原理,组装与损伤无关的形函数矩阵N、几何矩阵B、总体质量矩阵M、等效节点载荷向量F;
步骤2、运用有限元法建立基于损伤的陶瓷基复合材料刚度模型;
步骤2中,有限元法中单元刚度矩阵Ke的计算公式为
式中,Ve表示单元体积,V表示体积,D表示弹性矩阵,由陶瓷基复合材料非线性本构关系结合损伤情况实时确定;
运用高斯积分法在单元体积Ve内积分获得单元刚度矩阵Ke,进一步组装成总体刚度矩阵K,作为基于损伤的陶瓷基复合材料刚度模型;
步骤3、基于步骤1和步骤2求解广义特征值问题,获得固有频率和固有振型;
步骤4、由步骤2建立的刚度模型对Rayleigh阻尼进行同步更新变化,最终建立基于损伤的陶瓷基复合材料阻尼模型;
步骤5、将步骤3获得的模态结果代入步骤2建立的刚度模型和步骤4建立的阻尼模型,运用振型叠加法计算陶瓷基复合材料在振动载荷下的位移、应力、应变响应;
步骤6、基于陶瓷基复合材料的非线性应力-应变曲线,判断材料的损伤情况,若损伤饱和,输出模态和应力、应变,否则更新弹性参数,重复执行步骤2至步骤6;
步骤6中,基于陶瓷基复合材料的非线性应力-应变曲线,判断材料的损伤情况,计算割线刚度
E1=σx t/εx t (11)
式中,σx和εx分别是材料主方向的应力、应变;计算出的割线刚度E1用于更新下一迭代步中的弹性矩阵D;随着应变的变化,若当前应变大于历史最大应变值,损伤程度增大,采用当前应变计算割线刚度;若当前应变小于历史最大应变值,根据损伤不可逆性,采用历史最大应变值计算割线刚度;基于以上方法,实现当前损伤状态下的刚度折减;若损伤饱和,位移、应力、应变值趋于稳定,输出模态和应力、应变,否则重复执行步骤2-步骤6。
2.如权利要求1所述的一种陶瓷基复合材料模态的非线性计算方法,其特征在于:步骤3中,求解广义特征值问题,即
Kφ-ω2Mφ=0 (2)
式中,φ表示固有振型,ω表示固有频率;采用C++语言编程并调用Eigen矩阵库,求解式(2)得到各阶固有频率ωi及对应的固有振型向量φi。
4.如权利要求3所述的一种陶瓷基复合材料模态的非线性计算方法,其特征在于:步骤5中,已知振动系统运动控制方程为
式中,a(t)是系统的位移总响应;通过固有振型引入坐标变换
式中,Φ是由各阶固有振型向量φi组成的矩阵,x(t)是各阶广义位移值xi组成的向量;代入式(6)中解耦得到n个二阶非耦合微分方程
对结构施加正弦体积力载荷,即F=a cos ft,a和f分别是正弦载荷的幅值和频率,t表示时间;不考虑激励载荷开始时的瞬态响应,求解式(8)获得稳态响应为
xi(t)=Xicos(ft-θi) (9)
各阶位移响应的振幅Xi和相位角θi分别为
式中,ri=f/ωi是频率比;最终系统的位移总响应通过式(7)得到;基于位移响应,运用有限元同步获得应变、应力响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138360.7A CN110852015B (zh) | 2019-11-19 | 2019-11-19 | 一种陶瓷基复合材料模态的非线性计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911138360.7A CN110852015B (zh) | 2019-11-19 | 2019-11-19 | 一种陶瓷基复合材料模态的非线性计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110852015A CN110852015A (zh) | 2020-02-28 |
CN110852015B true CN110852015B (zh) | 2021-11-02 |
Family
ID=69602581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911138360.7A Active CN110852015B (zh) | 2019-11-19 | 2019-11-19 | 一种陶瓷基复合材料模态的非线性计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110852015B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005272192A (ja) * | 2004-03-24 | 2005-10-06 | National Institute Of Advanced Industrial & Technology | 多孔質構造体の製造方法 |
DE102006035247A1 (de) * | 2006-07-26 | 2008-02-07 | Grob-Werke Gmbh & Co. Kg | Thermomechanisches Stellelement |
CN103528901A (zh) * | 2013-10-23 | 2014-01-22 | 东北大学 | 一种测试硬涂层复合结构非线性刚度及阻尼的装置及方法 |
CN108287970A (zh) * | 2018-01-31 | 2018-07-17 | 东南大学 | 基于二维正交各向异性复合材料板的热模态对结构参数的灵敏度分析方法 |
CN108897213A (zh) * | 2018-05-28 | 2018-11-27 | 宝鸡文理学院 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
CN110006612A (zh) * | 2019-04-24 | 2019-07-12 | 哈尔滨工业大学 | 一种非接触式热模态试验系统及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416084B (zh) * | 2018-01-23 | 2022-02-18 | 南京理工大学 | 考虑复合材料弹塑性与损伤耦合的弹塑性损伤有限元方法 |
CN109920495B (zh) * | 2019-03-28 | 2020-03-13 | 南京航空航天大学 | 一种编织陶瓷基复合材料强度的多尺度预测方法 |
-
2019
- 2019-11-19 CN CN201911138360.7A patent/CN110852015B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005272192A (ja) * | 2004-03-24 | 2005-10-06 | National Institute Of Advanced Industrial & Technology | 多孔質構造体の製造方法 |
DE102006035247A1 (de) * | 2006-07-26 | 2008-02-07 | Grob-Werke Gmbh & Co. Kg | Thermomechanisches Stellelement |
CN103528901A (zh) * | 2013-10-23 | 2014-01-22 | 东北大学 | 一种测试硬涂层复合结构非线性刚度及阻尼的装置及方法 |
CN108287970A (zh) * | 2018-01-31 | 2018-07-17 | 东南大学 | 基于二维正交各向异性复合材料板的热模态对结构参数的灵敏度分析方法 |
CN108897213A (zh) * | 2018-05-28 | 2018-11-27 | 宝鸡文理学院 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
CN110006612A (zh) * | 2019-04-24 | 2019-07-12 | 哈尔滨工业大学 | 一种非接触式热模态试验系统及方法 |
Non-Patent Citations (4)
Title |
---|
ESTIMATION OF SONIC FATIGUE BY REDUCED-ORDER FINITE ELEMENT BASED ANALYSES;Rizzi Stephen A.,et al.;《International Conference on Recent Advances in Structural Dynamics》;20060719;129-145 * |
FRP约束钢筋混凝土柱非线性阻尼性能及动力响应研究;李笑然;《中国博士学位论文全文数据库(工程科技Ⅱ辑)》;20140115(第01期);C038-14 * |
Numerical Simulation of Dynamic Response of Fiber Reinforced Ceramic Matrix Composite Beam with Matrix Cracks Using Multiscale Modeling;Gao Xiguang, et al.;《Chinese Journal of Aeronautics》;20100831;第23卷(第5期);537-548 * |
陶瓷基复合材料本构模型及振动分析;陈晶;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20180315(第03期);C031-37 * |
Also Published As
Publication number | Publication date |
---|---|
CN110852015A (zh) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chuang et al. | Multi-scale design and optimization for solid-lattice hybrid structures and their application to aerospace vehicle components | |
CN107103138B (zh) | 一种激光喷丸变刚度轻量化方法 | |
Kan et al. | Nonlinear dynamic and deployment analysis of clustered tensegrity structures using a positional formulation FEM | |
Pagani et al. | Exact dynamic stiffness elements based on one-dimensional higher-order theories for free vibration analysis of solid and thin-walled structures | |
Loja et al. | Dynamic instability of variable stiffness composite plates | |
Bottasso et al. | Integrated aero-structural optimization of wind turbines | |
CN107563005B (zh) | 一种拉压不同刚度结构振动瞬时最优控制方法 | |
Devarajan | Free vibration analysis of curvilinearly stiffened composite plates with an arbitrarily shaped cutout using isogeometric analysis | |
CN110162826A (zh) | 薄壁结构热气动弹性动响应分析方法 | |
Doeva et al. | Exact solution for the deflection of composite beams under non-uniformly distributed loads | |
Battista et al. | Large deformations of 1D microstructured systems modeled as generalized Timoshenko beams | |
Wang et al. | Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems | |
Yan et al. | Structure/material concurrent optimization of lattice materials based on extended multiscale finite element method | |
Zhang et al. | An exact Fourier series method for vibration analysis of elastically connected laminated composite double-beam system with elastic constraints | |
CN110852015B (zh) | 一种陶瓷基复合材料模态的非线性计算方法 | |
Cea et al. | Nonlinear modal aeroelastic analysis from large industrial-scale models | |
WO2021077900A1 (zh) | 一种陶瓷基复合材料流固耦合响应计算方法 | |
Kirsch et al. | Accurate displacement derivatives for structural optimization using approximate reanalysis | |
CN111737908A (zh) | 一种基于动载荷静力等效的蒙皮桁条结构快速动态优化设计方法 | |
CN116187074A (zh) | 基于等几何的各向异性周期结构材料多尺度拓扑优化方法 | |
Behera et al. | Wave propagation analysis in laminated composite periodic frame structures using spectral finite element method | |
Vidal et al. | Modeling of composite and sandwich beams with a generic cross-section using a variable separation method | |
Ojo et al. | Geometrically nonlinear analysis of non-prismatic beam structures using strong Unified Formulation | |
Yu et al. | Orthogonal polynomials-ritz method for dynamic response of functionally graded porous plates using FSDT | |
Foutsitzi et al. | Modal shape control of smart composite beams using piezoelectric actuators |
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 |