CN105651217A - 一种大体积钢中非金属夹杂物尺寸的统计计算方法 - Google Patents

一种大体积钢中非金属夹杂物尺寸的统计计算方法 Download PDF

Info

Publication number
CN105651217A
CN105651217A CN201511025126.5A CN201511025126A CN105651217A CN 105651217 A CN105651217 A CN 105651217A CN 201511025126 A CN201511025126 A CN 201511025126A CN 105651217 A CN105651217 A CN 105651217A
Authority
CN
China
Prior art keywords
inclusions
size
steel
threshold value
distribution function
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
CN201511025126.5A
Other languages
English (en)
Other versions
CN105651217B (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.)
China National Petroleum Corp
CNPC Tubular Goods Research Institute
Original Assignee
China National Petroleum Corp
CNPC Tubular Goods Research Institute
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 China National Petroleum Corp, CNPC Tubular Goods Research Institute filed Critical China National Petroleum Corp
Priority to CN201511025126.5A priority Critical patent/CN105651217B/zh
Publication of CN105651217A publication Critical patent/CN105651217A/zh
Application granted granted Critical
Publication of CN105651217B publication Critical patent/CN105651217B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Abstract

本发明公开了一种大体积钢中非金属夹杂物尺寸的统计计算方法,包括以下步骤:1)采用金相检测方法测量若干样品金相试样中夹杂物的尺寸,建立非金属夹杂物统计模型,所述非金属夹杂物统计模型采用帕累托分布函数的形式,确定夹杂物尺寸的门限值,再根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值,然后将帕累托分布函数中的参数值代入非金属夹杂物统计模型中,得非金属夹杂物统计模型的表达式;2)确定待计算钢的体积,再将待计算钢的体积带入到非金属夹杂物统计模型的表达式中,得待计算钢中非金属夹杂物的尺寸信息。本发明能够准确的获取大体积钢中非金属夹杂物尺寸信息。

Description

一种大体积钢中非金属夹杂物尺寸的统计计算方法
技术领域
本发明属于大体积钢中非金属夹杂物水平的测量领域,涉及一种大体积钢中非金属夹杂物尺寸的统计计算方法。
背景技术
炼钢过程中,少量炉渣、耐火材料及钢水冶炼中的反应产物会进入钢液,在钢中形成非金属夹杂物。非金属夹杂物会降低钢的机械性能,特别是降低塑性、韧性及疲劳极限。严重时,非金属夹杂物会使钢在热加工与热处理时产生裂纹或发生突然脆断。
非金属夹杂物是钢中必然存在的产物,小尺寸非金属夹杂物对钢的性能是无害的,并且是大量存在的,而大尺寸夹杂物会严重恶化钢的机械性能,导致钢材突然断裂,发生灾难性后果。控制钢中夹杂物尺寸是钢铁冶炼和工程应用的重要目标之一。由于非金属夹杂物的尺寸基本在微米级,需要借助分析仪器才能测量。目前通用的非金属夹杂物的测量方法为金相法,金相法是首先利用机械研磨和抛光的方法进行样品制备,然后在金相显微镜上进行测量分析,根据夹杂物分类方法进行夹杂物评级。金相法进行非金属夹杂物评级时,由于采用小样品进行夹杂物观察,所以存在很大的偶然性;另外,金相小样品夹杂物检测评级只能反应取样样品的夹杂物水平,而不能反应大体积钢中夹杂物水平,特别是一些大体积钢结构件(如大锻件),并且现在的检测水平又无仪器装备法直接检测这类大锻件中夹杂物情况。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种大体积钢中非金属夹杂物尺寸的统计计算方法,该方法能够准确的获取大体积钢中非金属夹杂物尺寸信息。
为达到上述目的,本发明所述的大体积钢中非金属夹杂物尺寸的统计计算方法包括以下步骤:
1)采用金相检测方法测量若干样品金相试样中夹杂物的尺寸,建立非金属夹杂物统计模型,所述非金属夹杂物统计模型采用帕累托分布函数的形式,确定夹杂物尺寸的门限值,再根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值,然后将帕累托分布函数中的参数值代入非金属夹杂物统计模型中,得非金属夹杂物统计模型的表达式;
2)确定待计算钢的体积,再将待计算钢的体积带入到非金属夹杂物统计模型的表达式中,得待计算钢中非金属夹杂物的尺寸信息,其中,各样品金相试样的体积小于待计算钢的体积。
根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用极大似然法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用作图法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
所述渐近帕累托分布函数中参数值包括非金属夹杂物统计模型中夹杂物门槛值、分布函数的比例参数及夹杂物的形状参数。
样品金相试样的数量为5-10个。
设夹杂物尺寸的门槛值u,则有
F(x)=1-(1+ξ(x-u)/λ)-1/ξ(1)
其中,λ为比例参数,λ>0,ξ(-∞<ξ<∞)为形状参数,x为夹杂物的尺寸,F(x)为超过门槛值u的夹杂物数量,x-u为超过门槛值u的夹杂物,其中
0 < x - u < &infin; ...... i f &xi; &GreaterEqual; 0 0 < x - u < &lambda; ... i f &xi; < 0
随着减小,F(x)逐渐增大,则有
F(x)=1-exp(-(x-u)/λ)(2)。
设NV(u)为单位体积内尺寸超过门槛值u的夹杂物数量的期望值,xV
为V中估算的夹杂物尺寸,则有
NV(u)V(1-F(xV))=1(3)
由式(1)及式(3)得:
X v = u - &lambda; &xi; &lsqb; 1 - ( N v ( u ) V ) &xi; &rsqb; - - - ( 4 )
式中,NV(u)根据金相检测抛光横截面中观测到的夹杂数量运用Woodhead分析方法计算得到。
本发明具有以下有益效果:
本发明所述的大体积钢中非金属夹杂物尺寸的统计计算方法在计算过程中,通过传统的金相检测方法测量若干样品金相试样中夹杂物的尺寸,建立非金属夹杂物统计模型,非金属夹杂物统计模型采用帕累托分布函数形式,再确定帕累托分布函数中的参数值,然后将待计算钢的体积代入到非金属夹杂物统计模型,即可直接得到待计算钢中非金属夹杂物的尺寸信息,其中,待计算钢的体积大于样品金相试样的体积,本发明解决大体积钢或钢结构件中夹杂物水平的精确计算问题,为大体积钢件的冶炼和工程应用提供了技术参考,并且该方法便于操作和实施。
附图说明
图1为本发明的流程图;
图2为本发明中夹杂物的尺寸示意图;
图3为本发明中作图法确定帕累托分布函数中的参数值时的示意图;
图4为本发明中计算待计算钢中非金属夹杂物的尺寸信息时的示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
参考图1,本发明所述的大体积钢中非金属夹杂物尺寸的统计计算方法包括以下步骤:
1)采用金相检测方法测量若干样品金相试样中夹杂物的尺寸,建立非金属夹杂物统计模型,所述非金属夹杂物统计模型采用帕累托分布函数的形式,确定夹杂物尺寸的门限值,再根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值,然后将帕累托分布函数中的参数值代入非金属夹杂物统计模型中,得非金属夹杂物统计模型的表达式;
2)确定待计算钢的体积,再将待计算钢的体积带入到非金属夹杂物统计模型的表达式中,得待计算钢中非金属夹杂物的尺寸信息,其中,各样品金相试样的体积小于待计算钢的体积。
根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用极大似然法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用作图法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
所述渐近帕累托分布函数中参数值包括非金属夹杂物统计模型中夹杂物门槛值、分布函数的比例参数及夹杂物的形状参数。
样品金相试样的数量为5-10个。
设夹杂物尺寸的门槛值u,则有
F(x)=1-(1+ξ(x-u)/λ)-1/ξ(1)
其中,λ为比例参数,λ>0,ξ(-∞<ξ<∞)为形状参数,x为夹杂物的尺寸,F(x)为超过门槛值u的夹杂物数量,x-u为超过门槛值u的夹杂物,其中
0 < x - u < &infin; ...... i f &xi; &GreaterEqual; 0 0 < x - u < &lambda; ... i f &xi; < 0
随着减小,F(x)逐渐增大,则有
F(x)=1-exp(-(x-u)/λ)(2)。
设NV(u)为单位体积内尺寸超过门槛值u的夹杂物数量的期望值,xV为V中估算的夹杂物尺寸,则有
NV(u)V(1-F(xV))=1(3)
由式(1)及式(3)得:
X v = u - &lambda; &xi; &lsqb; 1 - ( N v ( u ) V ) &xi; &rsqb; - - - ( 4 )
式中,NV(u)根据金相检测抛光横截面中观测到的夹杂数量运用Woodhead分析方法计算得到,其中,
N V ( u ) = N A ( u ) / D i &OverBar; - - - ( 5 )
其中,为夹杂物尺寸大小的平均值。
通常门槛值u的选取采用作图法得到,即把样品金相试样检测获得夹杂物尺寸超过门槛值u的余数平均值对门槛值u作图,所作的曲线接近线性分布时,超过某一门槛值就为合理的线性分布的点,则该某一门槛值则为所要求的临界门槛值u,然后对线性曲线进行数据拟合,所得直线的纵轴截距为λ/(1-ξ),斜率为ξ/(1-ξ),从而计算得参数λ和ξ的值。
为了更精确估计参数λ和ξ的值,则采用最大似然函数法进行分析,其中,最大似然函数方程为:
L = &Pi; 1 k { 1 + &xi; ( X i - u ) &lambda; } - ( &xi; ) - 1 - - - ( 6 )
其中,λ和ξ的估计值为使L取最大值时对应的值。
实施例1
本实施例中使用材料为高强度低合金管线钢,具体操作如下所示:
1)试样制备:利用机械加工方法从钢板上沿横截面取样,试样尺寸为20mm×20mm×T,沿横截面进行机械研磨,然后利用抛光机进行表面抛光处理,抛光后用水冲洗,然后用电吹风吹干水迹,本实施例制备了8个样品;
2)金相检测:将步骤1)制备好的样品放在金相显微镜下进行夹杂物检测和测量,由于高强度低合金管线钢采用炉外精炼技术冶炼,夹杂物尺寸非常小,所以在观察时选择在400倍的放大倍数下进行,每个试样随机选取了10个视场,每个视场测量3个较大的夹杂物尺寸,参见图2,本实施例在夹杂物测量时统一按尺寸进行统计,共统计获得240个夹杂物;
3)概率函数参数值确定:采用作图法和最大似然函数法确定概率函数公式(4)中的参数值,作图法是根据检测统计的夹杂物尺寸,作夹杂物尺寸超过门槛值u的余数平均值对u作图,参考图2,所作的曲线接近线性分布时,超过门槛值u就为合理的线性分布的点,然后对线性曲线进行拟合,利用直线的纵轴截距和斜率计算λ和ξ的值,另外,根据式(6)对式(4)中参数进行计算,表1为采用两种方法获得的夹杂物门槛值u、参数λ和ξ的值;
表1
参数确定方法 门槛值(u),m 参数λ 参数ξ
作图法 3.0 1.88 -0.055
最大似然函数法 3.2 1.90 -0.056
4)大体积钢中夹杂物尺寸估算:把表1中确定的参数中代入式(4)中,给定估算夹杂物尺寸的钢体积,表2为获得的大体积钢中的夹杂物尺寸,参考图3,为了对比,把金相观察体累计积内的最大夹杂物尺寸也添加到了图3中,对比说明小体积内估算钢中夹杂物尺寸与检测结果非常接近。
表2
需要说明的是,随着估算钢体积的增加,其夹杂物尺寸并不是成线性增加,而是越来越接近极限值,这与钢中夹杂物的实际情况相吻合,因为当钢中夹杂物尺寸达到一定值时,由于夹杂物本身浮力和钢水流动,使大尺寸夹杂物会上浮到钢水表面而经过扒渣而清除掉。
以上实施例仅用以说明本材料的技术实施方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的原理和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (7)

1.一种大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,包括以下步骤:
1)采用金相检测方法测量若干样品金相试样中夹杂物的尺寸,建立非金属夹杂物统计模型,所述非金属夹杂物统计模型采用帕累托分布函数的形式,确定夹杂物尺寸的门限值,再根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值,然后将帕累托分布函数中的参数值代入非金属夹杂物统计模型中,得非金属夹杂物统计模型的表达式;
2)确定待计算钢的体积,再将待计算钢的体积带入到非金属夹杂物统计模型的表达式中,得待计算钢中非金属夹杂物的尺寸信息,其中,各样品金相试样的体积小于待计算钢的体积。
2.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用极大似然法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
3.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,根据夹杂物尺寸的门限值及各样品金相试样中夹杂物的尺寸得帕累托分布函数中的参数值的具体操作为:采用作图法根据夹杂物尺寸大于门限值的夹杂物确定渐近帕累托分布函数中参数值。
4.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,所述渐近帕累托分布函数中参数值包括非金属夹杂物统计模型中夹杂物门槛值、分布函数的比例参数及夹杂物的形状参数。
5.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,样品金相试样的数量为5-10个。
6.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,
设夹杂物尺寸的门槛值u,则有
F(x)=1-(1+ξ(x-u)/λ)-1/ξ(1)
其中,λ为比例参数,λ>0,ξ(-∞<ξ<∞)为形状参数,x为夹杂物的尺寸,F(x)为超过门槛值u的夹杂物数量,x-u为超过门槛值u的夹杂物,其中
0 < x - u < &infin; ... ... i f &xi; &GreaterEqual; 0 0 < x - u < &lambda; ... i f &xi; < 0
随着减小,F(x)逐渐增大,则有
F(x)=1-exp(-(x-u)/λ)(2)。
7.根据权利要求1所述的大体积钢中非金属夹杂物尺寸的统计计算方法,其特征在于,
设NV(u)为单位体积内尺寸超过门槛值u的夹杂物数量的期望值,xV为V中估算的夹杂物尺寸,则有
NV(u)V(1-F(xV))=1(3)
由式(1)及式(3)得:
X v = u - &lambda; &xi; &lsqb; 1 - ( N v ( u ) V ) &xi; &rsqb; - - - ( 4 )
式中,NV(u)根据金相检测抛光横截面中观测到的夹杂数量运用Woodhead分析方法计算得到。
CN201511025126.5A 2015-12-30 2015-12-30 一种大体积钢中非金属夹杂物尺寸的统计计算方法 Active CN105651217B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511025126.5A CN105651217B (zh) 2015-12-30 2015-12-30 一种大体积钢中非金属夹杂物尺寸的统计计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511025126.5A CN105651217B (zh) 2015-12-30 2015-12-30 一种大体积钢中非金属夹杂物尺寸的统计计算方法

Publications (2)

Publication Number Publication Date
CN105651217A true CN105651217A (zh) 2016-06-08
CN105651217B CN105651217B (zh) 2018-04-03

Family

ID=56490082

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511025126.5A Active CN105651217B (zh) 2015-12-30 2015-12-30 一种大体积钢中非金属夹杂物尺寸的统计计算方法

Country Status (1)

Country Link
CN (1) CN105651217B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110609042A (zh) * 2019-09-26 2019-12-24 江苏省沙钢钢铁研究院有限公司 一种钢中最大尺寸夹杂物的预测方法
CN111879784A (zh) * 2020-08-10 2020-11-03 耐世特凌云驱动系统(芜湖)有限公司 钢材洁净度评价新方法
CN112485118A (zh) * 2020-11-09 2021-03-12 中国科学院金属研究所 一种估算高强钢疲劳失效临界夹杂物尺寸的方法
CN115267119A (zh) * 2022-07-18 2022-11-01 中天钢铁集团有限公司 一种用于大尺寸非金属夹杂物检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001242163A (ja) * 2000-02-28 2001-09-07 Nisshin Steel Co Ltd シャドーマスク用熱延材の清浄度測定方法
JP2003294639A (ja) * 2002-04-03 2003-10-15 Jfe Steel Kk 非金属介在物の自動検査方法及び装置
CN1651905A (zh) * 2005-02-05 2005-08-10 石家庄钢铁股份有限公司 钢中非金属夹杂物定量分析方法
CN103123329A (zh) * 2012-12-31 2013-05-29 上海大学 金属中非金属夹杂物的快速检测方法及快速检测装置
CN104007245A (zh) * 2014-05-29 2014-08-27 西南交通大学 一种钢铁材料中最大非金属夹杂物尺寸的测试方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001242163A (ja) * 2000-02-28 2001-09-07 Nisshin Steel Co Ltd シャドーマスク用熱延材の清浄度測定方法
JP2003294639A (ja) * 2002-04-03 2003-10-15 Jfe Steel Kk 非金属介在物の自動検査方法及び装置
CN1651905A (zh) * 2005-02-05 2005-08-10 石家庄钢铁股份有限公司 钢中非金属夹杂物定量分析方法
CN103123329A (zh) * 2012-12-31 2013-05-29 上海大学 金属中非金属夹杂物的快速检测方法及快速检测装置
CN104007245A (zh) * 2014-05-29 2014-08-27 西南交通大学 一种钢铁材料中最大非金属夹杂物尺寸的测试方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110609042A (zh) * 2019-09-26 2019-12-24 江苏省沙钢钢铁研究院有限公司 一种钢中最大尺寸夹杂物的预测方法
CN111879784A (zh) * 2020-08-10 2020-11-03 耐世特凌云驱动系统(芜湖)有限公司 钢材洁净度评价新方法
CN112485118A (zh) * 2020-11-09 2021-03-12 中国科学院金属研究所 一种估算高强钢疲劳失效临界夹杂物尺寸的方法
CN112485118B (zh) * 2020-11-09 2022-03-08 中国科学院金属研究所 一种估算高强钢疲劳失效临界夹杂物尺寸的方法
CN115267119A (zh) * 2022-07-18 2022-11-01 中天钢铁集团有限公司 一种用于大尺寸非金属夹杂物检测方法

Also Published As

Publication number Publication date
CN105651217B (zh) 2018-04-03

Similar Documents

Publication Publication Date Title
CN105651217B (zh) 一种大体积钢中非金属夹杂物尺寸的统计计算方法
CN102507400A (zh) 一种t91钢管的剩余寿命定量分析方法
CN103439473B (zh) 一种12Cr1MoV钢受热面状态评估方法
CN107843510A (zh) 基于室温布氏硬度预测超临界机组t/p91耐热钢剩余持久寿命评估方法
CN102692425A (zh) 基于析出相面积分数的t/p91钢老化评级方法
He et al. Effect of constraint on creep crack initiation time in test specimens in ASTM-E1457 standard
Nguyen et al. Change of precipitate free zone during long-term creep in 2.25 Cr–1Mo steel
CN105891093B (zh) 一种铁磁性金属材料抗氢致开裂性能的检测方法
Li et al. Modeling on dynamic recrystallization of aluminium alloy 7050 during hot compression based on cellular automaton
Ma et al. Unified constraint parameter solutions for axial and circumferential surface cracks in pressurized pipes under creep condition
CN104614283B (zh) 一种金属材料热处理加工过程中的所对应物相变化的分析方法
Somerday et al. Effect of applied K level on the crack-arrest threshold in hydrogen environments: mechanics-based interpretation
CN107843509A (zh) 基于室温布氏硬度预测超临界机组t/p92耐热钢剩余持久寿命评估方法
Krechkovska et al. Substantiation of the critical structural and mechanical state of low-alloy heat-resistant steel from steam pipelines of thermal power plant
CN109187543A (zh) 一种在役乙烯裂解炉管脆化分级寿命评估方法
Kang et al. Initiation of stress corrosion cracks in X80 and X100 pipe steels in near-neutral pH environment
JP2018163148A (ja) 鋼材の硫化物応力腐食割れ性の評価方法および評価装置
Chen et al. Experimental investigation of various corroded metal materials used in building structures
JP6863152B2 (ja) 靱性予測装置、靱性予測方法、およびプログラム
CN109916737A (zh) 一种测试油井管抗应力腐蚀性能的方法
Kitade et al. Clarification of micromechanism on Brittle Fracture Initiation Condition of TMCP Steel with MA as the trigger point
Latypova et al. Hydrogen-induced cracking of 500 HBW steels studied using a novel tuning-fork test with integrated loadcell system
Das et al. Failure analysis of cold worked AISI 301 SS diaphragm of gas pump
Dzioba et al. Temperature dependency of fracture toughness of high-strength ferritic steel Hardox-400
CN108760501A (zh) 由紧凑拉伸试件确定岩石断裂韧度与拉伸强度的方法

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