CN106891896A - 一种商用车机械式自动变速器载荷识别计算方法 - Google Patents

一种商用车机械式自动变速器载荷识别计算方法 Download PDF

Info

Publication number
CN106891896A
CN106891896A CN201710186860.2A CN201710186860A CN106891896A CN 106891896 A CN106891896 A CN 106891896A CN 201710186860 A CN201710186860 A CN 201710186860A CN 106891896 A CN106891896 A CN 106891896A
Authority
CN
China
Prior art keywords
complete vehicle
gearshift
vehicle quality
acceleration
quality
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
CN201710186860.2A
Other languages
English (en)
Other versions
CN106891896B (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.)
FAW Group Corp
Original Assignee
FAW Group Corp
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 FAW Group Corp filed Critical FAW Group Corp
Priority to CN201710186860.2A priority Critical patent/CN106891896B/zh
Publication of CN106891896A publication Critical patent/CN106891896A/zh
Application granted granted Critical
Publication of CN106891896B publication Critical patent/CN106891896B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/12Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to parameters of the vehicle itself, e.g. tyre models
    • B60W40/13Load or weight

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Transmission Device (AREA)

Abstract

本发明公开了一种商用车机械式自动变速器载荷识别计算方法,包括:判断换挡状态;确定换挡后的加速度、驱动力、空气阻力和坡度值;计算整车质量;判断计算的整车质量是否满足存储条件;如果满足存储条件,则进行存储;存储满足要求的整车质量;判断存储的整车质量是否超过存储限值,如果超过,则对满足条件的整车质量取平均值。本发明通过在换挡之后计算加速度,并对加速度进行滤波处理,以及利用坡道传感器获得的坡度值来计算坡道阻力,从而保证了整车质量计算的准确性,并且将不满足存储条件的整车质量进行舍弃,进一步保证了整车质量计算值的准确性,从而能够准确地对车辆换挡过程中的换挡力进行控制。

Description

一种商用车机械式自动变速器载荷识别计算方法
技术领域
本发明属于商用车自动机械式自动变速器(AMT)领域,涉及一种载荷识别方法,具体涉及一种商用车机械式自动变速器载荷识别计算方法。
背景技术
在商用车自动变速器领域,特别在重型卡车自动变速器领域,AMT自动变速器逐渐成为商用车自动变速器的主要形式。在国外,商用车AMT变速器技术已经成熟,各大汽车厂商先后推出了自己的AMT变速器产品。在国内,各汽车厂商也纷纷加快AMT变速器产品研发步伐,力争突破商用车AMT产品研发关键技术。AMT变速器具有传统手动变速器的传动效率高、制造成本低、容易进行维修的特点。由于AMT变速器在换挡过程自动完成,可减小换挡冲击,提高换挡舒适性,进一步降低驾驶员劳动强度等优点,AMT变速器在商用车自动变速器领域具备较好的发展前景。
在商用车AMT自动变速器领域,车辆的载荷对于重型卡车有着至关重要的影响,对于起步控制中离合器结合深度,换挡过程中换挡力的控制,以及换挡时机的选择都有着很大的影响。比如在坡路环境下,重载起步时,平路的离合器控制就不能满足坡路起步控制,坡路换挡控制及换挡点的选择也不能满足需求。因此,需要知道载荷,在不同的环境下采取不同的控制方式。
发明内容
针对上述技术问题,本发明提供一种基于已有的车辆信息来进行商用车机械式自动变速器载荷识别的方法。
为解决上述技术问题,本发明采用的技术方案为:
本发明实施例提供了一种商用车机械式自动变速器载荷识别计算方法,包括:
S100:在车辆的行驶过程中,基于设定的时间间隔获取离合器位置、发动机输出轴和变速箱输入轴的转速,并基于获取的结果来判断变速器是否处于换挡状态,当判定为处于换挡状态时,进入步骤S200;
S200:在换挡之后,利用换挡前的车速和换挡后的车速确定车辆加速度,并确定当前车轮上的驱动力Ft、车辆受到的空气阻力Fw和车辆所处的坡度值α;
S300:基于确定的加速度、驱动力、空气阻力和坡度值来计算整车质量;
S400:判断计算的整车质量是否满足存储条件,如果满足,则执行步骤S500,如果不满足,则不存储,所述存储条件根据车辆的当前档位、油门开度和加速度来确定;
S500:存储满足存储条件的整车质量,并进入步骤S600;
S600:判断存储的整车质量的个数是否大于存储限值,如果不大于,则返回步骤S100,如果大于,则进入步骤S700;
S700:对满足整车质量最大最小值条件的所有整车质量值进行取平均值处理,以得到最终的整车质量估算值。
可选地,所述整车质量基于下述公式(1)确定:
其中,m为整车质量,Ft为车轮上的驱动力,Fw为车轮受到的空气阻力,α为坡度值,λ为旋转惯性系数,fR为滚动阻力系数,a为加速度。
可选地,所述加速度a通过下述公式(2)确定:
其中,VBef为换挡前的车速,VAft为换挡后的车速,t为换挡时间,单位为s。
可选地,所述驱动力Ft基于下述公式(3)确定:
其中,Ttq为发动机扭矩,iA为传动系传动比,η为传动系总效率,r为轮胎半径。
可选地,所述驱动力Fw基于下述公式(4)确定:
其中,CD为空气阻力系数,A为车辆的迎风面积,ρ为空气密度,V为车辆的当前行驶速度。
可选地,在计算整车质量时,对所述加速度进行滤波处理。
可选地,所述存储条件为落入基于车辆的当前的档位、油门开度和加速度通过查表确定的质量最大值和质量最小值之间。
与现有技术相比,本发明的有益效果为:
(1)利用离合器位置及发动机输出轴与变速器输入轴转速判断变速器是否处于换挡状态,在判定为换挡状态时,由变速箱控制系统准确获取换挡前后的车速用于加速度,并进行滤波处理,避免了加速度计算不准确对计算结果的影响;
(2)根据不同档位的牵引力及加速度限定质量最大值及最小值,保证质量计算值准确性;
(3)利用坡道传感器获取坡度值,增加了坡道阻力的计算准确性。
附图说明
图1为本发明的商用车机械式自动变速器载荷识别计算方法流程图。
图2为本发明的车辆在坡道上的受力分析图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明的商用车机械式自动变速器载荷识别计算通过安置在车辆中的换挡控制器进行计算,该换挡控制器可采用具有较强运算能力的控制芯片,例如瑞萨M12386,该换挡控制器与发动机控制单元、ABS控制单元通过CAN通信进行数据交互。换挡控制器中事先存储有与车辆相关的整车参数,包括车辆的传动系传动比iA,轮胎半径r,传动系总效率η,空气阻力系数CD,空气密度ρ,迎风面积A,滚动阻力系数fR,旋转惯性系数λ等。换挡控制器周期性获取离合器位置信息,例如以10ms为控制周期来周期性获取发动机输入轴转速信息、变速箱输入轴转速信息以及换挡前后车速信息等,并根据获取的离合器位置信息,发动机输入轴转速信息、变速箱输入轴转速信息判断变速器是否处于换挡状态,在判定为换挡状态后,执行整车质量计算。具体的计算过程如图1所示。
如图1所示,本发明实施例提供的商用车机械式自动变速器载荷识别计算方法包括:
S100:判断换挡状态;
S200:确定换挡后的加速度、驱动力、空气阻力和坡度值;
S300:计算整车质量;
S400:判断计算的整车质量是否满足存储条件;如果满足存储条件,则执行步骤S500,如果不满足,则不进行存储;
S500:存储满足要求的整车质量;
S600:判断存储的整车质量是否超过存储限值,如果没有超过,则执行步骤S100,如果超过,则执行步骤S700;
S700:对满足条件的整车质量取平均值。
具体地,上述各步骤具体为:
S100:判断换挡状态
在车辆的行驶过程中,换挡控制器基于设定的时间间隔,例如10ms获取离合器位置、发动机输出轴和变速箱输入轴的转速,并基于获取的结果来判断变速器是否处于换挡状态,当判定为处于换挡状态时,进入步骤S200,当判定为不处于换挡状态时,则继续采集数据进行判断。
S200:确定换挡后的加速度、驱动力、空气阻力和坡度值
具体地,在换挡之后,换挡控制器利用换挡前的车速和换挡后的车速确定车辆加速度,确定当前车轮上的驱动力Ft、车辆受到的空气阻力Fw和车辆所处的坡度值α。
其中,所述加速度a可通过下述公式(1)确定:
其中,VBef为换挡前的车速,具体为换挡开始时那一时刻的速度,VAft为换挡后的车速,具体为换挡结束后车速稳定时那一时刻的车速,t为换挡时间,单位为s,该换挡时间t是换挡变速器中存储的累积的换挡时间,即从开始换挡到换挡结束后车速稳定所经历的时间。
所述驱动力Ft可基于下述公式(2)确定:
其中,Ttq为换挡控制器从发动机控制单元获取的当前发动机扭矩,iA为传动系传动比,η为传动系总效率,r为轮胎半径。
所述驱动力Fw可基于下述公式(3)确定:
其中,CD为空气阻力系数,A为车辆的迎风面积,ρ为空气密度,V为换挡控制器获取的当前的车辆行驶速度,在本发明的应用中,为换挡后的车速。
S300:计算整车质量
基于步骤S200确定的加速度、驱动力、空气阻力和坡度值来计算整车质量。
图2为本发明的车辆在坡道上的受力分析图。如图2所示,整车质量可基于下述公式(4)确定:
其中,m为整车质量,Ft为车轮上的驱动力,Fw为车轮受到的空气阻力,α为由坡度传感器获取的车辆当前坡度值,λ为旋转惯性系数,fR为滚动阻力系数,a为加速度,为避免计算结果不准确,在计算过程中,对加速度进行滤波处理,例如利用卡尔曼滤波函数进行滤波处理。
S400:判断计算的整车质量是否满足存储条件
具体地,在计算得到整车质量时,换挡传感器会判断计算的整车质量是否满足存储条件,如果满足,则进行存储,如果不满足,则不存储。本发明中的存储条件为落入基于车辆的当前的档位、油门开度和加速度通过查表确定的质量最大值和质量最小值之间,对于特定的车型,整车质量的最大值和最小值与车辆的档位、油门开度和加速度之间的关系是确定的,存储在预先形成的关系表中,因此,当知道车辆的档位、油门开度和加速度信息时,便可通过查表法来确定整车质量的最大值和最小值。
S500:存储满足要求的整车质量
当计算的整车质量落入通过查表所确定的整车质量最大值和最小值之间时,则对该整车质量进行存储。
S600:判断存储的整车质量是否超过存储限值
换挡控制器会实时判断存储的整车质量的个数是否大于存储限值,如果不大于,则返回步骤S100继续计算整车质量,如果大于,则进入步骤S700;在本发明的一示例中,存储限值为15次,但并不局限于此,也可根据实际情况来设置其他的存储限值。
S700:整车质量平均值处理
当存储数量达到存储限值时,即已经存储了足够的整车质量数,则换挡控制器对满足整车质量最大最小值条件的所有整车质量值进行取平均值处理,以得到最终的整车质量估算值。
综上,本发明通过在换挡之后计算加速度,并对加速度进行滤波处理,以及利用坡道传感器获得的坡度值来计算坡道阻力,从而保证了整车质量计算的准确性,并且将不满足存储条件的整车质量进行舍弃,进一步保证了整车质量计算值的准确性,从而能够准确地对车辆换挡过程中的换挡力进行控制。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (7)

1.一种商用车机械式自动变速器载荷识别计算方法,其特征在于,包括:
S100:在车辆的行驶过程中,基于设定的时间间隔获取离合器位置、发动机输出轴和变速箱输入轴的转速,并基于获取的结果来判断变速器是否处于换挡状态,当判定为处于换挡状态时,进入步骤S200;
S200:在换挡之后,利用换挡前的车速和换挡后的车速确定车辆加速度,并确定当前车轮上的驱动力Ft、车辆受到的空气阻力Fw和车辆所处的坡度值α;
S300:基于确定的加速度、驱动力、空气阻力和坡度值来计算整车质量;
S400:判断计算的整车质量是否满足存储条件,如果满足,则执行步骤S500,如果不满足,则不存储,所述存储条件根据车辆的当前档位、油门开度和加速度来确定;
S500:存储满足存储条件的整车质量,并进入步骤S600;
S600:判断存储的整车质量的个数是否大于存储限值,如果不大于,则返回步骤S100,如果大于,则进入步骤S700;
S700:对满足整车质量最大最小值条件的所有整车质量值进行取平均值处理,以得到最终的整车质量估算值。
2.根据权利要求1所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,所述整车质量基于下述公式(1)确定:
m = F t - F w gf R c o s α + g sin α + λ a - - - ( 1 )
其中,m为整车质量,Ft为车轮上的驱动力,Fw为车轮受到的空气阻力,α为坡度值,λ为旋转惯性系数,fR为滚动阻力系数,a为加速度。
3.根据权利要求2所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,所述加速度a通过下述公式(2)确定:
a = V A f t - V B e f t - - - ( 2 )
其中,VBef为换挡前的车速,VAft为换挡后的车速,t为换挡时间,单位为s。
4.根据权利要求2所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,所述驱动力Ft基于下述公式(3)确定:
F t = T t q * i A * η r - - - ( 3 )
其中,Ttq为发动机扭矩,iA为传动系传动比,η为传动系总效率,r为轮胎半径。
5.根据权利要求2所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,所述驱动力Fw基于下述公式(4)确定:
F w = 1 2 * C D * A * ρ * V 2 - - - ( 4 )
其中,CD为空气阻力系数,A为车辆的迎风面积,ρ为空气密度,V为车辆的当前行驶速度。
6.根据权利要求3所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,在计算整车质量时,对所述加速度进行滤波处理。
7.根据权利要求1所述的商用车机械式自动变速器载荷识别计算方法,其特征在于,所述存储条件为落入基于车辆的当前的档位、油门开度和加速度通过查表确定的质量最大值和质量最小值之间。
CN201710186860.2A 2017-03-27 2017-03-27 一种商用车机械式自动变速器载荷识别计算方法 Active CN106891896B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710186860.2A CN106891896B (zh) 2017-03-27 2017-03-27 一种商用车机械式自动变速器载荷识别计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710186860.2A CN106891896B (zh) 2017-03-27 2017-03-27 一种商用车机械式自动变速器载荷识别计算方法

Publications (2)

Publication Number Publication Date
CN106891896A true CN106891896A (zh) 2017-06-27
CN106891896B CN106891896B (zh) 2019-06-07

Family

ID=59193244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710186860.2A Active CN106891896B (zh) 2017-03-27 2017-03-27 一种商用车机械式自动变速器载荷识别计算方法

Country Status (1)

Country Link
CN (1) CN106891896B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108128301A (zh) * 2017-12-21 2018-06-08 潍柴动力股份有限公司 一种发动机功率调节方法、装置及电子设备
CN110356409A (zh) * 2018-03-26 2019-10-22 北汽福田汽车股份有限公司 整车质量检测方法和装置、车辆
CN111267638A (zh) * 2020-01-21 2020-06-12 浙江吉利新能源商用车集团有限公司 一种商用车山路工况驾驶的控制方法、系统、设备及介质
CN111717214A (zh) * 2019-03-22 2020-09-29 长沙智能驾驶研究院有限公司 整车质量估算方法、装置、电子设备和存储介质
CN112441016A (zh) * 2019-09-05 2021-03-05 百度(美国)有限责任公司 基于齿轮的车辆载荷推断系统
CN112590804A (zh) * 2020-12-23 2021-04-02 苏州挚途科技有限公司 无人车辆的载荷轨迹跟踪方法和装置
CN112918553A (zh) * 2021-03-12 2021-06-08 合肥工业大学 纵向坡道车辆自适应转向控制方法及装置
CN113085876A (zh) * 2021-03-23 2021-07-09 浙江吉利控股集团有限公司 一种基于递推梯度矫正法的整车质量估算方法及装置
CN113119727A (zh) * 2021-05-07 2021-07-16 恒大新能源汽车投资控股集团有限公司 车辆超载处理方法及装置
CN113696967A (zh) * 2021-09-01 2021-11-26 国汽智控(北京)科技有限公司 车辆转向控制方法、装置、电子设备及存储介质
CN116872947A (zh) * 2023-09-08 2023-10-13 江西五十铃汽车有限公司 一种车辆载重测量方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004058960A (ja) * 2002-07-31 2004-02-26 Nissan Motor Co Ltd 接地荷重推定装置及び車両挙動制御装置
CN102627108A (zh) * 2012-04-11 2012-08-08 清华大学 一种基于高频信息提取的整车质量估算方法
CN103129560A (zh) * 2011-11-30 2013-06-05 通用汽车环球科技运作有限责任公司 用于估计车辆的质量的系统和方法
CN104848923A (zh) * 2014-05-05 2015-08-19 北汽福田汽车股份有限公司 车辆载荷的检测方法和检测装置
CN105209309A (zh) * 2013-05-24 2015-12-30 威伯科有限责任公司 用于确定车辆重量的方法和装置以及具有这种装置的车辆

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004058960A (ja) * 2002-07-31 2004-02-26 Nissan Motor Co Ltd 接地荷重推定装置及び車両挙動制御装置
CN103129560A (zh) * 2011-11-30 2013-06-05 通用汽车环球科技运作有限责任公司 用于估计车辆的质量的系统和方法
CN102627108A (zh) * 2012-04-11 2012-08-08 清华大学 一种基于高频信息提取的整车质量估算方法
CN105209309A (zh) * 2013-05-24 2015-12-30 威伯科有限责任公司 用于确定车辆重量的方法和装置以及具有这种装置的车辆
CN104848923A (zh) * 2014-05-05 2015-08-19 北汽福田汽车股份有限公司 车辆载荷的检测方法和检测装置

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108128301A (zh) * 2017-12-21 2018-06-08 潍柴动力股份有限公司 一种发动机功率调节方法、装置及电子设备
CN110356409A (zh) * 2018-03-26 2019-10-22 北汽福田汽车股份有限公司 整车质量检测方法和装置、车辆
CN111717214A (zh) * 2019-03-22 2020-09-29 长沙智能驾驶研究院有限公司 整车质量估算方法、装置、电子设备和存储介质
CN112441016A (zh) * 2019-09-05 2021-03-05 百度(美国)有限责任公司 基于齿轮的车辆载荷推断系统
CN111267638A (zh) * 2020-01-21 2020-06-12 浙江吉利新能源商用车集团有限公司 一种商用车山路工况驾驶的控制方法、系统、设备及介质
CN112590804A (zh) * 2020-12-23 2021-04-02 苏州挚途科技有限公司 无人车辆的载荷轨迹跟踪方法和装置
CN112918553A (zh) * 2021-03-12 2021-06-08 合肥工业大学 纵向坡道车辆自适应转向控制方法及装置
CN113085876A (zh) * 2021-03-23 2021-07-09 浙江吉利控股集团有限公司 一种基于递推梯度矫正法的整车质量估算方法及装置
CN113119727A (zh) * 2021-05-07 2021-07-16 恒大新能源汽车投资控股集团有限公司 车辆超载处理方法及装置
CN113696967A (zh) * 2021-09-01 2021-11-26 国汽智控(北京)科技有限公司 车辆转向控制方法、装置、电子设备及存储介质
CN116872947A (zh) * 2023-09-08 2023-10-13 江西五十铃汽车有限公司 一种车辆载重测量方法及系统
CN116872947B (zh) * 2023-09-08 2024-01-26 江西五十铃汽车有限公司 一种车辆载重测量方法及系统

Also Published As

Publication number Publication date
CN106891896B (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
CN106891896A (zh) 一种商用车机械式自动变速器载荷识别计算方法
CN109357003A (zh) 一种商用车机械式自动变速器起步档位输出控制方法
US7945365B2 (en) Method for influencing an automatic gearbox, taking into consideration the tractional resistance
US8872645B2 (en) Method for road grade estimation for enhancing the fuel economy index calculation
EP2236862B1 (en) Use of gradient road resistance strategies
US6633006B1 (en) Method and device for determining the mass of a vehicle
CN104937315B (zh) 车辆及用于提供车辆参数值的方法
EP0666435B1 (en) Method/system to determine gross combination weight of vehicles
US8370032B2 (en) Systems and methods for shift control for vehicular transmission
JP3619970B2 (ja) 車両の車両連結総重量を決定するための装置および方法
EP2359032B1 (en) Gear feedback system
US8676459B2 (en) Powertrain for a vehcile and system and method for controlling the powertrain
RU2487027C2 (ru) Устройство обеспечения маневрирования на склоне для транспортного средства
US10948074B2 (en) Methods for predictive shifting
CN104512404B (zh) 确定短期驾驶倾向的方法及利用该方法控制换挡的系统
EP2344786B1 (en) Method and system for controlling a vehicle powertrain
CN107117178A (zh) 考虑换挡和道路坡度因素的车辆质量估计方法
CN109723813B (zh) 一种目标挡位计算方法和装置
CN113339497B (zh) 自动变速箱的换挡规律的确定方法、相关装置及存储介质
EP2085657B1 (en) Method for controlling pickup on uphill for a motor-vehicle provided with an automatic or robotized gearbox
CN105046050B (zh) 一种基于汽车动力学的动态坡度计算方法
CN106143494A (zh) 使用重力加速度传感器估算道路坡度的装置和方法
CN106564504A (zh) 一种用于amt整车质量的计算方法
CN111252072A (zh) 针对纯电动汽车实现坡道检测处理的方法
JP2001208180A (ja) 自動車のドライブトレインに配置される自動化されたトランスミッションのための変速比をもとめるための方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant