CN109470734A - 一种计算马氏体组织转变量及相变点的方法 - Google Patents
一种计算马氏体组织转变量及相变点的方法 Download PDFInfo
- Publication number
- CN109470734A CN109470734A CN201811156239.2A CN201811156239A CN109470734A CN 109470734 A CN109470734 A CN 109470734A CN 201811156239 A CN201811156239 A CN 201811156239A CN 109470734 A CN109470734 A CN 109470734A
- Authority
- CN
- China
- Prior art keywords
- temperature
- transformation
- swell increment
- straight line
- curve
- 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
Links
- 230000009466 transformation Effects 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 28
- 229910000734 martensite Inorganic materials 0.000 title claims abstract description 23
- 238000001816 cooling Methods 0.000 claims abstract description 33
- 230000008961 swelling Effects 0.000 claims abstract description 7
- 238000000926 separation method Methods 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 241001269238 Data Species 0.000 claims description 6
- 229910052799 carbon Inorganic materials 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 abstract description 2
- 239000002184 metal Substances 0.000 abstract description 2
- 229910052751 metal Inorganic materials 0.000 abstract description 2
- 239000000284 extract Substances 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 7
- 230000007423 decrease Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910000746 Structural steel Inorganic materials 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 229910001563 bainite Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910001562 pearlite Inorganic materials 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/16—Investigating or analyzing materials by the use of thermal means by investigating thermal coefficient of expansion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/02—Investigating or analyzing materials by the use of thermal means by investigating changes of state or changes of phase; by investigating sintering
- G01N25/12—Investigating or analyzing materials by the use of thermal means by investigating changes of state or changes of phase; by investigating sintering of critical point; of other phase change
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开一种计算马氏体组织转变量及相变点的方法，属于金属相变的表征计算领域。本发明所述方法首先提取冷却部分的膨胀量‑温度数据，然后对冷却图像的数据在特定温度间隔下求平均值，并进行平滑处理，以最大膨胀量与最小膨胀量之差与该膨胀量下所对应的温度之差的比值作为拟合直线和拟合曲线的分离评判标准，以相变点以外的80±1℃温度区间步长内的数据进行迭代，采用最小二乘法回归拟合得到两条直线，最后根据杠杆原理求每个温度下对应的转变量。本发明提供的方法可以精确地计算出每个温度对应的转变量，克服了以金相图片来断定相变点及组织转变量的不准确性，大大节省了成本和时间。
Description
技术领域
本发明涉及一种计算马氏体组织转变量及相变点的方法，属于金属相变的表征计算领域。
背景技术
以往测定组织转变量及转变起止点一般根据经验公式或者金相图片，但是组织转变伴随着物性变化和力学性能变化，同时还要相变潜热释放，对于不同的材料，部分经验公式需要加以修正才可以应用，而根据金相图片判断组织转变量及相变点浪费大量的人力、物力、财力，所得结果往往也不够精确，具有较大的局限性。找到精确确定每个温度组织转变量及相变点的方法变得日益重要。
发明内容
本发明的目的在于提供一种计算马氏体组织转变量及相变点的方法，精确地确定每个温度下对应的转变量及转变起止点，具体包括以下步骤：
（1）将待测材料进行膨胀试验，最终得到相变膨胀曲线，将加热阶段与冷却阶段的膨胀量曲线分离，选取冷却部分的膨胀量曲线；
（2）以1~3℃温度区间划分为若干等距区间，将每个区间中的温度的平均值作为横坐标，将每个区间中的膨胀量的平均值作为纵坐标，对降温数据进行平均处理；
（3）对膨胀量-温度数据进行进行平滑处理，如果平滑效果不好，进行第二次平滑处理；
（4）将步骤（3）平滑处理后的高温段数据按照温度由高到低依次命名为（T1，D1）、（T2，D2）、（T3，D3）……（Tn，Dn），对（T1，D1）、（T2，D2）……（Tn，Dn）对应的n组数据进行线性回归得到拟合直线，其中n≥2，T表示温度，D表示膨胀量；若|D曲线-D直线|≥D总/T温差，则该点即为拟合直线与膨胀量-温度曲线的分离点，其中，D总为最大膨胀量与最小膨胀量之差，T温差为最大膨胀量对应的温度与最小膨胀量对应的温度之差；若|Tm-T1|=80±1℃，则Tm即为相变开始点，其中Tm为拟合直线与膨胀量-温度曲线的分离点对应的温度；此时（T1，D1）、（T2，D2）……（Tm，Dm）的拟合直线即为高温段切线；D=k×T+d，式中，k为高温段切线斜率，d为截距，D为膨胀量，T为温度；对于低温相变结束点和低温段切线，采用同样的方法处理；
（5）利用杠杆定律可计算相变开始点和相变结束点对应温度区间内的马氏体组织转变量，f=Xa/（Xa+Xb），式中：f为马氏体组织转变量，Xa为高温段数据拟合直线与平滑后膨胀曲线的竖直距离，Xb为低温段数据拟合直线与平滑后膨胀曲线的竖直距离。
本发明所述待测材料为发生相变的材料。
本发明所述步骤（3）中以30个膨胀量-温度数据为一个平滑组进行平滑，第二次平滑以20-40个膨胀量-温度数据为一个平滑组进行平滑。
本发明的有益效果：
（1）本发明可以精确地计算出不同冷速下对应的马氏体转变量及转变起止点，不受冷却速率的影响，适用于发生相变的任何材料（例如马氏体转变，贝氏体、珠光体的转变）；
（2）本发明可以有效保证实测数据的真实性和准确性，并能大大降低数据的繁杂性，降低人为因素的影响；
（3）本发明提供的方法操作方便，精确度较好，工艺简单，节约成本和人力。
附图说明
图1为40℃/s冷速下实验得到的膨胀量-温度曲线；
图2为127℃/s冷速下实验得到的膨胀量-温度曲线；
图3为实施例1在40℃/s的冷却速度下拟合计算得到的转变量图像；
图4为实施例1在127℃/s的冷却速度下拟合计算得到的转变量图像；
图5为实施例2在40℃/s的冷却速度下拟合计算得到的转变量图像。
具体实施方式
下面结合附图和具体实施例对本发明做进一步描述，但本发明的保护范围并不限于所述内容。
实施例1
一种计算马氏体组织转变量及相变点的方法，具体包括以下步骤：
（1）将待测材料进行膨胀试验，最终得到相变膨胀曲线，将加热阶段与冷却阶段的膨胀量曲线分离，选取冷却部分的膨胀量曲线；
（2）以1℃温度区间划分为若干等距区间，将每个区间中的温度的平均值作为横坐标，将每个区间中的膨胀量的平均值作为纵坐标，对降温数据进行平均处理；
（3）以30个膨胀量-温度数据为一个平滑组进行平滑处理，再以20个膨胀量-温度数据为一个平滑组进行平滑；
（4）将步骤（3）平滑处理后的高温段数据按照温度由高到低依次命名为（T1，D1）、（T2，D2）、（T3，D3）……（Tn，Dn），对（T1，D1）、（T2，D2）……（Tn，Dn）对应的n（n≥2，整数递增）组数据进行线性回归得到拟合直线，其中T表示温度，D表示膨胀量；若|D曲线-D直线|≥D总/T温差，则该点即为拟合直线与膨胀量-温度曲线的分离点，其中，D总为最大膨胀量与最小膨胀量之差，T温差为最大膨胀量对应的温度与最小膨胀量对应的温度之差；若|Tm-T1|=80±1℃，则Tm即为相变开始点，其中Tm为拟合直线与膨胀量-温度曲线的分离点对应的温度；此时（T1，D1）、（T2，D2）……（Tm，Dm）的拟合直线即为高温段切线，若|Tm-T1|≥81℃或者|Tm-T1|＜79℃，则缩小或扩大拟合直线所用的数据区域，采用前述方法重新进行拟合计算，直至|Tm-T1|=80±1℃为止，如此往复，得到相变开始点和高温段切线；D=k×T+d，式中，k为高温段切线斜率，d为截距，D为膨胀量，T为温度；对于低温相变结束点和低温段切线，采用同样的方法处理。
本实施例以Q235普通碳素结构钢为待测试样，进行膨胀实验，实验过程中冷速分别为为40℃/s和127℃/s。
当冷速为40℃/s的时候从膨胀量-温度曲线可以计算出：D总/T温差=0.02/190=0.00011，通过|D曲线-D直线|≥0.00011mm可以得到482.1℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到高温段切线D=0.000122×T-0.040。采用同样的方法处理，低温段186.5℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到低温段切线分别为D=0.0000579×T+0.0129。
当冷速为127℃/s的时候从膨胀量-温度曲线可以计算出：D总/T温差=0.17/239=0.00075mm，通过|D曲线-D直线|≥0.00075mm可以得到463.9℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到高温段切线D=0.000123×T-0.049。采用同样的方法处理，低温段172.3℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到低温段切线分别为D=0.0000553×T+0.00359。
（5）利用杠杆定律可计算相变开始点和相变结束点对应温度区间内的马氏体组织转变量，f=Xa/（Xa+Xb），式中：f为马氏体组织转变量，Xa为高温段数据拟合直线与平滑后膨胀曲线的竖直距离，Xb为低温段数据拟合直线与平滑后膨胀曲线的竖直距离。
本实施例中当冷速为40℃/s的时候，图3中虚线即对应马氏体组织转变量与温度的对应关系；当冷速为127℃/s时，如图4所示。
表1为不同冷速下计算得到的相变开始点和相变结束点与利用差示扫描量热仪实验得到的温度点的对比，由表可以看出，本实施例所述方法计算得到的相变点与实际实验得到的相变点比较接近，想变点伴随着冷却速度的下降而下降，。
表1
实施例2
一种计算马氏体组织转变量及相变点的方法，具体包括以下步骤：
（1）将待测材料进行膨胀试验，最终得到相变膨胀曲线，将加热阶段与冷却阶段的膨胀量曲线分离，选取冷却部分的膨胀量曲线；
（2）以3℃温度区间划分为若干等距区间，将每个区间中的温度的平均值作为横坐标，将每个区间中的膨胀量的平均值作为纵坐标，对降温数据进行平均处理；
（3）以30个膨胀量-温度数据为一个平滑组进行平滑处理，再以40个膨胀量-温度数据为一个平滑组进行平滑；
（4）将步骤（3）平滑处理后的高温段数据按照温度由高到低依次命名为（T1，D1）、（T2，D2）、（T3，D3）……（Tn，Dn），对（T1，D1）、（T2，D2）……（Tn，Dn）对应的n（n≥2，整数递增）组数据进行线性回归得到拟合直线，其中T表示温度，D表示膨胀量；若|D曲线-D直线|≥D总/T温差，则该点即为拟合直线与膨胀量-温度曲线的分离点，其中，D总为最大膨胀量与最小膨胀量之差，T温差为最大膨胀量对应的温度与最小膨胀量对应的温度之差；若|Tm-T1|=80±1℃，则Tm即为相变开始点，其中Tm为拟合直线与膨胀量-温度曲线的分离点对应的温度；此时（T1，D1）、（T2，D2）……（Tm，Dm）的拟合直线即为高温段切线，若|Tm-T1|≥81℃或者|Tm-T1|＜79℃，则缩小或扩大拟合直线所用的数据区域，采用前述方法重新进行拟合计算，直至|Tm-T1|=80±1℃为止，如此往复，得到相变开始点和高温段切线；D=k×T+d，式中， k为高温段切线斜率， d为截距，D为膨胀量，T为温度；对于低温相变结束点和低温段切线，采用同样的方法处理。
本实施例以Q235普通碳素结构钢为待测试样，进行膨胀实验，实验过程中冷速40℃/s。
当冷速为40℃/s的时候从膨胀量-温度曲线可以计算出：D总/T温差=0.035/251=0.00014，通过|D曲线-D直线|≥0.00014mm可以得到483.2℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到高温段切线D=0.000123×T-0.040。用同样的方法处理，低温段177.2℃对应的点即为拟合直线与膨胀量-温度曲线的分离点，拟合得到低温段切线分别为D=0.0000583×T+0.0131。
（5）利用杠杆定律可计算相变开始点和相变结束点对应温度区间内的马氏体组织转变量，f=Xa/（Xa+Xb），式中：f为马氏体组织转变量，Xa为高温段数据拟合直线与平滑后膨胀曲线的竖直距离，Xb为低温段数据拟合直线与平滑后膨胀曲线的竖直距离。
本实施例中当冷速为40℃/s的时候，图5中虚线即对应马氏体组织转变量与温度的对应关系。
表1为不同冷速下计算得到的相变开始点和相变结束点与利用差示扫描量热仪实验得到的温度点的对比，由表可以看出，本实施例所述方法计算得到的相变点与实际实验得到的相变点比较接近。
表2
实验和实施例2计算得出的的马氏体Ms和Mf点差值很小，且都伴随着冷却速度的下降而下降，说明该计算方法适用性较好，实施例2在40℃/s冷速下组织转变量如图5所示。
Claims (3)
1.一种计算马氏体组织转变量及相变点的方法，具体步骤如下：
（1）将待测材料进行膨胀试验，最终得到相变膨胀曲线，将加热阶段与冷却阶段的膨胀量曲线分离，选取冷却部分的膨胀量曲线；
（2）以1~3℃温度区间划分为若干等距区间，将每个区间中的温度的平均值作为横坐标，将每个区间中的膨胀量的平均值作为纵坐标，对降温数据进行平均处理；
（3）对膨胀量-温度数据进行进行平滑处理，如果平滑效果不好，进行第二次平滑处理；
（4）将步骤（3）平滑处理后的高温段数据按照温度由高到低依次命名为（T1，D1）、（T2，D2）、（T3，D3）……（Tn，Dn），对（T1，D1）、（T2，D2）……（Tn，Dn）对应的n组数据进行线性回归得到拟合直线，其中n≥2，T表示温度，D表示膨胀量；若|D曲线-D直线|≥D总/T温差，则该点即为拟合直线与膨胀量-温度曲线的分离点，其中，D总为最大膨胀量与最小膨胀量之差，T温差为最大膨胀量对应的温度与最小膨胀量对应的温度之差；若|Tm-T1|=80±1℃，则Tm即为相变开始点，其中Tm为拟合直线与膨胀量-温度曲线的分离点对应的温度；此时（T1，D1）、（T2，D2）……（Tm，Dm）的拟合直线即为高温段切线；D=k×T+d，式中，k为高温段切线斜率，d为截距，D为膨胀量，T为温度；对于低温相变结束点和低温段切线，采用同样的方法处理；
（5）利用杠杆定律可计算相变开始点和相变结束点对应温度区间内的马氏体组织转变量，f=Xa/（Xa+Xb），式中：f为马氏体组织转变量，Xa为高温段数据拟合直线与平滑后膨胀曲线的竖直距离，Xb为低温段数据拟合直线与平滑后膨胀曲线的竖直距离。
2.根据权利要求1所述的计算马氏体组织转变量及相变点的方法，其特征是：待测材料为发生相变的材料。
3.根据权利要求1所述的计算马氏体组织转变量及相变点的方法，其特征是：步骤（3）中以30个膨胀量-温度数据为一个平滑组进行平滑，第二次平滑以20-40个膨胀量-温度数据为一个平滑组进行平滑。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811156239.2A CN109470734B (zh) | 2018-09-30 | 2018-09-30 | 一种计算马氏体组织转变量及相变点的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811156239.2A CN109470734B (zh) | 2018-09-30 | 2018-09-30 | 一种计算马氏体组织转变量及相变点的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109470734A true CN109470734A (zh) | 2019-03-15 |
| CN109470734B CN109470734B (zh) | 2021-09-14 |
Family
ID=65663406
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811156239.2A Active CN109470734B (zh) | 2018-09-30 | 2018-09-30 | 一种计算马氏体组织转变量及相变点的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109470734B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113063813A (zh) * | 2021-03-15 | 2021-07-02 | 鞍钢股份有限公司 | 一种建立钢铁材料的连续冷却相变模型的方法 |
| CN114646284A (zh) * | 2022-03-16 | 2022-06-21 | 重庆邮电大学 | 一种精密齿轮几何参数的温度补偿方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6654713B1 (en) * | 1999-11-22 | 2003-11-25 | Hewlett-Packard Development Company, L.P. | Method to compress a piecewise linear waveform so compression error occurs on only one side of the waveform |
| CN105136842A (zh) * | 2015-07-29 | 2015-12-09 | 南京钢铁股份有限公司 | 一种适用于钢临界点测试的膨胀曲线分析方法及该方法的用途 |
-
2018
- 2018-09-30 CN CN201811156239.2A patent/CN109470734B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6654713B1 (en) * | 1999-11-22 | 2003-11-25 | Hewlett-Packard Development Company, L.P. | Method to compress a piecewise linear waveform so compression error occurs on only one side of the waveform |
| CN105136842A (zh) * | 2015-07-29 | 2015-12-09 | 南京钢铁股份有限公司 | 一种适用于钢临界点测试的膨胀曲线分析方法及该方法的用途 |
Non-Patent Citations (2)
| Title |
|---|
| Z.L.XIE,ETC: "stabilization of retained austenite due to partial martensitic transformations", <ACTA METALLURGICAET MATERIALLA> * |
| 操龙飞: "金属材料的热膨胀特征研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113063813A (zh) * | 2021-03-15 | 2021-07-02 | 鞍钢股份有限公司 | 一种建立钢铁材料的连续冷却相变模型的方法 |
| CN113063813B (zh) * | 2021-03-15 | 2023-11-10 | 鞍钢股份有限公司 | 一种建立钢铁材料的连续冷却相变模型的方法 |
| CN114646284A (zh) * | 2022-03-16 | 2022-06-21 | 重庆邮电大学 | 一种精密齿轮几何参数的温度补偿方法 |
| CN114646284B (zh) * | 2022-03-16 | 2023-08-11 | 重庆邮电大学 | 一种精密齿轮几何参数的温度补偿方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109470734B (zh) | 2021-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109470734A (zh) | 一种计算马氏体组织转变量及相变点的方法 | |
| CN109696788A (zh) | 一种基于显示面板的快速自动对焦方法 | |
| CN106296031A (zh) | 一种基于动态topsis的空中目标威胁评估方法 | |
| CN110443006B (zh) | 一种通过测量面积计算金属材料两相转变体积比例的方法 | |
| CN106304331A (zh) | 一种WiFi指纹室内定位方法 | |
| CN102331391B (zh) | 一种测算耐热钢中多相颗粒聚集粗化程度的方法 | |
| CN108615054B (zh) | 排水管网节点间相似性衡量的综合指标构建方法 | |
| CN105021533A (zh) | 一种快速检测x70屈服强度、抗拉强度和表面硬度的方法 | |
| Serigos et al. | Evaluation of 3D automated systems for the measurement of pavement surface cracking | |
| CN107462349A (zh) | 一种超声波热量表全温度全流量范围连续补偿的流量补偿方法 | |
| CN107702667B (zh) | 一种新型x射线测厚仪材质补偿系统的建立方法 | |
| CN110307910B (zh) | 一种步进梁加热炉钢坯温度获取方法及装置 | |
| CN107900114A (zh) | 一种对冷轧带钢板形质量进行评价的方法及装置 | |
| CN101201398B (zh) | 利用邻站测距辅助测角的非直达波定位方法 | |
| CN108170951B (zh) | 基于采样数据时空匹配示踪试验的纵向离散系数确定方法 | |
| Čelar et al. | Methodology for control delay estimation using new algorithm for critical points identification | |
| CN103196793B (zh) | 建立表征炼焦煤流变性能模型的方法 | |
| Nurulhuda et al. | Study of agricultural economic potential in West Kalimantan using Regional Analysis Techniques | |
| CN111157567A (zh) | 一种测量低合金高强钢连续冷却两相转变体积比例的方法 | |
| CN106680469A (zh) | 一种锂电池隔膜专用聚丙烯树脂的评价方法 | |
| CN107818240A (zh) | 基于伴随同化的纵向离散系数反演方法 | |
| CN110455662A (zh) | 铁素体钢冲击功与断裂韧性经验关系式确定方法 | |
| CN105572151B (zh) | 同一射线检测不同厚度材料透照参数确定方法 | |
| CN106932467B (zh) | 一种蒸汽发生器传热管管板缝隙区涡流精准测量方法 | |
| Kazakov et al. | Quantitative Description of Microstructural Banding in Steels |
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 |