CN104150301B - 电梯运行曲线修正方法 - Google Patents

电梯运行曲线修正方法 Download PDF

Info

Publication number
CN104150301B
CN104150301B CN201410360246.XA CN201410360246A CN104150301B CN 104150301 B CN104150301 B CN 104150301B CN 201410360246 A CN201410360246 A CN 201410360246A CN 104150301 B CN104150301 B CN 104150301B
Authority
CN
China
Prior art keywords
curve
delta
elevator
acceleration
zero
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
Application number
CN201410360246.XA
Other languages
English (en)
Other versions
CN104150301A (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.)
Yungtay Elevator Equipment China Co Ltd
Original Assignee
Yungtay Elevator Equipment China 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 Yungtay Elevator Equipment China Co Ltd filed Critical Yungtay Elevator Equipment China Co Ltd
Priority to CN201410360246.XA priority Critical patent/CN104150301B/zh
Publication of CN104150301A publication Critical patent/CN104150301A/zh
Application granted granted Critical
Publication of CN104150301B publication Critical patent/CN104150301B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Elevator Control (AREA)

Abstract

本发明涉及一种电梯运行曲线修正方法,该方法用于修正电梯进入减减速区后的预估运行曲线,从而解决由于位置的滑移造成电梯停靠不准确的问题。该方法包含以下步骤:测量得到电梯的实际剩余距离,与预估剩余距离进行比较得到偏差剩余距离;提供一种基于三角函数的电梯运行曲线修正方法,修正位移曲线、速度曲线、加速度曲线、加加速度曲线均平滑连续,且曲线方便求解,无需分段计算和存储;将修正曲线与理论预估曲线进行线性叠加,得到修正后的运行曲线。从而实现电梯的平滑停靠,同时最大程度地满足乘客的舒适度。

Description

电梯运行曲线修正方法
技术领域
本发明涉及电梯技术,具体涉及一种基于三角函数的电梯运行曲线修正方法。
技术背景
电梯在运行过程中,钢丝绳由于惯性可能会产生相对驱动轮的滑移,或由于变频器的响应速度而产生的运行速度偏差,造成电梯轿厢的实际走行距离与预估距离出现偏差,导致电梯停靠不准确的问题。
控制电梯运转的变频器输入信号由预估计的电梯运行速度曲线来决定,因此电梯速度曲线的生成和在线优化对电梯的运行性能具有极为重要的意义。通常的电梯运行曲线(如图1)包括加加速区(T1)、匀加速区(T2)、减加速区(T3)、恒速区(T4)、加减速区(T5)、匀减速区(T6)和减减速区(T7),其中加加速区(T1)、减加速区(T3)、加减速区(T5)和减减速区(T7)是其中必备的几个区段。因此,为了解决电梯停靠不准确的问题,则需要对减减速区(T7)区段的速度曲线进行修正,使得电梯能够精确停靠在预定停靠楼层,且停靠时的速度和加速度为0。
中国专利CN102234048A公开了一种电梯速度曲线修正方法。提供一种用于修正该减速后半行程的速度修正曲线;将该速度曲线与速度该修正曲线迭加,以获得修正后速度曲线。该修正速度曲线为平滑连续,但其加速度曲线被设计为4段直线,虽连续但不平滑。从实现的角度考虑,该曲线需要分段存储和计算,实现运行效率较低,无法较好满足利用电梯有限的计算资源实现高计算效率的需求。
然而,电梯的运转曲线与乘客搭乘电梯的舒适性具有密切联系。设计平滑的电梯运转曲线能够有效提高乘客搭乘电梯的舒适性,运转曲线包括位置对时间的曲线、速度对时间的曲线、加速度对时间的曲线、加加速度对时间的曲线。因此,为了最大程度提高乘客搭乘电梯的舒适性,有必要设计一种更加自然的电梯运行修正曲线,该修正的位移曲线、速度曲线、加速度曲线、加加速度曲线均需要平滑连续,同时,曲线应当方便求解,无需分段计算和存储,从而提高电梯有限的计算和存储效率。
考虑到三角函数曲线是工程领域中最常用且最自然的运动曲线之一,其具有较好的连续和平滑性,微分或积分后仍能够保留三角函数特性。与直线或抛物曲线相比,乘客搭乘采用三角函数运转曲线的电梯具有更好的舒适体验。同时,从实现的角度考虑,该类曲线计算简单、内存资源利用较少。
发明内容
本发明提出一种基于三角函数的电梯运行曲线修正方法,修正的位移曲线、速度曲线、加速度曲线、加加速度曲线均平滑连续,使得电梯轿厢能够精确且平顺地停靠在预定楼层,同时提高乘客搭乘电梯的舒适性。
本发明提出的电梯运行曲线修正方法用于修正包含减减速区(T7)段的运行曲线,在时间T7不变的情况下,当该段行程结束时,偏差位移被修正为0(准确停靠),同时修正速度、修正加速度也均为0。
为了实现上述目的,本发明所采用的技术方案如下:
电梯运行曲线修正方法,包括如下步骤:
(1)测量电梯距离终点的实际剩余距离;然后得到实际剩余距离与理论预估剩余距离的差值,即偏差剩余距离ΔS;
(2)当电梯在进入减减速区(T7)后,如果该偏差剩余距离ΔS不为零,则需要进行运行曲线的修正,得到修正如公式(1)的运行曲线ΔΩ(t),其中所述运行曲线ΔΩ(t)包含修正的位移曲线ΔS(t)、速度曲线Δv(t)、加速度曲线Δa(t)和加加速度曲线Δα(t):
ΔΩ(t)=[ΔS(t),Δv(t),Δa(t),Δα(t)]T(1)
当偏差剩余距离ΔS大于零时,运行曲线ΔΩ(t)按如下方法进行修正:
修正的位移曲线ΔS(t)为连续平滑的三角函数曲线,且电梯运行时间结束时,位移偏差为0,即ΔS(T7)=0;
修正的速度曲线Δv(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大再减小至零,即Δv(0)=Δv(T7)=0,当t=T7/2时,修正速度达到最大值Δv,即得到Δv(T7/2)=Δv;
修正的加速度曲线Δa(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大、再减小然后增大至零;当t=T7/2时,修正加速度为零,即Δa(0)=Δa(T7/2)=Δa(T7)=0;
修正的加加速度曲线Δα(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小再增大至零。即Δα(0)=Δα(T7)=0;当t=T7/2时,修正加加速度达到最小值-Δα,即得到Δα(T7/2)=-Δα;
当偏差剩余距离ΔS小于零时,运行曲线ΔΩ(t)的修正方法如下:
修正的位移曲线ΔS(t)为连续平滑的三角函数曲线,且电梯运行时间结束时,位移偏差为0,即ΔS(T7)=0;
修正的速度曲线Δv(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小再减小至零,即Δv(0)=Δv(T7)=0;当t=T7/2时,修正速度达到最小值-Δv,即得到Δv(T7/2)=-Δv;
修正的加速度曲线Δa(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小、再增大然后减小至零;当t=T7/2时,修正加速度为零,即Δa(0)=Δa(T7/2)=Δa(T7)=0;
修正的加加速度曲线Δα(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大再减小至零。即Δα(0)=Δα(T7)=0;当t=T7/2时,修正加加速度达到最大值Δα,即得到Δα(T7/2)=Δα;
(3)将修正曲线ΔΩ(t)与理论预估曲线Ω(t)进行线性叠加,得到修正后的公式(2)运行曲线Ω’(t)。
Ω'(t)=ΔΩ(t)+Ω(t)(2)
采用上述方法使电梯采用修正后的曲线Ω’(t)运行,可实现精确且平顺的停靠,同时提高乘客搭乘电梯的舒适性。
附图说明
图1为通常的电梯预估理想速度曲线示意图。
图2为本发明偏差剩余距离ΔS大于零时,修正前后的速度曲线和加速度曲线示意图。
图3为本发明偏差剩余距离ΔS小于零时,修正前后的速度曲线和加速度曲线示意图。
图4为本发明偏差剩余距离ΔS大于零时的修正曲线,包含修正位移曲线、速度曲线、加速度曲线、加加速度曲线示意图。
图5为本发明偏差剩余距离ΔS小于零时的修正曲线,包含修正位移曲线、速度曲线、加速度曲线、加加速度曲线。
具体实施方式
如下将结合附图和实施例,进一步阐述本发明。
当电梯的运行从匀减速区(T6)转换到减减速区(T7)的瞬间时刻时,此时电梯的速度设为v0,理论预估剩余距离为S0
由于理论预估速度曲线中减减速区(T7)的加加速度为恒定值,设该加加速度为α,则减减速区的初始加速度a0表达为:
a0=αT7(3)
减减速区的初始速度v0和理论预估剩余距离为S0分别表达为:
v 0 = 1 2 α T 7 2 - - - ( 4 )
S 0 = 1 6 α T 7 3 - - - ( 5 )
由公式(4)、(5)式得到减减速区的时间T7
T 7 = 3 S 0 v 0 - - - ( 6 )
由公式(3)、(4)式得到减减速区的初始加速度:
a 0 = 2 v 0 T 7 - - - ( 7 )
由公式(6)、(7)式得到:
a 0 = 2 v 0 2 3 S 0 - - - ( 8 )
因此,由公式(3)、(8)式得到加加速度α与v0、S0的关系:
α = 2 v 0 3 9 S 0 2 - - - ( 9 )
由公式(9)式得到减减速区理论预估位移曲线、速度曲线和加速度曲线分别为:
S ( t ) = 1 6 2 v 0 3 9 S 0 2 ( T 7 - t ) 3 = v 0 3 27 S 0 2 ( T 7 - t ) 3 - - - ( 10 )
v ( t ) = 1 2 2 v 0 3 9 S 0 2 ( T 7 - t ) 2 = v 0 3 9 S 0 2 ( T 7 - t ) 2 - - - ( 11 )
a ( t ) = 2 v 0 3 9 S 0 2 ( T 7 - t ) - - - ( 12 )
然而在电梯运行时由于出现滑移等问题,实际剩余距离与理论预估剩余距离为S0并不相等,此时电梯如果仍采用理论预估曲线运行,在停靠时就会产生偏差。因此,需要对该曲线进行修正,使得电梯能够准确停靠。
通过在每一停靠楼层设置的平层感应器能够测量估计得到实际的剩余距离。设实际剩余距离为S’,则偏差剩余距离ΔS表示为:
ΔS=S'-S0(13)
当偏差剩余距离ΔS大于零时,表明实际剩余距离大于理论预估剩余距离,则应当修正预估运行曲线,例如在运行过程中加速,使得预估剩余距离延长,从而在运行时间结束时,电梯能够准确停靠,并且速度、加速度参数也为零(参见图2、图4)。
当偏差剩余距离ΔS小于零时,表明实际剩余距离小于理论预估剩余距离,则应当修正预估运行曲线,例如在运行过程中减速,使得预估剩余距离缩短,从而在运行时间结束时,电梯能够准确停靠,并且速度、加速度参数也为零(图3、图5)。
因此,通过设计修正运行曲线ΔΩ(t),修正理论预估曲线Ω(t)为Ω’(t),从而使电梯能够准确停靠。
基于三角函数曲线设计的修正运行曲线ΔΩ(t)所包含的修正位移曲线、速度曲线、加速度曲线、加加速度曲线均平滑连续,同时,曲线方便求解,无需分段计算和存储。
图4、图5中修正的速度曲线Δv(t)为连续平滑的三角函数曲线,其中Δv(0)=Δv(T7)=0。由三角函数曲线的性质得到:当t=T7/2时,修正速度达到极值(当偏差剩余距离ΔS大于零时,该值为极大值;当偏差剩余距离ΔS小于零时,该值为极小值),设为Δv,即得到Δv(T7/2)=Δv,且当t=1T7/4=3T7/4时,Δv(1T7/4)=Δv(3T7/4)=Δv/2。由于行走距离为速度曲线所包含的面积,因此可以得到:
ΔS = 1 2 Δv T 7 - - - ( 14 )
得到:
Δv = 2 ΔS T 7 - - - ( 15 )
根据图4、5中所示的修正速度曲线的轨迹,假设曲线表达式为:
Δv(t)=-Acos(ωt)+b(16)
则得到:
A = Δv 2 , b = Δv 2 , ω = 2 π T 7 - - - ( 17 )
因此,修正速度曲线关于时间的函数为:
Δv ( t ) = - Δv 2 cos ( 2 πt T 7 ) + Δv 2 - - - ( 18 )
由Δa(t)=d(Δv(t))得到修正加速度曲线关于时间的函数:
Δa ( t ) = πΔv T 7 sin ( 2 πt T 7 ) - - - ( 19 )
由Δα(t)=d(Δa(t))得到修正加加速度曲线关于时间的函数:
Δα ( t ) = 2 π 2 Δv T 7 2 cos ( 2 πt T 7 ) - - - ( 20 )
通过(18)式计算偏差剩余距离ΔS关于时间的函数:
ΔS ( t ) = ΔS - ∫ Δv ( t ) dt = ΔS - ( Δvt 2 - Δv T 7 4 π sin ( 2 πt T 7 ) ) - - - ( 21 )
将修正曲线ΔΩ(t)与理论预估曲线Ω(t)进行线性叠加,得到修正后的运行曲线Ω’(t):
Ω ' ( t ) = S ' ( t ) v ' ( t ) a ' ( t ) α ' ( t ) = ΔΩ ( t ) + Ω ( t ) = ΔS ( t ) + S ( t ) Δv ( t ) + v ( t ) Δa ( t ) + a ( t ) Δα ( t ) + α ( t ) - - - ( 22 )
上式得到的修正运行曲线ΔΩ’(t)保证了电梯运行的平滑和连续,且该曲线方便求解,无需分段计算和存储。
上述实施例的修正运行曲线包含减减速区的[0,T7]时间段,但修正运行曲线也可以进一步包含[t1,t2]时间段,其中0≤t1<t2≤T7,此时时间间隔变量T7将为t2-t1
上述实施例进行了一次运行曲线的修正,但在电梯运行结束前,也可以进行多次的、不同时间段的运行曲线修正,例如电梯有多组位置校正装置的情况。此时基于三角函数的曲线修正方法具易于叠加等优点,通过相位的变化即能够实现相互叠加。
以上揭示和表述了本发明的基本原理、主要特征和优点。虽然通过实施例对本发明进行了阐述,但并不限制于所揭露的实施例,在不脱离本发明精神和范围的前提下仍有许多变化和改进。因此对于不偏离本发明的变化、改进和修饰均属本发明的涵盖范围。

Claims (3)

1.一种电梯运行曲线修正方法,其特征在于,该方法包含下列步骤:
1].测量电梯距离终点的实际剩余距离;然后得到实际剩余距离与理论预估剩余距离的差值,即偏差剩余距离ΔS;
2].当电梯在进入减减速区后,如果该偏差剩余距离不为零,则修正理论预估曲线中的位移曲线、速度曲线、加速度曲线、加加速度曲线并使位移曲线、速度曲线、加速度曲线、加加速度曲线均平滑连续,从而得到修正曲线;
3].将步骤2]得到的修正曲线与理论预估曲线进行线性叠加,得到修正后的运行曲线;
所述修正曲线包含修正速度曲线、修正加速度曲线、修正加加速度曲线、修正位移曲线,其修正速度曲线、修正加速度曲线、修正加加速度曲线、修正位移曲线关于时间的函数分别表示为如下公式:
&Delta; v ( t ) = - &Delta; v 2 c o s ( 2 &pi; t T 7 ) + &Delta; v 2
&Delta; a ( t ) = &pi; &Delta; v T 7 s i n ( 2 &pi; t T 7 )
&Delta; &alpha; ( t ) = 2 &pi; 2 &Delta; v T 7 2 c o s ( 2 &pi; t T 7 )
&Delta; S ( t ) = &Delta; S - &Integral; &Delta; v ( t ) d t = &Delta; S - ( &Delta; v t 2 - &Delta;vT 7 4 &pi; sin ( 2 &pi; t T 7 ) )
其中:
&Delta; v = 2 &Delta; S T 7 .
2.如权利要求1所述的电梯运行曲线修正方法,其特征在于,在所述步骤2]中,当偏差剩余距离ΔS大于零时,所述位移曲线、速度曲线、加速度曲线、加加速度曲线的修正方法分别如下:
修正的位移曲线ΔS(t)为连续平滑的三角函数曲线,且电梯运行时间结束时(t=T7),位移偏差为0,即ΔS(T7)=0;
修正的速度曲线Δv(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大再减小至零,即Δv(0)=Δv(T7)=0;当t=T7/2时,修正速度达到最大值Δv,即得到Δv(T7/2)=Δv;
修正的加速度曲线Δa(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大、再减小然后增大至零;当t=T7/2时,修正加速度为零,即Δa(0)=Δa(T7/2)=Δa(T7)=0;
修正的加加速度曲线Δα(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小再增大至零;即Δα(0)=Δα(T7)=0,当t=T7/2时,修正加加速度达到最小值-Δα,即得到Δα(T7/2)=-Δα;其中T7为减减速区的[0,T7]时间段。
3.如权利要求1所述的电梯运行曲线修正方法,其特征在于,在所述步骤2]中,当偏差剩余距离ΔS小于零时,所述位移曲线、速度曲线、加速度曲线、加加速度曲线的修正方法分别如下:
修正的位移曲线ΔS(t)为连续平滑的三角函数曲线,且电梯运行时间结束时,位移偏差为0,即ΔS(T7)=0;
修正的速度曲线Δv(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小再减小至零,即Δv(0)=Δv(T7)=0;当t=T7/2时,修正速度达到最小值-Δv,即得到Δv(T7/2)=-Δv;
修正的加速度曲线Δa(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先减小、再增大然后减小至零;当t=T7/2时,修正加速度为零,即Δa(0)=Δa(T7/2)=Δa(T7)=0;
修正的加加速度曲线Δα(t)为连续平滑的三角函数曲线,该曲线随着时间t的增加,从零开始先增大再减小至零;即Δα(0)=Δα(T7)=0;当t=T7/2时,修正加加速度达到最大值Δα,即得到Δα(T7/2)=Δα。
CN201410360246.XA 2014-07-25 2014-07-25 电梯运行曲线修正方法 Active CN104150301B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410360246.XA CN104150301B (zh) 2014-07-25 2014-07-25 电梯运行曲线修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410360246.XA CN104150301B (zh) 2014-07-25 2014-07-25 电梯运行曲线修正方法

Publications (2)

Publication Number Publication Date
CN104150301A CN104150301A (zh) 2014-11-19
CN104150301B true CN104150301B (zh) 2016-05-04

Family

ID=51875973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410360246.XA Active CN104150301B (zh) 2014-07-25 2014-07-25 电梯运行曲线修正方法

Country Status (1)

Country Link
CN (1) CN104150301B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104876079B (zh) * 2015-05-27 2017-01-04 上海贝思特控制技术有限公司 一种基于指数函数的电梯变频器高效电梯运行曲线修正方法
CN107070329B (zh) * 2017-02-20 2019-06-11 武汉科技大学 S型运行曲线计算方法以及步进电机旋转控制方法
CN108287527B (zh) * 2017-07-18 2019-02-01 山东大学 一种基于三角函数的改进s曲线加减速控制方法
CN107807625B (zh) * 2017-09-05 2020-05-12 百度在线网络技术(北京)有限公司 一种基于端到端的自动驾驶***舒适度的评估方法及装置
CN110498308B (zh) * 2019-08-30 2021-11-30 深圳市海浦蒙特科技有限公司 一种电梯运行控制方法和***
CN112850389B (zh) * 2020-12-31 2022-10-25 上海江菱机电有限公司 一种电梯运行速度的控制方法、***及存储介质
CN114212630B (zh) * 2021-11-04 2024-03-15 深圳市海浦蒙特科技有限公司 电梯运行控制方法、装置、电梯及计算机可读存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102234048B (zh) * 2010-04-22 2013-08-21 永大机电工业股份有限公司 电梯速度曲线修正方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3255839B2 (ja) * 1996-02-23 2002-02-12 フジテック株式会社 交流エレベータの制御装置
FI101780B (fi) * 1996-04-30 1998-08-31 Kone Corp Menetelmä ja laitteisto hissin hidastamiseksi
JP4874608B2 (ja) * 2005-09-14 2012-02-15 東芝エレベータ株式会社 エレベータ制御システム

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102234048B (zh) * 2010-04-22 2013-08-21 永大机电工业股份有限公司 电梯速度曲线修正方法

Also Published As

Publication number Publication date
CN104150301A (zh) 2014-11-19

Similar Documents

Publication Publication Date Title
CN104150301B (zh) 电梯运行曲线修正方法
JP4772445B2 (ja) 自動列車運転装置
WO2018126967A1 (zh) 穿梭车速度规划控制方法、装置及穿梭车
US20110046827A1 (en) Apparatus and method for controlling speed in automatic train operation
CN102649438B (zh) 轨道交通列车运行路径的制动点确定方法
CN102854885B (zh) 一种舰载机着舰纵向甲板运动补偿方法
CN110703757B (zh) 一种面向能耗优化的高铁列车速度规划方法
CN106143540A (zh) 轨道交通列车全速密贴运行控制曲线的自适应计算方法
CN102955884B (zh) 一种高速列车跟驰运行全速域安全车距标定方法
JP2016179812A (ja) 列車の移動を制御する方法およびシステム
CN103149937A (zh) 一种基于曲率补偿的横侧向曲线航迹跟踪方法
CN109070765A (zh) 列车控制装置、方法以及程序
CN109765887A (zh) 一种自动驾驶控制方法
US12030541B2 (en) Method and system for state feedback predictive control of high-speed train based on forecast error
CN112078631B (zh) 列车速度控制方法及***
CN107390517A (zh) 用于列车ato***的鲁棒自适应非奇异终端滑模控制方法
CN105005311A (zh) 一种飞行器俯仰通道攻角同步跟踪控制方法
CN102756954B (zh) 一种电梯速度补偿方法
CN104876079A (zh) 一种基于指数函数的电梯变频器高效电梯运行曲线修正方法
CN111580391B (zh) 基于模型预测的动车组牵引力矩控制方法
CN103662092A (zh) 一种衔接主减速和接近段的预测校正方法
CN103818393B (zh) 列车牵引制动***时间滞后特性的补偿方法
CN107515533A (zh) 一种用于列车ato***的鲁棒非奇异终端滑模控制方法
Albrecht et al. Optimal driving strategies for two successive trains on level track subject to a safe separation condition
Baranov et al. Energy-optimal control of vehicle traffic

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant