CN108152617B - 测量热电模块最大输出功率的方法、装置及终端设备 - Google Patents
测量热电模块最大输出功率的方法、装置及终端设备 Download PDFInfo
- Publication number
- CN108152617B CN108152617B CN201711344727.1A CN201711344727A CN108152617B CN 108152617 B CN108152617 B CN 108152617B CN 201711344727 A CN201711344727 A CN 201711344727A CN 108152617 B CN108152617 B CN 108152617B
- Authority
- CN
- China
- Prior art keywords
- thermoelectric module
- thermoelectric
- output power
- module
- determining
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明适用于热电转换技术领域,提供了一种测量热电模块最大输出功率的方法、装置及终端设备,根据热电模块输出电流的变化范围,以及热电模块的热端输入热能与热电模块的输出电流的变化关系,确定热端输入热能的变化范围,进而确定热电材料平均温度的变化范围和动态内阻的变化范围,在动态内阻的变化范围内,逐次调节测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率,确定热电模块的最大输出功率,缩小了测试电阻的调节范围,测量过程简单方便,节省了大量时间。
Description
技术领域
本发明属于热电转换技术领域,尤其涉及一种测量热电模块最大输出功率的方法、装置及终端设备。
背景技术
半导体温差发电主要利用半导体热电材料赛贝克效应将热能转换为电能。热电模块(Thermo electric Module,TEM)作为热能到电能转换的核心部件,最大输出功率是表征其热电转换性能的主要参数之一。目前,获得热电模块冷热端一定温差下最大输出功率的常用方法有三种:1.开路电压-短路电流法;2.伏-安特性法;3.负载匹配法。
但是,开路电压-短路电流法和伏-安特性法都是以理论值计算热电模块的最大输出功率,均未考虑热电模块中存在热电耦合效应,因此,按经典理论计算得到的热电模块最大输出功率准确度较低。负载匹配法为热电模块最大输出功率的经典测量法,基于热电模块冷热两端在一定温差下,调节负载阻值,找到负载电阻与热电模块内阻实际匹配工况,根据匹配工况下测量得到的热电模块输出端的电压(Vo)、电流(Io)计算得到热电模块最大输出功率(Pomax=VoIo),由于该测量方法完全模拟了热电模块实际工况,故得到的热电模块最大输出功率值比较准确,但是该测量方法需一个一个负载值改变,并且无法确定负载值调节范围是否覆盖指定温差下热点模块最大输出功率对应的内阻值,测量时均需要冷热端温差、电压/电流输出进入稳定状态,因此,测量过程复杂、耗时过长是目前负载匹配法测量热电模块最大输出功率存在的不足。
发明内容
有鉴于此,本发明实施例提供了一种测量热电模块最大输出功率的方法、装置及终端设备,以解决现有技术中采用负载匹配法测量热电模块最大输出功率,测量过程复杂、耗时过长的问题。
本发明实施例第一方面提供了一种测量热电模块最大输出功率的方法,包括以下步骤:
获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与热电模块的输出电流的变化关系,确定所述热端输入热能的变化范围。
根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围。
根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围。
在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
进一步地,所述测量热电模块最大输出功率的方法还包括:
根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻。
调节热电模块的外接测试电阻的初始阻值等于所述平均动态内阻。
进一步地,所述热电模块的热端输入热能与热电模块的输出电流的变化关系为:
其中,参数A=Rh-(1-η)Rc,参数B=1-0.5η+K[Rh+(1-η)Rc],Tc、Th分别为热电模块冷、热端温度,Qh为热电模块的热端输入热能,η为热电模块的热电转换效率,N为热电模块中的P-N电偶对数,α为热电模块的塞贝克系数,K为热电模块的导热系数,I为热电模块的输出电流。
所述热端输入热能与所述热电材料平均温度的变化关系为:
其中,参数A=Rh-(1-η)Rc,Tj为热电材料的平均温度,Tc、Th分别为热电模块冷、热端温度,Tcj、Thj分别为热电材料的冷、热端结面温度,Qh为热电模块的热端输入热能,Rc、Rh分别为热电模块冷、热端寄生电阻,η为热电模块的热电转换效率。
所述热电材料电导率与所述热电材料平均温度的变化关系为:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度。
所述热电材料电导率与所述动态内阻的变化关系为:
其中,N为热电模块中的P-N电偶对数,hn、hp分别为热电模块的N型、P型热电材料的臂高,ln、lp分别为热电模块的N型、P型热电材料横截面边长,σn、σp分别为热电模块的N型、P型热电材料在平均温度Tj下的电导率,Ri为热电模块在平均温度Tj下的动态内阻。
进一步地,所述在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率包括:
在所述动态内阻的变化范围内,确定所述测试电阻的初始调节方向。
在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
进一步地,所述根据测量的热电模块的输出功率确定热电模块的最大输出功率包括:
根据测量的热电模块的输出功率,确定热电模块的输出功率的峰值。
根据热电模块的输出功率的峰值,确定热电模块的最大输出功率。
本发明实施例第二方面提供了一种测量热电模块最大输出功率的装置,包括:
第一变化范围确定模块,用于获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与所述输出电流的变化关系,确定所述热端输入热能的变化范围;
第二变化范围确定模块,用于根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围;
第三变化范围确定模块,用于根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围;
输出功率测量模块,用于在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
最大输出功率确定模块,用于根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
进一步地,所述测量热电模块最大输出功率的装置还包括:
平均动态内阻确定模块,用于根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻;
初始阻值调节模块,用于调节所述测试电阻的初始阻值等于所述平均动态内阻。
进一步地,所述输出功率测量模块包括:
初始调节方向确定模块,用于在所述动态内阻的变化范围内,确定所述测试电阻的初始调节方向;
第一输出功率测量模块,用于在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
第二输出功率测量模块,用于在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
本发明实施例第三方面提供了一种测量热电模块最大输出功率的终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述实施例第一方面提供的测量热电模块最大输出功率的方法。
本发明实施例的第四方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述实施例第一方面提供的测量热电模块最大输出功率的方法。
本发明实施例与现有技术相比存在的有益效果是:本发明提供的一种测量热电模块最大输出功率的方法、装置及终端设备,根据热电模块输出电流的变化范围,以及热电模块的热端输入热能与热电模块的输出电流的变化关系,确定热端输入热能的变化范围,进而确定热电材料平均温度的变化范围和动态内阻的变化范围,在动态内阻的变化范围内,逐次调节测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率,确定热电模块的最大输出功率,缩小了测试电阻的调节范围,测量过程简单方便,节省了大量时间。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种测量热电模块最大输出功率的方法的流程示意图;
图2是本发明实施例提供的热电材料与温度的变化曲线图;
图3是本发明实施例提供的一种测量热电模块最大输出功率的装置的示意性框图;
图4是本发明实施例提供的一种测量热电模块最大输出功率的终端设备的结构示意图。
具体实施方式
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
为了说明本发明的技术方案,下面通过具体实施例进行说明。
参阅图1,本发明实施例提供了一种测量热电模块最大输出功率的方法,包括以下步骤:
S101,获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与所述输出电流的变化关系,确定所述热端输入热能的变化范围。
具体地,在热电模块输出开路的工况下,开路电流为Ioc,输出电流最小I=Ioc=0,在热电模块输出短路的工况下,短路电流为Isc,输出电流最大I=Isc,因此,输出电流的变化范围是0~Isc。
设热电模块冷、热端温度分别为Tc、Th;热电材料的冷、热结面温度分别为Tcj、Thj,热电模块冷、热端热能分别为Qc、Qh;热电模块的温差为ΔTj,热电材料平均温度为Tj,热电模块冷、热端寄生电阻分别为Rh、Rc,热电模块的输出功率为Po,热电模块的热电转换效率为η,热电模块的导热系数为K,热电模块的塞贝克系数为α,热电模块的动态内阻为Ri,热电模块的输出电流为I,测试电阻为RL,则:
Thj=Th-QhRh,Tcj=Tc+QcRc
Thj+Tcj=(Th+Tc)-QhRh+QcRc=(Th+Tc)-QhRh+(Qh-Po)Rc
=(Th+Tc)-QhRh+Qh(1-η)Rc=(Th+Tc)-[Rh-(1-η)Rc]Qh
ΔTj=Thj-Tcj=(Th-Tc)-QhRh-QcRc=(Th-Tc)-QhRh-(1-η)QhRc
=(Th-Tc)-[Rh+(1-η)Rc]Qh
令参数A=Rh-(1-η)Rc,Thj+Tcj=(Th+Tc)-AQh
Qh=NαIThj+K(Thj-Tcj)-0.5I2Ri=NαIThj+KΔTj-0.5I(NαΔTj-IRL)
=NαIThj-0.5NαIΔTj+KΔTj+0.5I2RL=0.5NαI(Thj+Tcj)+KΔTj+0.5I2RL
=0.5NαI(Thj+Tcj)+KΔTj+0.5Po
=0.5NαI(Thj+Tcj)+KΔTj+0.5ηQh
=0.5NαI(Th+Tc)-0.5NαIAQh+K{(Th-Tc)-[Rh+(1-η)Rc]Qh}+0.5ηQh
{1-0.5η+K[Rh+(1-η)Rc]+0.5NAαI}Qh=0.5Nα(Th+Tc)I+K(Th-Tc)
令参数B=1-0.5η+K[Rh+(1-η)Rc]
上式中,A=Rh-(1-η)Rc>0,并且
B(Th+Tc)-K(Th-Tc)A={1-0.5η+K[Rh+(1-η)Rc]}(Th+Tc)-K(Th-Tc)[Rh-(1-η)Rc]
=(1-0.5η)(Th+Tc)+KRhTh+KRhTc+K(1-η)RcTh+K(1-η)RcTc-KRhTh+KRhTc
+K(1-η)RcTh-K(1-η)RcTc
=(1-0.5η)(Th+Tc)+2K[RhTc+(1-η)RcTh]>0
因此,在热电模块的冷端温度Tc和热端温度Th一定时,热电模块所需的热端输入热能Qh与热电模块的输出电流I有关,并且热端输入热能Qh随输出电流I的增大而增大。
在热电模块输出开路的工况下,开路电流为Ioc,输出电流I=Ioc=0,此时输出电流为最小值,热端输入热能也取最小值,即Qh=Qhmin。
在热电模块输出短路的工况下,短路电流为Isc,输出电流I=Isc,此时输出电流为最大值,热端输入热能也取最大值,即Qh=Qhmax。
在热电模块输出功率最大的工况下,输出电流I在0~Isc之间,约为0.5Isc,此时输入热能Qh满足Qhmin<Qh<Qhmax,这样,根据热电模块的热端输入热能与热电模块的输出电流的变化关系,就确定热端输入热能的变化范围。
S102,根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围。
具体地,热电材料的平均温度满足:
即热电材料平均温度随所述热端输入热能的增大而降低。
在步骤S101中得出输入热能Qh的范围是Qhmin<Qh<Qhmax,结合上述热端输入热能与所述热电材料平均温度的变化关系可知:
在热电模块输出开路的工况下,开路电流为Ioc,输出电流I=Ioc=0,此时输出电流为最小值,热端输入热能也取最小值,即Qh=Qhmin,而热电材料平均温度Tj为最大值,设为Tjmax。
在热电模块输出短路的工况下,短路电流为Isc,输出电流I=Isc,此时输出电流为最大值,热端输入热能也取最大值,即Qh=Qhmax,而热电材料平均温度Tj为最小值,设为Tjmin。
因此,可以确定热电材料平均温度的变化范围是Tjmin<Tj<Tjmax。
S103,根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围。
具体地,热电材料的电导率与温度的变化特性如图2所示,可以看出无论半导体热电材料用于发电还是制冷,其电导率与材料的平均温度呈反向变化,即热电材料的电导率随热电材料平均温度的降低而增大。
以另一具体的热电材料为例,其电导率与温度的变化关系满足:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度。
热电材料电导率与热电模块动态内阻之间满足:其中,N为热电模块中的P-N电偶对数,hn、hp分别为热电模块的N型、P型热电材料的臂高,ln、lp分别为热电模块的N型、P型热电材料横截面边长,σn、σp分别为热电模块的N型、P型热电材料在平均温度Tj下的电导率,Ri为热电模块在平均温度Tj下的动态内阻。
通过上述分析,不难得出热电模块动态内阻随热电材料电导率的增大而减小,进而可知动态内阻随热电材料平均温度的降低而减小。
因此,在热电模块输出开路的工况下,热电材料平均温度Tj为最大值Tjmax,此时对应的动态内阻为最大值,设为Rioc。
在热电模块输出短路的工况下,热电材料平均温度Tj为最小值Tjmin,此时对应的动态内阻为最小值,设为Risc。
对于热电模块稳态开路工况,假设稳态开路电压为Voc1、瞬间短路电流为Isc1,由于热电模块温度场变化滞后于热电模块电场变化,当数据采样时间≤30ms,读取热电模块瞬间短路电流Isc1时,认为热电模块温度场近似没有改变,因此,
在热电模块的冷、热端温度分别为Tc、Th时,动态内阻Ri满足:Risc<Ri<Rioc。
S104,在所述动态内阻的变化范围内,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
具体地,当测试电阻的阻值与动态内阻相等时,热电模块的输出功率最大,在已知动态内阻的变化范围的情况下,只需在该变化范围内调节测试电阻的阻值即可,并测量各个测试电阻阻值对应的热电模块的输出功率。
S105,根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
具体地,从测量得到的热电模块输出功率中,可以发现输出功率存在峰值,并由此可以确定热电模块的最大输出功率。
本发明实施例提供的一种测量热电模块最大输出功率的方法、装置及终端设备,根据热电模块输出电流的变化范围,以及热电模块的热端输入热能与热电模块的输出电流的变化关系,确定所述热端输入热能的变化范围,进而确定热电材料平均温度的变化范围和动态内阻的变化范围,在动态内阻的变化范围内,逐次调节测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率,确定热电模块的最大输出功率,缩小了测试电阻的调节范围,测量过程简单方便,节省了大量时间。
进一步地,所述测量热电模块最大输出功率的方法还包括:
根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻。
调节所述测试电阻的初始阻值等于所述平均动态内阻。
具体地,在已知动态内阻的变化范围为Risc<Ri<Rioc时,可以根据Ri-mean=(Rioc+Risc)/2确定平均动态内阻,将平均动态内阻作为测试电阻的初始阻值,从能量输入角度分析,平均动态内阻更接近最大输出功率工况下的动态内阻,通过微调测试电阻的阻值,即可快速匹配热电模块输出功率最大的工况,并测量得到对应冷、热端温度下热电模块的最大输出功率,大大减少了测量时间。
进一步地,所述热电模块的热端输入热能与热电模块的输出电流的变化关系为:
其中,参数A=Rh-(1-η)Rc,参数B=1-0.5η+K[Rh+(1-η)Rc],Tc、Th分别为热电模块冷、热端温度,Qh为热电模块的热端输入热能,η为热电模块的热电转换效率,N为热电模块中的P-N电偶对数,α为热电模块的塞贝克系数,K为热电模块的导热系数,I为热电模块的输出电流。
所述热端输入热能与所述热电材料平均温度的变化关系为:
其中,参数A=Rh-(1-η)Rc,Tj为热电材料的平均温度,Tc、Th分别为热电模块冷、热端温度,Tcj、Thj分别为热电材料的冷、热端结面温度,Qh为热电模块的热端输入热能,Rc、Rh分别为热电模块冷、热端寄生电阻,η为热电模块的热电转换效率。
所述热电材料电导率与所述热电材料平均温度的变化关系为:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度。
所述热电材料电导率与所述动态内阻的变化关系为:
其中,N为热电模块中的P-N电偶对数,hn、hp分别为热电模块的N型、P型热电材料的臂高,ln、lp分别为热电模块的N型、P型热电材料横截面边长,σn、σp分别为热电模块的N型、P型热电材料在平均温度Tj下的电导率,Ri为热电模块在平均温度Tj下的动态内阻。
进一步地,所述在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率包括:
在所述动态内阻的变化范围内,确定所述测试电阻的初始调节方向;
在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
具体地,在动态内阻的阻值变化范围内,确定测试电阻的初始阻值和初始调节方向,这里的初始调节方向为测试电阻的阻值增大或减小的方向,可以以一定的差值逐次调节测试电阻的阻值,测量各个测试电阻阻值下的输出功率,在输出功率出现减小趋势时,表明测试电阻的阻值已经偏离了最大输出功率对应的动态内阻,此时,应该向相反的方向调节测试电阻的阻值,并测量热电模块的输出功率。
进一步地,所述根据测量的热电模块的输出功率确定热电模块的最大输出功率包括:
根据测量的热电模块的输出功率,确定热电模块的输出功率的峰值;
根据热电模块的输出功率的峰值,确定热电模块的最大输出功率。
具体地,热电模块的输出功率随输出电流的增大呈现先增大后减小的态势,同理,输出功率随动态内阻的变化也会呈现先增大后减小的态势,那么测量得到的输出功率的峰值,就是热电模块的最大输出功率。
参阅图3,本发明另一实施例提供了一种测量热电模块最大输出功率的装置,包括:
第一变化范围确定模块301,用于获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与所述输出电流的变化关系,确定所述热端输入热能的变化范围。
第二变化范围确定模块302,用于根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围。
第三变化范围确定模块303,用于根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围。
输出功率测量模块304,用于在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
最大输出功率确定模块305,用于根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
进一步地,所述测量热电模块最大输出功率的装置还包括:
平均动态内阻确定模块,用于根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻。
初始阻值调节模块,用于调节所述测试电阻的初始阻值等于所述平均动态内阻。
进一步地,所述输出功率测量模块304包括:
初始调节方向确定模块,用于在所述动态内阻的变化范围内,确定所述测试电阻的初始调节方向;
第一输出功率测量模块,用于在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
第二输出功率测量模块,用于在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
进一步地,所述测量热电模块最大输出功率的装置还包括:
关系确定模块,用于确定所述热电模块的热端输入热能与热电模块的输出电流的变化关系为:
其中,参数A=Rh-(1-η)Rc,参数B=1-0.5η+K[Rh+(1-η)Rc],Tc、Th分别为热电模块冷、热端温度,Qh为热电模块的热端输入热能,η为热电模块的热电转换效率,N为热电模块中的P-N电偶对数,α为热电模块的塞贝克系数,K为热电模块的导热系数,I为热电模块的输出电流。
确定所述热端输入热能与所述热电材料平均温度的变化关系为:
其中,参数A=Rh-(1-η)Rc,Tj为热电材料的平均温度,Tc、Th分别为热电模块冷、热端温度,Tcj、Thj分别为热电材料的冷、热端结面温度,Qh为热电模块的热端输入热能,Rc、Rh分别为热电模块冷、热端寄生电阻,η为热电模块的热电转换效率。
确定所述热电材料电导率与所述热电材料平均温度的变化关系为:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度。
确定所述热电材料电导率与所述动态内阻的变化关系为:
其中,N为热电模块中的P-N电偶对数,hn、hp分别为热电模块的N型、P型热电材料的臂高,ln、lp分别为热电模块的N型、P型热电材料横截面边长,σn、σp分别为热电模块的N型、P型热电材料在平均温度Tj下的电导率,Ri为热电模块在平均温度Tj下的动态内阻。
进一步地,所述最大输出功率确定模块305包括:
输出功率峰值确定模块,用于根据测量的热电模块的输出功率,确定热电模块的输出功率的峰值。
最大输出功率确定子模块,用于根据热电模块的输出功率的峰值,确定热电模块的最大输出功率。
如图4所示,本发明又一实施例提供了一种测量热电模块最大输出功率的终端设备,该实施例的测量热电模块最大输出功率的终端设备4包括:处理器40、存储器41以及存储在所述存储器41中并可在所述处理器40上运行的计算机程序42,例如测量热电模块最大输出功率的程序。所述处理器30执行所述计算机程序42时实现上述各个测量热电模块最大输出功率的方法实施例中的步骤,例如图1所示的步骤S101至S105。或者,所述处理器40执行所述计算机程序42时实现上述各装置实施例中各模块/单元的功能,例如图3所示模块301至305的功能。
示例性的,所述计算机程序42可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器41中,并由所述处理器40执行,以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序42在所述测量热电模块最大输出功率的终端设备4中的执行过程。例如,所述计算机程序42可以被分割成同步模块、汇总模块、获取模块、返回模块(虚拟装置中的模块),各模块的具体功能如下:
获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与热电模块的输出电流的变化关系,确定所述热端输入热能的变化范围。
根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围。
根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围。
在所述动态内阻的变化范围内,逐次调节热电模块的外接测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
进一步地,各个模块的具体功能还包括:
根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻;
调节所述测试电阻的初始阻值等于所述平均动态内阻。
进一步地,各个模块的具体功能还包括:
所述热电模块的热端输入热能与热电模块的输出电流的变化关系为:
其中,参数A=Rh-(1-η)Rc,参数B=1-0.5η+K[Rh+(1-η)Rc],Tc、Th分别为热电模块冷、热端温度,Qh为热电模块的热端输入热能,η为热电模块的热电转换效率,N为热电模块中的P-N电偶对数,α为热电模块的塞贝克系数,K为热电模块的导热系数,I为热电模块的输出电流。
所述热端输入热能与所述热电材料平均温度的变化关系为:
其中,参数A=Rh-(1-η)Rc,Tj为热电材料的平均温度,Tc、Th分别为热电模块冷、热端温度,Tcj、Thj分别为热电材料的冷、热端结面温度,Qh为热电模块的热端输入热能,Rc、Rh分别为热电模块冷、热端寄生电阻,η为热电模块的热电转换效率。
所述热电材料电导率与所述热电材料平均温度的变化关系为:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度。
所述热电材料电导率与所述动态内阻的变化关系为:
其中,N为热电模块中的P-N电偶对数,hn、hp分别为热电模块的N型、P型热电材料的臂高,ln、lp分别为热电模块的N型、P型热电材料横截面边长,σn、σp分别为热电模块的N型、P型热电材料在平均温度Tj下的电导率,Ri为热电模块在平均温度Tj下的动态内阻。
进一步地,各个模块的具体功能还包括:
在所述动态内阻的变化范围内,确定所述测试电阻的初始调节方向。
在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率。
进一步地,各个模块的具体功能还包括:
根据测量的热电模块的输出功率,确定热电模块的输出功率的峰值;
根据热电模块的输出功率的峰值,确定热电模块的最大输出功率。
所述测量热电模块最大输出功率的终端设备4可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述测量热电模块最大输出功率的终端设备可包括,但不仅限于,处理器40、存储器41。本领域技术人员可以理解,图4仅仅是测量热电模块最大输出功率的终端设备4的示例,并不构成对测量热电模块最大输出功率的终端设备4的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述测量热电模块最大输出功率的终端设备还可以包括输入输出设备、网络接入设备、总线等。
所称处理器40可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器41可以是所述工程项目竣工决算终端设备4的内部存储单元,例如工程项目竣工决算终端设备4的硬盘或内存。所述存储器41也可以是所述工程项目竣工决算终端设备4的外部存储设备,例如所述工程项目竣工决算终端设备4上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(FlashCard)等。进一步地,所述存储器41还可以既包括所述工程项目竣工决算终端设备4的内部存储单元也包括外部存储设备。所述存储器41用于存储所述计算机程序以及所述工程项目竣工决算终端设备所需的其他程序和数据。所述存储器41还可以用于暂时地存储已经输出或者将要输出的数据。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述***中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本发明所提供的实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的***实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括是电载波信号和电信信号。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (8)
1.一种测量热电模块最大输出功率的方法,其特征在于,包括以下步骤:
获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与所述输出电流的变化关系,确定所述热端输入热能的变化范围;
根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围;
根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围;
在所述动态内阻的变化范围内,确定热电模块的外接测试电阻的初始调节方向,在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率,在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
2.根据权利要求1所述的测量热电模块最大输出功率的方法,其特征在于,还包括:
根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻;
调节所述测试电阻的初始阻值等于所述平均动态内阻。
3.根据权利要求1所述的测量热电模块最大输出功率的方法,其特征在于,
所述热电模块的热端输入热能与热电模块的输出电流的变化关系为:
其中,参数A=Rh-(1-η)Rc,参数B=1-0.5η+K[Rh+(1-η)Rc],Tc、Th分别为热电模块冷、热端温度,Qh为热电模块的热端输入热能,η为热电模块的热电转换效率,N为热电模块中的P-N电偶对数,α为热电模块的塞贝克系数,K为热电模块的导热系数,I为热电模块的输出电流;
所述热端输入热能与所述热电材料平均温度的变化关系为:
其中,参数A=Rh-(1-η)Rc,Tj为热电材料的平均温度,Tc、Th分别为热电模块冷、热端温度,Tcj、Thj分别为热电材料的冷、热端结面温度,Qh为热电模块的热端输入热能,Rc、Rh分别为热电模块冷、热端寄生电阻,η为热电模块的热电转换效率;
所述热电材料电导率与所述热电材料平均温度的变化关系为:
其中,σn、σp分别为热电模块的N型、P型热电材料的电导率,Tj为热电材料平均温度;
所述热电材料电导率与所述动态内阻的变化关系为:
4.根据权利要求1所述的测量热电模块最大输出功率的方法,其特征在于,所述根据测量的热电模块的输出功率确定热电模块的最大输出功率包括:
根据测量的热电模块的输出功率,确定热电模块的输出功率的峰值;
根据热电模块的输出功率的峰值,确定热电模块的最大输出功率。
5.一种测量热电模块最大输出功率的装置,其特征在于,包括:
第一变化范围确定模块,用于获取热电模块的输出电流的变化范围,根据所述输出电流的变化范围,以及热电模块的热端输入热能与所述输出电流的变化关系,确定所述热端输入热能的变化范围;
第二变化范围确定模块,用于根据所述热端输入热能的变化范围,以及所述热端输入热能与热电模块的热电材料平均温度的变化关系,确定所述热电材料平均温度的变化范围;
第三变化范围确定模块,用于根据所述热电材料平均温度的变化范围,以及热电材料电导率与热电材料平均温度的变化关系,热电材料电导率与热电模块的动态内阻的变化关系,确定所述动态内阻的变化范围;
输出功率测量模块,用于在所述动态内阻的变化范围内,确定热电模块的外接测试电阻的初始调节方向,在所述初始调节方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率,在测量的热电模块的输出功率出现减小趋势时,在所述初始调节方向的相反方向上,逐次调节所述测试电阻的阻值,并测量各个测试电阻阻值下热电模块的输出功率;
最大输出功率确定模块,用于根据测量的热电模块的输出功率,确定热电模块的最大输出功率。
6.根据权利要求5所述的测量热电模块最大输出功率的装置,其特征在于,还包括:
平均动态内阻确定模块,用于根据所述动态内阻的变化范围,确定所述动态内阻的变化范围内的平均动态内阻;
初始阻值调节模块,用于调节测试电阻的初始阻值等于所述平均动态内阻。
7.一种测量热电模块最大输出功率的终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至4任意一项所述测量热电模块最大输出功率的方法的步骤。
8.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至4任意一项所述测量热电模块最大输出功率的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711344727.1A CN108152617B (zh) | 2017-12-15 | 2017-12-15 | 测量热电模块最大输出功率的方法、装置及终端设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711344727.1A CN108152617B (zh) | 2017-12-15 | 2017-12-15 | 测量热电模块最大输出功率的方法、装置及终端设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108152617A CN108152617A (zh) | 2018-06-12 |
CN108152617B true CN108152617B (zh) | 2020-05-26 |
Family
ID=62467177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711344727.1A Expired - Fee Related CN108152617B (zh) | 2017-12-15 | 2017-12-15 | 测量热电模块最大输出功率的方法、装置及终端设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108152617B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1890821A (zh) * | 2003-12-02 | 2007-01-03 | ***纪念研究所 | 热电装置及其应用 |
CN1975448A (zh) * | 2006-12-11 | 2007-06-06 | 中国科学院上海硅酸盐研究所 | 一种热电发电元器件性能测量装置及其方法 |
CN102540099A (zh) * | 2012-02-06 | 2012-07-04 | 天津大学 | 微型温差电池热电转换效率测试***及方法 |
CN202421371U (zh) * | 2011-12-22 | 2012-09-05 | 武汉理工大学 | 一种可控热电模块冷热端温度的功率测试实验装置 |
CN103323780A (zh) * | 2013-05-20 | 2013-09-25 | 南京航空航天大学 | 一种温差材料热电转换性能测试***及方法 |
CN104007139A (zh) * | 2014-06-10 | 2014-08-27 | 中国华能集团清洁能源技术研究院有限公司 | 热电模块的测试***及热电模块的测试方法 |
CN105203940A (zh) * | 2015-09-28 | 2015-12-30 | 中国科学院上海硅酸盐研究所 | 一种热电元件可靠性评价***及方法 |
CN106093654A (zh) * | 2016-07-11 | 2016-11-09 | 浙江大学 | 热电模组热电转换效率的测试装置及其测试方法 |
CN107290638A (zh) * | 2017-06-14 | 2017-10-24 | 复旦大学 | 一种多功能热电器件测试*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6856136B2 (en) * | 2002-05-22 | 2005-02-15 | Massachusetts Institute Of Technology | Thermoelectric device test structure |
-
2017
- 2017-12-15 CN CN201711344727.1A patent/CN108152617B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1890821A (zh) * | 2003-12-02 | 2007-01-03 | ***纪念研究所 | 热电装置及其应用 |
CN1975448A (zh) * | 2006-12-11 | 2007-06-06 | 中国科学院上海硅酸盐研究所 | 一种热电发电元器件性能测量装置及其方法 |
CN202421371U (zh) * | 2011-12-22 | 2012-09-05 | 武汉理工大学 | 一种可控热电模块冷热端温度的功率测试实验装置 |
CN102540099A (zh) * | 2012-02-06 | 2012-07-04 | 天津大学 | 微型温差电池热电转换效率测试***及方法 |
CN103323780A (zh) * | 2013-05-20 | 2013-09-25 | 南京航空航天大学 | 一种温差材料热电转换性能测试***及方法 |
CN104007139A (zh) * | 2014-06-10 | 2014-08-27 | 中国华能集团清洁能源技术研究院有限公司 | 热电模块的测试***及热电模块的测试方法 |
CN105203940A (zh) * | 2015-09-28 | 2015-12-30 | 中国科学院上海硅酸盐研究所 | 一种热电元件可靠性评价***及方法 |
CN106093654A (zh) * | 2016-07-11 | 2016-11-09 | 浙江大学 | 热电模组热电转换效率的测试装置及其测试方法 |
CN107290638A (zh) * | 2017-06-14 | 2017-10-24 | 复旦大学 | 一种多功能热电器件测试*** |
Non-Patent Citations (4)
Title |
---|
Beat the Deviations in Estimating Maximum Power of Thermoelectric Modules;Junling Gao et.al;《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》;20131031;第62卷(第10期);第2725-2729页 * |
Experimental analysis of peak power output of a thermoelectric liquid-to-liquid generator under an increasing electrical load resistance;Frédéric J. Lesage et.al;《Energy Conversion and Management》;20131231(第66期);第98-105页 * |
大功率热电转换关键技术及测试方法研究;高俊岭;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20160115(第01期);第C042-182页 * |
热电模块最大输出功率的估值方法;高俊岭 等;《华南理工大学学报(自然科学版)》;20130430;第41卷(第4期);第54-57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108152617A (zh) | 2018-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Afolabi et al. | Analysis of the load flow problem in power system planning studies | |
Zagirnyak et al. | Analysis of instantaneous power components of electric circuit with a semiconductor element | |
CN109342989B (zh) | 一种电能表的误差分析方法及装置 | |
CN108063456B (zh) | 分布式光伏发电并网规划方法及终端设备 | |
CN116031777A (zh) | 一种核电用配电柜温升智能控制方法及*** | |
CN108896803B (zh) | 基于温度补偿的电能表计量精度优化方法 | |
CN109000825B (zh) | 一种含谐波电流的电缆及其终端线芯温度计算方法 | |
CN115236392A (zh) | 多特征量的电能计量方法、装置、终端及存储介质 | |
Zhang et al. | A probabilistic assessment method for voltage stability considering large scale correlated stochastic variables | |
CN108152617B (zh) | 测量热电模块最大输出功率的方法、装置及终端设备 | |
US11169188B2 (en) | Low-observability matrix completion | |
CN103280413B (zh) | 获取晶圆的电阻温度系数的离散度的工业实现方法 | |
CN111611710A (zh) | 一种光伏电池等效建模方法 | |
Sun et al. | Static voltage stability analysis based on multi-dimensional holomorphic embedding method | |
CN116165547A (zh) | 锂电池的电压计算与修正方法及***、装置、存储介质 | |
CN113346842B (zh) | 一种太阳能电池串联电阻阻值的在线计算方法及装置 | |
Sarniak | Modeling of photovoltaic modules in Simulink and Simscape packages of Matlab software | |
CN109521349B (zh) | 半导体发电模块最大输出功率的测量方法及*** | |
CN113467546A (zh) | 温度控制方法、装置以及测温箱 | |
CN112082786A (zh) | 发电机组空气预热器正常运行性能的预测方法及装置 | |
CN114337215A (zh) | 功率变换电路的功率降额方法、装置、终端及存储介质 | |
CN111983411B (zh) | 多指栅型晶体管热阻测试方法、装置及终端设备 | |
CN114089144B (zh) | 二极管结参数测量方法及测量*** | |
CN114679132B (zh) | 一种光伏逆变器电热工况仿真方法、装置及存储介质 | |
KR101722826B1 (ko) | 태양광 패널의 모델링 방법 및 이를 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 기록매체 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20200526 Termination date: 20211215 |
|
CF01 | Termination of patent right due to non-payment of annual fee |