CN115618696A - 基于扩展有限元的三芯电缆接头温度场的计算方法及系统 - Google Patents
基于扩展有限元的三芯电缆接头温度场的计算方法及系统 Download PDFInfo
- Publication number
- CN115618696A CN115618696A CN202211631943.5A CN202211631943A CN115618696A CN 115618696 A CN115618696 A CN 115618696A CN 202211631943 A CN202211631943 A CN 202211631943A CN 115618696 A CN115618696 A CN 115618696A
- Authority
- CN
- China
- Prior art keywords
- cable
- finite element
- temperature
- formula
- crack
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims description 86
- 239000004020 conductor Substances 0.000 claims description 39
- 238000006073 displacement reaction Methods 0.000 claims description 37
- 238000010438 heat treatment Methods 0.000 claims description 27
- 239000012530 fluid Substances 0.000 claims description 26
- 238000012546 transfer Methods 0.000 claims description 22
- 238000004458 analytical method Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 19
- 230000004907 flux Effects 0.000 claims description 15
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 229910000831 Steel Inorganic materials 0.000 claims description 6
- 230000002500 effect on skin Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 239000002184 metal Substances 0.000 claims description 6
- 239000010959 steel Substances 0.000 claims description 6
- 230000007613 environmental effect Effects 0.000 claims description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000017525 heat dissipation Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000013021 overheating Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 229920003020 cross-linked polyethylene Polymers 0.000 description 1
- 239000004703 cross-linked polyethylene Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及电力电缆温度计算技术领域,公开了一种基于扩展有限元的三芯电缆接头温度场的计算方法及系统,其方法公开了通过利用有限元仿真软件构建电缆有限元仿真模型,并基于扩展有限元法对电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型,以提高温度场分布的求解精度,还设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场,从而能够准确地计算并反映直埋敷设下三芯电缆接头的温度场分布情况。
Description
技术领域
本发明涉及电力电缆温度计算技术领域,尤其涉及基于扩展有限元的三芯电缆接头温度场的计算方法及系统。
背景技术
随着我国经济和电力行业的发展,城市配电网的地下电缆系统已经成为我国电网改造工程中的重要部分,因此电力电缆温升的准确计算和及时预警对于电力系统的安全运行有着重要的意义。相较于电缆本体,电缆接头结构复杂、现场安装工艺欠妥、接触电阻的存在导致其发生过热故障的可能性增大,这在一定程度上影响电缆的正常运行。不同的电缆有不同的敷设条件,不同的敷设条件下电缆最外层的环境条件不同、散热条件不同,这对电缆本体及接头的温度产生较大的影响。三芯电缆常用于三相交流电网中的中亚配电线路和中性点接地的三相四线制系统,在城市配电网中应用广泛。相比于其他敷设条件,三芯电缆直埋敷设的优点为:敷设方便,节省材料和人工。但是由于土壤,电缆更不容易散热,长期工作在较高温的环境下会加速电缆接头的绝缘老化,影响其正常运行。因此有必要对三芯电缆接头温度场分布进行准确计算,把握其温度发展规律,做到温升过高及时预警,避免出现过热故障甚至火灾等安全事故。
公开号为CN 110083908A的中国发明专利申请公开了基于有限元分析的电缆缆芯温度预测方法,利用有限元分析与神经网络的有机结合,通过电缆的结构参数、材料热参数、敷设条件以及边界条件建立电缆的有限元模型,施加环境参数以及载荷进行有限元计算,得到电缆的温度场分布云图,选取电缆外护套及缆芯所对应的节点,查询相应的温度数据,获得样本数据,然后对获得的原始样本数据进行归一化处理,确定输入变量后,利用神经网络对数据进行训练,建立缆芯温度的预测模型,最后将预处理后的数据输入到训练好的神经网络中,就可以得到所求的缆芯温度预测值。通过提出的温度预测方法,可较准确的预测缆芯温度,为实现电网从传统预防性维护到主动预测性管理的转变提供一种有效途径。
而在现有技术中,多数是基于有限元分析法对于交联聚乙烯绝缘电缆的温度场进行了计算,而针对三芯电缆接头结构较为复杂的情况,尤其应用于出现裂纹或划痕现象的电缆接头,现有技术仍难以准确地计算并反映直埋敷设下三芯电缆接头的温度场分布情况。
发明内容
本发明提供了一种基于扩展有限元的三芯电缆接头温度场的计算方法及系统,解决了现有技术难以准确地计算并反映直埋敷设下三芯电缆接头的温度场分布情况的技术问题。
有鉴于此,本发明第一方面提供了一种基于扩展有限元的三芯电缆接头温度场的计算方法,包括以下步骤:
根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,所述电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
优选地,基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型的步骤具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
对所述电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
优选地,设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度的步骤具体包括:
利用二维泊松方程求解电缆接头的温度分布为:
根据电缆接头的温度分布设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
优选地,将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场的步骤具体包括:
通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
第二方面,本发明提供了一种基于扩展有限元的三芯电缆接头温度场的计算系统,包括:
建模模块,用于根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,所述电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
扩展有限元模块,用于基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
边界设置模块,用于设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
温度场分析模块,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
优选地,所述扩展有限元模块具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
网格剖分模块,用于对所述电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
优选地,所述边界设置模块具体包括:
温度分布求解模块,用于利用二维泊松方程求解电缆接头的温度分布为:
边界计算模块,用于根据电缆接头的温度分布设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
优选地,所述温度场分析模块具体包括:
介质损耗计算模块,用于通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
温度场分析模块,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
从以上技术方案可以看出,本发明具有以下优点:
本发明通过考虑三芯电缆散热的特点,利用有限元仿真软件构建电缆有限元仿真模型,并基于扩展有限元法对电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型,以提高温度场分布的求解精度,还设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场,从而能够准确地计算并反映直埋敷设下三芯电缆接头的温度场分布情况。
附图说明
图1为本发明实施例提供的一种基于扩展有限元的三芯电缆接头温度场的计算方法的流程图;
图2为本发明实施例提供的一种基于扩展有限元的三芯电缆接头温度场的计算系统的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了便于理解,请参阅图1,本发明提供的一种基于扩展有限元的三芯电缆接头温度场的计算方法,包括以下步骤:
S1、根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
可以理解的是,本实施例考虑到周围环境对电缆的温升影响,进而在进行有限元模型构建过程中,还对电缆直埋敷设下的环境结构进行建模。
其中,有限元仿真软件可以采用ANSYS,并设置三芯电缆结构的各种材料、结构尺寸和导热系数。
S2、基于扩展有限元法对电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
需要说明的是,对于带有裂纹的部位,采用扩展有限元法进行处理,同时进行网格剖分,对于不带裂痕或划痕的部分,可以采用普通有限元法处理以降低计算量。
S3、设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
S4、将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将热源加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
其中,在有限元仿真软件中,将热源加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,即可进行温度场分析,得到电缆接头的温度场。
需要说明的是,本实施例提供了一种于扩展有限元的三芯电缆接头温度场的计算方法,通过考虑三芯电缆散热的特点,利用有限元仿真软件构建电缆有限元仿真模型,并基于扩展有限元法对电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型,以提高温度场分布的求解精度,还设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场,从而能够准确地计算并反映直埋敷设下三芯电缆接头的温度场分布情况。
在一个具体实施例中,步骤S2具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
需要说明的是,扩展有限元法的基本思想是单位分解法,对于不连续问题,采用引入不连续函数来进行局部逼近,从而解决不连续问题。水平集方法(Level Set Method)是一种用于界面追踪和形状建模的数值技术,优点是可以在笛卡尔网格(Cartesian grid)上对演化中的曲线曲面进行数值计算而不必对曲线曲面参数化,同时可以方便地追踪单位的拓扑结构改变。
在理论研究中,电缆接头的裂纹或划痕与剖分的有限元网格完全独立,扩展有限元法可以通过上述理论实现裂纹的全自动扩展模拟,具有良好的逼近性能。单位分解法指的是任意函数都可以用一组局部函数和来表示。而连续的独立裂纹单元离散为各个节点位移,再通过有限元仿真软件实现连续性的计算。
同时,上述所涉及到的向量和公式中所涉及的其他函数是涉及裂纹的有限元计算中的原理性的内容,具体是通过成熟的有限元仿真软件计算实现的。
S202、对电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
在一个具体实施例中,步骤S3具体包括:
S301、利用二维泊松方程求解电缆接头的温度分布为:
S302、根据电缆接头的温度分布设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
在一个具体实施例中,步骤S4具体包括:
S401、通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
S402、将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将热源加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
以上为本发明提供的一种基于扩展有限元的三芯电缆接头温度场的计算方法的实施例的详细描述,以下为本发明提供的一种基于扩展有限元的三芯电缆接头温度场的计算系统的实施例的详细描述。
为了便于理解,请参阅图2,本发明提供了一种基于扩展有限元的三芯电缆接头温度场的计算系统,包括:
建模模块100,用于根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
扩展有限元模块200,用于基于扩展有限元法对电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
边界设置模块300,用于设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
温度场分析模块400,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将热源加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
在一个具体实施例中,扩展有限元模块具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
网格剖分模块,用于对电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
在一个具体实施例中,边界设置模块具体包括:
温度分布求解模块,用于利用二维泊松方程求解电缆接头的温度分布为:
边界计算模块,用于根据电缆接头的温度分布设置含有电缆裂纹部分的电缆有限元仿真模型的边界条件,边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
在一个具体实施例中,温度场分析模块具体包括:
介质损耗计算模块,用于通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
温度场分析模块,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将热源加载至含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本发明所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.基于扩展有限元的三芯电缆接头温度场的计算方法,其特征在于,包括以下步骤:
根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,所述电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
2.根据权利要求1所述的基于扩展有限元的三芯电缆接头温度场的计算方法,其特征在于,基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型的步骤具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
对所述电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
3.根据权利要求1所述的基于扩展有限元的三芯电缆接头温度场的计算方法,其特征在于,设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度的步骤具体包括:
利用二维泊松方程求解电缆接头的温度分布为:
根据电缆接头的温度分布设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
4.根据权利要求3所述的基于扩展有限元的三芯电缆接头温度场的计算方法,其特征在于,将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场的步骤具体包括:
通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
5.基于扩展有限元的三芯电缆接头温度场的计算系统,其特征在于,包括:
建模模块,用于根据三芯电缆的实际结构参数,采用有限元仿真软件构建三芯电缆的有限元仿真模型,得到电缆有限元仿真模型,所述电缆有限元仿真模型包括电缆直埋敷设下的环境结构、电缆缆芯、电缆接头、导体屏蔽、绝缘层、绝缘屏蔽、绕包带、钢铠套和外护套;
扩展有限元模块,用于基于扩展有限元法对所述电缆有限元仿真模型中的电缆裂纹部分进行定义,并进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型;
边界设置模块,用于设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度;
温度场分析模块,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场。
6.根据权利要求5所述的基于扩展有限元的三芯电缆接头温度场的计算系统,其特征在于,所述扩展有限元模块具体包括:
式1~5中,x表示独立裂纹单元中任一点的位置坐标,是任意函数,是表示函数的局部函数,表示形函数;表示节点间断附加自由度向量;表示奇异附加自由度向量;表示常规节点位移向量;表示裂纹位移不连续特性间断函数;表示裂纹位移不连续特性间断函数;表示裂纹在所选极坐标系下的半径;为裂纹在所选极坐标系下的夹角;
网格剖分模块,用于对所述电缆有限元仿真模型进行网格剖分,得到含有电缆裂纹部分的电缆有限元仿真模型。
7.根据权利要求5所述的基于扩展有限元的三芯电缆接头温度场的计算系统,其特征在于,所述边界设置模块具体包括:
温度分布求解模块,用于利用二维泊松方程求解电缆接头的温度分布为:
边界计算模块,用于根据电缆接头的温度分布设置所述含有电缆裂纹部分的电缆有限元仿真模型的边界条件,所述边界条件包括已知边界温度、已知边界法向热流密度、已知对流换热系数和流体温度分别为:
已知边界温度为,
已知边界法向热流密度为,
已知对流换热系数为,
流体温度为,
8.根据权利要求7所述的基于扩展有限元的三芯电缆接头温度场的计算系统,其特征在于,所述温度场分析模块具体包括:
介质损耗计算模块,用于通过下式计算电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗为:
式11~13中,为电缆内部缆芯导体发热产生的介质损耗,I为电流,R为交流电阻,为绝缘层介质损耗,f为频率,C为绝缘层电容,U为电压,δ为介质损耗角,λs为金属护套的损耗系数,为20℃时的电阻,α为温度系数,t为导体温度,ys为集合肤效应因数,yp为邻近效应因数;
温度场分析模块,用于将电缆内部缆芯导体发热产生的介质损耗和绝缘层介质损耗作为热源,将所述热源加载至所述含有电缆裂纹部分的电缆有限元仿真模型中的电缆接头中,并进行温度场分析,得到电缆接头的温度场分布云图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211631943.5A CN115618696A (zh) | 2022-12-19 | 2022-12-19 | 基于扩展有限元的三芯电缆接头温度场的计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211631943.5A CN115618696A (zh) | 2022-12-19 | 2022-12-19 | 基于扩展有限元的三芯电缆接头温度场的计算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115618696A true CN115618696A (zh) | 2023-01-17 |
Family
ID=84880350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211631943.5A Pending CN115618696A (zh) | 2022-12-19 | 2022-12-19 | 基于扩展有限元的三芯电缆接头温度场的计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115618696A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115795980A (zh) * | 2023-02-06 | 2023-03-14 | 广东电网有限责任公司中山供电局 | 基于多尺度有限元的电缆接头稳态温度场计算方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104731898A (zh) * | 2015-03-18 | 2015-06-24 | 华南理工大学 | 一种10kV三芯电缆有限元热分析方法 |
CN105205302A (zh) * | 2015-04-08 | 2015-12-30 | 辽宁达能电气股份有限公司 | 基于光纤测温主机的电缆动态流量计算方法 |
CN110083908A (zh) * | 2019-04-19 | 2019-08-02 | 陕西科技大学 | 基于有限元分析的电缆缆芯温度预测方法 |
CN115017772A (zh) * | 2022-06-14 | 2022-09-06 | 石家庄铁道大学 | 一种基于电场和温度场分布的三维电缆接头缺陷的分析方法和装置 |
-
2022
- 2022-12-19 CN CN202211631943.5A patent/CN115618696A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104731898A (zh) * | 2015-03-18 | 2015-06-24 | 华南理工大学 | 一种10kV三芯电缆有限元热分析方法 |
CN105205302A (zh) * | 2015-04-08 | 2015-12-30 | 辽宁达能电气股份有限公司 | 基于光纤测温主机的电缆动态流量计算方法 |
CN110083908A (zh) * | 2019-04-19 | 2019-08-02 | 陕西科技大学 | 基于有限元分析的电缆缆芯温度预测方法 |
CN115017772A (zh) * | 2022-06-14 | 2022-09-06 | 石家庄铁道大学 | 一种基于电场和温度场分布的三维电缆接头缺陷的分析方法和装置 |
Non-Patent Citations (3)
Title |
---|
张洪麟 等: ""基于有限元法的地下电缆群温度场及载流量的仿真计算"", 《高压电器》 * |
李录贤 等: ""扩展有限元法(XFEM)及其应用"", 《力学进展》 * |
陈梦曦: ""10kV三芯电缆接头温度场的仿真研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115795980A (zh) * | 2023-02-06 | 2023-03-14 | 广东电网有限责任公司中山供电局 | 基于多尺度有限元的电缆接头稳态温度场计算方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jullien et al. | Advanced modeling of crosstalk between an unshielded twisted pair cable and an unshielded wire above a ground plane | |
Patel et al. | Proximity-aware calculation of cable series impedance for systems of solid and hollow conductors | |
Dubitsky et al. | Comparison of finite element analysis to IEC-60287 for predicting underground cable ampacity | |
CN115618696A (zh) | 基于扩展有限元的三芯电缆接头温度场的计算方法及系统 | |
Tatematsu | A technique for representing lossy thin wires and coaxial cables for FDTD-based surge simulations | |
Shazly et al. | Thermal analysis of high‐voltage cables with several types of insulation for different configurations in the presence of harmonics | |
CN107315860A (zh) | 多根交联聚乙烯绝缘电缆共沟时温度场有限元计算方法 | |
Holyk et al. | Power cable rating calculations-a historical perspective [history] | |
Popoli et al. | Finite element analysis of mitigation measures for AC interference on buried pipelines | |
Dong et al. | Analysis method for the design of long submarine cables | |
Tsiamitros et al. | Earth return path impedances of underground cables for the two-layer earth case | |
Velásquez et al. | New methodology for design and failure analysis of underground transmission lines | |
CN115577601A (zh) | 穿管敷设下三芯电缆终端头稳态温度场分析方法和系统 | |
Che et al. | Improvement of cable current carrying capacity using COMSOL software | |
Dubyago et al. | Estimation of Insulating Materials Depreciation and Forecasting the Residual cable Resource considering the current core Temperature | |
Xu et al. | Grounding system analysis of 220 kV power cable lines installed underneath a bridge | |
Martins-Britto et al. | Inductive interferences between a 500 kV power line and a pipeline with a complex approximation layout and multilayered soil | |
Aras et al. | An assessment of the methods for calculating ampacity of underground power cables | |
Wang et al. | Steady-state electrical properties and loss reduction measures analysis for 220 kV tunnel cable system | |
Kropotin et al. | Mathematical model of XLPE insulated cable power line with underground installation | |
Szczegielniak et al. | Jabło nski, P | |
Ametani | Four-terminal parameter formulation of solving induced voltages and currents on a pipeline system | |
Wang et al. | Fault current calculation and grounding scattering optimization of 220 kV tunnel cable systems | |
CN115859744A (zh) | 电缆支架的涡流损耗计算方法及系统 | |
Patel et al. | Application of the MoM-SO method for accurate impedance calculation of single-core cables enclosed by a conducting pipe |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230117 |