CN110814053A - 基于尝试法的棒材水箱反向自决策控温智能计算方法 - Google Patents

基于尝试法的棒材水箱反向自决策控温智能计算方法 Download PDF

Info

Publication number
CN110814053A
CN110814053A CN201911125953.XA CN201911125953A CN110814053A CN 110814053 A CN110814053 A CN 110814053A CN 201911125953 A CN201911125953 A CN 201911125953A CN 110814053 A CN110814053 A CN 110814053A
Authority
CN
China
Prior art keywords
temperature
water
water tank
cooling
rolled piece
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
CN201911125953.XA
Other languages
English (en)
Other versions
CN110814053B (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.)
Huatian Engineering and Technology Corp MCC
Anhui University of Technology AHUT
Anhui Changjiang Steel Co Ltd
Original Assignee
Huatian Engineering and Technology Corp MCC
Anhui University of Technology AHUT
Anhui Changjiang Steel Co Ltd
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 Huatian Engineering and Technology Corp MCC, Anhui University of Technology AHUT, Anhui Changjiang Steel Co Ltd filed Critical Huatian Engineering and Technology Corp MCC
Priority to CN201911125953.XA priority Critical patent/CN110814053B/zh
Publication of CN110814053A publication Critical patent/CN110814053A/zh
Application granted granted Critical
Publication of CN110814053B publication Critical patent/CN110814053B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B37/00Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • B21B37/74Temperature control, e.g. by cooling or heating the rolls or the product
    • B21B37/76Cooling control on the run-out table
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B45/00Devices for surface or other treatment of work, specially combined with or arranged in, or specially adapted for use in connection with, metal-rolling mills
    • B21B45/02Devices for surface or other treatment of work, specially combined with or arranged in, or specially adapted for use in connection with, metal-rolling mills for lubricating, cooling, or cleaning
    • B21B45/0203Cooling
    • B21B45/0209Cooling devices, e.g. using gaseous coolants
    • B21B45/0215Cooling devices, e.g. using gaseous coolants using liquid coolants, e.g. for sections, for tubes
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B2201/00Special rolling modes
    • B21B2201/04Ferritic rolling
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B2201/00Special rolling modes
    • B21B2201/06Thermomechanical rolling

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Metal Rolling (AREA)
  • Control Of Heat Treatment Processes (AREA)

Abstract

本发明基于尝试法的棒材水箱反向自决策控温智能计算方法,包括下述步骤:1)计算轧件温度分布;2)计算得出不同水箱的组合下轧件上冷床时温度内外分布,及供水系统载荷;3)以表层温度值和温度内外差值为目标函数,查找出不越界的水箱开关组合;4)尝试按预定刻度值,以供水系统最小载荷为目标函数,逐一调节每个水箱水量阀大小;5)计算调低水量阀刻度后的轧件整体的时间‑温度演变曲线,与预期的cct冷却路径进行比对,满足要求则结束退出。本发明实现了快速自决策大水量快速冷却对各个水箱的控制参数,精细化按照水冷调控组织性能要求,用最快的冷却速度,在马氏体形成温度以上,确保轧件得到铁素体和珠光体,达标力学性能。

Description

基于尝试法的棒材水箱反向自决策控温智能计算方法
技术领域
本发明涉及领域为热轧棒材生产的控制冷却智能化算法,特别涉及一种基于尝试法的棒材水箱反向自决策控温智能计算方法。
背景技术
热轧棒材通常包含普通棒材(螺纹钢筋)和优棒特钢(如优质碳结钢、合金结构钢、冷镦钢等),二者合计产能长期占据我国钢铁行业约30%份额,在其制造全流程中,炼钢、连铸、轧钢工序段的管控,分别从成分、偏析、组织等方面对产品质量形成决定性影响。热轧棒材主要用作结构件,更关注屈服强度、拉伸强度、延伸率等力学性能,由于焊接条件的限制,以及地震/火灾等突发载荷服役条件的要求,通常热轧棒材追求以细化的10~11级铁素体+珠光体组织为目标,长期以来,主要基于控轧控冷工艺及配套喷嘴、水箱、水泵、控制予以实现。
一方面,轧材由于表面相较内部,头尾部相较中间,以及角部换热条件好,必然造成轧材温度不均匀;一方面由于棒材热轧过程中,前工序中的连铸、热送、加热、粗轧、中轧时存在工艺质量波动,引起中间轧件温度存在偏差;一方面,控冷设备,尤其喷嘴、导槽等属于消耗件,在对轧件快速冷却及支撑轧件运行中,形成的磨损将引起冷却不均,这三方面共同导致对传统冷却的控制,要实现高精度和稳定度较难实现,对模型和系统的可靠性要求高。此外,生产企业产线难以达到完全专业化,不可能只生产一种钢种及规格产品,必须随生产计划和任务安排,进行快节奏紧衔接的排产和换制度工艺,生产产品的变化,核心是钢种与规格对应的工艺制度参数变化,为保证高产能低成本,需要能迅速根据产线上新的变化形势情况,制定工艺参数,控制冷却用模型和算法成为关键。
目前,现场对热轧棒材控冷参数的制定仍以经验法为主,通常在积累出大量钢种和规格的控冷工艺值后,在经现场实际跑钢约3~5次调试后,得出适用值,固定使用一段时间,当最终检测值如上冷床温度、组织性能抽检不合格后,再更换参数,这种方法不仅费时费力,而且没有理论指导,难以固化,无法最大化发挥产线产能,此外也有基于傅里叶传热原理,通过差分法,形成轧件温度计算软件,结合经验或实验获得的水压、水流量等对换热系数的公式,预测温度变化,这种方法相较纯经验法有所进步,但随工况的变化,换热系数值并不稳定,需要更新,常在软件中未予体现。此外,实际应用中,最需要的是能在给定的新钢种和规格下,对于特定的产线水箱布置情况,基于上冷床温度值 (或要求值),快速方向给出具体的每一个水箱的开关与否,及开口下的开度阀值量,而不只是‘均匀’调控每个水箱,达到冷却温度值即可,按照材料冶金学调控螺纹钢棒材组织原理和实际控冷要求,应以大水量快速冷却,在满足冷却最终温度和温度不形成马氏体情况下,以开最少的水箱,用最快的冷却速度才最合理,这些要求已超出目前棒材控冷模型和系统的能力。
发明内容
为克服上述缺陷,本发明的目的是提出了一种更为先进合理、操作简易、投资较省的基于尝试法的棒材水箱反向自决策控温智能计算方法,以完成各个水箱工艺参数基于在线钢种和规格要求进行的调整。
为达到上述目的,本发明基于尝试法的棒材水箱反向自决策控温智能计算方法,包括下述步骤:
1)根据输入的钢种、规格,以及穿水前工艺参数和布置,计算出终轧机后轧件温度分布;
2)计算得出不同水箱的水量阀在满刻度或零刻度的组合下轧件上冷床时温度内外分布,及供水系统(含总管、支管)载荷;
3)以表层温度值和温度内外差值为目标函数,将计算得出的温度分布进行对比,查找出不越界的水箱开关组合;
4)尝试按预定刻度值,以供水系统最小载荷为目标函数,置入随机±5%~10%的开口度刻度值,逐一调节每个水箱水量阀大小;
5)计算每个水箱调低水量阀刻度后的轧件整体的时间-温度演变曲线,与预期的满足组织性能的cct冷却路径,进行比对,满足要求则结束退出,否则返回步骤4。
进一步的,所述的方法包括下述步骤:所述的轧件温度分布是基于傅里叶传热理论的有限差分算法计算得出的。
进一步的,所述的按预定刻度值逐一调整每个水箱水量阀的大小的步骤具体为:以表层温度值、温度内外差值、供水系统最小载荷为目标函数,所述的按预定刻度值逐一调整每个水箱水量阀的大小的步骤具体为:第1~第j-1个水箱置入随机±5%~10%的开口度刻度值,按预定刻度值调低第j个水箱调节水量阀的大小。
进一步的,还包括调低预定刻度值的步骤。
本发明基于尝试法的棒材水箱反向自决策控温智能计算方法,以傅里叶传热二维差分演化方程为基础,对于预设的(或者检测的)棒材上冷床表面温度和截面温差,先尝试轧件历经一个或多个水箱全开不同条件下,进行冷却后的上冷床温度情况,然后依据最大冷却和全空冷结果判定给定的预设要求是否越界,在不越界条件下,找出开水箱的总个数和具体对象,接着在确定开水箱个数后,按照控冷工艺要求,对可调的具体水箱,就开口度参数在范围内进行逐步寻优,找到最佳值,对输出的轧件冷却曲线与需求的满足产品组织性能的冷却路径进行比对,达到要求时则输出结果。实现了快速自决策大水量快速冷却对各个水箱的控制参数,精细化按照水冷调控组织性能要求,以开最少的水箱,用最快的冷却速度,在马氏体形成温度以上,确保轧件得到铁素体和珠光体,达标力学性能。
附图说明
图1一种基于尝试法的棒材水箱反向自决策控温智能计算方法程序框图
具体实施方式
典型热轧棒材轧制的生产线中,通常布置有18个机架,粗轧机组6架、中轧机组6架、精轧机组6架。开轧温度一般在950℃~1200℃,经18道次轧制后,出终轧机后的轧件温度约在900℃~1100℃,一方面,要将高温轧件温度快速降低到室温便于后续精整;
一方面,根据低碳钢材料冶金连续冷却转变原理,合理的冷却路径引起的冷却速度,可保证高温轧件高温奥氏体形成细小的铁素体,促进珠光体细化,得到良好的组织性能;
一方面,高冷却速度的冷却强度并不能太过,以至形成不利的马氏体组织。因此,在精轧机组结束后,会布置多组水箱,通过对每个水箱的水流量、水压控制,实现对穿行的轧件冷却强度控制,水箱数量一般在1~4组之间,同时两个相邻水箱间,还会放置一定长度距离的自然空冷距离,来形成心部到表面的恢复。
传统的棒材轧件温度调控中,通常只以最终的棒材表面温度为目标,忽略了本质上棒材所需的冷却路径控制对微观组织的需求,所开发的模型和系统,仅只以傅里叶热传导差分方程实现特定工艺参数下的温度计算,没有反算水箱最优参数的模型算法。
本发明提出了一种基于尝试法的棒材水箱反向自决策控温智能计算方法,首先根据输入的钢种、规格,以及穿水前工艺参数和布置,基于傅里叶传热理论的有限差分算法,得出出终轧机后轧件温度分布,以此作为后续计算的初始条件。随后根据穿水水箱现场工艺布置,计算得出不同水箱组合(仅考虑开关) 下轧件上冷床时的温度值和分布,与冷床上钢温度要求(检测值),及温度分布要求(一般<50℃)比对,如果越界,直接退出程序计算。
在不越界下,尝试获得最优水箱开关组合,以表层温度值、温度内外差值、供水系统最小载荷为目标函数,按预定刻度值逐一调整每个水箱水量阀的大小的步骤具体为:第1~第j-1个水箱置入随机±5%~10%的开口度刻度值,第j个水箱调节水量阀的顺序,间隔刻度先大后小的路径寻优冷却路径,得出后,对计算轧件整体的时间-温度演变曲线,并与预期的满足组织性能的cct冷却路径,进行比对,满足要求则结束退出程序,否则取一组新的水箱组合,再次进行尝试,直到满足误差内。
基于分间隔刻度的寻优,可快速缩小最优参数的所在范围,进而缩短查找时间,更能快速获得最佳结果。通过本发明的处理,可达到在最低投资下,基于稳定控冷工艺的冷却路径要求,提高轧件产品组织性能。
具体操作时,调节水量阀预定刻度值较大时,寻优速度较快,例如,25刻度时,但精度相应会比较低。而因此,调节水量阀预定刻度值较小时,寻优速度相对会慢一些,例如,5刻度时。但精度相应会比较高。因此根据实际需要调整预定刻度值,会有利于在精度和速度上得到保证。预定刻度值的范围可以为 1-30刻度。
综上可知,本发明提出的新型一种基于尝试法的棒材水箱反向自决策控温智能计算方法,利用傅里叶差分算法计算温度为基础,基于终轧前钢种、规格、工艺参数和布置,轧后水箱布置,针对满足产品组织性能的cct曲线冷却曲线目标,在计算中加入水箱组合开关和尝试从末个水箱反算开口度,形成计算冷却曲线,通过比对,实现反向子决策的智能计算,可在尽可能少的改建现有轧线条件下,达到尽可能满足轧件冷却要求,控制组织,达标力学性能的要求。

Claims (4)

1.一种基于尝试法的棒材水箱反向自决策控温智能计算方法,其特征在于,所述的方法包括下述步骤:
1)根据输入的钢种、规格,以及穿水前工艺参数和布置,计算出终轧机后轧件温度分布;
2)计算得出不同水箱的水量阀在满刻度或零刻度的组合下轧件上冷床时温度内外分布,及供水系统载荷;
3)以表层温度值和温度内外差值为目标函数,将计算得出的温度分布进行对比,查找出不越界的水箱开关组合;
4)尝试按预定刻度值,以供水系统最小载荷为目标函数,置入随机±5%~10%的开口度刻度值,逐一调节每个水箱水量阀大小;
5)计算每个水箱调低水量阀刻度后的轧件整体的时间-温度演变曲线,与预期的满足组织性能的cct冷却路径,进行比对,满足要求则结束退出,否则返回步骤4。
2.如权利要求1所述的基于尝试法的棒材水箱反向自决策控温智能计算方法,其特征在于,所述的方法包括下述步骤:所述的轧件温度分布是基于傅里叶传热理论的有限差分算法计算得出的。
3.如权利要求1所述的基于尝试法的棒材水箱反向自决策控温智能计算方法,其特征在于,以表层温度值、温度内外差值、供水系统最小载荷为目标函数,所述的按预定刻度值逐一调整每个水箱水量阀的大小的步骤具体为:第1~第j-1个水箱置入随机±5%~10%的开口度刻度值,按预定刻度值调低第j个水箱调节水量阀的大小。
4.如权利要求1所述的基于尝试法的棒材水箱反向自决策控温智能计算方法,其特征在于,还包括调节预定刻度值的步骤。
CN201911125953.XA 2019-11-15 2019-11-15 基于尝试法的棒材水箱反向自决策控温智能计算方法 Active CN110814053B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911125953.XA CN110814053B (zh) 2019-11-15 2019-11-15 基于尝试法的棒材水箱反向自决策控温智能计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911125953.XA CN110814053B (zh) 2019-11-15 2019-11-15 基于尝试法的棒材水箱反向自决策控温智能计算方法

Publications (2)

Publication Number Publication Date
CN110814053A true CN110814053A (zh) 2020-02-21
CN110814053B CN110814053B (zh) 2021-04-20

Family

ID=69556157

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911125953.XA Active CN110814053B (zh) 2019-11-15 2019-11-15 基于尝试法的棒材水箱反向自决策控温智能计算方法

Country Status (1)

Country Link
CN (1) CN110814053B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112090968A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090970A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090969A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制方法和系统
CN112329089A (zh) * 2020-09-29 2021-02-05 中冶华天工程技术有限公司 基于轧制程序表的热轧棒线材全流程温度场计算系统
CN113751512A (zh) * 2021-08-13 2021-12-07 福建三钢闽光股份有限公司 一种高线不冷段长度的测量方法及其装置
CN114951297A (zh) * 2022-06-29 2022-08-30 中冶华天工程技术有限公司 一种热轧棒材冷床物料跟踪及产品质量追溯方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101639691A (zh) * 2007-03-27 2010-02-03 宝山钢铁股份有限公司 钢板冷却控制模拟方法
CN101642780A (zh) * 2009-06-05 2010-02-10 中冶赛迪工程技术股份有限公司 一种中间坯冷却系统及冷却控制工艺
CN103464474A (zh) * 2013-08-12 2013-12-25 北京首钢自动化信息技术有限公司 棒材控冷系统的温度自动控制方法
CN105327949A (zh) * 2015-11-25 2016-02-17 北京金自天正智能控制股份有限公司 一种热轧带钢卷取温度的流量控制方法
CN107052061A (zh) * 2017-06-16 2017-08-18 安徽工程大学 一种考虑头尾温度偏差的中厚板及宽厚板层流冷却的水量动态调节方法
JP2017170488A (ja) * 2016-03-24 2017-09-28 株式会社神戸製鋼所 厚鋼板冷却方法及び厚鋼板冷却装置
CN107971345A (zh) * 2016-10-25 2018-05-01 宝山钢铁股份有限公司 钢管在线冷却的过程控制系统和控制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101639691A (zh) * 2007-03-27 2010-02-03 宝山钢铁股份有限公司 钢板冷却控制模拟方法
CN101642780A (zh) * 2009-06-05 2010-02-10 中冶赛迪工程技术股份有限公司 一种中间坯冷却系统及冷却控制工艺
CN103464474A (zh) * 2013-08-12 2013-12-25 北京首钢自动化信息技术有限公司 棒材控冷系统的温度自动控制方法
CN105327949A (zh) * 2015-11-25 2016-02-17 北京金自天正智能控制股份有限公司 一种热轧带钢卷取温度的流量控制方法
JP2017170488A (ja) * 2016-03-24 2017-09-28 株式会社神戸製鋼所 厚鋼板冷却方法及び厚鋼板冷却装置
CN107971345A (zh) * 2016-10-25 2018-05-01 宝山钢铁股份有限公司 钢管在线冷却的过程控制系统和控制方法
CN107052061A (zh) * 2017-06-16 2017-08-18 安徽工程大学 一种考虑头尾温度偏差的中厚板及宽厚板层流冷却的水量动态调节方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李振垒: "基于超快速冷却的热轧带钢轧后冷却控制系统与策略研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
李旭东: "超快冷条件下热轧板带钢轧后冷却控制方法研究与应用", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112090968A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090970A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090969A (zh) * 2020-08-28 2020-12-18 中冶华天工程技术有限公司 长材轧制穿水冷却控制方法和系统
CN112090968B (zh) * 2020-08-28 2022-03-25 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090970B (zh) * 2020-08-28 2022-04-22 中冶华天工程技术有限公司 长材轧制穿水冷却控制系统与控制方法
CN112090969B (zh) * 2020-08-28 2022-04-22 中冶华天工程技术有限公司 长材轧制穿水冷却控制方法和系统
CN112329089A (zh) * 2020-09-29 2021-02-05 中冶华天工程技术有限公司 基于轧制程序表的热轧棒线材全流程温度场计算系统
CN112329089B (zh) * 2020-09-29 2022-08-09 中冶华天工程技术有限公司 基于轧制程序表的热轧棒线材全流程温度场计算系统
CN113751512A (zh) * 2021-08-13 2021-12-07 福建三钢闽光股份有限公司 一种高线不冷段长度的测量方法及其装置
CN114951297A (zh) * 2022-06-29 2022-08-30 中冶华天工程技术有限公司 一种热轧棒材冷床物料跟踪及产品质量追溯方法

Also Published As

Publication number Publication date
CN110814053B (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN110814053B (zh) 基于尝试法的棒材水箱反向自决策控温智能计算方法
CN106391734B (zh) 一种用于热轧高强度钢筋的分段阶梯型冷却装置
CN109201748B (zh) 一种螺纹钢细晶轧制轧件冷却智能控制方法
CN109226327B (zh) 一种2.25Cr1Mo0.25V钢大型筒体锻件锻轧结合制造方法
CN105618491B (zh) 一种用于钢板在线固溶的钢板生产方法、装置及控制系统
CN111046457B (zh) 一种施工全周期混凝土拱坝温控曲线模型
CN102397875A (zh) 一种纵向变厚度钢板的生产方法
CN101168797B (zh) 一种优碳钢材在线预测与控制组织索氏体化方法
CN109108094B (zh) 一种螺纹钢细晶轧制智能控制方法
CN103611734A (zh) 一种层流冷却温度控制方法及系统
CN104988409A (zh) 免退火冷镦钢热轧盘条及其生产方法
CN110918642A (zh) 一种改善42CrMo棒材弯曲度的生产控制方法
CN103882221A (zh) 热轧加热炉模型中动态预测板坯最短在炉时间的方法
CN203991658U (zh) 一种高强度钢筋生产系统
CN109396759B (zh) 一种深潜用980钢大直径筒体锻轧结合制造方法
CN108907131B (zh) 一种降低板坯连铸头尾坯表面裂纹的二冷控制方法
CN107016509B (zh) 一种降低轧钢工序吨钢能耗的方法
CN110883105B (zh) 热轧螺纹钢轧后控冷寻优水箱参数实时调整的方法
CN113245376A (zh) 一种摩根五代轧机配套83m斯太尔摩辊道生产高品质焊丝钢的控轧控冷方法
CN101125343A (zh) 普通四辊热带钢连轧机带钢边部增厚综合控制方法
CN102817038B (zh) 一种在线调质、加热发黑涂装工艺生产高强螺纹钢筋的方法
CN104745790A (zh) 一种冷镦钢的软化轧制方法
Yuan et al. Control and application of cooling path after rolling for hot strip based on ultra fast cooling
RU2563911C2 (ru) Способ производства рулонного проката на непрерывном широкополосном стане
CN114798762B (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