(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 202011180585.1
(22)申请日 2020.10.29
(71)申请人 西安理工大学
地址 710048 陕西省西安市金花南路5号
(72)发明人 张嘉伟 冯凯 付庚
(74)专利代理机构 北京国昊天诚知识产权代理
有限公司 11315
代理人 杨洲
(51)Int.Cl.
G06F 30/20(2020.01)
G06F 113/16(2020.01)
G06F 119/04(2020.01)
(54)发明名称一种基于微元化物理模型的电缆寿命预测方法(57)摘要本发明提供了一种基于微元化物理模型的电缆寿命预测方法，将电缆试样分别在90℃、110℃和130℃下进行热老化实验，在不同老化时间下取样测量其断裂伸长率保留率数值，再将实验所得数据代入建模方程，得到不同的t0和v值，对t 0和v 值进行数据拟合，以确定最佳拟合值。对三组不同温度下对应的t 0和v值再进行数据拟合，分别得到t 0和v关于温度的函数关系，回带建模方程，得到断裂伸长保留率关于温度和时间的函数。本发明将失效标准的力学参数断裂伸长率保留率作为时间和温度相关的函数，一方面可以代入断裂伸长率保留率为50％时，求出电缆在额定工作温度下的电缆寿命，一方面也可直观显示绝
缘材料性能受温度和老化时间影响的下降趋势。权利要求书1页 说明书5页 附图2页CN 112446136 A 2021.03.05
C N 112446136
A
1.一种基于微元化物理模型的电缆寿命预测方法，其特征在于，包括以下步骤：
S1：利用微分法，将电缆试样分成若干个单位立方体块，再将每个单位立方体块分为m 个子立方体，m为正整数；
S2：假定在m个子立方体中，n个子立方体被降解，令每个单位立方体中退化部分体积比为V d ，则每个子立方体块边长为d ，可表示为d＝(V d /n)1/3；令k是沿x、y或z方向上的任意轴线退化的子立方体的数目，可表示为k＝n 1/3；
S3：设定电缆试件为一个单位长度，所以分母为1，具体如下：
1)试样老化前断裂伸长率EAB before ：
EAB before ＝α×1-1/1 (1)
其中：α-伸长因子，为大于1的实数；
2)试样老化后断裂伸长率EAB after ：
未退化部分其伸长因子a保持为原始值，退化部分伸长因子令为a d ，退化的个数为n 1/3，则退化部分总长度为∑d k ＝n 1/3×(V d /n)1/3，未退化部分长度为1-∑d k ；
EAB after ＝[α×(1-∑d k )+αd ×(∑d k )-1]/1
(2)3)断裂伸长率保留率σ：
σ＝EAB after /EAB before ≈1-(V d )1/3 (3)
4)v是由断裂伸长率的定义表征，可以用指数分布的累积分布函数表示：
σ＝1-(1-e -vt )1/3 (4)
其中，v表示降解速率；
假定t 0时刻开始发生降解反应，方程(4)可转化为如下(5)和(6)：
σ＝1，t≤t 0
(5)
S4：在不同老化时间下取样测量其断裂伸长率保留率数值，再将所得数据代入建模方程，得到不同的t 0和v值，对t 0和v值进行数据拟合，以确定最佳拟合值。
2.根据权利要求1所述的基于微元化物理模型的电缆寿命预测方法，其特征在于：所述步骤S4将电缆试样分别在90℃、110℃和130℃下进行热老化实验。
权 利 要 求 书1/1页CN 112446136 A
一种基于微元化物理模型的电缆寿命预测方法
技术领域
[0001]本发明属于电缆检测技术领域，具体涉及一种基于微元化物理模型的电缆寿命预测方法。
背景技术
[0002]随着电气设备越来越先进，越来越普遍，电力电缆的使用量也越来越多。电缆线路是城市电网的重要组成部分，随着电缆服役时间和数量的增长，其运行状态以及剩余使用寿命引起了人们的关注。
[0003]目前针对电缆的检测方法通常是巡视检查和绝缘电阻检测，采用这些手段可以大致分析电缆的运行状态，但是并不能确定其老化程度和使用寿命等信息。而橡胶绝缘电力电缆在运行中，热老化是导致其老化的关键因素，绝缘材料在受热条件下会发生降解反应，使得电缆绝缘性能逐渐降低。因此，开展电力电缆的热老化寿命分析工作尤为重要，评估电缆寿命，可以为确定电缆最佳检修时机与更换策略提供支持，提高电缆利用率，减少电力运维成本，降低经济损失。也可以预防电缆因绝缘老化失效而导致的安全隐患，保障电力系统的正常运行。
[0004]专利文献CN107590314 A公开的一种基于matlab的电缆寿命数据评估方法，本发明基于matlab软件，设定电缆老化的函数符合阿列纽斯曲线，将其对数化，得到不同温度下电缆使用寿命的对数与温度的倒数的线性关系，通过对输入的三组不同温度对应下的试验时间和断裂伸长率保留率数据进行拟合，求出不同温度对应下断裂伸长率保留率为50％时寿命时间，生成三组x,y值，再进行拟合，根据拟合直线得到热寿命与温度值的函数关系。但基于阿列纽斯曲线得出的最终结果只显示了电缆的热寿命与温度值
间的函数关系，未建立断裂伸长率保留率与老化时间和温度的函数关系，对受老化时间和温度影响下绝缘材料力学性能的下降趋势并未涉及，不够全面。
[0005]专利文献CN 106771891 A公开的一种电热联合老化的交联聚乙烯电缆剩余寿命评价方法中，通过对待评估交联聚乙烯电缆进行差式扫描量热实验，得到热老化化学反应活化能，利用电缆剩余寿命评估公式，结合电缆实际运行工况，计算得出电缆剩余寿命。但该方法无法从理论上建立常规法所具备的绝缘材料老化性能与老化程度之间的关系，并缺少重复性实验对其稳定性进行验证，具有一定的局限性。
[0006]需要注意的是，本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
发明内容
[0007]本发明目的在于提供了一种基于微元化物理模型的电缆寿命预测方法，建立一个基于物理的预测模型，将失效标准的力学参数断裂伸长率保留率作为时间和温度相关的函数，更能直观反映老化时间及温度对电缆绝缘材料性能的影响。
[0008]为实现上述目的本发明采用如下技术方案：
[0009]该基于微元化物理模型的电缆寿命预测方法，包括以下步骤：
[0010]S1：利用微分法，将电缆试样分成若干个单位立方体块，再将每个单位立方体块分为m个子立方体，m为正整数；
[0011]S2：假定在m个子立方体中，n个子立方体被降解，令每个单位立方体中退化部分体积比为V d，则每个子立方体块边长为d，可表示为d＝(V d/n)1/3；令k是沿x、y或z方向上的任意轴线退化的子立方体的数目，可表示为k＝n1/3；
[0012]S3：设定电缆试件为一个单位长度，所以分母为1，具体如下：
[0013]1)试样老化前断裂伸长率EAB before：
[0014]EAB before＝α×1-1/1 (1)
[0015]其中：α-伸长因子，为大于1的实数；
[0016]2)试样老化后断裂伸长率EAB after：
[0017]未退化部分其伸长因子a保持为原始值，退化部分伸长因子令为a d，退化的个数为n1/3，则退化部分总长度为∑d k＝n1/3×(V d/n)1/3，未退化部分长度为1-∑d k；
[0018]EAB after＝[α×(1-∑d k)+αd×(∑d k)-1]/1 (2)
[0019]3)断裂伸长率保留率σ：
[0020]σ＝EAB after/EAB before≈1-(V d)1/3 (3)
[0021]4)v是由断裂伸长率的定义表征，可以用指数分布的累积分布函数表示：[0022]σ＝1-(1-e-vt)1/3 (4)
[0023]其中，v表示降解速率；
[0024]假定t0时刻开始发生降解反应，方程(4)可转化为如下(5)和(6)：
[0025]σ＝1，t≤t0 (5)
[0026]
[0027]S4：在不同老化时间下取样测量其断裂伸长率保留率数值，再将所得数据代入建模方程，得到不同的t0和v值，对t0和v值进行数据拟合，以确定最佳拟合值。
[0028]进一步地，上述步骤S4将电缆试样分别在90℃、110℃和130℃下进行热老化实验。[0029]本发明的有益效果：
[0030]本发明利用微元法对电缆绝缘材料的断裂伸长保留率在不同温度下，受老化时间影响的下降曲线进行建模，建立一种基于物理的预测模型，将失效标准的力学参数断裂伸长率保留率作为时间和温度相关的函数，一方面可以代入断裂伸长率保留率为50％时，求出电缆在额定工作温度下的电缆寿命，一方面也可直观显示绝缘材料性能受温度和老化时间影响的下降趋势。
附图说明
[0031]图1是电缆试样微元化模型图；
[0032]图2是本发明的方法流程图。
[0033]附图标记说明：
[0034]1-电缆哑铃状试样微元化模型；2-单位子立方体块。
具体实施方式
[0035]现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的范例；相反，提供这些实施方式使得本发明将更加全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征或特性可以以任何合适的方式结合在一个或更多实施方式中。
[0036]下面结合附图1和2对本发明进行详细说明。
[0037]利用微分法，将电缆试样分成若干个单位立方体块，电缆试样采用电缆哑铃状试样微元化模型1，再将每个单位立方体块分为m个单位子立方体2，m为大的正整数。假定在m 个子立方体中，n个子立方体被降解。令每个单位立方体中退化部分体积比为V d，则每个子立方体块边长为d，可表示为d＝(V d/n)1/3。在单位立方体中，假定退化是均匀的，令k是沿x、y或z方向上的任意轴线退化的子立方体的数目，可表示为k＝n1/3。为了简化推导过程，本文对试件设定了一个单位长度，所以分母为1。
[0038]1)试样老化前断裂伸长率(EAB before)：
[0039]EAB before＝α×1-1/1 (1)
[0040]α：伸长因子，为大于1的实数
[0041]2)试样老化后断裂伸长率(EAB after)：
[0042]未退化部分其伸长因子a保持为原始值，退化部分伸长因子令为a d(退化部分无延展性，a d无限接近于1)，退化的个数为n1/3，则退化部分总长度为∑d k＝n1/3×(V d/n)1/3，未退化部分长度为1-∑d k。
[0043]EAB after＝[α×(1-∑d k)+αd×(∑d k)-1]/1 (2)
[0044]αd：接近于1
[0045]3)断裂伸长率保留率σ：
[0046]σ＝EAB after/EAB before≈1-(V d)1/3 (3)
[0047]表1、降解量的变化趋势
[0048]
[0049]在反应控制的降解中，V d遵循表1所示的下降速率v，其中v并不是由化学反应来测量，而是由断裂伸长率的定义表征，累计的虚拟退化部分为列(b)的总和，可以用指数分布的累积分布函数表示：
[0050]即σ＝1-(1-e-vt)1/3 (4)
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议