留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

PSM高压电源干式变压器的热分析计算

夏于洋 李青 毛晓惠

引用本文:
Citation:

PSM高压电源干式变压器的热分析计算

    作者简介: 夏于洋(1992—),男,硕士,从事波加热辅助高压电源技术工作;xiayy@swip.ac.cn.
  • 基金项目: 国家重点研发计划项目(2017YFE0300504)
  • 中图分类号: TM4

Thermal analysis calculation of dry-type transformer in PSM high voltage power supply

  • CLC number: TM4

  • 摘要: 多绕组干式整流变压器作为托卡马克装置加热系统中高压电源的一个重要设备,其性能优良与否直接关系到高压电源的输出品质。开展变压器的功率损耗与散热分析研究,减少变压器因温度上升造成的性能影响,对保证变压器的良好工作状态是十分重要的。通过MATLAB对干空气下的比热容、传热系数等热物理性质与温度的关系进行拟合分析,得到相关的关系方程,对变压器的对流传热与热辐射进行分析计算,得到仅在空气自然对流和热辐射的情况下,不能使变压器的温度控制在满足性能的温度范围之内。在强迫空气对流的情况下,实现了变压器良好的散热。进一步利用ANSYS对变压器的温度场分布进行分析,利用温度场的分布趋势图观察强迫风冷相对于自然冷却的优点,强迫风冷降低了整体的温度,使热量更快地散发到周围大气中,减少了变压器受高温的危害。
  • 图 1  简化的变压器图

    Fig. 1  Simplified transformer diagram

    图 2  干空气的空气密度、比热容、动力粘度、导热系数等与温度之间的关系

    Fig. 2  Relationship between air density, specific heat capacity, dynamic viscosity, thermal conductivity, and temperature of dry air

    图 3  强迫空气对流出口温度与空气流速的关系

    Fig. 3  Exit temperature changes with wind speed

    图 4  风冷通道中温度场分布等值线图

    Fig. 4  Contour map of temperature field distribution in air-cooled channels

    图 5  风冷通道中Y方向等值线图

    Fig. 5  Y-direction contour map in the air-cooled channel

    图 6  变压器自然冷却时温度变化图

    Fig. 6  Temperature change diagram of transformer during natural cooling

    表 1  变压器的基本参数

    Table 1  Basic parameters of the transformer

    R/mr/m${P_{{\rm{Cu}}}}$/W${P_{{\rm{Fe}}}}$/WH/m${T_{{\rm{fo}}}}$/${}^ \circ {\rm{C}}$D/mA/m2
    0.980.7439 82310 7701.8250.207.49
    下载: 导出CSV

    表 2  干空气的热物理性质

    Table 2  Thermophysical properties of dry air($p = 100 \;{\rm{kPa}}$

    $T/{}^\circ {\rm{C}}$$\rho /\left( {{\rm{kg}} \cdot {{\rm{m}}^{\rm{3}}}} \right)$${c_p}/\left( {{\rm{J}} \cdot {\rm{k}}{{\rm{g}}^{{\rm{ - 1}}}} \cdot {{\rm{K}}^{{\rm{ - 1}}}}} \right)$$\mu /\left( {{\rm{1}}{{\rm{0}}^{{\rm{ - 6}}}}{\rm{Pa}} \cdot {\rm{s}}} \right)$$k/\left( {{\rm{mW}} \cdot {{\rm{m}}^{ - 1}} \cdot {{\rm{K}}^{ - 1}}} \right)$
    201.1641.01318.2425.24
    401.0921.01319.2226.52
    601.0251.01720.4028.03
    800.9681.02220.992.931
    1000.9161.02221.7730.70
    1200.8701.02622.7531.98
    1400.8271.02623.5433.26
    1600.7891.03024.1234.42
    2000.7321.03425.8936.98
    2500.6531.04327.9539.77
    3000.5961.04729.7142.91
    下载: 导出CSV
  • [1] 孙雪静. 超高速永磁同步电机的热设计与大电流检测技术[D]. 西安: 西安电子科技大学, 2012: 43-46.

    Sun Xuejing. Thermal design and high current detection technology for ultra-high speed permanent magnet synchronous motor. Xi'an: Xidian University, 2012: 43-46
    [2] 庞营. 变压器发热与散热的数学研究[J]. 超星, 2014, 1(1):53-54. (Pang Ying. Mathematical study on heating and heat dissipation of transformers. Superstar, 2014, 1(1): 53-54
    [3] 张鲁, 李永光, 高向宗, 等. 变压器冷却技术的发展[J]. 能源与节能, 2016, 128(5):109-110. (Zhang Lu, Li Yongguang, Gao Xiangzong, et al. Development of transformer cooling technology. Energy and Energy Conservation, 2016, 128(5): 109-110 doi: 10.3969/j.issn.2095-0802.2016.05.054
    [4] 金立军. 变压器室对流换热的分析[J]. 电力建设, 2000(8):19-20. (Jin Lijun. Analysis of convective heat transfer in transformer room. Electric Power Construction, 2000(8): 19-20 doi: 10.3969/j.issn.1000-7229.2000.08.007
    [5] 许进峰. ANSYS Workbench 15.0完全自学一本通[M]. 北京: 电子工业出版社, 2014: 426-427.

    Xu Jinfeng. ANSYS Workbench 15.0 complete self-study guide. Beijing: Publishing House of Electronics Industry, 2014: 426-427
    [6] 苏玎, 于学莲. 110 kV级变压器冷却方式的选择[J]. 科技风, 2015(2):11. (Su Ding, Yu Xuelian. Selection of 110 kV transformer cooling method. Science and Technology, 2015(2): 11 doi: 10.3969/j.issn.1671-7341.2015.02.009
    [7] 星文彤, 李志尊, 胡仁喜, 等. ANSYS13.0热力学有限元分析从入门到精通[M]. 北京: 机械工业出版社, 2011: 11-16.

    Xing Wentong, Li Zhizun, Hu Renxi, et al. ANSYS13.0 thermodynamic finite element analysis from entry to proficiency. Beijing: China Machine Press, 2011: 11-16
    [8] 郭振华. 大型变压器绕组涡流损耗及热点温度的研究[D]. 河北: 河北工业大学, 2013: 5-8.

    Guo Zhenhua. Study on eddy current loss and hot spot temperature of large transformer windings. Hebei: Hebei University of Technology, 2013: 5-8
    [9] 张锋. 基于ANSYS的DC/DC电源模块热分析和热设计研究[D]. 重庆: 重庆大学, 2008: 27-29.

    Zhang Feng. Thermal analysis and thermal design of DC/DC power module based on ANSYS. Chongqing: Chongqing University, 2008: 27-29
    [10] 邢雅, 樊博, 康亚丽, 等. 树脂绝缘干式变压器内部温度场的计算[J]. 研究与设计, 2016, 40(1):34-38. (Xing Ya, Fan Bo, Kang Yali, et al. Calculation of internal temperature field of resin insulated dry type transformer. Research & Design, 2016, 40(1): 34-38
    [11] 谢龙汉. ANSYS FLOTRAN流体与热分析[M]. 北京: 电子工业出版社, 2012: 413-456.

    Xie Longhan. ANSYS FLOTRAN fluid and thermal analysis. Beijing: Publishing House of Electronics Industry, 2012: 413-456
  • [1] 郭静, 张国军, 王续博, 张文栋, . MEMS仿生矢量水听器微结构的有限元分析[J]. 强激光与粒子束, 2015, 27(02): 024125. doi: 10.11884/HPLPB201527.024125
    [2] 梁秀强, 袁杰红, 周仕明. 离子推力器栅极组件温度场仿真分析及试验研究[J]. 强激光与粒子束, 2018, 30(11): 114001. doi: 10.11884/HPLPB201830.180208
    [3] 姚列明, 杨中海, 李斌, 黄桃. 行波管收集极的热分析[J]. 强激光与粒子束, 2006, 18(01).
    [4] 牛婧杨, 王丽, 罗勇, 蒋伟. 回旋行波管电子枪阴极热分析[J]. 强激光与粒子束, 2013, 25(02): 446-450. doi: 10.3788/HPLPB20132502.0446
    [5] 焦俊科, 王新兵. CO2激光作用下运动石英玻璃的温度分布[J]. 强激光与粒子束, 2007, 19(01).
    [6] 徐斌, 刘景全, 江水东, 黄永华, 杨斌, 陈翔, 杨春生. 冷冻靶中硅冷却臂的温度场和力学分析[J]. 强激光与粒子束, 2015, 27(06): 062009. doi: 10.11884/HPLPB201527.062009
    [7] 申永明, 张耀飞, 侯天晋, 陈德章, 卿光弼. 固体热容激光器与热稳态固体激光器温度场和应力场的比较研究[J]. 强激光与粒子束, 2006, 18(11).
    [8] 雷文强, 杨中海, 肖礼. 高频腔本征值问题的三维数值模拟[J]. 强激光与粒子束, 2003, 15(02).
    [9] 陈辉, 王丽, 罗勇, 唐勇. Q波段回旋行波管新型盒型输出窗的设计[J]. 强激光与粒子束, 2015, 27(01): 013002. doi: 10.11884/HPLPB201527.013002
    [10] 陈庆红, 徐世珍, 黄进, 孙来喜, 郭德成, 程强, 刘红婕, 蒋晓东. 预应力对熔石英损伤的影响[J]. 强激光与粒子束, 2013, 25(10): 2531-2535. doi: 10.3788/HPLPB20132510.2531
    [11] 岱钦, 李新忠, 2, 王希军. LDA侧面泵浦Nd:YAG激光器的热效应分析[J]. 强激光与粒子束, 2007, 19(02).
    [12] 王伟平, 王翔. 激光辐照下金属/炸药温度场的计算[J]. 强激光与粒子束, 1999, 11(02).
    [13] 魏继锋, 卢飞, 蒋志雄, 常艳, 周文超, 周山, 胡晓阳. 阶梯反射锥型高能激光能量计温度场控制技术[J]. 强激光与粒子束, 2015, 27(12): 121001. doi: 10.11884/HPLPB201527.121001
    [14] 谭玉凤, 王继红, 任戈, 朱福音. 大口径轻量化主镜的温度场等效模型理论计算[J]. 强激光与粒子束, 2017, 29(06): 061001. doi: 10.11884/HPLPB201729.160555
    [15] 王伟平, 刘绪发, 孙承纬. 激光时空分布对钢靶温度场及热软化的影响[J]. 强激光与粒子束, 1997, 09(02).
    [16] 钟明, 罗大为, 朱柞金, 胡汉平, 何立群, 程曙霞. 激光辐照金属/炸药复合介质温度场的数值模拟[J]. 强激光与粒子束, 2000, 12(02).
    [17] 王艳茹, 李斌成, 刘明强. 激光照射有限尺寸高反射光学元件的温度场[J]. 强激光与粒子束, 2010, 22(02).
    [18] 杨光, 丁林林, 钦兰云, 卞宏友, 王维. TA15钛合金激光沉积温度场数值模拟与检测[J]. 强激光与粒子束, 2014, 26(11): 119002. doi: 10.11884/HPLPB201426.119002
    [19] 高晓锋, 韦建军, 周秀文, 安辛友, 吴卫东. 玻璃包覆纯铜丝快速冷却温度场数值模拟[J]. 强激光与粒子束, 2013, 25(08): 1970-1974. doi: 10.3788/HPLPB20132508.1970
    [20] 邵珠强, 胡增荣, 郭绍雄, 倪阳阳, 李悦, 张瑶, 陈长军, 王晓南. 激光烧结石墨烯铜复合材料温度场有限元模拟[J]. 强激光与粒子束, 2018, 30(03): 039001. doi: 10.11884/HPLPB201830.170366
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  73
  • HTML全文浏览量:  66
  • PDF下载量:  1
出版历程
  • 收稿日期:  2019-08-12
  • 录用日期:  2019-09-16
  • 网络出版日期:  2019-11-09

PSM高压电源干式变压器的热分析计算

    作者简介: 夏于洋(1992—),男,硕士,从事波加热辅助高压电源技术工作;xiayy@swip.ac.cn
  • 核工业西南物理研究院,成都 610225

摘要: 多绕组干式整流变压器作为托卡马克装置加热系统中高压电源的一个重要设备,其性能优良与否直接关系到高压电源的输出品质。开展变压器的功率损耗与散热分析研究,减少变压器因温度上升造成的性能影响,对保证变压器的良好工作状态是十分重要的。通过MATLAB对干空气下的比热容、传热系数等热物理性质与温度的关系进行拟合分析,得到相关的关系方程,对变压器的对流传热与热辐射进行分析计算,得到仅在空气自然对流和热辐射的情况下,不能使变压器的温度控制在满足性能的温度范围之内。在强迫空气对流的情况下,实现了变压器良好的散热。进一步利用ANSYS对变压器的温度场分布进行分析,利用温度场的分布趋势图观察强迫风冷相对于自然冷却的优点,强迫风冷降低了整体的温度,使热量更快地散发到周围大气中,减少了变压器受高温的危害。

English Abstract

  • 中国聚变工程实验堆(CFETR)是中国自主设计的下一代超导托卡马克。该装置中辅助加热系统包括NBI(中性束注入加热)、LHCD(低混杂波电流驱动与加热)、ICRF(离子回旋共振加热)、ECRH(电子回旋共振加热)。CFETR辅助加热系统所需的高压电源计划采用统一的拓扑结构设计,采用直流电源构成基本模块单元,通过串联迭加达到所需高压。电压调节控制采用直接级联与PWM调节相结合的方式,即采用PSM高压电源技术。

    LHCD高压电源系统是辅助波加热单套电源系统中输出电压、加热功率最高的电源系统,其主要组成单元为交流配电单元、变压器单元、直流模块单元、控制系统单元、保护单元等。多绕组干式变压器作为PSM高压电源系统中的一个重要设备,其变压器的性能对高压电源的输出品质影响较大,对此研究变压器在各种工作条件下对其自身的影响很有必要。温度对于电气设备的性能影响很大,一旦温度过高,变压器的负载损耗会增加,且长时间处于高温状态,变压器组件的寿命会缩短,其绝缘性能降低,影响变压器的运行安全和性能。所以有必要对变压器的发热与散热进行分析。变压器是一种不均匀的发热体,其内部热量流动和传递的过程是非常复杂的,这些都增加了变压器的发热计算的难度。变压器产生的各种损耗均以热量的形式从热源内部传递到热源的表面,然后通过对流和热辐射的形式传递到周围空间中[1-3]。因此需要讨论三种传热方式,传导、对流和辐射。在进行温度分析时需对每一实体的材料属性进行研究,包括热导率、辐射系数、密度、吸收率、对流换热系数、比热容等。

    • PSM高压电源系统中的干式变压器是一个常规的圆柱式变压器。在进行PSM干式变压器的热分析是采用简化的变压器模型,通过对流传热原理对变压器的结构温度与散热进行分析估计,简化的变压器模型如图1所示。图1(a)表示的是变压器的简化柱状图,其外面的红色部分表示的是风道,内部绿色表示的是变压器铁芯和绕组。图1(b)中的${T_{{\rm{ai}}}}$${T_{{\rm{ao}}}}$分别表示空气流进风道的进口温度和出口温度,${T_{\rm{o}}}$表示外壁温度,一般认为外壁温度和周围环境温度一样,在本文中,选择外壁温度为25 ℃。${T_{{\rm{fo}}}}$表示变压器的外壁温度,即铜导体与铁磁材料等各种变压器损耗的能量转化为热能进行散发的表面温度[4]

      图  1  简化的变压器图

      Figure 1.  Simplified transformer diagram

      LHCD辅助电源系统输出参数为−80 kV/100 A。辅助低杂波高压电源系统中,变压器的额定容量为6 000 kV,运行时间3 600 s,变压器的短路阻抗不大于8%。变压器的其他基本参数如表1所示,其中R表示变压器外半径,r表示变压器内半径,${P_{{\rm{Cu}}}}$为变压器铜损,${P_{{\rm{Fe}}}}$为变压器铁损,H为变压器的高度,${T_{{\rm{fo}}}}$为变压器外壁温度,D为风道宽度,A为变压器散热面积。干式变压器的损耗为PSM高压电源系统工作在输出直流电压−80 kV,电流100 A时其稳态过程中的损耗。

      表 1  变压器的基本参数

      Table 1.  Basic parameters of the transformer

      R/mr/m${P_{{\rm{Cu}}}}$/W${P_{{\rm{Fe}}}}$/WH/m${T_{{\rm{fo}}}}$/${}^ \circ {\rm{C}}$D/mA/m2
      0.980.7439 82310 7701.8250.207.49
    • 对流是由于温度不同的各部分流体之间发生相对运动所引起的热量传递方式,与高温部件相接触的空气部分,由于吸收传递出来的热量而膨胀,密度降低,从而向上运动,密度较大的冷空气将下降代替原来因受热而上升的热空气而引发空气对流现象。空气的对流情况有自然对流和强迫对流[5]

      热对流的换热能量公式

      $Q = hA\left( {{T_1} - {T_2}} \right)$

      (1)

      式中:$h$为平均换热系数,${\rm{W/}}\left( {{{\rm{m}}^{\rm{2}}} \cdot {\rm{K}}} \right)$$A$为换热表面积,${{\rm{m}}^{\rm{2}}}$$A = {\text{π}}H\left( {R + r} \right)$${T_1}$为固体表面温度,${\rm{K}}$${T_2}$为周围流体温度,${\rm{K}}$

      自然对流散热是最普遍发生的现象,在本文中,变压器周围空间空气产生的自然对流的换热系数相对于其他环境来说,其值较小,因几台变压器在实际运行之中都处于室内,这一条件限制了空气的自然对流[6-7]。所以取自然对流的换热系数$\tilde h = 10\;{\rm{W/}}\left( {{{\rm{m}}^{\rm{2}}} \cdot {\rm{K}}} \right)$

      平均换热系数

      $h = \frac{{{N_{u{\rm{L}}}}k}}{L}$

      (2)

      式中:${N_ {u{\rm{L}}}}$为平均努塞尔数;$k$为热传导系数,${\rm{W/}}\left( {{{\rm{m}}^{\rm{2}}} \cdot {\rm{K}}} \right)$$L$为变压器的竖直壁高度,本文采用$L = H$

      由干空气下的空气导热系数与温度的关系可以近似的表达为

      $100k = 2.319 \times {10^{ - 6}} \times {T^2} + 0.007 \times T + 2.385$

      (3)

      而平均换热系数与导热系数成正比,依据相关环境条件可以对强迫对流中的换热系数进行如下近似

      ${h_{{\rm{force}}}} = 0.125 \times T + 25$

      (4)
    • 热辐射是指物体发射电磁能,并被其他物体吸收转变为热的热量交换过程。热辐射过程中不需要任何介质,系统中的每个物体同时吸收和辐射热量。它们之间的热量传递可以用斯蒂芬-玻耳兹曼方程来计算

      ${Q_{\rm{r}}} = \varepsilon \sigma A\left( {{T_3} - {T_4}} \right)$

      (5)

      式中:$\varepsilon $为辐射率(黑度),本文取$\varepsilon = 0.7$$\sigma $为黑体辐射常数,$\sigma = 5.67 \times {10^{ - 8}}\;{\rm{W/}}\left( {{{\rm{m}}^{\rm{2}}} \cdot {\rm{K}}} \right)$${T_3}$为辐射面1的绝对温度,${\rm{K}}$${T_4}$为辐射面2的绝对温度,${\rm{K}}$

    • 对于干式变压器采用的是在自然条件下的干空气进行散热,因此把握干空气的热物理性质对PSM干式变压器的热分析处理将会具有十分重要的意义。干空气的基本参数见表2,其中$T$为温度,$ \; \rho $为空气密度,${c_p}$为标准大气压下的比热容,$\mu $为动力粘度,$k$为导热系数。

      表 2  干空气的热物理性质

      Table 2.  Thermophysical properties of dry air($p = 100 \;{\rm{kPa}}$

      $T/{}^\circ {\rm{C}}$$\rho /\left( {{\rm{kg}} \cdot {{\rm{m}}^{\rm{3}}}} \right)$${c_p}/\left( {{\rm{J}} \cdot {\rm{k}}{{\rm{g}}^{{\rm{ - 1}}}} \cdot {{\rm{K}}^{{\rm{ - 1}}}}} \right)$$\mu /\left( {{\rm{1}}{{\rm{0}}^{{\rm{ - 6}}}}{\rm{Pa}} \cdot {\rm{s}}} \right)$$k/\left( {{\rm{mW}} \cdot {{\rm{m}}^{ - 1}} \cdot {{\rm{K}}^{ - 1}}} \right)$
      201.1641.01318.2425.24
      401.0921.01319.2226.52
      601.0251.01720.4028.03
      800.9681.02220.992.931
      1000.9161.02221.7730.70
      1200.8701.02622.7531.98
      1400.8271.02623.5433.26
      1600.7891.03024.1234.42
      2000.7321.03425.8936.98
      2500.6531.04327.9539.77
      3000.5961.04729.7142.91

      根据干空气的热物理性质可以得到其空气密度、比热容、动力粘度、导热系数等与温度之间的关系。利用MATLAB的曲线拟合进行分析,可以得到这四个物理性质的大致曲线方程。干空气的空气密度、比热容、动力粘度、导热系数等与温度之间的关系图与拟合曲线如图2所示。

      图  2  干空气的空气密度、比热容、动力粘度、导热系数等与温度之间的关系

      Figure 2.  Relationship between air density, specific heat capacity, dynamic viscosity, thermal conductivity, and temperature of dry air

      在进行热设计分析的时候,由于温度的不同将引起空气的参数值的不同,利用图2中拟合得到的各参数与温度的关系式,代入热分析的计算公式中,可以简化计算过程,得到更加精确的结果。

    • 当采用空气自然散热时,变压器达到稳定温度时,对空气温升进行计算。

      绕组的热流密度

      ${q_{\rm{Q}}} = \frac{{\Delta {P_{{\rm{Cu}}}}}}{{3{F_{\rm{Q}}}}} = \frac{{39\;832}}{{3 \times 5.24}}\;{\rm{W/}}{{\rm{m}}^{\rm{2}}} = 2\;534\;{\rm{W/}}{{\rm{m}}^{\rm{2}}}$

      (6)

      绕组对周围空气的温升

      ${\theta _{\rm{Q}}} = 0.36q_{\rm{Q}}^{0.8} = 0.36 \times 2\;{534^{0.8}} \; {}^\circ {\rm{C}} = 190\;{}^\circ {\rm{C}}$

      (7)

      铁芯的热流密度

      ${q_{\rm{x}}} = \frac{{\Delta {P_{{\rm{Fe}}}}}}{{3{F_{\rm{x}}}}} = \frac{{10\;770}}{{5.9}}\;{\rm{W/}}{{\rm{m}}^{\rm{2}}} = 1\;825\;{\rm{W/}}{{\rm{m}}^{\rm{2}}}$

      (8)

      铁芯对周围空气的温升

      ${\theta _{\rm{x}}} = 0.36q_{\rm{x}}^{0.8} = 0.36 \times 1\;{825^{0.8}}\; {}^\circ {\rm{C}} = 146\;{}^\circ {\rm{C}}$

      (9)

      变压器整体的热流密度

      ${q_{\rm{b}}} = \frac{{{Q_{\rm{b}}}}}{A} = \frac{{31\;779}}{{7.49}}{\rm{W/}}{{\rm{m}}^{\rm{2}}} = 4\;242\;{\rm{W/}}{{\rm{m}}^{\rm{2}}}$

      (10)

      变压器对空气的温升

      ${\theta _{\rm{b}}} = 0.36q_{\rm{b}}^{0.8} = 287\;{}^\circ {\rm{C}}$

      (11)

      PSM干式变压器对空气的整体温升将达到287 ℃,因此仅在自然冷却的工作状态下,变压器的温升效果非常的明显。因此对于该PSM高压电源系统中的干式变压器只采用空气自然散热的方法,将会使变压器的环境温度达到一个影响其性能的值,必须采用其他措施来进行散热,保障变压器的良好性能。通过热力学可知,增加空气的流通速度将会增大换热系数,使热量更快的散发到周围空气中,减小变压器的温升。对此,分析PSM高压电源系统中干式变压器采用强迫风冷进行散热,研究得到其变压器的温升小于变压器规定的温度条件。

      根据热力学可知,对于变压器的总共热量除辐射散热的热量外,其余通过强迫风冷的方法将其热量散发掉,从而避免变压器的过高温升导致变压器的性能受到影响和造成更大的变压器损耗。强迫空气对流的平均温度取出口温度和入口温度的平均值进行计算[8-10]

      强迫对流空气平均温度

      ${T_{\rm{a}}} = \left( {{T_{\rm{ai} }} + {T_{{\rm{ao}}}}} \right)/2$

      (12)

      强迫对流换热能量

      ${Q_{\rm{k}}} = {h_{{\rm{force}}}}A\left( {{T_{{\rm{fo}}}} - {T_{\rm{a}}}} \right)$

      (13)

      强迫空气的温升热能

      ${Q_{\rm{T}}} = \rho AD{c_{\rm{p}}}\left( {{T_{\rm{a}}} - {T_{{\rm{ai}}}}} \right)$

      (14)

      在根据其强迫换流的能量等于空气温升所吸收的能量:${Q_{\rm{k}}}L/v = {Q_{\rm{T}}}$$v$表示空气速度。

      通过利用MATLAB编程计算可以得到热辐射能量为${Q_{\rm{r}}} = 6\;969\;{\rm{W}}$,此时的变压器剩余热量

      ${Q_{\rm{b}}} = {Q_{{\rm{Cu}}}} + {Q_{{\rm{Fe}}}} - {Q_{\rm{r}}} = 39\;832\;{\rm{W}} + 10\;770\;{\rm{W}} - 6\;969\;{\rm{W}} = 43 \; 633\;{\rm{W}}$

      (15)

      对此采用强迫风冷对变压器进行散热分析。采用强迫风冷时,当$v = 2\;{\rm{m/s}}$,可以得到其强迫空气的出口平均温度${T_{{\rm{ao}}}} $=73 ℃。变压器的出口温度与采用的强迫风冷其进口风速的变化如图3所示。可以看出当风速比较低的,由于空气的流通较慢,对此变压器所散发的绝大部分热量将会对风道内的空气进行加热,持续时间也会变长,因此导致变压器的出口温度将会很高,达到200 ℃以上,这对于变压器的稳态运行是不利的。在风速达到$v = 2\;{\rm{m/s}}$以后,其出口温度随着风速的增加下降的趋势变缓。风速的过快,使强迫风在风道的停留时间过短,吸收的热量少,因此转化为空气的内能变小,从而使得强迫风的温度上升偏低,这样可以使变压器维持在一个较低温度的工作环境。由于强迫风需要风机进行驱动,风速的增加会引起风机的功率增大,而风机的功率增大,会降低系统的有功功率,过快的强迫风不能充分利用空气吸收热量,所以干式变压器的风速需要控制在一个合理的范围之内。

      图  3  强迫空气对流出口温度与空气流速的关系

      Figure 3.  Exit temperature changes with wind speed

      PSM高压电源干式变压器在风速为$v = 2\;{\rm{m/s}}$时,利用MATLAB程序可以计算得到强迫对流可以换走热量${Q_{\rm{k}}} = 46\;858\;{\rm{W}}$,其总共散热能量为${Q_{\rm{k}}} + {Q_{\rm{r}}} = 46\;858\;{\rm{W}} + 6\;969\;{\rm{W}} = 53\;827\;{\rm{W}}$。考虑变压器在运行过程中,还存在一部分其他损耗将会变成热量,需要散发。当风速为$2\;{\rm{m/s}}$时,其散热能力足可以对变压器进行散热。

      通过利用ANSYS对变压器的散热进行分析,变压器采用空气冷却,压强为一个标准大气压。风道宽度采用20 mm,通过流体与热的分析,可以得到空气在风冷通道中的温度分布图与风速分布图[11],如图4图5所示。

      图  4  风冷通道中温度场分布等值线图

      Figure 4.  Contour map of temperature field distribution in air-cooled channels

      图  5  风冷通道中Y方向等值线图

      Figure 5.  Y-direction contour map in the air-cooled channel

      分析中将风道的比值变大20倍,以便于观测,同时将变压器进行倒向分析,通过变压器的温度场等值线分布图可以看到变压器壁的温度很高,这是由于空气与变压器壁接触中,空气的流速降低,阻碍变压器的散热性能,并且由于空气与壁的接触面加热密度很高,所以导致该处的温度较高,其最高温度在790 K。通过图5可以看到风速在整个通道中的具体流速,可以发现,温度低的部分其风速越快,温度高的部分风速越慢,整个通道的风速在$0 \sim 3\;{\rm{m/s}}$之间,并不是很大,所以对于驱动风冷的风机要求不会过高。

      当变压器停止工作之后,其变压器中因在工作期间导致的铜导线、矽钢片等的温度需要散发到空气当中,在停止工作之后其变压器的温度变化如图6所示。

      图  6  变压器自然冷却时温度变化图

      Figure 6.  Temperature change diagram of transformer during natural cooling

      从图中可以看出由于铜导体的导热系数很高,所以在铜导体的整个厚度区间之内,他们的温度相差很小,而在变压器停止使用之后,铜导体可以在20 min之内自然冷却到50 ℃左右。

    • 首先对变压器的换热过程进行理论分析,主要分析变压器在对流换热和热辐射换热两种情况下的热损耗,利用MATLAB的拟合工具对变压器的换热空气介质进行拟合,分析其主要参数与温度的关系,并得出相应的函数表达式,将此表达式代入热分析的各参数中进行计算,减少分析误差。通过ANSYS对变压器的对流换热的空气进行温度场和速度分析,得出强迫空气散热情况下,空气的温度分布情况。计算得到在风速为$2\;{\rm{m/s}}$的强迫对流中,其出口温度为73 ℃,并随着风速的增加可进一步降低出口温度,但引起的风机功率消耗将会变大,权衡二者的利弊,取风速值为$2\;{\rm{m/s}}$。在变压器停止运行后,绕组和铁芯在自然冷却条件下,20 min内变压的温度就可冷却到50 ℃左右。本文为PSM多绕组干式整流变压器的散热分析提供了一种研究分析方法,其结果对此类变压器的高性能和安全运行具有实用意义。

参考文献 (11)

目录

    /

    返回文章
    返回