CN105468815A - 一种防雷产品的压敏电阻热致损毁评估方法及系统 - Google Patents

一种防雷产品的压敏电阻热致损毁评估方法及系统 Download PDF

Info

Publication number
CN105468815A
CN105468815A CN201510742048.4A CN201510742048A CN105468815A CN 105468815 A CN105468815 A CN 105468815A CN 201510742048 A CN201510742048 A CN 201510742048A CN 105468815 A CN105468815 A CN 105468815A
Authority
CN
China
Prior art keywords
model
piezoresistor
lightning
surge protector
simplified
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
CN201510742048.4A
Other languages
English (en)
Other versions
CN105468815B (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.)
GUANGDONG MIG TECHNOLOGY Inc
Hefei University of Technology
Original Assignee
GUANGDONG MIG TECHNOLOGY Inc
Hefei University of Technology
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 GUANGDONG MIG TECHNOLOGY Inc, Hefei University of Technology filed Critical GUANGDONG MIG TECHNOLOGY Inc
Priority to CN201510742048.4A priority Critical patent/CN105468815B/zh
Publication of CN105468815A publication Critical patent/CN105468815A/zh
Application granted granted Critical
Publication of CN105468815B publication Critical patent/CN105468815B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明涉及防雷产品热致损毁领域,具体公开了一种防雷产品的压敏电阻热致损毁评估方法及系统。所述方法包括:S1、建立雷电流脉冲简化模型;S2、建立压敏电阻等效电路模型;S3、建立压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;S4、将雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与实体模型进行耦合,以获得多场耦合模型;S5、对多场耦合模型和有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性、安全性,以及对电涌保护器后置保护装置的灵敏性进行评估。

Description

一种防雷产品的压敏电阻热致损毁评估方法及系统
技术领域
本发明涉及防雷产品热致损毁领域,特别是涉及一种防雷产品的压敏电阻热致损毁评估方法及系统。
背景技术
雷电具有高电压,大电流和瞬间电磁冲击特性。雷击电涌(Surge)对各类微电子设备和电子信息系统(包括计算机系统、通信系统、电子控制系统、无线电系统以及电力电子装置等)安全运行带来极大的威胁。在现代防雷系统中,电涌保护器(SurgeProtectiveDevice,SPD)通过自身电子元器件的非线性伏安特性,有效地将瞬时过电压限制在安全范围之内,将大电流快速地泄放到大地,保证设备或系统不受冲击而损坏。
电涌保护器自身安全性和耐冲击性能受到关注,当受到多脉冲大电流重复冲击时,电涌保护器中的金属氧化物压敏电阻(MetalOxideVaristors,MOV)热吸收和热释放失去平衡,MOV热量迅速积聚,超过其极限耐受能量,导致电涌保护器热致损毁,严重影响设备或系统的安全,甚至引发火灾。通常,雷电脉冲冲击性试验只针对单个电流脉冲(如8/20μs)进行电涌保护器的耐冲击性检查,忽略重复电流脉冲及复合波形叠加电流脉冲的冲击次数及冲击间隔的影响,导致试验结果不够精确。而针对电涌保护器后置保护装置(如后置机械脱扣装置,熔断器等),现有的检测技术是基于破坏性试验,极大地浪费了生产资源,严重影响了生产效率。
发明内容
本发明的目的在于克服现有技术的不足,提供一种防雷产品的压敏电阻热致损毁评估方法,以可以非破坏性地评估不同类型的雷电流脉冲冲击电涌保护器的热积聚量以及电涌保护器后置保护装置精准脱扣动作的可靠性。
为了达到上述目的,本发明采用的技术方案是:
一种防雷产品的压敏电阻热致损毁评估方法,所述防雷产品具有电涌保护器与电涌保护器后置保护装置,所述方法包括:
S1、对雷电流脉冲进行简化建模,获得雷电流脉冲简化模型;
S2、对电涌保护器中的压敏电阻进行等效电路建模,获得压敏电阻等效电路模型;
S3、对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
S4、将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
S5、对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
优选地,所述步骤S1中,在对雷电流脉冲进行简化建模时,通过选取可积的双指数函数作为雷电流函数脉冲表达式建立雷电流脉冲简化模型,这是由于雷电流的电磁场表达式中既有与通道电流有关的项(感应场项)和与电流导数有关的项(辐射场项),又有与电流对时间积分有关的项(静电场项),因此选取可积的双指数函数作为雷电流函数脉冲表达式,所述可积的双指数函数为:
i ( t ) = I 0 ω ( e - α t - e - β t ) ;
式中,ω为峰值修正因子, ω = e - αt p - e - βt p ; 其中, t p = l n β α β - α , α,β为波形拟合常数,I0为短路电流峰值。
具体地,对于10/350μs的电流脉冲,选取下列函数作为雷电流函数脉冲表达式,建立雷电流脉冲简化模型:
i(t)=A1Im1(e-αt-e-βt);
式中,A1=1.025,α=2.05×10-3(μs)-1,β=0.564(μs)-1
对于8/20μs的电流脉冲,则选取下列作为雷电流函数脉冲表达式,建立雷电流脉冲简化模型:
i ( t ) = A 2 I m 2 e - t τ
式中,A2=0.01243,τ=3.11μs。
优选地,所述步骤S2中,在对电涌保护器中的压敏电阻进行等效电路建模时,根据压敏电阻的伏安特性曲线,建立分段等效电路模型;所述分段等效电路模型包括预击穿区等效电路模型,击穿区等效电路模型,以及回升区等效电路模型。
优选地,所述步骤S3中,利用CAD软件对压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模,并通过CAE软件赋予所述实体模型各个部分的材料特性及相关的性能参数。
为了简化建模,便于分析外加雷电电流脉冲产生的热效应对压敏电阻及电涌保护器后置保护装置的影响,在利用CAD软件对压敏电阻进行三维实体建模过程中,忽略次要影响因素对其进行建模,所述次要影响因素包括以下因素中的任意一个或任意多个:压敏电阻本身微观晶粒的各向异性、致密度不均匀性、气孔,以及压敏电阻老化、劣化系数。
优选地,所述步骤S4中,利用CAE软件导入所述雷电流脉冲简化模型、压敏电阻等效电路模型,并在CAE软件中添加电磁加热和传热这两个物理场,通过对电磁加热和传热这两个物理过程进行分析,进一步简化所述雷电流脉冲简化模型与压敏电阻等效电路模型。
具体地,在进一步简化雷电流脉冲简化模型时,根据雷电流的类型作进一步的简化建模,对于单一波形重复电流脉冲,选用以下函数作为雷电流函数脉冲表达式,建立雷电流单一波形重复多脉冲简化模型:
i ( t ) = 0.01243 I m 2 e - t 3.11 c o s ( 2 k π Δ t ) ;
其中,k=0,Δt,2Δt,...,(n-1)Δt;
式中,Im2为短路电流峰值,Δt为脉冲冲击时间间隔,n为脉冲冲击次数;
对于复合波形叠加电流脉冲,选用以下函数作为雷电流函数脉冲表达式建立雷电流复合波形叠加多脉冲简化模型:
i ( t ) = 1.025 I m 1 a ( e - 2.05 × 10 - 3 t - e - 0.564 t ) cos ( 2 k 1 π Δt 1 ) + 0.01243 I m 2 ( 1 - a ) e - t 3.11 cos ( 2 k 2 π Δt 2 ) ;
其中,k1=0,Δt1,2Δt1,...,(n-1)Δt1,k2=0,Δt2,2Δt2,...,(n-1)Δt2式中,Im1为其中一电流脉冲I1的短路电流峰值,a为电流脉冲I1所占百分比系数,Δt1为电流脉冲I1的脉冲冲击时间间隔,Im2为另一电流脉冲I2的短路电流峰值,Δt2为电流脉冲I2的脉冲冲击时间间隔,n为脉冲冲击次数;
在进一步简化压敏电阻等效电路模型过程中,通过计算出压敏电阻的吸收能量,并考虑压敏电阻的热弛豫现象来分析压敏电阻的受热情况来进行简化建模。
其中,在计算压敏电阻的吸收能量时,利用公式
W = ∫ 0 t C u ( t ) i ( t ) d t
来进行计算,式中,W为吸收能量,C为波形相关常数,u(t)为施加于压敏电阻上的瞬时电压,i(t)为施加于压敏电阻上的瞬时电流;即可以从施加于压敏电阻上的电流、电压的幅值和持续时间计算出压敏电阻的吸收能量。
由于压敏电阻在受到电流脉冲的冲击作用,产生热效应,其存在着热弛豫现象,即压敏电阻还未开始散热,下一个电流脉冲已开始作用于压敏电阻,因此,在对压敏电阻进行简化模型过程中将压敏电阻的热弛豫现象这一因素考虑在内。其中,热弛豫时间τ可表示成:
T = x 2 α t h ,
式中,x为压敏电阻元件厚度,αth为热扩散系数,其定义为:
α t h = K ρC p ,
式中,K为压敏电阻的热导率,ρ为压敏电阻的密度,Cp为压敏电阻的比热。
优选地,所述步骤S5中,利用CAE软件对步骤S4中建立的多场耦合模型和有限元数值模型进行求解,并对求解结果进行分析,借此来对防雷产品避免热致损毁的可靠性、安全性,以及对电涌保护器后置保护装置的灵敏性进行评估;
在评估过程中,通过对所述实体模型施加特定的单一波形重复电流脉冲和/或复合波形叠加电流脉冲,观察分析压敏电阻以及所述电涌保护器后置保护装置关键点的温度随时间的变化,并通过改变雷电流脉冲的参数,进一步分析在不同雷电流脉冲的作用下,分析压敏电阻的受热情况和电涌保护器后置保护装置关键点的温度变化。此外,还可以通过CAD改变所述实体模型的参数,观察特定规格压敏电阻在特定电流脉冲作用下所能耐受的冲击次数。
基于上述方法的发明构思,本发明还提供了一种防雷产品的压敏电阻热致损毁评估系统,其应用于具有电涌保护器与电涌保护器后置保护装置的防雷产品,所述系统包括:
雷电流脉冲简化模型建立模块,用于对雷电流脉冲进行简化建模,以获得雷电流脉冲简化模型;
压敏电阻等效电路模型建立模块,用于对电涌保护器中的压敏电阻进行等效电路建模,以获得压敏电阻等效电路模型;
实体模型与有限元数值模型建立模块,用于对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,以获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
多场耦合模型建立模块,用于将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
求解评估模块,用于对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
与现有技术相比,本发明具有的有益效果有:
(1)利用本发明提供的评估方法,可以在压敏电阻设计初期,进行粗略的压敏电阻雷电流脉冲模拟,来初步验证设计的合理性;
(2)本发明提供的评估方法,涉及参数可变的重复电流脉冲和复合电流脉冲冲击性实验的模拟,很好地弥补了传统单电流脉冲冲击性试验的局限性;
(3)本发明提供的评估方法,可以快速估计特定规格的压敏电阻所能耐受特定脉冲冲击的极限次数。
(4)本发明提供的评估方法,还可以验证电涌保护器后置保护装置(如熔断器、机械脱扣装置等)设计的合理性和灵敏性。
附图说明
图1为本发明实施例中所述防雷产品的压敏电阻热致损毁评估方法的流程示意图;
图2为本发明实施例中所述压敏电阻等效电路模型的示意图;
图3为本发明实施例中压敏电阻(氧化锌压敏电阻)与电涌保护器后置保护装置(低温焊锡机械脱扣装置)的实体模型的结构示意图;
图4为本发明实施例中对压敏电阻(氧化锌压敏电阻)及电涌保护器后置保护装置(低温焊锡机械脱扣装置)的有限元数值模型进行求解,得到在特定时间后的实体模型的温度分布云图;
图5为在传热过程中,低温焊锡机械脱扣装置中的低温焊锡的关键点的温度变化曲线图;
图6为本发明实施例中所述防雷产品的压敏电阻热致损毁评估系统的模块示意框图。
其中,图3中的附图标号为:1-脱扣簧片,2-低温焊锡,3-导电铜电极,4-氧化锌压敏电阻本体。
具体实施方式
下面结合实施例参照附图进行详细说明,以便对本发明的技术特征及优点进行更深入的诠释。
一种防雷产品的压敏电阻热致损毁评估方法,所述防雷产品具有电涌保护器与电涌保护器后置保护装置,如图1所示,所述方法包括:
S1、对雷电流脉冲进行简化建模,获得雷电流脉冲简化模型;
S2、对电涌保护器中的压敏电阻进行等效电路建模,获得压敏电阻等效电路模型;
S3、对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
S4、将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
S5、对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
在步骤S1中,考虑到雷电流脉冲的多样性和复杂性,因此对其进行简化建模。在对雷电流脉冲进行简化建模时,通过选取可积的双指数函数作为雷电流函数脉冲表达式建立雷电流脉冲简化模型,所述可积的双指数函数为:
i ( t ) = I 0 ω ( e - α t - e - β t ) - - - ( 1 ) ;
(1)式中,ω为峰值修正因子, ω = e - αt p - e - βt p ; 其中, t p = l n β α β - α , α,β为波形拟合常数,I0为短路电流峰值。
具体地,对于10/350μs的电流脉冲,(1)式可简化为:
i ( t ) = 1.025 I m 1 ( e - 2.05 × 10 - 3 t - e - 0.564 t ) - - - ( 2 ) ;
(2)式中,Im1为10/350μs电流脉冲的短路电流峰值,本实施例选取Im1=10kA。
对于8/20μs的电流脉冲,(1)式可简化为:
i ( t ) = 0.01243 I m 2 e - t 3.11 - - - ( 3 ) ;
(3)式中Im2为8/20μs电流脉冲的短路电流峰值,本实施例选取Im2=10kA。
本发明实施例中,选取氧化锌(ZnO)压敏陶瓷电阻(以下简称为氧化锌压敏电阻)作为电涌保护器中的压敏电阻(即将氧化锌压敏电阻作为限压泄放电流元件),选取低温焊锡机械脱扣装置作为电涌保护器后置保护装置,来对本发明作具体的描述。
在步骤S2中,在对压敏电阻进行等效电路建模时,根据压敏电阻的伏安特性曲线,并分析氧化锌压敏电阻中的氧化锌的非线性伏安特性进行建模,建立如附图2中(a)所示的简化的压敏电阻等效电路模型(其中,L为接线电感)。
其中,根据压敏电阻工作时处于预击穿区(小电流区)、击穿区(中电流区)、回升区(大电流区)的工作状态,可对建立分段等效电路模型:所述分段等效电路模型包括预击穿区等效电路模型,击穿区等效电路模型,以及回升区等效电路模型。
当压敏电阻处于预击穿区时,压敏电阻的非线性很低,可以忽略附图2(a)中的非线性电阻R2和串联电阻R3,得到如附图2(b)所示的简化电路模型(即预击穿区等效电路模型),在本实施例中,并联电阻R1=8MΩ,电极电容C=430pF;
当压敏电阻处于击穿区时,其击穿区等效电路模型如附图2(c)所示,其中非线性电阻R2随温度可变,高温条件下相当于短路;
当压敏电阻处于回升区时,非线性电阻R2相当于短路,得到如附图2(d)所示的简化电路模型(即预击穿区等效电路模型),本实施例中,串联电阻R3=5Ω。
在步骤S3中,利用CAD(ComputerAidedDesign,计算机辅助设计)软件对压敏电阻(本实施例优选为取氧化锌压敏电阻)以及所述电涌保护器后置保护装置(本实施例优选为低温焊锡机械脱扣装置)进行三维实体建模,建立压敏电阻及电涌保护器后置保护装置的实体模型,所述实体模型如附图3所示;同时,通过CAE(ComputerAidedEngineering,计算机辅助工程)软件赋予所述实体模型各个部分的材料特性及相关的性能参数。
本实施例中,主要通过分析对实体模型中的压敏电阻施加雷电流脉冲而产生的热效应对氧化锌压敏电阻的影响及低温焊锡机械脱扣装置的灵敏程度;因此,在建模过程中忽略氧化锌压敏电阻本身微观晶粒的各向异性、致密度不均匀性、气孔以及压敏电阻老化,劣化系数等次要影响因素中的一个或任意多个。
在所述步骤S4中,优选利用CAE软件导入所述雷电流脉冲简化模型、压敏电阻等效电路模型,并在CAE软件中添加电磁加热和传热这两个物理场,通过对电磁加热和传热这两个物理过程进行分析,进一步简化所述雷电流脉冲简化模型与压敏电阻等效电路模型。
具体到本实施例中,在通过分析对实体模型中的压敏电阻施加雷电流脉冲而产生的热效应而对氧化锌压敏电阻的影响及低温焊锡机械脱扣装置的灵敏程度时,分别通过对氧化锌压敏电阻施加8/20μs单电流脉冲,8/20μs单一波形重复电流脉冲,8/20μs与10/350μs复合波形叠加电流脉冲而产生的热效应。
因此,在进一步简化雷电流脉冲简化模型时,可根据雷电流的类型作进一步的简化建模。具体地,对于单一波形重复电流脉冲,以8/20μs的电流脉冲为例,选用以下函数作为雷电流函数脉冲表达式,建立雷电流单一波形重复多脉冲简化模型:
i ( t ) = 0.01243 I m 2 e - t 3.11 c o s ( 2 k π Δ t ) - - - ( 4 ) ;
其中,k=0,Δt,2Δt,...,(n-1)Δt;
(4)式中,Im2为短路电流峰值,Δt为脉冲冲击时间间隔,n为脉冲冲击次数。
对于复合波形叠加电流脉冲,以8/20μs和10/350μs的复合波形叠加电流脉冲为例,选用以下函数作为雷电流函数脉冲表达式建立雷电流复合波形叠加多脉冲简化模型:
i ( t ) = 1.025 I m 1 a ( e - 2.05 × 10 - 3 t - e - 0.564 t ) cos ( 2 k 1 π Δt 1 ) + 0.01243 I m 2 ( 1 - a ) e - t 3.11 cos ( 2 k 2 π Δt 2 ) - - - ( 5 ) ;
其中,k1=0,Δt1,2Δt1,...,(n-1)Δt1,k2=0,Δt2,2Δt2,...,(n-1)Δt2(2)式中,Im1为其中一电流脉冲I1(本实施例中电流脉冲I1为10/350μs电流脉冲)的短路电流峰值,a为电流脉冲I1所占百分比系数,Δt1为电流脉冲I1的脉冲冲击时间间隔,Im2为另一电流脉冲I2(本实施例中电流脉冲I2为8/20μs电流脉冲)的短路电流峰值,Δt2为电流脉冲I2的脉冲冲击时间间隔,n为脉冲冲击次数;
而在进一步简化压敏电阻等效电路模型过程中,通过计算出压敏电阻的吸收能量,并考虑压敏电阻的热弛豫现象来分析压敏电阻的受热情况来进行简化建模。
其中,在计算压敏电阻的吸收能量时,利用公式
W = ∫ 0 t C u ( t ) i ( t ) d t = C U I t - - - ( 6 )
来进行计算压敏电阻的吸收能量。(6)式中,W为吸收能量,C为波形相关常数,u(t)为施加于压敏电阻上的瞬时电压,i(t)为施加于压敏电阻上的瞬时电流;即可以从施加于压敏电阻上的电流、电压的幅值和持续时间计算出压敏电阻的吸收能量。
而结合上述步骤S2的分析,在小电流区(预击穿区),压敏电阻接近于绝缘体,根据公式(6),此阶段产生的热量可以忽略不计。
在中电流区(击穿区),压敏电阻呈现良好的导通性能,当电流不断增大时,压敏电阻两端的电压几乎不发生变化,从而很好地达到了限压的效果,因此,计算此阶段产生的热量时,公式(6)中的U可以设为定值。
大电流区(回升区)压敏电阻的伏安特性曲线又呈线性分布,利用欧姆定律,公式(6)可转化成:
W = ∫ 0 t C u ( t ) i ( t ) d t = ∫ 0 t i ( t ) 2 R d t - - - ( 7 ) ; 其中,R可以近似为固定值。
对氧化锌压敏电阻在不同电流脉冲区域的热量积累进行简化,将小电流区、中电流区、大电流区所积累的热量分别记作W1、W2、W3,则压敏电阻所产生的总热量W总=W1+W2+W3。
而在考虑压敏电阻的热弛豫现象时,热弛豫时间τ可表示成:
τ = x 2 α t h - - - ( 8 ) ,
(8)式中,x为压敏电阻元件厚度,αth为热扩散系数,其定义为:
α t h = K ρC p - - - ( 9 ) ,
(9)式中,K为压敏电阻的热导率,ρ为压敏电阻的密度,Cp为压敏电阻的比热。
具体到本实施例,氧化锌压敏电阻的热导率K=0.057J/(s·cm·K),密度ρ=5.4g/cm3,比热Cp=0.89J/(g·K);本实施例中,氧化锌压敏电阻的元件横截面积A=342=1156mm2,其元件厚度x=0.5cm;因此,可求得热扩散系数αth=0.0119cm2/s,氧化锌压敏陶瓷电阻热弛豫时间τ=21s。
在所述步骤S5中,优选利用CAE软件对步骤S4中建立的多场耦合模型和有限元数值模型进行求解,并对求解结果进行分析,借此来对防雷产品避免热致损毁的可靠性、安全性,以及对电涌保护器后置保护装置的灵敏性进行评估。当然,对多场耦合模型和有限元数值模型的求解,在实际求解过程中其不仅仅限于CAE软件进行模拟求解,还可以用其他可计算出求解结果的软件进行求解。
在评估过程中,通过对所述实体模型(主要为实体模型中的压敏电阻)施加特定的单一波形重复电流脉冲和/或复合波形叠加电流脉冲,观察分析压敏电阻以及所述电涌保护器后置保护装置关键点的温度随时间的变化,并通过改变雷电流脉冲的参数,进一步分析在不同雷电流脉冲的作用下,分析压敏电阻的受热情况和电涌保护器后置保护装置关键点的温度变化。此外,还可以通过CAD改变所述实体模型的参数,观察特定规格压敏电阻在特定电流脉冲作用下所能耐受的冲击次数。
具体到本实施例,在对氧化锌压敏电阻施加8/20μs单电流脉冲,其短路电流峰值为Im2=10kA,则氧化锌压敏电阻吸收能量为:W1≈0,W2=260J,W3=0(Im=10kA的脉冲几乎到达不了大电流区),因此W总1=260J;而单脉冲不考虑热弛豫现象,因此氧化锌压敏电阻的温升ΔT1=W总1/(Cp·ρ·V)=9.37K。
对于短路电流峰值Im1为10kA,冲击时间间隔为10s,冲击次数为10的8/20μs重复电流脉冲,由上述可知,热弛豫时间τ(τ=21s)约为冲击时间间隔Δt的两倍,即当电流脉冲开始第三次冲击时,电流脉冲第一次冲击产生的热量开始散失,W总2=2600J;而W散1=K·A·ΔT·Δt/x=1317J(其中K为热导率,A为元件横截面积(这里取A=1156mm2),ΔT为温差,x为元件厚度(这里取x=0.5cm)),ΔW总2=W总2-W散1=1283J。因此氧化锌压敏电阻的温升ΔT2=46.24K。
对于短路电流峰值为10kA,冲击时间间隔为10s,冲击次数为10的8/20μs、10/350μs复合波形叠加电流脉冲,其分析方法与重复电流脉冲类似,得到W总3=4000J,W散2=1317J,ΔW总3=W总3-W散2=2683J。因此氧化锌压敏电阻的温升ΔT3=96.70K。
考虑到本实施例选用的氧化锌压敏电阻的工作极限温度为403K,工作初始温度为25℃(298K),因此本实施例压敏电阻的极限温升为ΔT=105K;由此可知,本实施例中,氧化锌压敏电阻的三组温升均未超过其极限温升,因此三组雷电电流脉冲均在安全范围之内,即在此三种雷电流脉冲的作用下,防雷产品符合安全性的相关要求,防雷产品的安全性与可靠性好。而观察8/20μs、10/350μs复合波形电流脉冲模拟结果(ΔT3=96.70K)与极限温升(ΔT=105K)较接近,若再增加脉冲次数或者增大短路电流峰值,其温升很可能超过极限温升,因此氧化锌压敏电阻安全性不可控。
对于本实施例选取的低温焊锡机械脱扣装置,考虑到氧化锌压敏电阻的工作极限温度为403K,定义氧化锌压敏电阻和周围环境的温差,对步骤S3中所述的有限元模型进行求解,可分别得到压敏电阻开始散热特定时间后的实体模型的温度分布云图,如附图4所示;其中,图4中的(a)、(b)(c)分别为压敏电阻开始散热0s,0.5s,1s后,所述实体模型的温度分布云图。对实体模型中低温焊锡关键点(本实施例选取三个关键点)进行温度的信息收集,可得到如附图5所示的温度变化曲线图。从图5中可以看出,温度可以快速地传递到低温焊锡上,而所述防雷产品的脱扣基本温度设置在345K,因此,所述低温焊锡机械脱扣装置可以快速地,灵敏地完成压敏电阻的机械脱扣,即说明此低温焊锡机械脱扣装置的设计合理,其机械脱扣的灵敏性可符合防雷产品的安全性的相关要求。
本发明实施例提供的防雷产品的压敏电阻热致损毁评估系统,如附图6所示,所述系统包括:
雷电流脉冲简化模型建立模块,用于对雷电流脉冲进行简化建模,以获得雷电流脉冲简化模型;
压敏电阻等效电路模型建立模块,用于对电涌保护器中的压敏电阻进行等效电路建模,以获得压敏电阻等效电路模型;
实体模型与有限元数值模型建立模块,用于对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,以获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
多场耦合模型建立模块,用于将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
求解评估模块,用于对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
优选地,所述雷电流脉冲简化模型建立模块在对雷电流脉冲进行简化建模时,通过选取可积的双指数函数作为雷电流函数脉冲表达式建立雷电流脉冲简化模型;
所述可积的双指数函数为:
式中,ω为峰值修正因子, ω = e - αt p - e - βt p ; 其中, t p = l n β α β - α , α,β为波形拟合常数,I0为短路电流峰值。
本发明实施例提供的防雷产品的压敏电阻热致损毁评估系统,其工作方法/过程或应用原理如上一实施例所述的评估方法的工作方法/过程或应用原理相同,在此不再一一赘述。
上述实施例中提到的内容为本发明较佳的实施方式,并非是对本发明的限定,在不脱离本发明构思的前提下,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应视为等效的置换方式,其均在本发明的保护范围之内。

Claims (10)

1.一种防雷产品的压敏电阻热致损毁评估方法,所述防雷产品具有电涌保护器与电涌保护器后置保护装置,其特征在于,所述方法包括:
S1、对雷电流脉冲进行简化建模,获得雷电流脉冲简化模型;
S2、对电涌保护器中的压敏电阻进行等效电路建模,获得压敏电阻等效电路模型;
S3、对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
S4、将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
S5、对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
2.根据权利要求1所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:
所述步骤S1中,在对雷电流脉冲进行简化建模时,通过选取可积的双指数函数作为雷电流函数脉冲表达式建立雷电流脉冲简化模型;
所述可积的双指数函数为:
式中,ω为峰值修正因子,其中,α,β为波形拟合常数,I0为短路电流峰值。
3.根据权利要求1或2所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:所述步骤S3中,利用CAD软件对压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模,并通过CAE软件赋予所述实体模型各个部分的材料特性及相关的性能参数。
4.根据权利要求3所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:在利用CAD软件对压敏电阻进行三维实体建模过程中,忽略次要影响因素对其进行建模,所述次要影响因素包括以下因素中的任意一个或任意多个:压敏电阻本身微观晶粒的各向异性、致密度不均匀性、气孔,以及压敏电阻老化、劣化系数。
5.根据权利要求3所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:所述步骤S4中,利用CAE软件导入所述雷电流脉冲简化模型、压敏电阻等效电路模型,并在CAE软件中添加电磁加热和传热这两个物理场,通过对电磁加热和传热这两个物理过程进行分析,进一步简化所述雷电流脉冲简化模型与压敏电阻等效电路模型。
6.根据权利要求5所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:
所述步骤S4中,在进一步简化雷电流脉冲简化模型时,根据雷电流的类型作进一步的简化建模,对于单一波形重复电流脉冲,选用以下函数作为雷电流函数脉冲表达式,建立雷电流单一波形重复多脉冲简化模型:
i ( t ) = 0.01243 I m 2 e - t 3.11 c o s ( 2 k π Δ t ) , 其中,k=0,Δt,2Δt,...,(n-1)Δt;
式中,Im2为短路电流峰值,Δt为脉冲冲击时间间隔,n为脉冲冲击次数;
对于复合波形叠加电流脉冲,选用以下函数作为雷电流函数脉冲表达式建立雷电流复合波形叠加多脉冲简化模型:
i ( t ) = 1.025 I m 1 a ( e - 2.05 × 10 - 3 t - e - 0.564 t ) cos ( 2 k 1 π Δt 1 ) + 0.01243 I m 2 ( 1 - a ) e - t 3.11 c o s ( 2 k 2 π Δt 2 ) ;
其中,k1=0,Δt1,2Δt1,...,(n-1)Δt1,k2=0,Δt2,2Δt2,...,(n-1)Δt2式中,Im1为其中一电流脉冲I1的短路电流峰值,a为电流脉冲I1所占百分比系数,Δt1为电流脉冲I1的脉冲冲击时间间隔,Im2为另一电流脉冲I2的短路电流峰值,Δt2为电流脉冲I2的脉冲冲击时间间隔,n为脉冲冲击次数;
在进一步简化压敏电阻等效电路模型过程中,通过计算出压敏电阻的吸收能量,并考虑压敏电阻的热弛豫现象来分析压敏电阻的受热情况来进行简化建模;
其中,在计算压敏电阻的吸收能量时,利用公式来进行计算,式中,W为吸收能量,C为波形相关常数,u(t)为施加于压敏电阻上的瞬时电压,i(t)为施加于压敏电阻上的瞬时电流。
7.根据权利要求6所述的防雷产品MOV热致损毁评估模型,其特征在于:所述步骤S5中,利用CAE软件对步骤S4中建立的多场耦合模型和有限元数值模型进行求解,并对求解结果进行分析,借此来对防雷产品避免热致损毁的可靠性、安全性,以及对电涌保护器后置保护装置的灵敏性进行评估;
在评估过程中,通过对所述模型施加特定的单一波形重复电流脉冲和/或复合波形叠加电流脉冲,观察分析压敏电阻以及所述电涌保护器后置保护装置关键点的温度随时间的变化,并通过改变电流脉冲的参数,进一步分析在不同电流脉冲的作用下,分析压敏电阻的受热情况和电涌保护器后置保护装置关键点的温度变化。
8.根据权利要求1所述的防雷产品的压敏电阻热致损毁评估方法,其特征在于:
所述步骤S2中,在对电涌保护器中的压敏电阻进行等效电路建模时,根据压敏电阻的伏安特性曲线,建立分段等效电路模型;所述分段等效电路模型包括预击穿区等效电路模型,击穿区等效电路模型,以及回升区等效电路模型。
9.一种防雷产品的压敏电阻热致损毁评估系统,应用于具有电涌保护器与电涌保护器后置保护装置的防雷产品,其特征在于,所述系统包括:
雷电流脉冲简化模型建立模块,用于对雷电流脉冲进行简化建模,以获得雷电流脉冲简化模型;
压敏电阻等效电路模型建立模块,用于对电涌保护器中的压敏电阻进行等效电路建模,以获得压敏电阻等效电路模型;
实体模型与有限元数值模型建立模块,用于对所述压敏电阻以及所述电涌保护器后置保护装置进行三维实体建模与有限元数值建模,以获得压敏电阻及电涌保护器后置保护装置的实体模型与有限元数值模型;
多场耦合模型建立模块,用于将所述雷电流脉冲简化模型、压敏电阻等效电路模型中的物理场量与所述实体模型进行耦合,以获得压敏电阻及电涌保护器后置保护装置的感应电磁场和热力场的多场耦合模型;
求解评估模块,用于对所述多场耦合模型和所述有限元数值模型进行求解,在获得特定时间后压敏电阻及电涌保护器后置保护装置关键节点的温度变化及其温度峰值在时间和空间上的规律的基础上,对防雷产品避免热致损毁的可靠性与安全性进行评估,和/或对电涌保护器后置保护装置的灵敏性进行评估。
10.根据权利要求9所述的防雷产品的压敏电阻热致损毁评估系统,其特征在于:
所述雷电流脉冲简化模型建立模块在对雷电流脉冲进行简化建模时,通过选取可积的双指数函数作为雷电流函数脉冲表达式建立雷电流脉冲简化模型;
所述可积的双指数函数为:
式中,ω为峰值修正因子,其中,α,β为波形拟合常数,I0为短路电流峰值。
CN201510742048.4A 2015-11-02 2015-11-02 一种防雷产品的压敏电阻热致损毁评估方法及系统 Expired - Fee Related CN105468815B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510742048.4A CN105468815B (zh) 2015-11-02 2015-11-02 一种防雷产品的压敏电阻热致损毁评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510742048.4A CN105468815B (zh) 2015-11-02 2015-11-02 一种防雷产品的压敏电阻热致损毁评估方法及系统

Publications (2)

Publication Number Publication Date
CN105468815A true CN105468815A (zh) 2016-04-06
CN105468815B CN105468815B (zh) 2019-04-26

Family

ID=55606511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510742048.4A Expired - Fee Related CN105468815B (zh) 2015-11-02 2015-11-02 一种防雷产品的压敏电阻热致损毁评估方法及系统

Country Status (1)

Country Link
CN (1) CN105468815B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106707071A (zh) * 2017-02-28 2017-05-24 四川瑞云信通科技有限公司 用于计算电源防雷装置剩余使用寿命的系统及其方法
CN108957169A (zh) * 2018-05-15 2018-12-07 大连理工大学 基于极化去极化电荷差的氧化锌避雷器服役状态评价方法
CN109002586A (zh) * 2018-06-25 2018-12-14 国网湖南省电力有限公司 一种避雷器温度计算方法和系统
CN112001106A (zh) * 2020-08-26 2020-11-27 上海大学 一种基于有限元法的压敏电阻微观结构电热特性仿真方法
CN114184842A (zh) * 2021-11-15 2022-03-15 西南交通大学 一种基于能量吸收的多脉冲下ZnO电阻片性能评估方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2785116Y (zh) * 2005-02-02 2006-05-31 浙江德力西电器股份有限公司 电涌保护器的热脱离机构
CN101887790A (zh) * 2010-06-13 2010-11-17 深圳科安达电子科技股份有限公司 防雷用压敏电阻

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2785116Y (zh) * 2005-02-02 2006-05-31 浙江德力西电器股份有限公司 电涌保护器的热脱离机构
CN101887790A (zh) * 2010-06-13 2010-11-17 深圳科安达电子科技股份有限公司 防雷用压敏电阻

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
M GLASA: "The MOV computer models for thermal - Electric analysis", 《JOURNAL OF ELECTRICAL ENGINEERING》 *
于永站, 谭昌柏, 张金泉: "塑壳断路器热脱扣特性的仿真与试验分析", 《电器与能效管理技术》 *
孙建虎, 刘颖芳, 尹平: "雷电流数学模型分析", 《后勤工程学院学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106707071A (zh) * 2017-02-28 2017-05-24 四川瑞云信通科技有限公司 用于计算电源防雷装置剩余使用寿命的系统及其方法
CN108957169A (zh) * 2018-05-15 2018-12-07 大连理工大学 基于极化去极化电荷差的氧化锌避雷器服役状态评价方法
CN109002586A (zh) * 2018-06-25 2018-12-14 国网湖南省电力有限公司 一种避雷器温度计算方法和系统
CN109002586B (zh) * 2018-06-25 2023-04-18 国网湖南省电力有限公司 一种避雷器温度计算方法和系统
CN112001106A (zh) * 2020-08-26 2020-11-27 上海大学 一种基于有限元法的压敏电阻微观结构电热特性仿真方法
CN114184842A (zh) * 2021-11-15 2022-03-15 西南交通大学 一种基于能量吸收的多脉冲下ZnO电阻片性能评估方法
CN114184842B (zh) * 2021-11-15 2022-08-09 西南交通大学 一种基于能量吸收的多脉冲下ZnO电阻片性能评估方法

Also Published As

Publication number Publication date
CN105468815B (zh) 2019-04-26

Similar Documents

Publication Publication Date Title
CN105468815B (zh) 一种防雷产品的压敏电阻热致损毁评估方法及系统
WO2020114082A1 (zh) 多重连续雷电流分量作用下碳纤维复合材料损伤计算方法
CN111507034B (zh) 计算冲击负荷下避雷器温度场时变特性方法和系统
CN103808426B (zh) 电缆线芯温度的间接测量方法
CN103163421A (zh) 一种测试熔断器抗浪涌电流冲击能力的装置以及方法
Topcagic et al. Electrothermal and overload performance of metal-oxide varistors
CN109002586A (zh) 一种避雷器温度计算方法和系统
CN108254712A (zh) 一种电能表电源回路保护特性仿真分析方法及装置
CN110543717A (zh) 一种gis设备温升特性仿真模型的构建方法及装置
Zhou et al. Numerical modelling of MOV with Voronoi network and finite element method
Illias et al. Determination of surge arrester discharge energy using finite element analysis method
Drumea et al. Analysis of current carrying capacity of silver-based conductive pastes for PCB repair
CN105069258A (zh) 一种芯片设计可靠性的评估方法及装置
Zhou et al. Research on electrothermal characteristics of metal oxide varistor based on multi‐physical fields
Kim Failure prediction of metal oxide varistor using nonlinear surge look-up table based on experimental data
He et al. Experimental validation of MOA simulation models for energy absorption estimation under different impulse currents
Liu et al. Numerical inversion analysis on front‐face temperature rise of Al alloy suffered from long continuing current in lightning
Fu et al. Investigation on temperature behavior of CFRP during lightning strike using experiment and simulation
Bhuyan et al. Simulations of lightning impulse residual voltage test of surge arresters in matlab-simulink
Seyyedbarzegar et al. Application of finite element method for electro‐thermal modeling of metal oxide surge arrester
Bendakir et al. Towards the prospection of an optimal thermal response of ZnO surge arrester in HV power system
Marghescu et al. Investigation on current capabilities of Ni-based conductive pastes for PCB repair
Chen et al. Investigation on the accuracy of the junction temperature determined by the VCE (T) method under different load pulse durations in power cycling test
Meng et al. A thermal model for electrothermal simulation of power modules
Kuang et al. Temperature distribution in a ZnO arrester subjected to multiple current impulses

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

Termination date: 20201102

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