CN106017759A - 用转动惯量参数修正发动机转矩测量误差的方法 - Google Patents

用转动惯量参数修正发动机转矩测量误差的方法 Download PDF

Info

Publication number
CN106017759A
CN106017759A CN201610438784.5A CN201610438784A CN106017759A CN 106017759 A CN106017759 A CN 106017759A CN 201610438784 A CN201610438784 A CN 201610438784A CN 106017759 A CN106017759 A CN 106017759A
Authority
CN
China
Prior art keywords
flywheel
grades
inertial
level
rotary inertia
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
CN201610438784.5A
Other languages
English (en)
Other versions
CN106017759B (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.)
Guangxi Normal University
Original Assignee
Guangxi Normal University
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 Guangxi Normal University filed Critical Guangxi Normal University
Priority to CN201610438784.5A priority Critical patent/CN106017759B/zh
Publication of CN106017759A publication Critical patent/CN106017759A/zh
Application granted granted Critical
Publication of CN106017759B publication Critical patent/CN106017759B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L3/00Measuring torque, work, mechanical power, or mechanical efficiency, in general
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L3/00Measuring torque, work, mechanical power, or mechanical efficiency, in general
    • G01L3/24Devices for determining the value of power, e.g. by measuring and simultaneously multiplying the values of torque and revolutions per unit of time, by multiplying the values of tractive or propulsive force and velocity
    • G01L3/242Devices for determining the value of power, e.g. by measuring and simultaneously multiplying the values of torque and revolutions per unit of time, by multiplying the values of tractive or propulsive force and velocity by measuring and simultaneously multiplying torque and velocity

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明公开了用转动惯量参数修正发动机转矩测量误差的方法,包括的步骤有:在发动机曲轴后部串联设置一级和二级惯性飞轮;测量发动机、一级和二级惯性飞轮的角速度;计算发动机、一级和二级惯性飞轮的角加速度;利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩;根据测量工况需要设定发动机曲轴的转动惯量测量误差阀值,判断转动惯量的测定误差是否符合误差阀值要求;用优化后的转动惯量参数测定发动机的转矩。该方法具有原理可靠、工艺简单、耗油量少的优点,可以依据检测需求选择快速测量或精准测量。

Description

用转动惯量参数修正发动机转矩测量误差的方法
技术领域
本发明涉及汽车发动机转矩,具体是基于转动惯量原理用转动惯量参数修正发动机转矩测量误差的方法。
技术背景
发动机是汽车上重要的动力装置,其性能的好坏直接影响汽车的动力性、经济性、可靠性及排放等。发动机的有效功率是曲轴对外输出的功率,是一个综合性评价指标。通过该指标可以定性地确定发动机的技术状况,并定量地获得发动机的动力性。
现有检测发动机有效功率的方法,分为稳态测功和动态测功两种。稳态测功是指发动机在节气门开度一定、转速一定和其它参数保持不变的稳定状态下,在测功器上测定功率的一种方法。常见的测功器有水力测功器、电力测功器和电涡流测功器等。测功器可测出发动机的转速和转矩,然后通过计算得出功率。稳态测功时,不论发动机的工作行程数和形式如何,其有效功率Pe、有效转矩Ttq和转速n均具有下列关系
P e = T t q n 9550 - - - ( 1 ) (1)
动态测功是通过测试加速时间来测定平均功率,有的采用通过测角加速度以确定瞬时功率。
汽车的动力性系指汽车在良好路面上直线行驶时,由汽车受到的纵向外力决定的、所能达到的平均行驶速度。汽车是一种高效率的运输工具,运输效率 之高低在很大程度上取决于汽车的动力性。
F t = T t q i g i 0 η T r - - - ( 2 ) (2)
其中,Ft—驱动力;
r-车轮滚动半径;
Ttq—发动机转矩;
ig—变速器传动比;
i0—主减速器传动比;
ηT—传动系的机械效率。
从式(1)、(2)可以看出,只要准确测定发动机的转矩,就可进一步测算发动机的有效功率和汽车的驱动力。
发明内容
本发明的目的是针对现有技术的不足,而提供一种基于转动惯量原理用转动惯量参数修正发动机转矩测量误差的方法,该方法具有原理可靠、工艺简单、耗油量少、测量速度快、精准度高的优点。
实现本发明的目的的技术方案是:
用转动惯量参数修正发动机转矩测量误差的方法,包括如下步骤:
(1)在发动机曲轴后部串联设置一级和二级惯性飞轮;
(2)测量发动机、一级和二级惯性飞轮的角速度;
(3)计算发动机、一级和二级惯性飞轮的角加速度;
(4)利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩;
(5)根据测量工况需要设定发动机曲轴转动惯量的测量误差阀值,利用测定的二级惯性飞轮与一级惯性飞轮两者的转动惯量的差值,判断转动惯量的测定结果是否的在误差阀值范围内;如果分别利用二级惯性飞轮与一级惯性飞轮测出的转动惯量的差值的绝对值小于等于误差阀值,则满足测量要求进入下一步;如果利用二级惯性飞轮和一级惯性飞轮测出的转动惯量的差值大于等于误差阀值,给二级惯性飞轮减小一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;如果二级惯性飞轮和一级惯性飞轮测出的转动惯量的差值小于等于负误差阀值,给二级惯性飞轮增大一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;
(6)用步骤(5)修正后的转动惯量参数测量发动机转速,进一步计算发动机的转矩。
通过计算得出的发动机转矩及相应工况点的转速可以计算出发动机的功率。
步骤(1)所述在发动机曲轴后部串联设置一级和二级惯性飞轮,是利用发动机缸体内现有的曲轴及与曲轴配装的飞轮,串联一级惯量飞轮,通过一级飞轮电磁离合器与发动机曲轴飞轮连接,二级惯量飞轮通过二级飞轮电磁离合器与一级惯量飞轮连接。
在发明方法中,设步骤(2)中的测量的发动机的角速度为ω1,一级惯量飞轮的角速度为ω2,二级惯量飞轮的角速度为ω3
设步骤(3)中的计算得到发动机的角加速度为β1,一级惯量飞轮的角加速度为β2,二级惯量飞轮的角加速度为β3
步骤(4)所述利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩,先假设测定的发动机曲轴的转动惯量为J1,一级惯量飞轮的转动惯量为J2,二级惯量飞轮的转动惯量为J3
测定的发动机曲轴的转动惯量J1,发动机的转扭矩Ttq(t),具体方法是固定发动机油门在某一测试位置时,测出以下曲线:
发动机曲轴系统转速增长曲线:y1=ω1(t)
一级惯量飞轮转速增长曲线:y2=ω2(t)
二级惯量飞轮转速增长曲线:y3=ω3(t)
发动机的惯量系统的曲线为:
不带惯量飞轮的发动机曲轴转矩曲线:
带一级惯量飞轮的发动机曲轴转矩曲线:
带二级惯量飞轮的发动机曲轴转矩曲线:
因油门固定,发动机转矩稳定,令为角加速度,则:
Ttq(t)=J1β1,Ttq(t)=(J1+J22
J1β1=(J1+J22
J 1 = J 2 β 2 β 1 - β 2
T t q ( t ) = J 2 β 1 β 2 β 1 - β 2
进一步,Ttq(t)=(J1+J22=(J1+J2+J33
以J3=J1+J2为基准初次选用二级惯量飞轮,则可得:β2=2β3
用前面测定的一、二级惯量飞轮转速增长曲线:y2=ω2(t),y3=ω3(t);
进一步可求导出β2、β3
步骤(5)所述根据测量工况需要设定发动机曲轴转动惯量的误差阀值ε,判断转动惯量的测定误差是否在误差阀值范围内;具体是:如果满足测量要求进入步骤(6);如果给二级惯量飞轮J3减小一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;如果给二级惯量飞轮J3增大一个测试单位,返回第(2)步,再测量二级惯性飞轮的角速度,按流程进行。
步骤(6)所述用优化后的转动惯量参数及二级量飞轮J3测量发动机转速, 进一步计算发动机的转矩Ttq(t)=2J3β3,功率由于加载二级惯量飞轮载荷后,发动机转速上升平稳,所以测量结果更精准。
本发明方法是利用转动惯量与角加速的乘积等于转矩的基本原理测定发动机有效转矩,同时进一步测定发动机有效功率的新方法。利用发动机本身曲轴的等效转动惯量,在发动机曲轴后部串联一级和二级惯性飞轮;通过分别测量发动机转速的增长规律、串联一级惯性飞轮的加速时间曲线、串联二级惯性飞轮的加速时间曲线,完成如下四项工作:1、利用发动机自身及加载一级飞轮的两者的角加速度快速测量发动机曲轴系统的转动惯量;2、利用发动机自身及加载一级飞轮两者的角加速度快速测量发动机转矩;3、计算加载一级飞轮及二级惯量飞轮载荷测得的曲轴系统惯性误差大小,调整二级惯性飞轮惯量,进一步提高测量精准度;4、利用加载二级惯量飞轮载荷后,发动机加速运转平稳,准确测量发动机转矩。该方法具有原理可靠、工艺简单、耗油量少的优点,可以依据检测需求选择快速测量或精准测量。
附图说明
图1为本发明方法的流程图;
图2为在发动机曲轴后部串联设置一级和二级惯性飞轮的结构示意图。
图中,1.二级惯性飞轮 2.一级惯性飞轮 3.发动机曲轴飞轮 4.发动机曲轴 5.一级电磁离合器 6.二级电磁离合器。
具体实施方式
结合附图和实施例对本发明内容作进一步的说明,但不是对本发明的限定。
实施例
用转动惯量参数修正发动机转矩测量误差的方法,包括如下步骤:
(1)在发动机曲轴后部串联设置一级和二级惯性飞轮,如图2所示:是利用发动机缸体内现有的曲轴4及与曲轴配装的飞轮3,在发动机曲轴飞轮3后串联一级惯量飞轮2,通过一级飞轮电磁离合器5与发动机曲轴飞轮3连接,二级惯量飞轮1通过二级飞轮电磁离合器6与一级惯量飞轮2连接。
(2)测量发动机、一级和二级惯性飞轮的角速度,设测量的发动机的角速度为ω1,一级惯量飞轮的角速度为ω2,二级惯量飞轮的角速度为ω3
(3)计算发动机、一级和二级惯性飞轮的角加速度,设计算得到发动机的角加速度为β1,一级惯量飞轮的角加速度为β2,二级惯量飞轮的角加速度为β3
(4)利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩,先假设测定的发动机曲轴的转动惯量为J1,一级惯量飞轮的转动惯量为J2,二级惯量飞轮的转动惯量为J3
测定的发动机曲轴的转动惯量J1,发动机的转扭矩Ttq(t)。
具体方法是固定发动机油门在80%开度位置时,测出以下曲线:
发动机曲轴系统转速增长曲线:上标表示油门位置(80%开度),即工况点;下标表示不同惯量系统曲线类型;
一级惯量飞轮转速增长曲线:
二级惯量飞轮转速增长曲线:
发动机的惯量系统的惯量曲线为:
不带惯量飞轮的发动机曲轴转矩曲线: 带一级惯量飞轮的发动机曲轴转矩曲线:
带二级惯量飞轮的发动机曲轴转矩曲线:
因油门固定,发动机转矩稳定,令为角加速度,则:
T t q 8 ( t ) = J 1 β 1 8 , T t q 8 ( t ) = ( J 1 + J 2 ) β 2 8
J 1 β 1 8 = ( J 1 + J 2 ) β 2 8
J 1 = J 2 β 2 8 β 1 8 - β 2 8
T t q 8 ( t ) = J 2 β 1 8 β 2 8 β 1 8 - β 2 8
进一步,
以J3=J1+J2为基准初次选用二级惯量飞轮,则可得:
用前面测定的一、二级惯量飞轮转速增长曲线:
进一步可求导出
(5)根据测量工况需要设定发动机曲轴转动惯量的误差阀值ε,判断转动惯量的测定误差是否在误差阀值范围内;具体是:如果满足测量要求进入步骤(6);如果给二级惯量飞轮J3减小一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;如果 给二级惯量飞轮J3增大一个测试单位,返回第(2)步,再测量二级惯性飞轮的角速度,按流程进行。
(6)用步骤(5)优化后的转动惯量参数及二级量飞轮J3测量发动机转速,进一步计算发动机的转矩功率由于加载二级惯量飞轮载荷后,发动机转速上升平稳,所以测量结果更精准。

Claims (5)

1.用转动惯量参数修正发动机转矩测量误差的方法,其特征在于,包括如下步骤:
(1)在发动机曲轴后部串联设置一级和二级惯性飞轮;
(2)测量发动机、一级和二级惯性飞轮的角速度;
(3)计算发动机、一级和二级惯性飞轮的角加速度;
(4)利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩;
(5)根据测量工况需要设定发动机曲轴转动惯量的测量误差阀值,利用测定的二级惯性飞轮与一级惯性飞轮两者的转动惯量的差值,判断转动惯量的测定结果是否的在误差阀值范围内;如果分别利用二级惯性飞轮与一级惯性飞轮测出的转动惯量的差值的绝对值小于等于误差阀值,则满足测量要求进入下一步;如果利用二级惯性飞轮和一级惯性飞轮测出的转动惯量的差值大于等于误差阀值,给二级惯性飞轮减小一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;如果二级惯性飞轮和一级惯性飞轮测出的转动惯量的差值小于等于负误差阀值,给二级惯性飞轮增大一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;
(6)用步骤(5)修正后的转动惯量参数测量发动机转速,进一步计算发动机的转矩。
2.根据权利要求1所述转动惯量参数修正发动机转矩测量误差的方法,其特征在于:步骤(1)所述在发动机曲轴后部串联设置一级和二级惯性飞轮,是利用发动机缸体内现有的曲轴及与曲轴配装的飞轮,串联一级惯量飞轮,通过一级飞轮电磁离合器与发动机曲轴飞轮连接,二级惯量飞轮通过二级飞轮电磁离合器与一级惯量飞轮连接。
3.根据权利要求1所述用转动惯量参数修正发动机转矩测量误差的方法,其特征在于:设步骤(2)中的测量的发动机的角速度为ω1,一级惯量飞轮的角速度为ω2,二级惯量飞轮的角速度为ω3
设步骤(3)中的计算得到发动机的角加速度为β1,一级惯量飞轮的角加速度为β2,二级惯量飞轮的角加速度为β3
步骤(4)所述利用发动机和一级惯性飞轮两者的角加速度分别测定发动机曲轴的转动惯量和发动机的转动扭矩,先假设测定的发动机曲轴的转动惯量为J1,一级惯量飞轮的转动惯量为J2,二级惯量飞轮的转动惯量为J3
测定的发动机曲轴的转动惯量J1,发动机的转扭矩Ttq(t),具体方法是固定发动机油门在某一测试位置时,测出以下曲线:
发动机曲轴系统转速增长曲线:y1=ω1(t)
一级惯量飞轮转速增长曲线为:y2=ω2(t)
二级惯量飞轮转速增长曲线:y3=ω3(t)
发动机的惯量系统的曲线为:
不带惯量飞轮的发动机曲轴转矩曲线:
带一级惯量飞轮的发动机曲轴转矩曲线:
带二级惯量飞轮的发动机曲轴转矩曲线:
因油门固定,发动机转矩稳定,令为角加速度,则:
Ttq(t)=J1β1,Ttq(t)=(J1+J22
J1β1=(J1+J22
J 1 = J 2 β 2 β 1 - β 2
T t q ( t ) = J 2 β 1 β 2 β 1 - β 2
进一步,Ttq(t)=(J1+J22=(J1+J2+J33
以J3=J1+J2为基准初次选用二级惯量飞轮,则可得:β2=2β3
用前面测定的一、二级惯量飞轮转速增长曲线:y2=ω2(t),y3=ω3(t);进一步可求导出β2、β3
4.根据权利要求1所述用转动惯量参数修正发动机转矩测量误差的方法,其特征在于:步骤(5)所述根据测量工况需要设定发动机曲轴转动惯量的误差阀值ε,判断转动惯量的测定误差是否在误差阀值范围内;具体是:如果满足测量要求进入步骤(6);如果给二级惯量飞轮J3减小一个测试单位,返回步骤(2),再测量二级惯性飞轮的角速度,按流程进行;如果给二级惯量飞轮J3增大一个测试单位,返回第(2)步,再测量二级惯性飞轮的角速度,按流程进行。
5.根据权利要求1所述用转动惯量参数修正发动机转矩测量误差的方法,其特征在于:步骤(6)修正后的转动惯量参数及二级量飞轮J3测定发动机转速,进一步计算发动机的转矩Ttq(t)=2J3β3
CN201610438784.5A 2016-06-17 2016-06-17 用转动惯量参数修正发动机转矩测量误差的方法 Expired - Fee Related CN106017759B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610438784.5A CN106017759B (zh) 2016-06-17 2016-06-17 用转动惯量参数修正发动机转矩测量误差的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610438784.5A CN106017759B (zh) 2016-06-17 2016-06-17 用转动惯量参数修正发动机转矩测量误差的方法

Publications (2)

Publication Number Publication Date
CN106017759A true CN106017759A (zh) 2016-10-12
CN106017759B CN106017759B (zh) 2018-09-07

Family

ID=57088692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610438784.5A Expired - Fee Related CN106017759B (zh) 2016-06-17 2016-06-17 用转动惯量参数修正发动机转矩测量误差的方法

Country Status (1)

Country Link
CN (1) CN106017759B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112229570A (zh) * 2020-10-30 2021-01-15 广西师范大学 基于离心力法的转动惯量误差检测装置及检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1123424A (ja) * 1997-07-07 1999-01-29 Unisia Jecs Corp エンジンの発生トルク推定方法
CN1537221A (zh) * 2001-07-23 2004-10-13 �������֪ʶ��Ȩ�������ι�˾ 发动机转矩计算
CN103261530A (zh) * 2010-12-15 2013-08-21 沃尔沃建造设备有限公司 用于混合动力施工机械的转动控制系统
CN104122035A (zh) * 2014-07-02 2014-10-29 西南大学 一种直接设定型负载转矩和转动惯量模拟系统及其控制方法
CN104318083A (zh) * 2014-10-10 2015-01-28 北京交通大学 一种利用车辆发动机扭矩实时估计离合器传递扭矩的方法
CN104980078A (zh) * 2014-04-04 2015-10-14 广东美的制冷设备有限公司 电机的转动惯量的测量方法及测量装置和电机控制系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1123424A (ja) * 1997-07-07 1999-01-29 Unisia Jecs Corp エンジンの発生トルク推定方法
CN1537221A (zh) * 2001-07-23 2004-10-13 �������֪ʶ��Ȩ�������ι�˾ 发动机转矩计算
CN103261530A (zh) * 2010-12-15 2013-08-21 沃尔沃建造设备有限公司 用于混合动力施工机械的转动控制系统
CN104980078A (zh) * 2014-04-04 2015-10-14 广东美的制冷设备有限公司 电机的转动惯量的测量方法及测量装置和电机控制系统
CN104122035A (zh) * 2014-07-02 2014-10-29 西南大学 一种直接设定型负载转矩和转动惯量模拟系统及其控制方法
CN104318083A (zh) * 2014-10-10 2015-01-28 北京交通大学 一种利用车辆发动机扭矩实时估计离合器传递扭矩的方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112229570A (zh) * 2020-10-30 2021-01-15 广西师范大学 基于离心力法的转动惯量误差检测装置及检测方法

Also Published As

Publication number Publication date
CN106017759B (zh) 2018-09-07

Similar Documents

Publication Publication Date Title
CN102323206B (zh) 一种车辆滚动阻力系数的测量方法
CN105644565B (zh) 一种混合动力汽车负载的测量方法
CN105283749B (zh) 发动机台架系统
CN102305715A (zh) 一种用于汽车动力系统测试的动态负载模拟装置及方法
CN101706300B (zh) 汽车恒速全负荷燃油消耗率检测方法
CN110646223B (zh) 车辆驱动轮瞬时输出功率的测量系统
CN110646224B (zh) 车辆驱动轮瞬间输出功率的测量方法
CN104198196B (zh) 汽车和发动机转动件当量惯量检测方法
CN101936813B (zh) 柴油车加载工况废气排放检测方法
CN106872174A (zh) 汽车变速器台架敲击试验发动机瞬态周期扭矩模拟方法
CN111179465B (zh) 一种汽车油耗预测方法
CN103967631A (zh) 一种柴油机动力性能在线识别装置和识别方法
CN105157991B (zh) 驱动轮表面输出规定功率加载检测方法
CN106017759A (zh) 用转动惯量参数修正发动机转矩测量误差的方法
CN101963519B (zh) 台试模拟路试标准状态燃料消耗量检测方法
CN106052917A (zh) 用角加速度参数优化发动机转矩测量误差的方法
CN105258830B (zh) 汽车驱动轮表面输出功率检测方法
CN1195213C (zh) 增减惯量检测汽车惯量、阻力和功率方法
CN102538889B (zh) 在用车标准状态燃料消耗量检测加载方法
CN102589893B (zh) 柴油车额定功率部分负荷检测方法
CN105181351B (zh) 汽车发动机最大扭矩车速检测方法
Mamala et al. The estimation of the engine power with use of an accelerometer
CN201060085Y (zh) 汽车底盘测功机转动惯量测量装置
Murtagh et al. Development and validation of a forklift truck powertrain simulation
Duoba et al. Investigating Steady-State Road Load Determination Methods for Electrified Vehicles and Coordinated Driving (Platooning)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180907

CF01 Termination of patent right due to non-payment of annual fee