CN108416097A - 地线与预绞丝接触端口接触电阻的有限元热分析等效方法 - Google Patents
地线与预绞丝接触端口接触电阻的有限元热分析等效方法 Download PDFInfo
- Publication number
- CN108416097A CN108416097A CN201810102993.1A CN201810102993A CN108416097A CN 108416097 A CN108416097 A CN 108416097A CN 201810102993 A CN201810102993 A CN 201810102993A CN 108416097 A CN108416097 A CN 108416097A
- Authority
- CN
- China
- Prior art keywords
- contact
- ground wire
- preformed armor
- armor rods
- contact resistance
- 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
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]
Abstract
本发明公开了一种地线与预绞丝接触端口接触电阻的有限元热分析等效方法,首先通过实验获取架空地线与预绞丝接触端口的接触电阻的阻值;计算得到该接触电阻对应的接触面积;根据接触电阻和接触面积,计算出热流密度;建立地线与预绞丝的三维模型,并在接触端口的接触面积处加载计算得到的热流密度;最后,设计大电流温升实验,验证有限元仿真结果。本发明提供的方法,避免了在有限元仿真中对于复杂接触界面的建模要求,具有普遍推广使用的意义。
Description
技术领域
本发明涉及电网及过热分析研究技术领域,特别涉及一种地线与预绞丝接触端口的接触电阻有限元热分析等效方法。
背景技术
架空地线因过热导致断线的事故时有发生,致使大面积电网停电,严重影响电网的安全稳定运行。现有研究表明,当工频短路电流流经地线,在地线与预绞丝接触端口处,电流因集肤效应往预绞丝处扩散,地线与预绞丝接触端口的接触电阻是导致地线过热断线的一个重要因素。因此,对于地线与预绞丝接触端口过热分析研究,考虑接触电阻的存在是一个不可缺失的环节。
目前对于过热分析研究,有限元仿真计算是主要的手段之一。在有限元仿真中,对于复杂的接触界面进行几何建模难以实现,因此对于接触界面接触电阻进行等效设置是一个主要的解决途径,也是需要解决的一大难点。基于此,针对有限元仿真计算中接触电阻等处理问题,提出一种基于地线与预绞丝接触端口的接触电阻有限元热分析等效方法。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种地线与预绞丝接触端口接触电阻的有限元热分析等效方法,避免了对复杂接触界面的几何建模,直接在接触界面加载接触电阻对应的热流密度,同时不影响地线与预绞丝接触端口的传热情况,减少了仿真工作所需时间。
本发明的目的通过以下的技术方案实现:一种地线与预绞丝接触端口接触电阻的有限元热分析等效方法,包括以下步骤:
1、通过实验获取架空地线与预绞丝接触端口的接触电阻的阻值;
2、计算得到该接触电阻对应的接触面积;
3、根据实验得到的接触电阻和计算得到的接触面积,计算出热流密度;
4、建立地线与预绞丝的三维模型,并在接触端口的接触面积处加载计算得到的热流密度;
5、设计大电流温升实验,验证有限元仿真结果。
优选的,步骤1中,通过实验获取架空地线与预绞丝接触端口的接触电阻阻值,所采用的实验测量方法为交流压降法或者热电法;交流压降法可以通过欧姆定律直接得到接触电阻,热电法则是通过电阻的热效应间接求取。
优选的,步骤2中,先获得该接触电阻的轴向等效接触长度d和截面的环形接触长度l这两个参数,接触面积S=dl。
具体的,上述两个尺寸参数通过计算公式或者有限元仿真计算。
优选的,步骤2中接触面积通过Morgan的径向应力推导直接得到。
具体的,接触面积其中,d为接触电阻的轴向等效接触长度;Fr为单位长度地线表面径向所受的力;f为钢的屈服强度;Fr的可通过单位长度地线表面所受的总应力T与预绞丝螺旋角α2以及单根钢丝的直径dwire计算得到:
地线表面所受的总应力T和螺旋角α可由以下两式求解:
其中,T’为单位长度地线总应力;α1和m1分别为地线表面层的螺旋角和线股节径比,m2为预绞丝线股的节径比。
优选的,步骤3中,根据接触电阻以及接触面积计算接触电阻产生的热流密度,如下式:
式中,q为热流密度,I为加载的电流,Rc为接触电阻,S为接触面积。
优选的,步骤4中,收集地线和预绞丝的尺寸参数,利用有限元仿真软件ANSYS或者comsol建立地线与预绞丝的三维模型;添加激励时在预绞丝与地线接触部位加载上一步骤计算得到的接触电阻产生的热流密度。
优选的,步骤5中,利用大电流发生器设计地线与预绞丝的温升实验平台,加载一定电流得到地线与预绞丝接触端口的温升数据,结合仿真得到的数据进行验证分析。
本发明与现有技术相比,具有如下优点和有益效果:
本发明通过计算得到接触电阻的热流密度,加载至地线与预绞丝接触端口,以准确分析地线与预绞丝接触端口的温升情况。避免了对复杂接触界面的几何建模,直接在接触界面加载接触电阻对应的热流密度,同时不影响地线与预绞丝接触端口的传热情况,减少了仿真工作所需时间。
附图说明
图1为实施例基于地线与预绞丝接触端口的接触电阻有限元热分析等效方法流程图。
图2为实施例地线与预绞丝接触端口轴向电流密度分布示意图。
图3为实施例地线与预绞丝接触端口结构图。
图4为实施例地线与预绞丝接触端口径向示意图。
图5为实施例温升实验系统原理图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例1
一种地线与预绞丝接触端口的接触电阻有限元热分析等效方法,步骤流程图如图1所示,主要包括以下步骤:
步骤1,通过实验获取架空地线与预绞丝接触端口的接触电阻的阻值;
步骤2,在地线与预绞丝接触端口处,接触电阻产生的热效应因电流在此处扩散所致,利用有限元电磁仿真或者解析计算得到该接触电阻对应的接触面积;
步骤3,根据实验得到的接触电阻和计算得到的接触面积,计算出热流密度;
步骤4,利用有限元仿真软件ANSYS或者comsol建立地线与预绞丝的三维模型,并在接触端口的接触面积处加载计算得到的热流密度;
步骤5,设计大电流温升实验,验证有限元仿真结果。
所述步骤1,架空地线与预绞丝接触端口的接触电阻为交流接触电阻,难以用四线法直接测量。通过实验测量获取地线与预绞丝接触端口的接触电阻,测量方法有交流电压降法、热电法等。以交流电压降法为例,加载一定大小的交流电流I,测量得到地线与预绞丝接触端口两端的电压U以及功率因数并通过式(1)计算出接触电阻。
所述步骤2,一种计算接触面积的方法是:首先计算得到接触电阻的轴向等效接触长度和截面的环形接触长度。
对于接触电阻的轴向等效接触长度,如图2中的d所示。等效接触长度d可以通过有限元电磁仿真计算或者电磁解析计算得到。
地线与预绞丝接触端口的结构图以及径向接触情况如图3和图4所示。接触部位均出现在相邻两层线股凸起部位。径向接触长度l(即截面的环形接触长度)如图4中的红线所示。接触长度l的求解利用结构力学的计算公式或者有限元仿真计算获得。
接触面积由式(2)可以求解。
S=dl (2)
其次,也可以采用绞线径向应力推导接触面积的方法,如式(3)所示。
其中,d为等效接触长度,Fr为单位长度地线表面径向所受的力,单位N;f为钢的屈服强度,单位Mpa。Fr可通过单位长度地线表面所受的总应力T与预绞丝螺旋角α2以及单根钢丝的直径dwire计算得到。
地线表面所受的总应力T和螺旋角α可由以下两式求解。
其中,T’为单位长度地线总应力,单位Mpa,一般取为最小破拉断力的20%;α1和m1分别为地线表面层的螺旋角和线股节径比,m2为预绞丝线股的节径比。
所述步骤3,根据步骤1和步骤2得到的接触电阻和接触面积,通过式(6)求得接触电阻的热流密度。
所述步骤4,利用ANSYS或者comsol有限元仿真软件建立地线与预绞丝接触端口的三维几何模型。添加激励时,在接触端口处加载接触电阻等效的热流密度。最终基于此仿真模型进行传热分析。
所述步骤5,设计地线与预绞丝接触端口的温升实验,通过加载一定的交流电流,获取该部位离散点的温度分布,实验装置如图5所示。
综上所述,一种可靠的接触电阻有限元热分析的等效方法,避免在有限元仿真中对于复杂接触界面的建模要求,减少了仿真工作量,且不影响研究对象的传热情况,具有普遍推广使用的意义。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (9)
1.地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,包括以下步骤:
S1、通过实验获取架空地线与预绞丝接触端口的接触电阻的阻值;
S2、计算得到该接触电阻对应的接触面积;
S3、根据实验得到的接触电阻和计算得到的接触面积,计算出热流密度;
S4、建立地线与预绞丝的三维模型,并在接触端口的接触面积处加载计算得到的热流密度;
S5、设计大电流温升实验,验证有限元仿真结果。
2.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S1中,通过实验获取架空地线与预绞丝接触端口的接触电阻阻值,所采用的实验测量方法为交流压降法或者热电法;交流压降法可以通过欧姆定律直接得到接触电阻,热电法则是通过电阻的热效应间接求取。
3.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S2中,先获得该接触电阻的轴向等效接触长度d和截面的环形接触长度l这两个参数,接触面积S=dl。
4.根据权利要求3所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,其中两个尺寸参数通过计算公式或者有限元仿真计算。
5.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S2中接触面积通过Morgan的径向应力推导直接得到。
6.根据权利要求1或5所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,接触面积其中,d为接触电阻的轴向等效接触长度;Fr为单位长度地线表面径向所受的力;f为钢的屈服强度;Fr的可通过单位长度地线表面所受的总应力T与预绞丝螺旋角α2以及单根钢丝的直径dwire计算得到:
地线表面所受的总应力T和螺旋角α可由以下两式求解:
其中,T’为单位长度地线总应力;α1和m1分别为地线表面层的螺旋角和线股节径比,m2为预绞丝线股的节径比。
7.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S3中,根据接触电阻以及接触面积计算接触电阻产生的热流密度,如下式:
式中,q为热流密度,I为加载的电流,Rc为接触电阻,S为接触面积。
8.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S4中,收集地线和预绞丝的尺寸参数,利用有限元仿真软件ANSYS或者comsol建立地线与预绞丝的三维模型;添加激励时在预绞丝与地线接触部位加载上一步骤计算得到的接触电阻产生的热流密度。
9.根据权利要求1所述的地线与预绞丝接触端口接触电阻的有限元热分析等效方法,其特征在于,步骤S5中,利用大电流发生器设计地线与预绞丝的温升实验平台,加载一定电流得到地线与预绞丝接触端口的温升数据,结合仿真得到的数据进行验证分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810102993.1A CN108416097A (zh) | 2018-02-01 | 2018-02-01 | 地线与预绞丝接触端口接触电阻的有限元热分析等效方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810102993.1A CN108416097A (zh) | 2018-02-01 | 2018-02-01 | 地线与预绞丝接触端口接触电阻的有限元热分析等效方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108416097A true CN108416097A (zh) | 2018-08-17 |
Family
ID=63127663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810102993.1A Pending CN108416097A (zh) | 2018-02-01 | 2018-02-01 | 地线与预绞丝接触端口接触电阻的有限元热分析等效方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108416097A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109346984A (zh) * | 2018-11-12 | 2019-02-15 | 广东电网有限责任公司 | 一种减小地线与预绞丝接触端口接触电阻的方法和装置 |
CN109460616A (zh) * | 2018-11-12 | 2019-03-12 | 广东电网有限责任公司 | 接触电阻计算方法及装置 |
CN109508496A (zh) * | 2018-11-12 | 2019-03-22 | 广东电网有限责任公司 | 接触模型建立方法及装置 |
CN109683020A (zh) * | 2018-12-24 | 2019-04-26 | 广东电网有限责任公司 | 计算架空地线与预绞丝接触端口接触电阻的装置和方法 |
CN111707872A (zh) * | 2018-11-12 | 2020-09-25 | 广东电网有限责任公司 | 接触电阻测量方法及装置 |
CN112269060A (zh) * | 2020-09-09 | 2021-01-26 | 国联汽车动力电池研究院有限责任公司 | 一种接触电阻及接触电导率辨识方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11275723A (ja) * | 1998-03-23 | 1999-10-08 | Iimuro Dengyo Kk | スリーブプロテクター |
CN104731898A (zh) * | 2015-03-18 | 2015-06-24 | 华南理工大学 | 一种10kV三芯电缆有限元热分析方法 |
CN106934099A (zh) * | 2017-02-10 | 2017-07-07 | 华南理工大学 | 一种确定耐张线夹载流量与三维温度场分布的方法 |
CN107066719A (zh) * | 2017-04-06 | 2017-08-18 | 国网浙江省电力公司舟山供电公司 | 光纤复合海底电缆接地故障时温度分布仿真分析方法 |
CN107526872A (zh) * | 2017-07-24 | 2017-12-29 | 国网江苏省电力公司南京供电公司 | 一种500kV超高压电缆的热应力及形变量的计算方法 |
-
2018
- 2018-02-01 CN CN201810102993.1A patent/CN108416097A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11275723A (ja) * | 1998-03-23 | 1999-10-08 | Iimuro Dengyo Kk | スリーブプロテクター |
CN104731898A (zh) * | 2015-03-18 | 2015-06-24 | 华南理工大学 | 一种10kV三芯电缆有限元热分析方法 |
CN106934099A (zh) * | 2017-02-10 | 2017-07-07 | 华南理工大学 | 一种确定耐张线夹载流量与三维温度场分布的方法 |
CN107066719A (zh) * | 2017-04-06 | 2017-08-18 | 国网浙江省电力公司舟山供电公司 | 光纤复合海底电缆接地故障时温度分布仿真分析方法 |
CN107526872A (zh) * | 2017-07-24 | 2017-12-29 | 国网江苏省电力公司南京供电公司 | 一种500kV超高压电缆的热应力及形变量的计算方法 |
Non-Patent Citations (4)
Title |
---|
HANSHENG CAI 等: "Lightning performance of EHV and UHV overhead transmission Lines in China southern power grid", 《INTERNATIONAL CONFERENCE ON LIGHTNING PROTECTION》 * |
刘超 等: "考虑接触电阻下的电缆接头热点温度", 《高电压技术》 * |
郭德明 等: "缠绕预绞丝的地线悬垂串温度分布实验分析", 《广东电力》 * |
阮飞: "多因子作用下架空输电线路动态增容计算模型", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109346984A (zh) * | 2018-11-12 | 2019-02-15 | 广东电网有限责任公司 | 一种减小地线与预绞丝接触端口接触电阻的方法和装置 |
CN109460616A (zh) * | 2018-11-12 | 2019-03-12 | 广东电网有限责任公司 | 接触电阻计算方法及装置 |
CN109508496A (zh) * | 2018-11-12 | 2019-03-22 | 广东电网有限责任公司 | 接触模型建立方法及装置 |
CN109460616B (zh) * | 2018-11-12 | 2020-09-25 | 广东电网有限责任公司 | 接触电阻计算方法及装置 |
CN111707872A (zh) * | 2018-11-12 | 2020-09-25 | 广东电网有限责任公司 | 接触电阻测量方法及装置 |
CN109508496B (zh) * | 2018-11-12 | 2021-01-22 | 广东电网有限责任公司 | 接触模型建立方法及装置 |
CN111707872B (zh) * | 2018-11-12 | 2023-02-07 | 广东电网有限责任公司 | 接触电阻测量方法及装置 |
CN109683020A (zh) * | 2018-12-24 | 2019-04-26 | 广东电网有限责任公司 | 计算架空地线与预绞丝接触端口接触电阻的装置和方法 |
CN112269060A (zh) * | 2020-09-09 | 2021-01-26 | 国联汽车动力电池研究院有限责任公司 | 一种接触电阻及接触电导率辨识方法及系统 |
CN112269060B (zh) * | 2020-09-09 | 2024-01-02 | 国联汽车动力电池研究院有限责任公司 | 一种接触电阻及接触电导率辨识方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108416097A (zh) | 地线与预绞丝接触端口接触电阻的有限元热分析等效方法 | |
CN108334695B (zh) | 一种基于地线与预绞丝接触电阻的有限元设置方法 | |
CN103728539B (zh) | 基于分布式光纤测温的电缆电气故障仿真分析方法 | |
KR100934448B1 (ko) | 병렬로 연결된 복수의 초전도 선재에 흐르는 전류의비접촉식 측정방법 | |
CN104217061A (zh) | 低压配电柜的温度场仿真设计方法 | |
Chatzipetros et al. | Impact of proximity effects on sheath losses in trefoil cable arrangements | |
Holyk et al. | Power cable rating calculations-a historical perspective [history] | |
Martin et al. | High-frequency litz “air-gap” windings for high-power density electrical machines | |
Morgan | Effects of alternating and direct current, power frequency, temperature, and tension on the electrical parameters of ACSR conductors | |
Mitchell | Steady state analysis of non-uniform current distributions in cable-in-conduit conductors and comparison with experimental data | |
CN106960095A (zh) | 一种确定导线蠕变率的方法及系统 | |
Salmi et al. | Modeling heat transfer from quench protection heaters to superconducting cables in Nb3Sn magnets | |
Napieralski et al. | Nonuniform distribution of conductivity resulting from the stress exerted on a stranded cable during the manufacturing process | |
Rolando et al. | Performance assessment and optimization of the ITER toroidal field coil joints | |
Khodak et al. | Analysis of TF insert coil | |
CN105489519B (zh) | 高电阻结构的原位温度检测装置、芯片和失效检测方法 | |
Kajitani et al. | Analytical study of degradation of CIC conductor performance due to strand bending and buckling | |
Soldati et al. | AC loss investigation in MgB 2 multifilamentary wires: a numerical study. | |
Breschi et al. | Performance review of the joints for the ITER poloidal field coils | |
Breschi et al. | Electrodynamic losses of the ITER PF joints during the dynamic plasma scenario | |
CN110083905A (zh) | 一种隔离开关动态增容的热点温度分析方法 | |
Soppe et al. | Thermal modelling of fittings for high temperature low sag conductors | |
Chenying et al. | Simulation and analysis of power cable harmonic loss and temperature rise | |
Konstantopoulou et al. | Electro-magnetic performance of MgB 2 cables for the high current transmission lines at CERN | |
CN102394165B (zh) | 一种电阻板的绕制方法 |
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: 20180817 |