CN113838535B - 一种统一流变应力模型及计算方法 - Google Patents
一种统一流变应力模型及计算方法 Download PDFInfo
- Publication number
- CN113838535B CN113838535B CN202110951010.3A CN202110951010A CN113838535B CN 113838535 B CN113838535 B CN 113838535B CN 202110951010 A CN202110951010 A CN 202110951010A CN 113838535 B CN113838535 B CN 113838535B
- Authority
- CN
- China
- Prior art keywords
- stress
- strain
- dynamic recrystallization
- rheological
- model
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 14
- 238000001953 recrystallisation Methods 0.000 claims abstract description 44
- 239000000463 material Substances 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000005482 strain hardening Methods 0.000 claims abstract description 14
- 238000000611 regression analysis Methods 0.000 claims abstract description 4
- 230000006835 compression Effects 0.000 claims abstract description 3
- 238000007906 compression Methods 0.000 claims abstract description 3
- 238000002474 experimental method Methods 0.000 claims abstract description 3
- 238000004514 thermodynamic simulation Methods 0.000 claims abstract description 3
- 238000013211 curve analysis Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims 1
- 238000005096 rolling process Methods 0.000 description 11
- 229910000831 Steel Inorganic materials 0.000 description 6
- 239000010959 steel Substances 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910000851 Alloy steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012669 compression test Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000005098 hot rolling Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000003303 reheating Methods 0.000 description 1
- 238000003856 thermoforming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种统一流变应力模型及计算方法，通过单道次等温压缩热力模拟实验获得材料在不同温度和应变速率下的流变应力曲线；将曲线分为两类：未发生动态再结晶型曲线，该过程材料只发生加工硬化和回复软化；发生动态再结晶型，此时流变应力会出现峰值，并能够确定出发生动态再结晶的临界点，在该点之后继续变形即发生动态再结晶；获得加工硬化率θ与其对应的应力σ之间的关系曲线；建立变形温度、应变速率与特征点参数的关系，先找出特征点对应的应力和应变与参数Z的关系；采用数据分析软件进行回归分析，即可得到相应的流变应力模型。优点是：简单易行，能够快速地计算出流变应力。
Description
技术领域
本发明属于材料热加工领域，涉及一种统一流变应力模型及计算方法。
背景技术
在材料热成形过程中，流变应力作为变形温度、变形量和应变速率的函数，是压力加工工艺及设备的力、能量参数设计的重要参数之一。
轧钢生产在过去依赖经验进行操作，到近年来发展到生产流程的半自动化以及依靠计算机自动控制，流变应力模型作为加工过程中精确塑性成形控制的最基本数学模型，已成为从事轧钢技术工作者技术创新的必要工具，同时也是进行大规模高速连续化生产的保证。
因此，建立材料流变应力模型对于顺利实现轧制过程的计算机模拟以及精确控制，优化热轧工艺过程，改善产品的质量等方面都具有重要作用。
目前，轧制力预测普遍采用物理模拟实验手段，获得材料在不同变形条件下的流变应力曲线，基于曲线分析，得到材料在不同变形条件下的参量，进而基于经验模型建立起各个参量的关系，最终得到流变应力模型。
申请号为201210233146.1的专利一种用于轧机设计的材料变形抗力统一模型的构造方法，提出了一种轧机设计过程具有外推稳定性的材料变形抗力统一模型，该模型将材料的变形抗力表述为变形温度，应变速率以及变形量的非线性关系，从而得到材料变形抗力的统一模型。申请号为201310157302.5的专利一种板坯轧制过程中变形抗力的预测方法，将变形抗力分解为四个独立的参量，即基于参考屈服应力、变形温度影响系数、变形速率影响系数和变形程度影响系数这四个参量对变形抗力的影响关系，得到描述变形抗力的通用公式，以获得具有较高精度的结果。两者共同点在于所得的变形抗力模型都具有非线性特点，不同之处在于前者考虑了参数之间的交互作用，而后者将各个参量分立开来考虑。然而，两者都没有考虑材料在热加工过程中经历的组织变化阶段，且材料内部组织的变化对其变形抗力会产生直接影响，这必然会造成计算结果的偏差。申请号为201610802819.9的专利一种构建材料变形抗力模型的方法，考虑了材料再热加工过程中经历的组织变化，将模型进行了分段处理，但分段的模型不一致，建立和求解都有一定难度。
发明内容
为克服现有技术的不足，本发明的目的是提供一种统一流变应力模型及计算方法，该模型采用统一的形式表示材料在热加工过程中流变应力和应变关系，同时考虑了材料在热加工过程中发生的加工硬化、回复、动态再结晶等一系列变化对流变应力的影响，简化了流变应力模型结构及其计算。
为实现上述目的，本发明通过以下技术方案实现：
一种统一流变应力模型，包括：
确定统一流变应力模型为：
式(1)中，σ、ε分别为材料变形过程中的流变应力和应变；σT为特征应力；εT为特征应变；aT，bT，cT，dT为待定常数；
模型按照是否发生动态再结晶进行分类：一类为未发生动态再结晶情况；另一类为发生动态再结晶情况；
对于整个变形过程中未发生动态再结晶情况，模型表示为：
式(2)中，σ0为对应的屈服应力；ε0为对应的屈服应变；a0，b0，c0，d0为待定常数；
对于整个变形过程中发生动态再结晶情况，模型按照临界点进行分段，表示为：
式(3)中，σc为发生动态再结晶的临界应力，εc为发生动态再结晶的临界应变，a1，b1，c1，d1为待定常数；变形的应变范围为ε＜εc；
发生动态再结晶且变形量达到峰值时的流变应力模型表示为：
式(4)中，σp为流变应力中的峰值应力，εP为峰值应力对应的应变，a2，b2，c2，d2为待定常数；变形的应变范围为εc＜ε＜εp；
发生动态再结晶且变形量大于峰值时的流变应力模型表示为：
式(5)中，变形的应变范围为εp＜ε。
一种统一流变应力模型的计算方法，包括以下步骤：
1)获取基础数据
通过单道次等温压缩热力模拟实验获得材料在不同温度和应变速率下的流变应力曲线；
2)曲线分类
基于步骤1)中所得流变应力曲线分析，将曲线分为两类：
未发生动态再结晶型曲线，该过程材料只发生加工硬化和回复软化；
发生动态再结晶型，此时流变应力会出现峰值，并能够确定出发生动态再结晶的临界点，在该点之后继续变形即发生动态再结晶；
3)特征点获取
对步骤1)中的流变应力曲线进行分析，获得加工硬化率θ与其对应的应力σ之间的关系曲线；曲线中无拐点的则找出屈服点；曲线中有拐点的，则将拐点确定为发生动态再结晶的临界点，临界点对应临界应力σc，根据临界应力找到临界应变εc；流变应力曲线中有峰值点的，则找出峰值应力σp及其所对应的应变εp；
4)特征点对应的应力和应变的计算
建立变形温度、应变速率与特征点参数的关系，先找出特征点对应的应力和应变与参数Z的关系，关系由如下计算给出：
σ0＝aZb (6)
ε0＝cZd (7)
σc＝eZf (8)
εc＝gZh (9)
σp＝iZj (10)
εp＝kZl (11)
式(6)～(11)中，a,b,…,l为待定系数；
将式(6)～(11)分别与式(12)联立可以得到各个特征值与变形速率和温度的关系式；
式(12)中，Z为Z-hollomen参数，为应变速率，T为变形温度，Q为表观动态再结晶激活能，R为气体常数；
5)流变应力模型
将步骤5)计算得到的特征点对应的应力和应变对应地代入到式(2)～(5)中，再采用数据分析软件进行回归分析，即可得到相应的流变应力模型。
与现有技术相比，本发明的有益效果是：
本发明的模型采用统一的形式表示出了材料在热加工过程中流变应力和应变关系，并充分考虑了材料内部在轧制过程中经历的加工硬化、回复、动态再结晶等一系列的变化对变形抗力的影响，根据材料不同的组织变化进行分类，材料不同的组织变化阶段进行分段，模型更科学、准确，所提供的计算方法也简单易行，能够快速地计算出流变应力，为准确预测轧制力提供保障。
附图说明
图1是确定屈服应力应变的示意图。
图1中，1-确定的屈服点。
图2是加工硬化率与应力关系曲线图。
图2中，变形温度为900℃，应变速率为3s-1。
图3是加工硬化率与应力关系曲线图。
图3中，变形温度为1050℃，应变速率为0.1s-1。
图4是实验钢典型的流变应力曲线图。
图4中，变形温度为900℃，应变速率为3s-1。
图5是实验钢典型的流变应力曲线图。
图5中，变形温度为1050℃，应变速率为0.1s-1。
具体实施方式
下面结合说明书附图对本发明进行详细地描述，但是应该指出本发明的实施不限于以下的实施方式。
将本发明应用于一种合金钢，具体通过下列步骤实施：
1)获取基础数据
通过热力模拟试验机，对实验钢进行了单道次压缩实验，变形温度为900～1050℃，应变速率为0.1s-1，3s-1，总真应变量为70％，以获得给定变形温度和变形速率条件下的流变应力曲线，进而根据该曲线的不同特征进行分类，不同类型曲线，通过不同模型表示。
2)曲线分类
基于步骤1)中所得流变应力曲线分析，将曲线分为两类：
a一种为未发生动态再结晶型曲线，该过程材料只发生加工硬化和回复软化，变形温度为900，950℃条件下获得的实验曲线没有典型峰值，说明未发生动态再结晶，与该情况对应。
b另一种为发生动态再结晶型，此时流变应力会出现峰值，并能够确定出发生动态再结晶的临界点(在该点之后继续变形即发生动态再结晶)，变形温度为1000，1050℃条件下获得的实验曲线均出现典型峰值，说明发生了动态再结晶，与该情况对应。
3)特征点获取
对步骤1)中的流变应力曲线进行分析，获得加工硬化率θ与其对应的应力σ之间的关系曲线，曲线中无拐点的则找出屈服点，如图2所示，为变形温度为900℃，应变速率为3s-1条件下获得的曲线，曲线无拐点；曲线中有拐点的，则将拐点确定为发生动态再结晶的临界点，临界点对应临界应力σc，根据临界应力找到临界应变εc，图3给出了变形温度为1050℃，应变速率为0.1s-1条件下获得的曲线，曲线出现拐点。流变应力曲线中有峰值点的，则找出峰值应力σp及其所对应的应变εp。由步骤2)可知，变形温度为900，950℃条件下获得的实验曲线没有典型峰值，加工硬化率θ与其对应的应力σ之间的关系曲线无拐点，则由下一步骤找出屈服点，变形温度为1000℃，1050℃条件下获得的实验曲线有典型峰值，加工硬化率θ与其对应的应力σ之间的关系曲线有拐点，则找出相应的特征值点，见表1。
表1
4)屈服点的确定
对变形温度为900，950℃条件下获得的实验曲线进行分析，屈服点对应的屈服应力σ0与应变ε0按照图1确定。由应变为0.002处引一条与应力应变曲线初始段相平行的直线,该直线与应力应变曲线相交于一点,该点确定为屈服点，利用该点的坐标可确定出屈服应力(纵坐标)，不同条件下获取的屈服应力和应变，见表2。
表2
5)特征点参数(特征点对应的应力和应变)的计算
建立变形条件(变形温度、应变速率)与特征点参数的关系，先按照式(6)～(11)找出特征点对应的应力和应变与参数Z的关系，再分别与式(12)联立可以得到各个特征值与变形速率和温度的关系式；
6)流变应力模型
将步骤5)计算得到的特征点参数(特征点对应的应力和应变)对应地代入到式(2)～(5)中，再采用Origin数据分析软件进行回归分析，即可得到相应的流变应力模型，发生动态再结晶情况模型，见表3，未发生动态再结晶情况模型见表4。
表3
表4
将式(6)～(11)中确定的特征值根据不同的变形条件，分别代入表3、表4的流变应力模型中，即可以得到变形温度、应变速率与应变量对流变应力的影响关系，即流变应力模型。图4和图5中为实验钢典型的流变应力曲线，实线为模型计算曲线，数据点为实验曲线上的实际数据点，从图中可知，所建立的数据模型具有很高的计算精度，能够很好地预测轧制生产线上的轧机的轧制力。
Claims (1)
1.一种统一流变应力模型，其特征在于，包括：
确定统一流变应力模型为：
式(1)中，σ、ε分别为材料变形过程中的流变应力和应变；σT为特征应力；εT为特征应变；aT，bT，cT，dT为待定常数；
模型按照是否发生动态再结晶进行分类：一类为未发生动态再结晶情况；另一类为发生动态再结晶情况；
对于整个变形过程中未发生动态再结晶情况，模型表示为：
式(2)中，σ0为对应的屈服应力；ε0为对应的屈服应变；a0，b0，c0，d0为待定常数；
对于整个变形过程中发生动态再结晶情况，模型按照临界点进行分段，表示为：
式(3)中，σc为发生动态再结晶的临界应力，εc为发生动态再结晶的临界应变，a1，b1，c1，d1为待定常数；变形的应变范围为ε＜εc；
发生动态再结晶且变形量达到峰值时的流变应力模型表示为：
式(4)中，σp为流变应力中的峰值应力，εP为峰值应力对应的应变，a2，b2，c2，d2为待定常数；变形的应变范围为εc＜ε＜εp；
发生动态再结晶且变形量大于峰值时的流变应力模型表示为：
式(5)中，变形的应变范围为εp＜ε；
一种统一流变应力模型的计算方法，包括以下步骤：
1)获取基础数据
通过单道次等温压缩热力模拟实验获得材料在不同温度和应变速率下的流变应力曲线；
2)曲线分类
基于步骤1)中所得流变应力曲线分析，将曲线分为两类：
未发生动态再结晶型曲线，该过程材料只发生加工硬化和回复软化；
发生动态再结晶型，此时流变应力会出现峰值，并能够确定出发生动态再结晶的临界点，在该点之后继续变形即发生动态再结晶；
3)特征点获取
对步骤1)中的流变应力曲线进行分析，获得加工硬化率θ与其对应的应力σ之间的关系曲线；曲线中无拐点的则找出屈服点，屈服点对应的屈服应力σ0，根据屈服应力找到屈服应变ε0；曲线中有拐点的，则将拐点确定为发生动态再结晶的临界点，临界点对应临界应力σc，根据临界应力找到临界应变εc；流变应力曲线中有峰值点的，则找出峰值应力σp及其所对应的应变εp；
4)特征点对应的应力和应变的计算
建立变形温度、应变速率与特征点参数的关系，先找出特征点对应的应力和应变与参数Z的关系，关系由如下计算给出：
σ0＝aZb (6)
ε0＝cZd (7)
σc＝eZf (8)
εc＝gZh (9)
σp＝iZj (10)
εp＝kZl (11)
式(6)～(11)中，a,b,…,l为待定系数；
将式(6)～(11)分别与式(12)联立可以得到各个特征值与应变速率和温度的关系式；
式(12)中，Z为Z-hollomen参数，为应变速率，T为变形温度，Q为表观动态再结晶激活能，R为气体常数；
5)流变应力模型
将步骤4)计算得到的特征点对应的应力和应变对应地代入到式(2)～(5)中，再采用数据分析软件进行回归分析，即可得到相应的流变应力模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110951010.3A CN113838535B (zh) | 2021-08-18 | 2021-08-18 | 一种统一流变应力模型及计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110951010.3A CN113838535B (zh) | 2021-08-18 | 2021-08-18 | 一种统一流变应力模型及计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113838535A CN113838535A (zh) | 2021-12-24 |
| CN113838535B true CN113838535B (zh) | 2023-09-26 |
Family
ID=78960825
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110951010.3A Active CN113838535B (zh) | 2021-08-18 | 2021-08-18 | 一种统一流变应力模型及计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113838535B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651620A (zh) * | 2016-03-09 | 2016-06-08 | 中南大学 | 一种评估金属材料热加工性能的方法 |
| CN107818184A (zh) * | 2016-09-06 | 2018-03-20 | 鞍钢股份有限公司 | 一种构建材料变形抗力模型的方法 |
| CN108563880A (zh) * | 2018-04-20 | 2018-09-21 | 中国兵器工业第五九研究所 | 一种计及应变硬化和软化的金属材料本构模型建立方法 |
| JP2019155456A (ja) * | 2018-03-16 | 2019-09-19 | 日立金属株式会社 | 熱間鍛造条件の設定方法及び鍛造製品の製造方法 |
-
2021
- 2021-08-18 CN CN202110951010.3A patent/CN113838535B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651620A (zh) * | 2016-03-09 | 2016-06-08 | 中南大学 | 一种评估金属材料热加工性能的方法 |
| CN107818184A (zh) * | 2016-09-06 | 2018-03-20 | 鞍钢股份有限公司 | 一种构建材料变形抗力模型的方法 |
| JP2019155456A (ja) * | 2018-03-16 | 2019-09-19 | 日立金属株式会社 | 熱間鍛造条件の設定方法及び鍛造製品の製造方法 |
| CN108563880A (zh) * | 2018-04-20 | 2018-09-21 | 中国兵器工业第五九研究所 | 一种计及应变硬化和软化的金属材料本构模型建立方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113838535A (zh) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105631183B (zh) | 一种建立tc6钛合金锻造成形微观组织预测模型的方法 | |
| Qian et al. | 3D coupled macro–microscopic finite element modelling and simulation for combined blank-forging and rolling process of alloy steel large ring | |
| Lv et al. | 3D FEM simulation of the multi-stage forging process of a gas turbine compressor blade | |
| Ji et al. | Numerical analysis and experiment on cross wedge rolling and forging for engine valves | |
| CN102463268B (zh) | 一种锌合金棒材挤压的工艺设计与优化方法 | |
| CN111753453B (zh) | 一种高强度钢模锻件成形工艺的高精度模拟方法 | |
| CN107122502B (zh) | 一种优化合金挤压工艺的方法 | |
| CN112948986B (zh) | 结合有限元数值模拟程序的钛合金锻造工艺参数优化方法 | |
| Kukhar et al. | Experimental research of distribution of strains and stresses in work-piece at different modes of stretch-forging with rotation in combined dies | |
| CN111767665B (zh) | 一种高温合金盘锻件制坯用模具的型腔设计方法 | |
| CN113591341A (zh) | 一种基于数值模拟的钛合金锻造工艺优化方法 | |
| CN107423469B (zh) | 一种06Cr19Ni9NbN钢锻透的判定方法 | |
| Sun et al. | Thermo-mechanical coupled analysis of hot ring rolling process | |
| CN113838535B (zh) | 一种统一流变应力模型及计算方法 | |
| CN111889601B (zh) | 增大难变形合金模锻件表层变形的波浪形坯料设计方法 | |
| CN107818184B (zh) | 一种构建材料变形抗力模型的方法 | |
| Samołyk | Numerical investigation of producing a Ti6Al4V alloy jaw coupling sleeve-disk by orbital forging | |
| CN113343534B (zh) | 一种提高钛合金高筋锻件准β锻变形均匀性的方法 | |
| CN109033590B (zh) | 一种基于模锻模具磨损的工艺参数优化方法 | |
| CN107526877B (zh) | 一种消减残余应力的模拟方法 | |
| CN114283900B (zh) | 针对近β钛合金低倍粗晶组织分布的预测与调控方法 | |
| Liu et al. | Wear Analysis of Hot Extrusion Punch for Large Fan Shaft Based on Archard Wear Theory | |
| CN109909413B (zh) | 一种基于热加工图的锻模速度曲线迭代优化方法 | |
| Shahriari et al. | Optimal closed die finish forgings for nimonic80-A alloy using FEM method | |
| Wang et al. | Grain size control method for the nozzles of AP1000 primary coolant pipes |
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 |