Abstract：The element by element parallel finite element method(EBE-PFEM) has not been applied to engineering eddy current calculation till now. It is not easy to apply this efficient method directly due to the ill-conditioned system equations of engineering eddy current problems which consist of different conductive materials. This paper proposes a solution to the difficulty from both sides of the mathematical model and algorithm implementation. The modified magnetic vector potential and second order vector potential are adopted to simplify the mathematical model and improve the property of the system equations. The EBE-PFEM is introduced to reduce the large memory request and hardness of solution for eddy current analysis in large electrical equipments. Furthermore, the implementation process of EBE-PFEM on graphic processing units(GPU) is given.
唐任远, 吴东阳, 谢德馨. 单元级别并行有限元法求解工程涡流场的关键问题研究[J]. 电工技术学报, 2014, 29(5): 1-8.
Tang Renyuan, Wu Dongyang, Xie Dexin. Research on the Key Problem of Element by Element Parallel FEM Applied to Engineering Eddy Current Analysis. Transactions of China Electrotechnical Society, 2014, 29(5): 1-8.
[1] 历天威, 阮江军. 大规模电磁场数值计算中并行迭代方法的比较[J]. 电工技术学报, 2007, 22(8): 166-172. Li Tianwei, Ruan Jiangjun. Comparison of parallel iterative methods on large-scale electromagnetic field numerical computation[J]. Transactions of China Electrotechnical Society, 2007, 22(8): 166 -172. [2] 李月卉. 电磁场数值求解中迭代方法与预条件技术研究[D]. 成都: 电子科技大学, 2011. [3] Li Yuehui, Nie Zaiping, Meng Min, et al. An efficient MAINV preconditioned COCG method for FEM analysis of millimeter wave filters[J]. J. Infrared Milli Terahz Waves, 2011, 32(2): 216-224. [4] 张健飞. 采用PETSc的有限元并行计算实现与优化[J]. 计算机工程与应用, 2010, 46(10): 57-59. Zhang Jianfei. Parallel implementation and tuning of finite element computing with PETSc[J]. Computer Engineering and Applications, 2010, 46(10): 57-59. [5] Mifune T, Iwashita T, Shimasaki M. A fasst solver for FEM analysis using the parallelized algebraic multi- grid method[J]. IEEE Transactions on Magnetics, 2002, 38(2): 369-372. [6] Steve McFee, Wu Qingying, Mark Dorica, et al. Parallel and distributed processing for h-p adaptive finite-element analysis: a comparison of simulated and empirical studies[J]. IEEE Transactions on Magnetics, 2004, 40(2): 928-933. [7] 历天威, 阮江军, 吴田. 并行计算高压输电线路周围电场[J]. 电工技术学报, 2009, 24(7): 1-6. Li Tianwei, Ruan Jiangjun, Wu Tian. Parallel computation of electric field intensity nearby high voltage transmission lines[J]. Transactions of China Electrotechnical Society, 2009, 24(7): 1-6. [8] 历天威, 阮江军, 杜志叶, 等. 并行计算合成绝缘子串电压分布及金具表面电场强度[J]. 电工技术学报, 2010, 25(3): 6-13. Li Tianwei, Ruan Jiangjun, Du Zhiye, et al. Parallel computation of voltage distribustion along composite insulator strings and electric field intensity on the surface of hardware fittings[J]. Transactions of China Electrotechnical Society, 2010, 25(3): 6-13. [9] 谢开贵, 张怀勋, 胡博, 等. 大规模电力系统潮流计算的分布式GESP算法[J]. 电工技术学报, 2010, 25(6): 89-95. Xie Kaigui, Zhang Huaixun, Hu Bo, et al. Distributed algorithm for power flow of large-scale power system using the GESP technique[J]. Transactions of China Electrotechnical Society, 2010, 25(6): 89-95. [10] Hughus T J R, Levit I, Winget J. An element-by- element solution algorithm for problems of structural and solid mechanics[J]. Computer Methods in Applied Mechanics and Engineering 1983, 36: 241- 254. [11] 王顺绪, 孙柏国, 周树荃. 瞬态热传导问题的并行混EBE算法[J]. 淮海工学院学报, 1999, 8(3): 7-9. Wang Shunxu, Sun Boguo, Zhou Shuquan. A mixed EBE parallel algorithm for Transient heat conduction problems[J]. Journal of Huaihai Institute of Technology, 1999, 8(3): 7-9. [12] 李敏. 基于EBE技术的有限元并行算法的研究与实现[D]. 大连: 大连理工大学, 2007. [13] Liu Yaoru, Zhou Weiyuan, Yang Qiang. A distributed memory parallel element by element scheme based on Jacobi-conditioned conjugate gradient for 3D finite element analysis[J]. Finite Elements in Analysis and Design, 2007, 43: 494-503. [14] Bastos J P A. Is it possible to solve a FEM static case without assembling and storing a Ax = b matrix system?[J]. ICS Newslett, 2009, 16(1): 4-10. [15] Bastos J P A, Sadowski N. A method to solve FEM statics cases without assembling a matrix system: Application to 3Dedge elements[C]. Electric and Magnetic Fields Conference, Mondovi, Italy, 2009: 151-158. [16] Bastos J P A, Nelson Sadowski. A new method to solve 3D magnetodynamic problems without Ax = b system[J]. IEEE Transactions on Mganetics, 2010, 46(3): 3365-3368. [17] 李海良. 基于EBE策略的铝电解槽磁场有限元并行计算研究[D]. 长沙: 中南大学, 2011. [18] 史燕琨, 熊华强. 基于伪并行遗传算法的配电网电容器优化配置[J]. 电力系统保护与控制, 2009, 37(20): 57-60. Shi Yankun, Xiong Huaqiang. Optimal capacitor arrangement based on pseudo parallel genetic algorithm[J]. Power System Protection and Control, 2009, 37(20): 57-60. [19] 黄靖, 张晓锋, 叶志浩. 大型船舶电力系统分布式并行潮流计算方法[J]. 电力系统保护与控制, 2011, 39(7): 50-55. Huang Jing, Zhang Xiaofeng, Ye Zhihao. A distributed parallel load flow calculation method for large ship power system[J]. Power System Protection and Control, 2011, 39(7): 50-55. [20] 简金宝, 杨林峰, 全然, 等. 基于改进多中心校正解耦内点法的动态最优潮流并行算法[J]. 电工技术学报, 2012, 27(6): 232-241. Jian Jinbao, Yang Linfeng, Quan Ran. Parallel algorithm of dynamic optimal power flow based on improved multiple centrality corrections decoupling interior point method[J]. Transactions of China Electrotechnical Society, 2012, 27(6): 232- 241. [21] David M Fernández, Maryam Mehri Dehnavi, Warren J Gross. Alternate parallel processing approach for FEM[J]. IEEE Transactions on Magnetics, 2012, 48(2): 399-402. [22] 张舒, 褚艳利. GPU高性能运算之CUDA[M]. 北京: 中国水利水电出版社, 2009. [23] 夏俊峰, 杨帆, 李静, 等. 基于GPU的电力系统并行潮流计算的实现[J]. 电力系统保护与控制, 2010, 38(18): 100-103. Xia Junfeng, Yang Fan, Li Jing, et al. Implementation of parallel power flow calculation based on GPU[J]. Power System Protection and Control, 2010, 38(18): 100-103. [24] 江涵, 江全元. 基于GPU的交直流电力系统暂态稳定双层并行仿真[J]. 电力系统保护与控制, 2012, 40(21): 102-108. Jiang han, Jiang Quanyuan. A two-level parallel transient stability algorithm for AC/DC power system based on GPU platform[J]. Power System Protection and Control, 2012, 40(21): 102-108. [25] 柳有权, 尹康学, 吴恩华. 大规模稀疏线性方程组的GMRES-GPU快速求解算法[J]. 计算机辅助设计与图形学报, 2011, 23(4): 553-560. Liu Youquan, Yin Kangxue, Wu Enhua. Fast GMRES-GPU solver for large scale sparse linear systems[J]. Jounal of Computer-Aided Design & Computer Graphics, 2011, 23(4): 553-560. [26] 刘小虎. 大规模有限元系统的GPU加速计算研究 [J]. 计算力学学报, 2012, 2(1): 146-152. Liu Xiaohu. Solving large finite element system by GPU computation[J]. Chinese Journal of Computational Mechanics, 2012, 2(1): 146-152. [27] Barnat J, Bauch P. Employing multiple CUDA devices to accelerate LTL model checking[C]. IEEE 16th International Conference on Parallel and Distributed System, 2010: 259-266. [28] Thurley M J, Danell V. Fast morphological image processing open-source extensions for GPU processing with CUDA[J]. IEEE Journal of Selected Topics in Signal Processing, 2012, 6(7): 849-855. [29] 谢德馨. 三维涡流场的有限元分析[M]. 北京: 机械工业出版社, 2007. [30] Kershaw D S. The incomplete Cholesky-conjugate gradient method for iterative solution of system of linear equations[J]. Journal of Computational Physics, 1978, 26: 43-65. [31] Biro O, Preis K. On the use of the magnetic vector potential in the finite-element analysis of three- dimensional eddy currents[J]. IEEE Transactions on Magnetics, 1989, 25(4): 531-535. [32] 周克定. 工程电磁场数值计算的理论方法与应用[M]. 北京: 高等教育出版社, 1994. [33] 雷银照, 徐纪安. 时变电磁场唯一性定理的完整表述[J]. 电工技术学报, 2000, 15(1): 16-20. Lei Yinzhao, Xu Jian. New statement of uniqueness theorem for time-varying electromagnetic field[J]. Transactions of China Electrotechnical Society, 2000, 15(1): 16-20. [34] Lei Yinzhao, Qian Baoliang, Liao Xiaohui. On the divergence expression and the boundary value of modified magnetic vector potential[J]. Journal of Zhengzhou University of Technology, 2000, 21(3): 5-10. [35] 雷银照, 熊华俊, 王书彬. 线性瞬态涡流场定解问题中的法向边界条件与解的唯一性[J]. 中国电机工程学报, 2003, 23(4): 81-85. Lei Yinzhao, Xiong Huajun, Wang Shubin. The normal interface condition and the uniqueness of solution to the linear transient eddy current definite- solution problem[J]. Proceedings of the CSEE, 2003, 23(4): 81-85. [36] 雷银照. 时谐电磁场解析方法[M]. 北京: 科学出版社, 2000. [37] 毛雪飞, 雷银照. 金属管道外放置式偏心弯曲矩形线圈涡流场的解析解[J]. 电工技术学报, 2012, 27 (9): 153-159. Mao Xuefei, Lei Yinzhao. Analytical solution for the eddy current field due to air-cored eccentric curved rectangular probe coil on the outside of a metal pipe[J]. Transactions of China Electrotechnical Society, 2012, 27(9): 153-159. [38] Smythe W R. Static and dynamic electricity[M]. McGraW-Hill, 1968. [39] 周树荃, 梁维泰. 有限元结构分析并行计算[M]. 北京: 科学出版社, 1994. [40] 蔡大用, 白峰杉. 高等数值分析[M]. 北京: 清华大学出版社, 1997. [41] G H戈卢布, C F范洛恩. 矩阵计算[M]. 袁亚湘译.北京: 科学出版社, 2001. [42] 李路涛. GPU加速的大规模岩土工程有限元计算中的迭代求解[D]. 北京: 北京交通大学, 2012.
2014, Vol. 29 
