CN102830134B - 采用上下恒温参数辨识法测热界面材料性能的方法 - Google Patents

采用上下恒温参数辨识法测热界面材料性能的方法 Download PDF

Info

Publication number
CN102830134B
CN102830134B CN201210278771.8A CN201210278771A CN102830134B CN 102830134 B CN102830134 B CN 102830134B CN 201210278771 A CN201210278771 A CN 201210278771A CN 102830134 B CN102830134 B CN 102830134B
Authority
CN
China
Prior art keywords
thermal flow
flow meter
temperature
thermal
test point
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
CN201210278771.8A
Other languages
English (en)
Other versions
CN102830134A (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 Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201210278771.8A priority Critical patent/CN102830134B/zh
Publication of CN102830134A publication Critical patent/CN102830134A/zh
Application granted granted Critical
Publication of CN102830134B publication Critical patent/CN102830134B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种上下恒温参数辨识法测热界面材料性能的方法,包括以下步骤:第一步,测试标准材料热流量计的准备;第二步,两热流量计接触界面之间放置热界面材料,加载压应力,正向对热流量计加热;第三步,测试点测量温度的修正;第四步,接触热阻R的计算;第五步,热界面材料厚度L测量,热界面材料的当量导热系数的计算。本方法采用上下恒温参数辨识修正的对称测试结构进行测量可基本消去由于各温度传感器与热流量计的接触热阻不同或是温度不完全是线性的等原因引起的温度测量的不一致性误差,进而在保证热流量精度的前提下可极高精度的测得热界面材料的性能。

Description

采用上下恒温参数辨识法测热界面材料性能的方法
技术领域
本发明属于测试技术领域,具体涉及一种接触热阻测试方法及设备,适用于对常用材料的界面接触热阻的测试,尤其适用于对热界面材料的性能测试。
背景技术
接触热阻是一个受材料物性、机械特性、表面形貌、接触压力、温度、间隙材料等众多因素影响的参数。根据实验热流是否稳定,一般把接触热阻测量方法分为瞬态法和稳态法。瞬态法也是一种常用的接触热阻实验测量方法,其主要包括激光光热测量法、热成像法、“flash”闪光法、激光光声法等,其中激光光热测量法又包含调制光热法和热扫描法,调制光热法又有光热幅值法、光热相位法和脉冲法之分。虽然各种瞬态法虽宜于快速测量且可测量小到纳米数量级的薄膜,但其测量过程易受各种因素影响,且公式推导相对复杂,测量精度也较难保证。因此,界面接触热阻测量方法最常用的是稳态法:在两接触样品上维持一定的温差,测量两样品轴向上的温度值,再由傅里叶定律外推至接触界面处从而得到界面上的温差;热流量可由热流量计测量或由样品材料的热导率和温度梯度计算得到,从而R=|T1-T2|/Q。稳态接触热阻测试方法多是和美国国家标准ASTMD5470-06的测试标准设备相类似,但多有文献指出由于温度测量的不确定性误差和热损失误差很难保证对界面接触热阻有足够高的测量精度。
发明内容
本发明的目的在于提供一种上下恒温参数辨识法测热界面材料接触热阻,在保证热流量精度的前提下可极高精度的测得热界面材料的接触热阻和有效当量导热系数。
实现本发明目的的技术解决方案为:上下恒温参数辨识法测热界面材料性能,其特征在于如下步骤:
第一步,测试标准材料热流量计的准备。
加工出两个标准材料的热流量计,将热流量计竖直同轴夹装在两加热制冷套之间,在两加热制冷套上设置有应力加载装置,所述的热流量计上设置有温度传感器,温度传感器与数据采集***连接,用于测试试样的轴向温度;
热流量计上测试点之间的位置满足如下关系:以两个热流量计轴向方向上的接触界面截面位置为对称面,两个热流量计上的测试点位置完全对称,每个热流量计上的相邻两个测试点之间的轴向距离相等,每个热流量计从下端面到上端面之间设置n个测试点,测试点之间的距离为dx;
第二步,两热流量计接触界面之间放置热界面材料,加载压应力,正向对热流量计加热:
对两个热流量计轴向的其中的一端加热,另一端冷却,试样温度达到稳定后开始采集测试温度;所述的测试温度包括各热流量计上n个测试点的测量温度Ti,j,i=1,······n,n为热流量计上测试点数目并按对称面对称,j=1,2分别表示两不同的热流量计;
第三步,测试点测量温度的修正;
在做好充分的绝热条件下,对热流量计的两端同时设定一个恒定温度,试样温度达到稳定后开始采集测试温度;所述的测试温度包括热流量计上n个测试点的测量温度i=1,······n,n为热流量计上测试点数目;
对第二步中所采集的n个测试点的温度测量范围根据精度需要进行上述的多恒定温度点重复进行热流量计上n个测试点的测量温度采集,并把热流量计上各个测试点的测量温度与所设定的恒定温度进行参数辨识分析,进行线性拟合或多元拟合成相关函数;
第四步,接触热阻R的计算;
把第三步中相关函数对步骤二中采集的各个测试点的测量温度进行求解,得到一修正温度i=1,······n,n为热流量计上测试点数目;
进而在忽略热流损失的情况下,可较高精度的计算得到热流量计的接触热阻R;
第五步,热界面材料厚度L测量,热界面材料的当量导热系数的计算;
通过在两热流量计接触界面位置布置的原位测量***的参考点位置变化测得热界面材料的厚度L,计算表观接触热阻RA为:RA=A×R,其中A为接触面积,从而当量有效导热系数keff为:
k eff = L R A .
为保证温度梯度的一维性,热流量计的为圆柱体或长方体。
在正反双向测试时对接触界面温度Ts-1′、Ts-1′′和Ts-2′、Ts-2′′的计算也可采用最小平方法进行线性拟合求解或反问题方法求解。
在两个热流量计的接触界面位置装有原位测量***。
本发明与现有技术相比,本发明所述的一种采用上下恒温参数辨识法测热界面材料的方法采用上下恒温参数辨识修正的对称测试结构进行测量可基本消去由于各温度传感器与热流量计的接触热阻不同或是温度不完全是线性的等原因引起的温度测量的不一致性误差,进而在保证热流量精度的前提下可极高精度的测得热流量计的接触热阻。
附图说明
图1为本发明测试装置的主视图;
图2为本发明的***测试原理图;
图3为本发明中标准热流量计1的主视图;
图4为采用本发明方法在上下恒温测试时热流量计1和2上的温度偏差;
图5为采用本发明方法对热流量计1和2上的温度偏差进行相关函数拟合时的参数;
图6为采用本发明方法在压力2MPa时分别对上加热和下加热所得到的接触热阻进行恒温温度修正后得到的接触热阻随热流量的关系;
具体实施方式
本发明在美国国家标准ASTM D5470基础上提出了一种采用上下恒温参数辨识法测热界面材料的性能,所述的测试方法是采用可上下恒温的对称测试结构,通过对测试热流量计进行上下恒定温度参数辨识消去各温度传感器的不一致性误差和热损失误差,结合可控温热辐射防辐射屏来及辅助措施减小热流损失,来达到高精度测试热界面材料的热物性参数的目的,本方法可高精度的测量填充热界面材料后的界面接触热阻和有效当量导热系数。
下面结合附图对本发明作进一步详细描述。
结合图1,本发明公开了一种采用上下恒温参数辨识法测热界面材料的装置,该装置为上下正反双向热流测试的对称结构,包括控制***、支架3、第一滚珠套筒4-1、第二滚珠套筒4-2、滑动螺杆5、定向钢球和压力传感器6、辅助加热器7、真空罩9、试件测试区10、应力加载装置、真空抽放气口13、进出水口14、数据采集***、密封底盘16,支撑板17,水平调节杆20和加热丝21;其特征在于:应力加载装置由液压缸11和压力动力源12组成,液压缸11位于压力动力源12的上方;数据采集***由温度传感器、密封数据接头15组成,温度传感器通过导线与密封数据接头15相连;控制***由可控温防辐射屏2、加热制冷套1和控制防辐射屏加热丝R2组成;试样测试区10包括测试试件,其中定向钢球和压力传感器6、支架3、支撑板17和加热制冷套1上下对称,定向钢球和压力传感器6固定在支撑板17中心位置,应力加载装置通过支架3定位并和定向钢球以及压力传感器6接触,为试样加载应力,第一滚珠套筒4-1设置在滑动螺杆5的上下两端与支撑板17固定,第二滚珠套筒4-2设置在滑动螺杆5的底部并与支架3固定,辅助加热器7位于支撑板17和加热制冷套1之间,试样测试区10位于上下对称的两个加热制冷套1之间,两个可控温防辐射屏2位于试样测试区10的外部,真空罩9位于整个装置的外部固定于密封底盘16上,滑动螺杆5固定于密封底盘16的上部,真空抽放气口13、进出水口14和密封数据接头15均设置在密封底盘16上,液压缸11贯穿密封底盘16的中心,密封底盘上设置有四组水平调节杆20。
图2为本发明的测试原理示意图,在进行测试过程中,根据热流量计上的温度传感器的测量温度来调控防辐射屏上的加热装置使其达到和热流量计近似的温度梯度以此来减小热流损失。在上下加热制冷套的位置也相应布置有辅助加热器来调控和加热源近似的温度来减小热损失。
在图3中,本发明中插装有温度传感器的标准热流量计1的主视图,在该热流量计上插装有位置上下对称度有严格要求的3组温度传感器。标准热流量计可加工成圆柱体或者长方体,温度传感器的插装孔有严格的位置精度和形状精度要求,并保证有足够的上下对称度,在插装温度传感器前对标准热流量计(本案例选用的是Elkonite copper-tungsten alloy30W3材料,导热系数为216±2W/m K,硬度为276HV)进行酒精、丙酮、异丙酮和超声波清洗。温度传感器为对称等距排列,温度传感器的探头是通过焊接或导热膏粘结插装孔里。本实施例采用的温度传感器为热电阻。
本发明公开了一种上下恒温参数辨识法测定石墨垫片热界面材料的接触热阻的方法,其测试步骤如下:
第一步,测试标准材料热流量计的准备。
如图1和图2所示,上下热流量计之所以选择30W3钨铜合金考虑到以下原因:由于铜的良好的导热性能和钨的高强度使之在材料硬度和导热性能之间达到很好的平衡,所以可在长时间的使用过程中减小破坏接触端面表面形貌的可能性。制作出两标准材料(Elkonite copper-tungsten alloy30W3材料)热流量计,加工出两钨铜合金材料热流量计,将热流量计竖直安装在两个上下对称设置的加热制冷套之间,在两加热制冷套上设置有应力加载装置,所述的热流量计上设置有温度传感器,温度传感器与数据采集***连接,用于测试热流量计的轴向温度,若所采用的温度传感器为热电偶,则根据该热流量计轴向截面均匀布置的1-4个热电偶平均来求得该轴向点的温度;若所采用的温度传感器为热电阻,则对热流量计轴向截面均匀布置测温的1-4个热电阻采用4线制接法,对该均匀布置的1-4个热电阻的激励电流相同,对该均匀布置的1-4个热电阻信号采集端的引线可采用并联接法来平均求得该轴向点的温度并减小因热电阻的引线造成的热流损失。
热流量计上测试点之间的位置满足如下关系:以两热流量计纵轴方向上的接触界面截面位置为对称面,两热流量计上的测试点位置完全对称,每个热流量计从下端面到上端面之间均设置4个测试点,每个热流量计上相邻两个测试点之间的轴向距离相等,测试点之间的距离为dx=25mm,从接触界面到一个测试点的位置为2mm,如图2所示热流量计1的(T.x)3测试点到接触界面的距离为2mm,热流量计2同样从接触界面到一个测试点的位置为2mm。并按温度传感器尺寸在标准材料热流量计和热流量计上等距的加工出温度传感器的探头安装孔,所述的温度传感器的探头安装孔≤0.5mm,探头安装孔里通过焊接或导热膏粘结≤0.5mm的温度传感器探头,温度传感器通过真空腔壁的连接器与数据采集***连接,本发明的温度传感器采用四线制热电阻。
第二步,两热流量计接触界面之间放置热界面材料,加载压应力,正向对热流量计加热:
如图1和图2所示所示将布置有3组温度传感器的的热流量计1和2竖直夹装在两端都有热流量计、加热制冷套、辅助加热装置的真空腔中,为减小热损失在保温层外层加设一内嵌有加热装置的可控温防辐射屏,在抽真空后进行顶端加热底端制冷的正向热流测试,此时可控温防辐射屏模拟出近似热流量计的温度梯度,顶端布置的辅助加热器根据加热制冷套的温度控制其温度以减小纵轴向的热损失,达到稳态时进行温度数据采集,此时加载功率可通过上下对称布置的热流量计的来换算热流量。
例如在压力2MPa,在热流量为3.5W时,测得热流量计1上的(T.x)1=25.053、(T.x)2=24.658、(T.x)3=24.257,热流量计2上的(T.x)4=23.038、(T.x)5=22.64和(T.x)6=22.243。
第三步,测试点测量温度的修正:
在做好充分的绝热条件下,对步骤二中所采集的热流量计1上的(T.x)1、(T.x)2和(T.x)3以及热流量计2上的(T.x)4、(T.x)5和(T.x)6共6个测试点对上下对称布置的加热制冷套从22~28℃的进行7个恒定温度点的稳态采集数据。如图4所示,设定加热制冷套的恒定温度并记录(T.x)1、(T.x)2、(T.x)3、(T.x)4、(T.x)5和(T.x)6测试点的测量温度与该恒定温度的偏差。本案例按进行y=A0×x3+A1×x2+A2×x+A3的多元方程函数对该温度偏差与恒定温度进行数据拟合,拟合的参数A0、A1、A2和A3如图5所示。
第四步,接触热阻R的计算:
如图6所示,例如在热流量为3.5W时,把(T.x)1=25.053、(T.x)2=24.658、(T.x)3=24.257以及(T.x)4=23.038、(T.x)5=22.64和(T.x)6=22.243的温度值代入第三步的共6个测试点的拟合相关函数中得到温度偏差值,从而得到
( T ~ . x ) 1 = 25.039 , ( T ~ . x ) 2 = 24.625 , ( T ~ . x ) 3 = 24.198 , ( T ~ . x ) 4 = 23.061 , ( T ~ . x ) 5 = 22.642 , ( T ~ . x ) 6 = 22.241 .
在加正向热流测试时,在热流量为3.5W时,把修正后的根据热流量计1上的(T.x)1、(T.x)2和(T.x)3与3个测试点位置的修正后的温度梯度关系,以及热流量计2上的(T.x)4、(T.x)5和(T.x)6与3个测试点位置的修正后的温度梯度关系,通过数值方法外推得到的热流量计1的外推界面温度为热流量计2的外推界面温度为
两个热流量计的界面温差为:
继而接触热阻为: R ~ = Δ T ~ s Q = 0.307 K / W .
其中Q为热流量。
如图6所示,在压力2MPa,热流量在3.5W时,上加热由热流量计1和热流量计2上没有修正的温度(T.x)5=25.053、(T.x)6=24.658、(T.x)7=24.257、(T.x)9=23.038、(T.x)10=22.64和(T.x)11=22.243计算得到的接触热阻值为0.33K/W,而采用修正后的温度 ( T ~ . x ) 5 = 25.03866 , ( T ~ . x ) 6 = 24.6249 , ( T ~ . x ) 7 = 24.19786 , ( T ~ . x ) 9 = 23.061 , ( T ~ . x ) 10 = 22.642 ( T ~ . x ) 11 = 22.241 计算得到的接触热阻值为0.307 K/W。同理,对图6所示的采用上下不同热流方向计算得到的接触热阻值所测量的温度值进行修正,即得到图6所示的上下不同热流方向计算得到的接触热阻值与采用本方法温度修正后得到的接触热阻值同加载热流的关系。
第五步,热界面材料厚度L测量,热界面材料的当量导热系数的计算;
通过在两热流量计接触界面位置布置的原位测量***的参考点位置变化测得热界面材料的厚度L=0.2mm,计算表观接触热阻RA为:RA=A×R= 150.698mm2K/W,其中接触面积为A=490.874mm2,从而当量有效导热系数keff为:
k eff = L R A = 1.327 W / mK .
以上所述为本发明的较佳实施例的详细说明与图附,并非用来限制本发明,本发明的所有范围应以专利权利书所要求保护的范围为准,凡与本发明的设计思想及其类似变化的实施例、近似结构,都应包含于本发明的专利保护范围之中。

Claims (3)

1.一种采用上下恒温参数辨识法测热界面材料性能的方法,其特征在于如下步骤:
第一步,测试标准材料热流量计的准备;
加工出两个标准材料的热流量计,将热流量计竖直同轴夹装在两个加热制冷套之间,在两个加热制冷套上设置有应力加载装置,所述的热流量计上设置有温度传感器,温度传感器与数据采集***连接,用于测试试样的轴向温度;
热流量计上测试点之间的位置满足如下关系:以两个热流量计轴向方向上的接触界面截面位置为对称面,两个热流量计上的测试点位置完全对称,每个热流量计上的相邻两个测试点之间的轴向距离相等,每个热流量计从下端面到上端面之间设置n个测试点,测试点之间的距离为dx;
第二步,两热流量计接触界面之间放置热界面材料,加载压应力,正向对热流量计加热:
对两个热流量计轴向的其中的一端加热,另一端冷却,试样温度达到稳定后开始采集测试温度;所述的测试温度包括各热流量计上n个测试点的测量温度Ti,j,i=1,······n,n为热流量计上测试点数目并按对称面对称,j=1,2分别表示两不同的热流量计;
第三步,测试点测量温度的修正;
在做好充分的绝热条件下,对热流量计的两端同时设定一个恒定温度,试样温度达到稳定后开始采集测试温度;所述的测试温度包括热流量计上n个测试点的测量温度i=1,······n,n为热流量计上测试点数目;j=1,2分别表示两不同的热流量计;
对第二步中所采集的n个测试点的温度测量范围根据精度需要进行上述的多恒定温度点重复进行热流量计上n个测试点的测量温度采集,并把热流量计上各个测试点的测量温度与所设定的恒定温度进行参数辨识分析,进行线性拟合或多元拟合成相关函数;
第四步,接触热阻R的计算;
把第三步中相关函数对步骤二中采集的各个测试点的测量温度进行求解,得到一修正温度i=1,······n,n为热流量计上测试点数目;j=1,2分别表示两不同的热流量计;
进而在忽略热流损失的情况下,可较高精度的计算得到热流量计的接触热阻R;
第五步,热界面材料厚度L测量,热界面材料的当量导热系数的计算;
通过在两热流量计接触界面位置布置的原位测量***的参考点位置变化测得热界面材料的厚度L,计算表观接触热阻RA为:RA=A×R,其中A为接触面积,从而当量有效导热系数keff为:
k eff = L R A .
2.根据权利要求1所述的采用上下恒温参数辨识法测热界面材料性能的方法,其特征在于:热流量计为圆柱体或长方体。
3.根据权利要求1所述的采用上下恒温参数辨识法测热界面材料性能的方法,其特征在于:在正反双向测试时对接触界面温度Ts-1′、Ts-1″和Ts-2′、Ts-2″的计算采用最小平方法进行线性拟合求解或反问题方法求解。
CN201210278771.8A 2012-08-07 2012-08-07 采用上下恒温参数辨识法测热界面材料性能的方法 Active CN102830134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210278771.8A CN102830134B (zh) 2012-08-07 2012-08-07 采用上下恒温参数辨识法测热界面材料性能的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210278771.8A CN102830134B (zh) 2012-08-07 2012-08-07 采用上下恒温参数辨识法测热界面材料性能的方法

Publications (2)

Publication Number Publication Date
CN102830134A CN102830134A (zh) 2012-12-19
CN102830134B true CN102830134B (zh) 2014-07-16

Family

ID=47333345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210278771.8A Active CN102830134B (zh) 2012-08-07 2012-08-07 采用上下恒温参数辨识法测热界面材料性能的方法

Country Status (1)

Country Link
CN (1) CN102830134B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102980910B (zh) * 2012-11-30 2014-11-05 深圳市博恩实业有限公司 导热材料性能测试设备
CN104215661B (zh) * 2014-09-02 2017-02-22 兰州大学 基于超磁致伸缩智能材料的固体界面接触热阻测试装置
CN105388184A (zh) * 2015-12-17 2016-03-09 北京航空航天大学 用于接触热阻测试的试样安装夹具
CN106093108B (zh) * 2016-05-19 2018-10-16 南京航空航天大学 基于间隙缺陷识别的单向纤维增韧复合材料等效导热系数预估方法
CN109839406B (zh) * 2019-03-27 2022-08-09 桂林电子科技大学 一种界面接触热阻的高精度测试方法
CN110411397A (zh) * 2019-07-10 2019-11-05 南京航空航天大学 一种水平烟道底部积灰厚度的在线监测装置及方法
CN114062420B (zh) * 2020-08-10 2023-09-08 华为技术有限公司 一种导热参数的测试装置和方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6142662A (en) * 1998-06-16 2000-11-07 New Jersey Institute Of Technology Apparatus and method for simultaneously determining thermal conductivity and thermal contact resistance
CN101929970B (zh) * 2010-07-13 2012-05-16 北京航空航天大学 接触热阻测试方法
CN101929969B (zh) * 2010-07-13 2012-01-11 北京航空航天大学 带有冷却装置的接触热阻测试设备
CN101915780B (zh) * 2010-07-13 2012-03-28 北京航空航天大学 应用于高温高应力接触热阻测试的测试设备

Also Published As

Publication number Publication date
CN102830134A (zh) 2012-12-19

Similar Documents

Publication Publication Date Title
CN102768225B (zh) 一种高精度热界面材料测试方法
CN102830134B (zh) 采用上下恒温参数辨识法测热界面材料性能的方法
CN102768224B (zh) 正反双向热流法测固-固接触热阻的测试方法
CN102798645B (zh) 一种导热系数及接触热阻测试装置
CN102778475B (zh) 一种上下恒温参数辨识法测固-固接触热阻
CN101290299B (zh) 一种可变导热系数的测量装置及方法
CN101865865B (zh) 非接触式大坝混凝土早期热膨胀系数的测试方法及其设备
CN101343676A (zh) 高炉炉腔温度在线检测装置及检测方法
CN103148971B (zh) 一种测试超高压管式反应器热套端部结构局部应力场的方法
CN103217454B (zh) 圆筒结构热扩散率的光纤光栅测量方法
CN108333215B (zh) 一种整体式tps中气凝胶隔热层的热导率分析方法
CN102778474B (zh) 上下恒温参数辨识法测导热系数的方法
CN101021502A (zh) 低温电阻温度系数测试装置
CN105466633A (zh) 一种低温下压力传感器校准装置
CN108490237A (zh) 一种便携式金属管道热电势无损测量的装置及方法
CN108680284B (zh) 低温环境中光纤光栅温度传感器温度标定装置和标定方法
CN203849193U (zh) 一种室内岩样导热系数测试装置
CN203455295U (zh) 相变温度测试仪
CN102778476B (zh) 正反双向热流法测定导热系数的方法
CN102192922A (zh) 高温材料导热系数测量装置
CN112229871A (zh) 一种接触热阻测试装置及方法
CN109401956B (zh) 针对pcr仪的温度检测仪
CN100425977C (zh) 用位移传感器测量聚合物基泡沫材料线膨胀系数的方法
Sapozhnikov et al. The calibration of gradient heat flux sensors
CN102809445B (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
C14 Grant of patent or utility model
GR01 Patent grant