
应用研究
数事其本与成用
35kV绝缘子电位和电场的有限元求解研究
李文学
(三星数据系统(中国)有限公司北京100015)
摘要:本文简要介绍了趋嫁子在电气工程领域的作用和运行方式并利用基于有限元方法的原理,计算出了稳态时系数型编微分方程各项在一筛单元、三角形单元、四边形单元下的单元矩阵,并利用图标详细讲解了达一矩阵。本文还利用单元矩阵计算了3张V绝缘子的电位和电场分布,并给出了基于此原理的电气工程中的应用实际
关键词有限元系数型偏微分方程单元矩阵35kV绝缘子中图分类号:TM216
文献标识码:A
文章编号:1007-9416(2014)11-0094-01
绝缘子作为电力系统外绝缘的重要设备,起着紧固支撑和电气绝缘的作用,绝缘子表面的电位和电场分布对绝缘子伞形的设计、电气性能的改善有着重要作用,因此有必要结合工程实际,从理论上分析和计算绝缘子的电位和电场分布。有限元方法作为一种求解偏微分方程的数值方法,可以求解绝缘子的电位和电场分布,其基本核心是将求解域离散为单元,采用加权余量法形成单元矩阵,将不同的单元矩阵组合起来,构成整个求解域的内度矩阵并进行代数方程的求解。
不同的单元会形成不同的单元矩阵,本文以一维和二维单元为例,从原理上求解出了系数型偏微分方程不同项的单元矩阵,应用求解出来的二维单元矩阵,计算了35kV绝缘子的电位和电场分布,
得到了绝缘子表面的电场强度的最大值。 1有限元法的基本原理
1.1加权余量法
工程中的大部分问愿都能够以数理函数的形式表示成某个未知函数在求解域满足偏微分方程组,并且该未知函数应满足一定的边界条件,常用的办法是对未知函数进行假设并进行求解。常用的方法有配点法,子域法,最小二乘法、伽辽金法,本文基于伽辽金方法推导出不同维数下的各个单元的单元矩阵。
1.2系数型偏微分方程
考虑到工程的复杂情况,其偏微分方程形式各不相同,综合不同形式的偏微分方程,通过修改各项系数,则可获得不同的偏微分方程以满足工程需要
根据未知函数阶数的不同,可以将系数型偏微分方程分为四种
不同的形式,因此需要求出四种形式下的单元矩阵。 2一维单元的单元矩阵求解
-维单元的基本形式为U=C+cX,单元节点X处的值为
[, = c, + c,x
u,在节点处,的值为U2,则有:=G,+e
在一维单元下,采
用伽辽金法分别求偏微分方程各项的单元矩阵。
(a)电位分布云图
(b)电场强度分布云图
图1绝缘子的电位和电场强度分布云图
收福日期:20141016
3二维单元的单元矩阵求解
nene
二维情况下,未知函数u是X,的函数,因此V-Va
Ww-wid
e+
-会+号了,不同的单元形状会有不同的单元矩阵。三角形单元,对于源项f有了了s(x,ytuy=
[[s(5,)dedn=[11/21/2],四边形单元有
[u, =6, +,, +c,, +,, =, +c,X, + c,, +,,)
u, , +c,, +e,, +e,x,, 和u= (11), +
U =c, +c,X +c,Y+c,XY, 1
[u, = c +c,, + c,, + c,3,3
(1n)u, +, +(1),
四边形单元单元矩阵的求法和三角形单元类似,本文限于篇
幅,不再详细介绍。 4工程应用
本文利用上述推导出的单元矩阵,计算了35kV绝缘子的电位和电场分布,计算中假设绝缘子表面干净、空气湿度较低,可以当做静电场进行求解。
35kV绝缘子是一个轴对称模型,为节约计算资源,可采用二维轴对称的方法进行建模,然后对模型划分网格,计算出来的绝缘子的电位和电场强度分布如图1所示,电场强度的最大值为 957.85kV/m,小于空气的起晕场强3000kV/m,因此不需要在
高压端加装均压环来改善绝缘子表面的电场强度。 5结语
本文根据有限元方法的原理,推导出了常见单元下的单元矩阵,为后续工程计算提供了理论基础。根据求解的单元矩阵,计算了 35kV绝缘子的电位和电场分布,绝缘子表面的电场强度最大值为 957.85kV/m,小于空气的起晕场强,因此不需要加装均压环,该结论符合我国电力系统中220kV及以上电压等级的绝缘子才需要加装
均压环的事实。参考文献
[1]郭历伦.陈忠富,罗景润.陈刚.扩晨展有限元方法及应用综述[]].力学季刊,2011.04;612625
[2]李录贤,刘书静.张慧华.陈方方,王铁军.广义有限元方法研究进展[J].应用力学学报,2009,01:96-108+214.
[3]王建军.于长液.李其汉.工程中的随机有限元方法[J].应用力学学报.2009,02:297303+408
[4]肖静宇.儿类分数阶微分方程的数值方法研究[D].哈尔滨工业大学2013.
[5]蔚喜军,周铁.流体力学方程的间断有限元方法[J].计算物理, 2005,02:108116
作者简介:季文学(1977一),男,朝鲜族,吉林延边人,颈士,中级职称,研究方向为电气工程 94