CN109883823A - 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法 - Google Patents

金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法 Download PDF

Info

Publication number
CN109883823A
CN109883823A CN201910052688.0A CN201910052688A CN109883823A CN 109883823 A CN109883823 A CN 109883823A CN 201910052688 A CN201910052688 A CN 201910052688A CN 109883823 A CN109883823 A CN 109883823A
Authority
CN
China
Prior art keywords
constriction
point
moment
radius
hyperbola
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
CN201910052688.0A
Other languages
English (en)
Other versions
CN109883823B (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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN201910052688.0A priority Critical patent/CN109883823B/zh
Publication of CN109883823A publication Critical patent/CN109883823A/zh
Application granted granted Critical
Publication of CN109883823B publication Critical patent/CN109883823B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法,属于金属材料力学性能测试技术领域,该方法基于金属圆棒试样单轴拉伸失稳颈缩逐级扩展的双曲线旋转体模型,只需要测量圆棒试样拉伸过程中标距伸长量,并提取断后试棒外轮廓曲线坐标信息,通过反推计算每一时刻颈缩双曲线方程和颈缩最小截面半径,即可通过曲率半径公式求得此时刻颈缩处最小截面处的曲率半径,最后通过陈篪法计算该时刻真实应力和真实应变。本发明省去现有方法中颈缩最小截面半径和外轮廓曲率半径的试验测量环节,并且提高了大应变范围应力应变曲线的测量精度,对于金属材料力学性能测试具有重要意义。

Description

金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法
技术领域
本发明属于金属材料力学性能测试技术领域,具体涉及金属圆棒试样单轴拉伸试验大应变范围应力应变曲线的测量方法。
背景技术
单轴拉伸试验是金属材料力学性能测试的一种基本方法,通常由单轴拉伸试验获得材料的载荷-标距伸长量曲线,然后假定试样标距范围内均匀变形,将载荷-标距伸长量曲线转换为应力-应变曲线。但是试样在载荷最大点之后开始颈缩,试样标距范围内不均匀变形。因此,拉伸失稳后由载荷-瞬时标距长度曲线直接转换为应力-应变曲线是不准确的,而且试样产生颈缩后逐渐由单轴应力状态向三轴应力状态转变。由于在载荷最大点之前对应的应变范围相对较小(一般在0.2以内),为测量更大应变范围的应力应变曲线,Bridgman法假设圆棒拉伸颈缩轮廓呈圆弧形状,提出了颈缩后应力修正公式,并在实践应用。但实践证明圆棒拉伸颈缩轮廓呈双曲线形状,因此Bridgman法修正应力精度很难保证,一般最大误差都在10%。为此,陈篪基于颈缩轮廓双曲线假设提出了更为合理的应力修正公式。但是,针对颈缩最小截面半径和颈缩外轮廓曲线在最小截面处的曲率半径的实际测量环节非常繁琐,一般都采用光学测量,测量成本很高,且测量精度不易保证,限制了该方法的实际应用。
发明内容
本发明的目的,为了解决目前陈篪法中颈缩最小截面半径和颈缩外轮廓曲线在最小截面处的曲率半径测量繁琐且误差大的难题,提出一种金属圆棒试样单轴拉伸试验大应变范围应力应变曲线测量方法,无需直接测量颈缩最小截面半径和颈缩外轮廓曲线在最小截面处的曲率半径,而是通过提取断后试棒颈缩外轮廓曲线进行反推。首先将金属圆棒试样单轴拉伸颈缩最小截面部位简化为母线为双曲线的旋转体结构,根据韧性金属圆棒试样单轴拉伸失稳后颈缩扩展是材料颈缩最小截面部位的逐渐局部化变形过程,颈缩发展可视为前一时刻双曲线旋转体占一定百分比的材料变形发展到当前时刻的双曲线旋转体,由此逐渐积累构建金属圆棒试样单轴拉伸失稳颈缩逐级扩展的双曲线旋转体模型。只需要测量和记录圆棒试样拉伸标距伸长量和颈缩处最小截面半径,同时采用光学方法采集断后圆棒试样外轮廓曲线,基于颈缩双曲线旋转体模型反推前一时刻旋转体母线方程,计算母线对应的颈缩最小截面半径和颈缩外轮廓在最小截面处的曲率半径,最后通过陈篪法计算该时刻真实应力和真实应变。由此逐渐反推获得金属圆棒试样单轴拉伸试验不同时刻的应力应变计算结果,构成材料大应变范围内应力应变曲线。
本发明所采用的技术方案具体步骤如下:
(1)进行金属圆棒试样单轴拉伸试验,原始标距长度为l0,原始标距半径为a0,测量并记录开始拉伸直至断裂前的各时刻载荷Fi和标距伸长量dli,其中i=0~N,0和N分别对应拉伸开始和断裂点时刻,并用光学方法采集断后较长的一半试棒沿轴线方向的外轮廓信息,通过图像处理获得外轮廓曲线的各点坐标其中k=0~M,0和M分别对应试样颈缩断裂点和平行长度范围内接近圆弧过渡部分的某点,断裂点的确定是以yk值最小为准,以对应颈缩断裂点的轴线坐标为坐标原点,即对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数
(2)对断后试样颈缩外轮廓曲线坐标进行弹性修正,已知金属材料弹性模量E和泊松比λ,则修正后时刻i对应的各点坐标表示为
其中保留所有的点(数量为T+1),并将修正为(l0+dlN)/2,则对应断裂时刻标距点的坐标可表示为
对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数φ(x)。
(3)金属圆棒试样在颈缩失稳后始终在最小颈缩处扩散变形,而最小颈缩处之外不参与变形。将时刻i对应的颈缩最小截面部位简化为双曲线母线方程为fi(x)的旋转体,并与之前所有时刻的双曲线旋转体构成圆棒试样拉伸颈缩扩展的双曲线旋转体模型。时刻i母线fi(x)与颈缩外轮廓曲线函数φ(x)相切并交于Pi(xi,yi)点,其在最小截面处的坐标为Qi(x′i,y′i)点。
则有
其中ai和bi分别为i时刻颈缩最小截面半径和双曲线母线方程的待求参数。
(4)基于体积不变定律,时刻i对应的Pi(xi,yi)点至Qi(x′i,y′i)点的双曲线fi(x)旋转体体积Vi1和标距点至Pi(xi,yi)点的颈缩外轮廓曲线函数φ(x)旋转体体积Vi2的和等于初始标距范围内试样体积可表示为如下关系,
通过对上式求解,可计算出各时刻i对应未知量的xi、ai和bi值,并确定双曲线母线fi(x)方程。通过曲率半径公式
求得双曲线母线fi(x)在Qi(x′i,y′i)点处的曲率半径Ri,即颈缩处曲率半径。其中i=N时刻的颈缩曲率半径为,
(5)根据陈篪法,将i时刻的载荷Fi、颈缩处最小截面半径ai和曲率半径Ri,代入下式,
计算得到时刻i对应的真实应变εi和真实应力σi
(6)重复以上(3)~(5)步骤,分别计算i=0~N时刻的真实应变εi和真实应力σi,从而获取金属圆棒试样从拉伸开始至断裂结束范围内的应力应变曲线。
本发明的有益效果:
1)采用本发明的方案,只需要测量和记录圆棒试样标距伸长量,并提取断后试棒颈缩外轮廓曲线,与现有技术相比省去颈缩最小截面半径和外轮廓曲率半径的繁琐测量环节,节省了测量时间和测量成本;
2)本发明确定双曲线旋转体母线方程,通过曲率半径公式得到颈缩处最小截面处的曲率半径,计算精度远高于现有方法中曲率半径的试验测量精度,因此本发明测量精度远高于现有方法的测量精度。另外,本发明测量公式易于实现程序化,具有操作性。
附图说明
图1为圆棒试样单轴拉伸载荷位移曲线。
图2为单轴拉伸断后照片和外轮廓曲线坐标信息。
图3为圆棒试样颈缩双曲线旋转体反推模型。
图4为单轴拉伸大应变范围应力应变曲线。
具体实施方式
下面结合具体示例及附图对本发明作进一步说明。
实施例
(1)进行金属圆棒试样单轴拉伸试验,原始标距长度为l0为50mm,原始标距半径为a0为5mm,测量并记录开始拉伸直至断裂前的各时刻载荷Fi和标距伸长量dli,其中i=0~N,0和N分别对应拉伸开始和断裂点时刻,载荷-位移曲线如图1所示。如图2所示,用光学方法采集断后较长的一半试棒沿轴线方向的外轮廓信息,通过图像处理获得外轮廓曲线的各点坐标其中k=0~M,0和M分别对应试样颈缩断裂点和平行长度范围内接近圆弧过渡部分的某点,断裂点的确定是以yk值最小为准,以对应颈缩断裂点的轴线坐标为坐标原点,即对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数
(2)对断后试样颈缩外轮廓曲线坐标进行弹性修正,已知金属材料弹性模量E和泊松比λ,则修正后时刻i对应的各点坐标表示为
其中保留所有的点(数量为T+1),并将修正为(l0+dlN)/2,则对应断裂时刻标距点的坐标可表示为
对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数φ(x)。
(3)金属圆棒试样在颈缩失稳后始终在最小颈缩处扩散变形,而最小颈缩处之外不参与变形。将时刻i对应的颈缩最小截面部位简化为双曲线母线方程为fi(x)的旋转体,并与之前所有时刻的双曲线旋转体构成圆棒试样拉伸颈缩扩展的双曲线旋转体模型,如图3所示。时刻i母线fi(x)与颈缩外轮廓曲线函数φ(x)相切并交于Pi(xi,yi)点,其在最小截面处的坐标为Qi(x′i,y′i)点。则有
其中ai和bi分别为i时刻颈缩最小截面半径和双曲线母线方程的待求参数。
(4)基于体积不变定律,时刻i对应的Pi(xi,yi)点至Qi(x′i,y′i)点的双曲线fi(x)旋转体体积Vi1和标距点至Pi(xi,yi)点的颈缩外轮廓曲线函数φ(x)旋转体体积Vi2的和等于初始标距范围内试样体积可表示为如下关系,
通过对上式求解,可计算出各时刻i对应未知量的xi、ai和bi值,并确定双曲线母线fi(x)方程。通过曲率半径公式
求得双曲线母线fi(x)在Qi(x′i,y′i)点处的曲率半径Ri,即颈缩处曲率半径。其中i=N时刻的颈缩曲率半径为,
(5)根据陈篪法,将i时刻的载荷Fi、颈缩处最小截面半径ai和曲率半径Ri,代入下式,
计算得到时刻i对应的真实应变εi和真实应力σi
(6)重复以上(3)~(5)步骤,分别计算i=0~N时刻的真实应变εi和真实应力σi,从而获取金属圆棒试样从拉伸开始至断裂结束范围内的应力应变曲线,如图4所示。

Claims (2)

1.金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法,其特征在于,该测量方法的具体步骤如下:
(1)进行金属圆棒试样单轴拉伸试验,原始标距长度为l0,原始标距半径为a0,测量并记录开始拉伸直至断裂前的各时刻载荷Fi和标距伸长量dli,其中i=0~N,0和N分别对应拉伸开始和断裂点时刻,并用光学方法采集断后较长的一半试棒沿轴线方向的外轮廓信息,通过图像处理获得外轮廓曲线的各点坐标其中k=0~M,0和M分别对应试样颈缩断裂点和平行长度范围内接近圆弧过渡部分的某点,断裂点的确定是以yk值最小为准,以对应颈缩断裂点的轴线坐标为坐标原点,即对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数
(2)对断后试样颈缩外轮廓曲线坐标进行弹性修正,已知金属材料弹性模量E和泊松比λ,则修正后时刻i对应的各点坐标表示为
其中保留所有的点(数量为T+1),并将修正为(l0+dlN)/2,则对应断裂时刻标距点的坐标可表示为
对颈缩外轮廓曲线各点进行拟合,得到颈缩外轮廓曲线函数φ(x)。
(3)金属圆棒试样在颈缩失稳后始终在最小颈缩处扩散变形,而最小颈缩处之外不参与变形。将时刻i对应的颈缩最小截面部位简化为双曲线母线方程为fi(x)的旋转体,并与之前所有时刻的双曲线旋转体构成圆棒试样拉伸颈缩扩展的双曲线旋转体模型。时刻i母线fi(x)与颈缩外轮廓曲线函数φ(x)相切并交于Pi(xi,yi)点,其在最小截面处的坐标为Qi(x′i,y′i)点。则有
其中ai和bi分别为i时刻颈缩最小截面半径和双曲线母线方程的待求参数。
(4)基于体积不变定律,时刻i对应的Pi(xi,yi)点至Qi(x′i,y′i)点的双曲线fi(x)旋转体体积Vi1和标距点至Pi(xi,yi)点的颈缩外轮廓曲线函数φ(x)旋转体体积Vi2的和等于初始标距范围内试样体积可表示为如下关系,
通过对上式求解,可计算出各时刻i对应未知量的xi、ai和bi值,并确定双曲线母线fi(x)方程。通过曲率半径公式
求得双曲线母线fi(x)在Qi(x′i,y′i)点处的曲率半径Ri,即颈缩处曲率半径。其中i=N时刻的颈缩曲率半径为,
(5)根据陈篪法,将i时刻的载荷Fi、颈缩处最小截面半径ai和曲率半径Ri,代入下式,
计算得到时刻i对应的真实应变εi和真实应力σi
2.根据权利要求1所述的金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法,其特征在于,重复以上(3)~(5)步骤,分别计算i=0~N时刻的真实应变εi和真实应力σi,从而获取金属圆棒试样从拉伸开始至断裂结束范围内的应力应变曲线。
CN201910052688.0A 2019-01-21 2019-01-21 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法 Expired - Fee Related CN109883823B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910052688.0A CN109883823B (zh) 2019-01-21 2019-01-21 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910052688.0A CN109883823B (zh) 2019-01-21 2019-01-21 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法

Publications (2)

Publication Number Publication Date
CN109883823A true CN109883823A (zh) 2019-06-14
CN109883823B CN109883823B (zh) 2021-05-11

Family

ID=66926369

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910052688.0A Expired - Fee Related CN109883823B (zh) 2019-01-21 2019-01-21 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法

Country Status (1)

Country Link
CN (1) CN109883823B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110702513A (zh) * 2019-10-15 2020-01-17 吉林大学 一种金属棒材大应变范围硬化曲线的试验测量方法
CN111189701A (zh) * 2020-01-08 2020-05-22 吉林大学 一种金属双曲线试样的大应变压缩硬化曲线的测量方法
CN111272551A (zh) * 2020-02-17 2020-06-12 本钢板材股份有限公司 一种高速拉伸试验数据曲线处理方法
CN112729795A (zh) * 2020-12-25 2021-04-30 北京航空航天大学 一种叶根缘板过渡特征模拟试样
CN112927185A (zh) * 2021-01-19 2021-06-08 中国石油天然气集团有限公司 一种基于数字图像相关法的真应力-真应变曲线测试计算方法
CN113281171A (zh) * 2021-06-18 2021-08-20 中国船舶重工集团公司第七二五研究所 一种金属材料拉伸颈缩过程中应力分布测试方法
CN116202874A (zh) * 2023-05-05 2023-06-02 青岛宇通管业有限公司 一种排水管材柔韧性测试方法及系统
CN116858665A (zh) * 2023-06-09 2023-10-10 洛阳船舶材料研究所(中国船舶集团有限公司第七二五研究所) 一种圆棒试样单轴拉伸颈缩变形时外形轮廓的分析方法
CN116907985A (zh) * 2023-06-09 2023-10-20 洛阳船舶材料研究所(中国船舶集团有限公司第七二五研究所) 拉伸试样颈缩变形的最小截面中心应力与应变的分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010026356A1 (fr) * 2008-09-08 2010-03-11 Arkema France Procédé de détermination de la tenue à la fatigue d'une composition polymérique
CN102998175A (zh) * 2011-09-15 2013-03-27 北汽福田汽车股份有限公司 一种测试金属材料断裂伸长率的方法
CN103091167A (zh) * 2013-01-23 2013-05-08 西北工业大学 一种连续测定钛合金管材收缩应变比变化的方法
CN103175735A (zh) * 2012-12-06 2013-06-26 西南交通大学 材料拉伸真实本构曲线测试技术
JP5893542B2 (ja) * 2012-10-04 2016-03-23 公益財団法人鉄道総合技術研究所 長繊維強化金属基複合材料の引張強度試験方法及び引張試験片
CN107271279A (zh) * 2017-06-16 2017-10-20 西南石油大学 一种金属材料真应力真应变测试方法
CN108982223A (zh) * 2018-09-19 2018-12-11 吉林大学 金属圆棒试样单轴拉伸试验大应变范围应力应变的测量方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010026356A1 (fr) * 2008-09-08 2010-03-11 Arkema France Procédé de détermination de la tenue à la fatigue d'une composition polymérique
CN102998175A (zh) * 2011-09-15 2013-03-27 北汽福田汽车股份有限公司 一种测试金属材料断裂伸长率的方法
JP5893542B2 (ja) * 2012-10-04 2016-03-23 公益財団法人鉄道総合技術研究所 長繊維強化金属基複合材料の引張強度試験方法及び引張試験片
CN103175735A (zh) * 2012-12-06 2013-06-26 西南交通大学 材料拉伸真实本构曲线测试技术
CN103091167A (zh) * 2013-01-23 2013-05-08 西北工业大学 一种连续测定钛合金管材收缩应变比变化的方法
CN107271279A (zh) * 2017-06-16 2017-10-20 西南石油大学 一种金属材料真应力真应变测试方法
CN108982223A (zh) * 2018-09-19 2018-12-11 吉林大学 金属圆棒试样单轴拉伸试验大应变范围应力应变的测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈俊甫 等: "金属棒试样拉伸和扭转试验应变范围和力学特性对比", 《吉林大学学报(工学版)》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110702513B (zh) * 2019-10-15 2021-06-18 吉林大学 一种金属棒材大应变范围硬化曲线的试验测量方法
CN110702513A (zh) * 2019-10-15 2020-01-17 吉林大学 一种金属棒材大应变范围硬化曲线的试验测量方法
CN111189701A (zh) * 2020-01-08 2020-05-22 吉林大学 一种金属双曲线试样的大应变压缩硬化曲线的测量方法
CN111272551B (zh) * 2020-02-17 2022-08-05 本钢板材股份有限公司 一种高速拉伸试验数据曲线处理方法
CN111272551A (zh) * 2020-02-17 2020-06-12 本钢板材股份有限公司 一种高速拉伸试验数据曲线处理方法
CN112729795A (zh) * 2020-12-25 2021-04-30 北京航空航天大学 一种叶根缘板过渡特征模拟试样
CN112927185A (zh) * 2021-01-19 2021-06-08 中国石油天然气集团有限公司 一种基于数字图像相关法的真应力-真应变曲线测试计算方法
CN113281171A (zh) * 2021-06-18 2021-08-20 中国船舶重工集团公司第七二五研究所 一种金属材料拉伸颈缩过程中应力分布测试方法
CN113281171B (zh) * 2021-06-18 2023-10-20 中国船舶重工集团公司第七二五研究所 一种金属材料拉伸颈缩过程中应力分布测试方法
CN116202874A (zh) * 2023-05-05 2023-06-02 青岛宇通管业有限公司 一种排水管材柔韧性测试方法及系统
CN116858665A (zh) * 2023-06-09 2023-10-10 洛阳船舶材料研究所(中国船舶集团有限公司第七二五研究所) 一种圆棒试样单轴拉伸颈缩变形时外形轮廓的分析方法
CN116907985A (zh) * 2023-06-09 2023-10-20 洛阳船舶材料研究所(中国船舶集团有限公司第七二五研究所) 拉伸试样颈缩变形的最小截面中心应力与应变的分析方法
CN116907985B (zh) * 2023-06-09 2024-08-23 洛阳船舶材料研究所(中国船舶集团有限公司第七二五研究所) 拉伸试样颈缩变形的最小截面中心应力与应变的分析方法

Also Published As

Publication number Publication date
CN109883823B (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
CN109883823A (zh) 金属圆棒试样单轴拉伸应力应变的双曲线反推测量方法
CN109883825A (zh) 金属圆棒试样单轴拉伸大应变范围硬化曲线的测量方法
CN108982223B (zh) 金属圆棒试样单轴拉伸试验大应变范围应力应变的测量方法
CN110702513B (zh) 一种金属棒材大应变范围硬化曲线的试验测量方法
CN107991198A (zh) 拉弯扭三向复合引伸计
CN110411836B (zh) 一种测定金属材料塑性应变比的试验与计算方法
CN108982222A (zh) 金属板试样单轴拉伸试验大应变范围应力应变的测量方法
CN109883824A (zh) 金属圆棒试样单轴拉伸应力应变的圆弧反推测量方法
CN207600867U (zh) 一种塑料弯曲度测试装置
CN201096489Y (zh) 多用途测厚仪
CN107917773A (zh) 一种管棒形材料残余应力无压痕应变测定装置及方法
CN103335581B (zh) 内螺纹接头中径定位样板及中径测量方法
CN116907985B (zh) 拉伸试样颈缩变形的最小截面中心应力与应变的分析方法
CN107024401B (zh) 获取金属材料各向异性和拉压非对称性的方法及系统
CN109059734A (zh) 一种电气工程用电缆转弯半径测量仪
CN105371769A (zh) 一种测量动态裂纹长度的方法及装置
CN210426367U (zh) 发动机的曲轴轴瓦半径高度的测量装置
CN107238347B (zh) 应力应变测量系统中应变计连接导线的自动修正方法
CN207036559U (zh) 一种整o形橡胶圈拉伸试验断裂伸长量测量装置
CN111623909B (zh) 一种双叉臂二力杆的受力测试方法
CN202195789U (zh) 单向通止规
CN217338576U (zh) 一种子宫脱垂测量尺
CN209117485U (zh) 一种混凝土轴向拉伸变形装置
CN206095136U (zh) 一种混凝土膨胀率/收缩率自动测试装置
CN217717852U (zh) 测试棒辅助测量夹具

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210511

Termination date: 20220121

CF01 Termination of patent right due to non-payment of annual fee