CN104295387B - 一种基于指示转矩的质调节式发动机起动控制方法 - Google Patents

一种基于指示转矩的质调节式发动机起动控制方法 Download PDF

Info

Publication number
CN104295387B
CN104295387B CN201410401131.0A CN201410401131A CN104295387B CN 104295387 B CN104295387 B CN 104295387B CN 201410401131 A CN201410401131 A CN 201410401131A CN 104295387 B CN104295387 B CN 104295387B
Authority
CN
China
Prior art keywords
starting
trq
angular acceleration
electromotor
fuel injection
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.)
Expired - Fee Related
Application number
CN201410401131.0A
Other languages
English (en)
Other versions
CN104295387A (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201410401131.0A priority Critical patent/CN104295387B/zh
Publication of CN104295387A publication Critical patent/CN104295387A/zh
Application granted granted Critical
Publication of CN104295387B publication Critical patent/CN104295387B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)

Abstract

本发明公开了一种基于指示转矩的质调节式发动机起动控制方法,以解决现有技术难以实现发动机起动过程燃料喷射量精确控制而影响发动机起动性的问题。通过实验获取发动机不同曲轴转速及冷却液温度下的起动拖动转矩;根据发动机的起动控制要求,确定起动过程的目标角加速度函数,并计算加速阻力矩,实现对起动过程中发动机转速和起动时间的控制;根据发动机实际曲轴转速的变化计算发动机实际角加速度,发动机实际角加速度与目标角加速度的偏差确定补偿转矩;起动拖动转矩、加速阻力矩和补偿转矩相加得到起动需求的指示转矩;根据指示效率将需求指示转矩转换为燃料喷射量。本方法可保证各种起动条件下实现起动过程所需的最佳燃料喷射量。

Description

一种基于指示转矩的质调节式发动机起动控制方法
技术领域
本发明涉及一种燃烧发动机的控制,具体涉及基于需求指示转矩并且起动时间可控制的质调节式发动机的起动过程燃料喷射方式的控制策略
背景技术
对于质调节式发动机的起动过程中燃料喷射方式的控制,是通过相对进入气缸的空气量供给最佳燃料喷射量,使其在气缸内形成混合气并燃烧,推动活塞做功,克服发动机内部阻力,驱动曲轴转动,使发动机转速达到预期的目标值。
是否正常起动以及起动过程中发动机转速的变化与设定的燃料喷射量及其喷射方式直接相关。通常起动过程的燃油喷射量是通过发动机台架标定试验来确定的。由于起动过程控制策略不同,起动喷射量以及喷射方式有很大区别。如果燃料喷射量过小,喷射方式不合理,就不能完成正常的起动过程;否则燃料喷射量过多不仅燃料消耗量增加,不利于节能,而且造成碳氢、一氧化碳以及微粒等排放量增加。因此,根据起动工况的具体条件,设定合理的燃料喷射量且采用合理的喷射方式,不仅可以有效控制起动时间,而且还能有效控制起动过程的碳氢、一氧化碳以及微粒排放,以适应不断严格的车用发动机节能减排要求。
起动环境条件比较复杂,由于冷却液温度、喷射提前角、喷射压力、燃料喷射方式的不同,发动机指示效率差别较大,而且完成起动过程所需的燃料喷射量也不同。传统的起动控制方法,为了保证可靠起动,通常根据起动状态只简单地补偿一定的燃油喷射量,难以实现起动过程燃料喷射量的精确控制,因此直接影响发动机的起动性。随着汽车低碳化的发展,为完善车用发动机的起动性能,一种完善的起动控制策略是很有必要的。
发明内容
本发明为解决现有技术存在的难以实现发动机起动过程燃料喷射量的精确控制而影响发动机起动性的问题,提供了一种基于需求指示转矩和指示效率的起动燃料喷射量的控制方法。
为解决上述技术问题,本发明是采用如下技术方案实现的:
一种基于指示转矩的质调节式发动机起动控制方法,其特征在于包括如下步骤:
步骤一、通过实验获取发动机不同曲轴转速ncrank及冷却液温度Tw下的起动拖动转矩Trqd;根据发动机起动时的曲轴转速ncrank和冷却液温度Tw确定起动拖动转矩Trqd
步骤二、根据发动机的起动控制要求,确定起动过程的目标角加速度函数β0(t),基于起动过程的目标角加速度β0,计算加速阻力矩TrqI,实现对起动过程中发动机转速和起动时间的控制。
步骤三、根据发动机实际曲轴转速ncrank的变化计算发动机实际角加速度β,由发动机实际角加速度β与目标角加速度β0的偏差确定补偿转矩Trqc;实际角加速度β与目标角加速度β0的偏差为正,补偿转矩减小,偏差为负,则补偿转矩增加。
步骤四、将起动拖动转矩Trqd、加速阻力矩TrqI和补偿转矩Trqc相加得到起动需求的指示转矩Trqi,即Trqi=Trqd+TrqI+Trqc
步骤五、根据指示效率将需求指示转矩Trqi转换为燃料喷射量q,燃料喷射量q由基础燃料喷射量q0和燃料喷射修正量△q两部分组成:
进一步的技术方案如下:
步骤二的具体过程为:
a.根据起动时间的要求设定的预期起动时间t0,根据起动电机和发动机特性设定起动开始喷射转速n1、起动过程完成转速n2;获得目标发动机的转动惯量I。
b.根据发动机的起动控制要求,确定目标角加速度函数β0(t),β0(t)满足在预期起动时间t0内的积分为起动过程完成转速n2与起动开始喷射转速n1的差值,即 ∫ 0 t 0 β 0 ( t ) dt = n 2 - n 1 ;
c.利用公式TrqI=I·β0计算加速阻力矩TrqI,其中:β0为目标角加速度,I为发动机转动惯量。
步骤五的具体过程为:
a.根据标准常温起动条件下发动机的指示效率ηi0计算需求指示转矩Trqi所需的基础燃料喷射量q0其中π为圆周率,τ为冲程数,i为气缸数,Hu为燃料低热值;
b.起动条件改变时,根据指示效率的变化,修正燃料喷射量△q,即k为指示效率修正因子,通过实验获得。
c.将基础燃料喷射量q0和燃料喷射修正量△q相加获得起动燃料喷射量q,q=q0+△q。
与现有技术相比,本发明具有如下有益效果:起动过程严格遵循预期起动时间计算得出的角加速度进行,并附有补偿环节,保证各种环境条件下实现起动过程所需的最佳指示转矩,同时保证目标起动时间,起动过程可靠性高;通过设定不同的预期起动时间,从而改变发动机的实际起动时间;起动过程的燃料喷射量由需求指示转矩转换而得,在更改燃料喷射参数或燃烧条件变更时只需调整指示效率修正因子即可,保证各种起动条件下实现起动过程所需的最佳燃料喷射量;采用不同的目标角加速度的算法,可以调整起动过程曲线,实现起动向怠速及其他工况的平滑过渡。
附图说明
下面结合附图对本发明作进一步的说明:
图1为本发明专利的基本控制流程图
图2为本发明专利的实施例一的效果图
图3为本发明专利的实施例二的效果图
具体实施方式
为了详细地介绍本发明的内容,下面定义一些相关的概念:
定义一:指示转矩,指气缸一个工作循环热功转换所获得的指示功等效成的曲轴输出转矩。在起动阶段,所需求的指示转矩是用来克服发动机起动过程中阻力矩和惯性阻力矩,用Trqi表示。
定义二:起动拖动转矩,指发动机在某稳定冷却液温度时停油,利用电动机拖动发动机使其维持在某一起动转速所提供的转矩,用Trqd表示。
定义三:加速阻力矩(或惯性阻力矩),指起动过程中发动机转速升高时,用来克服发动机转动惯量所需要的转矩。加速阻力矩为发动机转动惯量与曲轴角加速度之积,用TrqI表示。
下面结合具体的实施例对本发明作详细的描述:
实施例一:
本发明所述的基于指示转矩的质调节式发动机起动控制包括以下步骤:
步骤1,通过实验获取发动机起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系,并结合实验获得的起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系,根据发动机当前曲轴转速ncrank和冷却液温度Tw确定起动拖动转矩Trqd。具体地,在本实施例中,在冷却液温度Tw下,停止燃料喷射,并将电力测功机转换为电动机倒拖发动机到给定的曲轴转速ncrank,并保持冷却液温度Tw不变,电力测功机所需提供的转矩即为起动拖动转矩Trqd;调整冷却液温度Tw和曲轴转速ncrank为其他值,得到对应的起动拖动转矩Trqd,从而建立起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系。
步骤2:根据发动机的起动控制要求,确定起动过程的目标角加速度函数β0(t),基于起动过程的目标角加速度β0,计算加速阻力矩TrqI。包括以下三个具体步骤:
a:根据起动时间的要求设定预期起动时间t0,根据起动电机和发动机特性设定起动开始喷射转速n1、起动过程完成转速n2;获得目标发动机的转动惯量I。起动开始转速n1是由起动电机带动发动机到达的开始燃料喷射的转速;起动结束转速n2是设定的发动机脱离起动状态的转速。本实施例中,t0=3s,n1=250r/min;n2=750r/min,I=0.5kg·m2
b:根据发动机的起动控制要求,确定目标角加速度函数β0(t),目标角加速度β0满足其在预期起动时间t0内的积分为起动过程完成转速n2与起动开始喷射转速n1的差值,即在本实施例中,n2=750r/min,n1=250r/min,I=0.5kg·m2,采用起动过程中目标角加速递减的方法,π为圆周率,t为当前相对于开始进入起动阶段的相对时刻。目标角加速度β0随着相对时刻t的增加而减少,当相对时刻t=t0时,目标角加速度β0为零。
c:基于目标角加速度β0和发动机转动惯量I,利用公式TrqI=I·β0计算加速阻力矩TrqI。本实施例中,TrqI=5.8(3-t)。
步骤3:根据发动机实际曲轴转速ncrank计算发动机实际角加速度β,根据发动机实际角加速度β与目标角加速度β0的偏差确定补偿转矩Trqc。实际角加速度β由曲轴转速ncrank微分计算后获得,即π为圆周率,t为当前相对于开始进入起动阶段的相对时刻。若实际角加速度β大于步骤4计算的目标角加速度β0,则补偿转矩Trqc减小,若实际际角加速度β小于步骤4计算的目标角加速度β0,则补偿转矩Trqc增加。本实施例中,补偿转矩I为发动机转动惯量,Fc为补偿系数,t为当前相对于开始进入起动阶段的相对时刻,t0为设定的预期起动时间,为3s,I=0.5kg·m2,Fc=1,所以有补偿系数Fc可根据实际发动机调整。
步骤4:获取起动需求的指示转矩Trqi。需求的指示转矩为起动拖动转矩Trqd、加速阻力矩TrqI、补偿转矩Trqc三者之和,即Trqi=Trqd+TrqI+Trqc
步骤5:根据指示效率将需求指示转矩Trqi转换为燃料喷射量q,燃料喷射量q由基础燃料喷射量q0和燃料喷射修正量△q两部分组成。包括以下具体三个步骤:
a:根据标准常温起动条件下(GB/T18297-2001《汽车发动机性能试验方法》)发动机的指示效率ηi0计算需求指示转矩Trqi所需的基础燃料喷射量q0其中π为圆周率,τ为冲程数,i为气缸数,Hu为燃料低热值。本实施例中,冲程数τ=4,气缸数i=4。
b:起动条件改变时,根据指示效率的变化计算燃料喷射修正量△q。其中kη为指示效率修正因子,通过实验获得。
c:将计算的基础燃料喷射量q0和燃料喷射修正量△q相加获得起动燃料喷射量q=q0+△q。
实施例二:
本发明所述的基于指示转矩的质调节式发动机起动控制包括以下步骤:
步骤1,通过实验获取发动机起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系,并结合实验获得的起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系,根据发动机当前曲轴转速ncrank和冷却液温度Tw确定起动拖动转矩Trqd。具体地,在本实施例中,在冷却液温度Tw下,停止燃料喷射,并将电力测功机转换为电动机倒拖发动机到给定的曲轴转速ncrank,并保持冷却液温度Tw不变,电力测功机所需提供的转矩即为起动拖动转矩Trqd;调整冷却液温度Tw和曲轴转速ncrank为其他值,得到对应的起动拖动转矩Trqd,从而建立起动拖动转矩Trqd与曲轴转速ncrank及冷却液温度Tw的关系。
步骤2:根据发动机的起动控制要求,确定起动过程的目标角加速度函数β0(t),基于起动过程的目标角加速度β0,计算加速阻力矩TrqI。包括以下三个具体步骤:
a:根据起动时间的要求设定预期起动时间t0,根据起动电机和发动机特性设定起动开始喷射转速n1、起动过程完成转速n2;获得目标发动机的转动惯量I。起动开始转速n1是由起动电机带动发动机到达的开始燃料喷射的转速;起动结束转速n2是设定的发动机脱离起动状态的转速。本实施例中,t0=3s,n1=250r/min;n2=750r/min,I=0.5kg·m2
b:根据发动机的起动控制要求,确定目标角加速度函数β0(t)。目标角加速度函数β0(t)满足其在预期起动时间t0内的积分为起动过程完成转速n2与起动开始喷射转速n1的差值,即在本实施例中,n2=750r/min,n1=250r/min,I=0.5kg·m2,t0=2s,采用起动过程中目标角加速度恒定的方法,即π为圆周率,t为当前相对于开始进入起动阶段的相对时刻。
c:利用公式TrqI=I·β0计算加速阻力矩TrqI,其中:β0为目标角加速度,I为发动机转动惯量。本实施例中:TrqI=13.09(kg·m)。
步骤3:根据发动机实际曲轴转速ncrank计算发动机实际角加速度β,根据发动机实际角加速度β与目标角加速度β0的偏差确定补偿转矩Trqc。实际角加速度β由曲轴转速ncrank微分计算后获得,即π为圆周率。本实施例中,补偿转矩I为发动机转动惯量,为0.5kg·m2,Fc为补偿系数,取1,t为当前相对于开始进入起动阶段的相对时刻,t0为设定的预期起动时间,为2s,所以有补偿系数Fc可根据实际发动机调整。
步骤4:获取起动需求的指示转矩Trqi。需求的指示转矩为起动拖动转矩Trqd、加速阻力矩TrqI、补偿转矩Trqc三者之和,即Trqi=Trqd+TrqI+Trqc
步骤5:根据指示效率将需求指示转矩Trqi转换为燃料喷射量q,燃料喷射量q由基础燃料喷射量q0和燃料喷射修正量△q两部分组成。包括以下具体三个步骤:
a:根据标准常温起动条件下(GB/T18297-2001《汽车发动机性能试验方法》)发动机的指示效率ηi0计算需求指示转矩Trqi所需的基础燃料喷射量q0其中π为圆周率,τ为冲程数,i为气缸数,Hu为燃料低热值。本实施例中,冲程数τ=4,气缸数i=4。
b:起动条件改变时,计算燃料喷射修正量△q。kη为指示效率修正因子,通过实验获得。
c:将计算的基础燃料喷射量q0和燃料喷射修正量△q相加获得起动燃料喷射量q=q0+△q。
实施例一和实施例二采用了不同的预期起动时间和不同的目标角加速度函数β0(t),起动过程的目标角加速规律不同,加速阻力矩的变化规律也不相同,起动过程的速度变化也不相同,起动效果分别如图2和图3所示。

Claims (3)

1.一种基于指示转矩的质调节式发动机起动控制方法,其特征在于包括如下步骤:
步骤一、通过实验获取发动机不同曲轴转速ncrank及冷却液温度Tw下的起动拖动转矩Trqd;根据发动机起动时的曲轴转速ncrank和冷却液温度Tw确定起动拖动转矩Trqd
步骤二、根据发动机的起动控制要求,确定起动过程的目标角加速度函数β0(t),基于起动过程的目标角加速度β0,计算加速阻力矩TrqI,实现对起动过程中发动机转速和起动时间的控制;
步骤三、根据发动机实际曲轴转速ncrank的变化计算发动机实际角加速度β,由发动机实际角加速度β与目标角加速度β0的偏差确定补偿转矩Trqc;实际角加速度β与目标角加速度β0的偏差为正,补偿转矩减小,偏差为负,则补偿转矩增加;
步骤四、将起动拖动转矩Trqd、加速阻力矩TrqI和补偿转矩Trqc相加得到起动需求的指示转矩Trqi,即Trqi=Trqd+TrqI+Trqc
步骤五、根据指示效率将需求指示转矩Trqi转换为燃料喷射量q,燃料喷射量q由基础燃料喷射量q0和燃料喷射修正量Δq两部分组成。
2.按照权利要求1所述的一种基于指示转矩的质调节式发动机起动控制方法,其特征在于步骤二的具体过程如下:
a.根据起动时间的要求设定的预期起动时间t0,根据起动电机和发动机特性设定起动开始喷射转速n1、起动过程完成转速n2;获得目标发动机的转动惯量I;
b.根据发动机的起动控制要求,确定目标角加速度函数β0(t),β0(t)满足在预期起动时间t0内的积分为起动过程完成转速n2与起动开始喷射转速n1的差值,即
c.利用公式TrqI=I·β0计算加速阻力矩TrqI,其中:β0为目标角加速度,I为发动机转动惯量。
3.按照权利要求1所述的一种基于指示转矩的质调节式发动机起动控制方法,其特征在于步骤五的具体过程如下:
a.根据标准常温起动条件下发动机的指示效率ηi0计算需求指示转矩Trqi所需的基础燃料喷射量q0其中π为圆周率,τ为冲程数,i为气缸数,Hu为燃料低热值;
b.起动条件改变时,根据指示效率的变化,修正燃料喷射量Δq,即k为指示效率修正因子,通过实验获得;
c.将基础燃料喷射量q0和燃料喷射修正量Δq相加获得起动燃料喷射量q,q=q0+Δq。
CN201410401131.0A 2014-08-14 2014-08-14 一种基于指示转矩的质调节式发动机起动控制方法 Expired - Fee Related CN104295387B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410401131.0A CN104295387B (zh) 2014-08-14 2014-08-14 一种基于指示转矩的质调节式发动机起动控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410401131.0A CN104295387B (zh) 2014-08-14 2014-08-14 一种基于指示转矩的质调节式发动机起动控制方法

Publications (2)

Publication Number Publication Date
CN104295387A CN104295387A (zh) 2015-01-21
CN104295387B true CN104295387B (zh) 2016-06-08

Family

ID=52315287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410401131.0A Expired - Fee Related CN104295387B (zh) 2014-08-14 2014-08-14 一种基于指示转矩的质调节式发动机起动控制方法

Country Status (1)

Country Link
CN (1) CN104295387B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6414125B2 (ja) * 2016-04-04 2018-10-31 トヨタ自動車株式会社 車両制御装置
CN106052917B (zh) * 2016-06-17 2018-09-07 广西师范大学 用角加速度参数优化发动机转矩测量误差的方法
CN106274514B (zh) * 2016-08-31 2018-09-14 北京新能源汽车股份有限公司 一种增程器控制方法、装置及增程式电动汽车
US10597020B2 (en) 2017-12-08 2020-03-24 GM Global Technology Operations LLC Powertrain with engine start function using resettable engine speed profile
CN108638915B (zh) * 2018-05-16 2021-04-09 江铃汽车股份有限公司 电动汽车行驶到蠕行速度前人为加油时的扭矩控制方法
CN111636970B (zh) * 2020-06-05 2021-08-24 吉林大学 一种基于复合喷射发动机的冷起动喷油点火控制方法
CN112282943B (zh) * 2020-10-30 2021-08-06 吉林大学 一种基于有效热效率的质调节式发动机的压缩比控制方法
CN113482784B (zh) * 2021-08-25 2022-08-23 潍柴动力股份有限公司 一种发动机起动扭矩确定方法及相关设备

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63248945A (ja) * 1987-04-06 1988-10-17 Toyota Motor Corp 内燃機関の燃料噴射制御装置
JP2004162623A (ja) * 2002-11-14 2004-06-10 Nissan Motor Co Ltd エンジンの燃料供給装置
JP4123244B2 (ja) * 2005-03-30 2008-07-23 トヨタ自動車株式会社 内燃機関の燃料噴射制御装置
JP2007327406A (ja) * 2006-06-07 2007-12-20 Toyota Motor Corp 内燃機関の制御装置及び方法
US7822536B2 (en) * 2006-06-13 2010-10-26 Toyota Jidosha Kabushiki Kaisha Start-up control device for internal combustion engine
JP2009121307A (ja) * 2007-11-14 2009-06-04 Nikki Co Ltd エンジンの燃料噴射制御方法及び燃料噴射制御装置
CN101492046A (zh) * 2008-11-25 2009-07-29 吉林大学 并联式混合动力车辆整车能量管理与总成协调控制方法
CN103414421A (zh) * 2013-07-23 2013-11-27 同济大学 基于Bang-bang控制的无刷直流电机换相转矩波动的抑制方法

Also Published As

Publication number Publication date
CN104295387A (zh) 2015-01-21

Similar Documents

Publication Publication Date Title
CN104295387B (zh) 一种基于指示转矩的质调节式发动机起动控制方法
CN100410514C (zh) 发动机的起动装置
CN102733962A (zh) 可变压缩比发动机的控制装置
CN102060014A (zh) 传动系刚度控制***和方法
CN104005861B (zh) 双燃料发动机的控制方法、***及车辆
CN104302895B (zh) 可变压缩比内燃机的控制装置
RU2699150C2 (ru) Способ регулирования расхода воздуха в двигателе (варианты)
CN103806998B (zh) 交通工具的动力系的能量控制***和方法
CN102400807A (zh) 发动机扭矩估计***和方法
CN102213150B (zh) 在低怠速时使用曲轴转速传感器使能气缸平衡的方法和***
WO2017138279A1 (ja) エンジンの制御装置
CN106246437B (zh) 基于发动机循环的发动机点火控制控制方法
CN108167084A (zh) 一种点燃与压燃方式结合的旋转活塞发动机电喷控制方法
CN104295388A (zh) 一种基于指示转矩的发动机各缸不均匀性补偿控制方法
EP2765291B1 (en) Method for operating a compression ignition internal combustion engine with direct fuel injection
CN204984644U (zh) 一种小型汽油机数字智能油门控制***
JP5273547B2 (ja) エンジン制御装置
CN101555840B (zh) 利用冷却剂和油温预测发动机增量摩擦转矩的***及方法
CN106870182A (zh) 一种发动机目标怠速转速控制方法
CN105626284A (zh) 双燃料发动机燃气流量控制方法及装置
CN109247022B (zh) 发动机的控制方法和发动机***
CN101871400B (zh) 用于协调式转矩控制的转矩储备和排放控制***
CN103498733A (zh) 一种自由活塞内燃发电机的运动控制方法
CN204532559U (zh) 降低柴油机瞬态工况烟度的控制***
CN203347928U (zh) 一种带点火功能的小型发动机电子控制单元

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160608

Termination date: 20190814

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