CN101957426A - 集成电路综合环境试验方法 - Google Patents

集成电路综合环境试验方法 Download PDF

Info

Publication number
CN101957426A
CN101957426A CN2009100882400A CN200910088240A CN101957426A CN 101957426 A CN101957426 A CN 101957426A CN 2009100882400 A CN2009100882400 A CN 2009100882400A CN 200910088240 A CN200910088240 A CN 200910088240A CN 101957426 A CN101957426 A CN 101957426A
Authority
CN
China
Prior art keywords
test
integrated circuit
temperature
environment
minimum
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
CN2009100882400A
Other languages
English (en)
Other versions
CN101957426B (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.)
BEIJING SAN-TALKING TESTING ENGINEERING ACADEMY Co Ltd
Original Assignee
BEIJING SAN-TALKING TESTING ENGINEERING ACADEMY 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 BEIJING SAN-TALKING TESTING ENGINEERING ACADEMY Co Ltd filed Critical BEIJING SAN-TALKING TESTING ENGINEERING ACADEMY Co Ltd
Priority to CN2009100882400A priority Critical patent/CN101957426B/zh
Publication of CN101957426A publication Critical patent/CN101957426A/zh
Application granted granted Critical
Publication of CN101957426B publication Critical patent/CN101957426B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种集成电路综合环境试验方法,包括:S1,根据集成电路的产品基本信息和使用的环境信息,获取试验最高温度TTH、试验最低温度TTL、试验最高温度变化率ΔTE、试验功率密度谱、试验湿度和试验电应力;S2,根据上述试验数据,分别在给所述集成电路施加额定工作电压、最高工作电压、最低工作电压时,同时进行温度循环试验、随机振动试验和湿度试验;S3,上述试验后,对所述集成电路进行功能测试,如果所述集成电路功能丧失,或任一性能参数值超出设计范围,则所述集成电路失效;否则,所述集成电路完好。本发明能真实模拟集成电路的实际使用环境,更快更全面的激发产品的故障缺陷,为改进集成电路设计、制造缺陷提供参考。

Description

集成电路综合环境试验方法
技术领域
本发明涉及集成电路技术领域,尤其涉及一种集成电路综合环境试验方法。
背景技术
集成电路在实际使用中会经受到各种环境应力的作用,产品的各薄弱环节在这种综合环境应力的作用下极容易暴露出来。近年来,大规模集成电路的研制单位和用户都在不同场合提出,迫切需求能快速有效的评价大规模集成电路的质量与可靠性的试验方法。资料表明,在众多环境因素中对产品质量与可靠性影响最大的是温度(40%),振动(28%)和湿度(18%,湿度是指相对湿度,即空气中实际所含水蒸汽密度和同温度下饱和水蒸汽密度的百分比值。)。开展集成电路综合环境极限应力试验就是研究温度-湿度-振动的综合极限应力试验的方法和程序,使用该方法和程序可以更快更全面的激发产品的故障缺陷,为研制方提供一种快速、准确评价大规模集成电路质量与可靠性的手段,试验结果可为改进集成电路设计、工艺上的缺陷提供参考资料,为集成电路用户验收产品提供试验方法。现有技术中还没有形成集成电路综合极限应力试验的试验方法。集成电路的设计研制周期短,产品的质量与可靠性要求高,现有的成熟的试验方法不能满足快速评价集成电路的质量与可靠性以及失效率的需求。
发明内容
本发明的目的是提供一种集成电路综合环境试验方法,以解决现有技术中不能真实模拟集成电路实际使用环境的缺陷。
为了达到上述目的,本发明的技术方案提出一种集成电路综合环境试验方法,所述方法包括以下步骤:
S1,根据集成电路的产品基本信息和使用的环境信息,获取试验数据,所述试验数据包括:试验最高温度TTH、试验最低温度TTL、试验最高温度变化率ΔTE、试验功率密度谱、试验湿度和试验电应力;
S2,根据上述试验数据,分别在给所述集成电路施加额定工作电压、最高工作电压、最低工作电压时,同时进行温度循环试验、随机振动试验和湿度试验;
S3,上述试验后,对所述集成电路进行功能测试,如果所述集成电路功能丧失,或任一性能参数值超出设计范围,则所述集成电路失效;否则,所述集成电路完好。
其中,所述集成电路的产品基本信息包括:最高工作温度TWH、最低工作温度TWL、额定工作电压V、最高工作电压VH和最低工作电压VL
使用的环境信息包括:最高环境温度TEH、最低环境温度TEL、最高功率谱密度W、最高温度变化率ΔTE
其中,所述试验最高温度TTH取TWH与TEH+Tq的最小值,试验最低温度TTL取TWL与TEL-Tq的最大值,其中,Tq为加严温度;试验最高温度变化率ΔTE为集成电路使用环境信息的最高温度变化率ΔTE;试验功率密度谱的值为集成电路使用环境信息的最高功率谱密度W。
其中,同时进行温度循环试验、随机振动试验和湿度试验的试验步骤为:
S2-1,试验温度为TTL,所述集成电路的开关应至少通断两次;
S2-2,温度由TTL上升到TTH,所述集成电路的开关处于接通状态;
S2-3,温度维持在TTH,所述集成电路的开关至少通断两次,并且控制湿度为80%-90%;
S2-4,温度由TTH降至TTL,所述集成电路的开关处于接通状态;
随机振动试验在上述过程中的开关处于稳定状态时施加。
其中,所述步骤S2-1中,试验最低温度TTL持续时间为40min。
其中,所述步骤S2-2中,温度由TTL上升到TTH的持续时间为t=(TTH-TTL)/ΔTE
其中,所述步骤S2-3中,温度TTH持续时间为40min。
其中,所述步骤S2-4中,温度由TTH下降到TTL的持续时间为t′=(TTH-TTL)/ΔTE
本发明的技术方案结合温度、湿度和振动三个方面的影响对集成电路进行综合环境试验,更能真实的模拟集成电路的实际使用环境,可以更快更全面的激发产品的故障缺陷,为研制方提供一种快速、准确评价大规模集成电路质量与可靠性的方法,并且试验结果可对改进集成电路设计、工艺上的缺陷提供参考。
附图说明
图1是本发明实施例的一种集成电路综合环境试验方法的流程示意图;
图2是本发明实施例的一种集成电路综合环境试验方法的随机振动图谱;
图3是本发明实施例的一种集成电路综合环境试验方法的时间对照图。
具体实施方式
以下实施例用于说明本发明,但不用来限制本发明的范围。
图1为本发明集成电路综合环境试验方法的实施例流程图,本实施例的试验方法包括以下步骤:
S101、集成电路提供商给出产品的基本信息,用户提供器件(集成电路)使用环境的信息,做为制定综合试验条件的输入。
集成电路提供商应提供的产品基本信息,包括:器件最高工作温度、器件最低工作温度、器件额定工作电压、器件最高工作电压、器件最低工作电压等。
用户提供器件使用的环境信息,包括:最高环境温度、最低环境温度、最高振动量值、最高温度变化率等。
表1集成电路提供商应提供的产品基本信息
Figure B2009100882400D0000041
表2用户提供器件使用的环境信息
S102、根据器件的基本信息和使用环境信息分项确定试验条件;
首先根据输入的信息分项确定试验条件。需要确定的试验条件包括:试验最高温度、试验最低温度、试验温度变化率、试验湿度、试验振动量、试验电应力。
试验最高温度TTH根据以下公式确认:TTH=min(TWH,TEH+Tq);即取TWH与TEH+Tq的最小值。
试验最低温度TTL根据以下公式确认:TTL=max(TWL,TEL-Tq);即取TWL与TEL-Tq的最大值。
其中,Tq为加严温度,可由用户与研制单位协商,一般推荐值为11℃。
试验最高温度变化率ΔTE为器件使用的最高温度变化率ΔTE;试验振动条件按图2所示的随机振动图谱进行试验,试验功率密度谱的值为器件使用的最高功率谱密度W。
试验湿度:在试验最高温度时注入85%~90%湿度,其余时间不加控制。
试验电应力:每次试验循环过程依次以额定工作电压V、最高工作电压VH、最低工作电压VL的顺序进行循环。
S103、制定综合环境试验方案;
综合环境试验需用温度应力、振动应力、湿度应力、电应力同时变化模拟综合环境,根据分项试验条件来确定综合环境的试验条件,从而制定综合环境的试验方案。
整个试验应至少重复十个子循环,一个子循环由一个完整的电应力(VL、VH、V)循环构成,每个电应力(VL、VH、V)包括A、B、C、D四个阶段(参考图3和表3)。
下面对一个子循环的试验步骤进行详细描述:
1)电应力为额定电压V:A阶段:施加温度应力TTL,此阶段开关应至少通断两次。B阶段:温度由TTL上升到TTH,开关处于接通状态。C阶段:温度维持在TTH,开关至少通断两次,加入80%-90%的湿度应力。D阶段:温度由TTH降至TTL,开关处于接通状态。振动应力在整个过程的开关处于稳定状态时施加。
2)电应力为最高工作电压VH,重复步骤1)中的A、B、C、D四个阶段的步骤。
3)电应力为最低工作电压VL。重复步骤1)中的A、B、C、D四个阶段的步骤。
对图3的详细描述如下:
1.电应力的施加按以下要求变化:第一循环的输入电压为标称电压V;第二循环的输入电压为VH;第三循环的输入电压为VL。三个试验循环的输入电压变化构成一个完整的电应力循环,整个试验期间重复这一电应力循环。
2.温度应力的施加,在每一循环的开始即A段施加试验最低温度TTL,持续时间为40min;然后由TTL上升到TTH为B段,持续时间为t=(TTH-TTL)/ΔTE;C段施加试验最高温度TTH,持续时间为40min;最后进入D段,试验温度由TTH下降到TTL,此阶段持续时间为t′=(TTH-TTL)/ΔTE
3.湿度应力的施加,在每一循环的C段即施加试验最高温度TTH时加入湿度应力,其余阶段不加控制。
4.振动应力的施加,振动应力应在开关处于稳定状态时施加,振动类型为随机振动,振动方向为X、Y、Z各向,振动量值为W。
5.开关状态,开关在冷浸和热浸阶段需至少启动2次,在整个试验循环的其它阶段,开关均应处于接通状态。
表3综合环境试验剖面时序表
Figure B2009100882400D0000061
S104、试验结果判据分析、试验夹具的要求。
1.试验结果判据
器件在断电后突然加电能否正常工作,是考核开关在极端温度下的启动能力。若在此条件下无法实现应有功能,即视为失效。
在综合试验结束后进行的电性能测试中,若器件在加电状态下不能正常工作,要求的任一性能参数值超出器件的设计范围,即视为失效。
若是由于试验夹具和试验箱在试验过程中失效导致的故障不视为器件失效。
2.振动试验夹具设计要求
试验夹具是用于将试验件安装在振动台台面上的过渡段,它的一端与振动台运动部分机械连接;而另一端与试验件机械连接,从而将振动台的机械能传递到试验件上。为了将振动台的机械振动能保真地传递到试验件上,夹具应满足下列要求:
为了保证试验的准确性,夹具应满足下列要求:
a.在试验频带内夹具的传递特性无共振峰;
b.在试验频带内不产生与主振动方向正交的横向运动;
c.为防止试验中出现温度不均匀、升降温缓慢等现象,夹具应选取热容小的材料制作,且要考虑受热面对温度影响等因素。
d.夹具与试验件接口面运动均匀。
本发明可用于集成电路但不仅限于集成电路。其它类型电子元器件如分立器件、阻容元件、电感、继电器、连接器、光电器件、电子模块、电路板等产品也可使用本方法。以上为本发明的最佳实施方式,依据本发明公开的内容,本领域的普通技术人员能够显而易见地想到一些雷同、替代方案,均应落入本发明保护的范围。

Claims (8)

1.一种集成电路综合环境试验方法,其特征在于,所述方法包括以下步骤:
S1,根据集成电路的产品基本信息和使用的环境信息,获取试验数据,所述试验数据包括:试验最高温度TTH、试验最低温度TTL、试验最高温度变化率ΔTE、试验功率密度谱、试验湿度和试验电应力;
S2,根据上述试验数据,分别在给所述集成电路施加额定工作电压、最高工作电压、最低工作电压时,同时进行温度循环试验、随机振动试验和湿度试验;
S3,上述试验后,对所述集成电路进行功能测试,如果所述集成电路功能丧失,或任一性能参数值超出设计范围,则所述集成电路失效;否则,所述集成电路完好。
2.如权利要求1所述的集成电路综合环境试验方法,其特征在于,所述集成电路的产品基本信息包括:最高工作温度TWH、最低工作温度TWL、额定工作电压V、最高工作电压VH和最低工作电压VL
使用的环境信息包括:最高环境温度TEH、最低环境温度TEL、最高功率谱密度W、最高温度变化率ΔTE
3.如权利要求2所述的集成电路综合环境试验方法,其特征在于,所述试验最高温度TTH取TWH与TEH+Tq的最小值,试验最低温度TTL取TWL与TEL-Tq的最大值,其中,Tq为加严温度;试验最高温度变化率ΔTE为集成电路使用环境信息的最高温度变化率ΔTE;试验功率密度谱的值为集成电路使用环境信息的最高功率谱密度W。
4.如权利要求3所述的集成电路综合环境试验方法,其特征在于,同时进行温度循环试验、随机振动试验和湿度试验的试验步骤为:
S2-1,试验温度为TTL,所述集成电路的开关应至少通断两次;
S2-2,温度由TTL上升到TTH,所述集成电路的开关处于接通状态;
S2-3,温度维持在TTH,所述集成电路的开关至少通断两次,并且控制湿度为80%-90%;
S2-4,温度由TTH降至TTL,所述集成电路的开关处于接通状态;
随机振动试验在上述过程中的开关处于稳定状态时施加。
5.如权利要求4所述的集成电路综合环境试验方法,其特征在于,所述步骤S2-1中,试验最低温度TTL持续时间为40min。
6.如权利要求4所述的集成电路综合环境试验方法,其特征在于,所述步骤S2-2中,温度由TTL上升到TTH的持续时间为t=(TTH-TTL)/ΔTE
7.如权利要求4所述的集成电路综合环境试验方法,其特征在于,所述步骤S2-3中,温度TTH持续时间为40min。
8.如权利要求4所述的集成电路综合环境试验方法,其特征在于,所述步骤S2-4中,温度由TTH下降到TTL的持续时间为t′=(TTH-TTL)/ΔTE
CN2009100882400A 2009-07-14 2009-07-14 集成电路综合环境试验方法 Expired - Fee Related CN101957426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100882400A CN101957426B (zh) 2009-07-14 2009-07-14 集成电路综合环境试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100882400A CN101957426B (zh) 2009-07-14 2009-07-14 集成电路综合环境试验方法

Publications (2)

Publication Number Publication Date
CN101957426A true CN101957426A (zh) 2011-01-26
CN101957426B CN101957426B (zh) 2012-07-18

Family

ID=43484870

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100882400A Expired - Fee Related CN101957426B (zh) 2009-07-14 2009-07-14 集成电路综合环境试验方法

Country Status (1)

Country Link
CN (1) CN101957426B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102628901A (zh) * 2012-03-30 2012-08-08 株洲南车时代电气股份有限公司 一种轨道交通电气设备可靠性试验方法
CN103852672A (zh) * 2014-03-31 2014-06-11 中国航天科工运载技术研究院北京分院 一种以步进应力对电子设备实施环境应力筛选的方法
CN104155200A (zh) * 2014-07-23 2014-11-19 西安空间无线电技术研究所 一种快速温度变化试验抗结霜凝露的方法
CN104181456A (zh) * 2013-05-27 2014-12-03 海洋王(东莞)照明科技有限公司 印刷电路板测试方法
CN104596719A (zh) * 2014-12-25 2015-05-06 北京自动测试技术研究所 一种航天用集成电路极限应力强度的快速评价方法
CN104932328A (zh) * 2015-04-30 2015-09-23 中国电子科技集团公司第四十一研究所 一种综合环境试验处理方法
CN106546834A (zh) * 2015-09-17 2017-03-29 大唐移动通信设备有限公司 一种环境应力筛选试验方法及装置
CN106777730A (zh) * 2016-12-26 2017-05-31 中国航空工业集团公司西安飞机设计研究所 一种可靠性综合应力条件的制定方法
CN106841847A (zh) * 2016-12-21 2017-06-13 吴中区穹窿山倪源交通器材经营部 一种轨道交通电气设备可靠性试验方法
CN106918772A (zh) * 2017-02-27 2017-07-04 中国人民解放军军械工程学院 雷达电路板电应力加速退化方法和装置
CN107064665A (zh) * 2016-12-01 2017-08-18 中国北方车辆研究所 电子产品综合应力强化试验方法
CN107560788A (zh) * 2016-07-01 2018-01-09 无锡华润上华科技有限公司 压力传感器芯片的测试方法
CN108427126A (zh) * 2018-01-31 2018-08-21 交通运输部水运科学研究所 基带芯片高低温特性的测试方法和装置
CN110389296A (zh) * 2019-06-21 2019-10-29 深圳米飞泰克科技有限公司 集成电路成品处理方法、装置及终端设备
CN112268781A (zh) * 2020-10-15 2021-01-26 中国电子科技集团公司第四十九研究所 一种高分子湿敏电容极限环境应力的确定方法
CN112684324A (zh) * 2020-12-30 2021-04-20 无锡市同步电子科技有限公司 一种机载电子控制器用pcb故障快速激发与验证的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5471877A (en) * 1993-07-16 1995-12-05 At&T Ipm Corp. Environmental stress screening process with liquid coupled vibration
US5744975A (en) * 1996-06-06 1998-04-28 International Business Machines Corporation Enhanced defect elimination process for electronic assemblies via application of sequentially combined multiple stress processes
US6446508B1 (en) * 2001-01-17 2002-09-10 Venturedyne, Ltd. Vibration compartment environmental control
TWM310340U (en) * 2006-11-13 2007-04-21 Princeton Technology Corp Test system for testing chips
CN101339302B (zh) * 2007-07-05 2010-06-16 比亚迪股份有限公司 一种ito测试板和测试方法

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102628901A (zh) * 2012-03-30 2012-08-08 株洲南车时代电气股份有限公司 一种轨道交通电气设备可靠性试验方法
CN102628901B (zh) * 2012-03-30 2014-11-12 株洲南车时代电气股份有限公司 一种轨道交通电气设备可靠性试验方法
CN104181456A (zh) * 2013-05-27 2014-12-03 海洋王(东莞)照明科技有限公司 印刷电路板测试方法
CN103852672A (zh) * 2014-03-31 2014-06-11 中国航天科工运载技术研究院北京分院 一种以步进应力对电子设备实施环境应力筛选的方法
CN104155200A (zh) * 2014-07-23 2014-11-19 西安空间无线电技术研究所 一种快速温度变化试验抗结霜凝露的方法
CN104596719A (zh) * 2014-12-25 2015-05-06 北京自动测试技术研究所 一种航天用集成电路极限应力强度的快速评价方法
CN104932328A (zh) * 2015-04-30 2015-09-23 中国电子科技集团公司第四十一研究所 一种综合环境试验处理方法
CN104932328B (zh) * 2015-04-30 2017-11-17 中国电子科技集团公司第四十一研究所 一种综合环境试验处理方法
CN106546834A (zh) * 2015-09-17 2017-03-29 大唐移动通信设备有限公司 一种环境应力筛选试验方法及装置
CN106546834B (zh) * 2015-09-17 2019-08-09 大唐移动通信设备有限公司 一种环境应力筛选试验方法及装置
CN107560788A (zh) * 2016-07-01 2018-01-09 无锡华润上华科技有限公司 压力传感器芯片的测试方法
CN107064665B (zh) * 2016-12-01 2019-08-09 中国北方车辆研究所 电子产品综合应力强化试验方法
CN107064665A (zh) * 2016-12-01 2017-08-18 中国北方车辆研究所 电子产品综合应力强化试验方法
CN106841847A (zh) * 2016-12-21 2017-06-13 吴中区穹窿山倪源交通器材经营部 一种轨道交通电气设备可靠性试验方法
CN106777730A (zh) * 2016-12-26 2017-05-31 中国航空工业集团公司西安飞机设计研究所 一种可靠性综合应力条件的制定方法
CN106918772A (zh) * 2017-02-27 2017-07-04 中国人民解放军军械工程学院 雷达电路板电应力加速退化方法和装置
CN106918772B (zh) * 2017-02-27 2019-08-06 中国人民解放军军械工程学院 雷达电路板电应力加速退化方法和装置
CN108427126A (zh) * 2018-01-31 2018-08-21 交通运输部水运科学研究所 基带芯片高低温特性的测试方法和装置
CN110389296A (zh) * 2019-06-21 2019-10-29 深圳米飞泰克科技有限公司 集成电路成品处理方法、装置及终端设备
CN112268781A (zh) * 2020-10-15 2021-01-26 中国电子科技集团公司第四十九研究所 一种高分子湿敏电容极限环境应力的确定方法
CN112684324A (zh) * 2020-12-30 2021-04-20 无锡市同步电子科技有限公司 一种机载电子控制器用pcb故障快速激发与验证的方法

Also Published As

Publication number Publication date
CN101957426B (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN101957426B (zh) 集成电路综合环境试验方法
CN103778292B (zh) 一种热振联合载荷下bga焊点疲劳寿命预测方法
Xia et al. Fatigue life prediction of package-on-package stacking assembly under random vibration loading
CN109614684B (zh) 一种多应力综合下的卫星电子产品寿命预测方法
CN102663201B (zh) 一种考虑电气互连结构可靠性的电子产品可靠性预计方法
Catelani et al. Critical components test and reliability issues for photovoltaic inverter
Ye et al. Reliability evaluation of CSP soldered joints based on FEM and Taguchi method
Al-Yafawi et al. Random vibration test for electronic assemblies fatigue life estimation
Wyrwas et al. Accurate quantitative physics-of-failure approach to integrated circuit reliability
JP2005026250A (ja) はんだ接合部の亀裂発生寿命予測方法及び装置
CN116973726A (zh) Ic老化测试方法、装置、设备及存储介质
Qin et al. Comparing and modifying estimation methods of fatigue life for PCBA under random vibration loading by finite element analysis
Xu et al. Numerical studies of the mechanical response of solder joint to drop/impact load
Wang et al. Prediction of chip solder fatigue in IGBTs
CN112287582B (zh) 一种陶瓷焊柱阵列的热疲劳寿命优化方法
Núñez et al. Experimental Characterization of Time-Dependent Variability in Ring Oscillators
CN202770931U (zh) 具有样品接线端子的安全型环境试验箱
Bosco et al. Simulation and experiment of thermal fatigue in the CPV die attach
Wang et al. Thermal stress analysis of BGA solder joint power load based on comsol
CN117849592B (zh) 一种复杂电磁环境下电力芯片的加速老化测试方法
Krishna et al. Bivariate zero-inflated power series distribution
Wang et al. SMD reliability research base on Stolkarts fatigue model
Li et al. Profile simplification of reliability prediction based on failure physics
Latha et al. Critical Component Identification in Analog Electronic Circuits using Sensitivity Analysis
Kang et al. Design analysis and optimization of metal core bgas for warpage

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: 20120718

Termination date: 20180714

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