Abstract：The development of digital twin technology for electrical equipment requires extensive calculation. For the electrical equipment with large physical dimension or complicated structure, a large number of meshing grids are needed to analyze its physical field accurately using finite element method(FEM), which leads to a huge order of algebraic equations and makes the solution very difficult, and sometimes it cannot be solved. For sparse symmetric equations obtained by finite element analysis, this paper researches effective order-reduced algorithms. Firstly, the block iterative algorithm is derived and explored in detail, in which the block diagonalization method is proposed to transform the finite element stiffness matrix into the diagonal matrix. The diagonal matrix is the matrix with unit matrix at diagonal, and the matrix is divided into sub-matrix blocks according to the unit matrix. Next, the block left multiplication method is proposed to realize the movement of sub-matrix blocks and block reorganization, and the permutation matrix method is also derived to realize the process. The block reorganization can reduce the order of the equations to half of its original order, and can reduce the order repeatedly. Taking the finite element analysis of the magnetic field of the open slot of the motor as an example, the block iterative algorithm is used twice to reduce the order of the equations, and the original 10th-order stiffness matrix is reduced to 3rd order. Next, the mathematical model of orthogonal order-reduced decomposition algorithm(OORDA) is also derived. OORDA is to orthogonalize any column of the matrix with other column vectors, removes elements from corresponding row and column, and the matrix will be reduced by one order. This operation is repeated until the matrix is reduced to first-order, and then substitute the result back to complete the solution. Aiming at the characteristics of sparsity of finite element stiffness matrix, the taboo search algorithm is introduced into the Gaussian elimination method to form an improved Gaussian elimination method(IGEM).The information of non-zero elements is extracted and recorded in the taboo information table, so as to eliminate the memory occupation and computation increase caused by zero elements during elimination and back-substitution processes, which will greatly reduce the number of operations and computational cost of Gaussian Elimination method. Due to the notable memory requirements of both OORDA and Gaussian elimination method, and the fact that the block iterative algorithm can substantially reduce the order of the equations, the Hybrid Algorithm 1 (HA1) that combine the block iterative algorithm with OORDA and Hybrid Algorithm 2 (HA2) that combine the block iterative algorithm with IGEM are proposed in this paper. The four algorithms, OORDA, IGEM, HA1 and HA2 are applied to the finite element analysis of the magnetic field produced by a long straight wire carrying current, the correctness of the algorithms have been verified by comparing the calculation results with the analytical solution. The four algorithms are applied to calculate the magnetic field distribution of a single-phase transformer under different meshes. The results show that using the block iterative algorithm to reduce the order of the equations, and then employing the OORDA or IGEM to solve them can greatly improve the computational efficiency. The larger the calculation scale, the more significant the effect. The HA2 shows the highest computational efficiency. The OORDA is not suitable for solving high-order equations alone, but it is a feasible choice for matrices that are inconvenient to inverse. All the algorithms in this paper are realized by programming in C++, and the order reduction algorithms can provide a solution approach and reference for large-scale numerical calculation.
金军, 阎秀恪, 钟立国, 张艳丽, 任自艳. 方程组降阶算法及其在电工装备数值计算中的应用[J]. 电工技术学报, 2025, 40(7): 2020-2032.
Jin Jun, Yan Xiuke, Zhong Liguo, Zhang Yanli, Ren Ziyan. Equation System Order-Reduced Algorithms and Their Application in Numerical Computation for Electrical Equipment. Transactions of China Electrotechnical Society, 2025, 40(7): 2020-2032.
[1] 董华军, 温超阳, 孙鹏, 等. 基于正交实验新型真空灭弧室触头磁场仿真与参数优化设计[J]. 电工技术学报, 2022, 37(21): 5598-5606. Dong Huajun, Wen Chaoyang, Sun Peng, et al.Simulation and optimization of the contact magnetic field of a new type of vacuum interrupter based on orthogonal experiment[J]. Transactions of China Electrotechnical Society, 2022, 37(21): 5598-5606. [2] 刘云飞, 张炳义, 宗鸣, 等. 基于非线性混合模型的模块组合式永磁电机磁场解析[J]. 电工技术学报, 2022, 37(18): 4593-4603. Liu Yunfei, Zhang Bingyi, Zong Ming, et al.Analytical prediction of magnetic field in modular combined permanent magnet motor by a nonlinear hybrid model[J]. Transactions of China Electrotechnical Society, 2022, 37(18): 4593-4603. [3] 邱浩, 王曙鸿, 孙凤举, 等. 基于时域有限积分技术的四级串联快脉冲直线型变压器驱动源电磁特性[J]. 电工技术学报, 2022, 37(4): 816-826. Qiu Hao, Wang Shuhong, Sun Fengju, et al.The electromagnetic characteristics of the four-stage series-connected fast linear transformer driver based on time-domain finite integration technique[J]. Transactions of China Electrotechnical Society, 2022, 37(4): 816-826. [4] 黄天超, 王泽忠, 李宇妍. 换流变压器直流偏磁对油箱涡流损耗的影响[J]. 电工技术学报, 2023, 38(8): 2004-2014. Huang Tianchao, Wang Zezhong, Li Yuyan.The influence of converter transformer DC bias on eddy current loss of tank[J]. Transactions of China Electrotechnical Society, 2023, 38(8): 2004-2014. [5] Guo Ze, Wang Shuhong, Tang Zuqi, et al.Fast time-domain solution of dynamic electromagnetic problems based on sinc interpolation[J]. IEEE Transactions on Magnetics, 2021, 57(6): 7200704. [6] 赵志刚, 王丽美, 陈天缘, 等. 基于等效复数磁导率的利兹线绕组损耗计算模型[J]. 电工技术学报, 2024, 39(4): 948-955. Zhao Zhigang, Wang Limei, Chen Tianyuan, et al.Calculation model of winding loss of Litz-wire based on equivalent complex permeability[J]. Transactions of China Electrotechnical Society, 2024, 39(4): 948-955. [7] 龚天勇, 马光同, 任刚. 跑道型高温超导磁体的磁场解析计算模型[J]. 电工技术学报, 2023, 38(8): 1991-2003. Gong Tianyong, Ma Guangtong, Ren Gang.Analytical models for the magnetic field calculation of racetrack high-temperature superconducting magnet[J]. Transac-tions of China Electrotechnical Society, 2023, 38(8): 1991-2003. [8] 刘刚, 胡万君, 郝世缘, 等. 油浸式变压器绕组瞬态温升降阶快速计算方法[J]. 电工技术学报, 2024, 39(3): 643-657. Liu Gang, Hu Wanjun, Hao Shiyuan, et al.Reduced order calculation method of steady temperature rise of oil immersed power transformer[J]. Transactions of China Electrotechnical Society, 2024, 39(3): 643-657. [9] 刘刚, 郝世缘, 朱章宸, 等. 基于DMD-ATS油浸式电力变压器绕组瞬态温升快速计算方法研究[J]. 电工技术学报, 1-12[2025-04-02]. https://doi.org/10.19595/j.cnki.1000-6753.tces.230465. Liu Gang, Hao Shiyuan, Zhu Zhangchen, et al. Research on rapid calculation method of transient temperature rise of winding of DMD-ATS oil-immersed power transformer[J]. Transactions of China Electro-technical Society, 1-12[2025-04-02].https://doi.org/10.19595/j.cnki.1000-6753.tces.230465 [10] 康军, 马凡, 胡健, 等. 一种多时间尺度线性系统模型降阶的误差预测方法[J]. 电工技术学报, 2017, 32(3): 56-64. Kang Jun, Ma Fan, Hu Jian, et al.An error prediction method of model order reduction for multi-time scale linear system[J]. Transactions of China Electro-technical Society, 2017, 32(3): 56-64. [11] 许建成, 孙建军, 钟佩军, 等. 基于平衡实现理论的变流器并网系统降阶模型[J]. 电工技术学报, 2021, 36(增刊1): 255-264. Xu Jiancheng, Sun Jianjun, Zhong Peijun, et al.Reduced-order model of grid-connected converter system based on balanced realization theory[J]. Transactions of China Electrotechnical Society, 2021, 36(S1): 255-264. [12] 陈鹏伟, 戚陈陈, 陈新, 等. 附加频率控制双馈风电场频率响应特性建模与参数辨识[J]. 电工技术学报, 2021, 36(15): 3293-3307. Chen Pengwei, Qi Chenchen, Chen Xin, et al.Frequency response modeling and parameter identification of doubly-fed wind farm with additional frequency control[J]. Transactions of China Electrotechnical Society, 2021, 36(15): 3293-3307. [13] 顾理成, 陈前, 赵文祥, 等. 五相永磁容错电机的相间短路容错控制[J]. 电工技术学报, 2022, 37(8): 1972-1981. Gu Licheng, Chen Qian, Zhao Wenxiang, et al.Inter-phase short-circuit fault-tolerant control for five-phase permanent magnet fault-tolerant motors[J]. Transactions of China Electrotechnical Society, 2022, 37(8): 1972-1981. [14] 马益平, 姚艳, 汪雅静, 等. 配电网电压稳定的分布式评估方法[J]. 电测与仪表, 2021, 58(3): 28-33. Ma Yiping, Yao Yan, Wang Yajing, et al.Distributed evolution method of voltage stability in radial distribution networks[J]. Electrical Measurement & Instrumentation, 2021, 58(3): 28-33. [15] Liu Zhaocheng, Cui Xiang, Li Xuebao, et al.Interface finite element method with lower calculation amount for step size-varying electro-quasistatic field prob-lem[J]. CSEE Journal of Power and Energy Systems, 2023, 9(5): 1710-1719. [16] Qiu Hao, Wang Shuhong, Zhang Naming, et al.Field-circuit coupling and electromagnetic-thermal-mechanical coupling analysis of the single-stage fast linear transformer driver using time-domain finite integration technique[J]. IEEE Transactions on Magnetics, 2021, 57(6): 8401205. [17] 李胜利, 王川龙, 温瑞萍. 求解非奇异线性方程组的正交降阶法[J]. 中北大学学报(自然科学版), 2014, 35(3): 248-251, 257. Li Shengli, Wang Chuanlong, Wen Ruiping.Ortho-gonal reduced-order method for solving nonsingular linear systems[J]. Journal of North University of China (Natural Science Edition), 2014, 35(3): 248-251, 257. [18] 王芙蓉, 杨帆, 张亚, 等. 基于奇异值分解的矩阵低秩近似量子算法[J]. 物理学报, 2021, 70(15): 7-14. Wang Furong, Yang Fan, Zhang Ya, et al.Matrix low-rank approximate quantum algorithm based on singular value decomposition[J]. Acta Physica Sinica, 2021, 70(15): 7-14. [19] 李文强, 刘晓. 循环三对角Toeplitz线性方程组的分组降阶算法[J]. 科技导报, 2012, 30(5): 43-48. Li Wenqiang, Liu Xiao.Grouping and order reducing algorithm for solving Toeplitz type circular tri-diagonal liner algebraic equation systems[J]. Science & Technology Review, 2012, 30(5): 43-48. [20] Sabariego R V, Vanbroekhoven B, Gyselinck J, et al.RL-ladder circuit models for eddy-current problems with translational movement[J]. IEEE Transactions on Magnetics, 2022, 58(9): 6301304. [21] 尹松柱, 谭领. 解某些线代数方程组的降阶方法[J]. 东北师大学报(自然科学版), 1965(2): 17-24. Ing Shungzhu, Tan Ling.A lowering order method for solving certain system of linear equations[J]. Journal of Northeast Normal University, 1965(2): 17-24. [22] 刘长河. 大型稀疏线性方程组的数值解法[J]. 北京建筑大学学报, 2023, 39(1): 103-108. Liu Changhe.Numerical methods for solving the large-scale sparse linear equations[J]. Journal of Beijing University of Civil Engineering and Archi-tecture, 2023, 39(1): 103-108. [23] 胡之光. 电机电磁场的分析与计算[M]. 北京: 机械工业出版社, 1982.
2025, Vol. 40 
