CN109583112A - 一种基于abaqus有限元软件的架空导线找形方法 - Google Patents
一种基于abaqus有限元软件的架空导线找形方法 Download PDFInfo
- Publication number
- CN109583112A CN109583112A CN201811487109.7A CN201811487109A CN109583112A CN 109583112 A CN109583112 A CN 109583112A CN 201811487109 A CN201811487109 A CN 201811487109A CN 109583112 A CN109583112 A CN 109583112A
- Authority
- CN
- China
- Prior art keywords
- wire
- lead
- calculating
- length
- stress
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000005484 gravity Effects 0.000 claims abstract description 66
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 239000004020 conductor Substances 0.000 claims abstract description 27
- 238000004088 simulation Methods 0.000 claims abstract description 14
- 230000001133 acceleration Effects 0.000 claims description 11
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 description 8
- 238000005452 bending Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于ABAQUS有限元软件的架空导线找形方法,包括以下步骤:S1、计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度即导线原长;S2、计算导线挂点直线距离,将导线挂点直线距离和S1中计算得到的导线原长的大小作对比,从而判定导线的初始形状;S3、按照导线的初始形状建模并预设初始装配应力,设置分析步中的衰减系数后进行有限元分析求取解算结果,获得模拟弧垂和模拟最大应力;S4、计算解算结果的相对误差值并判定建模的正确性,若误差值符合预定阈值条件,结束解算,若误差值不符合预定阈值条件,返回S3校正分析步中的衰减系数后继续执行,直到符合预定阈值条件。本发明的方法操作简单,计算结果稳定,精确度高。
Description
技术领域
本发明属于电力架空导线找形技术领域,具体涉及一种基于ABAQUS有限元软件的架空导线找形方法。
背景技术
对于输电线路安全,分析它的各种受力状况是研究线路安全的基础。输电线作为一种悬索结构,只能受拉而不能受压,而且不承受弯矩,导线的初始构型关系在导线在动力学分析中有着重要地位,确定导线的初始构型是分析导线动力学响应的起点,直接关系到了以后的计算精度。
目前基于有限元软件的找形方法有以下几种:
1、直接迭代法,该方法见于张旺海的“基于ANSYS的架空输电导线找形研究”以及李力的“基于ANSYS的输电导线找形方法研究”;
2、小弹性模量法,该方法见于张旺海的“基于ANSYS的架空输电导线找形研究”;
3、预制模型更新法:此方法见于孔伟的“架空输电线路导线在ANSYS中的找形分析”以及黄新波的“均匀覆冰下的架空线有限元找形分析”。
其中,直接迭代法和小弹性模量迭代法都有着连贯性差的缺点,小弹性模量法还需要设置外置循环,从而不断的更新模型,可操作性较差,预制模型更新法使用有局限性,不能对应更加复杂的情况,而且计算精度不高。
发明内容
有鉴于此,本发明提供一种基于ABAQUS有限元软件的架空导线找形方法,以解决现有技术的不足。
本发明的技术方案是:
一种基于ABAQUS有限元软件的架空导线找形方法,包括以下步骤:
S1、计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度即导线原长;
S2、计算导线挂点直线距离,将导线挂点直线距离和S1中计算得到的导线原长的大小作对比,从而判定导线的初始形状;
S3、按照导线的初始形状建模并预设初始装配应力,设置分析步中的衰减系数后进行有限元分析求取解算结果,获得模拟弧垂和模拟最大应力;
S4、计算解算结果的相对误差值并判定建模的正确性,若误差值符合预定阈值条件,结束解算,若误差值不符合预定阈值条件,返回S3校正分析步中的衰减系数后继续执行,直到符合预定阈值条件。
优选的,S1中计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度的方法包括以下步骤:
S11、利用公式(1)计算所需参数h0
其中,h0为由导线属性确定的悬链线方程参数,无物理意义,σ0为导线重力下水平应力,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S12、利用公式(2)计算曲线左端点的横坐标la
其中,la为导线左端点的横坐标,h0为悬链线方程参数,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S13、利用公式(3)计算曲线右端点横坐标lb
lb=l+la (3)
其中,lb导线右端点横坐标;
S14、利用公式(4)计算重力作用下导线长度L
其中,L为重力作用下导线长度;
S15、利用公式(5)计算重力作用下的导线伸长量ΔL
其中,ΔL为重力引起的导线伸长量,E为弹性模量;
S16、利用公式(6)计算无应力情况下导线长度L0
L0=L-ΔL (6)
其中,L0为无应力情况下导线长度即导线原长,ΔL为重力引起的导线伸长量,L为重力作用下导线长度。
优选的,所述S2中判定导线的初始形状的方法包括以下步骤:
S21、利用公式(7)计算导线挂点直线距离Ls
其中,Ls为导线挂点直线距离,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S22、利用公式(8)计算无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值C
C=L0-Ls (8)
其中,C为无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值,Ls为导线挂点直线距离,L0为无应力情况下导线长度,即导线原长;
S23、判定导线的初始形状,判定方法是:
若C≤0,导线的初始形状为直线,
若C>0,导线的初始形状为v形折线。
优选的,所述S4中计算解算结果的相对误差值并判定建模的正确性的方法包括以下步骤:
S41、利用公式(9)计算导线弧垂处横坐标xmid
其中,xmid为导线弧垂处横坐标,h0为悬链线方程参数,h为导线挂点竖直距离,即高差,l为导线挂点水平距离,即档距;
S42、利用公式(10)计算理论弧垂f
其中,f为理论弧垂,la为导线左端点横坐标,lb导线右端点横坐标;
S43、利用公式(11)计算理论最大应力σmax
其中,σmax为理论最大应力,m取la和lb中绝对值较大的数,h0为悬链线方程参数,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S44、利用公式(12)计算导线弧垂相对误差δ1
其中,δ1为导线弧垂相对误差,f1为模拟弧垂,f为理论弧垂;
S45、利用公式(13)计算导线最大应力相对误差δ2
其中,δ2为导线最大应力相对误差,σmax1为模拟最大应力,σmax为理论最大应力;
S46、设定预定阈值为5%,若δ1≤5%且δ2≤5%,则建模条件正确,符合预定阈值条件,模拟成功;若δ1>5%或δ2>5%,则相对误差值不符合预定阈值条件。
与现有技术相比,本发明提供的一种基于ABAQUS有限元软件的架空导线找形方法操作简单,可操作性好,计算结果稳定,精确度高,实用性好,值得推广。
附图说明
图1为本发明的找形方法的流程示意图;
图2为本发明的导线带塔示意图;
图3为本发明满足导线拐点坐标条件的椭圆。
具体实施方式
本发明提供了一种基于ABAQUS有限元软件的架空导线找形方法,下面结合图1到图3,对本发明进行说明。
如图1所示,本发明提供的一种基于ABAQUS有限元软件的架空导线找形方法,包括以下步骤:
S1、计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度即导线原长;
S2、计算导线挂点直线距离,将导线挂点直线距离和S1中计算得到的导线原长的大小作对比,从而判定导线的初始形状;
S3、按照导线的初始形状建模并预设初始装配应力,设置分析步中的衰减系数后进行有限元分析求取解算结果,获得模拟弧垂和模拟最大应力;
S4、计算解算结果的相对误差值并判定建模的正确性,若误差值符合预定阈值条件,结束解算,若误差值不符合预定阈值条件,返回S3校正分析步中的衰减系数后继续执行,直到符合预定阈值条件。
进一步的,S1中计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度的方法包括以下步骤:
S11、利用公式(1)计算所需参数h0
其中,h0为由导线属性确定的悬链线方程参数,无物理意义,σ0为导线重力下水平应力,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S12、利用公式(2)计算曲线左端点的横坐标la
其中,la为导线左端点的横坐标,h0为悬链线方程参数,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S13、利用公式(3)计算曲线右端点横坐标lb
lb=l+la (3)
其中,lb导线右端点横坐标;
S14、利用公式(4)计算重力作用下导线长度L
其中,L为重力作用下导线长度;
S15、利用公式(5)计算重力作用下的导线伸长量ΔL
其中,ΔL为重力引起的导线伸长量,E为弹性模量;
S16、利用公式(6)计算无应力情况下导线长度L0
L0=L-ΔL (6)
其中,L0为无应力情况下导线长度即导线原长,ΔL为重力引起的导线伸长量,L为重力作用下导线长度。
进一步的,所述S2中判定导线的初始形状的方法包括以下步骤:
S21、利用公式(7)计算导线挂点直线距离Ls
其中,Ls为导线挂点直线距离,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S22、利用公式(8)计算无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值C
C=L0-Ls (8)
其中,C为无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值,Ls为导线挂点直线距离,L0为无应力情况下导线长度,即导线原长;
S23、判定导线的初始形状,判定方法是:
若C≤0,导线的初始形状为直线,
若C>0,导线的初始形状为v形折线。
进一步的,S4中计算解算结果的相对误差值并判定建模的正确性的方法包括以下步骤:
S41、利用公式(9)计算导线弧垂处横坐标xmid
其中,xmid为导线弧垂处横坐标,h0为悬链线方程参数,h为导线挂点竖直距离,即高差,l为导线挂点水平距离,即档距;
S42、利用公式(10)计算理论弧垂f
其中,f为理论弧垂,la为导线左端点横坐标,lb导线右端点横坐标;
S43、利用公式(11)计算理论最大应力σmax
其中,σmax为理论最大应力,m取la和lb中绝对值较大的数,h0为悬链线方程参数,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S44、利用公式(12)计算导线弧垂相对误差δ1
其中,δ1为导线弧垂相对误差,f1为模拟弧垂,f为理论弧垂;
S45、利用公式(13)计算导线最大应力相对误差δ2
其中,δ2为导线最大应力相对误差,σmax1为模拟最大应力,σmax为理论最大应力;
S46、设定预定阈值为5%,若δ1≤5%且δ2≤5%,则建模条件正确,符合预定阈值条件,模拟成功;若δ1>5%或δ2>5%,则相对误差值不符合预定阈值条件。
其中,当导线初始形状为直线时,在ABAQUS中有限元建模仿真的方法包括以下步骤:
(1)计算导线挂点坐标
取导线挂点连线中点为原点,则左挂点为右挂点为
(2)利用公式(14)计算建模所需导线初始装配应力σL
其中,σL为导线初始装配应力;
(3)进入有限元软件进行建模,在软件中以两挂点B1,B2直接建立直线;
(4)模型施加约束:
在初始分析步中,将建好的导线模型两端简支,即固定两端x,y,z三个个方向的位移;
(5)模型施加重力:在随后的分析步中,在模型中施加重力加速度g,
设置衰减系数,使用系统的默认值即可;
(6)划分网格后,提交并进行计算。
其中,当导线初始形状为v形时,在ABAQUS中有限元建模仿真的方法包括以下步骤:
(1)计算导线挂点坐标以及导线拐点坐标;
取导线挂点连线中点为原点,则左挂点为右挂点为
需计算拐点坐标,拐点Q(x1,y1)满足以下公式(15),见图3:
利用公式(16)计算当拐点位于中间位置时的坐标点D(x2,y2)
其中,为导线左挂点坐标,为导线右挂点坐标,Q(x1,y1)为折线拐点一般情况下的坐标,D(x2,y2)为折线拐点中间位置时的坐标。
(2)得到初始装配应力
这里导线初始装配应力为零,但是为了有限元软件运算的收敛这里需要使导线处于紧绷,所以设置σL=1000Mpa,小于应力Mpa级别的1%,对结果影响不大。
(3)进入有限元软件进行建模,在软件中以两挂点B1,B2和拐点D建立折线。
(4)模型施加约束:在初始分析步中,将建好的导线模型两端简支,即固定两端x,y,z三个方向的位移。
(5)模型施加重力:在随后的分析步中,在模型中施加重力加速度g,设置衰减系数,使用默认即可。
(6)划分网格后,提交并进行计算。
其中,建模完成后,在软件的后处理结果中即可获取模拟弧垂f1和模拟最大应力σmax1。
本发明的找形方法首先根据设定的导线最低点水平应力得到无应力情况下导线长度,然后利用有限元软件进行非线性有限元分析完成导线的找形,计算不考虑重力作用时的导线初始长度是本方法的关键之处。
计算原理:
导线为不可压缩且不可承受弯矩的索,导线形状为悬链线,曲线为平面曲线,简化为平面问题。
根据悬链线的相关方程,可以求得导线在受重力作用下的导线弧长,以及导线各点应力的理论值。
根据导线各点应力值计算得到重力引起的导线伸长量,而无应力情况下导线长度是受重力作用下的导线长度与重力引起的导线伸长量之差,利用无应力情况下导线长度,导线挂点距离以及胡克定律可得导线初始装配应力,进而得到导线初始形状。
在有限元软件中根据导线初始形状建模,并施加初始装配应力,然后施加重力进行计算,就可以得到导线在重力下的初始构型。
实施例1:
如图2所示,现有一导线模型,其具体数据为,弹性模量E=55GPa,直径26.81mm,横截面积为426.28mm2,自重为1.18kg/m,现将其安装在档距为300m,无高差的两座塔之间,安装时水平应力为30Mpa,分析并计算导线的形状以及内力。
详细过程:
计算所需参数
其中,h0为由导线属性确定的悬链线方程参数,σ0为导线重力下水平应力,
A为导线的横截面积,p为导线单位长度的质量,g为重力加速度;
计算曲线左端点,右端点
lb=l+la=150
其中,la为导线左端点横坐标,lb为导线右端点横坐标;l为导线档距,h为导线高差。
计算重力作用下导线长度
其中,L为重力作用下导线长度;
计算重力引起的导线伸长量
其中,ΔL为重力引起的导线伸长量,E为弹性模量;
计算导线初始形状和初始装配应力:
其中,L0为无应力情况下导线长度,即导线原长,Ls为导线挂点的直线距离,C为导线长度和导线挂点直线距离之间的差值。
这里C>0,设σL=1000Mpa,导线的初始形状为v形折线,
计算拐点坐标
其中,σL为导线无重力情况下的导线内的应力,D(x2,y2)为折线拐点中间位置时的坐标。
在ABAQUS软件中建模,将模型两端简支,导线上施加预应力重力,然后模拟计算,求取解算结果后得到模拟弧垂f1=10.1899m,模拟最大应力σmax1=30.25Mpa。
理论计算:f=10.1885m,σmax=30.276Mpa。
计算相对误差值:
其中,δ1为导线弧垂相对误差,δ2为导线最大应力相对误差,f1为模拟弧垂,σmax1为模拟最大应力,f为理论弧垂,σmax为理论导线最大应力。
误差满足阈值条件,模拟成功。
实施例2:
现有一导线模型,其具体数据为,弹性模量E=55GPa,直径26.81mm,横截面积为426.28mm2,自重为1.18kg/m,现将其安装在档距为300m,高差为50m的两座塔之间,安装时水平应力为55.1Mpa,分析并计算导线的形状以及内力。
详细过程:
计算所需参数
其中,h0为由导线属性确定的悬链线方程参数,σ0为导线重力下水平应力,A为导线的横截面积,p为导线单位长度的质量,g为重力加速度;
计算曲线左端点,右端点
lb=l+la=-186.6710
其中,la为导线左端点横坐标,lb为导线右端点横坐标,l为导线档距,h为导线高差。
计算重力作用下导线长度
其中,L为重力作用下导线长度;
计算重力引起的导线伸长量
其中,ΔL为由重力引起的导线伸长量,E为弹性模量,
计算导线初始构型和初始装配应力:
其中,L0为导线无重力情况下的长度,Ls为导线挂点的直线距离,C为导线长度和导线挂点直线距离之间的差值。
这里C<0,因此
导线的初始形状为直线;
其中σL为导线无重力情况下的导线内的应力。
在ABAQUS软件中建模,将模型两端简支,导线上施加预应力,然后模拟计算求取解算结果后得到模拟弧垂f1=5.619m,模拟最大应力σmax1=56.689Mpa。
理论计算:f=5.6176m,σmax=56.689Mpa。
计算相对误差值:
其中,δ1为导线弧垂相对误差,δ2为导线最大应力相对误差,f1为模拟弧垂,σmax1为模拟最大应力,f为理论弧垂,σmax为理论导线最大应力。
误差满足,模拟成功。
实施例3:
现有一导线模型,其具体数据为,弹性模量E=55GPa,直径26.81mm,横截面积为426.28mm2,自重为1.18kg/m,现将其安装在档距为300m,无高差的两座塔之间,安装时水平应力为30Mpa。分析并计算导线的形状以及内力,本实施例3与实施例1相同,但更改了拐点坐标。
详细的计算过程:
计算所需参数
其中,h0为由导线属性确定的悬链线方程参数,σ0为导线重力下水平应力,A为导线的横截面积,p为导线单位长度的质量,g为重力加速度;
计算曲线左端点,右端点:
lb=l+la=150
其中,la为导线左端点横坐标,lb为导线右端点横坐标,l为导线档距,h为导线高差。
计算导线现长
其中,L为导线重力作用下长度。
计算重力引起的导线伸长量:
其中,ΔL为重力引起的导线伸长量,E为弹性模量。
计算导线初始构型和初始装配应力:
其中,L0为导线无重力情况下的长度,Ls为导线挂点的直线距离,C为导线长度和导线挂点直线距离之间的差值。
这里C>0,设σL=1000Mpa,导线的初始形状为v形折线,
计算过程中,为了验证拐点的选取可以取任意点,如图3所示的满足导线拐点坐标条件的椭圆,这里设置拐点为1/3处,令x1=100,带入椭圆方程得到y1=-10.0649,其中σL为导线无重力情况下导线内的应力,(x1,y1)为折线拐点为1/4处时的坐标,其中f为导线和挂点连线的最大值,即弧垂;
在ABAQUS软件中建模,将模型两端简支,导线上施加预应力,然后模拟计算得到结果,模拟弧垂f1=10.186m,模拟最大应力σmax1=30.28Mpa。
理论值:f=10.1885m,σmax=30.276Mpa。
计算相对误差值:
其中,δ1为导线弧垂相对误差,δ2为导线最大应力相对误差,f1为模拟弧垂,σmax1为模拟最大应力,f为理论弧垂,σmax为理论导线最大应力。
误差满足,模拟成功。
与现有技术相比,本发明提供的一种基于ABAQUS有限元软件的架空导线找形方法,建模简单,减少了计算步骤,减轻了操作负担,计算结果稳定,能够收敛,适用范围广,实用性好,值得推广。
以上公开的仅为本发明的较佳的具体实施例,但是,本发明实施例并非局限于此,任何本领域技术人员能思之的变化都应落入本发明的保护范围。
Claims (4)
1.一种基于ABAQUS有限元软件的架空导线找形方法,其特征在于,包括以下步骤:
S1、计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度即导线原长;
S2、计算导线挂点直线距离,将导线挂点直线距离和S1中计算得到的导线原长的大小作对比,从而判定导线的初始形状;
S3、按照导线的初始形状建模并预设初始装配应力,设置分析步中的衰减系数后进行有限元分析求取解算结果,获得模拟弧垂和模拟最大应力;
S4、计算解算结果的相对误差值并判定建模的正确性,若误差值符合预定阈值条件,结束解算,若误差值不符合预定阈值条件,返回S3校正分析步中的衰减系数后继续执行,直到符合预定阈值条件。
2.如权利要求1所述的一种基于ABAQUS有限元软件的架空导线找形方法,其特征在于,所述S1中计算重力作用下的导线长度和导线伸长量,求取无应力情况下导线长度的方法包括以下步骤:
S11、利用公式(1)计算所需参数h0
其中,h0为由导线属性确定的悬链线方程参数,无物理意义,σ0为导线重力下水平应力,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S12、利用公式(2)计算曲线左端点的横坐标la
其中,la为导线左端点的横坐标,h0为悬链线方程参数,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S13、利用公式(3)计算曲线右端点横坐标lb
lb=l+la (3)
其中,lb导线右端点横坐标;
S14、利用公式(4)计算重力作用下导线长度L
其中,L为重力作用下导线长度;
S15、利用公式(5)计算重力作用下的导线伸长量ΔL
其中,ΔL为重力引起的导线伸长量,E为弹性模量;
S16、利用公式(6)计算无应力情况下导线长度L0
L0=L-ΔL (6)
其中,L0为无应力情况下导线长度即导线原长,ΔL为重力引起的导线伸长量,L为重力作用下导线长度。
3.如权利要求1所述的一种基于ABAQUS有限元软件的架空导线找形方法,其特征在于,所述S2中判定导线的初始形状的方法包括以下步骤:
S21、利用公式(7)计算导线挂点直线距离Ls
其中,Ls为导线挂点直线距离,l为导线挂点水平距离,即档距,h为导线挂点竖直距离,即高差;
S22、利用公式(8)计算无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值C
C=L0-Ls (8)
其中,C为无应力情况下导线长度L0和导线挂点直线距离Ls之间的差值,Ls为导线挂点直线距离,L0为无应力情况下导线长度,即导线原长;
S23、判定导线的初始形状,判定方法是:
若C≤0,导线的初始形状为直线,
若C>0,导线的初始形状为v形折线。
4.如权利要求1所述的一种基于ABAQUS有限元软件的架空导线找形方法,其特征在于,所述S4中计算解算结果的相对误差值并判定建模的正确性的方法包括以下步骤:
S41、利用公式(9)计算导线弧垂处横坐标xmid
其中,xmid为导线弧垂处横坐标,h0为悬链线方程参数,h为导线挂点竖直距离,即高差,l为导线挂点水平距离,即档距;
S42、利用公式(10)计算理论弧垂f
其中,f为理论弧垂,la为导线左端点横坐标,lb导线右端点横坐标;
S43、利用公式(11)计算理论最大应力σmax
其中,σmax为理论最大应力,m取la和lb中绝对值较大的数,h0为悬链线方程参数,A为导线的横截面积,p为导线单位长度质量,g为重力加速度;
S44、利用公式(12)计算导线弧垂相对误差δ1
其中,δ1为导线弧垂相对误差,f1为模拟弧垂,f为理论弧垂;
S45、利用公式(13)计算导线最大应力相对误差δ2
其中,δ2为导线最大应力相对误差,σmax1为模拟最大应力,σmax为理论最大应力;
S46、设定预定阈值为5%,若δ1≤5%且δ2≤5%,则建模条件正确,符合预定阈值条件,模拟成功;若δ1>5%或δ2>5%,则相对误差值不符合预定阈值条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811487109.7A CN109583112B (zh) | 2018-12-06 | 2018-12-06 | 一种基于abaqus有限元软件的架空导线找形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811487109.7A CN109583112B (zh) | 2018-12-06 | 2018-12-06 | 一种基于abaqus有限元软件的架空导线找形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109583112A true CN109583112A (zh) | 2019-04-05 |
CN109583112B CN109583112B (zh) | 2022-03-29 |
Family
ID=65926096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811487109.7A Expired - Fee Related CN109583112B (zh) | 2018-12-06 | 2018-12-06 | 一种基于abaqus有限元软件的架空导线找形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109583112B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363677A (zh) * | 2019-03-13 | 2019-10-22 | 国网山西省电力公司经济技术研究院 | 一种覆冰输电线路找形方法 |
CN114818420A (zh) * | 2022-04-18 | 2022-07-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008077414A (ja) * | 2006-09-21 | 2008-04-03 | Olympus Corp | 撚り線ワイヤの変形解析シミュレーションモデル、撚り線ワイヤの加工方法、及び撚り線ワイヤの変形解析シミュレーション方法 |
CN104166776A (zh) * | 2014-08-29 | 2014-11-26 | 南京理工大学 | 一种基于ansys的输电线路导线找形方法 |
CN104318130A (zh) * | 2014-11-25 | 2015-01-28 | 华北电力大学(保定) | 一种拉线塔拉线的非线性计算方法 |
CN104504284A (zh) * | 2015-01-04 | 2015-04-08 | 西安电子科技大学 | 一种基于悬链线单元的松弛索网找形方法 |
CN105426597A (zh) * | 2015-11-09 | 2016-03-23 | 中国电力科学研究院 | 一种架空输电线路大截面导线的找形确定方法 |
CN107194120A (zh) * | 2017-06-20 | 2017-09-22 | 南京理工大学 | 一种基于有限质点法的覆冰输电线路找形方法 |
CN108710763A (zh) * | 2018-05-22 | 2018-10-26 | 国网江西省电力有限公司经济技术研究院 | 220kV输电线路塔线耦联体系覆冰仿真方法 |
-
2018
- 2018-12-06 CN CN201811487109.7A patent/CN109583112B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008077414A (ja) * | 2006-09-21 | 2008-04-03 | Olympus Corp | 撚り線ワイヤの変形解析シミュレーションモデル、撚り線ワイヤの加工方法、及び撚り線ワイヤの変形解析シミュレーション方法 |
CN104166776A (zh) * | 2014-08-29 | 2014-11-26 | 南京理工大学 | 一种基于ansys的输电线路导线找形方法 |
CN104318130A (zh) * | 2014-11-25 | 2015-01-28 | 华北电力大学(保定) | 一种拉线塔拉线的非线性计算方法 |
CN104504284A (zh) * | 2015-01-04 | 2015-04-08 | 西安电子科技大学 | 一种基于悬链线单元的松弛索网找形方法 |
CN105426597A (zh) * | 2015-11-09 | 2016-03-23 | 中国电力科学研究院 | 一种架空输电线路大截面导线的找形确定方法 |
CN107194120A (zh) * | 2017-06-20 | 2017-09-22 | 南京理工大学 | 一种基于有限质点法的覆冰输电线路找形方法 |
CN108710763A (zh) * | 2018-05-22 | 2018-10-26 | 国网江西省电力有限公司经济技术研究院 | 220kV输电线路塔线耦联体系覆冰仿真方法 |
Non-Patent Citations (3)
Title |
---|
YANG WEIGUO: "Second Shape Finding Analysis of Membrane Structures", 《TSINGHUA SCIENCE AND TECHNOLOGY》 * |
林晖: "一种弹性边界单层平面索网快速找形方法", 《钢结构》 * |
贾玉琢: "覆冰输电架空导线初始构形研究", 《水电能源科学》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363677A (zh) * | 2019-03-13 | 2019-10-22 | 国网山西省电力公司经济技术研究院 | 一种覆冰输电线路找形方法 |
CN114818420A (zh) * | 2022-04-18 | 2022-07-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
CN114818420B (zh) * | 2022-04-18 | 2024-03-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109583112B (zh) | 2022-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Song et al. | Nonlinear modelling of high-speed catenary based on analytical expressions of cable and truss elements | |
CN109583112B (zh) | 一种基于abaqus有限元软件的架空导线找形方法 | |
CN108959742B (zh) | 大跨越输电塔线体系气动弹性模型设计方法 | |
CN105335569B (zh) | 特高压输电线路多联v型绝缘子串力学特性仿真模拟方法 | |
CN113177343A (zh) | 一种输电线路弧垂计算方法、设备及计算机介质 | |
CN106096105A (zh) | 输电线路风偏瞬态响应计算方法 | |
CN102721875A (zh) | 一种基于输电导线表面等效电荷曲线积分的三维工频电场的测定方法 | |
CN111783201A (zh) | 一种三跨自锚式悬索桥动力特性的快速分析方法 | |
CN109858117B (zh) | 一种单跨覆冰输电线路脱冰最大跳跃高度的算法 | |
JP6366477B2 (ja) | 鉄塔装柱設計支援システム | |
CN110363677B (zh) | 一种覆冰输电线路找形方法 | |
CN111506944A (zh) | 考虑塔线耦合效应的杆塔风荷载脉动折减系数计算方法 | |
CN110968935A (zh) | 一种任意高差覆冰输电线路脱冰最大跳跃高度的算法 | |
CN104850746A (zh) | 一种基于四阶龙格-库塔和模拟退火的等值盐密预测方法 | |
WO2020245892A1 (ja) | 設備状態解析装置、設備状態解析方法、及びプログラム | |
CN108199792B (zh) | 一种wcdma基站电磁辐射预测方法 | |
CN116933602A (zh) | 输电塔线耦联体系覆冰舞动响应的数值模拟方法及系统 | |
Wu et al. | Characteristics of multi-span transmission lines following ice-shedding | |
CN111241731B (zh) | 一种依据水平张力的输电导线找形方法 | |
CN110645939A (zh) | 一种两条导线之间的相对垂直距离的动态预测系统及方法 | |
CN116189831A (zh) | 一种铝合金低温断裂延伸率的测量方法及系统 | |
CN110224335A (zh) | 一种考虑蠕变量的架空导线架设方法和系统 | |
CN111125626B (zh) | 一种基于s型函数随机子空间识别的模型定阶方法 | |
CN111737863B (zh) | 基于能量法的多跨架空输电线覆冰找形计算方法 | |
CN111737862B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220329 |
|
CF01 | Termination of patent right due to non-payment of annual fee |