Loading [MathJax]/jax/element/mml/optable/SuppMathOperators.js

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

300 MV·A脉冲发电机组向HL-2M装置环向场供电仿真研究

徐玉麟 彭建飞 卜明南 李华俊 王英翘 李维斌 李波

陶雪峰, 刘昆. 补偿脉冲发电机放电波形优化方法[J]. 强激光与粒子束, 2018, 30: 095001. doi: 10.11884/HPLPB201830.170325
引用本文: 徐玉麟, 彭建飞, 卜明南, 等. 300 MV·A脉冲发电机组向HL-2M装置环向场供电仿真研究[J]. 强激光与粒子束, 2022, 34: 065002. doi: 10.11884/HPLPB202234.210477
Tao Xuefeng, Liu Kun. Pulse shaping method for compulsator[J]. High Power Laser and Particle Beams, 2018, 30: 095001. doi: 10.11884/HPLPB201830.170325
Citation: Xu Yulin, Peng Jianfei, Bu Mingnan, et al. Simulation of 300 MV·A motor generator feeding the toroidal field power supply[J]. High Power Laser and Particle Beams, 2022, 34: 065002. doi: 10.11884/HPLPB202234.210477

300 MV·A脉冲发电机组向HL-2M装置环向场供电仿真研究

doi: 10.11884/HPLPB202234.210477
基金项目: 四川省青年科技创新研究团队项目(21CXTD108)
详细信息
    作者简介:

    徐玉麟,xuyulin@swip.ac.cn

  • 中图分类号: TM921.5

Simulation of 300 MV·A motor generator feeding the toroidal field power supply

  • 摘要: 为实现对环向场电流的高效准确控制,提高中国环流器二号M托卡马克装置的实验效率,开展了针对300 MV·A六相脉冲发电机组向环向场脉冲放电的仿真研究。从300 MV·A六相脉冲发电机组向聚变装置环向场供电系统的拓扑出发,以等效处理分析为基础建立了仿真模型并进行脉冲放电动态过程仿真。通过对比分析仿真波形与300 MV·A六相脉冲发电机组向HL-2A装置供电的实际波形,验证了仿真模型的可靠性和合理性。在此基础上进行了300 MV·A六相脉冲发电机组向HL-2M环向场供电过程仿真研究,根据不同的环向场电流平顶值,得到了其对应的300 MV·A六相脉冲发电机组励磁电压控制波形,可以作为下阶段装置高参数实验实际控制的重要参考和先进控制策略研究平台。
  • 随着现代科技的发展,越来越多的先进技术被应用到国防领域中。近年来,各国在新概念武器方面加大了投入,发展迅速。新概念武器主要包括电磁炮、高功率微波武器、高能激光武器、电热化学炮等。虽然这些高能武器的工作机理不同,但它们都需要高功率脉冲电源(PPS)的激励或驱动。脉冲电源的种类很多,它们的储能方式及放电机理不同,都有各自的特点[1]。就目前的发展水平来看,实际使用的脉冲电源系统(主要为脉冲电容器)体积大、质量大,往往超过武器系统本身,因此发展小型化、轻型化的高功率脉冲电源是新概念武器走向实战的重要保证[2]。补偿脉冲发电机(CPA)是一种特殊的同步电机,通过补偿屏蔽筒、补偿线圈等结构实现磁通压缩,能够极大地降低放电时电机的瞬态电感,从而实现大电流、高功率输出[3]。CPA采用飞轮惯性储能方式,通过机械能与电能之间的转换,能够达到极高的储能密度,且能量转换效率高,是非常有发展前景的脉冲电源方案[4]。补偿脉冲电机类型多种多样,也可驱动不同类型的武器系统。由于不同负载的阻抗特性不同,对放电波形的需求差异也很大,CPA的脉冲波形灵活调节的能力十分重要[5]。电机结构确定后,单相电枢的放电波形一般是固定的。通常将补偿脉冲电机设计成多相结构,通过各相点火角的控制及相电流组合实现总脉冲的调节[6]。以往的研究通常集中在调波方法的定性分析上,较少涉及最优化问题。文献[6]以一台四相电机为研究对象,使用遗传算法对放点波形进行优化,但未能得到很好的结果。通过建立CPA放电数学模型,在放电仿真中可对电枢转子转速、励磁电流、点火角等参数进行灵活更改,结合智能优化算法,可以找到最优控制量组合,得到优化脉冲波形。本文以一台空心补偿脉冲发电机为对象,开展两相电机脉冲波形调节的最优化问题研究。

    CPA的放电仿真方法多样,各有优缺点:场路耦合仿真方法将电机有限元模型和放电控制电路模型结合,开展联合仿真[7]。这种方法计算准确,且能够得到电磁场分布情况,但是需要消耗大量计算资源,效率较低;将电机模型进行简化,例如将被动补偿CPA简化为电压源加电感、电阻的形式,在电路仿真软件中分析其放电特性[8]。这种方法有一定局限性,对于等效电感随转子位置变化的情况不能有效模拟;建立放电过程的数学模型,通过解微分方程组的方法求解各物理量的变化。这种方法在建模阶段工作量较大,但精确的数学模型一旦建立,就能通过数值方法进行快速仿真,具有很强的灵活性[6]。本文研究的波形调节最优化问题需要大量计算,故采取建立数学模型的方法实现放电仿真。

    本文的两相四极空心CPA剖面如图 1所示,a,b两相电枢绕组正交放置,消除电磁耦合。不使用补偿屏蔽筒或补偿绕组等专门的补偿结构,放电工作时由励磁绕组实现补偿作用。

    图  1  补偿脉冲发电机剖面图
    Figure  1.  Cross section of compulsator

    由于空心电机不使用铁磁材料,需要很强的励磁电流激励。为减轻系统质量,通常采用自励磁方法。自励磁及放电原理如图 2所示。先由一个电容器为励磁绕组充电,得到种子电流,随即将电枢绕组与励磁绕组连接,形成正反馈过程,励磁电流呈指数增大。当励磁电流达到要求时,打开放电主开关为负载供电。

    图  2  自励磁及放电等效电路图
    Figure  2.  The equivalent circuit of excitation and discharge process

    电流正方向以图 2标注为准,根据基尔霍夫定律可列写方程组

    {dψadt=iaRa+(ia+ib)Rload +dLload (ia+ib)dtdψbdt=ibRb+(ia+ib)Rload +dLload (ia+ib)dtdψfdt+ifRf=0iload =|ia+ib| if: sign(iaib)>0,{dψadt=iaRa+(iaib)Rload +dLload (iaib)dtdψbdt=ibRb(iaib)Rload dLload (iaib)dtdψfdt+ifRf=0iload =|iaib| if: :sign(iaib)<0 (1)

    式中:ψ为磁链;iRL分别为电流、电阻、电感;下标a,b,f,load分别表示a, b两相电枢绕组、励磁绕组和负载;sign(·)为符号函数。

    通过磁共能Wm可计算电磁转矩

    Wm=12[iaibif][La0maf0LbmbfmafmbfLf][iaibif] (2)

    式中:mafmbf分别为a,b两相电枢绕组与励磁绕组的互感;Lf为励磁绕组的自感。

    对四极电机,电磁转矩

    T=2Wm/θ=2if(Mafiasinθ+Mbfibcosθ) (3)

    式中:q为转子的电角度,本文中为机械角度的2倍;MafMbf分别为互感mafmbf在电机旋转过程中的最大值。

    综合式(1)~(3),经整理可写出完整的状态方程

    {p[iaibif]=[(La+Lload )SLload mafSLload (Lb+Lload )mbfmafmbfLf]1[Ra+Rload +˙Lload S(Rload +˙Lload )˙maf S(Rload +˙Lload )Rb+Rload +˙Lload ˙mbf˙maf˙mbfRf][iaibif]˙θ=ω,˙ω=2TJ (4)

    式中:p表示微分算子;S=sign(ia·ib);˙maf,˙mbf为互感的变化率;˙Lload 为负载电感变化率;J为转子的转动惯量;˙θ=ω,为电角速度,˙ω为电角加速度。

    对于结构已经固定的CPA,只能通过外部控制调波,且控制量有限。但是电枢点火角、放电周期选择对单相脉冲影响很大,通过点火角控制以及多相脉冲组合是可以得到理想波形的。

    本文讨论三种典型负载:电磁轨道炮、脉冲激光器和电热化学炮。这三种负载均可用CPA驱动,但它们的负载特性以及对脉冲波形的需求差异很大。因此,有必要对负载特性进行深入分析,进而研究CPA与负载的匹配问题。

    本文电磁轨道炮属于阻感负载,其电阻、电感值与弹丸位移成正比,弹丸所受电磁推力F与负载电流的平方成正比[9],即

    Rload =Rx,Lload =Lx,F=Li2load /2 (5)

    式中:L′,R′分别为电磁炮导轨的电感梯度和电阻梯度;x为弹丸位移。

    电磁轨道炮的弹丸能够承受的最大过载(加速度)是有限制的。因此,在放电过程中尽量保持加速度均匀能够最大限度地加速弹丸,提高发射效率,达到最大炮口速度[10]。为量化脉冲波形是否适用于电磁轨道炮,通常会计算放电过程平均加速度与峰值加速度的比值[11],本文简称为“加速度比”。设弹丸加速度为a,加速过程中最大加速度为amax,放电起始和终止时间分别为t0tf,则加速度比AR的表达式为

    AR=tft0a(t)dt/[amax(tft0)] (6)

    需要说明的是,加速度比与电流幅值无关,仅仅表征电流形状。加速度比越大,则说明波形越接近平顶波。设a,b两相电枢点火角分别为θaθb,其组合决定了放电脉冲波形,加速度比也就随之确定,亦即加速度比是点火角的函数

    AR=AR(θa,θb) (7)

    这样,CPA驱动电磁轨道炮的波形优化问题可转化为如下函数优化问题

    {maxAR(θa,θb) s. t. θLj (8)

    式中:θjLθjU分别为点火角取值的上下界,由预先设定的放电周期决定。

    找到使加速度比最大的点火角组合,也就找到了最接近平顶脉冲、加速过程最平缓的放电波形。使用差分进化算法对该函数求最大值,仿真基本条件设置见表 1。令两相电枢各放电一个周期,优化后的波形如图 3所示,优化结果为:[qa, qb]=[1.546, 1.571],AR=0.513。根据优化结果可知,当两相点火角在π/2附近时,合成后的波形最接近平顶脉冲。

    表  1  仿真参数设置
    Table  1.  Simulation parameters
    initial rotor speed/(r·min-1) initial field current/kA rotor’s moment of inertia/(kg·m2) mass of projectile/g inductance gradient/(μH·m-1) resistance gradient/(mΩ·m-1) discharge period
    18 000 8 0.177 5 5 1 0.5 1
    下载: 导出CSV 
    | 显示表格
    图  3  电磁轨道炮优化放电波形
    Figure  3.  Optimized discharge pulse for EM gun

    本文电磁轨脉冲激光器是一种非线性电阻性负载,阻值在整个放电过程中呈先下降再上升的U型曲线[12],稳定放电时,有[13]

    R_{\text {load }}=K_0 / \sqrt{i_{\text {load }}} (9)

    式中:K0是脉冲氙灯的阻抗系数,与其物理结构有关。

    脉冲激光器需要尖顶脉冲驱动。虽然没有弹丸加速过程,但同样可使用“加速度比”的概念来量化波形效果,求解加速度比最小值即可。设阻抗系数K0=5时,优化结果为:[qa, qb]=[4.015, 3.142],AR=0.293。优化后的波形如图 4所示,可见,将加速度比的计算方法应用于脉冲激光器的脉冲优化是可行的,能够得到尖顶脉冲。

    图  4  脉冲激光器优化放电波形
    Figure  4.  Optimized discharge pulse for flashlamp

    电热化学炮的核心部件为毛细管,也可以看作非线性电阻。本文使用文献[14]建立的物理模型表示

    R_{\text {load }}=0.22 l \beta^{8 / 11} \alpha^{-13 / 11} i_{\text {load }}^{-6 / 11} (10)

    式中:l为毛细管长度,β为比例常数,α为毛细管截面半径。

    电热化学炮所需波形是分段上升的,很难找出一个简单的计算量来表征其优劣。文献[15]电热炮对电流的要求,即

    i_{\text {desire }}(t)=\frac{1}{2} p_0 A \sqrt{\frac{1}{m(\gamma-1)}} \sqrt{\frac{t}{R_{\text {load }}}} (11)

    式中:idesire表示需求电流,p0为炮膛加速压力,A为炮膛截面积,m为弹丸质量,γ为气体绝热常数。

    电流幅值可以通过励磁电流大小、转速等调节。为消除幅值影响,表征实际电流与需要电流波形的匹配程度,本文提出“形状方差”的概念,计算方法如下。

    (1) 设定点火角,计算实际放电电流iload(t)= iload(qa, qb),Rload(t)=Rload(qa, qb)。根据式(10)~(11)计算需求放电电流idesire(t),其中,t1tk分别对应放电初始时刻和结束时刻。

    (2) 将放电时间适当离散化,取k个时间点:t′=[t1, t2, t3, …, tk]。其中,t1tk分别对应放电初始时刻和结束时刻。

    (3) 求出各时间点处实际电流iload(tk)与需求电流idesire(tk)的比值,虚拟的统计量

    R_{\text {ld }}(k)=i_{\text {load }}\left(t_k\right) / i_{\text {desire }}\left(t_k\right) (12)

    (4) 计算电流比值的方差,定义为“形状方差”SV,可表示为

    S_V=\operatorname{var}\left(R_{\text {ld }}\right) (13)

    这样,得到的形状方差能够描述实际电流与需求电流形状的匹配程度。在相同的初始励磁电流、转速等条件下,形状方差越小,说明波形越符合要求。优化结果为:[qa, qb]=[2.005, 4.712]。优化的放电波形如图 5所示,图 6显示了需求电流与实际电流的对比情况。可见,实际电流与需求电流有一定差别,但基本符合电热化学炮的要求。

    图  5  电热化学炮优化放电波形
    Figure  5.  Optimized discharge pulse for ETC
    图  6  需求电流与实际电流对比
    Figure  6.  Comparison of desired and actual current

    本文建立了两相空心补偿脉冲发电机放电过程的数学模型,在此基础上开展了放电模型最优化研究。针对三种不同负载,分析负载特性及需求波形,提出相应的优化指标。优化结果表明,采用加速度比作为优化指标,可以有效得到适合电磁轨道炮和脉冲激光器的放电脉冲。对于电热化学炮,本文定义的“形状方差”可较好地对放电脉冲进行优化。基于电机数学模型的仿真具有很高的灵活性。但是,文中建立的数学模型是否准确,提出的波形优化方法是否合理,则需要进一步的实验验证。根据实验结果对仿真模型进行修正,可以得到更加精确的模型。两相电机具有较强的波形调节能力。可以推论,随着相数进一步增加,电机放电波形调节的灵活性将进一步增强。

  • 图  1  机组向装置环向场供电原理图

    Figure  1.  Schematic diagram of the MG feeding the TF power supply of fusion device

    图  2  励磁控制系统工作拓扑图

    Figure  2.  Working topology diagram of excitation control system

    图  3  300MV·A脉冲发电机组向环向场供电仿真模型

    Figure  3.  Simulation model of a 300 MV·A MG feeding the TF power supply

    图  4  脉冲发电机组励磁电压波形

    Figure  4.  MG excitation voltage waveform

    图  5  机组向环向场供电的对比图

    Figure  5.  Comparison between simulation results and experimental waveforms

    图  6  环向场电流的对比图

    Figure  6.  Contrast diagram of the TF current

    图  7  环向场电流仿真波形

    Figure  7.  Simulation waveforms of the TF current

    图  8  励磁电压控制波形(初始转速400 r/min)

    Figure  8.  Excitation voltage control waveforms (initial speed 400 r/min)

    图  9  励磁电压控制波形(初始转速500 r/min)

    Figure  9.  Excitation voltage control waveforms (initial speed 500 r/min)

    图  1  4.5 kW窄线宽线偏振光纤激光器结构示意图

    图  2  实验结果图

    表  1  同步发电机的主要参数

    Table  1.   Main parameters of synchronous generator

    rated
    voltage/V
    rated
    current/A
    connecting
    mode
    rated speed/
    (r/min)
    rated
    frequency/Hz
    rotational
    inertia/(kg·m2)
    rated excitation
    voltage/V
    rated excitation
    current/A
    300029 000dual Y shift 30°335~50067~1007 300 0003492~4093301~353
    下载: 导出CSV

    表  2  环向场的主要参数

    Table  2.   Main parameters of the toroidal field

    maximum magnetic
    induction/T
    inductance/mHresistance/mΩtime constant/smaximum flat-top
    time/s
    TF of HL-2A2.82305845
    TF of HL-2M332487
    下载: 导出CSV
  • [1] 刘永, 李强, HL-2M研制团队. 中国环流器二号M(HL-2M)托卡马克主机研制进展[J]. 中国核电, 2020, 13(6):747-752. (Liu Yong, Li Qiang, HL-2M Development Team. The development progress of the tokamak machine for HL-2M[J]. China Nuclear Power, 2020, 13(6): 747-752
    [2] 宣伟民, 彭建飞, 李华俊, 等. 300MVA脉冲发电机系统研制进展[C]//2013核工业西南物理研究院年报. 2014: 98-100

    Xuan Weimin, Peng Jianfei, Li Huajun, et al. Development progress of 300MVA motor generator[C]//Southwestern Institute of Physics, Annual Report 2013 of Southwestern Institute of Physics. 2013: 98-100
    [3] 彭建飞, 王树锦, 李华俊, 等. 基于Matlab/Simulink的环向磁场脉冲电源建模与仿真[J]. 核聚变与等离子体物理, 2004, 24(4):254-258. (Peng Jianfei, Wang Shujin, Li Huajun, et al. Modeling and simulation of toroidal field power supply system of the HL-2A tokamak based on Matlab/Simulink[J]. Nuclear Fusion and Plasma Physics, 2004, 24(4): 254-258 doi: 10.3969/j.issn.0254-6086.2004.04.003
    [4] 李华俊, 徐丽荣, 刘学梅, 等. HL-2A装置环向场电源系统及其控制[J]. 兵工自动化, 2004, 23(1):55-57. (Li Huajun, Xu Lirong, Liu Xuemei, et al. Power supply system and its control of HL-2A toroidal field[J]. Ordnance Industry Automation, 2004, 23(1): 55-57 doi: 10.3969/j.issn.1006-1576.2004.01.022
    [5] 陈宇红, 王尚柄. HL-2A装置环向场电源系统设计[J]. 中国核科技报告, 1999(s2):52-55. (Chen Yuhong, Wang Shangbing. Design of HL-2A toroidal field power supply system[J]. China Nuclear Science and Technology Report, 1999(s2): 52-55
    [6] 康莉, 王树锦. HL-2A装置环向场电源系统的EMTP仿真[J]. 电工技术, 2003(6):14-15. (Kang Li, Wang Shujing. EMTP simulation of HL-2A unit circular field power supply system[J]. Electric Engineering, 2003(6): 14-15 doi: 10.3969/j.issn.1002-1388.2003.06.006
    [7] Peng Jianfei, Xuan Weimin, Wang Haibing, et al. Study on matching a 300 MVA motor generator with an ohmic heating power supply in HL-2M[J]. Plasma Science and Technology, 2013, 15(3): 300-302. doi: 10.1088/1009-0630/15/3/22
    [8] 彭建飞, 宣伟民, 将宝钢, 等. HL-2M装置300MVA脉冲发电机组的设计[J]. 核聚变与等离子体物理, 2021, 41(s2):443-448. (Peng Jianfei, Xuan Weimin, Jiang Baogang, et al. Design of 300MVA motor generator on HL-2M tokamak[J]. Nuclear Fusion and Plasma Physics, 2021, 41(s2): 443-448 doi: 10.16568/j.0254-6086.2021s2002
    [9] 袁保山, 姜韶风, 陆志鸿. 托卡马克装置工程基础[M]. 北京: 原子能出版社, 2011

    Yuan Baoshan, Jiang Shaofeng, Lu Zhihong. Engineering foundaments of tokamak devices[M]. Beijing: Atomic Energy Press, 2011
    [10] 宣伟民, 王英翘, 李华俊, 等. HL-2M装置供电系统研制[J]. 核聚变与等离子体物理, 2021, 41(s2):437-442. (Xuan Weimin, Wang Yingqiao, Li Huajun, et al. Development of HL-2M power supply system[J]. Nuclear Fusion and Plasma Physics, 2021, 41(s2): 437-442 doi: 10.16568/j.0254-6086.2021s2001
  • 期刊类型引用(5)

    1. 杨颜冰,杨文杰,吴巍,冯文天,欧贤金,李大为,杨静,张翔,陈广全,张旭,姚庆高,马力祯. 张力线磁场测量系统的平台搭建与相关测试实验研究. 原子核物理评论. 2023(03): 377-384 . 百度学术
    2. 马娜,李波,董岚,王小龙,门玲鸰,王铜,梁静,何振强,柯志勇,卢尚,韩圆颖,闫路平,张露彦,闫皓月. CSNS四极铁标定方案分析. 强激光与粒子束. 2022(08): 97-102 . 本站查看
    3. 龙虎,彭志勇. 海量激光超声数据采集与处理系统设计. 激光杂志. 2020(05): 153-157 . 百度学术
    4. 张晓芳. 基于机器学习的光纤网络激光器异常功率数据采集系统设计. 激光杂志. 2020(09): 111-115 . 百度学术
    5. 刘博,门玲鸰,辛洪兵,董岚,李波,王小龙,罗涛,梁静,王铜,马娜,何振强,柯志勇,张晓辉,沈建新. HEPS预准直单元BQ磁铁精密高稳垫锁位置调整装置设计. 核技术. 2019(11): 37-41 . 百度学术

    其他类型引用(2)

  • 加载中
图(11) / 表(2)
计量
  • 文章访问数:  665
  • HTML全文浏览量:  195
  • PDF下载量:  42
  • 被引次数: 7
出版历程
  • 收稿日期:  2021-12-13
  • 修回日期:  2022-03-24
  • 录用日期:  2022-04-01
  • 网络出版日期:  2022-04-09
  • 刊出日期:  2022-06-15

目录

/

返回文章
返回