高压SF6断路器开断零区电弧特性仿真分析

Simulation Analysis of Arc Characteristics on Breaking Zero Zone in High Voltage SFCircuit Breaker

王飞鸣1,2, 赵义松1,2, 史可鉴1,3, 李斌1,3, 张彬1,3, 徐振乙1,2

1. 国网辽宁省电力有限公司电力科学研究院,沈阳110006

2. 辽宁东科电力也有限公司,沈阳110006

3. 国网高电压强电流实验室,沈阳110006

WANG Feiming1,2, ZHAO Yisong1,2, SHI Kejian1,3, LI Bin1,3, ZHANG Bin1,3, XU Zhenyi1,2

1. Research Institute of State Grid Liaoning Electric Power Company, Shenyang 110006, China
2. Liaoning Dongke Power Co., Ltd., Shenyang 110006, China
3. High Voltage and Large Current Laboratory of State Grid Corporation, Shenyang 110006, China

  • 王飞鸣(通信作者)1986—,男,博士,主要从事高电压绝缘技术和高压断路器瞬态特性的研究工作E-mail: famingehao@126.com

摘要

为了分析高压SF6断路器灭弧室内零区电弧等离子体与气流分子间能量运输物理过程,寻求最佳的小电流电弧开断控制策略,需要开展断路器开断零区电弧特性的仿真计算分析。建立了断路器开断电弧能量多物理场耦合计算模型,并计算了零区电弧等离子体密度及热力学特性参数的空间分布;然后采用多物理状态参数综合评判的方法,分析了不同开断电流、燃弧时间、开断速度等因素对零区电弧与气流分子间能量运输过程的影响,并求得了零区电弧等离子体物理特性变化规律及参数范围。计算结果表明:气流与电弧间的对流效应明显;当弧道出现温度小于4 000 K、密度大于3 kg/m3、焓值小于1.6×107 J/kg、电导率小于800 S/m的绝缘断点时,电弧熄灭。研究结果可以为高压SF6断路器灭弧室设计、操作机械特性选择、电弧熄灭过程物理特性分析提供重要的理论基础。

关键词 : 断路器; SF6; 零区; 等离子体; 电弧特性;

DOI:10.13336/j.1003-6520.hve.20190430009 2019年5月31日第45卷May

ABSTRACT

In order to analyze the physical process between zero zone arc plasma and airflow intermolecular energy transport in the high voltage SF6 circuit breaker, and find the optimal control strategy of small current arc switching, the simulation analysis of arc characteristics on breaking zero zone is presented. The coupled field calculation physics of the breaking arc energy in the circuit breaker was established to calculate the spatial distribution of the zero zone arc plasma density and thermodynamic characteristics parameters. Using the comprehensive evaluation method of multiple physical state parameters, the effects of different breaking current, arcing time, breaking speed and other factors on the energy transport process between zero zone arc and gas flow were analyzed, the variation law and parameter range of the zero zone arc plasma physical properties were obtained. The calculation results show that there is a significant convection effect between the airflow and the arc. When the insulation break point with a temperature less than 4 000 K, a density greater than 3 kg/m3, enthalpy value less than 1.6×107 J/kg, and a conductivity less than 800 S/m, the arc extinguishes. The research results can provide important theoretical basis for the design of the arc-extinguishing chamber of the high-voltage SF6 circuit breaker, the selection of the mechanical characteristics of the operation, and the physical characteristics of the arc-extinguishing process.

KEY WORDS : circuit breaker; SF6; zero zone; plasma; arc characteristics;

0 引言

高压SF6断路器是电力系统重要的保护及控制单元,断路器能否可靠地开断故障电流直接关系着电力系统的稳定运行[1-2]。断路器开断过程主要分为两个阶段,燃弧阶段和绝缘介质恢复阶段[3-4]。断路器电弧零区即故障电流过零区间,是电弧熄灭和介质恢复过程的过渡过程。在弧道电流过零区间,系统电流为开断电弧提供的功率逐渐减小,在断路器强烈气吹的作用下,压气缸内大量的冷态气体分子与电弧等离子体进行能量交换,带走了足够多的粒子热动能。随着等离子体内电子和重粒子的热动能下降,电子和重粒子温度下降,粒子间的碰撞频率和电离率下降,直到电子密度下降到不足以维持弧道的导通,电弧熄灭[5-7]。电弧熄灭后残余电子与正离子复合为原子、分子,带电粒子密度下降,中性粒子密度上升,介质绝缘强度快速恢复,弧道进入介质恢复阶段。此时触头间的系统瞬态恢复电压(transient recovery voltage,TRV)快速升高,当弧后介质绝缘强度恢复速度大于触头间TRV上升速度时,不会发生弧后重击穿现象,断路器成功开断。

断路器零区电弧与气吹分子间的热动能交换是开断电弧熄灭的关键因素,这主要取决于电弧过零时刻的能量分布、气吹速度、弧道等离子体与气流分子间的热力学特性参数空间变化等。断路器开断工况不同,使得电弧熄灭的物理过程不同。电力系统短路故障工况下断路器开断电流最高可达80 kA,断路器必须具备强烈的气吹条件及很高的开断速度才能在电流零区熄灭电弧[8-10]。但是针对容性、感性负载电流开断工况,过于强烈的气吹条件和快速的开断速度会造成电弧提前熄灭,直接导致熄弧开距小,介质绝缘距离不够,极易产生弧后重击穿现象。目前,断路器多采用弹簧、液压、永磁操动机构[11-12],只能用一条操动速度特性曲线应对不同开断条件的需求,在保证短路故障开断的前提下,很难保证容性、感性负荷开断的可靠性。随着断路器智能操动机构的不断发展,电机操动机构以其高可控性使得断路器智能控制变为可能[13-16],但是针对不同的开断工况选取合适的操动速度特性,保证最佳的熄弧特性及介质恢复特性成为断路器智能控制策略亟需解决的核心问题。

本文基于文献[17-18]计算得出的SF6气体电弧等离子体热力学参数,针对容性负载电流开断工况,建立了断路器开断电弧能量、等离子体密度、灭弧室压强和气流速度的耦合计算物理模型,得出零区电弧等离子体密度及热力学特性参数的空间分布,分析了不同开断电流、相角、速度对零区电弧与气流分子间能量输运过程的影响,给出零区电弧等离子体热力学特性变化规律及参数范围,最后提出了符合容性小电流开断工况下的断路器控制策略。研究结果是高压SF6断路器灭弧室设计、操作机械特性选择、电弧熄灭过程物理特性分析的重要理论基础。

1 仿真计算模型

本文将文献[17-18]中计算得出的等离子体物理特性参数结果带入到气流场k-ε和Navier-Stokes方程中,结合SF6断路器灭弧室的结构,应用FLUENT气流场计算软件中的自定义函数单元,自动计算电弧等离子体在不同温度及压强下电导率和电弧能量,实现电弧能量根据灭弧室内流场变化的空间分布加载[19-21]。针对特高压无功补偿电容器组专用126 kV SF6断路器开断容性负荷AC 1 600 A电流的工况,建立断路器弧触头、灭弧室、压气缸的物理计算模型,断路器灭弧室模型图如图1所示。

设置断路器开断距离为150 mm,弧触头运行行程为200 mm,开断速度为9.6 m/s,灭弧室内气流初始速度为0,初始温度为300 K,初始压强为0.7 MPa。灭弧室结构采用三角形剖分网格结构,压气缸采用四面体网格结构,通过控制不同位置单元网格大小降低网格畸变率,提高气流场的计算精度,断路器网格结构如图2所示。采用动态层铺的动网格计算模型表征断路器压气缸运动和压气物理过程,分别设置固体壁面、对称轴、压力入口、压力出口4种边界条件,具体设置条件及参数在文献[20-24]中详细给出。

断路器电弧能量和开断速度是影响断路器压气过程、气流特性和电弧熄灭特性的重要因素。电弧能量通过弧道电流值和开断电流相角来改变,通过分析不同开断速度和电弧能量下灭弧室内零区电弧温度、压强、气流、密度等参数(电弧参数取弧道中心区域平均值)的变化分析开断电弧的熄弧特性。

图1 灭弧室结构图 Fig.1 Schematic diagram of circuit breaker interrupter

图2 灭弧室剖分网格图 Fig.2 Circuit breaker interrupter meshed structure diagram

2 不同电流电弧熄弧特性分析

计算电弧电流开断相角为0.5π,电流值分别为交流6 400、3 200、1 600和800 A。燃弧时间、熄弧开距与开断电流的关系如图3所示,不同开断电流下的气流参数和灭弧室内参数分布情况分别如表1和图4所示。

图3中电弧燃弧时间及熄弧时刻的开距值随着开断电流值的增大而增大。由于开断800 A到6 400 A电流在电弧过零区间存在提前熄弧现象,且随着弧道电流值的减小,过零前提前熄弧现象越明显。因此,当开断速度和相位相同时,开断电流越大,燃弧时间相对越长,熄弧开距值越大。

结合表1和图4电弧及气流参数的空间分布可以看出,0.5π的开断相角下,燃弧时间并没有达到5 ms。由于弧道电流只有几百到几千A,弧道能量较小,在电流第一个过零区间内,小喷口打开后,电弧等离子体进入绝缘介质恢复过程。在电弧熄灭时刻,弧道内最低温度下降到4 000 K以下,此时弧道内带电粒子快速与电子附着为原子、分子等中性粒子,电子密度下降,介质密度快速上升,弧道内的介质密度最低值上升到3 kg/m3以上,弧触头间介质绝缘强度快速恢复。在此过程中,小喷口打开后断路器的压气过程对电弧等离子体能量的扩散作用明显。小喷口打开前,在压气缸的运动及电弧能量的共同作用下,小喷口上下游区域建立起明显的压强差,且开断电流值越大,压强差值越大(Δp=0.189~0.295 MPa)。小喷口打开后,在压强差的作用下,弧触头间产生明显的气流,不同开断电流下气流速度为0.42~0.69 Ma,电弧等离子体内粒子运动的热动能被气流快速带走。以开断电流1 600 A为例可以看出:(1)弧道内等离子体热导率为0.84 W·(m·K)-1,气流内热导率为1.32 W·(m·K)-1,气流的导热能力明显大于弧道内的等离子体,说明在气流与电弧间的能量交换明显;(2)弧道内等离子体粘性系数为0.068 g·(m·s)-1,气流内粘性系数为0.24 g·(m·s)-1,气流的动量传导能力明显大于弧道内的等离子体,说明在气流与电弧间的动量交换明显;(3)弧道内等离子体焓值为1.56×107 J/kg,气流内焓值为2.46×107 J/kg,气流带走等离子体内大量的热能,降低弧道温度,弧道内电导率下降到672 S/m,介质绝缘性能快速恢复,出现绝缘断点,弧道断裂,电弧熄灭。

图3 不同开断电流燃弧时间与熄弧开距 Fig.3 Arcing time and open distance arc extinguishing under different breaking currents

表1 不同开断电流下的气流参数 Table 1 Airflow parameters under different breaking currents

3 不同开断相角电弧熄弧特性分析

电弧电流开断相角θ分别为0.8π、0.7π、0.6π和0.5π,开断电流为I=1.414×1 600×sin(ωt+θ),开断速度为9.6 m/s时,开断电流、熄弧开距和燃弧时间如图5所示。

图5可以看出,开断相角θ是影响熄弧开距lx和燃弧时间Ts的关键因素,直接影响开断电弧能量,针对容性小电流开断工况,合理选择开断相角是保证开断可靠性的关键。随着开断相角θ由0.8π到0.5π变化,电弧燃弧时间变大,熄弧时刻弧触头开距值变大,弧触头刚分时刻电流值变大,因此燃弧能量值越大。

不同开断相角下的气流参数和灭弧室内参数分布情况分别如表2和图6所示。可知,燃弧时间由1.34 ms到4.27 ms变化,气吹速度速度由0.3 Ma增加到0.62 Ma,气流焓值由2.11×107 J/kg增加到2.46×107 J/kg,气流热导率由1.08 W/(m·K)增加到1.32 W/(m·K),气流粘性系数由0.19 g·(m·s)-1增加到0.24 g·(m·s)-1。说明燃弧时间越大,电弧上下游

图4 不同开断电流下灭弧室内参数分布 Fig.4 Arc parameter distribution of interrupter under different breaking currents

图5 不同开断相角下开断电流、熄弧开距与燃弧时间 Fig.5 Breaking current, arcing time and open distance arc extinguishing under different phase angles

区域压强差值越大,气吹速度越大,气流与电弧等离子体间的能量交换越明显。燃弧时间越长,需要断路器以更强烈的气吹强度带走更多的能量,直到熄弧时刻弧道温度<4 000 K,弧道密度>3 kg/m3,弧道焓值<1.6×107 J/kg,弧道电导率<800 S/m,弧道出现绝缘断点,电弧熄灭,介质绝缘特性开始恢复。

4 不同开断速度电弧熄弧特性分析

本文设定断路器开断电流为1 600 A,开断相角为0.5π,开断速度分别为9.6、8.64、7.68、6.72、5.76和4.8 m/s。给出不同开断速度下熄弧时刻灭弧室内电弧参数分布,具体参数如表3和图7—图10所示。

通过分析可以看出,随着开断速度由4.8 m/s增大到9.6 m/s,弧触头刚分时间由33.07 ms下降到16.53 ms,熄弧时刻的开距值由17 mm增大到41 mm,燃弧时间维持在4.3 ms左右变化。由此可知,开断速度直接影响熄弧开距值,对燃弧时间影响不大。开断速度增大,导致熄弧时刻气流速度由0.46马赫增加到0.65 Ma,气流焓值由2.21×107 J/kg增大到2.46×107 J/kg,气流热导率由1.03 W/(m·K)-1增大到1.323 W/(m·K)-1,气流粘性系数由0.215 g/(m·s)-1增大到0.24 g/(m·s)-1。气流速度越快,气流分子与等离子体间的碰撞频率越大,增大了气流与等离子体间的能量传导,因此气流的热导率、粘性系数及焓值与开断速度成正比关系。

针对断路器开断容性负载的工况,由于负荷电流仅为1 600 A,断路器在开断电流第一个过零点即可将电弧熄灭。短燃弧时间下,熄弧开距值小,弧触头间的介质绝缘距离小,灭弧室内电场值较大。在弧后剩余能量的共同作用下,弧触头间介质绝缘

表2 不同开断相角下的气流参数 Table 2 Airflow parameters under different breaking phase angles

表3 不同开断速度下关键参数 Table 3 Key parameters under different breaking speeds

强度恢复速度较慢。当介质绝缘强度恢复电压小于弧触头间瞬态恢复电压时,发生弧后重击穿现象,导致断路器开断失败。所以合理地控制熄弧开距及气流速度是决定容性负荷电流电弧熄灭及弧后介质恢复速度的主要因素。通过本文的计算分析可知,断路器开断相角决定燃弧时间,燃弧时间与开断速度决定熄弧开距值,电弧能量与开断速度决定吹弧气流速度,因此,合理地选择开断相角及开断速度是保证开断特性的关键因素。本文结合计算结果,

图6 不同相角开断灭弧室内参数分布 Fig.6 Arc parameter distribution of interrupter under different breaking phase angles

图7 不同速度开断灭弧室内温度和压强分布 Fig.7 Temperature and pressure distribution of interrupter under different speeds

图8 不同速度开断灭弧室内焓值和电导率分布 Fig.8 Enthalpy and conductivity distribution of interrupter under different speeds

建议1 600 A容性电流开断燃弧时间在2 ms以上,熄弧开距值在40 mm以上。因此,选取开断速度为9.6 m/s,开断相角在0~0.8π范围内。

5 结论

1)本文通过仿真计算的方法给出高压SF6断路器零区电弧熄弧过程物理特性参数空间分布及变化规律,定量分析开断电流、相角及速度等宏观参量对电弧热导率、粘性系数、电导率、焓值等微观参数的影响权重,给出气流参数对电弧熄灭物理过程的宏观描述。

2)熄弧过程等离子体参数的定量变化规律为:当弧道温度<4 000 K,弧道密度>3 kg/m3,弧道焓值<1.6×107 J/kg,弧道电导率<800 S/m时,弧道出现绝缘断点,电弧熄灭,介质绝缘特性开始恢复。

图9 不同速度开断灭弧室内密度和马赫数分布 Fig.9 Density and Mach number distribution of interrupter under different speeds

图10 不同速度开断灭弧室内热导率和粘性系数分布 Fig.10 Thermal conductivity and viscosity coefficient distribution of interrupter under different speeds

3)针对126 kV SF6断路器容性负载1 600 A电流开断工况,给出合理的燃弧时间及熄弧开距控制策略,即选取最佳的开断相角及速度,具体表现为:1 600 A容性电流开断燃弧时间>2 ms,熄弧开距值>40 mm,则选取开断速度为9.6 m/s,开断相角在0~0.8π范围内。

编辑 曹昭君 余洋洋

参考文献

[1] 林莘. 现代高压电器技术[M]. 2版. 北京:机械工业出版社,2006:323-334. LIN Xin.The modern high voltage electrical technology[M]. 2nd ed. Beijing, China: Mechanical Industry Press, 2006: 323-334.

[2] 刘振亚. 特高压交流输电系统过电压与绝缘配合[M]. 北京:中国电力出版社,2008:33-37. LIU Zhenya.Over-voltage and insulation coordination of UHV AC transmission system[M]. Beijing, China: China Electric Power Publishing House, 2008: 33-37.

[3] 傅中,陈维江,李志兵,. SF6断路器动态接触电阻测量方法及影响因素[J]. 高电压技术,2017,43(5):1535-1542. FU Zhong, CHEN Weijiang, LI Zhibing, et al.Measurement method and influencing factors of dynamic contact resistance of SF6 circuit breaker[J]. High Voltage Engineering, 2017, 43(5): 1535-1542.

[4] 宋冬冬,程林,林志法,. 电弧热等离子体建模、仿真及应用综述[J]. 高电压技术,2018,44(3):932-943. SONG Dongdong, CHENG Lin, LIN Zhifa, et al.Review of arc thermal plasma modeling, simulation and application[J]. High Voltage Engineering, 2018, 44(3): 932-943.

[5] 王飞鸣. 高压SF6断路器非平衡态电弧机理及开断特性研究[D]. 沈阳:沈阳工业大学,2016. WANG Feiming.Research on mechanism of non-equilibrium plasma arc and breaking characteristics in high-voltage SF6 circuit breakers[D]. Shenyang, China: Shenyang University of Technology, 2016.

[6] 荣命哲,刘定新,李美,. 非平衡态等离子体的仿真研究现状与新进展[J]. 电工技术学报,2014,29(6):271-282. RONG Mingzhe, LIU Dingxin, LI Mei, et al.Research status and new progress on the numerical simulation of non-equilibrium plasmas[J]. Transactions of China Electrotechnical Society, 2014, 29(6): 271-282

[7] 李和平,于达仁,孙文廷,. 大气压放电等离子体研究进展综述[J]. 高电压技术,2016,42(12):3697-3727. LI Heping, YU Daren, SUN Wenting, et al.State-of-the-art of atmospheric discharge plasmas[J]. High Voltage Engineering, 2016, 42(12): 3697-3727.

[8] 王瑞雪,海彬,田思理,. 绝缘材料表面电荷测量优化及等离子体处理对其表面电特性的影响[J]. 高电压技术,2017,43(6):1808-1815. WANG Ruixue, HAI Bin, TIAN Sili, et al.Optimization of dielectric material surface charge measurement and impact of plasma treatment on their surface electrical characteristics[J]. High Voltage Engineering, 2017, 43(6): 1808-1815.

[9] 王伟宗,荣命哲,YAN J D,. 高压断路器SF6电弧电流零区动态特征和衰减行为的研究综述[J]. 中国电机工程学报,2015,35(8):2059-2072. WANG Weizong, RONG Mingzhe, YAN J D, et al.Investigation of the dynamic characteristics and decaying behavior of SF6 arcs in high voltage circuit breakers during current-zero period: a review[J]. Proceedings of the CSEE, 2015, 35(8): 2059-2072.

[10] 林莘,王飞鸣,刘卫东,. 高压SF6断路器关合预击穿特性计算与试验研究[J]. 中国电机工程学报,2016,36(5):1445-1452. LIN Xin, WANG Feiming, LIU Weidong, et al.Research on calculation and experiment of closing pre-strike characteristic in high voltage SF6 circuit breaker[J]. Proceedings of the CSEE, 2016, 36(5): 1445-1452.

[11] 王飞鸣,林莘,谢寅志,. 特高压半封闭气体绝缘组合电器的暂态电压仿真分析[J]. 高电压技术,2012,38(2):407-413. WANG Feiming, LIN Xin, XIE Yinzhi, et al.Transient voltage simulation analysis of UHV semi-enclosed gas insulated composites[J]. High Voltage Engineering, 2012, 38(2): 407-413.

[12] 史可鉴. GIS中隔离开关电机操动机构及控制方法研究[D]. 沈阳:沈阳工业大学,2016. SHI Kejian.Research on the operating mechanism and control method of isolating switch motor in GIS[D]. Shenyang, China: Shenyang University of Technology, 2016.

[13] 林莘,王飞鸣,冯凯,. 高压SF6断路器介质绝缘强度恢复机理研究[J]. 中国电机工程学报,2017,37(20):6118-6125. LIN Xin, WANG Feiming, FENG Kai, et al.Study on recovery mechanism of dielectric insulation strength of high voltage SF6 circuit breaker[J]. Proceedings of the CSEE, 2017, 37(20): 6118-6125.

[14] 陈坚,向金秋,郭恒,. 基于一维电子平衡流体模型的平行板离子引出数值模拟[J]. 高电压技术,2017,43(6):1830-1836. CHEN Jian, XIANG Jinqiu, GUO Heng, et al.One-dimension modeling on parallel-plate ion extraction process based on electron-equilibrium fluid model[J]. High Voltage Engineering, 2017, 43(6): 1830-1836.

[15] WANG W Z, WU Y, RONG M Z.Influence of ablated PTFE vapor entrainment on critical dielectric strength of hot SF6 gas[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2014, 21(4): 1478-1485.

[16] WANG W Z.Investigation of the dynamic characteristics and decaying behavior of SF6 arcs in switching applications[D]. Liverpool, UK: The University of Liverpool, 2013.

[17] LESLIE S F, RICHARD W L.Composition and transport properties of SF6 and their use in a simplified enthalpy flow arc model[J]. Proceedings of the IEEE, 1971, 59(4): 474-485.

[18] CHERVY B, GLEIZES A, RAZAFINIMANANA M.Thermodynamic properties and transport coefficients in SF6-Cu mixtures at temperatures of 300-30000 K and pressures of 0.1-1 MPa[J]. Journal of Physics D: Applied Physics, 1994, 27(6): 1193.

[19] 林莘,王亮,王飞鸣,. 非平衡态双温度SF6等离子体电弧数学模型研究[J]. 高压电器,2015,51(3):1-7. LIN Xin, WANG Liang, WANG Feiming, et al.Non-equilibrium dual temperature SF6 plasma arc mathematical model research[J]. High Voltage Apparatus, 2015, 51(3): 1-7.

[20] 王飞鸣,林莘,徐建源. 高压SF6断路器非平衡态等离子体电弧的熄弧特性[J]. 高电压技术,2014,40(10):3073-3084. WANG Feiming, LIN Xin, XU Jianyuan.Extinction characteristics of non-equilibrium plasma arc in high-voltage SF6 circuit breaker[J]. High Voltage Engineering, 2014, 40(10): 3073-3084.

[21] 王亮. SF6断路器开断容性小电流介质恢复特性研究[D]. 沈阳:沈阳工业大学,2014. WANG Liang.Research on dielectric recovery characteristic of small capacitive current in SF6 circuit breaker[D]. Shenyang, China: Shenyang University of Technology, 2014.

[22] 林莘,王飞鸣,刘卫东,. 高压SF6断路器冷态介质恢复特性试验研究[J]. 中国电机工程学报,2016,36(14):3973-3981. LIN Xin, WANG Feiming, LIU Weidong, et al.Experiment researches on cold state dielectric recovery characteristics in high voltage SF6 circuit breakers[J]. Proceedings of the CSEE, 2016, 36(14): 3973-3981.

[23] LIN X, WANG F M, XU J Y, et al.Study on the mathematical model of dielectric recovery characteristics in high voltage SF6 circuit breaker[J]. Plasma Science and Technology, 2016, 18(3): 223-229.

[24] 徐建源,刘振祥,夏亚龙,. 1 100 kV GIL快速接地开关分闸过程C4F7N/CO2混合气体介质绝缘强度计算分析[J]. 高电压技术,2018,44(10):3125-3134. XU Jianyuan, LIU Zhenxiang, XIA Yalong, et al.Calculation and analysis on dielectric strength of C4F7N/CO2 in 1 100 kV GIL HSGS during the opening process[J]. High Voltage Engineering, 2018, 44(10): 3125-3134.

  • 目录

    图1