CN103455686A - 架空输电塔-线耦合体系有限元模型的建模方法 - Google Patents
架空输电塔-线耦合体系有限元模型的建模方法 Download PDFInfo
- Publication number
- CN103455686A CN103455686A CN2013104211878A CN201310421187A CN103455686A CN 103455686 A CN103455686 A CN 103455686A CN 2013104211878 A CN2013104211878 A CN 2013104211878A CN 201310421187 A CN201310421187 A CN 201310421187A CN 103455686 A CN103455686 A CN 103455686A
- Authority
- CN
- China
- Prior art keywords
- finite element
- element model
- tower
- power transmission
- overhead power
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种架空输电塔-线耦合体系有限元模型的建模方法,属于架空输电线路抗冰能力评估技术领域。本发明利用计算机,通过程序,先建立架空输电杆塔、导线及地线有限元模型库,后组建输电塔-线耦合体系有限元模型基本结构,再用几何非线性有限元找形分析方法建立自重力作用下准确的架空输电塔-线耦合体系有限元模型,从而实现架空输电塔-线耦合体系有限元模型的快速高效建立。本发明具有方法简单、可操作性强;建模工作量少、效率高;模型准确真实靠性、精确度高,适用性强等优点。本发明可广泛应用于架空输电塔-线耦合体系有限元模型的建立,特别适用于大面积多条线路架空输电塔-线耦合体系的建立,为抗冰能力的评估提供准确而可靠的基础。
Description
技术领域
本发明属于架空输电线路抗冰能力评估技术领域,具体涉及架空输电线路的塔-线耦合体系有限元模型的建模方法。
背景技术
电力***的架空输电线路长期暴露在自然界中,容易受覆冰等自然灾害的影响。我国架空输电线路的覆冰灾害较为严重,大面积覆冰灾害事故时有发生,给国民经济及人民生活造成了巨大的灾难和损失。随着我国经济建设的飞速发展和人民生活水平的提高,防止架空输电线路覆冰灾害的发生,保障电力***的安全稳定运行已经是亟待解决的问题。架空输电线路塔-线耦合体系的抗冰能力评估,能有效指导覆冰在线监测设备的安装和线路抗冰改造工作,从而提高电力***架空输电线路抵御覆冰灾害天气的能力,确保电力***的安全稳定运行。建立架空输电塔-线耦合体系有限元模型是架空输电线路塔-线耦合体系抗冰能力评估的基础,因此研究架空输电塔-线耦合体系有限元模型的建模方法具有极其重要的意义。
现有架空输电塔-线耦合体系有限元模型的建模方法,如2011年11月23日公布的公布号为CN102254061A的“输电塔线塔耦合体系的有限元建模及受力方法”专利,公开的建模方法是:利用计算机,通过程序,先通过对输电塔整体结构件的三维CAD模型主动简化并提取框体特征结构,后利用BEAM188梁单元建立特征结构的ANSYS有限元模型库,再按照输电塔结构的实际形状,通过变换工作平面,由单面向多面逐渐过渡,动态的由单面结构逐渐“拼接”形成多面结构,同时通过解析与数值的综合方法对到导线进行参数化建模,并与输电塔模型进行组合,最后确定输电塔线塔耦合体系有限元模型。该建模方法的主要缺点是:①建立杆塔有限元模型时,需要建立三维CAD模型并导入到ANSYS中,人工简化提取框体特征结构并需变换工作平面,导致建模的速度较慢并且耗时较长;②未考虑绝缘子串对架空输电塔-线耦合体系的影响,并且建立的是简化的杆塔模型,并非实际截面形状,因此建模的准确度较低;③建立的导、地线有限元模型未经初始形态找形分析,导致建立的塔线塔有限元模型与实际模型有差异,模型的可靠度和精确度较低,从而导致用此模型进行受力分析进行抗冰能力评估得到结果不准确;④建立杆塔模型时,需要人工进行简化提取框体特征结构,需要有CAD与ANSYS软件的接口;用BEAM188梁单元建立了特征结构的有限元模型库,每建立一个框体结构都需要调用相应特征结构的模型库;同时还要变换工作平面,单面向多面逐渐过渡,动态的由单面结构逐渐“拼接”形成多面结构,导致该建模方法比较复杂,技术难度大,可操作性较低,适用性差,不适合用此法建立的模型来进行大面积多条架空输电线路的抗冰能力评估。
发明内容
本发明针的目的是针对现有架空输电塔-线耦合体系有限元模型的建模方法的不足,提供一种架空输电塔-线耦合体系有限元模型的建模方法,具有方法简单、可操作性强;建模工作量少、效率高;模型准确真实靠性、精确度高,适用性强等优点。能为架空输电塔-线耦合体系抗冰能力的评估,特别是针对大面积多条架空输电塔-线耦合体系抗冰能力的评估,提供准确而可靠的模型基础。
本发明的基本原理:
①“杆塔定位模块”的基本原理
当在迪科尔坐标系中,该模型的第一基杆塔的塔身顶部中心点为坐标原点时(本发明杆塔有限元模型建立时也是以该点为坐标原点,为了简化,本发明的杆塔坐标均指该点坐标),坐标为(xi,yi,zi)=(0,0,0),其中:x轴为杆塔横担方向、y轴为杆塔高度方向、z轴为线路延伸方向。并以该坐标原点作为全局坐标系来定位后面所有杆塔的位置。又设第i基杆塔相对于第i-1基杆塔的档距、高差和转角分别为li、hi和θi。第i基杆塔的坐标为(xi,yi,zi)则:
第i基杆塔的x轴坐标xi为:
式中li、θi分别表示第i基杆塔相对于第i-1基杆塔的档距和转角。
第i基杆塔的y轴坐标yi为:
式中hi表示第i基杆塔相对于第i-1基杆塔的高差。
第i基杆塔的z轴坐标zi为:
式中li、θi分别表示第i基杆塔相对于第i-1基杆塔的档距和转角。
i即杆塔相对于第一基杆塔的转角θ为:
式中θi分别表示第i基杆塔相对于第i-1基杆塔的转角。
其中li×sinθi和li×cosθi分别表示第i基杆塔相对上一基杆塔在x轴和z轴方向上的位移。因此基于第i基杆塔的坐标为(xi,yi,zi)及i即杆塔相对于第1基杆塔的距离,即可得到第i基杆塔在以第1基杆塔为坐标原点的全局坐标系中的位置坐标。
②“初始形态分析”的基本原理
根据导、地线在正常工作状态下的受力变形均在弹性形变范围内,不涉及到材料非线性问题,属于几何非线性范畴的大变形小应变问题,因此,本发明采用几何非线性有限元找形方法对导、地线进行找形分析。其基本原理是:将索膜结构离散为若干单元,然后针对索膜结构的大位移小应变特点,应用几何非线性理论,建立以节点位移为基本未知量的非线性有限元方程组,从一个假定的初始几何形状开始,通过设定索膜结构的初应力,进行迭代求解,得到相应的几何形状。在ANSYS非线性分析中,通过找形结果的具体表现为:找形结束后施加自重力进行求解,得到导线最低点水平应力在误差范围内和设置的预应力相当,而同时只有较小的位移,理论上进位移为零时则误差也为零。
实现本发明目的的技术方案是:一种架空输电塔-线耦合体系有限元模型的建模方法,利用计算机,通过程序,先建立架空输电杆塔、导线及地线有限元模型库,后组建架空输电塔-线耦合体系有限元模型基本结构,再利用几何非线性有限元找形分析方法,建立自重力作用下准确的架空输电塔-线耦合体系有限元模型,从而实现架空输电塔-线耦合体系有限元模型的快速高效建立。所述方法的具体步骤如下:
(1)建立杆塔、导线和地线有限元模型库
首先,利用现有的ANSYS软件(根据有限元理论,以高质量高效率完成复杂结构的分析计算的大型通用软件)中现有的BEAM188梁单元来模拟架空输电杆塔(简称杆塔)构件。即先对BEAM188梁单元进行设置:分别设置钢材Q235、Q345和Q420的杨氏模量、泊松比、密度、截面形状、屈服强度;后根据杆塔的实际形状和实际尺寸,利用现有的ANSYS软件建立实际杆塔有限元模型;再将已建立的实际杆塔有限元模型保存成为APDL文件,并存在ANSYS目录下,形成架空输电线路常用的不同种杆塔有限元模型库,以备用。
然后,用现有的ANSYS软件中现有的LINK10单元模拟导线和地线。即:输入架空输电线路常用的架空输电导、地线(简称导、地线)的基本参数:等效截面积、杨氏模量、密度、泊松比及初始应变,并分别保存成为各导、地线型号的APDL文件,从而形成架空输电线路常用的不同型号的导、地线有限元模型库,以备用。
(2)组建架空输电塔-线耦合体系有限元模型基本结构
第(1)步完成后,通过程序,在第(1)步中建立的架空输电线路的杆塔、导线和地线有限元模型库中,根据实际架空输电线路调取杆塔、导线和地线的有限元模型,组建架空输电塔-线耦合体系的有限元模型基本结构,其具体步骤如下:
①输入模型实际参数
利用现有的ANSYS软件,先在编写的程序中分别输入实际架空输电塔-线耦合体系中的杆塔、绝缘子串、导线和地线的实际参数。其中,杆塔的实际参数为:模型中杆塔总数、杆塔型号、杆塔呼称高度、塔间档距、塔间高差;导线和地线的实际参数为:导、地线型号及其***数;绝缘子串的实际参数为:绝缘子串长度、绝缘子串质量和绝缘子串直径。
②检查实际参数完整准确性
第(2)-①步完成后,检查输入实际参数的完整性与准确性:当输入的实际参数不完整或不正确时,则返回(2)-①步,继续输入未完的实际参数或更改输入实际参数,以保证输入实际参数的完整性与准确性;当输入的实际参数完整准确时,则可进入下一步调用第(1)步建立的有限元模型库开始建模。
③建立杆塔有限元模型
第(2)-②步完成后,将输有实际参数的程序输入至现有的ANSYS软件中,在第(1)步建立的杆塔有限元模型库中,先调出实际架空输电线路的各基杆塔有限元模型;后将“杆塔定位模块”程序输入到现有的ANSYS中,定位各基杆塔。从而实现了整个架空输电塔-线耦合体系有限元模型中杆塔模型的建立及定位。
所述“杆塔定位模块”:以第一基杆塔的塔身顶点中间位置为坐标原点,以第i基杆塔与i-1基杆塔之间的档距、高差和转角为基本参数来实现第i基杆塔坐标位置的确定。
④建立绝缘子串有限元模型
第(2)-③步完成后,先将输有实际参数的程序输入至现有的ANSYS软件中,后调用“绝缘子串生成模块”程序。所述“绝缘子串生成模块”:由第i基杆塔转角θi来判断此杆塔悬挂绝缘子串的类型。当转角θi为0时,第i基杆塔为直线塔,则在此基杆塔绝缘子悬挂点处生成悬垂绝缘子串有限元模型;当转角θi不为0时,第i基杆塔为耐张塔(特殊情况:杆塔为直线转角塔,则将转角θi的值改成一个小数值,如0.001),则在此基杆塔绝缘子悬挂点处生成耐张绝缘子串有限元模型。从而实现了整个架空输电塔-线耦合体系有限元模型中绝缘子串有限元模型的建立。
⑤建立导、地线有限元模型
第(2)-④步完成后,先将“导、地线生成模块”程序输入现有的ANSYS软件中,在第(1)步建立的导线和地线有限元模型库中,调出实际型号的导、地线有限元模型;后在导、地线弦线位置处生成导、地线有限元模型,建立架空输电塔-线耦合体系有限元模型基本结构。
所述“导、地线生成模块”采用直接生成有限元的方式,并同时综合考虑计算量和计算精度来剖分单元,将每个单元长度定为0.9-1.1米,根据两悬挂点间的直线距离取整数即为单元数目,从而便于通过循环命令流直接生成有限元模型,同时也能保证具有足够的计算精度。
(3)建立准确架空输电塔-线耦合体系有限元模型
第(2)步完成后,对第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,即用非线性有限元分析方法,确定导、地线在自重力作用下的初始形态。对架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型,其具体步骤如下:
①设置初始应变、边界条件及求解精度
通过程序,调用第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构(其导线和地线的有限元模型用编写的程序自动生成直线形状,避免了采用悬链线或抛物线建模带来的诸多困难),为初始形态找形分析做好准备;后设置模型初始应变,根据应力和应变关系得到的真实初始应变,通过实常数设置LINK10单元的属性;再设置边界条件,设置每基杆塔塔腿及边界塔的导、地线悬挂点处均为全自由度约束;然后设置模型求解精度ERR0,求解精度是根据实际的精度要求,在综合考虑计算量的情况下给出的。对第(2)步建立的架空输电塔-线耦合体系有限元模型设置初始应变、边界条件及求解精度,以确保架空输电塔-线耦合体系有限元模型的精确性和准确性。
②进行非线性求解
第(3)-①步完成后,先设置非线性求解选项,后进行非线性求解。其中,非线性选项的设置包括大变形(Large Deformation Effects)、应力刚化(Stress Stiffening Effects),同时设置相应的子步数。求解中可采用Newton-Raphson方法进行迭代收敛控制,以帮助得到较好的收敛结果。Newton-Raphson法可有效克服逐步荷载增量法中因误差积累导致求解结果发散的问题,在某个容限范围内使求解达到平衡收敛,采用Newton-Raphson法计算输电导、地线找形时需不断更新体系位形,进行多次非线性静力计算,直至导、地线的最大位移矢量差接近于零且方向一致时,找形计算结束,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型。
③计算最低点水平应力和求解误差
第(3)-②步完成后,利用现有的ANSYS软件,通过程序,先提取中间点应力,后求解其水平分量,再计算求解误差ERR1。
④建立准确架空输电塔-线耦合体系有限元模型
第(3)-③步完成后,判断是否满足非线性求解终止条件:当ERR1≥ERR0,则返回(3)-②步继续进行非线性求解迭代,直到满足求解终止条件为止;当ERR1<ERR0时,则找形分析结束,即可得到自重力作用下准确的导、地线有限元模型,从而建立了自重力作用下的准确架空输电塔-线耦合体系有限元模型。
本发明采用以上的技术方案后主要有以下效果:
(1)本发明利用现有的ANSYS软件,结合输电线路杆塔结构、材料特性等,建立ANSYS有限元模型库,通过编制程序,调用模型库,避免了ANSYS与CAD软件的接口问题,能准确、快速高效、省时省力的建立架空输电塔-线耦合体系的有限元模型,大大减少建立模型的工作量,这为实现基于有限元的输电线路塔-线耦合体系抗冰能力评估提供了基础;
(2)本发明开发的架空输电塔-线耦合体系有限元模型的建模方法,建模时考虑了绝缘子串的影响,用实际截面形状和实际尺寸在ANSYS软件中建立杆塔有限元模型,未对其进行简化,更加真实的模拟了实际运行的输电杆塔;建立了导、地线基本模型之后,对其进行了初始形态找形分析,建立的模型与实际吻合程度高,更加真实可靠。用此方法建立的模型进行分析计算,结果更精确,适用性更强,更加有利于大面积多条架空输电塔-线耦合体系的建立并对其抗冰能力进行评估,从而确保架空输电线路的安全可靠运行;
(3)本发明开发的架空输电塔-线耦合体系有限元模型的建模方法,运用了几何非线性有限元方法编制程序,对导、地线模型基本结构进行初始形态找形分析,建立的导、地线模型相较直接建立悬链线导、地线结构,更加接近实际运行状况,有更高的精确度和可靠度,有更高的建模效率,从而实现了架空输电塔-线耦合体系有限元模型高效快速的建立;
(4)本发明开发的架空输电塔-线耦合体系有限元模型的建模方法,简单易行,可操作性强。此方法是在编制的程序中输入架空输电塔-线耦合体系的实际参数(输电杆塔数量和型号、档距、高差、转角,导线地线型号等),自动调用杆塔、导线及地线的有限元模型库、初始形态找形分析后,建立架空输电塔-线耦合体系有限元模型。人员进行操作时,只需修改各段架空输电线路的实际参数,输入程序至ANSYS软件就可以自动建立对应的架空输电塔-线耦合体系有限元模型,方法简单易行,可操作性强。
本发明可广泛应用于架空输电塔-线耦合体系有限元模型的建立,特别适用于大面积多条架空输电塔-线耦合体系的建立,并为抗冰能力的评估提供准确而可靠的基础。
附图说明
图1为本发明的架空输电塔-线耦合体系有限元模型建模的***流程框图;
图2为图1中建立绝缘子串有限元模型中“绝缘子串生成模块”子流程图;
图3为图1中建立准确架空输电塔-线耦合体系有限元模型中“初始形态找形分析”子流程图;
图4为实施例的架空输电塔-线耦合体系有限元模型基本结构;
图5为图4中的模型经过初始形态找形分析后的准确架空输电塔-线耦合体系有限元模型;
图6为图5中的黑框中架空输电塔-线耦合体系有限元模型局部放大图;
图中:1为杆塔有限元模型,2为绝缘子串有限元模型,3为地线有限元模型,4为导线有限元模型。
具体实施方式
下面结合具体实施方式,进一步说明本发明。
实施例:
一种架空输电塔-线耦合体系有限元模型的建模方法,以某500kV输电线路的一耐张段为例,建立其架空输电塔-线耦合体系有限元模型的具体步骤如下:
(1)建立杆塔、导线和地线有限元模型库
首先,利用现有的ANSYS软件中现有的BEAM188梁单元来模拟架空输电杆塔构件,即先对BEAM188梁单元进行设置:所使用的Q235、Q345和Q420三种钢材具有相同的杨氏模量、泊松比及密度,分别为206000N/mm2、0.3、及7850kg/m3,屈服强度分别为235MPa、345MPa、420MPa,截面形状为L40×3到L200×24不同的角钢;后根据杆塔的实际形状和实际尺寸,利用现有的ANSYS软件建立实际杆塔有限元模型;再将已建立的实际杆塔有限元模型保存成为APDL文件,并存在ANSYS目录下,形成架空输电线路常用的不同种杆塔有限元模型库。
然后,用现有的ANSYS软件中现有的LINK10单元模拟导线和地线。500kV架空输电线路常用导线的型号:LGJ-400/50、LGJ-400/35、LGJ-630/45、A3/S3A-465/60、A3/S1A-732/92、JLHA1/G1A-400/50、JLHA2/G1A-630/45等;500kV架空输电线路常用地线的型号:GJ-80、GJ-100、GJ-120、GJ-150、GJ-210、GJ-240、LBGJ-100/20AC、LBGJ-100-40AC、LBGJ-120-40AC、LBGJ-150-40AC、JLB20A-100、JLB20A-120、JLB14A-150等。分别输入架空输电线路各型号的导、地线基本参数:等效截面积、杨氏模量、密度、泊松比及初始应变,并分别保存成为各导、地线型号的APDL文件,从而形成架空输电线路常用的不同型号的导、地线有限元模型库。导线LGJ-400/50和地线JLB20A-100特性参数如表1所示:
表1导线LGJ-400/50和地线JLB20A-100特性参数
导线型号 | LGJ-400/50 | 导线型号 | JLB20A-100 |
等效截面积(mm2) | 451.55 | 等效截面积(mm2) | 100.88 |
密度(kg/m3) | 3346 | 密度(kg/m3) | 6681 |
杨氏模量(N/mm2) | 69000 | 杨氏模量(N/mm2) | 147200 |
破坏应力(N/mm2) | 259.6 | 破坏应力(N/mm2) | 1206 |
初始应变 | 1e-6 | 初始应变 | 1e-5 |
泊松比 | 0.3 | 泊松比 | 0.3 |
(2)组建架空输电塔-线耦合体系有限元模型基本结构
第(1)步完成后,通过程序,调用在第(1)步中建立的架空输电线路的杆塔、导线和地线有限元模型库,限元模型库中,根据实际架空输电线路调取杆塔、导线和地线的有限元模型,组建架空输电塔-线耦合体系的有限元模型基本结构,其具体步骤如下:①输入模型实际参数
利用现有的ANSYS软件,先在编写的程序中分别输入实际架空塔-线耦合体系中的杆塔、绝缘子串、导线和地线的实际参数。其中,杆塔的实际参数为:模型中杆塔总数7、各基杆塔型号、杆塔呼称高度、塔间档距、塔间高差;导线和地线的实际参数为:导线型号LGJ-400/50、***数为4,地线型号JLB20A-100、***数为1;绝缘子串的实际参数为:绝缘子串长度、绝缘子串质量和绝缘子串直径。
表2实际架空输电线路某耐张区间杆塔明细表
杆塔运行号 | 杆塔型式 | 呼称高度(m) | 水平档距(m) | 高差(m) | 转角(°) |
308# | JG11 | 27 | — | — | — |
309# | ZB61A | 27 | 123 | 15.8 | — |
310# | ZB62A | 43.5 | 240 | -2.4 | — |
311# | ZB62A | 42 | 356 | 34.9 | — |
312# | ZB61A | 31.5 | 225 | 16.8 | — |
313# | ZB62A | 31.5 | 351 | -30.9 | — |
314# | JG22 | 27 | 692 | -6.2 | 左21°03' |
②检查实际参数完整准确性
第(2)-①步完成后,检查输入实际参数的完整性与准确性:经检查输入的实际参数完整准确,则可进入下一步调用有限元模型库开始建模。
③建立杆塔有限元模型
第(2)-②步完成后,将输有实际参数的程序输入至现有的ANSYS软件中,在第(1)步建立的杆塔有限元模型库中,先调出实际架空输电线路的JG11、ZB61A、ZB62A、ZB62A、ZB61A、ZB62A、JG22杆塔有限元模型;后将“杆塔定位模块”程序输入到现有的ANSYS中,定位各基杆塔。从而实现了整个架空输电塔-线耦合体系有限元模型中杆塔模型的建立及定位。
④建立绝缘子串有限元模型
第(2)-③步完成后,先将输有实际参数的程序输入至现有的ANSYS软件中,后调用“绝缘子串生成模块”程序。经判断,为第1和第7基杆塔生成耐张绝缘子串有限元模型,为第2-6基杆塔生成悬垂绝缘子串有限元模型。从而实现了整个架空输电塔-线耦合体系有限元模型中绝缘子串有限元模型的建立。
⑤建立导、地线有限元模型
第(2)-④步完成后,先将“导、地线生成模块”程序输入现有的ANSYS软件中,在第(1)步建立的导线和地线有限元模型库中,调出实际型号为LGJ-400/50的导线和实际型号为JLB20A-100的地线有限元模型;后在导、地线弦线位置处生成导、地线有限元模型,建立架空输电塔-线耦合体系有限元模型基本结构,如图4所示。
(3)建立准确架空输电塔-线耦合体系有限元模型
第(2)步完成后,对第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,即用非线性有限元分析方法,确定导、地线在自重力作用下的初始形态。对架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型,其具体步骤如下:
①设置初始应变、边界条件及求解精度
通过程序,调用第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构,为初始形态找形分析做好准备;后设置导线初始应变为1e-6,地线初始应变为1e-5,通过实常数设置LINK10单元横截面积;再设置边界条件,设置每基杆塔塔腿及边界塔的导、地线悬挂点处均为全自由度约束;然后设置模型求解精度ERR0为0.001。
②进行非线性求解
第(3)-①步完成后,先设置非线性求解选项,后进行非线性求解。其中,非线性选项的设置包括大变形、应力刚化,设置相应的子步数为20。采用Newton-Raphson方法,经多次非线性迭代计算,直至导、地线的最大位移矢量差接近于零且方向一致时,找形计算结束,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型。
③计算最低点水平应力和求解误差
第(3)-②步完成后,利用现有的ANSYS软件,通过程序,先提取中间点应力,后求解其水平分量,再计算求解误差ERR1。
④建立准确架空输电塔-线耦合体系有限元模型
第(3)-③步完成后,判断是否满足非线性求解终止条件:当ERR1≥ERR0,则返回(3)-②步继续进行非线性求解迭代,直到满足求解终止条件为止;当ERR1<ERR0时,则找形分析结束,即可得到自重力作用下准确的导、地线有限元模型,从而建立了自重力作用下的准确架空输电塔-线耦合体系有限元模型,如图5所示。
Claims (2)
1.一种架空输电塔-线耦合体系有限元模型的建模方法,利用计算机,通过程序,建立架空输电塔-线耦合体系有限元模型,其特征在于所述方法的具体步骤如下:
(1)建立杆塔、导线和地线有限元模型库
首先,利用现有的ANSYS软件中现有的BEAM188梁单元来模拟架空输电杆塔,即杆塔构件,即先对BEAM188梁单元进行设置:分别设置钢材Q235、Q345和Q420的杨氏模量、泊松比、密度、截面形状、屈服强度;后根据杆塔的实际形状和实际尺寸,利用现有的ANSYS软件建立实际杆塔有限元模型;再将已建立的实际杆塔有限元模型保存成为APDL文件,并存在ANSYS目录下,形成架空输电线路常用的不同种杆塔有限元模型库,以备用;
然后,用现有的ANSYS软件中现有的LINK10单元模拟导线和地线,即:输入架空输电线路常用的架空输电导、地线,即导、地线的基本参数:等效截面积、杨氏模量、密度、泊松比及初始应变,并分别保存成为各导、地线型号的APDL文件,从而形成架空输电线路常用的不同型号的导、地线有限元模型库,以备用;
(2)组建架空输电塔-线耦合体系有限元模型基本结构
第(1)步完成后,通过程序,在第(1)步中建立的架空输电线路的杆塔、导线和地线有限元模型库中,根据实际架空输电线路调取杆塔、导线和地线的有限元模型,组建架空输电塔-线耦合体系的有限元模型基本结构,其具体步骤如下:
①输入模型实际参数
利用现有的ANSYS软件,先在编写的程序中分别输入实际架空输电塔-线耦合体系中的杆塔、绝缘子串、导线和地线的实际参数,其中,杆塔的实际参数为:模型中杆塔总数、杆塔型号、杆塔呼称高度、塔间档距、塔间高差;导线和地线的实际参数为:导、地线型号及其***数;绝缘子串的实际参数为:绝缘子串长度、绝缘子串质量和绝缘子串直径;
②检查实际参数完整准确性
第(2)-①步完成后,检查输入实际参数的完整性与准确性:当输入的实际参数不完整或不正确时,则返回(2)-①步,继续输入未完的实际参数或更改输入实际参数,以保证输入实际参数的完整性与准确性;当输入的实际参数完整准确时,则可进入下一步调用第(1)步建立的有限元模型库开始建模;
③建立杆塔有限元模型
第(2)-②步完成后,将输有实际参数的程序输入至现有的ANSYS软件中,在第(1)步建立的杆塔有限元模型库中,先调出实际架空输电线路的各基杆塔有限元模型;后将“杆塔定位模块”程序输入到现有的ANSYS中,定位各基杆塔,从而实现了整个架空输电塔-线耦合体系有限元模型中杆塔模型的建立及定位;
④建立绝缘子串有限元模型
第(2)-③步完成后,先将输有实际参数的程序输入至现有的ANSYS软件中,后调用“绝缘子串生成模块”程序,所述“绝缘子串生成模块”:由第i基杆塔转角θi来判断此杆塔悬挂绝缘子串的类型,当转角θi为0时,第i基杆塔为直线塔,则在此基杆塔绝缘子悬挂点处生成悬垂绝缘子串有限元模型;当转角θi不为0时,第i基杆塔为耐张塔,特殊情况:杆塔为直线转角塔,则将转角θi的值改成一个小数值,如0.001,则在此基杆塔绝缘子悬挂点处生成耐张绝缘子串有限元模型,从而实现了整个架空输电塔-线耦合体系有限元模型中绝缘子串有限元模型的建立;
⑤建立导、地线有限元模型
第(2)-④步完成后,先将“导、地线生成模块”程序输入现有的ANSYS软件中,在第(1)步建立的导线和地线有限元模型库中,调出实际型号的导、地线有限元模型;后在导、地线弦线位置处生成导、地线有限元模型,建立架空输电塔-线耦合体系有限元模型基本结构;
(3)建立准确架空输电塔-线耦合体系有限元模型
第(2)步完成后,对第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,即用非线性有限元分析方法,确定导、地线在自重力作用下的初始形态,对架空输电塔-线耦合体系有限元模型基本结构进行初始形态找形分析,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型,其具体步骤如下:
①设置初始应变、边界条件及求解精度
通过程序,调用第(2)步建立的架空输电塔-线耦合体系有限元模型基本结构,为初始形态找形分析做好准备;后设置模型初始应变,根据应力和应变关系得到的真实初始应变,通过实常数设置LINK10单元的属性;再设置边界条件,设置每基杆塔塔腿及边界塔的导、地线悬挂点处均为全自由度约束;然后设置模型求解精度ERR0,求解精度是根据实际的精度要求,在综合考虑计算量的情况下给出的,对第(2)步建立的架空输电塔-线耦合体系有限元模型设置初始应变、边界条件及求解精度,以确保架空输电塔-线耦合体系有限元模型的精确性和准确性;
②进行非线性求解
第(3)-①步完成后,先设置非线性求解选项,后进行非线性求解,其中,非线性选项的设置包括大变形、应力刚化,同时设置相应的子步数,求解中可采用Newton-Raphson方法进行迭代收敛控制,以帮助得到较好的收敛结果,Newton-Raphson法可有效克服逐步荷载增量法中因误差积累导致求解结果发散的问题,在某个容限范围内使求解达到平衡收敛,采用Newton-Raphson法计算输电导、地线找形时需不断更新体系位形,进行多次非线性静力计算,直至导、地线的最大位移矢量差接近于零且方向一致时,找形计算结束,建立自重力作用下的准确架空输电塔-线耦合体系有限元模型;
③计算最低点水平应力和求解误差
第(3)-②步完成后,利用现有的ANSYS软件,通过“求解”程序,先提取中间点应力,后求解其水平分量,再计算求解误差ERR1;
④建立准确架空输电塔-线耦合体系有限元模型
第(3)-③步完成后,判断是否满足非线性求解终止条件:当ERR1≥ERR0,则返回(3)-②步继续进行非线性求解迭代,直到满足求解终止条件为止;当ERR1<ERR0时,则找形分析结束,即可得到自重力作用下准确的导、地线有限元模型,从而建立了自重力作用下的准确架空输电塔-线耦合体系有限元模型。
2.按照权利要求1所述的架空输电塔-线耦合体系有限元模型的建模方法,其特征在于:
所述“杆塔定位模块”:以第一基杆塔的塔身顶点中间位置为坐标原点,以第i基杆塔与i-1基杆塔之间的档距、高差和转角为基本参数来实现第i基杆塔坐标位置的确定;
所述“导、地线生成模块”:采用直接生成有限元的方式,并同时综合考虑计算量和计算精度来剖分单元,将每个单元长度定为0.9-1.1米,根据两悬挂点间的直线距离取整数即为单元数目,从而便于通过循环命令流直接生成有限元模型,同时也能保证具有足够的计算精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013104211878A CN103455686A (zh) | 2013-09-17 | 2013-09-17 | 架空输电塔-线耦合体系有限元模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013104211878A CN103455686A (zh) | 2013-09-17 | 2013-09-17 | 架空输电塔-线耦合体系有限元模型的建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103455686A true CN103455686A (zh) | 2013-12-18 |
Family
ID=49738042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013104211878A Pending CN103455686A (zh) | 2013-09-17 | 2013-09-17 | 架空输电塔-线耦合体系有限元模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103455686A (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955867A (zh) * | 2014-04-25 | 2014-07-30 | 常州东瑞电力软件有限公司 | 一种配电网接线图自动生成装置 |
CN104166776A (zh) * | 2014-08-29 | 2014-11-26 | 南京理工大学 | 一种基于ansys的输电线路导线找形方法 |
CN104298808A (zh) * | 2014-08-26 | 2015-01-21 | 国家电网公司 | 输电铁塔非线性柔性构件的应力计算方法 |
CN104331566A (zh) * | 2014-11-12 | 2015-02-04 | 国家电网公司 | 500kV交流直线转角塔电场仿真*** |
CN104504167A (zh) * | 2014-08-29 | 2015-04-08 | 广西电网公司电力科学研究院 | 一种输电杆塔结构静态力学分析***及其仿真方法 |
CN105335569A (zh) * | 2015-11-16 | 2016-02-17 | 国家电网公司 | 特高压输电线路多联v型绝缘子串力学特性仿真模拟方法 |
CN105354382A (zh) * | 2015-11-05 | 2016-02-24 | 中国电力科学研究院 | 一种建立输电线路塔线体系有限元模型的方法 |
CN105740549A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 220kV输电线路塔线耦联体系风振仿真方法 |
CN105740548A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 一种随机风载荷下输电线路风振计算方法 |
CN107092766A (zh) * | 2017-06-27 | 2017-08-25 | 国网江苏省电力公司电力科学研究院 | 一种基于有限元的导线风偏绝缘校核方法 |
CN107133379A (zh) * | 2017-03-29 | 2017-09-05 | 河北省电力勘测设计研究院 | 特高压塔架柱建模***及方法 |
CN107220396A (zh) * | 2014-08-26 | 2017-09-29 | 江苏省电力公司南通供电公司 | 解决奇异性矩阵奇异性问题的输电铁塔杆件应力计算方法 |
CN107784142A (zh) * | 2016-08-31 | 2018-03-09 | 株洲中车时代电气股份有限公司 | 半导体功率组件的热仿真方法 |
CN108197353A (zh) * | 2017-12-17 | 2018-06-22 | 内蒙航天动力机械测试所 | 一种基于ansys的apdl语言的固体火箭发动机工装设计方法 |
CN109726498A (zh) * | 2019-01-09 | 2019-05-07 | 华南理工大学 | 一种结合仿真精确计算架空地线预绞丝接触端口受力方法 |
CN109918772A (zh) * | 2019-03-05 | 2019-06-21 | 中国科学院近代物理研究所 | 线圈建模方法、装置、电子设备及存储介质 |
CN110334459A (zh) * | 2019-07-11 | 2019-10-15 | 广东电网有限责任公司 | 一种输电塔塔线体系快速精细化建模***及方法 |
CN111125945A (zh) * | 2019-11-28 | 2020-05-08 | 广东电网有限责任公司 | 一种输电线路极限抗风能力评估方法和装置 |
CN111159925A (zh) * | 2019-11-22 | 2020-05-15 | 国网江苏省电力有限公司电力科学研究院 | 瓷绝缘子运行荷载条件下机械性能评估的有限元分析方法、***和存储介质 |
CN112036061A (zh) * | 2020-08-07 | 2020-12-04 | 国网浙江省电力有限公司电力科学研究院 | 全方位长短腿塔线体系有限元建模及动力响应分析方法 |
CN112269397A (zh) * | 2020-10-30 | 2021-01-26 | 国家电网有限公司 | 一种基于高精度定位和视觉跟踪的无人机巡检方法 |
CN113553747A (zh) * | 2021-08-27 | 2021-10-26 | 贵州电网有限责任公司 | 一种覆冰输电塔力学特性分析及状态评估方法 |
CN113642222A (zh) * | 2021-08-27 | 2021-11-12 | 贵州电网有限责任公司 | 综合考虑微地形与荷载的基于有限元分析的杆塔选型方法 |
CN114818420A (zh) * | 2022-04-18 | 2022-07-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409439A (zh) * | 2008-11-19 | 2009-04-15 | 中国电力工程顾问集团公司 | 一种确定重冰区输电线路杆塔导地线布置的方法 |
CN102663215A (zh) * | 2012-05-14 | 2012-09-12 | 重庆大学 | 架空输电线路塔-线耦合体系抗冰能力的评估方法 |
-
2013
- 2013-09-17 CN CN2013104211878A patent/CN103455686A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409439A (zh) * | 2008-11-19 | 2009-04-15 | 中国电力工程顾问集团公司 | 一种确定重冰区输电线路杆塔导地线布置的方法 |
CN102663215A (zh) * | 2012-05-14 | 2012-09-12 | 重庆大学 | 架空输电线路塔-线耦合体系抗冰能力的评估方法 |
Non-Patent Citations (3)
Title |
---|
YANG FENGLI等: "Numerical Simulation on the HV Transmission Tower-Line System under Ice Shedding", 《TRANSMISSION & DISTRIBUTION CONFERENCE & EXPOSITION: ASIA AND PACIFIC》 * |
张旺海 等: "基于ANSYS的架空输电导线找形研究", 《电力建设》 * |
毛峰: "基于有限元方法的不同覆冰状态输电塔线体系力学特性研究", 《万方学术论文数据库》 * |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955867A (zh) * | 2014-04-25 | 2014-07-30 | 常州东瑞电力软件有限公司 | 一种配电网接线图自动生成装置 |
CN104298808A (zh) * | 2014-08-26 | 2015-01-21 | 国家电网公司 | 输电铁塔非线性柔性构件的应力计算方法 |
CN105005660A (zh) * | 2014-08-26 | 2015-10-28 | 国家电网公司 | 贴近实际运行状态的非线性柔性构件的应力计算方法 |
CN105069189A (zh) * | 2014-08-26 | 2015-11-18 | 国家电网公司 | 贴近实际输电铁塔构造的非线性柔性构件的应力计算方法 |
CN104298808B (zh) * | 2014-08-26 | 2017-09-22 | 国家电网公司 | 输电铁塔非线性柔性构件的应力计算方法 |
CN107220396A (zh) * | 2014-08-26 | 2017-09-29 | 江苏省电力公司南通供电公司 | 解决奇异性矩阵奇异性问题的输电铁塔杆件应力计算方法 |
CN104504167B (zh) * | 2014-08-29 | 2018-11-06 | 广西电网公司电力科学研究院 | 一种输电杆塔结构静态力学分析***及其仿真方法 |
CN104166776A (zh) * | 2014-08-29 | 2014-11-26 | 南京理工大学 | 一种基于ansys的输电线路导线找形方法 |
CN104504167A (zh) * | 2014-08-29 | 2015-04-08 | 广西电网公司电力科学研究院 | 一种输电杆塔结构静态力学分析***及其仿真方法 |
CN104166776B (zh) * | 2014-08-29 | 2017-04-05 | 南京理工大学 | 一种基于ansys的输电线路导线找形方法 |
CN104331566A (zh) * | 2014-11-12 | 2015-02-04 | 国家电网公司 | 500kV交流直线转角塔电场仿真*** |
CN104331566B (zh) * | 2014-11-12 | 2017-09-15 | 国家电网公司 | 500kV交流直线转角塔电场仿真*** |
CN105354382A (zh) * | 2015-11-05 | 2016-02-24 | 中国电力科学研究院 | 一种建立输电线路塔线体系有限元模型的方法 |
CN105354382B (zh) * | 2015-11-05 | 2020-06-26 | 中国电力科学研究院 | 一种建立输电线路塔线体系有限元模型的方法 |
CN105335569B (zh) * | 2015-11-16 | 2019-03-26 | 国家电网公司 | 特高压输电线路多联v型绝缘子串力学特性仿真模拟方法 |
CN105335569A (zh) * | 2015-11-16 | 2016-02-17 | 国家电网公司 | 特高压输电线路多联v型绝缘子串力学特性仿真模拟方法 |
CN105740549A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 220kV输电线路塔线耦联体系风振仿真方法 |
CN105740548A (zh) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | 一种随机风载荷下输电线路风振计算方法 |
CN105740548B (zh) * | 2016-02-01 | 2019-05-24 | 西安交通大学 | 一种随机风载荷下输电线路风振计算方法 |
CN107784142B (zh) * | 2016-08-31 | 2021-10-19 | 株洲中车时代半导体有限公司 | 半导体功率组件的热仿真方法 |
CN107784142A (zh) * | 2016-08-31 | 2018-03-09 | 株洲中车时代电气股份有限公司 | 半导体功率组件的热仿真方法 |
CN107133379A (zh) * | 2017-03-29 | 2017-09-05 | 河北省电力勘测设计研究院 | 特高压塔架柱建模***及方法 |
CN107092766A (zh) * | 2017-06-27 | 2017-08-25 | 国网江苏省电力公司电力科学研究院 | 一种基于有限元的导线风偏绝缘校核方法 |
CN108197353A (zh) * | 2017-12-17 | 2018-06-22 | 内蒙航天动力机械测试所 | 一种基于ansys的apdl语言的固体火箭发动机工装设计方法 |
CN109726498A (zh) * | 2019-01-09 | 2019-05-07 | 华南理工大学 | 一种结合仿真精确计算架空地线预绞丝接触端口受力方法 |
CN109726498B (zh) * | 2019-01-09 | 2022-10-25 | 华南理工大学 | 一种结合仿真精确计算架空地线预绞丝接触端口受力方法 |
CN109918772A (zh) * | 2019-03-05 | 2019-06-21 | 中国科学院近代物理研究所 | 线圈建模方法、装置、电子设备及存储介质 |
CN109918772B (zh) * | 2019-03-05 | 2023-04-18 | 中国科学院近代物理研究所 | 线圈建模方法、装置、电子设备及存储介质 |
CN110334459A (zh) * | 2019-07-11 | 2019-10-15 | 广东电网有限责任公司 | 一种输电塔塔线体系快速精细化建模***及方法 |
CN111159925B (zh) * | 2019-11-22 | 2022-08-23 | 国网江苏省电力有限公司电力科学研究院 | 瓷绝缘子机械性能评估的有限元分析方法和*** |
CN111159925A (zh) * | 2019-11-22 | 2020-05-15 | 国网江苏省电力有限公司电力科学研究院 | 瓷绝缘子运行荷载条件下机械性能评估的有限元分析方法、***和存储介质 |
CN111125945A (zh) * | 2019-11-28 | 2020-05-08 | 广东电网有限责任公司 | 一种输电线路极限抗风能力评估方法和装置 |
CN112036061A (zh) * | 2020-08-07 | 2020-12-04 | 国网浙江省电力有限公司电力科学研究院 | 全方位长短腿塔线体系有限元建模及动力响应分析方法 |
CN112036061B (zh) * | 2020-08-07 | 2024-05-10 | 国网浙江省电力有限公司电力科学研究院 | 全方位长短腿塔线体系有限元建模及动力响应分析方法 |
CN112269397A (zh) * | 2020-10-30 | 2021-01-26 | 国家电网有限公司 | 一种基于高精度定位和视觉跟踪的无人机巡检方法 |
CN113553747A (zh) * | 2021-08-27 | 2021-10-26 | 贵州电网有限责任公司 | 一种覆冰输电塔力学特性分析及状态评估方法 |
CN113642222A (zh) * | 2021-08-27 | 2021-11-12 | 贵州电网有限责任公司 | 综合考虑微地形与荷载的基于有限元分析的杆塔选型方法 |
CN114818420A (zh) * | 2022-04-18 | 2022-07-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
CN114818420B (zh) * | 2022-04-18 | 2024-03-29 | 内蒙古电力(集团)有限责任公司内蒙古电力经济技术研究院分公司 | 输电导线初始构型快速建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103455686A (zh) | 架空输电塔-线耦合体系有限元模型的建模方法 | |
Meng et al. | Sensitivity analysis of wind pressure coefficients on CAARC standard tall buildings in CFD simulations | |
CN105095576B (zh) | 一种输电铁塔杆件应力计算方法 | |
CN107346357B (zh) | 一种基于整体耦合模型的海上风机疲劳分析*** | |
CN103676667B (zh) | 变电站工程软导线装配三维仿真测控*** | |
Ju et al. | Fatigue design of offshore wind turbine jacket-type structures using a parallel scheme | |
CN103106305A (zh) | 一种基于实测模态的空间网格结构模型分步修正方法 | |
CN105224741A (zh) | 电动汽车驱动***电磁辐射测试规划方法 | |
CN106250590A (zh) | 一种基于cfd数值模拟的高空露天平台行人风环境评估方法 | |
CN102682175B (zh) | 基于屈曲模态组合的网格结构施工误差可靠性分析方法 | |
CN104361157A (zh) | 一种建筑物间风环境的评价方法 | |
CN104967114B (zh) | 电网负荷实时数字建模方法及*** | |
Alipour et al. | Analytical approach to characterize tornado-induced loads on lattice structures | |
CN109726435B (zh) | 一种接头钉群载荷计算的方法 | |
Li et al. | Wind-induced collapse mechanism and failure criteria of super-large cooling tower based on layered shell element model | |
CN112214843B (zh) | 一种用于风洞试验机翼模型有限元刚度修正方法及装置 | |
Chen et al. | Mechanical behavior of welded hollow spherical joints with diameter exceeding 1.0 m | |
CN105302979A (zh) | 两相流体网络模型中阀门组的建模方法和*** | |
CN212375727U (zh) | 一种用于钢-混组合桁架桥梁桁架节点位置调控的*** | |
CN110851904A (zh) | 一种外框内筒电视塔结构参数灵敏度快速分析方法及*** | |
CN113553747A (zh) | 一种覆冰输电塔力学特性分析及状态评估方法 | |
CN103806423B (zh) | 自升式平台三角桁架式桩腿等效的方法 | |
CN104732039A (zh) | 一种批量计算式桩基础设计方法及数据处理装置 | |
CN103745068A (zh) | 一种套筒铆接结构的有限元建模方法 | |
CN109238627A (zh) | 一种应用于风场特性风洞试验的多点同步排管测压装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20131218 |
|
RJ01 | Rejection of invention patent application after publication |