球坐标系下的面积微元
油浸式变压器的热点温度是其运行过程中的主要参数,在一定程度上决定着变压器的绝缘老化水平。
变压器内部结构极其复杂,且多为导电导磁的材料,这就使得漏磁场极易穿过内部结构件形成回路从而产生涡流损耗,最终导致变压器局部温度升高。
当热点温度高于基准值时,每当温度增加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":"变压器热点温度建模与仿真,对多物理场耦合及温升特性的影响
卜翠旭1308求半径为a的球的表面积(要求用直角坐标系) -
包田昌17826786270 ______ 我们已经知道的就是xyz的坐标系了,也叫空间直角坐标系.那么我们以A点为 那么球的表面积就是S=4πr^2=28π 体积是V=(4πr^3)/3=7根号14
卜翠旭1308球坐标系φ的范围
包田昌17826786270 ______ 三重积分有3个独立变量,在直角坐标系中分别是dux、zhiy、z,在球坐标中是r,ψ,θ.它们的取值范围,关键是取决于“积分区域”,对于“积分区域”要有一个直观形象地概念,实际上就是一个三维立体图形所占的空间区域.球坐标(r,θ,φ),φ是与正x轴之间的方位角,θ是与正z轴的夹角,r是与原点距离.球坐标系是三维坐标系的一种,用以确定三维空间中点、线、面以及体的位置,它以坐标原点为参考点,由方位角、仰角和距离构成.球坐标系在地理学、天文学中都有着广泛应用.
卜翠旭1308高数: 如何理解拉梅系数? -
包田昌17826786270 ______ (1).就是说x1=p, x2=q, x3任取 x1=p, x2任取, x3=r x1=p, x2=q, x3任取 p,q,r为任意常数 形成的曲线网,在空间中的(p,q,r)点是正交的. 曲线网本身可以是弯曲的. 你可以想象一下球坐标系,有助于理解. (2).我理解x1, x2, x3是原先那点,在...
卜翠旭1308怎样用球坐标求球体面积和体积RT, -
包田昌17826786270 ______[答案] 高数的东西忘记的差不多了,球坐标知道之后,可否计算出球体半径(重点以及难点),如果能算出半径,其他的都迎刃而解了.还需要你自己结合你的问题进一步研究.
卜翠旭1308经纬度坐标下的球面多边形面积计算公式 -
包田昌17826786270 ______ 对于球面n边形,所有角的和为S,球的半径为R,那么其面积就是 R^2*(S-(n-2)*Pi)
卜翠旭1308二重积分转换成极坐标计算的面积元素,三重积分转换成柱坐标、球面坐标计算的体积元素是怎么得出来的? -
包田昌17826786270 ______ 球面坐标计算的体积公式=∫∫∫_V dV 此处是球体,那么利用球坐标 =∫∫∫ ρ^2 sin φ dρdφdθ =∫dθ ∫sin φdφ ∫ ρ^2dρ =2π*[-cosφ |]*[ρ^3/3 |] =2π*2*r^3/3 =4πr^3/3 扩展资料 球面坐标系是三大常用的坐标系之一,其它二个常用的坐标系是标准的欧氏坐标系、柱面坐标系.球面坐标变换公式描述了空间中一点P在欧氏坐标系下的坐标 与球面坐标系下的坐标 之间的变换关系.该变换关系如下述公式给出 : 或者,将表达成的形式: 参考资料来源:百度百科—球面坐标变换
卜翠旭1308利用球面坐标计算三重积分球面坐标系中的体积元素:dv=r^2sinkdrdkdm纬线方向的宽为rsinkdm 是怎么得出来的? -
包田昌17826786270 ______[答案] 球面坐标系 x=rsinkcosm y=rsinksinm z=rcosk 然后是rsink是x,y,z的关于r,k,m雅克比(JOCOBI行列式)的值
卜翠旭1308大学高数中如何用微元法求证球的体积 -
包田昌17826786270 ______ 球坐标系 sinordrdo o到牌 0到2牌 0到r
卜翠旭1308二重积分直角坐标转化成极坐标后为什么多了一个r -
包田昌17826786270 ______ dxdy=rdrdθ 根据极坐标和直角坐标的转化公式,代人d的不等式中即可,极坐标的基本公式x=rcosθ,y=rsinθ,由此可知x²+y²=r^2,代人x²+y²≦x+y中有r^2≤rcosθ+rsinθ,由于r≥0,所以0≦r≦sinθ+cosθ. 面积微元从直角坐标系转化为极坐标...
卜翠旭1308球坐标知识求解已知球半径为R,球面两点A(α,β)B(φ,θ),
包田昌17826786270 ______ 首先,你没有说明一个非本质问题:哪个是经度,哪个是纬度,纬度是以北极为0度(高等数学上通常这样讲),还是以赤道为0度(数学分析通常这样讲). 我下面按高等数学习惯给你讲,经度:α,0°