-
The pulsed high current linear driving device operates under extreme conditions, and various forms of metal damages will reduce the service life of the device. At present, the multi-physics coupling mechanism of pulsed high current linear driving device is still unclear, and the multi-parameter diagnosis method in the laboratory environment is limited. Therefore, it is urgent to clarify the evolution process of multiple physical parameters through numerical modeling methods, in order to guide the optimization of the overall performance and improve the service life of the device. In this work, mathematical and physical models of electromagnetic field, temperature field and structural field under dynamic conditions are established. The local solution is carried out by using the characteristics of rail reverse motion and the invariant physical quantities at the distal end of the contact. The constraint equations of the non-equipotential surface of the rail entrance and the armature-rail interface conditions under the technical framework are derived. The constraint equations are applied by the penalty function method. The model also takes into account the practical factors such as the temperature dependence of the material properties, thermal stresses, and the frictional heat of the contact surface. The finite element discrete format of the electromagnetic field and the temperature field is solved in the form of Euler’s backward differentiation, and the structural field is solved by the Newmark method. The reliability of the model is verified by comparing the calculation results with the numerical tools EMAP3D and Comsol under the same configuration and input conditions, as well as related experiments. Through the numerical simulation of the C-type armature, the typical evolution process of the corresponding multi-parameter is obtained. During sliding electrical contact, the velocity skin effect becomes more pronounced with velocity increasing. The current is gradually concentrated on the surface of the rail, and the highest current density is found at the rear edge of the contact surface and at the edge of the outer arm of the armature. Moreover, the magnetic induction intensity at the tail of the contact surface continues to shrink over time. The heat-concentrated region appears at the top edge of the contact surface, and with time going by, it extends along the sliding direction and bottom direction of the armature. In addition, there is peak stress in the front of the rail contact and significant stress at the armature throat. When the local stress at the throat of the armature exceeds the corresponding yield strength, it can cause serious deformation or even fracture of the armature.
-
Keywords:
- finite element numerical simulation /
- sliding electrical contact /
- friction heat /
- local modeling
[1] Fair H D 2001 IEEE Trans. Magn. 37 25
Google Scholar
[2] 马伟明, 鲁军勇 2023 电工技术学报 38 3943
Google Scholar
Ma W M, Lu J Y 2023 Trans. Chin. Electrotech. Soc. 38 3943
Google Scholar
[3] Sun J, Cheng J, Wang Q 2021 IEEE Trans. Plasma Sci. 49 3988
Google Scholar
[4] Stefani F, Parker J V 1999 IEEE Trans. Magn. 35 312
Google Scholar
[5] Li S, Li J, Xia S, Zhang Q, Liu P 2019 IEEE Trans. Plasma Sci. 47 2399
Google Scholar
[6] Sun J, Cheng J, Wang Q 2022 IEEE Trans. Plasma Sci. 50 1032
Google Scholar
[7] 殷强, 张合, 李豪杰, 史云雷 2016 强激光与粒子束 28 025008
Google Scholar
Yin Q, Zhang H, Li H J, Shi Y L 2016 High Power Laser Part. Beams 28 025008
Google Scholar
[8] 李昕, 翁春生 2009 火炮发射与控制学报 2009 1
Google Scholar
Li X, Weng C S 2009 J. Gun Launch Control 2009 1
Google Scholar
[9] Li X, Weng C S 2008 Prog. Nat. Sci. 18 1565
Google Scholar
[10] Tang B, Xu Y, Lin Q, Li B 2017 IEEE Trans. Plasma Sci. 45 1361
Google Scholar
[11] 郑杜成, 徐蓉, 成文凭, 赵伟康, 袁伟群, 严萍 2019 电工电能新技术 38 33
Google Scholar
Zheng D C, Xu R, Cheng W P, Zhao W K, Yuan W Q, Yan P 2019 Adv. Technol. Electr. Eng. Energy 38 33
Google Scholar
[12] 翟小飞, 杨帆, 张晓, 刘华 2021 海军工程大学学报 33 19
Google Scholar
Zhai X F, Yang F, Zhang X, Liu H 2021 J. Naval Univ. Eng. 33 19
Google Scholar
[13] Hsieh K T 1995 IEEE Trans. Magn. 31 604
Google Scholar
[14] Hsieh K T, Kim B K 1999 IEEE Trans. Magn. 35 166
Google Scholar
[15] Hsieh K T 2007 IEEE Trans. Magn. 43 1131
Google Scholar
[16] Lin Q H, Li B M 2020 IEEE Trans. Plasma Sci. 48 2287
Google Scholar
[17] Lin Q H, Li B M 2016 Defence Technol. 12 101
Google Scholar
[18] 林庆华, 栗保明 2020 兵工学报 41 1697
Google Scholar
Lin Q H, Li B M 2020 Acta Armamentarii 41 1697
Google Scholar
[19] Shatoff H, Pearson D A, Kull A E 2005 IEEE Pulsed Power Conference Monterey, CA, USA, June 13–17, 2005 p253
[20] Wang G H, Xie L, He Y, Song S Y, Gao J J 2016 IEEE Trans. Plasma Sci. 44 1424
Google Scholar
[21] 王刚华, 谢龙, 赵海龙 2021 爆炸与冲击 41 111
Google Scholar
Wang G H, Xie L, Zhao H L 2021 Explosion and Shock Waves 41 111
Google Scholar
-
图 1 多物理场耦合模型
Figure 1. Multi physical field coupling model.
图 2 基于轨道反向运动的滑动电接触时间演化示意图
Figure 2. Schematic of time evolution of sliding electrical contact based on rail reverse motion.
图 3 块状电枢下电流波形的比较
Figure 3. Comparison of current waveforms under block armature.
图 4 网格及结构场边界条件
Figure 4. Mesh and structural field boundary conditions.
图 5 X, Y, Z三方向位移三维对比 (a) X方向; (b) Y方向; (c) Z方向
Figure 5. 3D comparison of displacement in X, Y and Z directions: (a) X-direction; (b) Y-direction; (c) Z-direction.
图 6 磁场对比
Figure 6. Magnetic field comparison.
图 7 电流对比
Figure 7. Current comparison.
图 8 网格剖分示意图(1/4)
Figure 8. Schematic diagram of mesh (1/4).
图 9 激励电流波形
Figure 9. Excitation current waveform.
图 10 速度及位移波形
Figure 10. Velocity and displacement waveforms.
图 11 电流对比图
Figure 11. Current comparison diagram.
图 12 电流密度的演化 (a) 0.1 ms; (b) 0.2 ms; (c) 0.6 ms
Figure 12. Evolution of current density: (a) 0.1 ms; (b) 0.2 ms; (c) 0.6 ms.
图 13 磁感应强度的演化 (a) 0.1 ms; (b) 0.2 ms; (c) 0.6 ms
Figure 13. Evolution of the magnetic induction intensity: (a) 0.1 ms; (b) 0.2 ms; (c) 0.6 ms.
图 14 经验公式及动态计算的电磁推力比较
Figure 14. Comparison of electromagnetic thrust based on empirical formula and dynamic calculation.
图 15 0.18 ms时的导体温度分布
Figure 15. Temperature distribution of conductor at 0.18 ms.
图 16 电枢接触面温度演化 (a) 0.1 ms; (b) 0.3 ms
Figure 16. Temperature evolution of armature contact surface: (a) 0.1 ms; (b) 0.3 ms.
图 17 0.6 ms时结构场位移分布 (a) X方向; (b) Y方向; (c) Z方向
Figure 17. Displacement distribution of structural field at 0.6 ms: (a) X-direction; (b) Y-direction; (c) Z-direction.
图 18 0.5 ms时的von Mises应力分布
Figure 18. Von Mises stress distribution at 0.5 ms.
图 19 电磁法向力与接触压力的比较
Figure 19. Comparison of electromagnetic normal force and contact force.
图 20 t = 0.6 ms时电枢接触面温度 (a) 无接触压力; (b) 电磁法向力近似接触压力; (c) 结构场接触压力
Figure 20. Temperature of the armature contact surface at t = 0.6 ms: (a) No contact pressure; (b) electromagnetic normal force approximates contact pressure; (c) structural field contact pressure.
表 1 仿真模型电磁场参数
Table 1. Electromagnetic field parameters of simulation model.
参数 初始电导率
σ0/(MS·m–1)磁导率
μ/(10–6 H·m–1)质量密度
ρ/(kg·m–3)比热容
c/(J·kg–1·K–1)热导率
λ/(W·m–1·K–1)电阻率温度系数
α/K熔化温度
Tm/K铜 58.8 1.2566 8960 385 385 0.0039 1356.15 铝 25 1.2566 2700 896 167 0.0041 924.85 
DownLoad: CSV
表 2 仿真模型热场及结构场参数
Table 2. Thermal field and structural field parameters of simulation model.
参数 弹性模量
E/GPa线膨胀系数
β/(µm·m–1·K–1)泊松比 μPoisson 铜 110 16.5 0.34 铝 68.9 23.6 0.33
DownLoad: CSV
-
[1] Fair H D 2001 IEEE Trans. Magn. 37 25
Google Scholar
[2] 马伟明, 鲁军勇 2023 电工技术学报 38 3943
Google Scholar
Ma W M, Lu J Y 2023 Trans. Chin. Electrotech. Soc. 38 3943
Google Scholar
[3] Sun J, Cheng J, Wang Q 2021 IEEE Trans. Plasma Sci. 49 3988
Google Scholar
[4] Stefani F, Parker J V 1999 IEEE Trans. Magn. 35 312
Google Scholar
[5] Li S, Li J, Xia S, Zhang Q, Liu P 2019 IEEE Trans. Plasma Sci. 47 2399
Google Scholar
[6] Sun J, Cheng J, Wang Q 2022 IEEE Trans. Plasma Sci. 50 1032
Google Scholar
[7] 殷强, 张合, 李豪杰, 史云雷 2016 强激光与粒子束 28 025008
Google Scholar
Yin Q, Zhang H, Li H J, Shi Y L 2016 High Power Laser Part. Beams 28 025008
Google Scholar
[8] 李昕, 翁春生 2009 火炮发射与控制学报 2009 1
Google Scholar
Li X, Weng C S 2009 J. Gun Launch Control 2009 1
Google Scholar
[9] Li X, Weng C S 2008 Prog. Nat. Sci. 18 1565
Google Scholar
[10] Tang B, Xu Y, Lin Q, Li B 2017 IEEE Trans. Plasma Sci. 45 1361
Google Scholar
[11] 郑杜成, 徐蓉, 成文凭, 赵伟康, 袁伟群, 严萍 2019 电工电能新技术 38 33
Google Scholar
Zheng D C, Xu R, Cheng W P, Zhao W K, Yuan W Q, Yan P 2019 Adv. Technol. Electr. Eng. Energy 38 33
Google Scholar
[12] 翟小飞, 杨帆, 张晓, 刘华 2021 海军工程大学学报 33 19
Google Scholar
Zhai X F, Yang F, Zhang X, Liu H 2021 J. Naval Univ. Eng. 33 19
Google Scholar
[13] Hsieh K T 1995 IEEE Trans. Magn. 31 604
Google Scholar
[14] Hsieh K T, Kim B K 1999 IEEE Trans. Magn. 35 166
Google Scholar
[15] Hsieh K T 2007 IEEE Trans. Magn. 43 1131
Google Scholar
[16] Lin Q H, Li B M 2020 IEEE Trans. Plasma Sci. 48 2287
Google Scholar
[17] Lin Q H, Li B M 2016 Defence Technol. 12 101
Google Scholar
[18] 林庆华, 栗保明 2020 兵工学报 41 1697
Google Scholar
Lin Q H, Li B M 2020 Acta Armamentarii 41 1697
Google Scholar
[19] Shatoff H, Pearson D A, Kull A E 2005 IEEE Pulsed Power Conference Monterey, CA, USA, June 13–17, 2005 p253
[20] Wang G H, Xie L, He Y, Song S Y, Gao J J 2016 IEEE Trans. Plasma Sci. 44 1424
Google Scholar
[21] 王刚华, 谢龙, 赵海龙 2021 爆炸与冲击 41 111
Google Scholar
Wang G H, Xie L, Zhao H L 2021 Explosion and Shock Waves 41 111
Google Scholar
DownLoad:
Catalog
Metrics
- Abstract views: 1868
- PDF Downloads: 44
- Cited By: 0
