球坐标系拉普拉斯算子
油浸式变压器的热点温度是其运行过程中的主要参数,在一定程度上决定着变压器的绝缘老化水平。
变压器内部结构极其复杂,且多为导电导磁的材料,这就使得漏磁场极易穿过内部结构件形成回路从而产生涡流损耗,最终导致变压器局部温度升高。
当热点温度高于基准值时,每当温度增加6℃,变压器的老化率就增加一倍,这是经典的六度法则。
![](https://p0.ssl.img.360kuai.com/t01174767da6e16be6e.jpg?size=756x486)
局部过热是危害变压器运行的主要因素,因此,测算变压器的热点温度有利于及时掌握变压器的运行情况,从而尽量避免热故障的发生。
一方面,变压器长期处于运行状态,并承受交直流的共同作用,高热环境使其绝缘材料的老化加快,而其运行时的温度又决定了绝缘的老化程度,也就是变压器的寿命;另一方面,研究热点温度可以对冷却系统的运行状况进行判断,指导变压器设计。
![](https://p0.ssl.img.360kuai.com/t017a25d4af2f2d15b6.jpg?size=788x578)
因此,研究基于多物理场耦合计算的变压器温升特性对变压器的绝缘、散热、老化等问题均有重要意义。
目前,国内外关于油浸式变压器的热点温度计算方法主要分为4类:多物理场建模计算方法、热路模型法、人工智能算法和经验公式法。
2002年建立了油浸式电力变压器的模型,通过计算分析了变压器负载系数、环境温度等因素对绝缘热点温度的影响。
![](https://p0.ssl.img.360kuai.com/t0183e9ae12779d6b22.jpg?size=728x596)
2014年,提出了用户定义函数(UDF的算法,将热场与电磁场和流体场耦合,讨论了油浸变压器温度特性的多物理场分析方法,以预测温升。
2016年,采用有限元法和Ansys软件进行三维模型开发,对160KVN油浸式配电变压器进行温度场分布研究,得到了随温度变化的热场分布.
2019年,采用有限元法及有限体积法分别就典型变压器在出厂温升试验和现场实际运行两种工况下的绕组损耗及热点温升进行了仿真计算,并对计算结果进行了对比分析。
![](https://p0.ssl.img.360kuai.com/t01f36de5cacdc11901.jpg?size=692x572)
通过对变压器热点温度计算和多物理场耦合研究的国内外现状分析,本文拟以变压器热点温度为中心,从理论上分析变压器内部产热机理和内部热量传递过程,探讨油流在变压器内部的对流换热行为。
充分考虑变压器内部多物理场耦合情况,建立一种绕组热点温度的改进模型,并应用流固热耦合温度场分析方法对油浸式变压器热点温度进行仿真分析,具有较高模型精度的同时节省了计算量。
![](https://p0.ssl.img.360kuai.com/t01433e3648674fa895.jpg?size=706x596)
![](https://p0.ssl.img.360kuai.com/t0180240e38cd928b83.jpg?size=1080x137)
●—≺ 变压器的热特性分析 ≻—●
对变压器温度分布计算的研究是基于能量守恒原理。变压器在实际运行的过程中,电磁能量的一部分要转化为热量,这些热量是铁芯、绕组及相关结构件中产生的损耗所引起的。
当这些损耗引起的热量分散到各部位周围的介质中便会导致变压器内部发热及温度升高。
![](https://p0.ssl.img.360kuai.com/t0181c51ec23d2120c6.jpg?size=1298x780)
变压器的冷却措施即是用来避免变压器因温度升高而造成损坏的系统。变压器在工作时会有负载电流流过,会有磁通在绕组周围产生,这种由负载电流产生的磁通即是漏磁通,由漏磁通引起的损耗是变压器温升的主要热源。
变压器中产生的热量向周围介质中传递的方式(即散热)分为3种类型,分别是热传导、热对流和热辐射。
热传导是指由于温度差引起的热能传递现象。传导散热的效率取决于接触两个机体的接触面积、温度差以及两者的材料属性。
假设两个物体的接触面积为S,距离为x,温度分别为θ1、θ2,则其通过传导散热方式传递的热量表示为:
![](https://p0.ssl.img.360kuai.com/t0102a6dcf7a2c92a41.png?size=163x66)
式中:Q是传递的热量;λ是传导散热系数。
热对流是热量通过流动介质传递的过程。在变压器中,高温物体附近的流体受热膨胀,密度降低并向上流动,同时密度较大的冷流体下降并代替原来的受热流体,使得高温物体表面常常发生对流现象。可以表述为:
![](https://p0.ssl.img.360kuai.com/t010ca8c51a36e40ba7.png?size=150x26)
式中:h是对流换热系数;tw是固体表面平均温度;tf是流体平均温度;A是对流散热面积。
热辐射是物体由于具有温度而辐射电磁波的现象。辐射散热是双向的,它通过电磁波的方式将热量由热源向直线方向放射热量,在其散热的同时也在吸收热量,其传递出去的热量减去吸收的热量即为净传递热量,可表述为:
![](https://p0.ssl.img.360kuai.com/t015d045418ff5a691a.png?size=253x58)
式中:C0是辐射系数;T1为发射温度;T2为接收温度;F为固体表面积。变压器内部散热形式示意图大致如图1所示。
铁芯、绕组与变压器油之间主要呈现热对流和热传导的散热方式,整个油箱与油箱外界的空气主要呈现的是热对流与热辐射的散热方式。
![](https://p0.ssl.img.360kuai.com/t01b4b133bfa79b4adb.jpg?size=750x680)
![](https://p0.ssl.img.360kuai.com/t015ca2a79c8bff5c3e.jpg?size=1080x137)
●—≺ 变压器温度场与流场数学模型 ≻—●
在进行变压器温升特性研究时,需要考虑流体、固体以及传热等多方面的问题,还需要满足质量守恒、能量守恒以及动能守恒。
质量守恒定律在油浸式变压器温度场和流场的耦合场中的数学表达。单位时间t内变压器油中体积为V的部分通过控制面S流入变压器内部的质量之和等于单位时间内油中质量的增量。其数学表达式为:
![](https://p0.ssl.img.360kuai.com/t01f5cb7f9f6fc91cc3.png?size=270x53)
式中:ρ为密度;n指的是微元面积矢量dS外法线的单位向量;U指的是微元表面dS上流体的流速。该方程的微分形式为:
![](https://p0.ssl.img.360kuai.com/t019615e781e363ef5a.png?size=248x57)
式中:u、v、w分别为油流速在x、y、z方向的分量。
动量守恒定律在变压器油流流速计算中的数学表达。在变压器内部,动量方程可以表示在变压器油中取体积为V的控制体,对于该控制体,单位时间内变压器油的动量增量相当于这段时间内通过控制面S流入的动量和控制体的外力合力之和。所以,变压器油流流速计算的数学表达式为:
![](https://p0.ssl.img.360kuai.com/t0158cd15eb1bfeb620.png?size=287x95)
式中:U表示为油流速;f为单位质量变压器油所受到的质量力;Π为微元面积矢量dS的应力张量。X、Y、Z三个方向的运动方程为:
![](https://p0.ssl.img.360kuai.com/t014d190ecd74367431.png?size=268x51)
式中:px为流体内热源及由于粘性作用流体机械能转换为热能的部分。Y、Z方向的运动方程与X方向的方程类似。
能量守恒定律在变压器内部的热量传递过程中的数学表达单位时间内由外界传递给控制体的热量、功以及通过控制面流入控制体的能量之和等于在该时间内控制体中油流能量的增量。其数学表达式为:
![](https://p0.ssl.img.360kuai.com/t016a50e55bd78d7920.png?size=290x52)
式中:T为控制体的温度;c为变压器油的比热容;k为导热系数;Q为对应微元的生热量;V为调和算子。
热传递方程式的建立。在运用热力学第一定律来建立热传递方程式时,整个求解结构中热量的传导速率应等于机械能和热能进入结构的能量速率与离开结构的能量速率之差。其在直角坐标系中的微分方程为:
![](https://p0.ssl.img.360kuai.com/t015eccaf8eab5aed49.png?size=275x100)
式中:τ为温度梯度;qV为单位体积的传热体在单位时间内散出的热量;c为体积比热容;ρ为介质密度;λ为传热系数;
![](https://p0.ssl.img.360kuai.com/t016beee191447b08db.png?size=168x47)
为拉普拉斯算子。
在分析变压器温度场分布时,平面或轴对称圆柱模型更适用于计算变压器这类常物性、有内热源的二维稳态温度场,所以我们将直角坐标系的热传递微分方程变换为圆柱坐标系的热传导微分方程,推导如下:
![](https://p0.ssl.img.360kuai.com/t017c1f852717850a3e.jpg?size=640x820)
轴对称结构的变压器在柱坐标系中温度分布T(r、z),可以通过下式进行求解:
![](https://p0.ssl.img.360kuai.com/t014da0960797ee3b1c.png?size=370x70)
![](https://p0.ssl.img.360kuai.com/t01e83f4008b30507b2.jpg?size=1080x137)
●—≺ 变压器温度场与流场数学模型 ≻—●
对变压器内部变压器油建立流-固耦合的变压器物理模型,如图 2。
![](https://p0.ssl.img.360kuai.com/t01cb87122986b1f9c1.jpg?size=676x528)
变压器内部绕组及铁芯铁轭损耗所产生的热量,经过其热传导、热对流的导热方式传递给油箱内部的变压器油,再通过这两种方式传递给油箱,通过油箱将热量散发出去,从而达到热平衡,进而降低变压器内部的温度,完成散热的目的。
物理参数与材料特性。对变压器流-固-热耦合场分析时,需要对变压器中材料的比热容、导热系数、密度、粘度以及热扩散系数等热性能参数进行设定。根据变压器型号及实际情况给出内部金属材料的参数设定,如表1所示。
![](https://p0.ssl.img.360kuai.com/t01a8560264175adaf5.jpg?size=808x432)
油浸式变压器的散热过程是:铁芯和线圈把热量首先传给在其附近的油,使油的温度升高。
温度高的油体积增加,比重减小,就向油箱的上部运动。冷油将自然运动补充到热油原来的位置。
而热油沿箱壁或散热器管将热量放出,经箱壁被周围的空气带走,温度降低后又回到油箱下部参加循环。这样,因油温的差别,产生了油的自然循环流动,进而达到散热的目的。变压器油的密度、热导率、比热、动力粘度等参数均会随着温度变化而变化。变压器油的材料特性参数与温度T的函数关系,如表2所示。
![](https://p0.ssl.img.360kuai.com/t018c3432cd338a4f3e.jpg?size=798x460)
边界条件的设定。本文采用有限元方法计算换流变压器温度和油流分布,计算前需确定初始条件及边界条件。
在变压器外部热量传导过程中,油箱壁与空气之间的热交换属于对流传热,该传热过程满足第三类边界条件即
![](https://p0.ssl.img.360kuai.com/t018516c7d6b5ae4eb3.png?size=216x58)
式中:Tf为油箱表面空气的温度;α为对流换热系数;T为待求量油箱壁上的温度。
(2)初始温度设置与环境温度相等为300K。
(3)变压器流场计算设置初始条件为:变压器油的初始速度设置为0m/s;变压器油的重力加速度为98m/s2,方向为沿Z轴负方向。
本文应用流-固-热耦合温度场计算方法对变压器进行温度场仿真计算,通过设定1000次迭代次数,得出在环境温度300K、额定负载情况下,4500kV油浸式变压器温度场及流场的分布情况。
![](https://p0.ssl.img.360kuai.com/t0190911721d5b6eb3a.jpg?size=1282x982)
变压器油箱表面的温度场分布如图3所示。
图3(a)为主视角,图3(b)为后视角。可以看出,变压器满载运行时,油箱上表面油温最大约为3478K。变压器热点温度主要集中在油箱顶部偏中间的位置。油箱体上表面的温度远大于下表面的温度,温度差最大约可达55K,造成这一现象的主要原因是散热条件的差异。
![](https://p0.ssl.img.360kuai.com/t011f40fbd88d6eb4b9.png?size=363x728)
为了得出高准确度的变压器热点温度动态数值,本文给出了一种考虑变压器温升特性和多物理场耦合情况的变压器热点温度计算模型,并利用有限元分析软件对变压器各结构的温度场及变压器油的流场进行了分析和研究,通过计算与分析得到如下结论:
变压器的热点温度主要出现在变压器的阀侧绕组上,为360.8K,铁芯及铁轭上的温度主要集中在铁芯的中间部分,最高温度为3565K;变压器油的温度与其流速有一定的联系,温度高的地方流速越大,最高温度出现在绕组之间油道的出口处和变压器顶端偏下的位置,为3607K。
本文所述方法得到的变压器热点温度值与经验公式计算得到的变压器绕组热点温度值相比,其变化趋势基本一致,但准确度提升了4%。
","gnid":"90b619a2985211c95","img_data":[{"flag":2,"img":[{"desc":"","height":"486","title":"","url":"https://p0.ssl.img.360kuai.com/t01174767da6e16be6e.jpg","width":"756"},{"desc":"","height":"578","title":"","url":"https://p0.ssl.img.360kuai.com/t017a25d4af2f2d15b6.jpg","width":"788"},{"desc":"","height":"596","title":"","url":"https://p0.ssl.img.360kuai.com/t0183e9ae12779d6b22.jpg","width":"728"},{"desc":"","height":"572","title":"","url":"https://p0.ssl.img.360kuai.com/t01f36de5cacdc11901.jpg","width":"692"},{"desc":"","height":"596","title":"","url":"https://p0.ssl.img.360kuai.com/t01433e3648674fa895.jpg","width":"706"},{"desc":"","height":"137","title":"","url":"https://p0.ssl.img.360kuai.com/t0180240e38cd928b83.jpg","width":"1080"},{"desc":"","height":"780","title":"","url":"https://p0.ssl.img.360kuai.com/t0181c51ec23d2120c6.jpg","width":"1298"},{"desc":"","height":"66","title":"","url":"https://p0.ssl.img.360kuai.com/t0102a6dcf7a2c92a41.png","width":"163"},{"desc":"","height":"26","title":"","url":"https://p0.ssl.img.360kuai.com/t010ca8c51a36e40ba7.png","width":"150"},{"desc":"","height":"58","title":"","url":"https://p0.ssl.img.360kuai.com/t015d045418ff5a691a.png","width":"253"},{"desc":"","height":"680","title":"","url":"https://p0.ssl.img.360kuai.com/t01b4b133bfa79b4adb.jpg","width":"750"},{"desc":"","height":"137","title":"","url":"https://p0.ssl.img.360kuai.com/t015ca2a79c8bff5c3e.jpg","width":"1080"},{"desc":"","height":"53","title":"","url":"https://p0.ssl.img.360kuai.com/t01f5cb7f9f6fc91cc3.png","width":"270"},{"desc":"","height":"57","title":"","url":"https://p0.ssl.img.360kuai.com/t019615e781e363ef5a.png","width":"248"},{"desc":"","height":"95","title":"","url":"https://p0.ssl.img.360kuai.com/t0158cd15eb1bfeb620.png","width":"287"},{"desc":"","height":"51","title":"","url":"https://p0.ssl.img.360kuai.com/t014d190ecd74367431.png","width":"268"},{"desc":"","height":"52","title":"","url":"https://p0.ssl.img.360kuai.com/t016a50e55bd78d7920.png","width":"290"},{"desc":"","height":"100","title":"","url":"https://p0.ssl.img.360kuai.com/t015eccaf8eab5aed49.png","width":"275"},{"desc":"","height":"47","title":"","url":"https://p0.ssl.img.360kuai.com/t016beee191447b08db.png","width":"168"},{"desc":"","height":"820","title":"","url":"https://p0.ssl.img.360kuai.com/t017c1f852717850a3e.jpg","width":"640"},{"desc":"","height":"70","title":"","url":"https://p0.ssl.img.360kuai.com/t014da0960797ee3b1c.png","width":"370"},{"desc":"","height":"137","title":"","url":"https://p0.ssl.img.360kuai.com/t01e83f4008b30507b2.jpg","width":"1080"},{"desc":"","height":"528","title":"","url":"https://p0.ssl.img.360kuai.com/t01cb87122986b1f9c1.jpg","width":"676"},{"desc":"","height":"432","title":"","url":"https://p0.ssl.img.360kuai.com/t01a8560264175adaf5.jpg","width":"808"},{"desc":"","height":"460","title":"","url":"https://p0.ssl.img.360kuai.com/t018c3432cd338a4f3e.jpg","width":"798"},{"desc":"","height":"58","title":"","url":"https://p0.ssl.img.360kuai.com/t018516c7d6b5ae4eb3.png","width":"216"},{"desc":"","height":"982","title":"","url":"https://p0.ssl.img.360kuai.com/t0190911721d5b6eb3a.jpg","width":"1282"},{"desc":"","height":"728","title":"","url":"https://p0.ssl.img.360kuai.com/t011f40fbd88d6eb4b9.png","width":"363"}]}],"original":0,"pat":"art_src_0,sexf,sex4,sexc,fts0,sts0","powerby":"pika","pub_time":1703260920000,"pure":"","rawurl":"http://zm.news.so.com/358ac55766af76b56c24dfa46e3ae91c","redirect":0,"rptid":"f64c14cb1c585550","rss_ext":[],"s":"t","src":"观史通鉴","tag":[],"title":"变压器热点温度建模与仿真,对多物理场耦合及温升特性的影响
苍应菡3322什么是拉普拉斯算子? -
曾池胡19660874319 ______ 首先介绍hamilton算子,埃,怎么说呢,太难打出来了.hamilton算子就是偏X,偏Y,偏Z,laplace算子就是偏偏X,偏偏Y,偏偏Z,举个例子,有一个二阶可偏导函数U,用laplace 算子就是Uxx+Uyy+Uzz.
苍应菡3322微分算子并不是变换矩阵或者坐标系,仅仅是一个算子,给出了坐标系的...
曾池胡19660874319 ______ 圆球坐标系: 圆柱坐标系:
苍应菡3322拉普拉斯算子————————————————————————————————————————纯工科生,高数没好好学~现在做渗流力学碰见拉普... -
曾池胡19660874319 ______[答案] 你说的倒三角叫nabla,是哈密尔顿引入的一个算符,和四元数有关,讲出来会让你更糊涂.总之,如你理解是个简写的符号. 拉普拉斯算子作用在某个函数f(x,y,z)上(拿三维举个例子),就是将这个函数对每个变量求二阶偏导数,然后求和,仅此而已. ...
苍应菡3322拉普拉斯算子的实际应用是什么?
曾池胡19660874319 ______ 拉普拉斯算子有许多用途,是椭圆型算子中的一个重要例子. 在物理中,常用于波方程的数学模型、热传导方程以及亥姆霍兹方程. 在静电学中,拉普拉斯方程和泊松方程的应用随处可见.在量子力学中,其代表薛定谔方程式中的动能项. 在数学中,经拉普拉斯算子运算为零的函数称为调和函数;拉普拉斯算子是霍奇理论的核心,并且是德拉姆上同调的结果.
苍应菡3322拉普拉斯算子的物理意义是什么? -
曾池胡19660874319 ______ 拉普拉斯算子从形式上看表示,一个场变量的梯度的散度.散度的概念是很清晰的,从高斯方程应用到静电场领域可以知道,散度可以表示一个矢量在单位空间内产生通量的强度,静电场中因为一个封闭的曲面内部有静电荷,那么这个封闭曲面...
苍应菡3322拉普拉斯算子的实际应用是什么? -
曾池胡19660874319 ______ 在物理中,常用于波方程的数学模型、热传导方程以及亥姆霍兹方程.在静电学中,拉普拉斯方程和泊松方程的应用随处可见.在量子力学中,其代表薛定谔方程式中的动能项.在数学中,经拉普拉斯算子运算为零的函数称为调和函数;拉普拉斯算子是霍奇理论的核心,并且是德拉姆上同调的结果.
苍应菡3322哈密尔顿算符,拉普拉斯算子,梯度和散度 -
曾池胡19660874319 ______ 哈密顿算符的谱为测量系统总能时所有可能结果的集合.如同其他自伴算符,哈密顿算符的谱可以透过谱测度被分解,成为纯点、绝对连续、奇点三种部分. 拉普拉斯算子是n维欧几里德空间中的一个二阶微分算子,定义为梯度(▽f)的散...