CN111968220A - 一种基于响应曲面法的真空烧结炉结构参数优化方法 - Google Patents
一种基于响应曲面法的真空烧结炉结构参数优化方法 Download PDFInfo
- Publication number
- CN111968220A CN111968220A CN202010753286.6A CN202010753286A CN111968220A CN 111968220 A CN111968220 A CN 111968220A CN 202010753286 A CN202010753286 A CN 202010753286A CN 111968220 A CN111968220 A CN 111968220A
- Authority
- CN
- China
- Prior art keywords
- vacuum sintering
- sintering furnace
- workpiece
- furnace
- dimensional 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22F—WORKING METALLIC POWDER; MANUFACTURE OF ARTICLES FROM METALLIC POWDER; MAKING METALLIC POWDER; APPARATUS OR DEVICES SPECIALLY ADAPTED FOR METALLIC POWDER
- B22F3/00—Manufacture of workpieces or articles from metallic powder characterised by the manner of compacting or sintering; Apparatus specially adapted therefor ; Presses and furnaces
- B22F3/10—Sintering only
- B22F3/1003—Use of special medium during sintering, e.g. sintering aid
- B22F3/1007—Atmosphere
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/25—Process efficiency
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Mechanical Engineering (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Powder Metallurgy (AREA)
Abstract
本发明提供一种基于响应曲面法的真空烧结炉结构参数优化方法,包括以下步骤:S1:建立真空烧结炉三维模型,并在真空烧结炉三维模型中进行数值仿真采集数据;S2:建立真空烧结炉的评价指标;S3:基于真空烧结炉三维模型,采用响应曲面法设计结构参数优化试验,计算评价指标;S4:根据计算得到的评价指标,建立结构参数优化的回归模型;S5:对回归模型进行方差分析,得到方差分析结果;S6:对回归模型进行回归分析,得到回归分析结果;S7:根据方差分析和回归分析的结果确定优化的真空烧结炉结构参数。本发明提供一种基于响应曲面法的真空烧结炉结构参数优化方法,解决了目前我国的真空烧结炉热效率水平不够高的问题。
Description
技术领域
本发明涉及结构参数优化技术领域,更具体的,涉及一种基于响应曲面法的真空烧结炉结构参数优化方法。
背景技术
真空烧结炉是粉末冶金行业中不可缺少的设备之一,属于工业炉中的一种高能耗设备。在真空烧结炉的工作过程中,金属粉末颗粒在高温下扩散形成烧结颈,从而连接金属颗粒使粒子间的空隙减少,让零件获得较高的力学强度。为了保证这一现象的正常发生,需要在真空烧结炉工作时中充入具有还原性的气氛,将颗粒表面的的氧化层进行还原,同时保证其他烧结材料的表面不被氧化。我国有着大量的真空烧结炉制造企业,满足了粉末冶金产业的烧结炉需求,同时因为激烈的价格竞争,使得国内的真空烧结炉价格和进口设备对比有着明显的优势,但是气氛和温度的控制技术,尤其是能耗方面存在很大差距,国际上的真空烧结炉的热效率水平达到了60%,而我国的真空烧结炉热效率水平不到35%。
目前我国对于真空烧结炉的研究大部分集中于炉内温度场的研究或者烧结炉的温控系统研究,而缺少对节能降耗的研究,导致我国的真空烧结炉热效率水平不够高。
现有技术中,还没有结合响应曲面法对烧结炉结构参数进行优化的方案,如2015年5月13日授权的中国专利,一种基于响应曲面分析的声屏障优化设计方法,公告号为CN102663199B,采用响应曲面分析优化声屏障设计参数,有效控制声屏障的降噪效果和工程造价,但不能直接用于优化真空烧结炉结构参数。
发明内容
本发明为克服目前我国的真空烧结炉热效率水平不够高的技术缺陷,提供一种基于响应曲面法的真空烧结炉结构参数优化方法。
为解决上述技术问题,本发明的技术方案如下:
一种基于响应曲面法的真空烧结炉结构参数优化方法,包括以下步骤:
S1:建立真空烧结炉三维模型,并在真空烧结炉三维模型中进行数值仿真采集数据;
S2:建立真空烧结炉的评价指标,包括炉最大能效指标和工件温度均匀性指标;
S3:基于真空烧结炉三维模型,采用响应曲面法设计结构参数优化试验,计算评价指标;
S4:根据计算得到的评价指标,建立结构参数优化的回归模型;
S5:对回归模型进行方差分析,得到方差分析结果;
S6:对回归模型进行回归分析,得到回归分析结果;
S7:根据方差分析和回归分析的结果确定最终优化结果,最终优化结果即为优化的真空烧结炉结构参数。
优选的,真空烧结炉三维模型包括真空烧结炉本体、工件和加热体;采集的数据包括工件从常温至最终保温阶段过程中所吸收的热量Qg,真空烧结炉三维模型中电源输入的功率P,真空烧结炉三维模型中电源输入的电能Qq,工件的比热容cg,工件的密度ρg,加热时间t5,工件整体的平均温度T,工件的初始温度T0,工件温度的标准差σ,工件上所取温度点的数量N,工件上第i个点的温度Ti。
优选的,在步骤S3中,设计L(4≤L≤8)次不同加热体结构下的两因素多水平的试验;其中,每次试验均有M(20≤M≤30)个试验组,两因素分别为真空烧结炉三维模型内部前后加热体距炉内壁距离参数d1和左右加热体距炉内部距离参数d2,响应值为炉最大能效指标和工件温度均匀性指标,且各因素在不同加热体结构下的取值一致。
优选的,在步骤S3中,具体包括以下步骤:
S3.1:在真空烧结炉三维模型中,修改加热体的结构组成和其在真空烧结炉本体内的位置参数,并去除工件,形成真空烧结炉的空炉三维模型;
S3.2:对空炉三维模型进行有效加热空间的测定,得到有效加工空间Vgmax;
S3.3:在空炉三维模型的中心添置工件,形成含工件的真空烧结炉三维模型;
S3.4:求解含工件的真空烧结炉三维模型的炉最大能效指标和工件温度均匀性指标;
S3.5:重复步骤S3.1至步骤S3.4,直到完成L次不同加热体结构下的两因素多水平的试验。
优选的,在步骤S3.3中,工件为六面体的氮化硅工件。
优选的,采用炉最大效能比值η来评价炉最大能效指标:
采用工件温度变异性系数ε来评价工件温度均匀性指标:
优选的,结构参数优化的回归模型包括:加热体在真空烧结炉本体内的位置参数与炉最大能效指标的多元二次回归模型和加热体在真空烧结炉本体内的位置参数与工件温度均匀性指标的多元二次回归模型。
优选的,在步骤S5中,方差分析包括P值分析和复相关系数R2分析,分别得到显著性分析结果和可信度分析结果。
优选的,在步骤S6中,回归分析包括残差正态分布图分析和预测值与实验值对比图分析,得到准确性分析结果。
优选的,根据方差分析和回归分析的结果选取炉最大能效指标最高的结构参数作为最终优化结果,得到优化的真空烧结炉结构参数。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供了一种基于响应曲面法的真空烧结炉结构参数优化方法,根据真空烧结炉三维模型和评价指标设计基于响应曲面法的结构参数优化试验,得到结构参数优化的回归模型,并通过对回归模型进行方差分析和回归分析,从而确定最终优化结果,得到优化的真空烧结炉结构参数,提高了真空烧结炉的热效率,实现了节能降耗。
附图说明
图1为本发明的技术方案实施步骤流程图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,一种基于响应曲面法的真空烧结炉结构参数优化方法,包括以下步骤:
S1:建立真空烧结炉三维模型,并在真空烧结炉三维模型中进行数值仿真采集数据;
S2:建立真空烧结炉的评价指标,包括炉最大能效指标和工件温度均匀性指标;
S3:基于真空烧结炉三维模型,采用响应曲面法设计结构参数优化试验,计算评价指标;
S4:根据计算得到的评价指标,建立结构参数优化的回归模型;
S5:对回归模型进行方差分析,得到方差分析结果;
S6:对回归模型进行回归分析,得到回归分析结果;
S7:根据方差分析和回归分析的结果确定最终优化结果,最终优化结果即为优化的真空烧结炉结构参数。
在实施过程中,采用ANSYS软件对真空烧结炉进行仿真建模,建立真空烧结炉三维模型和评价指标,并设计基于响应曲面法的结构参数优化试验,得到结构参数优化的回归模型,然后通过对回归模型进行方差分析和回归分析,从而确定最终优化结果,得到优化的真空烧结炉结构参数,提高了真空烧结炉的热效率,实现了节能降耗。
更具体的,真空烧结炉三维模型包括真空烧结炉本体、工件和加热体;采集的数据包括工件从常温至最终保温阶段过程中所吸收的热量Qg,真空烧结炉三维模型中电源输入的功率P,真空烧结炉三维模型中电源输入的电能Qq,工件的比热容cg,工件的密度ρg,加热时间t5,工件整体的平均温度T,工件的初始温度T0,工件温度的标准差σ,工件上所取温度点的数量N,工件上第i个点的温度Ti。
在实施过程中,结合传热学知识建立真空烧结炉三维模型并采集数据,且对于真空烧结炉三维模型的传热方式只考虑内部热传导和辐射换热,不考虑热对流的影响。
更具体的,在步骤S3中,设计四次不同加热体结构下的两因素多水平的试验;其中,每次试验均有25个试验组,两因素分别为真空烧结炉三维模型内部前后加热体距炉内壁距离参数d1和左右加热体距炉内部距离参数d2,响应值为炉最大能效指标和工件温度均匀性指标,且各因素在不同加热体结构下的取值一致。
在实施过程中,采用Design-Expert软件基于响应曲面法设计四次不同加热体结构组成下加热体位置参数与炉最大能效指标和工件温度均匀性指标的多目标优化模型进行试验,每次均有25个试验组。
更具体的,在步骤S3中,具体包括以下步骤:
S3.1:在真空烧结炉三维模型中,修改加热体的结构组成和其在真空烧结炉本体内的位置参数,并去除工件,形成真空烧结炉的空炉三维模型;
S3.2:将空炉三维模型导入ANSYS软件中的Fluent仿真模块中对空炉三维模型进行有效加热空间的测定,得到有效加工空间Vgmax;
S3.3:在空炉三维模型的中心添置工件,形成含工件的真空烧结炉三维模型;
S3.4:将含工件的真空烧结炉三维模型导入Fluent仿真模块,求解炉最大能效指标和工件温度均匀性指标;
S3.5:重复步骤S3.1至步骤S3.4,直到完成四次不同加热体结构下的两因素多水平的试验。
在实施过程中,步骤S3.1至步骤S3.4为一次两因素多水平试验,通过修改加热体的结构组成和其在真空烧结炉本体内的位置参数,采用响应曲面法建立真空烧结炉结构参数与两项评价指标之间的函数模型,并对真空烧结炉最大能效指标和工件温度均匀性指标进行多目标优化求解,实现四次不同加热体结构下的两因素多水平的试验设计。
更具体的,在步骤S3.3中,工件为六面体的氮化硅工件。
更具体的,采用炉最大效能比值η来评价炉最大能效指标:
采用工件温度变异性系数ε来评价工件温度均匀性指标:
在实施过程中,结合热平衡原理建立炉最大能效比指标和工件温度均匀性的质量指标。在初始结构参数不变的情况下,原始的氮化硅真空烧结炉的最大能效比值η=11.165%。
更具体的,结构参数优化的回归模型包括:加热体在真空烧结炉本体内的位置参数与炉最大能效指标的多元二次回归模型和加热体在真空烧结炉本体内的位置参数与工件温度均匀性指标的多元二次回归模型。
更具体的,在步骤S5中,方差分析包括P值分析和复相关系数R2分析,分别得到显著性分析结果和可信度分析结果。
在实施过程中,若P≥0.05,则因素对响应值的影响不显著;若0.01≤P<0.05,则因素对响应值的影响显著;若0.0001≤P<0.01,则因素对响应值的影响高度显著;若P<0.0001,则因素对响应值的影响极高度显著。若R2的值越接近1,则回归模型的可信度越高。
更具体的,在步骤S6中,回归分析包括残差正态分布图分析和预测值与实验值对比图分析,得到准确性分析结果。
在实施过程中,若残差正态分布图中各点的位置分布越靠近图对角线,则该试验组别的残差值呈正态分布且响应面拟合的回归模型的准确性越高;若预测值与实验值的对比图中各点的位置分布越靠近图对角线,则该回归模型的准确性越高。
更具体的,根据方差分析和回归分析的结果选取炉最大能效指标最高的结构参数作为最终优化结果,得到优化的真空烧结炉结构参数。
在实施过程中,优化的真空烧结炉结构参数为:d1=26mm、d2=75mm、Vgmax=275mm*275mm*406mm,炉最大能效比值为η=12.3105%,与原始的真空烧结炉的η=11.165%相比,炉最大能效比值提高了10.26%,提高了真空烧结炉的热效率,实现了节能降耗。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,包括以下步骤:
S1:建立真空烧结炉三维模型,并在真空烧结炉三维模型中进行数值仿真采集数据;
S2:建立真空烧结炉的评价指标,包括炉最大能效指标和工件温度均匀性指标;
S3:基于真空烧结炉三维模型,采用响应曲面法设计结构参数优化试验,计算评价指标;
S4:根据计算得到的评价指标,建立结构参数优化的回归模型;
S5:对回归模型进行方差分析,得到方差分析结果;
S6:对回归模型进行回归分析,得到回归分析结果;
S7:根据方差分析和回归分析的结果确定最终优化结果,最终优化结果即为优化的真空烧结炉结构参数。
3.根据权利要求2所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,在步骤S3中,设计L(4≤L≤8)次不同加热体结构下的两因素多水平的试验;其中,每次试验均有M(20≤M≤30)个试验组,两因素分别为真空烧结炉三维模型内部前后加热体距炉内壁距离参数d1和左右加热体距炉内部距离参数d2,响应值为炉最大能效指标和工件温度均匀性指标,且各因素在不同加热体结构下的取值一致。
4.根据权利要求3所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,在步骤S3中,具体包括以下步骤:
S3.1:在真空烧结炉三维模型中,修改加热体的结构组成和其在真空烧结炉本体内的位置参数,并去除工件,形成真空烧结炉的空炉三维模型;
S3.2:对空炉三维模型进行有效加热空间的测定,得到有效加工空间Vgmax;
S3.3:在空炉三维模型的中心添置工件,形成含工件的真空烧结炉三维模型;
S3.4:求解含工件的真空烧结炉三维模型的炉最大能效指标和工件温度均匀性指标;
S3.5:重复步骤S3.1至步骤S3.4,直到完成L次不同加热体结构下的两因素多水平的试验。
5.根据权利要求4所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,在步骤S3.3中,工件为六面体的氮化硅工件。
7.根据权利要求2所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,结构参数优化的回归模型包括:加热体在真空烧结炉本体内的位置参数与炉最大能效指标的多元二次回归模型和加热体在真空烧结炉本体内的位置参数与工件温度均匀性指标的多元二次回归模型。
8.根据权利要求7所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,在步骤S5中,方差分析包括P值分析和复相关系数R2分析,分别得到显著性分析结果和可信度分析结果。
9.根据权利要求7所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,在步骤S6中,回归分析包括残差正态分布图分析和预测值与实验值对比图分析,得到准确性分析结果。
10.根据权利要求1所述的一种基于响应曲面法的真空烧结炉结构参数优化方法,其特征在于,根据方差分析和回归分析的结果选取炉最大能效指标最高的结构参数作为最终优化结果,得到优化的真空烧结炉结构参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010753286.6A CN111968220B (zh) | 2020-07-30 | 2020-07-30 | 一种基于响应曲面法的真空烧结炉结构参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010753286.6A CN111968220B (zh) | 2020-07-30 | 2020-07-30 | 一种基于响应曲面法的真空烧结炉结构参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111968220A true CN111968220A (zh) | 2020-11-20 |
CN111968220B CN111968220B (zh) | 2023-06-23 |
Family
ID=73363670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010753286.6A Active CN111968220B (zh) | 2020-07-30 | 2020-07-30 | 一种基于响应曲面法的真空烧结炉结构参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111968220B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112893778A (zh) * | 2021-02-18 | 2021-06-04 | 福州大学 | 基于计算机辅助的消失模内花键叉头铸造工艺设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040216558A1 (en) * | 2003-04-25 | 2004-11-04 | Robert Mariani | Method of forming sintered valve metal material |
CN107220400A (zh) * | 2017-04-05 | 2017-09-29 | 中国石油化工股份有限公司 | 一种基于响应面法的高温空气燃烧组织优化方法 |
CN108804789A (zh) * | 2018-05-29 | 2018-11-13 | 广东工业大学 | 一种陶瓷喷雾干燥过程能耗建模方法 |
CN110580377A (zh) * | 2019-08-06 | 2019-12-17 | 天津大学 | 一种基于响应曲面法的植草沟设计参数优化方法 |
-
2020
- 2020-07-30 CN CN202010753286.6A patent/CN111968220B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040216558A1 (en) * | 2003-04-25 | 2004-11-04 | Robert Mariani | Method of forming sintered valve metal material |
CN107220400A (zh) * | 2017-04-05 | 2017-09-29 | 中国石油化工股份有限公司 | 一种基于响应面法的高温空气燃烧组织优化方法 |
CN108804789A (zh) * | 2018-05-29 | 2018-11-13 | 广东工业大学 | 一种陶瓷喷雾干燥过程能耗建模方法 |
CN110580377A (zh) * | 2019-08-06 | 2019-12-17 | 天津大学 | 一种基于响应曲面法的植草沟设计参数优化方法 |
Non-Patent Citations (1)
Title |
---|
张永学 等: "基于响应面法的油田蓄热式加热炉燃烧器优化设计", 中国石油大学学报(自然科学版), vol. 41, no. 2, pages 156 - 161 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112893778A (zh) * | 2021-02-18 | 2021-06-04 | 福州大学 | 基于计算机辅助的消失模内花键叉头铸造工艺设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111968220B (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111400958B (zh) | 一种基于ansys的高温碳化炉氧含量分布的模拟方法 | |
US8437991B2 (en) | Systems and methods for predicting heat transfer coefficients during quenching | |
CN111651908B (zh) | 一种基于ansys的高温碳化炉多场耦合应力分布模拟方法 | |
CN111222685B (zh) | 基于模型迁移的锅炉宽负荷NOx排放浓度预测方法 | |
CN111627503B (zh) | 一种氧化铝陶瓷基复合材料涂层激光熔覆制造过程中应力场的预测方法 | |
CN106897520B (zh) | 一种含有模糊参数的传热系统可靠性分析方法 | |
CN110579201B (zh) | 一种基于差分进化算法的平面度评定方法 | |
CN111709096A (zh) | 一种强化自然对流换热的异型翅片结构设计方法 | |
CN111968220A (zh) | 一种基于响应曲面法的真空烧结炉结构参数优化方法 | |
Zhang et al. | Automatic optimization design of a feeder extrusion die with response surface methodology and mesh deformation technique | |
CN112528569A (zh) | 基于数字孪生的工业加热炉温度场预测方法 | |
CN111881630A (zh) | 一种核设施近场环境模拟数值风洞的建立系统及方法 | |
CN110362861B (zh) | 一种考虑能效的模具结构参数多目标优化方法 | |
CN116738891A (zh) | 一种增强飞行器流场模拟稳定性的lu-sgs改进方法 | |
CN110029198B (zh) | 一种高炉冷却系统的冷却效果的计算机标定方法 | |
CN109918765B (zh) | 一种高速飞行器舱内热环境耦合精细计算方法 | |
Faulkner et al. | Fast prediction of indoor airflow distribution inspired by synthetic image generation artificial intelligence | |
CN112632856A (zh) | 回焊炉的传送带速度及温度控制方法 | |
Tajik et al. | Multi-objective optimization of aluminum anode baking process employing a response surface methodology | |
CN116796546A (zh) | 基于小样本学习的高温合金增材制造工艺优化预测方法 | |
CN106909747B (zh) | 一种热对流扩散系统中模糊参数隶属度函数识别方法 | |
CN111368402B (zh) | 一种管材挤压工艺优化方法 | |
CN103699748B (zh) | 撬装化工厂框架吊装点位置的确定方法 | |
WO2022011724A1 (zh) | 一种基于workbench的低温碳化炉多耦合场应力分布的模拟方法 | |
CN111859671A (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 |