基于改进传输线法的电器温度场并行有限元计算

doi:10.16628/j.cnki.2095-8188.2022.03.005

摘要/Abstract

摘要：

大功率电器产品的优化设计,以及电器多物理场耦合仿真的发展,对电器温度场仿真的求解精度以及求解效率提出了更高的需求。使用传输线法(TLM)加速三维非线性热传导有限元(FEM)的求解,为进一步提高计算效率,提出了传输线导纳确定原则和松弛迭代方法,并且使用并行计算技术来实现TLM的加速求解。将改进的TLM应用于直动式继电器温度场的求解中,通过与商用软件COMSOL以及Newton-Raphson (NR)迭代法的求解结果对比,验证了算法的精度。相比于传统的NR迭代法,TLM能够取得稳定的求解效率提升。

关键词: 有限元分析, 传输线法, 并行计算, 继电器, 温度场

Abstract:

Higher solution accuracy and solution efficiency of the temperature field simulation are required,to satisfy the optimization design and the development of multi-physical field coupling simulation of high-power electrical apparatus.The transmission line iteration method (TLM) is used to accelerate the finite element analysis of three-dimensional nonlinear heat conduction problem.In order to further improve the efficiency of TLM,the transmission line admittance determination principle and relaxation iteration method are proposed.The parallel computing technology is used to realize the acceleration of TLM.The improved TLM is applied to solve the temperature field of the direct-acting relay,and the accuracy of the algorithm is verified by comparing the solution results with commercial software COMSOL and Newton-Raphson iteration.Compared with the traditional Newton-Raphson iteration,TLM can achieve a stable solution efficiency improvement.

Key words: finite element analysis, transmission line method, parallel computing, relay, temperature field

中图分类号:

TM764

邱子澜, 杨文英, 彭飞, 翟国富. 基于改进传输线法的电器温度场并行有限元计算[J]. 电器与能效管理技术, 2022, 0(3): 30-38.

QIU Zilan, YANG Wenying, PENG Fei, ZHAI Guofu. Improved Transmission Line Method for Parallel Finite Element Calculation in Temperature Field of Electrical Apparatus[J]. LOW VOLTAGE APPARATUS, 2022, 0(3): 30-38.

图/表 12

图1

图2

图3

图4

图5

图6

表1

图7

表2

表3

表4

图8

参考文献 29

[1]	李丹, 张代润, 杨林. 电磁继电器温升特性的仿真研究[J]. 电气技术, 2018, 19(3):26-30.
[2]	姜宁, 任秀红, 李合, 等.50 kA大电流温升试验系统的建立与研究[J]. 电器与能效管理技术, 2018(13):59-63.
[3]	翟国富, 崔行磊, 杨文英. 电磁继电器产品及研究技术发展综述[J]. 电器与能效管理技术, 2016(2):1-8.
[4]	薄凯, 周学, 翟国富, 等. 大电流条件下真空接触器温升特性的试验测量与仿真分析[J]. 电工技术学报, 2019, 34(24):81-89.
[5]	李茹瑶, 杨文英, 肖斌, 等. 高压直流继电器线圈温度场数值算法[J]. 电工技术学报, 2018, 33(增刊2):194-202.
[6]	杨文英, 郭久威, 王茹, 等. 继电器电磁机构电磁-热耦合模型建立与计算方法[J]. 电工技术学报, 2017, 32(13):169-177.
[7]	屈建宇, 徐宏宇. 万能式断路器温升特性仿真与实验研究[J]. 电器与能效管理技术, 2021(4):40-45.
[8]	周金发, 李茹瑶, 周文博, 等. 功率型电磁继电器热特性计算方法研究[J]. 电器与能效管理技术, 2017(4):32-34,70.
[9]	杨文强, 张蓬鹤, 张保亮. 小型断路器温升预测与仿真分析[J]. 电器与能效管理技术, 2020(6):25-29.
[10]	黄世泽, 郭其一, 李凡璋, 等. 基于Ansys的电磁系统温度场仿真技术研究[J]. 电器与能效管理技术, 2015(2):32-37.
[11]	蔡文智. 基于有限元法的继电器三维稳态温度场分析[J]. 电器与能效管理技术, 2015(10):32-36.
[12]	纽春萍, 陈德桂, 刘颖异, 等. 计及主回路和电磁系统发热的交流接触器数值热分析[J]. 中国电机工程学报, 2007, 27(15):53-58.
[13]	苏秀苹, 陆俭国, 刘帼巾, 等. 小型直流电磁继电器温度场仿真分析[J]. 电工技术学报, 2011, 26(8):185-189.
[14]	DLALA E, ARKKIO A. General formulation for the Newton-Raphson method and the fixed-point method in finite-element programs[C]// XIX International Conference on Electrical Machines-ICEM 2010.IEEE, 2010:1-5.
[15]	DOMMEL H W. Digital computer solution of electromagnetic transients in single-and multiphase networks[J]. IEEE Transactions on Power Apparatus and Systems, 1969, PAS-88(4):388-399. doi: 10.1109/TPAS.1969.292459
[16]	LOBRY J, TRECAT J, BROCHE C. Symmetry and TLM method in nonlinear magnetostatics[J]. IEEE Transactions on Magnetics, 1996, 32(3):702-705. doi: 10.1109/20.497336
[17]	FLACK T J, KNIGHT R E. On the domain decomposition and transmission line modelling finite element method for time-domain induction motor analysis[J]. IEEE Transactions on Magnetics, 2002, 35(3):1290-1293. doi: 10.1109/20.767187
[18]	KNIGHT A M. Time-stepped eddy-current analysis of induction machines with transmission line modeling and domain decomposition[J]. IEEE Transactions on Magnetics, 2003, 39(4):2030-2035. doi: 10.1109/TMAG.2003.812707
[19]	YANG W, PENG F, DINAVAHI V. Nonlinear axisymmetric magnetostatic analysis for electromagnetic device using TLM-based finite-element method[J]. IEEE Transactions on Magnetics, 2017, 53(4):1-9.
[20]	彭飞, 杨文英, 翟国富. 基于优化的传输线法的并行静磁场有限元方法[J]. 电工技术学报, 2019, 34(13):2716-2725.
[21]	彭飞, 杨文英, 杨春玲. 基于传输线法的非线性静磁场四边形有限元加速求解方法研究[J]. 中国电机工程学报, 2019, 39(13):3947-3955.
[22]	YANG W, PENG F, DINAVAHI V,et al. A generalized parallel transmission line iteration for finite element analysis of permanent magnet axisymmetrical actuator[J]. IEEE Transactions on Magnetics, 2019, 55(3):1-10.
[23]	LIU P, DINAVAHI V. Real-time finite-element simulation of electromagnetic transients of transformer on FPGA[J]. IEEE Transactions on Power Delivery, 2018, 33(4):1991-2001. doi: 10.1109/TPWRD.2018.2812753
[24]	LIU P, LI J, DINAVAHI V. Matrix-free nonlinear finite-element solver using transmission-line modeling on GPU[J]. IEEE Transactions on Magnetics, 2019(7):1-5.
[25]	LOBRY J, TRECAT J, BROCHE C. The transmission line modeling (TLM) method as a new iterative technique in nonlinear 2-D magnetostatics[J]. IEEE transactions on magnetics, 1996, 32(2):559-566. doi: 10.1109/20.486548
[26]	杨世铭, 陶文铨. 传热学[M]. 4版. 北京: 高等教育出版社, 2006.
[27]	RAO S S. The finite element method in engineering[M]. Fourth Edition.xFord:Butterworth-Heinemann, 2005.
[28]	XIAOYE S.Li. An Overview of SuperLU:Algorithms,Implementation,and user interface[J]. Acm Transactions on Mathematical Software, 2005, 31(3):302-325. doi: 10.1145/1089014.1089017
[29]	胡方, 杨文英, 赵瑞平, 等. 大功率直流接触器温度场仿真及影响因素分析[J]. 低压电器, 2011(23):1-6.

材料名称	组件	热导率(数据来源于COMSOL材料库)
黄铜	静触头	λ=555.4-2.116 905T+0.008 971 429 T²-1.266 667×10^-5T³
	动触头
	线圈
	轴套
DT4E	轭铁	λ=155.635 5-0.447 094 3T+ 8.992 2×10^-4T²-9.859 391×10^-7 T³+4.862 049×10^-10T⁴-7.562 752× 10^-14T⁵
	孔轭铁
	衔铁
	外壳
不锈钢	连杆	λ=6.742 253+0.028 649 15T
陶瓷	陶瓷座	λ=16.7
聚酰亚胺	线圈骨架	λ=-0.007 707 532+0.005 769 13T+ 5.622 796×10^-7T²-4.329 984× 10^-10T³
空气	其余部分	λ=-0.002 275 8+1.154 80×10^-4T -7.902 52×10^-8T²+4.117 02× 10^-11T³-7.438 64×10^-15T⁴

节点编号	坐标/mm	COMSOL 温度/K	NR温度/K	TLM温度/K
1	(0,0,0)	323.026	323.091	323.091
2	(-1.27,5.86,-6.86)	323.068	323.129	323.129
3	(-8.50,1.06,14.90)	321.430	321.521	321.542
4	(-4.55,6.65,20.20)	320.113	320.272	320.273
5	(-13.45,18.78,4.39)	323.183	323.215	323.315
6	(2.56,-19.01,15.77)	321.390	321.490	321.490
7	(17.37,-15.14,13.39)	321.909	322.029	322.029
8	(10.01,13.89,23.30)	320.458	320.594	320.595
9	(20.55,-2.48,26.30)	320.891	321.013	321.013
10	(-23.92,-3.30,-32.50)	322.971	323.030	323.030

节点编号	NR相对误差/%	TLM相对误差/%
1	0.020	0.020
2	0.018	0.018
3	0.028	0.035
4	0.049	0.050
5	0.010	0.041
6	0.031	0.031
7	0.037	0.037
8	0.042	0.042
9	0.038	0.038
10	0.018	0.018

编号	单元数目	NR		TLM
		迭代次数/次	单步迭代时间/s	迭代次数/次	单步迭代时间/s
		1	540 580	4	85.82	235	0.181
2	492 177	4	63.31	250	0.161
3	344 274	4	23.62	285	0.123
4	174 825	4	6.48	269	0.052