CN113866008B - 一种基于阈值应力和抗拉强度的蠕变寿命预测方法 - Google Patents

一种基于阈值应力和抗拉强度的蠕变寿命预测方法 Download PDF

Info

Publication number
CN113866008B
CN113866008B CN202111040892.4A CN202111040892A CN113866008B CN 113866008 B CN113866008 B CN 113866008B CN 202111040892 A CN202111040892 A CN 202111040892A CN 113866008 B CN113866008 B CN 113866008B
Authority
CN
China
Prior art keywords
stress
creep
temperature
threshold stress
tensile strength
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
Application number
CN202111040892.4A
Other languages
English (en)
Other versions
CN113866008A (zh
Inventor
宋迎东
赵旭
孙志刚
牛序铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202111040892.4A priority Critical patent/CN113866008B/zh
Publication of CN113866008A publication Critical patent/CN113866008A/zh
Application granted granted Critical
Publication of CN113866008B publication Critical patent/CN113866008B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G01N3/18Performing tests at high or low temperatures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001Type of application of the stress
    • G01N2203/0003Steady
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0017Tensile
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0071Creep
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0073Fatigue
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/022Environment of the test
    • G01N2203/0222Temperature
    • G01N2203/0226High temperature; Heating means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0676Force, weight, load, energy, speed or acceleration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0694Temperature

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)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

本发明的公开了一种基于阈值应力和抗拉强度的蠕变寿命预测方法,包括以下步骤:步骤1,进行不同温度T下的材料高温拉伸试验,得到对应温度下的抗拉强度σb;步骤2,在不同温度下分别开展多组不同应力条件下的高温蠕变试验,获得对应温度应力条件下的最小蠕变速率以及蠕变寿命tf;步骤3,根据最小蠕变速率得到不同温度对应的阈值应力σth;步骤4,基于不同温度下的阈值应力σth,建立阈值应力与温度之间的函数关系;步骤5,在阈值应力、抗拉强度基础上,建立蠕变寿命预测公式;通过蠕变寿命预测公式,对不同应力温度条件下的蠕变寿命进行预测。本发明考虑了阈值应力和抗拉强度的共同影响,寿命预测精度高,并具有较强的外推能力。

Description

一种基于阈值应力和抗拉强度的蠕变寿命预测方法
技术领域
本发明涉及一种蠕变寿命预测方法,尤其涉及一种基于阈值应力和抗拉强度的蠕变寿命预测方法。
背景技术
随着低碳排放以及高燃烧效率要求的不断提高,高温部件承受着日益严重的蠕变载荷影响,准确预测其蠕变寿命至关重要。由于高温合金优异的抗蠕变性能,使得高温合金被广泛应用于高温部件。高温合金优异的抗蠕变性能主要来源于大量的沉淀强化相与位错的相互作用,通过强化相阻碍位错的运动,以限制蠕变变形的进一步发展。强化相的存在使得沉淀强化类以及固溶强化类等合金存在阈值应力特征,当外加应力小于该阈值应力时,位错难以克服强化相的阻碍继续运动,所以难以产生蠕变变形,此时可以认为蠕变寿命为无穷大。而当应力趋向于抗拉强度时,蠕变寿命趋近于0。可以看出,蠕变寿命与阈值应力与抗拉强度之间存在的固有的联系。然而,目前却没有研究将阈值应力和抗拉强度相结合以实现对蠕变寿命的准确预测。
发明内容
本发明的目的是提供一种基于阈值应力和抗拉强度的蠕变寿命预测方法,以实现对蠕变寿命的准确预测。
为实现上述目的,本发明采用如下技术方案:
一种基于阈值应力和抗拉强度的蠕变寿命预测方法,包括以下步骤:
步骤1,进行不同温度T下的材料高温拉伸试验,得到对应温度下的抗拉强度σb
步骤2,在不同温度下分别开展多组不同应力条件下的高温蠕变试验,获得对应温度应力条件下的最小蠕变速率以及蠕变寿命tf
步骤3,根据步骤2获得的最小蠕变速率得到不同温度对应的阈值应力σth
步骤4,基于步骤3获得的不同温度下的阈值应力σth,建立阈值应力与温度之间的函数关系;
步骤5,在步骤4得到的阈值应力、步骤1得到的抗拉强度基础上,建立蠕变寿命预测公式;通过蠕变寿命预测公式,对不同应力温度条件下的蠕变寿命进行预测。
所述步骤3中,根据步骤2中高温蠕变试验获得的最小蠕变速率利用公式建立同一温度下最小蠕变速率/>应力σ与阈值应力σth之间的关系,从而求解该温度下对应的阈值应力σth,式中,Am为常数;对不同温度进行同样操作,进而得到不同温度对应的阈值应力水平。
所述步骤4中,基于步骤3获得的不同温度下的阈值应力水平σth,利用多项式形式进行拟合,从而建立阈值应力与温度之间的函数关系,即:式中,m为多项式项数,ai为拟合参数,i=0,1,2···,m,一般可取m≤3。
所述步骤5中,建立的蠕变寿命预测公式为:
式中,A、n为常数,σth为阈值应力,σb为抗拉强度,σ为施加的应力,T为施加的温度,R为气体常数,为表观激活能。
所述步骤5中,表观激活能由以下方法得到:在同一/>值下,由试验蠕变速率的对数/>与温度的倒数1/T之间的线性拟合直线斜率确定表观激活能/>
所述步骤5中,常数A、n由以下方法得到:对蠕变寿命预测公式进行转换,两边取对数,有:
因此,常数A、n通过的线性拟合直线的斜率与截距获得。
所述步骤5中,通过蠕变寿命预测公式,对不同应力温度条件下的蠕变寿命进行预测时,应当保证当应力趋向于阈值应力σth时,蠕变寿命趋近于无穷大,当应力趋向于抗拉强度σb时,蠕变寿命趋近于0。
有益效果:与现有技术相比,本发明具有以下技术效果:
1、本发明提出的一种基于阈值应力和抗拉强度的蠕变寿命预测方法,只需进行高温拉伸试验与高温蠕变试验,获得相应温度下的抗拉强度以及阈值应力,所需参数少,试验简单,成本低;
2、本发明建立了阈值应力、抗拉强度以及蠕变寿命之间的关系,将蠕变寿命限制在阈值应力和抗拉强度所对应的应力区间之内,当应力趋向于阈值应力时,蠕变寿命趋近于无穷大,当应力趋向于抗拉强度时,蠕变寿命趋近于0,该方法对蠕变寿命预测精度高,更具有理论依据;
3、本发明可适用于工程中常见的沉淀强化合金、固溶强化合金等材料,具有广泛的应用能力。
附图说明
图1是阈值应力计算方法图;
图2是对数最小蠕变速率与温度的倒数关系图;
图3是线性拟合求解常数A、n图;
图4是蠕变寿命预测值随应力温度的变化关系图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
本发明的一种基于阈值应力和抗拉强度的蠕变寿命预测方法,包括以下步骤:
步骤1,首先开展在不同温度T下的材料高温拉伸试验,得到对应温度下的抗拉强度σb
步骤2,在不同温度下进行多组应力下的高温蠕变试验,一般可开展2-4个温度,每个温度下5-7组不同应力的高温蠕变试验;获得对应温度应力条件下的最小蠕变速率以及蠕变寿命tf
步骤3,基于上述步骤2获得的最小蠕变速率数据,在相同温度下,利用最小蠕变速率/>施加的应力σ以及阈值应力σth之间的关系公式:/>将上式方程两边同乘以1/5次幂,得到:/>对同一温度下获得的试验数据/>进行线性拟合,拟合直线与X轴截距即为该温度下对应的阈值应力;对不同的温度条件重复上述过程,即可得到不同温度对应的阈值应力。
步骤4,基于步骤3获得的不同温度下的阈值应力水平σth,利用多项式形式进行拟合,从而建立阈值应力与温度之间的函数关系,即:式中,m为多项式项数,ai为拟合参数,i=0,1,2···,m,一般可取m≤3。
步骤5,通过上述步骤得到了阈值应力、抗拉强度之后,建立蠕变寿命预测公式:
式中,A、n为常数,σth为阈值应力,σb为抗拉强度,σ为施加的应力,T为施加的温度,单位为开尔文K,R为气体常数(R=8.314J/(mol·K)),为表观激活能,可以在相同的值下由试验最小蠕变速率的对数与温度的倒数之间的关系确定。具体为当/>值相同时,/>与1/T数据的线性拟合直线的斜率即为/>进而可以获得表观激活能/>数值。
对蠕变寿命预测公式两边进行相应数学变换,并同时取对数,得到
A与n可通过试验数据的最佳线性拟合直线的斜率与截距结果获得。
因此,可以通过步骤5建立的蠕变寿命预测公式,以对任意不同应力温度条件下的蠕变寿命进行准确的预测。保证当应力趋向于阈值应力σth时,蠕变寿命趋近于无穷大,当应力趋向于抗拉强度σb时,蠕变寿命趋近于0。
下面,结合具体实施例对本发明做进一步说明,根据实施例可以更好地理解本发明。实施例所描述的具体物料配比、工艺条件及其结果仅用于说明本发明,而不应当也不会限制权利要求书中所详细描述的本发明。
实施例
在本实施例中,将本发明的蠕变寿命预测方法应用到镍基高温合金GH4169材料的蠕变寿命预测中,包括以下步骤:
(1)在600℃和650℃下进行GH4169材料高温拉伸试验,得到对应的抗拉强度分别为1440MPa和1255MPa。
(2)在600℃和650℃下分别开展6个不同应力值下的GH4169材料高温蠕变试验,得到对应的最小蠕变速率以及蠕变寿命tf。具体试验方案与获得的试验数据如表1所示。
表1 GH4169材料蠕变试验方案与数据
(3)利用公式分别在600℃和650℃下线性拟合/>数据,拟合直线与X轴交点所对应的应力值即为该温度下的阈值应力。所计算的阈值应力如图1所示,进而得到600℃下的阈值应力为593MPa,650℃下的阈值应力为309MPa。利用该两个温度下的阈值应力,其他温度下的阈值应力可通过线性内插或者外推的方法计算得到。
(4)基于上述获得的600℃和650℃下的阈值应力水平σth,可利用多项式形式进行拟合,从而建立阈值应力与温度之间的函数关系。由于试验仅开展了两个温度,因此采用了线性拟合方式,即取多项式形式的前两项。得到阈值应力与温度之间的函数关系式为:σth=-5.68*T+5551.64,其中T为开尔文温度。
(5)在相同的值下,利用函数/>线性拟合试验数据,斜率即为/>进而求得/>如图2所示。随后,线性拟合/>试验数据,相应拟合直线斜率与截距即为未知参数n与lnA值。拟合直线如图3所示,确定系数为0.9296。拟合得到n=-2.3563,A=4.2769。于是,得到蠕变寿命预测公式为:
因此,对于任意温度应力下的蠕变寿命预测,可通过试验得到的抗拉强度与阈值应力,利用多项式拟合方法得到相应温度下的抗拉强度以及阈值应力后,将温度应力条件代入上式即可获得对应的蠕变寿命预测结果。在600℃和650℃下所预测的蠕变寿命与应力之间的关系如图4所示。从图中可以看出,蠕变寿命预测结果与试验结果吻合程度较高,展现出了令人满意的预测精度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.一种基于阈值应力和抗拉强度的蠕变寿命预测方法,其特征在于:包括以下步骤:
步骤1,进行不同温度T下的材料高温拉伸试验,得到对应温度下的抗拉强度σb
步骤2,在不同温度下分别开展多组不同应力条件下的高温蠕变试验,获得对应温度应力条件下的最小蠕变速率以及蠕变寿命tf
步骤3,根据步骤2获得的最小蠕变速率得到不同温度对应的阈值应力σth
所述步骤3中,根据步骤2中高温蠕变试验获得的最小蠕变速率利用公式建立同一温度下最小蠕变速率/>应力σ与阈值应力σth之间的关系,从而求解该温度下对应的阈值应力σth,式中,Am为常数;对不同温度进行同样操作,进而得到不同温度对应的阈值应力水平;
步骤4,基于步骤3获得的不同温度下的阈值应力σth,建立阈值应力与温度之间的函数关系;
所述步骤4中,基于步骤3获得的不同温度下的阈值应力水平σth,利用多项式形式进行拟合,从而建立阈值应力与温度之间的函数关系,即:式中,m为多项式项数,ai为拟合参数,i=0,1,2···,m,m≤3;
步骤5,在步骤4得到的阈值应力、步骤1得到的抗拉强度基础上,建立蠕变寿命预测公式;通过蠕变寿命预测公式,对不同应力温度条件下的蠕变寿命进行预测;
所述步骤5中,建立的蠕变寿命预测公式为:
式中,A、n为常数,σth为阈值应力,σb为抗拉强度,σ为施加的应力,T为施加的温度,R为气体常数,为表观激活能。
2.根据权利要求1所述的基于阈值应力和抗拉强度的蠕变寿命预测方法,其特征在于:所述步骤5中,表观激活能由以下方法得到:在同一/>值下,由试验蠕变速率的对数与温度的倒数1/T之间的线性拟合直线斜率确定表观激活能/>
3.根据权利要求1所述的基于阈值应力和抗拉强度的蠕变寿命预测方法,其特征在于:所述步骤5中,常数A、n由以下方法得到:对蠕变寿命预测公式进行转换,两边取对数,有:
因此,常数A、n通过的线性拟合直线的斜率与截距获得。
4.根据权利要求1所述的基于阈值应力和抗拉强度的蠕变寿命预测方法,其特征在于:所述步骤5中,通过蠕变寿命预测公式,对不同应力温度条件下的蠕变寿命进行预测时,应当保证当应力趋向于阈值应力σth时,蠕变寿命趋近于无穷大,当应力趋向于抗拉强度σb时,蠕变寿命趋近于0。
CN202111040892.4A 2021-09-07 2021-09-07 一种基于阈值应力和抗拉强度的蠕变寿命预测方法 Active CN113866008B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111040892.4A CN113866008B (zh) 2021-09-07 2021-09-07 一种基于阈值应力和抗拉强度的蠕变寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111040892.4A CN113866008B (zh) 2021-09-07 2021-09-07 一种基于阈值应力和抗拉强度的蠕变寿命预测方法

Publications (2)

Publication Number Publication Date
CN113866008A CN113866008A (zh) 2021-12-31
CN113866008B true CN113866008B (zh) 2024-06-11

Family

ID=78989762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111040892.4A Active CN113866008B (zh) 2021-09-07 2021-09-07 一种基于阈值应力和抗拉强度的蠕变寿命预测方法

Country Status (1)

Country Link
CN (1) CN113866008B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1477383A (zh) * 2003-06-23 2004-02-25 国电热工研究院 高温部件蠕变寿命的测试方法
JP2012202908A (ja) * 2011-03-28 2012-10-22 Hiroyuki Sato クリープ曲線およびクリープ寿命の予測方法
CN106568655A (zh) * 2016-10-28 2017-04-19 沈阳工业大学 一种预测耐热合金蠕变寿命的方法
CN110411864A (zh) * 2018-04-26 2019-11-05 天津大学 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法
CN110688788A (zh) * 2019-08-28 2020-01-14 南京航空航天大学 一种高温材料蠕变变形和寿命预测方法及模型
CN112730061A (zh) * 2020-12-24 2021-04-30 南京航空航天大学 一种多级变温变载蠕变寿命评价方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1477383A (zh) * 2003-06-23 2004-02-25 国电热工研究院 高温部件蠕变寿命的测试方法
JP2012202908A (ja) * 2011-03-28 2012-10-22 Hiroyuki Sato クリープ曲線およびクリープ寿命の予測方法
CN106568655A (zh) * 2016-10-28 2017-04-19 沈阳工业大学 一种预测耐热合金蠕变寿命的方法
CN110411864A (zh) * 2018-04-26 2019-11-05 天津大学 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法
CN110688788A (zh) * 2019-08-28 2020-01-14 南京航空航天大学 一种高温材料蠕变变形和寿命预测方法及模型
CN112730061A (zh) * 2020-12-24 2021-04-30 南京航空航天大学 一种多级变温变载蠕变寿命评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"A novel damage constitutive model for creep deformation and damage evolution prediction";Xu Zhao等;《Fatigue Fract Eng Mater Struct》;20221103;全文 *
"Improvement of Creep Behavior Prediction Using Threshold Stress and Tensile Properties: Introduction of the TTC Relations";XU ZHAO等;《METALLURGICAL AND MATERIALS TRANSACTIONS A》;20220724;全文 *
岩石蠕变及疲劳损伤特性的研究进展;王亚松;马林建;刘新宇;陈艳;;工业建筑;20160420(第04期);全文 *

Also Published As

Publication number Publication date
CN113866008A (zh) 2021-12-31

Similar Documents

Publication Publication Date Title
CN108170905B (zh) 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法
CN109142083B (zh) 一种变载荷历程下的蠕变损伤计算方法
CN114295491B (zh) 一种蠕变损伤及变形随时间演化行为的预测方法
CN112730061B (zh) 一种多级变温变载蠕变寿命评价方法
CN110220805A (zh) 一种基于蠕变疲劳损伤累积的变幅多轴热机疲劳寿命预测方法
CN113866008B (zh) 一种基于阈值应力和抗拉强度的蠕变寿命预测方法
CN112417580B (zh) 一种考虑保载效应的高温低周疲劳寿命预测方法
CN115415542A (zh) 基于响应面法的双相不锈钢3d打印件性能的预测方法
CN108460183A (zh) 一种小样本量测定航空发动机材料高周疲劳p-s-n曲线测试方法
Ishihara et al. Acceleration of Grain Boundary Cracking in Ni-Base Alloy 617 Under Creep-Fatigue Loading at 800° C
CN113008677B (zh) 一种镍基高温合金的蠕变持久性能预测方法
CN110940572A (zh) 一种高Cr铁素体耐热钢的蠕变寿命预测方法
CN112329270A (zh) 一种基于加速因子模型的步降应力加速试验方法
Kikuchi et al. Molecular Dynamics Analysis of the Acceleration of Intergranular Cracking of Ni-Base Superalloy Caused by Accumulation of Vacancies and Dislocations Around Grain Boundaries
CN112214875B (zh) 一种通过沉淀颗粒相尺寸评估工件真实服役温度的方法
CN114566236A (zh) 含涂层高温合金高温持久/蠕变寿命的预测方法及系统
Frangini Corrosion rate and anodic dissolution behavior of a B2-iron aluminide alloy in sulfuric acid
CN113702204B (zh) 一种高参数电站用先进高温合金持久强度的预测方法及系统
CN114936474B (zh) 一种通过迟滞回线宽度预测纤维增强陶瓷基复合材料界面剪应力的方法
Gray et al. An investigation of creep under non-constant stress and temperature conditions in Waspaloy
CN102383024A (zh) 一种镍铁基合金和形成镍铁基合金的方法
LU501526B1 (en) A small sample data life prediction method
CN114329952A (zh) 一种高温合金的长程蠕变寿命预测方法
Gishvarov et al. Methods of accelerated tests of turbine blades in power plants operated under multicomponent loading conditions
Haque Modification of the MPC omega model to predict primary and tertiary creep

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