1、1连续发射脉冲功率电源用水冷模拟负载研究 戴宇峰 1，鲁军勇 1，张晓 1，王杰 1（1.海军工程大学 舰船综合电力技术国防科技重点实验室，湖北 武汉 430033）摘要：针对脉冲功率电源连续循环放电吸能需求，本文提出了一种 1.9MA 级循环脉冲功率水冷模拟负载方案。首先，考虑水冷负载与实际电磁发射负载的相似性要求，提出了 48 钢管阵列组成的电阻网络，可方便的实现电源不同组合方式的放电考核；其次，针对 1.9MA 级电流产生的脉冲电磁力可能引起水冷模拟负载损坏的问题，建立了水冷模拟负载三维有限元分析模型，进行了电磁力计算和结构分析，保证了水冷模拟负载稳定性；最后，利用热网络法对模拟负载在循
2、环脉冲模式下的强迫风冷、自然冷却和去离子水冷却等方式进行了温升分析，结果表明：模拟负载采用钢管内通去离子水冷却效果最好，循环放电时温度可以恢复到初始状态，最高温度 62.5，满足连续放电实验的需求。利用提出的方法设计了一台水冷模拟负载样机并进行了连续 2 次 1.9MA 放电研究，试验结果与理论分析吻合较好，负载结构运行稳定，验证了本文理论分析的正确性。关键字：电磁发射；循环脉冲功率；水冷负载；有限元法中图分类号：TM832 文献标志码：A 文章编号：Research on Water-Cooling Simulated Resistance for Continuously Emission
3、 Pulsed Power SupplyDAI Yufeng1, LU Junyong1, ZHANG Xiao1, WANG Jie1（1. National Key Laboratory of Science and Technology on Vessel Integrated Power System, Naval University of Engineering, Wuhan 430033, China）Abstract: For the demand of continuous energy absorption of pulsed power supply, the propo
4、sal about cycle pulse power water-cooling simulated resistance with 1.9 MA current levels is carried out. First, the similarities between resistance and actual electromagnetic launcher are given full consideration. Resistance network composed of 48 array of steel pipes is designed, which is convenie
5、nt for realizing discharging assessment with different power combination. Next, for the problem that 1.9MA level current can cause too large electromagnetic force, the three-dimension finite element model of resistance is constructed. Electromagnetic force calculation and structure analysis of resis
6、tance are implemented to ensure the stability of resistance. Finally, analysis on resistance temperature rise with the ways of natural, gale and deionized water cooling using thermal network is implemented. The results show that the way of deionized water cooling reaches the best performance. The re
7、sistance temperature can reach initial state when pulsed power supply continuously discharges, and the maximum temperature can reach 62.5, which can satisfy the requirement of continuous discharging. A prototype of resistance is manufactured, and is used for two times discharging with 1.9MA current.
8、 Experimental result has good agreement with theoretical analysis. Resistance structure runs well, which verifies the correctness of theoretical analysis.Keywords: electromagnetic emission; cycle pulsed power; water-cooling resistance; finite element method收稿日期：基金项目：国家自然科学基金资助项目（51522706,51407191,51
9、307176）；实验室基金（9140C840409140C84026, 9140C840409150C84358）；国家部委资助项目作者简介：戴宇峰（1988），男，安徽宿州人，博士研究生，E-mail：；鲁军勇（通信作者），男，研究员，博士，博士生导师，E-mail：电磁发射系统是利用电磁能对负载进行超高速发射的新型高效能发射系统，主要由初级储能装置、脉冲功率电源系统及发射装置构成。脉冲功率电源系统作为发射装置能量来源，为发射装置提供瞬时高功率能量 1。考虑动态发射的成本高、样本次数有限等制约因素，脉冲功率电源在进行全系统联调前，需要设计一种模拟负载，进行循环脉冲功率释放，以验证其功率输出性
10、能和运行可靠性。目前，国内外尚未见到电流为 1.9MA 级的快速连续发射脉冲电2源用模拟负载的相关报道。该模拟负载具有以下特点：一是要能保证放电特性与实际放电特性具有可比性，因此所设计的负载特性与实际发射装置特性在电气性能上如电流幅值、上升时间、脉宽方面应具有相似性；二是脉冲功率电源系统 1.9MA 级脉冲大电流对水冷模拟负载会产生强电磁力，因此必须对水冷模拟负载进行电磁力和结构分析，确保水冷模拟负载在大电流下工作稳定；三是模拟负载必须能吸收连续发射脉冲功率电源的能量，因而必须设计高效的冷却结构，确保其温升性能。本文提出了一种新结构水冷模拟负载，采用钢管作为电阻本体，钢管中通去离子水确保钢管不
11、过温，采用环氧夹板固定钢管，确保结构稳定性。通过样机的试验研究，验证了水冷模拟负载性能。1 水冷模拟负载电气性能分析电磁发射装置需要多个脉冲功率电源模块构成脉冲功率电源系统并联同步对负载放电 2-6，在设计水冷模拟负载时需要考虑电源系统同步放电时每个模块对输出电流的承受能力，防止负载不匹配造成模块输出电流过大损坏器件，因此在设计水冷模拟负载时首先要确保负载与脉冲功率电源系统的匹配性。通过前期探索，单个脉冲功率电源模块可释放 60kA 电流，负载阻值 16m，为保证每个模块同时放电时电流达到 60kA，负载阻值需要作相应调整。为满足不同规模脉冲功率电源系统放电性能考核，水冷模拟负载阻值应具有可调
12、节性，并至少能满足 2 电源模块同步放电需求。图 1 为典型 2电源模块放电示意图，通过增加并联回路单元个数，可以灵活配置脉冲功率电源系统能级，通过调节回路单元阻值，可以控制水冷模拟负载尺寸大小。本文设计的水冷模拟负载用于 32个电源模块同步放电实验，输出总电流达1.9MA，模拟负载由 16 个回路单元并联构成，单个回路单元采用 2 根钢管串联构成，回路单元阻值为 8m。CcRLd0R8mloadCcRLd0图 1 脉冲功率电源放电模型Fig.1 Discharge module of pulsed power supply0 1 2 3 4 5 6 700.511.52(MA)(ms) 图
13、2 输出电流波形对比Fig.2 Current comparison between different loads采用 SIMPLORER 仿真软件对 32 个电源模块对水冷负载装置及电磁发射装置同步放电过程进行仿真，图 2 为仿真电流波形，由图可得通过水冷负载装置与电磁发射装置电流峰值分别为 1.92MA 及 1.91MA，电流上升时间分别为 1.58ms 及 1.50ms，从能量的角度定义等效脉冲宽度为I( t)2dt/Imax2，得到不同负载装置下等效脉冲宽度分别为 1.21ms 及 1.14ms。对比发现脉冲功率电源系统在不同负载装置条件下的输出电流在电流幅值、上升时间及脉宽方面具有
14、相似性，从而确保脉冲功率电源在电磁力大小、器件过流能力及热量生成方面在两种负载装置下具有相似性。2 水冷模拟负载电磁力及结构分析水冷模拟负载由 32 根空心钢管构成 48 电阻阵列，相邻 2 根钢管在端部通过铜排连接构成电回路单元，总计 16 个回路单元。2.1 水冷模拟负载电磁力计算对于钢管阵列构成的水冷模拟负载，不同馈电方式钢管所受电磁力大小不同，因此为尽可能减小钢管所受电磁力大小，确保水冷模拟负载稳定性，需要对馈电方式进行优化分析。图 3 显示了两种水冷模拟负载馈电方式，分别为均匀馈电和交叉馈电，针对两种不同馈电方式采用解析方法计算钢管受力情况，计算时每个回路电流峰值为 120kA，相当
15、于两个电源模块并联对回路单元放电时的峰值电流。. . .(a) 均匀馈电方式(a) Model in uniform feeding way3. . .(b) 交叉馈电方式(b) Model in cross feeding way图 3 不同馈电方式示意图Fig.3 Schematic diagram of different feeding models解析法计算钢管电磁力时假设回路电流为直流，并且忽略钢管截面大小，将其等效成细线模型，即假设全部电流从钢管轴心流过 7。任意两通流导线之间相互作用力为：F=0I1I2/(2R)，其中 0 为真空磁导率，I 1、I 2为导线电流，R 为导线之间
16、的中心距离，采用叠加方法计算每根钢管所受合力。-0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6-0.6-0.5-0.4-0.3-0.2-0.100.10.2X(m)Y(m)(a) 均匀馈电方式(a) Model of uniform feeding-0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6-0.6-0.5-0.4-0.3-0.2-0.100.10.2X(m)Y(m)(b) 交叉馈电方式(b) Model of cross feeding图 4 不同馈电方式电磁力分布Fig.4 Electromagnetic force distributi
17、on with different feeding models表 1 不同馈电及计算方式电磁力对比表Tab.1 Electromagnetic force comparison between different feeding and computation ways馈电/计算方式 (0,0) (0.2,-0.15) (0.6,-0.15)均匀/解析 30.9kN 18.7kN 13.7kN均匀/有限元 30.0kN 19.2kN 14.7kN交叉/解析 13.1kN 1.33kN 0.59kN交叉/有限元 13.4kN 1.42kN 0.61kN由表 1 可以看出，解析和有限元两种计算结
18、果误差最大 2.6%，验证两种计算方法的准确性。由图 4 可以看出，两种馈电方式下钢管电磁力分布均具有对称性，钢管受力对称可以避免水冷负载装置向侧边倾斜，减小水冷模拟负载底座受力；同时，在相同位置处，钢管所受电磁力在两种馈电方式下相差加大，交叉馈电钢管最大电磁力只有均匀馈电的 42%，而最小电磁力只有均匀馈电的 4.3%，因此水冷模拟负载采用交叉馈电方式；另外，均匀馈电方式钢管主要电磁力分布在 X 轴方向，钢管受力整体均匀，而交叉馈电方式钢管主要电磁力分布在Y 轴方向，顶层和底层钢管电磁力远大于中间层钢管电磁力，而钢管受力方向不同决定着环氧夹板固定及螺栓安装方向，因此对于交叉馈电方式，环氧夹板
19、采用水平放置，同时螺栓自上而下穿过环氧夹板并与底座固定；在 Z 轴方向上，为防止钢管出现大变形，减小钢管应力，环氧夹板沿 Z 轴等间距依次排列，构成环氧夹板阵列。2.2 水冷模拟负载结构分析采用有限元仿真对水冷模拟负载进行结构分析 8-9。首先对钢管电磁场分析，获取钢管节点电磁力，然后将节点电磁力导入结构模型中进行结构分析。考虑到水冷模拟负载在 Z 轴方向尺寸较大及端部连接母排结构复杂性，导致网格剖分单元总数较大，消耗较多仿真时间及计算机资源，因此对负载结构模型进行简化。首先，忽略钢管端部连接母排对钢管形变影响，即假设钢管端部为自由端，用于模拟突出钢管形变更加严重的工况；其次，环氧夹板等间隔排
20、列，可假设任意两列环氧夹板之间钢管形变相同，因此只对两列环氧夹板之间的钢管进行结构分析，图 5 和图 6 分别为简化前后水冷负载示意图。图 5 简化前水冷负载示意图Fig.5 Graph of water-cooling load before simplificationXYZ4图 6 简化后水冷负载示意图Fig.6 Graph of water-cooling load after simplification仿真时，钢管和环氧夹板均采用六面体单元剖分，钢管截面尺寸相对环氧夹板尺寸较小，分别设置单元尺寸为 1mm 和 10mm，钢管和环氧弹性模量分别设置为 205GPa 和 25GPa，泊
21、松比均为 0.3。激励电流 120kA、钢管内外径30mm 和 38mm、钢管阵列横向和竖向间隔200mm 和 150mm、环氧夹板阵列间隔300mm，突出钢管长度 100mm。图 7、图 8 为有限元方法计算得到的钢管变形量及应力云图，其中最大变形量 0.02mm，分布在钢管最前端，最大应力 20MPa，分布在环氧夹板边缘与钢管突出部分相交处，最大应力值远远小于钢管屈服强度 200MPa，且变形量较小，满足工程应用需求。图 7 钢管变形量云图Fig.7 Deformation contour graph of load device图 8 钢管应力云图Fig.8 Stress contour
22、 graph of load device3 水冷模拟负载温度特性分析水冷模拟负载在循环脉冲大电流激励下，温度不断升高，钢管温度过高会影响钢管的结构特性和电磁特性，以及减少相邻环氧夹板的寿命，因此有必要对水冷模拟负载冷却降温。水冷模拟负载可采用自然冷却、强风冷却及去离子水冷却方式，通过对比水冷模拟负载在不同冷却方式下的冷却效果确定冷却方案。3.1 水冷模拟负载热网络法分析水冷模拟负载温度分布采用热网络法进行求解，为此作以下假设：（1）在循环脉冲电流激励下每根钢管温度分布相同，因此仅仅分析单根钢管的温度特性；（2）钢管径向尺寸较大，因此认为钢管径向温度分布相同；（3）脉冲电源放电时间为毫秒级，相
23、对放电周期较短，因此假设脉冲电源放电时刻钢管温度发生跃变；（4）钢管导热率大，钢管截面温差较小，可假设钢管截面为等温面；（5）强风冷却时钢管外表面冷却效果明显，钢管内部空气对流受限，冷却效果较差，可假设钢管内壁绝热。热网络法中采用热容、热阻及温度等效电路中的电容、电阻及电压，并采用电路求解方法求解温度分布 10-14。在强风和去离子水冷却方式下，钢管温升采用式(1)、式(2)和式(3) 计算，将式(1)中的 RW替换成 RN即为自然冷却方式下钢管温度计算公式。(1)1d()0)TWATtCRiCVte(2)()(d)FNTFNFATAiRtCNFttVe(3)(0)1,itti其中，C T表示
24、钢管热容、 RF表示冷却流体热阻、R N表示自然冷却方式下空气热阻、 RW表示强风冷却方式下空气热阻；V F表示钢管管芯冷却流体温度、V A表示钢管外侧空气温度、VXYZ5表示钢管温度、V i(0-)表示第 i 次放电后钢管温度初始值；N 表示脉冲功率电源系统循环放电次数，T 为电源模块放电周期； V 表示钢管瞬时温升，且 V=2Q/(cm)，其中， 表示负载热量产生率，Q 表示电源模块初始储能， c表示钢管比热容，m 表示回路单元钢管质量。3.2 水冷模拟负载温度变化趋势分析在不同冷却方式下，计算相应热阻和热容值分别为：CT=866.7F、R F=1.4m、R W=83.8m、R N=837
25、.7m，且 VF=25、V A=25、V(0) =25、N=25、T=10s、 =40%、Q=176kJ，计算得到钢管瞬时温升 V=37.6。0 50 100 150 200 25002004006008001000()(s)(a) 强风和自然冷却(a) Natural and gale-cooling0 50 100 150 200 250203040506070()(s)(b) 去离子水冷却(b) Deionized water-cooling图 9 不同冷却方式钢管温度Fig.9 Temperature of steel pipe in different cooling ways由图
26、9 可以看出，在强风冷却方式下，钢管温度不断上升，最终在 307左右达到稳定，对于环氧夹板 150极限工作温度，可进行连续 3 次放电；在自然冷却方式下，钢管温度不断升高，连续 25 次放电后钢管温度达到 823，在环氧夹板可承受的温度范围内，可进行连续2 次放电；在去离子水冷却方式下，在冷却水温度及流速不变的前提下，单个放电周期内，钢管温度可以恢复到初始状态，最高温度62.5，满足连续放电实验的需求，因此水冷模拟负载采用去离子水冷却方式。4 水冷模拟负载实验结果分析按照设计方案制造了包含馈电母排、钢管电阻及水冷机组的水冷模拟负载，运用该水冷模拟负载对脉冲功率电源系统放电性能考核，进行了不同能
27、级脉冲功率电源系统放电实验，输出最大峰值电流达 1.9MA，实现了对脉冲功率电源系统的充分考核，实验过程中水冷模拟负载工作稳定，输出总电流达到设计峰值。同时进行了高能级连续 2 次同步放电实验，实验过程中水冷模拟负载采用自然冷却方式，图 10为计算和实际测量钢管温度变化曲线，可以看出测量温度变化曲线与计算得到变化曲线较吻合，验证计算方法的正确性。0 200 400 600 800 1000 1200 1400 1600 1800 2000020406080100120()(s)图 10 钢管温度变化趋势Fig.10 Steel pipe temperature variety trend5 结
28、论为验证脉冲功率电源系统放电性能，设计并制造了与脉冲功率电源系统相匹配的循环功率水冷模拟负载，从电磁、结构及温度等方面分析水冷模拟负载的性能。仿真结果显示，对水冷模拟负载放电时，电源系统输出总电流幅值、上升时间及脉宽达到设计指标。电源系统高能级同步放电时水冷模拟负载结构稳定，无大变形及裂纹出现，电源系统高能级连续 2 次放电后钢管温度变化趋势与计算结果相吻合，说明计算方法的正确性，该计算方法可用于评估水冷负载装置在去离子水冷却方式下温升特性。实验结果表明，该水冷模拟负载可以配合脉冲功率电源系统进行高能级循环放电性能考核。下一步将进行电源系统高能级连续 25 次放电实验，进一步验证水冷模拟负载的
29、稳定性。参考文献（References ）1 Akiyama H, Sakugawa T, Namihira T. Industrial applications of pulsed power technologyJ. IEEE Transactions on Dielectrics and Electrical Insulation, 2007, 14(5): 1051-1064.0 50 10207012062 Bernardes J S, Stumborg M F, Jean T E. Analysis of a capacitor-based pulsed-power system
30、for driving long-range electromagnetic gunsJ. IEEE Transactions on Magnetics, 2003, 39(1): 486-489.3 Liu Y, Lin F C, Dai L, et al. Development of a compact 450-kJ pulsed-power-supply system for electromagnetic launcherJ. IEEE Transactions on Plasma Science, 2011, 39(1): 304-309.4 Jin Y S, Kim Y B, K
31、im J S, et al. Fabrication and testing of a 600-kJ pulsed power systemJ. IEEE Transactions on Plasma Science, 2013, 41(10): 2671-2674.5 Lee B, An S, Kim S H, et al. Operation of a 2.4-MJ pulsed power system for railgunJ. IEEE Transactions on Plasma Science, 2014, 42(10): 2886-2890.6 Yang K S, Kim S
32、H, Lee B, et al. Electromagnetic launch experiments using a 4.8-MJ pulsed power supplyJ. IEEE Transactions on Plasma Science, 2015, 43(5): 1358-1361.7 任人, 董志强, 国伟, 等. MA 级盘式汇流排的设计J. 高电压技术, 2014, 40(4): 1148-1152.REN Ren, DONG Zhiqiang, GUO Wei, et al. Design of disc-shaped bus bar for MA-level appli
33、cationsJ. High Voltage Engineering, 2014, 40(4): 1148-1152. (in Chinese)8 俞小华, 董健年, 王浩, 等. 电磁发射用脉冲电抗器端面电磁力计算J. 弹道学报, 2013, 25(2): 106-110.YU Xiaohua, DONG Jiannian, WANG Hao, et al. Calculation of electromagnetic force on end face of pulsed reactor for electromagnetic launchJ. Journal of Ballistics,
34、 2013, 25(2): 106-110. (in Chinese)9 杨玉东, 王建新, 薛文. 大电流汇流排电磁参数的分析与仿真J. 南京理工大学学报: 自然科学版, 2011, 35(1): 102-107.YANG Yudong, WANG Jianxin, XUE Wen. Analysis and simulation on heavy current bus bar electromagnetic parametersJ. Journal of Nanjing University of Science and Technology. 2011, 35(1): 102-107.
35、 (in Chinese)10 刘刚, 雷成华, 刘毅刚. 根据电缆表面温度推算导体温度的热路简化模型暂态误差分析J. 电网技术, 2011, 35(4): 212217.LIU Gang, LEI Chenghua, LIU Yigang. Analysis on transient error of simplified thermal circuit model for calculating conductor temperature by cable surface temperatureJ. Power System Technology, 2011, 35(4): 212217.
36、 (in Chinese)11 张尧, 周鑫, 牛海清, 等. 单芯电缆热时间常数的理论计算与试验研究J. 高电压技术, 2009, 35(11): 2801-2806.ZHANG Yao, ZHOU Xin, NIU Haiqing, et al. Theoretical calculation and experimental research on thermal time constant of single-core cablesJ. High Voltage Engineering, 2009, 35(11): 2801-2806. (in Chinese)12 雷鸣, 刘刚, 赖
37、育庭, 等. 采用 Laplace 方法的单芯电缆线芯温度动态计算J. 高电压技术 , 2010, 36(5): 1150-1154.LEI Ming, LIU Gang, LAI Yuting, et al. Dynamic calculation of core temperature of single core cables using Laplace methodJ. High Voltage Engineering, 2010, 36(5): 1150-1154. (in Chinese)13 李道奎, 段静波, 雷勇军. 金属热防护系统瞬态热分析的并联一维模型J. 国防科技大学学
38、报, 2009, 31(4): 126-130.LI Daokui, DUAN Jingbo, LEI Yongjun. Parallel one-dimensional model for transient thermal analysis of metallic thermal protection systemJ. Journal of National University of Defense Technology, 2009, 31(4): 126-130. (in Chinese)14 Liu S B. Calculation of the steady-state and transient temperature rises of round cable bundlesJ. IEEE Transactions on Power Delivery, 2010, 25(3): 1229-1235.
