CN104834781B - 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 - Google Patents
一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 Download PDFInfo
- Publication number
- CN104834781B CN104834781B CN201510232180.0A CN201510232180A CN104834781B CN 104834781 B CN104834781 B CN 104834781B CN 201510232180 A CN201510232180 A CN 201510232180A CN 104834781 B CN104834781 B CN 104834781B
- Authority
- CN
- China
- Prior art keywords
- heat
- smoothing reactor
- temperature
- oil
- mrow
- 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.)
- Active
Links
- 238000009499 grossing Methods 0.000 title claims abstract description 95
- 238000004458 analytical method Methods 0.000 title claims abstract description 45
- 230000001052 transient effect Effects 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000009826 distribution Methods 0.000 claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 239000007787 solid Substances 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000004519 manufacturing process Methods 0.000 claims abstract description 7
- 239000003921 oil Substances 0.000 claims description 90
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 64
- 238000012546 transfer Methods 0.000 claims description 45
- 238000004804 winding Methods 0.000 claims description 35
- 239000012530 fluid Substances 0.000 claims description 32
- 230000005855 radiation Effects 0.000 claims description 21
- 230000004907 flux Effects 0.000 claims description 18
- 229910052742 iron Inorganic materials 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 11
- 239000002828 fuel tank Substances 0.000 claims description 7
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 6
- 238000001816 cooling Methods 0.000 claims description 6
- 229910052802 copper Inorganic materials 0.000 claims description 6
- 239000010949 copper Substances 0.000 claims description 6
- 230000017525 heat dissipation Effects 0.000 claims description 5
- 229910000976 Electrical steel Inorganic materials 0.000 claims description 4
- 239000012080 ambient air Substances 0.000 claims description 3
- 239000000295 fuel oil Substances 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 3
- 230000005389 magnetism Effects 0.000 claims description 3
- 239000003758 nuclear fuel Substances 0.000 claims description 3
- 230000007257 malfunction Effects 0.000 abstract description 6
- 239000003570 air Substances 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 230000005611 electricity Effects 0.000 description 5
- 239000007788 liquid Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 206010037660 Pyrexia Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 208000032368 Device malfunction Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000005619 thermoelectricity Effects 0.000 description 1
- 230000032258 transport Effects 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Housings And Mounting Of Transformers (AREA)
Abstract
本发明公开了一种基于平波电抗器多次换向失败时的暂态温度场仿真方法，包括以下步骤：A：基于ANSYS有限元软件，采用实体建模的方式建立油浸式平波电抗器三维几何模型；B：对油浸式平波电抗器进行产热分析；C：对油浸式平波电抗器进行散热分析；D：确定温度场仿真计算的边界条件；E：计算正常工况下平波电抗器稳态温度场；F：将求得稳态温度场结果作为故障暂态模型的初始值，然后加载故障电流，采用分步长的方法，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。本发明能够为查明平波电抗器瓦斯继电器误动的原因和为平波电抗器的瓦斯校验整定值设置提供依据。
Description
技术领域
本发明涉及电力系统仿真领域，尤其涉及一种基于平波电抗器多次换向失败时的暂态温度场仿真方法。
背景技术
目前,我国已建成多项大型直流输电工程，随着这些直流输电工程的投运，直流输电的运行稳定工作显得日益重要。统计数据表明，平波电抗器瓦斯继电器误动是故障主要因素之一，其原因主要是平波电抗器瓦斯继电器无法避开连续发生换向失败而导致瓦斯动作，而事后对油质的分析结果却都是正常的，并无大量油气产生。引发上述故障的具体动作机理无法确定，是因电流产生的热效应导致油体积膨胀，还是因大电流突变产生的电动力使带电体变形造成对变压器油的挤出，都无法考证。另一方面，在进行平波电抗器的瓦斯校验整定值设置时，目前的方法还是借鉴常规油浸式变压器，且依靠于经验值设置，效果不佳。
平波电抗器结构与变压器类似，对变压器的研究也适用于平波电抗器，变压器温度场描述方法也可用于平波电抗器。现有变压器故障的研究主要集中在通过油中溶解气体在线监控与色谱分析来进行故障识别，很少有基于油流速来进行故障识别。由于油流速是故障发生时的一个典型特征，影响油流速的产气和油膨胀都是由故障发热引起的，但平波电抗器故障发热是一个暂态过程，其故障发热过程复杂。因此故障后暂态温度场的仿真研究对查明平波电抗器气体继电器误动作原因具有重要意义。现有研究都只考虑变压器正常运行时的稳态温度场描述，并没有涉及对故障情况下变压器暂态温度场进行计算。
综上，考虑到多次换向失败引起的平波电抗器气体继电器误动作的现实问题，有必要对故障情况下平波电抗器内部暂态温度场进行仿真计算，为查明平波电抗器瓦斯继电器误动的原因和为平波电抗器的瓦斯校验整定值设置提供依据。
发明内容
本发明的目的是提供一种基于平波电抗器多次换向失败时的暂态温度场仿真方法，能够为查明平波电抗器瓦斯继电器误动的原因和为平波电抗器的瓦斯校验整定值设置提供依据。
本发明采用下述技术方案：
一种基于平波电抗器多次换向失败时的暂态温度场仿真方法，包括以下步骤：
A：基于ANSYS有限元软件，采用实体建模的方式建立油浸式平波电抗器三维几何模型；然后对油浸式平波电抗器三维几何模型进行材料物理参数确定、单元类型确定以及网格划分；
B：对油浸式平波电抗器进行产热分析，分别确定油浸式平波电抗器中绕组电流产生的铜损、铁芯产生的铁损和由于漏磁通在电抗器结构件中产生的杂散损耗；
C：对油浸式平波电抗器进行散热分析，分别确定油浸式电抗器中的热传导、热辐射和热对流；
D：确定温度场仿真计算的边界条件；
E：计算正常工况下平波电抗器稳态温度场；
F：将求得稳态温度场结果作为故障暂态模型的初始值，然后加载故障电流，采用分步长的方法，使用ANSYS有限元软件的瞬态分析功能，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。
所述的步骤A中，油浸式平波电抗器三维几何模型包括绕组、铁芯、油箱、油、套管和冷却系统；在对油浸式平波电抗器三维几何模型中材料物理参数进行确定时，依次输入绕组、铁芯、油箱和油的电阻率及热传导系数，并输入铁芯非线性磁导率。
所述的步骤A中，在ANSYS有限元软件中，以铁芯的B-H曲线的形式输入非线性磁导率，其中铁芯的B-H曲线以二维数组的形式拟合。
所述的步骤B中，绕组电流产生的铜损p的计算式为：，I1为绕组的相电流；r1为绕组的总电阻；铁芯产生的铁损P0的计算式为：P0＝K0GPc，K0为铁芯损耗的工艺系数，G为铁芯的重量，Pc为单位重量铁耗，P1/50为铁耗系数，为Bm＝1T且f＝50Hz时每千克硅钢片的铁耗，β为频率指数，Bm为最大磁通密度，B指代磁通密度，下角标m指代max；F为频率；漏磁通在电抗器结构件中产生的杂散损耗在总损耗中所占的比重较小，在此不再考虑。
所述的步骤C中，对油浸式平波电抗器进行散热分析，分别确定油浸式电抗器中的热传导、热辐射和热对流；
(1)油浸式电抗器中的热传导：
热传导在导热介质两侧温度差(△θ＝θ1-θ2)可以用以下公式计算：
△θ＝Q×RT，
其中△θ为热传导在导热介质两侧温度差，θ1为导热介质一侧温度，θ2为导热介质另一侧温度，Q为绝缘层或油箱壁的热通量，RT为绝缘层或油箱壁热阻，下角标T代表绝缘层或油箱壁，β1为导热介质厚度，A为导热介质的横截面积；k为热导率；
定义q为单位面积的热通量即热流密度，则上式变为
油浸式电抗器油纸绝缘的热导率k为定值；
(2)油浸式电抗器中的热辐射：
辐射的热转移斯特凡波尔兹曼定律表示为：
其中，PR单位时间内物体表面辐射出的总能量，下角标R仅用于区分并非为变量，斯特凡波尔兹曼常数η＝5.67×10-8w/m2·k4，E为表面发射率，AR为辐射表面积，Ts为辐射表面平均温度，Tα为周围空气温度，下角标s和下角标α表示仅用于区分并非为变量；
(3)油浸式电抗器中的热对流
对流换热公式如下：
Q＝α(tw-tf)S；
其中，Q为对流换热量，单位为W；tw和tf为壁面和流体平均温度，单位为℃，w表示油箱壁壁面，f表示流体，下角标w、f非变量，仅用于区分两个温度；S为对流换热面积，单位为㎡；α为对流换热系数，单位为W/㎡·℃。
所述的步骤D中，温度场仿真计算的边界条件分别为：
(1)第一类边界是已知任何时刻物体边界面的温度分布，即
T|S1＝T0；
其中，T为物体温度，S1为物体的边界面，T0在稳态热传导过程中表示给定的温度，在非稳态导热过程中表示温度值；
(2)第二类边界是已知任何时刻物体边界面上的热流密度分布，即
其中，λ为导热系数，T为物体温度，表示温度沿n方向的梯度，n为热流传递方向；S2为物体的边界面；q0为通过边界面S2的热流密度，对于稳态导热过程，q0为常量；对于非稳态过程，q0是随时间变化的量。
(3)第三类边界条件是已知物体与周围流体间的对流换热系数和流体温度，即
其中，λ为导热系数；T为物体温度；表示温度沿n方向的梯度，n为热流传递方向；S3为物体的边界面；、α为物质的散热系数，Tf为流体温度，T为物体温度。
所述的步骤E中，使用ANSYS有限元分析软件，利用上述步骤中确定的网格划分后的油浸式平波电抗器三维几何模型，输入材料物理参数和边界条件，然后在平波电抗器绕组上加载稳态电流，并在平波电抗器油箱外侧远场模型上加载外界环境温度；外界环境温度加载后，在产热的同时考虑散热对温度的影响，根据步骤C中确定的热传导公式、热辐射公式和对流换热公式，计算出平波电抗器内部热量交换以及油箱壁的散热；最后通过ANSYS有限元分析软件计算出稳态温度场的结果，即平波电抗器三维模型内各个点的温度值，并以三维图像的形式显示出来；其中，稳态电流指油浸式平波电抗器正常运行时绕组上通过的电流。
所述的步骤F中，将求得稳态温度场结果作为故障暂态模型的初始值，然后加载故障电流，采用分步长的方法，使用ANSYS有限元软件的瞬态分析功能，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。
本发明通过使用实体建模的方式建立油浸式平波电抗器三维几何模型，然后经材料物理参数确定、单元类型确定、网格划分和边界条件确定，通过在平波电抗器绕组上加载稳态电流，利用用ANSYS有限元分析软件计算出稳态温度场的结果；再将求得稳态温度场结果作为故障暂态模型的初始值，通过加载故障电流并采用分步长的方法，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。本发明为故障情况下平波电抗器温度场计算提供了一种方法，为通过油流速整定平波电抗器气体继电器并防止其误动作提供了相应的基础模型，同时也为查明平波电抗器瓦斯继电器误动的原因和为平波电抗器的瓦斯校验整定值设置提供依据。
附图说明
图1为本发明的流程示意图。
具体实施方式
以下结合附图和实施例对本发明作以详细的描述：
如图1所示，本发明所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于，包括以下步骤：
A：基于ANSYS有限元软件，采用实体建模的方式建立油浸式平波电抗器三维几何模型；然后对油浸式平波电抗器三维几何模型进行材料物理参数确定、单元类型确定以及网格划分；
油浸式平波电抗器的结构与变压器相似，油浸式平波电抗器三维几何模型主要包括绕组、铁芯、油箱、油、套管和冷却系统；
在对油浸式平波电抗器三维几何模型中材料物理参数进行确定时，依次输入绕组、铁芯、油箱和油的电阻率及热传导系数，并输入铁芯非线性磁导率，在ANSYS有限元软件中，以铁芯的B-H曲线的形式输入非线性磁导率，其中铁芯的B-H曲线以二维数组的形式拟合。
在确定油浸式平波电抗器三维几何模型的材料物理参数后，还需确定油浸式平波电抗器三维几何模型的单元类型并进行网格划分。本发明中，油浸式平波电抗器三维几何模型的单元类型采用三维热电耦合单元。在进行网格划分过程中，由于油浸式平波电抗器三维几何模型中绕组近似于圆环形体，可对绕组采用先进行面剖分，再通过体扫掠完成体的剖分；由于油的形状极不规则，因此先将其拆分成规则形状，然后再分别进行剖分控制。在进行面剖分的时候，通过设置边线的长度或分段数，保证模型结构之间的网格能够合理的过渡。进行旋转体剖分的时候，通过设置沿旋转方向生成单元个数控制网格剖分的精度。一般网格越密，计算的精度会更高，但同时所需的计算时间也会增加。本发明中，由于在油箱附近温度变化较小，铁芯跟绕组作为主要发热源，其附近温度变化较大，铁芯跟绕组附近需增加网格密度，离铁芯跟绕组较远的部分可降低网格密度，网格的密度主要通过设置边线的分段数或者生成单元个数来控制。通过上述一系列的剖分控制，保证了模型结构之间的网格平滑过渡，得到质量较好的有限元模型。此外还需在油箱外侧加一层远场单元，用以模拟外界环境温度对平波电抗器的影响。
B：对油浸式平波电抗器进行产热分析，分别确定油浸式平波电抗器中绕组电流产生的铜损、铁芯产生的铁损和由于漏磁通在电抗器结构件中产生的杂散损耗；
其中，绕组电流产生的铜损p的计算式为：I1为绕组的相电流；r1为绕组的总电阻；
铁芯产生的铁损P0的计算式为：P0＝K0GPc，K0为铁芯损耗的工艺系数，G为铁芯的重量，Pc为单位重量铁耗(通过查阅硅钢片型号的铁损曲线即可得到单位重量铁损值Pc)，P1/50为铁耗系数，为Bm＝1T且f＝50Hz时每千克硅钢片的铁耗，β为频率指数，Bm为最大磁通密度，B指代磁通密度，下角标m指代max；F为频率。
油浸式平波电抗器运行时，绕组中的电流产生的磁通不全通过铁芯闭合，相当一部分通过电抗器内部介质和结构件(如油、拉板、夹件、油箱壁等)闭合，在油浸式平波电抗器一定负载运行下，其绕组中的电流产生的漏磁场在导电介质中产生一定的杂散损耗，通常这些损耗在总损耗中所占的比重较小，在此不再考虑。
C：对油浸式平波电抗器进行散热分析，分别确定油浸式电抗器中的热传导、热辐射和热对流；
(1)油浸式电抗器中的热传导：
在油浸式电抗器中热量通过绝缘层流出，热量通常流过几层绝缘层到达冷却介质，油箱壁两侧的热量传递同样通过热传导进行的；根据热电类比,热传导在导热介质两侧温度差(△θ＝θ1-θ2)可以用以下公式计算：
△θ＝Q×RT，其中△θ为热传导在导热介质两侧温度差，θ1为导热介质一侧温度，θ2为导热介质另一侧温度，Q为绝缘层或油箱壁的热通量，RT为绝缘层或油箱壁热阻，下角标T代表绝缘层或油箱壁，β1为导热介质厚度，A为导热介质的横截面积；k为热导率。
定义q为单位面积的热通量即热流密度，则上式变为
油浸式电抗器油纸绝缘的热导率k是随温度变化的,在一般工程计算中，可取合理的定值计算。
(2)油浸式电抗器中的热辐射：
油浸式电抗器内部绕组和油、铁芯和油、铁芯和绕组、绕组之间以及绕组与油箱壁都存在温度差，都会发生热辐射，但相对于其间的热传导和热对流小很多，可忽略不计；油浸式电抗器油箱壁和散热器也会通过热辐射向周围空气散失热量，辐射物体表面的辐射率影响辐射程度，辐射的热转移斯特凡波尔兹曼定律表示为：
其中，PR单位时间内物体表面辐射出的总能量，下角标R仅用于区分并非为变量，斯特凡波尔兹曼常数η＝5.67×10-8w/m2·k4，E为表面发射率，AR为辐射表面积，Ts为辐射表面平均温度，Tα为周围空气温度，下角标s和下角标α表示仅用于区分并非为变量。
(3)油浸式电抗器中的热对流
油浸式电抗器内部绕组和油、铁芯和油、油与油箱壁以及油箱壁和外界空气之间的换热都是通过热对流进行的；油浸式电抗器油箱壁和散热器的热量转移到临近空气中，使周围的空气受热，空气是可压缩气体，受热体积会膨胀，产生运动。但是，加热的空气冷却后并不回到原来的状态，而是由新的同样密度的空气补充。
对流换热公式如下：
Q＝α(tw-tf)S
其中，Q为对流换热量，单位为W；tw和tf为壁面和流体平均温度，单位为℃，w表示油箱壁壁面，f表示流体，下角标w、f非变量，仅用于区分两个温度；S为对流换热面积，单位为㎡；α为对流换热系数，单位为W/㎡·℃。
变压器中有三种传热散热方式:热传导、热对流以及热辐射。在油浸式变压器中，对流和传导起着相当重要的作用，热辐射的作用相对的比较小。在强油循环的变压器中，热辐射的作用可以忽略不计；而自然油循环变压器中，传导起到主要作用，对流次之，辐射起到一定的作用。
D：确定温度场仿真计算的边界条件，分别为：
温度场仿真计算的相关背景为：
导热的基本定律及微分方程
为了研究物体的导热问题，本发明引入连续介质假定，即认为物体内部介质是连续的，运用连续函数描述物体的温度分布，物体内部各点的温度可以表述为时间和空间坐标的函数，即
T＝T(x,y,z,t)；
其中，T表示温度，x，y，z为物体温度的几何空间坐标，t为时间；
导热的基本定律为在单位时间内,通过一定截面的导热量和这个截面的温度梯度呈正比例，而热量的传递方向则与温度梯度方向相反,即傅里叶定律：
其中，q为热流密度；Φ为传导热流量；A为与传热方向垂直的横截面积；λ为导热系数；表示温度沿n方向的梯度，T为温度；n为热流传递方向。
固液耦合界面
在固液交界面上,固体和液体的热耦合是靠对流换热进行的，对流换热的基本计算公式一般采用牛顿冷却公式。在电抗器的铁芯和绕组与油流之间，及油箱壁和空气之间，流体温度低于固体温度，流体被加热：
q＝h(tw-tf)；
在油流与油箱壁之间，流体温度高于固体温度，固体被加热：
q＝h(tf-tw)；
其中，q为热流密度；h为表面对流换热系数；tw为壁面温度；tf为流体温度。
表面对流换热系数h的大小与换热过程的许多因素有关,它不仅取决于流体的性质以及换热表面的形状、大小与布置,而且还与流速有密切的关系,可用下式表示：
h＝f(vf,tw,tf,λ,c,ρ,μ,β2,x1,y1,z1,φ1)；
其中，vf为流体的速度；tw为壁面温度；tf为流体温度；λ、C和ρ分别为流体的导热系数、比热容和密度；μ为流体的黏度；β2为流体的体积膨胀系数；x1、y1、z1分别为直角坐标下被求点的坐标；φ1为固体表面形状系数。
由于对流传热系数是一个受多变量控制的复杂函数，很难准确的进行理论计算，所以一般在工程上多采用经验公式。
上述内容为温度场仿真计算中导热微分方程的确定，由于导热微分方程是描述导热过程共性的通用表达式，从数学角度来看，求解导热微分方程式可获得方程的通解。然而，每一个具体的导热过程总是在特定条件进行的。因此，为获得唯一的温度分布函数，还必须给出表达该过程特点的补充说明条件，即确定微分方程特解的边界条件。
确定温度场仿真计算的边界条件：
(1)第一类边界是已知任何时刻物体边界面的温度分布，即
T|S1＝T0；
其中，T为物体温度，S1为物体的边界面，T0在稳态热传导过程中表示给定的温度，在非稳态导热过程中表示温度值；
(2)第二类边界是已知任何时刻物体边界面上的热流密度分布，即
其中，λ为导热系数，T为物体温度，表示温度沿n方向的梯度，n为热流传递方向；S2为物体的边界面；q0为通过边界面S2的热流密度，对于稳态导热过程，q0为常量；对于非稳态过程，q0是随时间变化的量。
(3)第三类边界条件是已知物体与周围流体间的对流换热系数和流体温度，即
其中，λ为导热系数；T为物体温度；n为热流传递方向；S3为物体的边界面；、α为物质的散热系数，Tf为流体温度，T为物体温度。上述内容中，S1、S2和S3中的下角标1、2和3，以及Tf的下角标f，均是仅用于区分并非为变量。
E：计算正常工况下平波电抗器稳态温度场；
使用ANSYS有限元分析软件，利用上述步骤中确定的网格划分后的油浸式平波电抗器三维几何模型，输入材料物理参数和边界条件，然后在平波电抗器绕组上加载稳态电流，并在平波电抗器油箱外侧远场模型上加载外界环境温度，本实施例中外界环境温度为20度；外界环境温度加载后，在产热的同时考虑散热对温度的影响，根据步骤C中确定的热传导公式、热辐射公式和对流换热公式，计算出平波电抗器内部热量交换以及油箱壁的散热；最后通过ANSYS有限元分析软件计算出稳态温度场的结果，即平波电抗器三维模型内各个点的温度值，并以三维图像的形式显示出来。其中，稳态电流指油浸式平波电抗器正常运行时绕组上通过的电流。
F：将求得稳态温度场结果作为故障暂态模型的初始值，然后加载故障电流，采用分步长的方法，使用ANSYS有限元软件的瞬态分析功能，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。
本发明中，以典型换相失败故障电流作为故障模型中使用的故障电流。连续换相失败电流即假设多个故障时间点，在每个故障时间点加载一个典型故障电流。
换相失败故障电流的分析原理及确定过程：
由于关断角γ与触发超前角β3、换流变压器变比k、逆变侧直流电流Id、等效换相电抗Xc、逆变器交流侧母线电压有效值UL的相互关系为：
从γ的表达式可以得出，当故障导致Id急剧增大或者UL急剧下降时，会导致双极直流电压迅速大幅下降，故障初期直流电流急剧增加，产生过冲电流，过冲电流峰值主要受直流平波电抗器和换流变压器阻抗的限制，直流电流过冲到约为额定值的2倍；换相重叠角μ也因换相电压瞬时跌落而增大，逆变侧的控制系统还来不及将超前触发角β3加大，使得关断角γ从先前的稳定值瞬时跌落到零，导致换相失败的发生。为抑制急剧增大的短路电流，VDCL环节起作用，整流侧提高触发角α以减小直流电流减小；故障消除后，VDCL确保直流电流逐渐爬升至故障前水平。
在ANSYS有限元软件中可以以表格的形式加载不同电流，这里将电流-时间图以二维数组的形式制成表格，作为计算中加载的电流，便能在计算中加载故障电流。
由于换相是一个几分钟左右的暂态过程，分析中需要使用ANSYS有限元软件中求解模块的暂态分析功能，采用分步长的方法，即依据以故障后数秒或者数毫秒作为一个步长，在ANSYS有限元软件中只需设置好每个步长的长度，并设置一个计算时间终点，就可以分步计算出不同时间点的温度场分布，一般步长越小计算精度越高，但同时所需要的计算时间也更长，可在故障发生时依据故障波形以几十到一百毫秒为步长，非故障时步长可设置为0.5-5秒不等。
利用ANSYS，在设置好加载的电流跟步长后就可以进行故障后暂态温度场分布的计算，ANSYS有限元软件中的后处理模块能保存并查看每步中平波电抗器三维模型内各点的温度值。
Claims (6)
1.一种基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于，包括以下步骤：
A：基于ANSYS有限元软件，采用实体建模的方式建立油浸式平波电抗器三维几何模型；然后对油浸式平波电抗器三维几何模型进行材料物理参数确定、单元类型确定以及网格划分；
B：对油浸式平波电抗器进行产热分析，分别确定油浸式平波电抗器中绕组电流产生的铜损、铁芯产生的铁损和由于漏磁通在电抗器结构件中产生的杂散损耗；绕组电流产生的铜损p的计算式为：I1为绕组的相电流；r1为绕组的总电阻；铁芯产生的铁损P0的计算式为：P0＝K0GPc，K0为铁芯损耗的工艺系数，G为铁芯的重量，Pc为单位重量铁耗，P1/50为铁耗系数，为Bm＝1T且f＝50Hz时每千克硅钢片的铁耗，β为频率指数，Bm为最大磁通密度，B指代磁通密度，下角标m指代max；F为频率；漏磁通在电抗器结构件中产生的杂散损耗在总损耗中所占的比重较小，在此不再考虑；
C：对油浸式平波电抗器进行散热分析，分别确定油浸式电抗器中的热传导、热辐射和热对流；
D：确定温度场仿真计算的边界条件；
E：计算正常工况下平波电抗器稳态温度场；
F：将求得稳态温度场结果作为故障暂态模型的初始值，然后加载故障电流，采用分步长的方法，使用ANSYS有限元软件的瞬态分析功能，计算出故障后温度场随时间的分布，完成基于平波电抗器多次换向失败时的暂态温度场仿真。
2.根据权利要求1所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于：所述的步骤A中，油浸式平波电抗器三维几何模型包括绕组、铁芯、油箱、油、套管和冷却系统；在对油浸式平波电抗器三维几何模型中材料物理参数进行确定时，依次输入绕组、铁芯、油箱和油的电阻率及热传导系数，并输入铁芯非线性磁导率。
3.根据权利要求2所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于：所述的步骤A中，在ANSYS有限元软件中，以铁芯的B-H曲线的形式输入非线性磁导率，其中铁芯的B-H曲线以二维数组的形式拟合。
4.根据权利要求1所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于：所述的步骤C中，对油浸式平波电抗器进行散热分析，分别确定油浸式电抗器中的热传导、热辐射和热对流；
(1)油浸式电抗器中的热传导：
热传导在导热介质两侧温度差可以用以下公式计算：
Δθ＝θ1-θ2＝Q×RT，
其中Δθ为热传导在导热介质两侧温度差，θ1为导热介质一侧温度，θ2为导热介质另一侧温度，Q为绝缘层或油箱壁的热通量，RT为绝缘层或油箱壁热阻，下角标T代表绝缘层或油箱壁，β1为导热介质厚度，A为导热介质的横截面积；k为热导率；
定义q为单位面积的热通量即热流密度，则上式变为
<mrow>
<mi>&Delta;</mi>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>q</mi>
<mo>&times;</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>k</mi>
</mfrac>
<mo>;</mo>
</mrow>
油浸式电抗器油纸绝缘的热导率k为定值；
(2)油浸式电抗器中的热辐射：
辐射的热转移斯特凡波尔兹曼定律表示为：
<mrow>
<msub>
<mi>P</mi>
<mi>R</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&eta;EA</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mi>S</mi>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>&alpha;</mi>
<mn>4</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，PR单位时间内物体表面辐射出的总能量，下角标R仅用于区分并非为变量，斯特凡波尔兹曼常数η＝5.67×10-8w/m2·k4，E为表面发射率，AR为辐射表面积，Ts为辐射表面平均温度，Tα为周围空气温度，下角标s和下角标α表示仅用于区分并非为变量；
(3)油浸式电抗器中的热对流
对流换热公式如下：
Q＝α(tw-tf)S；
其中，Q为对流换热量，单位为W；tw和tf为壁面和流体平均温度，单位为℃，w表示油箱壁壁面，f表示流体，下角标w、f非变量，仅用于区分两个温度；S为对流换热面积，单位为㎡；α为对流换热系数，单位为W/㎡·℃。
5.根据权利要求1所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于：所述的步骤D中，温度场仿真计算的边界条件分别为：
(1)第一类边界是已知任何时刻物体边界面的温度分布，即
T|S1＝T0；
其中，T为物体温度，S1为物体的边界面，T0在稳态热传导过程中表示给定的温度，在非稳态导热过程中表示温度值；
(2)第二类边界是已知任何时刻物体边界面上的热流密度分布，即
<mrow>
<mo>-</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mi>n</mi>
<msub>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
其中，λ为导热系数，T为物体温度，表示温度沿n方向的梯度，n为热流传递方向；S2为物体的边界面；q0为通过边界面S2的热流密度，对于稳态导热过程，q0为常量；对于非稳态过程，q0是随时间变化的量；
(3)第三类边界条件是已知物体与周围流体间的对流换热系数和流体温度，即
<mrow>
<mo>-</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mi>n</mi>
<msub>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，λ为导热系数；T为物体温度；表示温度沿n方向的梯度，n为热流传递方向；S3为物体的边界面；α为物质的散热系数，Tf为流体温度，T为物体温度。
6.根据权利要求5所述的基于平波电抗器多次换向失败时的暂态温度场仿真方法，其特征在于：所述的步骤E中，使用ANSYS有限元分析软件，利用上述步骤中确定的网格划分后的油浸式平波电抗器三维几何模型，输入材料物理参数和边界条件，然后在平波电抗器绕组上加载稳态电流，并在平波电抗器油箱外侧远场模型上加载外界环境温度；外界环境温度加载后，在产热的同时考虑散热对温度的影响，根据步骤C中确定的热传导公式、热辐射公式和对流换热公式，计算出平波电抗器内部热量交换以及油箱壁的散热；最后通过ANSYS有限元分析软件计算出稳态温度场的结果，即平波电抗器三维模型内各个点的温度值，并以三维图像的形式显示出来；其中，稳态电流指油浸式平波电抗器正常运行时绕组上通过的电流。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510232180.0A CN104834781B (zh) | 2015-05-08 | 2015-05-08 | 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510232180.0A CN104834781B (zh) | 2015-05-08 | 2015-05-08 | 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104834781A CN104834781A (zh) | 2015-08-12 |
| CN104834781B true CN104834781B (zh) | 2018-04-27 |
Family
ID=53812665
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510232180.0A Active CN104834781B (zh) | 2015-05-08 | 2015-05-08 | 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104834781B (zh) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106021663B (zh) * | 2016-05-09 | 2019-11-29 | 广东电网有限责任公司电力科学研究院 | 一种变压器重瓦斯保护的分析方法和装置 |
| CN106529096A (zh) * | 2016-12-19 | 2017-03-22 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种整体油箱温度场分析方法 |
| CN107273628A (zh) * | 2017-06-26 | 2017-10-20 | 浙江群力电气有限公司 | 一种地下配电室的温度场计算方法及系统 |
| CN107391835B (zh) * | 2017-07-18 | 2020-09-01 | 广东电科院能源技术有限责任公司 | 一种高压变压器数值模拟网格划分方法及装置 |
| CN107451381B (zh) * | 2017-09-18 | 2020-12-18 | 海南电网有限责任公司电力科学研究院 | 一种基于磁-热-流耦合的变压器热稳定计算方法 |
| CN108280278B (zh) * | 2018-01-12 | 2024-07-02 | 中国电力科学研究院有限公司 | 一种高压直流平抗瓦斯动态特性的仿真方法及系统 |
| CN108683164B (zh) * | 2018-04-17 | 2020-10-30 | 全球能源互联网研究院有限公司 | 饱和电抗器及其参数整定方法、装置及仿真设备 |
| CN109086514B (zh) * | 2018-07-25 | 2023-07-25 | 国网江苏省电力有限公司检修分公司 | 一种基于多物理场耦合的电抗器温度场计算方法 |
| CN109085452B (zh) * | 2018-09-17 | 2021-07-06 | 广东电网有限责任公司 | 短路故障检测装置 |
| CN109508502B (zh) * | 2018-11-22 | 2023-06-30 | 国家电网有限公司 | 基于电磁-温湿度耦合的三维开关柜温湿度分布计算方法 |
| CN111368385B (zh) * | 2018-12-25 | 2023-06-27 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电热耦合场中套管运行特性仿真计算方法 |
| CN110765685A (zh) * | 2019-10-18 | 2020-02-07 | 南方电网科学研究院有限责任公司 | 一种电抗器多物理场耦合的仿真方法、装置及存储介质 |
| CN110889209B (zh) * | 2019-11-18 | 2023-04-28 | 中国北方车辆研究所 | 一种润滑油加温仿真方法 |
| CN111597748B (zh) * | 2020-05-22 | 2024-05-17 | 中国矿业大学 | 一种基于gil热特性实现故障判定的方法 |
| CN111985075B (zh) * | 2020-07-03 | 2024-03-26 | 国网山东省电力公司电力科学研究院 | 适用于氧化锌避雷器的温度分布计算方法及系统 |
| CN112380753B (zh) * | 2020-11-25 | 2022-04-22 | 三峡大学 | 一种提高油浸式电抗器绕组散热性能的方法 |
| EP4266067A4 (en) * | 2020-12-17 | 2024-02-14 | Nissin Electric Co., Ltd. | METHOD FOR TEMPERATURE INCREASE TESTING OF REACTORS |
| CN112926240B (zh) * | 2021-02-01 | 2022-07-01 | 中国人民解放军海军工程大学 | 短时工作制大容量水冷电抗器设计方法 |
| CN113128025B (zh) * | 2021-03-19 | 2022-09-16 | 广西电网有限责任公司电力科学研究院 | 一种变压器绕组流体温度场仿真模型的优化方法 |
| CN113237920B (zh) * | 2021-05-17 | 2022-04-22 | 西南交通大学 | 一种特高压换流变压器阀侧套管故障热源检测方法 |
| CN114896920B (zh) * | 2022-07-12 | 2022-09-30 | 中南民族大学 | 一种三维堆叠芯片热仿真模型建立及热点温度预测方法 |
| CN115659659A (zh) * | 2022-10-31 | 2023-01-31 | 云南电网有限责任公司电力科学研究院 | 一种计算变压器自然对流散热下绕组温度场的方法及系统 |
| CN117725780B (zh) * | 2023-12-01 | 2024-07-30 | 国网宁夏电力有限公司电力科学研究院 | 一种平波电抗器的绝缘故障概率评估方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819239A (zh) * | 2010-04-06 | 2010-09-01 | 嘉兴学院 | 基于三维温度场快速构建的变压器故障诊断系统 |
| CN103728539A (zh) * | 2014-01-23 | 2014-04-16 | 华北电力大学（保定） | 基于分布式光纤测温的电缆电气故障仿真分析方法 |
| CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
-
2015
- 2015-05-08 CN CN201510232180.0A patent/CN104834781B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819239A (zh) * | 2010-04-06 | 2010-09-01 | 嘉兴学院 | 基于三维温度场快速构建的变压器故障诊断系统 |
| CN103728539A (zh) * | 2014-01-23 | 2014-04-16 | 华北电力大学（保定） | 基于分布式光纤测温的电缆电气故障仿真分析方法 |
| CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| "电子式电流互感器暂态传变延时测试技术研究";赵勇 等;《电力系统保护与控制》;20140901;第42卷(第17期);第125-130页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104834781A (zh) | 2015-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104834781B (zh) | 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 | |
| Lv et al. | A novel thermal management structure using serpentine phase change material coupled with forced air convection for cylindrical battery modules | |
| Radwan et al. | Uniform cooling for concentrator photovoltaic cells and electronic chips by forced convective boiling in 3D-printed monolithic double-layer microchannel heat sink | |
| Zhao et al. | Thermal performance of mini-channel liquid cooled cylinder based battery thermal management for cylindrical lithium-ion power battery | |
| Vijayan et al. | Steady state and stability characteristics of single-phase natural circulation in a rectangular loop with different heater and cooler orientations | |
| CN109086514B (zh) | 一种基于多物理场耦合的电抗器温度场计算方法 | |
| CN112115628B (zh) | 基于油浸式变压器温度场分布计算的热点温度检测方法 | |
| CN105045966A (zh) | 一种变压器内部温度场的混合计算方法 | |
| CN107451381B (zh) | 一种基于磁-热-流耦合的变压器热稳定计算方法 | |
| CN111753449A (zh) | 一种获取电力变压器不同工作情况下热点温度的仿真方法 | |
| Eslami-Nejad et al. | Heat transfer in double U-tube boreholes with two independent circuits | |
| CN107871037A (zh) | 一种35kV干式电抗器内部温度估算方法 | |
| Lamaison et al. | Two-phase flow control of electronics cooling with pseudo-CPUs in parallel flow circuits: dynamic modeling and experimental evaluation | |
| Liu et al. | Pressure drop and void fraction in steam-water two-phase flow at high pressure | |
| Wan | Thermal performance of heat pipe array in battery thermal management | |
| Wang et al. | Heat transfer characteristic of water at near critical pressure in circumferentially non-uniformly heated vertical tubes | |
| Jiang et al. | An immersion flow boiling heat dissipation strategy for efficient battery thermal management in non-steady conditions | |
| Santisteban et al. | Thermal modelling of electrical insulation system in power transformers | |
| Ebenezer et al. | Determination of temperature distribution and reduction of life of a dry type three phase transformer | |
| CN118228669B (zh) | 油浸式配电变压器绕组匝间绝缘劣化的温度场分析方法 | |
| Wan | A numerical investigation of the thermal performances of an array of heat pipes for battery thermal management | |
| Mirahmadi et al. | Study on preheating techniques for start-up of tubular solid oxide fuel cells | |
| Zeng et al. | Transient temperature field simulation based on ANSYS for oil-immersed smoothing reactor when repeatedly Commutation Failure | |
| Dong et al. | Multi-Physical Field Coupled Method for Temperature Rise of Winding in Oil-immersed Power Transformer | |
| Mariani et al. | Qualitative thermal characterization and cooling of lithium batteries for electric vehicles |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yang Haijing Inventor after: Gong Liang Inventor after: Han Wei Inventor after: Shi Guang Inventor after: Li Bin Inventor after: Chen Kai Inventor after: Li Chaohui Inventor after: Dang Yiqi Inventor after: Ma Rui Inventor after: Zeng Kai Inventor before: Yang Haijing Inventor before: Shi Guang Inventor before: Han Wei Inventor before: Li Bin Inventor before: Ma Rui Inventor before: Zeng Kai Inventor before: Gong Liang |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |