CN105628511A - 一种高温合金蠕变寿命预测方法 - Google Patents

一种高温合金蠕变寿命预测方法 Download PDF

Info

Publication number
CN105628511A
CN105628511A CN201510997412.1A CN201510997412A CN105628511A CN 105628511 A CN105628511 A CN 105628511A CN 201510997412 A CN201510997412 A CN 201510997412A CN 105628511 A CN105628511 A CN 105628511A
Authority
CN
China
Prior art keywords
creep
curve
exp
formula
stress
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.)
Granted
Application number
CN201510997412.1A
Other languages
English (en)
Other versions
CN105628511B (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201510997412.1A priority Critical patent/CN105628511B/zh
Publication of CN105628511A publication Critical patent/CN105628511A/zh
Application granted granted Critical
Publication of CN105628511B publication Critical patent/CN105628511B/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

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

本发明提供了一种高温合金的蠕变寿命预测方法,该方法是在θ投影法的基础上,引入蠕变过程中材料形状变化导致真应力变化这一因素,得到修正θ投影法,具体形式为:ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)。该方法可用于利用蠕变曲线的第一阶段至第三阶段初期数据预测整个均匀变形阶段的蠕变曲线,有利于缩短获得蠕变数据的时间。同时还可用于预测其他蠕变条件下的蠕变曲线,包括预测蠕变某一蠕变条件下的蠕变断裂时间和某一蠕变应变对应的蠕变中断时间。该方法简单可靠,适用于不同类型的高温合金,适合工程应用,在高温合金蠕变寿命管理方面具有广阔前景。

Description

一种高温合金蠕变寿命预测方法
技术领域
本发明属于材料科学与工程技术应用,具体涉及高温合金蠕变寿命预测方法。
背景技术
自20世纪30年代开始,世界各国开始高温合金的研究,并在接下来的半个多世纪中迅猛发展。高温合金以其优异的高温强度,良好的抗氧化和抗热腐蚀性能,良好的疲劳性能、断裂韧性等综合性能,被广泛应用于制造航空、舰艇和工业用燃气轮机的涡轮叶片、导向叶片、涡轮盘、高压压气机盘和燃烧室等高温部件,已成为军民用燃气涡轮发动机热端部件不可替代的关键材料。
高温合金在服役过程中通常会因为蠕变变形导致失效,从而造成巨大的损失。因此若能够准确地预测高温合金蠕变时间,尤其是某一蠕变应变条件对应的蠕变中断时间,则可对高温合金服役时间进行有效的管理,在避免造成合金蠕变失效的情况下最大程度地对其进行利用。
从20世纪50年代开始,各种蠕变寿命预测模型被提出。1982年R.M.Evans等人提出通过对蠕变曲线进行描述,从而预测蠕变寿命的方法,称为“θ投影法”。θ投影法认为蠕变过程可分为材料软化过程和材料硬化过程,且将材料软化过程和材料硬化过程均以指数方程的形式表示,二者叠加可得θ投影法方程。由于其不仅能对不同条件下的蠕变断裂时间进行预测,同时还可预测某一蠕变应变所对应的蠕变中断时间,θ投影法在提出后被广泛应用于各种类型的合金。但是,θ投影法的提出基于理想蠕变条件,即温度和应力均为恒定;而在真实蠕变过程中,蠕变工程应力恒定,随着应变量的增加,蠕变样品的横截面逐渐缩小,从而导致蠕变真应力逐渐增加,二者互相矛盾。由于蠕变过程中形状变化导致的应力变化不可避免,且会导致蠕变寿命低于恒定应力下的蠕变寿命,从而使得利用模拟蠕变中断曲线所得的θ投影法方程预测同一条件下更高应变对应的蠕变时间时结果偏大;更有甚者,对于蠕变应变较大的蠕变曲线,其曲线不符合θ投影法的规律。因此需要针对此现象对θ投影法进行修正,以求得到更高精度的蠕变寿命预测模型。
发明内容
本发明的目的是提出一种可准确描述蠕变曲线的蠕变寿命预测方法,该方法通过引入蠕变实验中由于蠕变样品形状变化导致的应力变化这一因素,对θ投影法进行修正,减小预测的蠕变曲线与实验所得的蠕变曲线之间的差别,该方法简单可靠,适用于不同类型的高温合金。
一种高温合金的蠕变寿命预测方法,其核心要点在于:引入实际蠕变过程中,蠕变真应力随应变增加而增加,从而导致蠕变寿命降低这一因素,即考虑蠕变过程中形状变化这一损伤参量。本发明的实现步骤如下:
(1)获取至少5条不同温度和应力条件下合金的高温蠕变曲线,每条蠕变曲线包括温度、应力、蠕变应变和蠕变时间。
(2)将蠕变曲线第一阶段至第三阶段初期的数据按公式
ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)进行拟合,求出公式中的参数θ1、θ2、θ3、θ4、θ5,得到修正后的θ投影法方程。
(3)利用(2)中结果可通过延长θ投影法曲线,获得均匀变形阶段的蠕变曲线。
(4)将不同温度和应力条件下的蠕变曲线拟合得到的θ1、θ2、θ3、θ4、θ5按公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ分别计算参数a、b、c、d、e、f。
(5)利用(2)和(4)所得结果可预测不同温度和应力下的蠕变曲线,从而获得对应温度和应力下的蠕变断裂时间或某一应变所对应的蠕变中断时间。
所述公式ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)中,ε表示蠕变应变,t表示蠕变时间。
所述公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ中,T表示蠕变温度,σ表示蠕变工程应力,a、b、c、d、e、f为与材料相关的常数。
所述参数θ1、θ2、θ3、θ4、θ5是利用数学分析软件,按最小二乘法回归,将蠕变曲线时间-应变数值输入,按公式ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)求得待定系数θ1、θ2、θ3、θ4、θ5
所述参数a、b、c、d、e、f是利用数学分析软件,按最小二乘法回归,将蠕变温度、蠕变应力及对应的θ1、θ2、θ3、θ4、θ5输入,按公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ求得待定系数a、b、c、d、e、f。
本发明的优点在于:能够准确地对不同类型的高温合金蠕变曲线进行描述,能够提高高温合金蠕变断裂时间和某一应变量对应的蠕变中断时间的预测精度。该方法简单可靠,利用常规蠕变试验机获得蠕变数据。利用此方法可对高温合金的服役时间进行有效估计,减小危害,降低成本,适用于不同类型的高温合金。
附图说明
图1为K465合金不同蠕变条件下的蠕变时间-应变曲线。
图2为K465合金蠕变实验曲线与修正后的θ投影法预测曲线对比图。
图3为K465合金蠕变实验曲线与经典θ投影法预测曲线对比图(a)900℃/300MPa,(b)900℃/320MPa,(c)950℃/300MPa,(d)975℃/200MPa,(e)1000℃/137MPa。
图4为修正后的θ投影法预测不同条件下的K465合金蠕变曲线。
具体实施方式
以下实例将对本发明予以进一步的说明,以便本领域人员更好地理解本发明的优点和特征。
首先,以K465合金为研究对象,利用高温蠕变试验机,获得其在五个不同温度和应力条件下的蠕变曲线,实验条件分别为1000℃/137MPa、975℃/200MPa、950℃/300MPa、900℃/320MPa、900℃/300MPa,实验时每隔30min取一个数据点,实验结果见图1。
其次,分别利用公式ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)及经典θ投影法方程对蠕变曲线ε<1.3%的蠕变数据进行拟合,获得待定系数θ1、θ2、θ3、θ4、θ5(见表1),得到修正后的θ投影法方程。具体拟合方法为利用Originlab软件的非线性拟合工具,将方程写入并选择合适的初始值,选择需要拟合的数据进行拟合。
图2为利用修正后的θ投影法方程作出的完整蠕变曲线与实验所得蠕变曲线的对比。图3为经典θ投影法方程作出的完整蠕变曲线与实验所得蠕变曲线的对比。由图2和图3可以看出,修正后的θ投影法大大提高了预测精度。利用不同蠕变条件下的系数θ1、θ2、θ3、θ4、θ5及公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ拟合得出待定系数a、b、c、d、e、f,见表2。具体拟合过程由Originlab软件的线性拟合工具实现。
利用上述结果,通过计算不同温度和应力下的系数θ1、θ2、θ3、θ4、θ5,可得到对应的修正后的θ投影法方程,从而达到预测不同条件蠕变曲线的目的,见图4。
表1修正后的θ投影法拟合K465合金不同蠕变条件下的系数
900℃/300MPa 900℃/320MPa 950℃/300MPa 975℃/200MPa 1000℃/137MPa
θ1 0.09795 0.15025 0.10966 0.18653 0.22057
θ2 0.05956 0.06252 0.89837 0.11496 0.07638
θ3 1.19714 1.47450 2.72696 1.43088 0.59260
θ4 0.00287 0.00380 0.01257 0.00462 0.00402
θ5 0.27149 0.26709 0.26504 0.31151 0.33348
表2K465合金各θ值对应的蠕变材料常数
a b c d e f
θ1 -7.768 0.02253 0.007390 -2.431E-5 -- --
θ2 -22.68 -0.003030 0.01959 1.712E-5 -- --
θ3 7.069 -0.04362 -0.008130 4.980E-5 -- --
θ4 -20.08 0.01924 0.01671 -1.221E-5 -- --
θ5 -- -- -- -- -0.3931 -5.908E-4

Claims (5)

1.一种高温合金的蠕变寿命预测方法,其特征在于,具体包括如下步骤:
1)获取至少5条不同温度和应力条件下合金的高温蠕变曲线,每条蠕变曲线包括温度、应力、蠕变应变和蠕变时间;
2)将蠕变曲线第一阶段至第三阶段初期的数据按公式
ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)进行拟合,求出公式中的待定系数θ1、θ2、θ3、θ4、θ5,得到修正后的θ投影法方程;
3)利用2)中结果通过延长θ投影法曲线,获得均匀变形阶段的蠕变曲线;
4)将不同温度和应力条件下的蠕变曲线拟合得到的θ1、θ2、θ3、θ4、θ5按公式logθi=a+b*σ+c*T+d*σ*T,i=1,2,3,4和公式logθ5=e+f*σ分别计算参数a、b、c、d、e、f;
5)利用2)和4)所得结果预测不同温度和应力下的蠕变曲线,从而获得对应温度和应力下的蠕变断裂时间或蠕变中断时间。
2.按照权利要求1所述的一种高温合金蠕变寿命预测方法,其特征在于:公式ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)中,ε表示蠕变应变,t表示蠕变时间。
3.按照权利要求1所述的一种高温合金蠕变寿命预测方法,其特征在于:公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ中,T表示蠕变温度,σ表示蠕变工程应力,a、b、c、d、e、f为与材料相关的常数。
4.按照权利要求1所述的一种高温合金蠕变寿命预测方法,其特征在于:按最小二乘法回归,将蠕变曲线时间-应变数值输入,按公式ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)求得待定系数θ1、θ2、θ3、θ4、θ5
5.按照权利要求1所述的一种高温合金蠕变寿命预测方法,其特征在于:按最小二乘法回归,将蠕变温度、蠕变应力及对应的θ1、θ2、θ3、θ4、θ5输入,按公式logθi=a+b*σ+c*T+d*σ*T(i=1,2,3,4)和公式logθ5=e+f*σ求得待定系数a、b、c、d、e、f。
CN201510997412.1A 2015-12-25 2015-12-25 一种高温合金蠕变寿命预测方法 Active CN105628511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510997412.1A CN105628511B (zh) 2015-12-25 2015-12-25 一种高温合金蠕变寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510997412.1A CN105628511B (zh) 2015-12-25 2015-12-25 一种高温合金蠕变寿命预测方法

Publications (2)

Publication Number Publication Date
CN105628511A true CN105628511A (zh) 2016-06-01
CN105628511B CN105628511B (zh) 2018-09-25

Family

ID=56043650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510997412.1A Active CN105628511B (zh) 2015-12-25 2015-12-25 一种高温合金蠕变寿命预测方法

Country Status (1)

Country Link
CN (1) CN105628511B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106568655A (zh) * 2016-10-28 2017-04-19 沈阳工业大学 一种预测耐热合金蠕变寿命的方法
CN108170905A (zh) * 2017-12-08 2018-06-15 南昌航空大学 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法
CN108931448A (zh) * 2018-05-07 2018-12-04 华南理工大学 一种高铬钢材料热力学响应及疲劳-蠕变损伤的预测方法
CN109142362A (zh) * 2018-09-10 2019-01-04 中国大唐集团科学技术研究院有限公司华中分公司 一种基于bp神经网络的电厂用金属材料寿命评估方法
CN109253930A (zh) * 2017-07-12 2019-01-22 天津大学 高温管道中含椭圆形轴向外表面裂纹的孕育期预测方法
CN109932241A (zh) * 2017-12-15 2019-06-25 天津大学 塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法
CN110411851A (zh) * 2019-07-22 2019-11-05 北京科技大学 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法
CN110411850A (zh) * 2019-07-22 2019-11-05 北京科技大学 一种高温合金涡轮叶片服役条件的评估方法
CN110967245A (zh) * 2018-09-28 2020-04-07 中国航发商用航空发动机有限责任公司 材料蠕变时间及寿命实验方法以及实验系统
CN113125275A (zh) * 2021-04-06 2021-07-16 西北工业大学 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法
CN113252465A (zh) * 2021-05-20 2021-08-13 天津理工大学 一种基于m-h法的耐热钢蠕变寿命预测方法
CN113420502A (zh) * 2021-06-11 2021-09-21 北京科技大学 一种含Re、Ru镍基单晶高温合金蠕变寿命预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710053A (zh) * 2009-11-06 2010-05-19 上海师范大学 高温材料的蠕变寿命预测方法
CN103267683A (zh) * 2013-04-28 2013-08-28 扬州大学 一种确定耐热金属材料剩余使用寿命的方法
CN103712865A (zh) * 2012-10-01 2014-04-09 株式会社日立制作所 高温部件的蠕变损伤评价方法以及损伤评价系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710053A (zh) * 2009-11-06 2010-05-19 上海师范大学 高温材料的蠕变寿命预测方法
CN103712865A (zh) * 2012-10-01 2014-04-09 株式会社日立制作所 高温部件的蠕变损伤评价方法以及损伤评价系统
CN103267683A (zh) * 2013-04-28 2013-08-28 扬州大学 一种确定耐热金属材料剩余使用寿命的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
江冯: "采用不同方法预测10Cr9MoW2VNbBN钢蠕变曲线的比较", 《机械工程材料》 *
耿鲁阳: "基于修正θ投影法的炉管用Cr25Ni35Nb钢蠕变寿命预测", 《机械工程材料》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106568655A (zh) * 2016-10-28 2017-04-19 沈阳工业大学 一种预测耐热合金蠕变寿命的方法
CN106568655B (zh) * 2016-10-28 2019-04-12 沈阳工业大学 一种预测耐热合金蠕变寿命的方法
CN109253930B (zh) * 2017-07-12 2020-09-08 天津大学 高温管道中含椭圆形轴向外表面裂纹的孕育期预测方法
CN109253930A (zh) * 2017-07-12 2019-01-22 天津大学 高温管道中含椭圆形轴向外表面裂纹的孕育期预测方法
CN108170905B (zh) * 2017-12-08 2020-02-04 南昌航空大学 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法
CN108170905A (zh) * 2017-12-08 2018-06-15 南昌航空大学 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法
CN109932241A (zh) * 2017-12-15 2019-06-25 天津大学 塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法
CN108931448B (zh) * 2018-05-07 2021-08-10 华南理工大学 一种高铬钢材料热力学响应及疲劳-蠕变损伤的预测方法
CN108931448A (zh) * 2018-05-07 2018-12-04 华南理工大学 一种高铬钢材料热力学响应及疲劳-蠕变损伤的预测方法
CN109142362A (zh) * 2018-09-10 2019-01-04 中国大唐集团科学技术研究院有限公司华中分公司 一种基于bp神经网络的电厂用金属材料寿命评估方法
CN109142362B (zh) * 2018-09-10 2021-01-26 中国大唐集团科学技术研究院有限公司华中分公司 一种基于bp神经网络的电厂用金属材料寿命评估方法
CN110967245A (zh) * 2018-09-28 2020-04-07 中国航发商用航空发动机有限责任公司 材料蠕变时间及寿命实验方法以及实验系统
CN110411851A (zh) * 2019-07-22 2019-11-05 北京科技大学 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法
CN110411851B (zh) * 2019-07-22 2020-06-26 北京科技大学 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法
CN110411850A (zh) * 2019-07-22 2019-11-05 北京科技大学 一种高温合金涡轮叶片服役条件的评估方法
CN113125275A (zh) * 2021-04-06 2021-07-16 西北工业大学 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法
CN113125275B (zh) * 2021-04-06 2024-05-03 西北工业大学 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法
CN113252465A (zh) * 2021-05-20 2021-08-13 天津理工大学 一种基于m-h法的耐热钢蠕变寿命预测方法
CN113252465B (zh) * 2021-05-20 2022-07-01 天津理工大学 一种基于m-h法的耐热钢蠕变寿命预测方法
CN113420502A (zh) * 2021-06-11 2021-09-21 北京科技大学 一种含Re、Ru镍基单晶高温合金蠕变寿命预测方法
CN113420502B (zh) * 2021-06-11 2023-08-11 北京科技大学 一种含Re、Ru镍基单晶高温合金蠕变寿命预测方法

Also Published As

Publication number Publication date
CN105628511B (zh) 2018-09-25

Similar Documents

Publication Publication Date Title
CN105628511A (zh) 一种高温合金蠕变寿命预测方法
Gallo et al. Generalized approach to estimation of strains and stresses at blunt V‐notches under non‐localized creep
CN103267683B (zh) 一种确定耐热金属材料剩余使用寿命的方法
CN108170905A (zh) 一种用于镍基高温合金叶片热机械疲劳载荷下的寿命预测方法
CN103886125B (zh) 一种钛合金热复合成形数值模拟方法
Campagnolo et al. Rapid finite element evaluation of the averaged strain energy density of mixed‐mode (I+ II) crack tip fields including the T‐stress contribution
JP2007256042A (ja) ガスタービン高温部品のき裂進展予測方法及びこの方法を用いたき裂進展予測装置
CN111079276A (zh) 多失效模式下电子产品可靠性指标加速验证方法
Vakili-Tahami et al. Investigating the possibility of replacing IN 738LC gas turbine blades with IN 718
CN108875266A (zh) 一种涡轮叶片区间损伤容限分析方法
Bhachu et al. Application of 3D fracture mechanics for improved crack growth predictions of gas turbine components
Zhang et al. Creep‐fatigue behavior of thin‐walled plate with holes: stress state characterization and life estimation
Zhou et al. Model and application of product conflict problem with integrated TRIZ and Extenics for low-carbon design
Zhou et al. An evaluation method based on mechanical parts structural characteristics for proactive remanufacturing
AU2022201697B2 (en) Method and system of multi-objective and multi-dimensional online joint monitoring for nuclear turbine
CN115979804A (zh) 一种增材制造材料寿命预测方法
CN114462788B (zh) 一种基于半马尔科夫过程的多状态元件可靠性评估方法
Moneta et al. Influence of manufacturing tolerances on vibration frequencies of turbine blade
CN115455594A (zh) 基于涡轮叶片冷却设计数据库的内流通道数据修正方法
CN117057082A (zh) 维修极限尺寸设计方法、维修检查方法和存储介质
CN115204013A (zh) 一种多轴应力状态下的材料的寿命预测方法
CN114139276A (zh) 一种盘轴一体式整体叶盘结构疲劳寿命分析方法
CN113656970A (zh) 一种航空发动机产品工艺提升方法、系统、设备及介质
CN109959555B (zh) 考虑厚度和裂纹深度影响的中心裂纹拉伸试样的孕育期预测方法
Chen et al. Numerical simulation and experiment study on the nuclear fuel spacer grid stamping of Inconel 718

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant