CN110489914A - 一种基于应力损伤等效的耐久性计算方法 - Google Patents

一种基于应力损伤等效的耐久性计算方法 Download PDF

Info

Publication number
CN110489914A
CN110489914A CN201910797068.XA CN201910797068A CN110489914A CN 110489914 A CN110489914 A CN 110489914A CN 201910797068 A CN201910797068 A CN 201910797068A CN 110489914 A CN110489914 A CN 110489914A
Authority
CN
China
Prior art keywords
equivalent
stress
damage
loading spectrum
durability
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
CN201910797068.XA
Other languages
English (en)
Other versions
CN110489914B (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.)
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Original Assignee
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
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 Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC filed Critical Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Priority to CN201910797068.XA priority Critical patent/CN110489914B/zh
Publication of CN110489914A publication Critical patent/CN110489914A/zh
Application granted granted Critical
Publication of CN110489914B publication Critical patent/CN110489914B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本申请属于飞机结构强度设计领域,特别涉及一种基于应力损伤等效的耐久性计算方法。包括:步骤一:计算第一载荷谱的当量过载损伤;步骤二:计算第二载荷谱下结构的当量应力损伤;步骤三:根据所述当量过载损伤和所述当量应力损伤等效折算所述结构在所述第一载荷谱下的等效应力;步骤四:根据所述第一载荷谱下数据集以及所述等效应力,对所述结构进行耐久性计算评估。本申请的基于应力损伤等效的耐久性计算方法,可以构建不同载荷谱下结构耐久性的关系,从而通过已有的(裂纹长度a,飞行小时t)数据集在新飞机结构设计早期进行耐久性计算评估。

Description

一种基于应力损伤等效的耐久性计算方法
技术领域
本申请属于飞机结构强度设计领域,特别涉及一种基于应力损伤等效的耐久性计算方法。
背景技术
随着时代的前进,对飞机飞行的安全性和使用的经济性提出了越来越高的要求。飞机结构的耐久性是一种结构具有在使用寿命期间承受重复载荷谱作用下而不产生功能性损坏或引起不经济维修等问题的特性,是飞机结构设计必须满足的结构特性。
在飞机结构设计的早期,因缺乏耐久性计算所需要的(裂纹长度a,飞行小时t)数据集,耐久性计算评估工作难以开展。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种基于应力损伤等效的耐久性计算方法,以解决现有技术存在的至少一个问题。
本申请的技术方案是:
一种基于应力损伤等效的耐久性计算方法,包括:
步骤一:计算第一载荷谱的当量过载损伤;
步骤二:计算第二载荷谱下结构的当量应力损伤;
步骤三:根据所述当量过载损伤和所述当量应力损伤等效折算所述结构在所述第一载荷谱下的等效应力;
步骤四:根据所述第一载荷谱下数据集以及所述等效应力,对所述结构进行耐久性计算评估。
可选地,步骤一中,所述计算第一载荷谱的当量过载损伤为:
GA=∑(n当量 i)mNi
n当量 i=[2n幅值(n幅值+|n均值|)]1/2
n幅值=(nmax-nmin)/2
n均值=(nmax+nmin)/2
其中,GA为第一载荷谱A的当量过载损伤;Ni为n当量 i的作用次数;m为应力-循环曲线的指数,取m=4;nmax为过载循环的峰值;nmin为过载循环的谷值。
可选地,步骤二中,所述计算第二载荷谱下结构的当量应力损伤为:
LB=∑(σ当量 i)mNi
σ当量 i=[2σ幅值幅值+|σ均值|)]1/2
σ幅值=(σmaxmin)/2
σ均值=(σmaxmin)/2
其中,LB为第二载荷谱B下结构C的当量应力损伤;Ni为当量应力σ当量 i的作用次数;m为应力-循环曲线的指数,取m=4;σmax为应力循环的峰值;σmin为应力循环的谷值。
可选地,步骤三中,所述结构在所述第一载荷谱下的等效应力为:
σ=(LB/GA)1/4
其中,σ为结构C在第一载荷谱A下的等效应力;GA为第一载荷谱A的当量过载损伤;LB为第二载荷谱B下结构C的当量应力损伤。
可选地,步骤四中,所述第一载荷谱下数据集包括裂纹长度a和飞行小时t。
可选地,步骤四中,对所述结构进行耐久性计算评估的方法为概率断裂力学方法。
发明至少存在以下有益技术效果:
本申请的基于应力损伤等效的耐久性计算方法,可以构建不同载荷谱下结构耐久性的关系,从而通过已有的数据集在新飞机结构设计早期进行耐久性计算评估。
附图说明
图1是本申请一个实施方式的基于应力损伤等效的耐久性计算方法流程图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1对本申请做进一步详细说明。
本申请提供了一种基于应力损伤等效的耐久性计算方法,方法包括:
步骤一:计算第一载荷谱的当量过载损伤;
步骤二:计算第二载荷谱下结构的当量应力损伤;
步骤三:根据当量过载损伤和当量应力损伤等效折算结构在第一载荷谱下的等效应力;
步骤四:根据第一载荷谱下数据集以及等效应力,对结构进行耐久性计算评估。
具体的,步骤一中,第一载荷谱A的当量过载损伤为GA
GA=∑(n当量 i)mNi
n当量 i=[2n幅值(n幅值+|n均值|)]1/2
n幅值=(nmax-nmin)/2
n均值=(nmax+nmin)/2
其中,GA为第一载荷谱A的当量过载损伤;Ni为n当量 i的作用次数;m为应力-循环曲线(S-N曲线)的指数,取m=4;nmax为过载循环的峰值;nmin为过载循环的谷值。
步骤二中,第二载荷谱B下结构的当量应力损伤为:
LB=∑(σ当量 i)mNi
σ当量 i=[2σ幅值幅值+|σ均值|)]1/2
σ幅值=(σmaxmin)/2
σ均值=(σmaxmin)/2
其中,LB为第二载荷谱B下结构C的当量应力损伤;Ni为当量应力σ当量 i的作用次数;m为应力-循环曲线的指数,取m=4;σmax为应力循环的峰值;σmin为应力循环的谷值。
步骤三中,结构C在第一载荷谱A下的等效应力σ为:
σ=(LB/GA)1/4
其中,σ为结构C在第一载荷谱A下的等效应力;GA为第一载荷谱A的当量过载损伤;LB为第二载荷谱B下结构C的当量应力损伤。
步骤四中,根据第一载荷谱下(裂纹长度a,飞行小时t)数据集以及等效应力,采用概率断裂力学方法对结构C进行耐久性计算评估
在本申请的一个实施方式中,第一载荷谱A的过载峰谷值如表1所示:
表1
根据步骤一中的公式计算可得n幅值、n均值及n当量 i,如表2所示:
表2
序号 n<sub>幅值</sub> n<sub>均值</sub> n<sup>当量</sup><sub>i</sub>
1 1.10 2.90 2.97
2 1.60 3.10 3.88
3 1.80 3.20 4.24
4 1.40 3.00 3.51
5 1.05 3.05 2.93
6 1.25 2.95 3.24
7 1.50 2.90 3.63
8 1.25 2.75 3.16
计算可得GA=∑(n当量 i)mNi=2623.391。
本实施例中,结构C耐久性关键部位分为3个应力区,各应力区在第二载荷谱B下的应力峰谷值如表3所示:
表3
根据步骤二中的公式计算可得σ幅值、σ均值及σ当量 i,如表4所示:
表4
计算可得各应力区的当量应力损伤LB,如表5所示。
表5
进一步,本实施例步骤三中,通过当量过载损伤和当量应力损伤等效折算结构C在第一载荷谱A下的等效应力σ,计算可得各应力区的等效应力,如表6所示。
表6
应力区 1 2 3
等效应力σ(MPa) 43.21 40.12 37.04
本实施例中,第一载荷谱A下(裂纹长度a,飞行小时t)数据集如表7所示:
表7
各应力区对应的细节数量如表8所示:
应力区 1 2 3
等效应力σ(MPa) 43.21 40.12 37.04
细节数量(个) 20 30 40
取可靠度R=0.95,经济修理极限ae=0.8mm,采用概率断裂力学方法即可计算得到在第二载荷谱B下结构C的耐久性评估结果,如表9所示。
表9
损伤度 使用时间t(小时)
1 1825
2 2295
3 2612
4 2863
5 3067
6 3245
7 3402
8 3545
9 3674
10 3794
本申请的基于应力损伤等效的耐久性计算方法,可以构建不同载荷谱下结构耐久性的关系,从而通过已有的(裂纹长度a,飞行小时t)数据集在新飞机结构设计早期进行耐久性计算评估。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。

Claims (6)

1.一种基于应力损伤等效的耐久性计算方法,其特征在于,包括:
步骤一:计算第一载荷谱的当量过载损伤;
步骤二:计算第二载荷谱下结构的当量应力损伤;
步骤三:根据所述当量过载损伤和所述当量应力损伤等效折算所述结构在所述第一载荷谱下的等效应力;
步骤四:根据所述第一载荷谱下数据集以及所述等效应力,对所述结构进行耐久性计算评估。
2.根据权利要求1所述的基于应力损伤等效的耐久性计算方法,其特征在于,步骤一中,所述计算第一载荷谱的当量过载损伤为:
GA=∑(n当量 i)mNi
n当量 i=[2n幅值(n幅值+|n均值|)]1/2
n幅值=(nmax-nmin)/2
n均值=(nmax+nmin)/2
其中,GA为第一载荷谱A的当量过载损伤;Ni为n当量 i的作用次数;m为应力-循环曲线的指数,取m=4;nmax为过载循环的峰值;nmin为过载循环的谷值。
3.根据权利要求2所述的基于应力损伤等效的耐久性计算方法,其特征在于,步骤二中,所述计算第二载荷谱下结构的当量应力损伤为:
LB=∑(σ当量 i)mNi
σ当量 i=[2σ幅值幅值+|σ均值|)]1/2
σ幅值=(σmaxmin)/2
σ均值=(σmaxmin)/2
其中,LB为第二载荷谱B下结构C的当量应力损伤;Ni为当量应力σ当量 i的作用次数;m为应力-循环曲线的指数,取m=4;σmax为应力循环的峰值;σmin为应力循环的谷值。
4.根据权利要求3所述的基于应力损伤等效的耐久性计算方法,其特征在于,步骤三中,所述结构在所述第一载荷谱下的等效应力为:
σ=(LB/GA)1/4
其中,σ为结构C在第一载荷谱A下的等效应力;GA为第一载荷谱A的当量过载损伤;LB为第二载荷谱B下结构C的当量应力损伤。
5.根据权利要求4所述的基于应力损伤等效的耐久性计算方法,其特征在于,步骤四中,所述第一载荷谱下数据集包括裂纹长度a和飞行小时t。
6.根据权利要求1所述的基于应力损伤等效的耐久性计算方法,其特征在于,步骤四中,对所述结构进行耐久性计算评估的方法为概率断裂力学方法。
CN201910797068.XA 2019-08-27 2019-08-27 一种基于应力损伤等效的耐久性计算方法 Active CN110489914B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910797068.XA CN110489914B (zh) 2019-08-27 2019-08-27 一种基于应力损伤等效的耐久性计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910797068.XA CN110489914B (zh) 2019-08-27 2019-08-27 一种基于应力损伤等效的耐久性计算方法

Publications (2)

Publication Number Publication Date
CN110489914A true CN110489914A (zh) 2019-11-22
CN110489914B CN110489914B (zh) 2023-01-17

Family

ID=68554495

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910797068.XA Active CN110489914B (zh) 2019-08-27 2019-08-27 一种基于应力损伤等效的耐久性计算方法

Country Status (1)

Country Link
CN (1) CN110489914B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111400923A (zh) * 2020-03-23 2020-07-10 南京航空航天大学 单向陶瓷基复合材料任意应力加卸载本构行为快速计算方法
CN112213090A (zh) * 2020-09-25 2021-01-12 中国直升机设计研究所 一种直升机动部件损伤容限简化谱编制方法
CN113704874A (zh) * 2021-08-03 2021-11-26 中国航空工业集团公司沈阳飞机设计研究所 一种复杂受载进气道结构设计载荷分析方法

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1983455A2 (de) * 2007-04-12 2008-10-22 Autoform Engineering Gmbh Grenzdehnungsanalyse
US20110288790A1 (en) * 2009-01-08 2011-11-24 Battelle Memorial Institute Path-Dependent Cycle Counting And Multi-Axial Fatigue Evaluation Of Engineering Structures
CN103258105A (zh) * 2012-02-17 2013-08-21 中国航空工业集团公司沈阳飞机设计研究所 一种飞机起落架寿命消耗计算方法
CN103984815A (zh) * 2014-05-12 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种多信息代表起落载荷谱编制方法
CN104792633A (zh) * 2015-04-17 2015-07-22 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种飞机机体裂纹扩展寿命预测方法
CN104809321A (zh) * 2013-12-31 2015-07-29 中国航空工业集团公司沈阳飞机设计研究所 一种基于细节疲劳额定值的高机动飞机寿命分析方法
JP2016161524A (ja) * 2015-03-05 2016-09-05 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
CN106529094A (zh) * 2016-11-30 2017-03-22 中国航空工业集团公司沈阳飞机设计研究所 一种基于典型飞行动作的飞机严重载荷谱编制方法
CN106557630A (zh) * 2016-11-21 2017-04-05 中国石油大学(华东) 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法
CN107742025A (zh) * 2017-10-16 2018-02-27 北京航空航天大学 一种预测飞机金属蒙皮冲击后疲劳寿命的方法
CN107944160A (zh) * 2017-11-29 2018-04-20 中国航空工业集团公司沈阳飞机设计研究所 一种用于发动机固定装置的疲劳载荷谱计算方法
CN108052717A (zh) * 2017-12-01 2018-05-18 中国航空工业集团公司沈阳飞机设计研究所 一种基于局部应力-应变法的疲劳寿命校准方法
US20190054573A1 (en) * 2017-08-18 2019-02-21 The Regents Of The University Of Michigan Unified Fatigue Life Evaluation Method For Welded Structures
CN109374449A (zh) * 2018-09-25 2019-02-22 南京航空航天大学 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法
CN109684697A (zh) * 2018-12-14 2019-04-26 中国航空工业集团公司西安飞机设计研究所 一种当量损伤模型的确定方法
CN109918854A (zh) * 2019-04-09 2019-06-21 国网江苏省电力有限公司电力科学研究院 一种复合绝缘横担耐久性能评估方法及装置

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1983455A2 (de) * 2007-04-12 2008-10-22 Autoform Engineering Gmbh Grenzdehnungsanalyse
US20110288790A1 (en) * 2009-01-08 2011-11-24 Battelle Memorial Institute Path-Dependent Cycle Counting And Multi-Axial Fatigue Evaluation Of Engineering Structures
CN103258105A (zh) * 2012-02-17 2013-08-21 中国航空工业集团公司沈阳飞机设计研究所 一种飞机起落架寿命消耗计算方法
CN104809321A (zh) * 2013-12-31 2015-07-29 中国航空工业集团公司沈阳飞机设计研究所 一种基于细节疲劳额定值的高机动飞机寿命分析方法
CN103984815A (zh) * 2014-05-12 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种多信息代表起落载荷谱编制方法
JP2016161524A (ja) * 2015-03-05 2016-09-05 株式会社神戸製鋼所 残留応力推定方法及び残留応力推定装置
CN104792633A (zh) * 2015-04-17 2015-07-22 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种飞机机体裂纹扩展寿命预测方法
CN106557630A (zh) * 2016-11-21 2017-04-05 中国石油大学(华东) 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法
CN106529094A (zh) * 2016-11-30 2017-03-22 中国航空工业集团公司沈阳飞机设计研究所 一种基于典型飞行动作的飞机严重载荷谱编制方法
US20190054573A1 (en) * 2017-08-18 2019-02-21 The Regents Of The University Of Michigan Unified Fatigue Life Evaluation Method For Welded Structures
CN107742025A (zh) * 2017-10-16 2018-02-27 北京航空航天大学 一种预测飞机金属蒙皮冲击后疲劳寿命的方法
CN107944160A (zh) * 2017-11-29 2018-04-20 中国航空工业集团公司沈阳飞机设计研究所 一种用于发动机固定装置的疲劳载荷谱计算方法
CN108052717A (zh) * 2017-12-01 2018-05-18 中国航空工业集团公司沈阳飞机设计研究所 一种基于局部应力-应变法的疲劳寿命校准方法
CN109374449A (zh) * 2018-09-25 2019-02-22 南京航空航天大学 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法
CN109684697A (zh) * 2018-12-14 2019-04-26 中国航空工业集团公司西安飞机设计研究所 一种当量损伤模型的确定方法
CN109918854A (zh) * 2019-04-09 2019-06-21 国网江苏省电力有限公司电力科学研究院 一种复合绝缘横担耐久性能评估方法及装置

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LIU CHENCHEN等: "Analysis for In-situ Corrosion-Multiaxial Fatigue Behaviors of 7075 Aerospace Aluminum Alloys", 《CHINA MECHANICAL ENGINEERING》 *
LOK, S. KRISHNA等: "Prescience Life of Landing Gear using Multiaxial Fatigue Numerical Analysis", 《STRUCTURAL INTEGRITY》 *
吴德锋等: "先进飞机典型结构关键部位损伤容限分析研究", 《飞机设计》 *
孟宪红等: "谱载荷下三维构件疲劳分析的损伤力学方法", 《北京航空航天大学学报》 *
王长江: "民机载荷谱分散性及其等效方法研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *
董彦民等: "变更载荷谱情况下的耐久性分析技术与应用", 《北京航空航天大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111400923A (zh) * 2020-03-23 2020-07-10 南京航空航天大学 单向陶瓷基复合材料任意应力加卸载本构行为快速计算方法
CN111400923B (zh) * 2020-03-23 2021-10-19 南京航空航天大学 单向陶瓷基复合材料任意应力加卸载本构行为快速计算方法
CN112213090A (zh) * 2020-09-25 2021-01-12 中国直升机设计研究所 一种直升机动部件损伤容限简化谱编制方法
CN112213090B (zh) * 2020-09-25 2022-11-18 中国直升机设计研究所 一种直升机动部件损伤容限简化谱编制方法
CN113704874A (zh) * 2021-08-03 2021-11-26 中国航空工业集团公司沈阳飞机设计研究所 一种复杂受载进气道结构设计载荷分析方法
CN113704874B (zh) * 2021-08-03 2024-02-02 中国航空工业集团公司沈阳飞机设计研究所 一种复杂受载进气道结构设计载荷分析方法

Also Published As

Publication number Publication date
CN110489914B (zh) 2023-01-17

Similar Documents

Publication Publication Date Title
CN110489914A (zh) 一种基于应力损伤等效的耐久性计算方法
Ahmed The perils of unearned foreign income: Aid, remittances, and government survival
EP2128025A3 (en) Anti-collision system between aircraft boarding bridges
CN106779619B (zh) 一种完善业务审批的审核加签方法及系统
CN103926152A (zh) 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法
CN103344514A (zh) 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法
CN103226714A (zh) 基于较大编码系数强化的稀疏编码方法
CN103914623B (zh) 用于飞机结构腐蚀疲劳关键件延寿的寿命包线扩展方法
Crighton et al. Influenza and pneumonia hospitalizations in Ontario: a time-series analysis
CN103544402A (zh) 一种疲劳开裂结构等效分析谱的构造方法
CN103674588B (zh) 一种飞机金属壁板结构疲劳临界部位筛选方法
CN106599405A (zh) 一种计算主减速器与机身任意连接点位置实际载荷的方法
CN109507040A (zh) 一种蜂窝夹层结构面板压缩应力评估方法
CN208395690U (zh) 一种定位方便的桥台背墙裂缝加固结构
CN109086471B (zh) 一种嵌入式轨道简支梁桥上无缝线路的简便设计方法
CN109376476A (zh) 装配应力工程计算方法
CN206369581U (zh) 一种使用拉力传感器的动态称重系统
CN109625312A (zh) 飞机结构蒙皮损伤修复装置
CN204137048U (zh) 档案专用车
CN208166313U (zh) 电梯轿底结构
CN112478130B (zh) 一种加强框及其座舱盖结构
Majery Work ability of a population of 40+ in Luxembourg
Tsukimori et al. Study on behaviours of multi-ply bellows subjected to pressure and displacement loads
CN205675829U (zh) 长钢轨卸车用一体化锁轨设备
CN213479940U (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