CN106126941A - 一种冷轧机组功率建模方法 - Google Patents

一种冷轧机组功率建模方法 Download PDF

Info

Publication number
CN106126941A
CN106126941A CN201610487049.3A CN201610487049A CN106126941A CN 106126941 A CN106126941 A CN 106126941A CN 201610487049 A CN201610487049 A CN 201610487049A CN 106126941 A CN106126941 A CN 106126941A
Authority
CN
China
Prior art keywords
power
rolling
strip steel
roll
computing formula
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
CN201610487049.3A
Other languages
English (en)
Other versions
CN106126941B (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.)
China University of Geosciences
Original Assignee
China University of Geosciences
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 China University of Geosciences filed Critical China University of Geosciences
Priority to CN201610487049.3A priority Critical patent/CN106126941B/zh
Publication of CN106126941A publication Critical patent/CN106126941A/zh
Application granted granted Critical
Publication of CN106126941B publication Critical patent/CN106126941B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Metal Rolling (AREA)

Abstract

一种冷轧机组功率建模方法。本发明针对酸洗‑连轧机组功率特点,提供了一种酸洗‑连轧产线机组功率建模的方法,对酸洗‑连轧产线中由变频电机控制的大功率机组建立了电功率模型,对酸洗槽建立了热功率模型,并将电功率模型和热功率模型耦合成为统一的功率模型。本发明所述基于电功率的机组以轧机机架、矫直机、卷曲机为例,其他机组的建模方法与之类似;热功率建模方法与电功率建模方法类似。本发明所建立的功率模型可以为后续冷轧工艺能量效率研究提供基础,也可以进行加工能耗评估和预测。建立的功率模型物理意义明确,合理、简单而准确,解决了酸洗‑连轧产线机组能效评估的问题。

Description

一种冷轧机组功率建模方法
技术领域
本发明提供了一种酸洗-连轧机组的功率建模方法,属于冷轧机组的能耗优化领域。
背景技术
2010年中国钢铁业工序能耗(吨钢可比综合能耗差7.47%,焦化工序差4.96%,炼铁差17.45%,轧钢差49.85%)中轧制工艺能耗落后于国际先进水平49.85%,是各种工序能耗差距最大的,节能降耗潜力很大。因此,提高冷轧工序机组能源消耗的效率,具有重要的意义。
建立酸洗-连轧产线机组输出功率的模型,不仅可以描述和评估酸洗-连轧产线机组的能源使用情况,也可以指导后续能量效率优化的研究。在当前的能效评估的研究中,“能量流(Energy Flow)”已经成为一种评估制造业能效的有效方法,并且在机床领域的研究已经达到比较成熟的层次,黄拯滔提出的一种基于能量守恒原理的数控车削能耗模型的精度已经达到99%以上(HUANG Zhengtao,ZHANG Chaoyong,LUO Min.An AssessmentModel of Energy Consumption for NC Turning Process Based on Principle ofConservation ofEnergy[J].Chinese Mechanical Engineer,2015,18(26):2419-2421)。冷轧工序作为多源能耗生产流程与传统机床能耗构成完全不同,将机床的能量流模型运用到冷轧机组,从电、热等多能量源进行综合能量流分析,对于冷轧工序的能耗控制和能效优化具有实际意义和广阔的应用前景。
目前钢铁冶金的功率模型的分析研究一般都集中在炼钢以及连铸等工序。赵业清(Zhao Yeqing.Research on Energy Flow Network Model in Iron and SteelEnterprises Based on Hybird Petri Net[J].Metallurgical IndustryAutomation,2014,38(5):27-28)应用能量流和物质流的方法对钢铁行业从原料到热轧的过程进行了分析。朱里红和王俊飞(Zhu Lihong,Wang Junfei,Zhang Feng.Development of EnergyMonitoring System Reheating Furnace in Hot Rolling Process[J].EngineeringControl Computer,2011,24(5):72-72)提出了一种在热轧过程的能耗监测加热炉实时能耗的***。李姗姗(Li Shanshan.Research on Energy Efficiency Assessment MethodofIron and Steel Production Process Based on Exergy Analysis[D].Jinan:Shandong University,2013)提出了一种基于分析的能耗分析方法,主要针对钢铁生产的炼钢工序。而关于冷轧(Cold Rolling)工序,尤其是生产流程复杂,涉及到的机组、参数众多的酸洗-连轧产线,机组的功率模型未见有文献报道。
冷轧作为钢铁生产的高能耗工序,节能潜力巨大。由于产线机组和工艺参数众多,生产过程复杂,关于冷轧能效优化的研究比较少,也缺少一个精确的冷轧机组功率模型,来分析冷轧机组功率与工艺参数之间关系。
发明内容
本发明提供了一种酸洗-连轧机组功率建模方法,用于评估和预测机组的能耗,解决了现有技术中的不足。
实现本发明上述目的所采用的技术方案为:
一种冷轧机组功率建模方法,包括以下步骤:(1)、对酸洗-连轧产线中由变频电机控制的大功率机组建立电功率模型,所述电功率模型包括轧机机架输入功率Pi、矫直机输入功率Pin、卷取机输入功率Pim;轧机机架输入功率Pi=Pz+Pf1+Pf2,其中Pz为轧制工艺所需功率,Pf1为辊颈与轴承摩擦损失的功率,Pf2为轧辊与带钢摩擦损失的功率;矫直机输入功率Pin=Pj+Pm,其中Pj为用于矫直工艺的功率,Pm为矫直辊与带钢摩擦损失的功率;卷取机输入功率Pim=Pq+Pqm,其中Pq为用于带钢卷取工艺的功率,Pqm为卷筒轴颈与轴承之间摩擦损失的功率;
(2)、对酸洗-连轧产线中的酸洗槽建立热功率模型,所述酸洗槽的输入功率为Qin=Q1+Q2+Q3+Q4;其中Q1为用于带钢加热功率,Q2为新酸加热功率,Q3为酸液蒸发散失功率,Q4为设备散热散失功率;
(3)、将电功率模型和热功率模型耦合成为统一的功率模型。
步骤(1)中Pz的计算公式如下:
P z = 2 p ‾ b ‾ ψ ( R Δ h + ( c p R ‾ ) 2 + c p ‾ R ) 2 × n 0 h 0 h i η ,
其中为轧制单位压力,为轧制前后带钢的平均宽度,ψ为轧制力臂系数,Δh为压下量,R为轧辊的半径,c为与工艺有关的系数,n0为机架轧辊转速,hi为带钢出i轧机时的厚度,η为电机到轧机的传动效率;
Pf1的计算公式为:
其中p为轧制力,d1为辊颈直径,f为轴颈和轴承之间的摩擦系数;
Pf2的计算公式为:
其中Mf2、Mf1为摩擦力矩,i为传动机构的传动比。
步骤(1)中Pj的计算公式为:
其中v为带钢运行的速度,σs为钢种的屈服强度,r0为带钢的初始曲率,kn为带钢经过第n个矫直辊时的残余曲率,E为带钢的弹性模量;
Pm的计算公式为:
其中μ和μ1分别为工作辊和支承辊轴承摩擦系数,d和d1分别为工作辊和支承辊的轴辊直径,f1为工作辊与支承辊间的滚动摩擦系数,为相邻上下两辊连线与垂直线夹角余弦值,t为相邻下辊圆心距离,Dg为下辊直径,Dz为上辊直径。
步骤(1)中Pq的计算公式为:
其中T为卷曲时的张力,P为卷曲支点作用力;
Pqm的计算公式为:Pqm=Ff2d2v5
其中,F为卷重与张力的合力,f2为轴颈的摩擦系数,d2为卷曲机轴颈直径,v5为卷曲机滚筒转速。
步骤(2)中Q1的计算公式为:
其中b为带钢的宽度,h为带钢厚度,v为带钢在酸洗工艺段运行速度,ρ1为带钢密度,C1为带钢的比热容,t1为带钢进酸洗槽时候的温度,t2为带钢的最终温度;
Q2的计算公式为:Q2=qC2(t2-t1),
其中q为新酸的流量,C2为新酸的比热容,t1为新酸温度,t2为酸洗槽内酸液温度;
Q3的计算公式为:Q3=kS1
其中K为蒸发散热系数,S1为蒸发液体表面积;
Q4的计算公式为:
其中λ2为设备导热系数,S2为设备表面积,Δt为设备内外壁温差,d2为设备厚度。
本发明中将电功率和电热功率耦合成为一个统一的功率模型后,表示如下:
P z = 2 p b ‾ ψ ( R Δ h + ( c p ‾ R ) 2 + c p ‾ R ) 2 η × n 0 h 0 h 1 P f = [ ( 1 η - 1 ) M Z + pd 1 f i + pd 1 f + ( 0.03 ~ 0.06 ) M H ] × n 0 h 0 h 1 . . . Q 1 = bhvρ 1 C 1 ( t 2 - t 1 ) 1000 Q 2 = qC 2 ( t s 2 - t s 1 ) Q 3 = kS 1 Q 4 = λ 2 S 2 Δ t d 2
与现有技术相比,本发明具有以下优点:
1.本发明主要针对冷轧工序的酸洗-连轧产线,通过分析机组输入输出功率的特性和机理,建立了电功率模型和热功率模型。不同于传统的功率模型,本发明将电功率模型和热功率模型耦合成为统一的功率模型。解决了当前酸洗-连轧产线缺少一种能效评估方法的问题,同时为后续的能效优化研究提供了理论依据。
2.本发明提出的方法基于能量流的理论方法,依据输入到机组的功率与机组输出的功率相等的原理,对每个机组的每个功率输出进行分解,并分别建立与工艺参数有关的平衡方程,获得的功率模型较为精确,如图1所示。
3.采用本发明方法建立的酸洗-连轧机组功率模型,只需给出合同中带钢的尺寸、弹性模量等参数,即可快速获得较为准确的轧制能耗,在实际生产中具有一定的预测功能。
附图说明
图1为本发明提供的酸洗-轧机组功率模型总示意图;
图2为迭代过程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实例,对本发明进行进一步详细说明。此处所描述的具体实例,仅用以解释本发明,并不用于限定本发明。
本发明所建立的功率模型可用于对不同品种的带钢加工时所需的功率进行预测,也能为后续的能效优化研究提供基础。以功率预测为例,机组的工艺能耗都与带钢的工艺参数存在一种映射的关系,所建立的功率模型实际上就是关于工艺参数的方程,方程的自变量为带钢自身的参数,因变量为输出功率的值。将产品的参数代入到模型方程,就可以实现加工所需的功率的预测和评估。模型自身应该具有一定的可靠性和准确性,表1中的数据来自某冷轧厂提供的数据,作为对照。
表1酸洗-连轧产线机组电机型号
机器名称 电机型号 容量(Kw)
七辊矫直机 1PQ8355-6PM80-Z 200
1号机架 TDZBS 4174-4 4174
2号机架 TDZBS 5250-4 5250
3号机架 TDZBS 5250-4 5250
4号机架 TDZBS 5250-4 5250
5号机架 TDZBS 5250-4 5250
卷取机 TDZBS 2850-4 2850
由于酸洗-连轧产线属于流程式生产,具有连续性的特点,一旦启动将不能随意停机,所以模型中的参数无法用实验的方法直接测的。本发明基于某冷轧厂的酸洗-连轧产线建立的功率模型,涉及到机组参数从该冷轧厂的提供的数据中获得;与工艺有关的参数将通过文献资料获取。
下面以单个的轧机机架为例对本发明的实施方法做详细具体的说明,矫直机输入功率、卷取机输入功率以及酸洗槽的输入功率的建模均可参照下述方法,本实施例中不进行一一累述。
本实施例中功率模型中出现的机组工艺参数如表2中所示:
表2模型中参数值
轧制工艺的功率方程为:
P z = 2 p ‾ b ‾ ψ ( R Δ h + ( c p R ‾ ) 2 + c p ‾ R ) 2 × n 0 h 0 h i η
表2中的参数代入后得到如下方程
P z = 1.05 p ‾ b ‾ ( 300 Δ h + ( 3.22 × 10 - 3 p ‾ ) 2 + 3.22 × 10 - 3 p ‾ ) 2 × n 0 h 0 h 1
上述方程中只剩是不确定的,需要通过公式(a)、(b)进行迭代求得,迭代过程如图2所示:
p ‾ = 95 e 0.05 l ′ h - 1 0.05 l ′ h ‾ - - - ( a )
l ′ = 300 Δ h + ( 3.22 × 10 - 3 p ‾ ) 2 + 3.22 × 10 - 3 p ‾ - - - ( b )
按照上述迭代法求出单位轧制压力带入轧制工艺功率模型中,得到一个以带钢自身工艺参数为自变量的方程。即可根据将要加工的带钢自身的工艺参数求出轧制工艺功率Pz
以1号机架为例进行计算,从某冷轧厂提供的合同明细中随机取5种带钢作为样本,将其尺寸参数代入到方程,得出的轧制不同钢种对应的轧机所需功率列在表3所示:
表3机组所需功率预测
表3中5种带钢样本的入口厚度、出口厚度以及宽度均来自某冷轧厂某年提供的一份合同明细。根据表中计算获得的数据,5种样本在1号机架中加工所需的轧制功率分1982.8Kw、2304.1Kw、2577.6Kw、3046.5Kw、4028.9Kw,表1中1号机架的电机容量为4174Kw,模型所预测的功率小于机架电机的容量,而且轧制功率随带钢厚度的增大而增大。由上可知本发明所建立的模型具有一定的可靠性,在实际工程应用中具有一定的价值。
本领域的技术人员容易理解,以上所述仅为本发明的实施例子而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种冷轧机组功率建模方法,其特征在于包括以下步骤:
(1)、对酸洗-连轧产线中由变频电机控制的大功率机组建立电功率模型,所述电功率模型包括轧机机架输入功率Pi、矫直机输入功率Pin、卷取机输入功率Pim;轧机机架输入功率Pi=Pz+Pf1+Pf2,其中Pz为轧制工艺所需功率,Pf1为辊颈与轴承摩擦损失的功率,Pf2为轧辊与带钢摩擦损失的功率;矫直机输入功率Pin=Pj+Pm,其中Pj为用于矫直工艺的功率,Pm为矫直辊与带钢摩擦损失的功率;卷取机输入功率Pim=Pq+Pqm,其中Pq为用于带钢卷取工艺的功率,Pqm为卷筒轴颈与轴承之间摩擦损失的功率;
(2)、对酸洗-连轧产线中的酸洗槽建立热功率模型,所述酸洗槽的输入功率为Qin=Q1+Q2+Q3+Q4;其中Q1为用于带钢加热功率,Q2为新酸加热功率,Q3为酸液蒸发散失功率,Q4为设备散热散失功率;
(3)、将电功率模型和热功率模型耦合成为统一的功率模型。
2.根据权利要求1所述的冷轧机组功率建模方法,其特征在于:步骤(1)中Pz的计算公式如下:
P z = 2 p ‾ b ‾ ψ ( R Δ h + ( c p R ‾ ) 2 + c p ‾ R ) 2 × n 0 h 0 h i η ,
其中为轧制单位压力,为轧制前后带钢的平均宽度,ψ为轧制力臂系数,Δh为压下量,R为轧辊的半径,c为与工艺有关的系数,n0为机架轧辊转速,hi为带钢出i轧机时的厚度,η为电机到轧机的传动效率;
Pf1的计算公式为:
其中p为轧制力,d1为辊颈直径,f为轴颈和轴承之间的摩擦系数;
Pf2的计算公式为:
其中Mf2、Mf1为摩擦力矩,i为传动机构的传动比。
3.根据权利要求1所述的冷轧机组功率建模方法,其特征在于:步骤(1)中Pj的计算公式为:
其中v为带钢运行的速度,σs为钢种的屈服强度,r0为带钢的初始曲率,kn为带钢经过第n个矫直辊时的残余曲率,E为带钢的弹性模量;
Pm的计算公式为:
其中μ和μ1分别为工作辊和支承辊轴承摩擦系数,d和d1分别为工作辊和支承辊的轴辊直径,f1为工作辊与支承辊间的滚动摩擦系数,为相邻上下两辊连线与垂直线夹角余弦值,t为相邻下辊圆心距离,Dg为下辊直径,Dz为上辊直径。
4.根据权利要求1所述的冷轧机组功率建模方法,其特征在于:步骤(1)中Pq的计算公式为:
其中T为卷曲时的张力,P为卷曲支点作用力;
Pqm的计算公式为:Pqm=Ff2d2v5
其中,F为卷重与张力的合力,f2为轴颈的摩擦系数,d2为卷曲机轴颈直径,v5为卷曲机滚筒转速。
5.根据权利要求1所述的冷轧机组功率建模方法,其特征在于:步骤(2)中Q1的计算公式为:
其中b为带钢的宽度,h为带钢厚度,v为带钢在酸洗工艺段运行速度,ρ1为带钢密度,C1为带钢的比热容,t1为带钢进酸洗槽时候的温度,t2为带钢的最终温度;
Q2的计算公式为:Q2=qC2(t2-t1),
其中q为新酸的流量,C2为新酸的比热容,t1为新酸温度,t2为酸洗槽内酸液温度;
Q3的计算公式为:Q3=kS1
其中K为蒸发散热系数,S1为蒸发液体表面积;
Q4的计算公式为:
其中λ2为设备导热系数,S2为设备表面积,Δt为设备内外壁温差,d2为设备厚度。
CN201610487049.3A 2016-06-28 2016-06-28 一种冷轧机组功率建模方法 Expired - Fee Related CN106126941B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610487049.3A CN106126941B (zh) 2016-06-28 2016-06-28 一种冷轧机组功率建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610487049.3A CN106126941B (zh) 2016-06-28 2016-06-28 一种冷轧机组功率建模方法

Publications (2)

Publication Number Publication Date
CN106126941A true CN106126941A (zh) 2016-11-16
CN106126941B CN106126941B (zh) 2018-06-26

Family

ID=57266837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610487049.3A Expired - Fee Related CN106126941B (zh) 2016-06-28 2016-06-28 一种冷轧机组功率建模方法

Country Status (1)

Country Link
CN (1) CN106126941B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107392491A (zh) * 2017-08-02 2017-11-24 中国地质大学(武汉) 一种评价冷轧带钢生产综合能量效率的方法
CN109829219A (zh) * 2019-01-22 2019-05-31 中国地质大学(武汉) 连轧机轧辊变速能量流模型和连轧机能量流模型建模方法
CN112462699A (zh) * 2020-11-02 2021-03-09 湖南华菱涟源钢铁有限公司 热轧生产电耗管理***及方法
CN113010980A (zh) * 2021-03-04 2021-06-22 西南石油大学 一种基于㶲分析的数控机床综合能效测试方法
CN113983721A (zh) * 2021-09-14 2022-01-28 山东澳信供热有限公司 空气源热泵维保***
CN114406006A (zh) * 2022-04-01 2022-04-29 承德建龙特殊钢有限公司 一种无缝钢管逐支跟踪生产***

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101934290B (zh) * 2009-06-30 2014-04-16 上海宝信软件股份有限公司 不锈钢冷连轧负荷分配调整方法
CN102266865B (zh) * 2011-05-30 2013-04-17 山东大学 热冷轧负荷分配方法
CN104998913B (zh) * 2015-06-29 2017-03-15 东北大学 一种冷轧轧制过程中电机功率的预测方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107392491A (zh) * 2017-08-02 2017-11-24 中国地质大学(武汉) 一种评价冷轧带钢生产综合能量效率的方法
CN109829219A (zh) * 2019-01-22 2019-05-31 中国地质大学(武汉) 连轧机轧辊变速能量流模型和连轧机能量流模型建模方法
CN112462699A (zh) * 2020-11-02 2021-03-09 湖南华菱涟源钢铁有限公司 热轧生产电耗管理***及方法
CN113010980A (zh) * 2021-03-04 2021-06-22 西南石油大学 一种基于㶲分析的数控机床综合能效测试方法
CN113983721A (zh) * 2021-09-14 2022-01-28 山东澳信供热有限公司 空气源热泵维保***
CN113983721B (zh) * 2021-09-14 2023-03-24 山东澳信供热有限公司 空气源热泵维保***
CN114406006A (zh) * 2022-04-01 2022-04-29 承德建龙特殊钢有限公司 一种无缝钢管逐支跟踪生产***
CN114406006B (zh) * 2022-04-01 2022-06-03 承德建龙特殊钢有限公司 一种无缝钢管逐支跟踪生产***

Also Published As

Publication number Publication date
CN106126941B (zh) 2018-06-26

Similar Documents

Publication Publication Date Title
CN106126941B (zh) 一种冷轧机组功率建模方法
CN103962390B (zh) Vc辊平整机湿平整过程中张力与轧制压力综合设定方法
CN102941232B (zh) 一种热连轧精轧过程控制方法
CN100550029C (zh) 冷连轧高速轧制过程中摩擦系数预报与设定方法
Wang et al. Research and development trend of shape control for cold rolling strip
CN103544340B (zh) 五机架冷连轧机组极薄带轧制中乳化液浓度的设定方法
CN105512804B (zh) 冷连轧过程以成本综合控制为目标的乳化液流量设定方法
CN113333474B (zh) 一种基于数字孪生的带钢热轧板形控制方法和***
CN107377634A (zh) 一种热轧带钢出口凸度预报方法
CN101507978B (zh) 冷连轧机基于机理与工况相结合的工作辊磨损预报方法
CN109241619A (zh) 3d模拟软件优化铝硅合金热轧工艺的方法
CN102513351A (zh) 一种带钢冷连轧轧制方法及装置
CN102266865A (zh) 热冷轧负荷分配方法
CN102847721A (zh) 一种热轧带钢轧辊热凸度确定方法
CN104942019A (zh) 一种带钢冷轧过程宽度自动控制方法
CN103028615A (zh) 一种预测带钢热连轧过程温度演变的方法
CN104998913A (zh) 一种冷轧轧制过程中电机功率的预测方法
CN109711101B (zh) 一种求解厚钢板同速异径蛇形轧制力能参数的方法
CN101733291A (zh) 中厚板轧后冷却辊道速度的控制方法
CN102172629A (zh) 一种用于cvc-6型平整机组的辊型控制方法
CN106345817A (zh) 工作辊长径比大且直径绝对值小的六辊轧机板形预报方法
CN102161054A (zh) 基于影响矩阵自学习的板形闭环控制方法
CN105234188A (zh) 冷连轧过程中以效益控制为目标的轧制速度优化方法
CN106825068B (zh) 一种轧制过程带钢表面粗糙度的预测方法
CN102688896A (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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180626

Termination date: 20190628

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