CN111551434B - 一种碳纤维复合材料轴向压缩性能测试方法与本构模型 - Google Patents
一种碳纤维复合材料轴向压缩性能测试方法与本构模型 Download PDFInfo
- Publication number
- CN111551434B CN111551434B CN202010366881.4A CN202010366881A CN111551434B CN 111551434 B CN111551434 B CN 111551434B CN 202010366881 A CN202010366881 A CN 202010366881A CN 111551434 B CN111551434 B CN 111551434B
- Authority
- CN
- China
- Prior art keywords
- composite material
- carbon fiber
- strain
- fiber composite
- compression
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/30—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight
- G01N3/307—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight generated by a compressed or tensile-stressed spring; generated by pneumatic or hydraulic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/003—Generation of the force
- G01N2203/0042—Pneumatic or hydraulic means
- G01N2203/0044—Pneumatic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0067—Fracture or rupture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种碳纤维复合材料轴向压缩性能测试方法与本构模型:首先,采用MTS万能试验机对碳纤维复合材料试件进行轴向的准静态压缩,得到低应变率压缩性能;然后,通过SHPB实验装置对碳纤维复合材料试件进行轴向的动态压缩实验,得到高应变率压缩性能;接着,采用复合材料混合定律验证准静态实验的合理性,采用波形变化关系验证动态压缩实验的准确性,提出联合型屈曲强度计算公式;最后,引入应变率项对实验结果进行拟合,得到碳纤维复合材料轴向压缩的应力应变本构模型。本发明得到的本构模型能较为准确地描述碳纤维复合材料轴向压缩的力学性能,且能为数值仿真和其他类似材料的轴向压缩性能测试提供参考。
Description
技术领域
本发明属于压缩性能研究领域,具体涉及一种碳纤维复合材料轴向压缩性能测试方法与本构模型。
背景技术
环氧树脂是碳纤维增强聚合物复合材料中最常用的热固性聚合物,将其作为基体制成的碳纤维增强环氧树脂复合材料具有模量高、强度高、抗拉强度高、耐疲劳、耐高温、耐腐蚀、耐电击、密度小、可设计强等优点,广泛地应用于汽车、船舶、航空、航天、建筑等领域,被称为“黑色黄金”。碳纤维的发展非常迅速,目前常见的有T300、T700、T800等系列。近年来,更是出现了含有铺层的、多向编织的、多种碳纤维的复合材料。学者们的研究转向于含铺层的本构研究,并且加入了连续损伤效应,这样使得研究通常局限于现象描述和更新结构,然而对基础的单向碳纤维的压缩本构模型建立尚有不足。
曹翠微等[曹翠微,李照谦,李贺军,等.轴棒法编织C/C复合材料的热物理及烧蚀性能[J].固体火箭技术,2011,34(1):113-118.]对轴棒法编织三维四向C/C复合材料压缩及弯曲性能做了研究,观测了材料的压缩应力-应变曲线,给出了材料的失效过程,但是并未对本构模型进行研究。李彪等[李彪,杨勇新,岳清瑞,等.碳纤维复合材料棒材力学性能试验方法研究[J].工业建筑,2013,43(6):5-8.]通过实验研究了碳纤维长棒的拉伸性能,指出长碳纤维棒的拉伸破坏模式为脆性断裂和爆炸式断裂,并绘制了拉伸应力应变曲线;C.SOUTIS[Soutis C.Compression testing of pultruded carbon fibre-epoxycylindrical rods[J].Journal of Materials Science,2000,35(14):3441-3446.]开展了挤压碳纤维-环氧圆柱杆的压缩试验,材料为T300/828,直径4mm,有效长度13mm,压缩速率1mm/min;这两个文献的相同之处在于,试件是长尺寸试件且只进行了静态的拉伸和压缩实验,未对本构模型进行研究。
本构模型是描述一类材料的力学性能的基本方程,是材料力学特性研究的重中之重,是数值计算发展的核心要点。为了弥补和丰富单向碳纤维复合材料的压缩力学领域,需要提出一种真实应力应变状态下的本构模型。本发明基于MTS和SHPB实验设备,开展了碳纤维复合材料轴向压缩性能测试,提出了一种本构模型。
发明内容
本发明的目的在于提供一种碳纤维复合材料轴向压缩性能测试方法与本构模型,以建立碳纤维复合材料轴向压缩的本构模型,为数值仿真和其他类似材料的轴向压缩性能测试提供参考。
实现本发明目的的技术解决方案为:一种碳纤维复合材料轴向压缩性能测试方法与本构模型,包括以下步骤:
步骤1、采用MTS万能试验机对碳纤维复合材料试件(直径10mm×长度10mm的圆柱)进行轴向的准静态压缩,得到低应变率压缩性能;
步骤2、通过SHPB实验装置对碳纤维复合材料试件(直径10mm×长度5mm的圆柱)进行轴向的动态压缩实验,得到高应变率压缩性能;
步骤3、采用复合材料混合定律验证准静态实验的合理性,采用波形变化关系验证动态压缩实验的准确性,提出联合型屈曲强度计算公式;
步骤4、引入应变率项对实验结果进行拟合,得到碳纤维复合材料轴向压缩的应力应变本构模型。
本发明与现有技术相比,其显著优点是:
(1)本发明采用准静态和动压压缩实验方法,获得了碳纤维复合材料在低应变率和高应变率下的轴向压缩性能。
(2)本发明通过复合材料混合定律和波形变化关系,验证了准静态实验和动态压缩实验的合理性和准确性;
(3)本发明通过实验观察了试件的失效模式,综合剪切型屈曲和拉压型屈曲的特点,提出了联合型屈曲强度计算公式;
(4)本发明通过引入应变率项,结合实验数据,拟合得到了一种描述碳纤维复合材料轴向真实应力应变关系的本构模型。
附图说明
图1为本发明方法的流程图。
图2为实施例准静态和动态真实应力应变曲线图,其中图(a)为准静态应力应变曲线图,图(b)为动态真实应力应变曲线图。
图3为实施例动态压缩波形验证图。
图4为实施例准静态压缩本构方程验证图,其中图(a)为应变率为0.1s-1的准静态应力应变曲线图,图(b)为应变率为0.01s-1的准静态应力应变曲线图,图(c)为应变率为0.001s-1的准静态应力应变曲线图。
图5为实施例动态压缩本构方程验证图,其中图(a)为应变率为537s-1的动态应力应变曲线图,图(b)为应变率为1158s-1的动态应力应变曲线图,图(c)为应变率为1351s-1的动态应力应变曲线图。
具体实施方式
下面结合附图及具体实施例对本发明做进一步的介绍。
结合图1,本发明的一种碳纤维复合材料轴向压缩性能测试方法与本构模型,包括以下步骤:
步骤1、采用MTS万能试验机对碳纤维复合材料试件进行轴向的准静态压缩,得到低应变率压缩性能:
通过MTS电子万能试验机进行低应变率实验,实验时加装引伸计,能获得更准确的压缩变形。通过改变压头速度v,可以得到材料在应变率为0.001s-1、0.01s-1、0.1s-1三中工况下的压缩特性,实验前进行0.1KN的预压。v通过工程应变率的计算公式(1)反推得到。处理实验获得的力-位移曲线,根据式(2)可得到工程应力σ和应变ε,再根据关系式(3)可转换成真实的应力σT和应变εT。
式中,l0为试件的初始长度,F为轴向载荷,A0为试件的初始横截面积,Δl为压缩距离。
步骤2、通过SHPB实验装置对碳纤维复合材料试件进行轴向的动态压缩实验,得到高应变率压缩性能:
通过分离式霍普金森压杆(split Hopkinson pressure bar)SHPB实验装置进行高应变率实验。SHPB实验的建立是基于一维应力波假设和均匀性假设,将试样的压缩状态视为一维应力状态,并认为应力应变沿长度方向均匀分布。控制应变率的方法是控制子弹的弹射气压,从而控制子弹的撞击速度。实验前,需要对齐子弹与入射杆、入射杆与透射杆,减小弥散效应;在试件两端涂抹凡士林,减小摩擦效应。SHPB实验装置获取应力应变关系的方法有三种:一波法、两波法和三波法。本文采用反射波和透射波组成的两波法,计算公式为
式中:σ为试件的一维应力(工程应力),A为当前试件横截面积,E为试件的弹性模量,A0为原始试件横截面积,ε为工程应变(工程应变),εt和εr为分别为透射波和反射波独立传播过程中的应变,c0为压杆中应力波的传播速度,l0为试件的原始长度,l为试件的当前长度,τ为时间,为应变率。
步骤3、采用复合材料混合定律验证准静态实验的合理性,采用波形变化关系验证动态压缩实验的准确性:
根据混合定律,取下标f代表纤维(fiber)的参数,m代表基体(matrix)的参数,v为体积分数,ρ为密度,E为弹性模量,则单向碳纤维复合材料的密度、弹性模量可表示为
ρ=ρfvf+ρmvm (5)
对反射波和透射波进行加法计算,与入射波进行对比,以验证应力均匀分布假设是否成立,参见公式(7)。应变值来自于采集信号的换算,因此可以直接采用实验采集得到的信号直接计算。
εt+εr=εi (7)
式中,εi为入射波独立传播过程中的应变。
根据实验过程中观察到的屈曲情况,建立包含拉压型屈曲强度σcr,T-C和剪切型屈曲强度σcr,S的联合型屈曲强度σcr,(T-C)+S计算公式
σcr,(T-C)+S=(vfσcr,T-C+vmσcr,S) (9)
式中,Gm为基体的剪切模量,η为匹配系数。
步骤4、引入应变率项对实验结果进行拟合,得到碳纤维复合材料轴向压缩的应力应变本构模型:
低应变情况下,材料的压缩过程主要包括线弹性阶段和失效后的三个阶段;高应变率情况下,材料的压缩过程主要分为线弹性阶段和波动阶段。低应变实验得出应力应变结果与动态压缩得到的结果一致性不是很好,需要将低应变率与高应变率情况分开讨论。
式中,a,b,c,a0,a1,a2,a3,a4,a5,a6,a7,a8,a9均为实验常数。以应变率0.1s-1为参考进行拟合,得出实验常数,再对应变率为0.01s-1和0.001s-1的实验数据进行验证。
4.2、动态(高应变率)情况:假设弹性阶段为线性关系,发生屈服后的流动应力为多项式曲线关系,采用分段分析方法作进一步的研究。弹性阶段采用广义Hook定律描述,弹性模量通过屈服时的应力σYD和应变εYD计算得到,取应变率为537s-1的屈服应变为参考应变参考εref,得到动态应力σD、动态应变εD、动态应变率相关的本构方程为
式中,ε′D为转换应变,A0,A1,A2,A3,A4,A,B均为实验常数。以应变率为537s-1的实验数据作为参考进行拟合,得出实验常数,再对应变率为1158s-1和1351s-1的实验数据进行验证。
实施例
步骤1和步骤2中所用的材料为碳纤维复合材料,是通过75%的碳纤维束与25%的环氧树脂经高温高压固化而成。纤维是3K/T300聚丙烯腈基碳纤维,拉伸强度为3530MPa,弹性模量为230GPa,直径为7μm,密度1.76g/cm3。树脂为E-51环氧树脂,弹性模量为1.99GPa,密度1.16g/cm3。制成的棒材属于单向碳纤维增强环氧树脂复合材料,密度1.55g/cm3。压缩试样通过截取碳纤维棒得到,MTS准静态实验:圆截面直径d=10mm,长度l=10mm;SHPB动压压缩实验:圆截面直径d=10mm,长度l=5mm。
准静态压缩和SHPB动态压缩结果参见图2。(1)准静态轴向压缩性能:从图中可以看出,三种应变率下碳纤维复合材料的真实应力应变曲线走势基本相同,可大致划分为线弹性阶段、屈服阶段、平台应力阶段、密实压缩阶段。在应变率为0.01s-1和0.001s-1时,具有双线性阶段和双屈服阶段,且流动应力有阶跃下降的特点。随着应变率的增大,屈服强度和流动应力增大,屈服点滞后,具有应变率效应。(2)动态轴向压缩性能:随着应变率的增大,材料的屈服点出现滞后,屈服极限增大,流动应力也增大。结合准静态实验现象,可以确定该材料具有应变率效应。四种高应变率下的应力应变曲线走势相同,可分为两个阶段,即线性上升阶段和流动应力波动阶段。与准静态情况相同,波动阶段出现的原因是基体失去承载能力,载荷转而由纤维承受,但纤维的屈曲、断裂、剪切等复杂不稳定情况使得流动应力发生波动。
根据步骤3,计算得到的密度为1.61g/cm3,与实际密度1.55g/cm3的相对误差仅为3.87%,验证了混合定律的合理性;计算得到的理论弹性模量为8.57GPa,准静态实验测得的第一弹性模量均值为8.32GPa,实验相对理论值的误差在2.92%以内,证明了实验的准确性;此外,由于工艺过程中不可避免的空隙,导致纤维与基体的体积分数存在一定差异,使得理论计算结果偏大。动态验证结果参见图3,从计算结果可以得知,入射波幅值h1与反射波和透射波幅值之和h2基本相同,证明了应力沿长度方向是均匀分布的。在联合型屈曲强度计算过程中,材料的剪切模量G=1GPa,匹配系数η=0.01;单独采用剪切型屈曲强度计算方法得到的结果为40MPa,单独采用拉压型屈曲强度计算方法得到的结果为321MPa,采用联合区屈曲强度计算方法得到的结果为250.75MPa;实验得到的屈曲强度为246.41MPa,与联合型屈曲强度预测结果相近,两者的相对误差仅为1.76%,从而极大地提高了预测精度。
根据步骤4,拟合得出的准静态实验常数为:a=1.4,b=13.6,c=2.9×107,a0=0,a1=143.5,a2=-2105.2,a3=13709.7,a4=-49178.4,a5=105311.4,a6=-137963.8,a7=108379,a8=-46868.5,a9=8575.9。准静态压缩的真实应力应变计算结果参见图4,可以看出,对于应变率为0.01s-1和0.001s-1的情况,这个本构方程得出的曲线与实验曲线走势基本吻合,特别是在初始上升阶段和平台阶段效果更好。
拟合得出的动态实验常数为:A0=-7.3×107,A1=1.2×1011,A2=-1.2×1013,A3=4.6×1014,A4=-6×1015,A=0.3,B=159.1。动态压缩的真实应力应变曲线参见图5,可以看出,对于应变率为1153s-1和1358s-1的情况,线性段的的结果是吻合的非常好的,流动应力波动段也具有一定的匹配性,基本能描述单向碳纤维复合材料的动态压缩应力应变行为。
Claims (3)
1.一种碳纤维复合材料轴向压缩性能测试方法,其特征在于,包括以下步骤:
步骤1、采用MTS万能试验机对碳纤维复合材料试件进行轴向的准静态压缩,得到碳纤维复合材料试件的低应变率压缩性能;
步骤2、通过SHPB实验装置对碳纤维复合材料试件进行轴向的动态压缩实验,得到碳纤维复合材料试件的高应变率压缩性能;
步骤3、采用复合材料混合定律验证准静态实验的合理性,采用波形变化关系验证动态压缩实验的准确性,提出联合型屈曲强度计算公式:
式中,σcr(T-C)+S为联合型屈曲强度,η为匹配系数,vm、vf分别为基体和纤维的体积分数,Em、Ef分别为基体和纤维的弹性模量,Gm为基体的剪切模量;
步骤4、引入应变率项对实验结果进行拟合,得到碳纤维复合材料轴向压缩的应力应变本构模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010366881.4A CN111551434B (zh) | 2020-04-30 | 2020-04-30 | 一种碳纤维复合材料轴向压缩性能测试方法与本构模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010366881.4A CN111551434B (zh) | 2020-04-30 | 2020-04-30 | 一种碳纤维复合材料轴向压缩性能测试方法与本构模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111551434A CN111551434A (zh) | 2020-08-18 |
CN111551434B true CN111551434B (zh) | 2022-11-22 |
Family
ID=72000173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010366881.4A Active CN111551434B (zh) | 2020-04-30 | 2020-04-30 | 一种碳纤维复合材料轴向压缩性能测试方法与本构模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111551434B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112304756B (zh) * | 2020-10-16 | 2023-07-21 | 中国航发四川燃气涡轮研究院 | 一种纤维增强复合材料环状结构周向拉伸性能表征方法 |
CN113125259B (zh) * | 2021-04-12 | 2022-10-28 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种点阵结构准静态压缩试验件及试验方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106053222A (zh) * | 2016-06-07 | 2016-10-26 | 大连理工大学 | 一种铝硅合金ADC12材料Johnson‑Cook本构模型的拟合方法 |
CN106202683A (zh) * | 2016-07-04 | 2016-12-07 | 南京理工大学 | 一种考虑相关效应的橡胶材料粘超弹本构模型的建模方法 |
CN106525612B (zh) * | 2016-09-23 | 2019-03-26 | 西北工业大学 | 基于拉伸及压缩实验的聚碳酸酯统一本构模型的构建方法 |
CN107633120B (zh) * | 2017-09-07 | 2018-07-17 | 东南大学 | 一种纤维增强复合材料动态剪切本构模型的构建方法 |
CN108549743B (zh) * | 2018-03-13 | 2019-03-12 | 东南大学 | 一种纤维增强复合材料动态拉伸失效评估方法 |
CN109902362B (zh) * | 2019-02-03 | 2023-08-25 | 南京理工大学 | 考虑纤维增强复合材料复杂非线性行为的弹粘塑性本构模型构建方法 |
-
2020
- 2020-04-30 CN CN202010366881.4A patent/CN111551434B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111551434A (zh) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sockalingam et al. | Fiber-matrix interface characterization through the microbond test | |
Daoud et al. | Characterization of the vibrational behaviour of flax fibre reinforced composites with an interleaved natural viscoelastic layer | |
CN111551434B (zh) | 一种碳纤维复合材料轴向压缩性能测试方法与本构模型 | |
Zarei et al. | On the free vibrations of joined grid-stiffened composite conical-cylindrical shells | |
Rashedi et al. | Fracture characterization of glass fiber composite laminate under experimental biaxial loading | |
Jumahat et al. | A graphical method predicting the compressive strength of toughened unidirectional composite laminates | |
Daniel et al. | Characterization and constitutive modeling of composite materials under static and dynamic loading | |
Lu et al. | A novel experiment-based approach on nonlinear bending analysis of thermoplastic composite laminates | |
Kroupa et al. | Improved nonlinear stress-strain relation for carbon-epoxy composites and identification of material parameters | |
Han et al. | Behaviors of composite laminates under low-energy impact using a novel analytical framework | |
Catangiu et al. | Experimental results for bending fatigue behaviour of glass-epoxy composite materials | |
Wang et al. | Numerical modeling of tensile behavior of fiber-reinforced polymer composites | |
Samal et al. | Effect of fiber orientation on vibration response of glass epoxy composite beam | |
Makeev et al. | A method for measurement of three-dimensional constitutive properties for composite materials | |
Ji et al. | Study of in-plane fatigue failure and life prediction of weave composites under constant and variable amplitude loading | |
Böhm et al. | Experimental investigation of the strain rate dependent behaviour of 2D biaxially and triaxially reinforced braided composites | |
CN111143975B (zh) | 一种锻造的碳纤维复合材料压缩力学本构计算方法 | |
Prabhakaran et al. | Modelling and simulation of natural fibre/epoxy composites-prediction of stress state and deformations | |
Changfang et al. | INVESTIGATION OF COMPRESSION MECHANICS OF STRAIN RATE-DEPENDENT FORGED/LAMINATED CARBON FIBER− EPOXY RESIN COMPOSITES | |
Njim et al. | Flexural Bending and Fatigue Analysis of Functionally Graded Viscoelastic Materials: Experimental and Numerical Approaches | |
Gowtham et al. | Experimental investigation of vibration characteristics of polymer composites subjected to edge crack | |
Deshpande | Characterization of CFRP and GFRP composite materials at high strain rate tensile loading | |
Azhaguvel et al. | Optimization of mechanical properties of e-glass woven fabric composite | |
Makris et al. | Biaxial failure envelopes for glass fibre reinforced composite laminates | |
Bao et al. | A study on the torsional wave of fiber-reinforced composite materials |
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 |