CN103171559B - 分模式最优化混联式混合动力汽车能量管理方法 - Google Patents

分模式最优化混联式混合动力汽车能量管理方法 Download PDF

Info

Publication number
CN103171559B
CN103171559B CN201310090705.2A CN201310090705A CN103171559B CN 103171559 B CN103171559 B CN 103171559B CN 201310090705 A CN201310090705 A CN 201310090705A CN 103171559 B CN103171559 B CN 103171559B
Authority
CN
China
Prior art keywords
mode
speed
motor
control algorithm
state
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
CN201310090705.2A
Other languages
English (en)
Other versions
CN103171559A (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.)
Tongji University
Original Assignee
Tongji 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 Tongji University filed Critical Tongji University
Priority to CN201310090705.2A priority Critical patent/CN103171559B/zh
Publication of CN103171559A publication Critical patent/CN103171559A/zh
Application granted granted Critical
Publication of CN103171559B publication Critical patent/CN103171559B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Electric Propulsion And Braking For Vehicles (AREA)
  • Hybrid Electric Vehicles (AREA)

Abstract

本发明涉及一种分模式最优化混联式混合动力汽车能量管理方法,具体包括以下步骤:1)状态机根据车速、电池充电状态值和车辆需求功率确定当前汽车所处的工作模式;2)模式选择控制器获取状态机确定的工作模式,选取与当前工作模式相对应的控制算法及所需传感器信号;3)由选定的控制算法得到各驱动源转矩信号和制动转矩信号,并将驱动源转矩信号传送给发动机控制器和电机控制器,制动转矩信号传送给制动控制器。与现有技术相比,本发明具有实用性高、计算负荷低、优化性能好、适用范围广等优点。

Description

分模式最优化混联式混合动力汽车能量管理方法
技术领域
本发明涉及一种混合动力汽车的能量管理方法,尤其是涉及一种分模式最优化混联式混合动力汽车能量管理方法。
背景技术
混合动力汽车是当前清洁能源汽车的发展主流,它可弥补纯电动汽车续驶里程短的不足,同时相对传统汽车又可大幅度地提高能源利用率,减少尾气排放。混合动力汽车有内燃机、电机等多个动力源,结构复杂,如何实时有效地对多个动力源进行能量管理仍是混合动力汽车控制领域的难题之一。
混联式混合动力汽车的传动系统主要有两种结构,即行星齿轮结构和双转子电机结构,其中行星齿轮结构被广泛采用。混联式混合动力汽车的能量管理策略主要有三种方法,即基于规则方法、最优化方法和智能化方法,而实际被采用的只有基于规则的能量管理策略。基于规则的能量管理策略主要以发动机为优化对象,由发动机的万有特性曲线确定发动机的最优工作区域,根据最优工作区域来制定控制策略,目前,该方法主要采用发动机恒定工作点控制策略或发动机最优工作曲线控制策略。基于规则的方法虽然控制简单,实用性强,但无法保证整车在运行过程中始终处于最优状态。基于最优化方法的能量管理策略是以车辆所有动力部件为优化对象,根据车辆当前的运行状态以及将来的运行信息,对车辆各个动力部件进行实时优化控制,保证整车在运行过程中始终处于最优状态。目前,正在研究的主要方法有动态规划法、等效消耗最小法和预测控制法。基于最优化的方法理论上可以保证车辆驱动系统实时处于最优运行状态,但由于控制变量多、计算负荷大,且需获取将来的行驶路况信息,因此现在还无法实际应用到普通车辆上。智能化能量管理策略大都以车辆所有动力部件为优化对象,以人的思维方式对车辆驱动系统在运行过程中进行优化控制,如模糊控制和神经网络等方法,智能化方法虽然鲁棒性强,但优化效果无法得到保障。
在专利《一种混合动力车辆能量管理方法及能量管理系统》(申请公布号为CN102717797A)中,公布了一种基于远程服务器进行最优化控制的能量管理策略,把采集到的相关数据通过GPRS实时得传递给远程服务器,建立驾驶员需求功率转移概率矩阵,并基于随机动态规划建立能量管理状态转移矩阵,通过迭代计算得到最优解,再传输给车辆控制器进行控制。这种方法计算负荷大,且仅适用于路线固定的公交汽车。在专利《一种混联式混合动力车辆的动力控制方法》(申请公布号CN102815295A)中,公布了一种基于需求功率进行分模式控制的能量管理策略,但未能保证每个模式下车辆驱动系统都处于最优状态。在文献″Afour-stepmethodtodesignanenergymanagementstrategyforhybridvehicles″(ZhuYuan,ChenYaobin,TianGuangyu,eta1.ProceedingsofAmericanControlConference,pp.156-161,2004)中,针对一种混联式混合动力电动汽车,提出过分模式控制方法,但每个模式内均采用模糊控制方法,优化方法只用来调整相关参数,所以控制效果的优化性不能保证。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种实用性高、计算负荷低、优化性能好、适用范围广的分模式最优化混联式混合动力汽车能量管理方法。
本发明的目的可以通过以下技术方案来实现:
一种分模式最优化混联式混合动力汽车能量管理方法,所述的混合动力汽车中设有行星齿轮耦合装置以及与行星齿轮耦合装置连接的发动机、第一电机和第二电机,发动机有发动机控制器,两电机有各自的控制器,所述的控制方法具体包括以下步骤:
1)状态机根据车速、电池充电状态值和车辆需求功率对汽车的工作模式进行切换,以确定当前汽车所处的工作模式;
2)模式选择控制器获取状态机确定的工作模式,选取与当前工作模式相对应的控制算法及所需传感器信号;
3)由选定的控制算法得到各驱动源转矩信号和制动转矩信号,并将驱动源转矩信号传送给发动机控制器和电机控制器,制动转矩信号传送给制动控制器。
所述的状态机对汽车的工作模式进行切换具体为:
101)状态机判断是否满足当前车速小于下限值、电池SOC值大于下限值且车辆需求功率小于电池额定功率,若满足,则将工作模式切换为低速模式,并执行步骤103),否则执行步骤102);
102)判断当前车速是否小于上限值,若是,则将工作模式切换为常规模式,并执行步骤103),若否,则将工作模式切换为高速模式;
103)判断车辆需求功率是否小于等于零,若是,则进入再生制动模式,并执行步骤104),若否,则执行步骤101);
104)判断车辆需求功率是否大于零,若是,则执行步骤101),若否,则保持再生制动模式。
当汽车的工作模式为低速模式时,所述的控制算法为PID。控制器采用PID控制算法根据电机需求转速与电机实际转速的差值输出低速模式下电机需求转矩;
当汽车的工作模式为常规模式、高速模式或再生制动模式时,所述的控制算法为模型预测控制算法,模式选择控制器获取当前工作模式,由当前工作模式下的状态方程和目标函数采用模型预测控制算法获得最优控制量。
当汽车的工作模式为常规模式时,所述的模型预测控制算法的具体控制过程为:
201)确定常规模式下的状态方程及目标函数,常规模式下的状态方程为
x · = Ax + B u u + B v v y = Cx
其中状态量x=(weng,wMG2,SOC)T,weng为发动机转速,wMG2为第二电机转速,SOC为电池充电状态值;控制量u=(Teng,TMG2,TMG1T,Teng为发动机转矩,TMG2为第二电机转矩,TMG1第一电机转矩;v为扰动量;y为输出量,A、Bu、Bv和C为关系矩阵,u为输入量;
常规模式下的目标函数是
J=ayfy(y,u,Δu)+aflffl(u,Δu)+aesfes(u,△u)+aSOCfSOC(u,△u)
其中fy(y,u,Δu)为驱动性能函数,fyl(u,Δu)为燃油性函数,fes(u,Δu)为排放性函数,fSOC(u,Δu)为电池充电状态函数,ay、afi、aes和aSOC为各函数的权值;
202)根据当前时刻的各状态量和控制量,结合驾驶需求功率,在满足各状态量及控制量的约束条件下,使目标函数最小,求得下一时刻的最优控制量值。
与现有技术相比,本发明具有以下优点:
1、综合了基于规则和基于最优化方法的能量管理策略的优点;
2、实用性高,与已有最优化法能量管理策略相比,由于采用分模式最优化方法,控制变量个数减少,状态空间方程得以精简,计算负荷降低,实用性提高;
3、优化性能好,与已实用的能量管理策略相比,由于各模式下都采用最优化控制策略,车辆驱动系统始终处于最优状态;
4、适用范围广,由于无需知晓将来路况信息,适用于各种路线行驶的车辆;
5、模式划分条件与模式切换条件简单,操作性好,由于模式划分主要依据车速,而电池充电状态和车辆需求功率只做补充的切换条件。
附图说明
图1为本发明的工作模式切换流程图;
图2为本发明结构示意图;
图3为动力耦合装置连接示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
一种分模式最优化混联式混合动力汽车能量管理方法,所述的混合动力汽车中设有行星齿轮耦合装置以及与行星齿轮耦合装置连接的发动机和电机,如图3所示,所述的行星齿轮耦合装置包括相连的发动机、第一电机和第二电机,发动机与齿架相连,即A2端口;第一电机与太阳轮相连,即A1端口;第二电机与齿圈相连,即A3端口;A4端口与驱动轴连接,发动机与发动机控制器连接,第一电机和第二电机均有各自的电机控制器。混合动力汽车中还设有相连接的状态机、模式选择控制器及控制信号选择切换模块,控制信号选择切换模块与发动机控制器和电机控制器相连接。
如图1-图2所示,所述的控制方法具体包括以下步骤:
1)状态机根据车速、电池充电状态值和车辆需求功率对汽车的工作模式进行切换,以确定当前汽车所处的工作模式;
2)模式选择控制器获取状态机确定的工作模式,选取与当前工作模式相对应的控制算法及所需传感器信号;
3)由选定的控制算法得到各驱动源转矩信号和制动转矩信号,并将驱动源转矩信号传送给发动机控制器和电机控制器,制动转矩信号传送给制动控制器。
如图1所示,汽车起动后,结合图1对控制流程进行详细描述:
在步骤401中,判断是否满足当前车速小于常规模式车速的下限值、电池SOC值大于下限值且车辆需求功率小于电池额定功率,若满足,则进入步骤402,否则进入步骤403。
在步骤402中,运行低速模式,由第二电机单独驱动,动力电池提供功率。同步进入步骤406。
在步骤403中,判断当前车速是否小于常规模式车速的上限值,若是,则进入步骤404,否则进入步骤405。
在步骤404中,运行常规模式,发动机和电机混合驱动,无滑行制动和刹车制动,采用模型预测控制算法,通过传感器信息和驾驶员信息实时更新状态方程参数,得到最优控制信息,保证车辆动力性、高效性,维持电池充电状态。同步进入步骤406。
在步骤405中,运行高速模式,除状态方程不同外,其余同常规模式。
在步骤406中,判断驾驶需求功率是否小于等于零,若是,则进入步骤407,否则进入步骤401。
在步骤407中,运行再生制动模式,发动机空转或停止,电机回馈制动,同刹车制动协调控制,得到最大回收能量。同步进入步骤408。
在步骤408中,判断驾驶需求功率是否大于零,若是,则进入步骤401,否则进入步骤407。
如图2所示,模式选择控制器中设有低速模式控制单元、常规模式控制单元、高速模式控制单元和再生制动模式控制单元,各模式控制单元采用各自的控制算法输出控制信号。
当汽车的工作模式为低速模式时,所述的控制算法为PID,低速模式控制单元采用PID控制算法根据电机需求转速与电机实际转速的差值输出低速模式下电机需求转矩;
当汽车的工作模式为常规模式、高速模式或再生制动模式时,所述的控制算法为模型预测控制算法,各模式控制单元根据当前工作模式下的状态方程和目标函数,采用模型预测控制算法获得最优控制量。
采用模型预测控制算法时,目标函数考虑驱动性、油耗、排放及电池充电状态,各部件协调控制,保证车辆驱动系统处于最优状态。以常规模式为例,所述的模型预测控制算法的具体控制过程为:
确定常规模式下的状态方程及目标函数,常规模式下的状态方程为
x · = Ax + B u u + B v v y = Cx
其中状态量x=(weng,wMG2,SOC)T,weng为发动机转速,wMG2为第二电机转速,SOC为电池充电状态值;控制量u=(Teng,TMG2,TMG1)T,Teng为发动机转矩,TMG2为第二电机转矩,TMG1第一电机转矩;v为扰动量;y为输出量,A、Bu、Bv和C为关系矩阵,u为输入量;
常规模式下的目标函数是
J=ayfy(y,u,△u)+aflffl(u,Δu)+aesfes(u,Δu)+asocfsoc(u,Δu)
其中fy(y,u,Δu)为驱动性能函数,ffl(u,Δu)为燃油性函数,fes(u,Δu)为排放性函数,fSOC(u,Δu)为电池充电状态函数,ay、afl、aes和aSOC为各函数的权值;
根据当前时刻的各状态量和控制量,结合驾驶需求功率,在满足各状态量及控制量的约束条件下,使目标函数最小,求得下一时刻的最优控制量值。

Claims (3)

1.一种分模式最优化混联式混合动力汽车能量管理方法,所述的混合动力汽车中设有行星齿轮耦合装置以及与行星齿轮耦合装置连接的发动机、第一电机和第二电机,发动机连接有发动机控制器,两电机有各自的控制器,其特征在于,所述的能量管理方法具体包括以下步骤:
1)状态机根据车速、电池充电状态值和车辆需求功率对汽车的工作模式进行切换,以确定当前汽车所处的工作模式;
所述的状态机对汽车的工作模式进行切换具体为:
101)状态机判断是否满足当前车速小于下限值、电池SOC值大于下限值且车辆需求功率小于电池额定功率,若满足,则将工作模式切换为低速模式,并执行步骤103),否则执行步骤102);
102)判断当前车速是否小于上限值,若是,则将工作模式切换为常规模式,并执行步骤103),若否,则将工作模式切换为高速模式;
103)判断车辆需求功率是否小于等于零,若是,则进入再生制动模式,并执行步骤104),若否,则执行步骤101);
104)判断车辆需求功率是否大于零,若是,则执行步骤101),若否,则保持再生制动模式;
2)模式选择控制器获取状态机确定的工作模式,选取与当前工作模式相对应的控制算法及所需传感器信号;
3)由选定的控制算法得到各驱动源转矩信号和制动转矩信号,并将驱动源转矩信号传送给发动机控制器和电机控制器,制动转矩信号传送给制动控制器。
2.根据权利要求1所述的一种分模式最优化混联式混合动力汽车能量管理方法,其特征在于,当汽车的工作模式为低速模式时,所述的控制算法为PID,控制器采用PID控制算法根据电机需求转速与电机实际转速的差值输出低速模式下电机需求转矩;
当汽车的工作模式为常规模式、高速模式或再生制动模式时,所述的控制算法为模型预测控制算法,模式选择控制器获取当前工作模式,由当前工作模式下的状态方程和目标函数采用模型预测控制算法获得最优控制量。
3.根据权利要求2所述的一种分模式最优化混联式混合动力汽车能量管理方法,其特征在于,当汽车的工作模式为常规模式时,所述的模型预测控制算法的具体控制过程为:
201)确定常规模式下的状态方程及目标函数,常规模式下的状态方程为
x · = A x + B u u + B v v y = C x
其中状态量x=(weng,wMG2,SOC)T,weng为发动机转速,wMG2为第二电机转速,SOC为电池充电状态值;控制量u=(Teng,TMG2,TMG1)T,Teng为发动机转矩,TMG2为第二电机转矩,TMG1第一电机转矩;v为扰动量;y为输出量,A、Bu、Bv和C为关系矩阵,u为输入量;
常规模式下的目标函数是
J=ayfy(y,u,Δu)+aflffl(u,Δu)+aesfes(u,Δu)+aSOCfSOC(u,Δu)
其中fy(y,u,Δu)为驱动性能函数,ffl(u,Δu)为燃油性函数,fes(u,Δu)为排放性函数,fSOC(u,Δu)为电池充电状态函数,ay、afl、aes和aSOC为各函数的权值;
202)根据当前时刻的各状态量和控制量,结合驾驶需求功率,在满足各状态量及控制量的约束条件下,使目标函数最小,求得下一时刻的最优控制量值。
CN201310090705.2A 2013-03-20 2013-03-20 分模式最优化混联式混合动力汽车能量管理方法 Expired - Fee Related CN103171559B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310090705.2A CN103171559B (zh) 2013-03-20 2013-03-20 分模式最优化混联式混合动力汽车能量管理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310090705.2A CN103171559B (zh) 2013-03-20 2013-03-20 分模式最优化混联式混合动力汽车能量管理方法

Publications (2)

Publication Number Publication Date
CN103171559A CN103171559A (zh) 2013-06-26
CN103171559B true CN103171559B (zh) 2016-04-20

Family

ID=48631881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310090705.2A Expired - Fee Related CN103171559B (zh) 2013-03-20 2013-03-20 分模式最优化混联式混合动力汽车能量管理方法

Country Status (1)

Country Link
CN (1) CN103171559B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6060839B2 (ja) * 2013-07-17 2017-01-18 トヨタ自動車株式会社 ハイブリッド車両の制御装置
DE102015215546A1 (de) * 2015-08-14 2017-02-16 Robert Bosch Gmbh Verfahren und Vorrichtung zum Feststellen, ob in einem Kraftfahrzeug ein Fehlerzustand vorliegt oder nicht
CN106080585B (zh) * 2016-06-22 2020-04-03 江苏大学 一种双行星排式混合动力汽车非线性模型预测控制方法
JP6725880B2 (ja) * 2016-09-29 2020-07-22 三菱自動車工業株式会社 ハイブリッド車両の制御装置
CN106950862A (zh) * 2017-04-09 2017-07-14 怀宁鑫橙信息技术有限公司 一种快递车动力模式智能远程控制方法
CN108394403B (zh) * 2018-04-28 2019-07-05 吉林大学 一种功率分流式混合动力汽车模式切换控制方法
US11117567B2 (en) * 2018-06-26 2021-09-14 Toyota Motor Engineering & Manufacturing North America, Inc. Real time trajectory optimization for hybrid energy management utilizing connected information technologies
CN109495050B (zh) * 2018-11-27 2020-09-18 浙江大学 基于二次型价值函数的双电机转矩同步模型预测控制方法
CN109849894B (zh) * 2019-03-15 2020-09-01 重庆长安新能源汽车科技有限公司 一种用于四驱混合动力汽车的能量管理方法、系统及设备
CN110001620B (zh) * 2019-04-24 2021-04-13 吉林大学 一种液压轮毂混合动力车辆多模式切换控制方法
CN111746260A (zh) * 2020-06-29 2020-10-09 一汽解放汽车有限公司 一种48v系统及48v系统控制方法
JP7280902B2 (ja) * 2021-01-28 2023-05-24 本田技研工業株式会社 車両制御装置
CN113173153B (zh) * 2021-06-10 2024-03-08 北京理工大学 一种无轨胶轮车的控制方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101492046A (zh) * 2008-11-25 2009-07-29 吉林大学 并联式混合动力车辆整车能量管理与总成协调控制方法
CN102069798A (zh) * 2010-12-29 2011-05-25 浙江钱江摩托股份有限公司 一种并联式轻度混合动力摩托车的控制装置及其控制方法
CN202499132U (zh) * 2012-03-05 2012-10-24 浙江大学城市学院 新型Plug_in混合动力汽车能量管理控制器
CN102815295A (zh) * 2012-08-23 2012-12-12 浙江吉利汽车研究院有限公司杭州分公司 一种混联式混合动力车辆的动力控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7076350B2 (en) * 2003-12-19 2006-07-11 Lear Corporation Vehicle energy management system using prognostics
US8560153B2 (en) * 2010-02-15 2013-10-15 Ford Global Technologies, Llc Vehicle battery state of charge hold function and energy management

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101492046A (zh) * 2008-11-25 2009-07-29 吉林大学 并联式混合动力车辆整车能量管理与总成协调控制方法
CN102069798A (zh) * 2010-12-29 2011-05-25 浙江钱江摩托股份有限公司 一种并联式轻度混合动力摩托车的控制装置及其控制方法
CN202499132U (zh) * 2012-03-05 2012-10-24 浙江大学城市学院 新型Plug_in混合动力汽车能量管理控制器
CN102815295A (zh) * 2012-08-23 2012-12-12 浙江吉利汽车研究院有限公司杭州分公司 一种混联式混合动力车辆的动力控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
混联型混合动力汽车控制策略优化研究;刘文杰;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20071115(第05期);第19-24、35-41页以及图3.1、4.1-4.4 *

Also Published As

Publication number Publication date
CN103171559A (zh) 2013-06-26

Similar Documents

Publication Publication Date Title
CN103171559B (zh) 分模式最优化混联式混合动力汽车能量管理方法
CN108382186B (zh) 混联式混合动力系统及车辆工作模式决策方法
CN108394403B (zh) 一种功率分流式混合动力汽车模式切换控制方法
CN109606348B (zh) 一种插电式行星混联汽车能量管理控制方法
Banvait et al. A rule-based energy management strategy for plug-in hybrid electric vehicle (PHEV)
KR100949260B1 (ko) 전기자동차용 전지 충전 시스템
CN101557974B (zh) 混合动力车辆及混合动力车辆的控制方法
CN102729987B (zh) 一种混合动力公交车能量管理方法
CN103770778B (zh) 一种基于最低使用成本考虑的增程式电车的控制方法
JP4292131B2 (ja) ハイブリッド自動車のエンジン運転停止を要求する方法及びシステム
CN103587522B (zh) 一种混合动力汽车动力总成智能控制方法
CN104787036B (zh) 基于预测的电池能量使用控制电池功率的系统和方法
CN101254744B (zh) 一种混合动力多能源控制系统及其控制方法
CN106080585B (zh) 一种双行星排式混合动力汽车非线性模型预测控制方法
CN106080579B (zh) 一种基于悬架振动能量回收的混合动力汽车整车控制方法
CN103112450B (zh) 一种插电式并联混合动力汽车的实时优化控制方法
CN107458369B (zh) 一种同轴并联式混合动力电动汽车能量管理方法
CN101492046A (zh) 并联式混合动力车辆整车能量管理与总成协调控制方法
CN101708722A (zh) 基于模糊逻辑的串联式混合动力电动汽车整车控制方法
CN104760591B (zh) 混合动力综合控制系统
CN108327512A (zh) 混合动力驱动系统及车辆
CN105599755A (zh) 一种插电式混合动力客车的驱动轴转矩控制方法
CN105015543B (zh) 混合动力汽车的扭矩分配方法
CN104260720A (zh) 一种混联混合动力控制系统及采用该系统实现的控制方法
CN105128855A (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: 20160420

Termination date: 20190320

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