系统级HEMP耦合分析方法研究进展.pdf

系统级H EMP耦合分析方法研究进展谢海燕(西北核技术研究所,西安 7 1 0 0 2 4)摘 要:高空电磁脉冲(h i g h a l t i t u d e e l e c t r o m a g n e t i c p u l s e,HEMP)幅值强、频谱宽、覆盖范围广,可对大面积范围内的绝大多数未加固的电子电力系统造成损伤。系统级HEMP耦合分析是HEMP评估和加固防御的基础。本文介绍了作者团队近1 0年来在系统级HEMP耦合分析方法方面的研究进展,主要包括两方面:一是针对局域多尺度多元复杂系统HEMP耦合分析的瞬态电磁拓扑方法;二是针对广域分布式系统HEMP耦合分析的,基于有效耦合长度的网络耦合分析方法。关键词:高空电磁脉冲;电磁拓扑;有效耦合长度;耦合分析;局域系统;广域系统中图分类号:O 4 4 1.4;T L 9 1文献标志码:A D O I:1 0.1 2 0 6 1/j.i s s n.2 0 9 5 6 2 2 3.2 0 2 3.0 2 0 1 0 2收稿日期:2 0 2 2 0 9 2 1;修回日期:2 0 2 3 0 4 2 3基金项目:国家重点研发计划资助项目(2 0 2 0 Y F A 0 7 0 9 8 0 0)作者简介:谢海燕(1 9 8 4-),女,安徽黄山人,研究员,博士,主要从事强电磁脉冲耦合效应研究。E-m a i l:x i e h a i y a n n i n i t.a c.c nR e s e a r c h P r o g r e s s o f S y s t e m L e v e l H EMP C o u p l i n g A n a l y s i s M e t h o d sX I E H a i y a n N o r t h w e s t I n s t i t u t e o f N u c l e a r T e c h n o l o g y X i a n 7 1 0 0 2 4 C h i n a A b s t r a c t H i g h-a l t i t u d e e l e c t r o m a g n e t i c p u l s e H E M P w i t h h i g h m a g n i t u d e w i d e b a n d w i d t h a n d l a r g e c o v e r a g e c a n d a m a g e m o s t u n h a r d e n e d e l e c t r o n i c a n d e l e c t r i c a l s y s t e m s a n d d r a w m a n y a t t e n t i o n s w o r l d w i d e T h e a n a l y s i s o f H E M P c o u p l i n g i n t o s y s t e m s i s t h e b a s i s f o r t h e H E M P e v a l u a t i o n a n d p r o t e c t i o n T h i s p a p e r i n t r o d u c e s t h e r e s e a r c h p r o g r e s s o f s y s t e m l e v e l H E M P c o u p l i n g a n a l y s i s m e t h o d s w h i c h h a s b e e n m a d e b y t h e a u t h o r s t e a m T h e p r o g r e s s i n c l u d e s t w o a s p e c t s o n e i s t h e t r a n s i e n t e l e c t r o m a g n e t i c t o p o l o g y m e t h o d w h i c h i s a p p l i e d t o a n a l y z e H E M P c o u p l i n g i n t o l o c a l m u l t i-s c a l e m u l t i v a r i a t e c o m p l e x s y s t e m s t h e o t h e r i s t h e m e t h o d b a s e d o n e f f e c t i v e c o u p l i n g l e n g t h w h i c h i s e m p l o y e d t o a n a l y z e t h e H E M P c o u p l i n g i n t o w i d e a r e a d i s t r i b u t e d s y s t e m s K e y w o r d s h i g h-a l t i t u d e e l e c t r o m a g n e t i c p u l s e e l e c t r o m a g n e t i c t o p o l o g y e f f e c t i v e c o u p l i n g l e n g t h c o u p l i n g a n a l y s i s l o c a l a r e a s y s t e m w i d e a r e a s y s t e m 高空电磁脉冲(h i g h-a l t i t u d e e l e c t r o m a g n e t i c p u l s e,HEMP)通常是指爆高大于3 0 k m的核爆炸在地球表面产生的电磁脉冲,按照产生机制和时间可分成早期(E 1)、中期(E 2)和晚期(E 3)3部分14。由于幅值强,频谱宽和覆盖范围广,HEMP可对大面积范围内的绝大多数电子电力系统造成损伤,引起国内外广泛关注。HEMP可通过电子系统天线的“前门耦合”或屏蔽腔体上孔缝和各类线缆的“后门耦合”进入电子系统内部57,在关键易损设备端口产生高电压或大电流,从而对内部电路和器件造成损伤。根据作用目标的特点,作用目标可分成两大类:一类是飞机等局域目标;另一类是电力系统等广域目标。飞机等局域目标包含腔体、线缆、电路及器件等1-201020第1 4卷 第2期2 0 2 3年6月现 代 应 用 物 理MO D E R N A P P L I E D P HY S I C SV o l.1 4,N o.2J u n.2 0 2 3结构,其HEMP耦合分析问题是一个多尺度多元复杂系统耦合分析问题。针对该问题,B a u m等8和T e s c h e等9提出了电磁拓扑的概念,并建立传输线HEMP耦合计算的B a u m-L i u-T e s c h e(B L T)方程。P a l e t t a等1 0和P a r m a n t i e r等1 1进一步发展了电磁拓扑方法,将全波数值模拟和B L T方程相结合,用于飞 机 等 局 域 目 标 的HEMP耦 合 分 析。由 于HEMP幅值高,与电子系统相互作用存在非线性,一方面是由于电子系统内部限幅器等非线性保护器件在高电压或大电流情况下进入工作状态,另一方面是由于产生的高电压或大电流对内部器件产生损伤等非线性效应。而B L T方程为频域方程,不能用于非线性分析。T e s c h e1 2建立了B L T方程的时域形式,然而该方法涉及大量卷积和迭代计算,同时不便于与器件电路模拟相结合。电力系统等广域目标包含输电网及配电网等网络结构,其HEMP耦合分析问题是一个广域分布式系统 耦 合 分 析 问 题。国 外 多 次 开 展 电 力 系 统HEMP耦合效应 评估,最有影 响力 的 为2 0 1 0年M e t a t e c h公司的评估1 3。然而这些评估基本是针对单根架空线开展耦合分析,虽指出耦合电压不随着架空线线性增加,但对架空线互连的影响及如何选择架空线长度未开展研究。近1 0年 来,本 文 作 者 所 在 团 队 围 绕 系 统 级HEMP耦合分析开展深入研究,分别针对局域目标和广域目标的耦合特点建立了瞬态电磁拓扑方法和基于有效耦合长度(e f f e c t i v e c o u p l i n g l e n g t h,E C L)的网络耦合分析方法,取得了多项进展。1瞬态电磁拓扑方法1.1主要思想瞬态电磁 拓 扑 方 法 的 整 体 思 想 与 电 磁 拓 扑方法的 整 体 思 想 类 似,都 是 基 于 良 好 屏 蔽 近 似(g o o d s h i e l d i n g a p p r o x i m a t i o n,G S A)91 1。在G S A下,只 考 虑 外 部 场 对 系 统 内 线 缆 的 耦 合 作用,而忽略系统内部线缆电流辐射对场的影响。一 般 而 言,对 电 子 系 统 内 的 大 部 分 线 缆,基 于G S A,当频率小于5 0 0 MH z时,可得到较好的结果1 1。在G S A下,外 部HEMP与 系 统 耦 合 问题,可分解 成HEMP与 系 统 屏 蔽 腔 等 耦 合 的 耦合场,及 耦 合 场 与 线 缆 的 耦 合 作 用。图1为HME P与电子系统耦合示意图。图1 H EM P与电子系统耦合示意图F i g.1 S c h e m a t i c d i a g r a m o f H EM P c o u p l i n g w i t h e l e c t r o n i c s y s t e m s时域有限差分法和有限积分方法等3维电磁场数值模拟方法可有效分析HEMP与系统屏蔽腔等耦合的耦合场。对于耦合场与线缆的耦合作用,原有电磁拓扑方法采用B L T方程求解场激励下的传输线网络响应。B L T方程可表示为9V=(I+S)(-S)-1V S(1)其中:V为电压向量;I为单位矩阵;S为传输线网络节点或终端的散射矩阵;为传播矩阵,V S为外部场对传输线网络的激励作用。S参数把终端看成黑匣子,只能用于线性负载,同时不易与电路器件模拟相结合。瞬态电磁拓扑方法提出采用外场激励下的传输线的S P I C E等效电路模型来分析系统内线缆的响应1 4。传输线的S P I C E等效电路模型把场线耦合问题转换成电路分析问题,具有以下优点:(1)可用于非线性负载;(2)可易与电路、器件模拟相结合;(3)模块性好,传输线网络的S P I C E等效电路可很方便地由传输线的S P I C E等效电路组合而成。1.2外场激励下传输线的S P I C E等效电路模型传输线的S P I C E等效电路模型基于传输线理论。由于电子系统内传输线的损耗较小,在HEMP频段可忽略传输线的损耗1 5。为简单起见,本文以无损单导线为例,其时域传输线方程可表示为zV(z,t)+LtI(z,t)=Vs(z,t)zI(z,t)+CtV(z,t)=Is(z,t)(2)其中:V和I为传输线上的电压和电流;L和C为传输线的电感和电容;Vs和Is为外部激励场的等效分布电压源和电流源;t为时间。式(2)为波动方程,通过系列变换可表示为2-201020 谢海燕:系统级HEMP耦合分析方法研究进展第2期V(0,t)-ZCI(0,t)=V(l,t-T)-ZCI(l,t-T)+V0(t)V(l,t)-ZCI(l,t)=V(l,t-T)+ZCI(l,t-T)+Vl(t)(3)其中:V0和Vl为外场在传输线两端的等效集总电压源;ZC为传输线的阻抗。式(3)可用电路实现,无损单导传输线S P I C E等效电路的一般形式如图2所示。图2无损单导传输线S P I C E等效电路的一般形式F i g.2 G e n e r a l f o r m o f S P I C E c i r c u i t f o r l o s s l e s s t r a n s m i s s i o n l i n e 在平面波激励下多导体传输线的S P I C E等效电路模型基础上,本文作者团队建立了非平面波激励下多导 传 输 线 的F D T D-S P I C E混 合 模 型1 6,其 中F D T D方法用于求解腔内传输线的非均匀激励场。地面上水平架设的屏蔽电缆可看成两个传输线系统:一是屏蔽层外部与地构成外部传输线系统;二是屏蔽层内部与内部芯线构成内部传输线系统。内外两个传输线系统通过转移阻抗和转移导纳相联系。基于内外传输线系统的传输线方程,文献1 7和文献1 8 分别建立平面波激励下的同轴和多芯屏蔽电缆的S P I C E等效电路模型。结合F D T D数值模拟,文献1 9 建立了非均匀场激励下的同轴和多芯屏蔽电缆的S P I C E等效电路模型。除多导线及屏蔽电缆等线型外,本文作者团队针对常见的双绞线进一步建立了平面波激励下双绞线的S P I C E等效电路模型2 0。上述各类模型中,传输线的布放假设都是均匀且确定的。然而实际系统中,传输线的布放具有一定的不确定性且沿线方向具有非均匀性。文献2 1 针对非均匀传输线建立了级联传输线的模型,即非均匀传输线看成多段均匀传输线的级联。文献2 2 和文献2 3 研究了传输线随机平移和随机转动对耦合信号的影响,并在此基础上建立了平面波激励 下随机平 移 和 随 机 转 动 传 输 线 的S P I C E等效电路模型。同时,基于电磁拓扑思想和外场激励下传输线的S P I C E等效电路模型,研发了瞬态电磁拓扑软件T E MT,并通过实验初步验证了其在系统级H E M P耦合分析方面的能力2 42 5。1.3 T EMT及典型应用T EMT软件采用C+语言编制,可支持多种传输线类型,如表1所列。T EMT可根据定义的线缆结构和激励场生成支持S P I C E 2和S P I C E 3格式的传输线等效电路,能用于P s p i c e,n g s p i c e等软件。为引入激励场,T EMT软件设有与C S T等3维电磁场数值模拟软件的接口。表1 T EM T支持的传输线类型T a b.1 T y p e s o f t r a n s m i s s i o n l i n e s u p p o r t e d b y T EMTL i n e t y p eN o e x c i t a t i o nf i e l dP l a n e w a v ef i e l dN o n-u n i f o r mf i e l dS i n g l e c o n d u c t o r/m u l t i p l e c o n d u c t o r sC o a x i a l/m u l t i-c o r es h i e l d e dP C B m i c r o s t r i p l i n eR a n d o m s i n g l e c o n d u c t o r/m u l t i p l e c o n d u c t o r s图3为外场激励下含有复杂线缆网络的小车。布线近端为三线传输线,分成两组传输线与不同终端相连。传输线在布放过程具有一定的随机性,图3中给出了三线传输线和双线传输线截面的可能变化。图4为车内负载上的平均耦合电压随时间的变化关系2 6。图3外场激励下含有复杂线缆网络的小车F i g.3 A c a r w i t h c o m p l i c a t e d c a b l e n e t w o r ke x c i t e d b y a p l a n e w a v e f i e l d 3-201020第1 4卷现 代 应 用 物 理(a)V o l t a g e o f n e a r l o a d(b)V o l t a g e o f f a r l o a d图4车内负载上的平均耦合电压随时间的变化关系2 6F i g.4 A v e r a g e c o u p l e d v o l t a g e s v s.t i m e a t t h e l o a d s i n s i d e t h e c a r2 62基于E C L的网络耦合分析方法2.1E C L的形成机制针对局域范围内的传输线网络,B L T方程和传输线的S P I C E等效电路模型都可很好地建模计算。然而输电网、配电网及高铁控制网等都是广域网络,结构复杂、未知量多,无论是采用B L T方程还是传输线的S P I C E等效电路模型都存在计算量大及建模困难等问题。针对HEMP与架空线网络耦合问题,本文作者团队提出了长传输线的E C L的概念,揭示了其形成机制,并针对架空线HEMP耦合给出了E C L的经验表征公式2 72 9。一般而言,传输线越长,HEMP的耦合信号越大;当传输线增加到一定长度时,HEMP的耦合信号不再增加,此时对应长度称为E C L。HEMP与传输线耦合存在E C L的内在原因为HEMP的脉冲特性。图5为HEMP与长线耦合模型及E C L形成机制示意图。假设在t0时刻,HEMP到达传输线的近端并耦合产生一个信号,那么在t0+t时刻耦合信号主要由两部分组成,一部分是t0时刻产生的信号并经过t时间的衰减,另一部分为t0+t时刻产生的耦合信号。假如t0+nt时刻,在传输线上新耦合产生信号与t0时刻产生的信号不再重叠,此时耦合信号不再随着传输线的长度增加而增加,对应的传输线长度为E C L。(a)L o n g t r a n s m i s s i o n l i n e m o d e l(b)F o r m a t i o n m e c h a n i s m o f e f f e c t i v e c o u p l i n g l e n g t h(E C L)图5 H EM P与长线耦合模型及E C L形成机制示意图F i g.5 S c h e m a t i c d i a g r a m o f l o n g t r a n s m i s s i o n l i n e a n d i t s m e c h a n i s m f o r E C L E C L与传输线的距地高度h及HEMP的入射仰角、方位角和极化角等参数有关。对1 0 m高的架空线,E C L主要分布在几十米到几百米,一般不超过2 k m2 7。2.2基于E C L的网络耦合分析方法由于HEMP与传输线耦合存在E C L,可得出其终端耦合信号主要取决于与该终端连接且长度等于E C L的那一段传输线与HEMP的耦合2 8。基于此规律,长传输线网络与HEMP的耦合可分解成多根独立传输线与HEMP的耦合,其中独立传输线的长度应该大于等于E C L。图6为传输线网络及其耦合电压与E C L理论模型及其耦合电压2 8。4-201020 谢海燕:系统级HEMP耦合分析方法研究进展第2期图6传输线网络及其耦合电压与E C L理论模型及其耦合电压F i g.6 T r a n s m i s s i o n l i n e n e t w o r k w i t h c o u p l i n g v o l t a g e s a n d i t s E C L m o d e l w i t h c o u p l i n g v o l t a g e s 该传输线网络中传输线T1、T2和T3的长度分别为5 0,6 0,4 0 m,距地高度都为8 m。该场景下其有效耦合长度不超过4 0 m。由图6计算结果可见,基本可忽略E C L长度以外的线缆段耦合对终端耦合电压的影响。值得注意的是,若与终端直接相连那段线缆与HEMP的耦合为弱耦合,则从别处反射回来的信号可能占主要地位,但此时终端信号一般不会 很 大。通 过E C L理 论,可 使 传 输 线 网 络HEMP耦合计算的维度大大降低。2.3架空线网络H EMP耦合分析及快速预测在E C L理论的基础上,HEMP与架空线网络的耦合可分解成长度等于E C L的独立架空线与HEMP的耦合。然而由于架空线网络不同位置处面临的HEMP环境的幅值、入射方向和极化方向不同3 0,开展架空线网络HEMP耦合分析需考虑不同位置处HEMP环境的差异;同时由于架空线网络的错综复杂性,尤其是配网架空线网络,耦合分析也需考 虑 架 空 线 走 向 的 随 机 性。虽 然 预 测 单 一HEMP环境下某一架空线的计算量不大,但预测大面积范围下架空线网络HEMP耦合电压分布的计算量很大。本文作者团队在研究架空线HEMP耦合规律的基础上,基于人工神经网络建立了架空线HEMP耦合的快速预测模型3 1。图7为水平极化HEMP不同入射角时架空线匹配负载耦合电压的累积概率分布。图7中快速预测模型计算所需时间为1.5 s,是数值模拟所需时间的1/3 8 0。图7快速预测模型得到的累积概率与数值模拟结果的对比F i g.7 C u m u l a t i v e p r o p a b i l i t y o b t a i n e d b y p r e d i c t i o n m o d e lc o m p a r e d w i t h n u m e r i c a l r e s u l t在综合考虑不同位置处HEMP幅值、极化方向和入射方向差异性和架空线走向的随机性的基础上,本文作者团队基于上面的快速预测技术在分钟量级内给出了某一场景下1 0 k V架空线网络匹配负载时HEMP耦合电压峰值分布,如图8所示。与HEMP辐射环境的“笑脸”分布图1 3不同,HEMP架空线传导环境分布图更像一个“双肺”图,这主要是因为HEMP辐射环境在南北方向主要是水平极化,而水平极化HEMP与长线缆的耦合信号比垂直极化HEMP的耦合信号小。5-201020第1 4卷现 代 应 用 物 理图8某一场景下,1 0 k V架空线网络匹配负载时H EM P耦合相对峰值电压F i g.8 D i s t r i b u t i o n o f r e l a t i v e c o u p l e d v o l t a g e p e a k s g e n e r a t e db y H EMP o n 1 0 k V o v e r h e a d t r a n s m i s s i o n l i n e n e t w o r k s3结论及展望近1 0年来,本文作者团队围绕系统级HEMP耦合分析开展深入研究,取得了系列成果:(1)针对HEMP与多尺度多元复杂系统耦合分析问题,建立了瞬 态 电 磁 拓 扑 方 法,包 括 提 出 采 用 传 输 线 的S P I C E等效电路模型计算系统内传输线网络的耦合响应,建立了屏蔽电缆、双绞线、随机和非均匀线等多种线缆外场激励下的S P I C E等效电路模型,完成实验验证并给出典型应用示例。与传统频域电磁拓扑方法相比,瞬态电磁拓扑方法不仅可用于非线性情况,且便于与电路器件效应分析相结合。(2)针对HEMP与广域分布式系统耦合分析问题,发现长传输线HEMP耦合存在E C L的机制,建立了架空线E C L的经验表征和基于E C L理论的传输线网络耦合分析方法,并基于人工神经网络建立了架空线HEMP耦合快速预测模型,实现了分钟量级内给出典型场景下架空线网络耦合信号的空间分布。2 0 2 0年,本文作者团队提出的“信息化条件下国家关键基础设施如何防范重大电磁威胁”入选2 0 2 0年科协十大工程难题3 2,后续将进一步深入研究国家关键基础设施HEMP耦合效应评估。参考文献 1 I E C E l e c t r o m a g n e t i c c o m p a t i b i l i t y EM C -P a r t 2 E n v i r o n m e n t-S e c t i o n 9 D e s c r i p t i o n o f HEMP e n v i r o n m e n t-r a d i a t e d d i s t u r b a n c e b a s i c EMC p u b l i c a t i o n f i r s t e d i t i o n I E C 6 1 0 0 0-2-9 S 1 9 9 6 2 王建国 牛胜利 张 殿辉 等 高空核爆 炸效应参数手册 M 北京 原子能出版社 2 0 1 0 WANG J i a n-g u o N I U S h e n g-l i Z HAN G D i a n-h u i e t a l P a r a m e t e r h a n d b o o k o f h i g h a l t i t u d e n u c l e a r d e t o n a t i o n e f f e c t M B e i j i n g A t o m i c E n e r g y P r e s s 2 0 1 0 3 L I Y L I U L WAN G J G e t a l N u m e r i c a l s i m u l a t i o n o f t h e i n t e r m e d i a t e-t i m e h i g h-a l t i t u d e e l e c t r o m a g n e t i c p u l s e J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 2 2 6 4 5 1 4 2 31 4 3 0 4 王建国 刘 利 牛胜利 等 高空核爆炸环境数值模拟 J 现代应 用 物 理 2 0 2 3 1 4 1 0 1 0 1 0 1 WAN G J i a n-g u o L I U L i N I U S h e n g-l i e t a l N u m e r i c a l s i m u l a t i o n s o f e n v i-r o n m e n t a l p a r a m e t e r s o f h i g h-a l t i t u d e n u c l e a r e x p l o s i o n J M o d e r n A p p l i e d P h y s i c s 2 0 2 3 1 4 1 0 1 0 1 0 1 5 王建国 刘国治 周金山 微波孔缝线性耦合函数研究 J 强激光与粒子束 2 0 0 3 1 5 1 1 1 0 9 3 1 0 9 9 WAN G J i a n-g u o L I U G u o-z h i Z HOU J i n-s h a n I n v e s t i g a t i o n s o n f u n c-t i o n f o r l i n e a r c o u p l i n g o f m i c r o w a v e s i n t o s l o t s J H i g h P o w e r L a s e r a n d P a r t i c l e B e a m s 2 0 0 3 1 5 1 1 1 0 9 31 0 9 9 6 C HE N J WAN G J G A T h r e e-d i m e n s i o n a l s e m i-i m p l i c i t F D T D s c h e m e f o r c a l c u l a t i o n o f s h i e l d i n g e f f e c t i v e n e s s o f e n c l o s u r e w i t h t h i n s l o t s J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 0 7 4 9 2 3 5 4 3 6 0 7 C HE N J WAN G J G T h r e e-d i m e n s i o n a l d i s p e r s i v e h y b r i d i m p l i c i t-e x p l i c i t f i n i t e-d i f f e r e n c e t i m e-d o m a i n m e t h o d f o r s i m u l a t i o n s o f g r a p h e n e J C o m p u t e r P h y s i c s C o mm u n i c a t i o n s 2 0 1 6 2 0 7 2 1 1 2 1 6 8 B AUM C E L I U T K T E S C HE F M O n t h e a n a l y s i s o f g e n e r a l m u l t i c o n d u c t o r t r a n s m i s s i o n-l i n e n e t w o r k s I n t e r a c t i o n N o t e s 3 5 0 R 1 9 7 8 9 T E S C HE F M I ANO Z M V KA R L S S ON T EMC a n a l y s i s m e t h o d s a n d c o m p u t a t i o n a l m o d e l s M N e w Y o r k W i l e y 1 9 9 7 1 0 P A L E T T A L P A RMA N T I E R J P I S S A C F e t a l S u s c e p t i b i l i t y a n a l y s i s o f w i r i n g i n a c o m p l e x s y s t e m c o m b i n i n g a 3-D s o l v e r a n d a t r a n s m i s s i o n-l i n e n e t w o r k s i m u l a t i o n J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 0 2 4 4 2 3 0 93 1 7 1 1 P A RMAN T I E R J P N u m e r i c a l c o u p l i n g m o d e l s f o r c o m p l e x s y s t e m s a n d r e s u l t s J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 0 4 4 6 3 3 5 9 3 6 7 1 2 T E S C HE F M D e v e l o p m e n t a n d u s e o f t h e B L T e q u a t i o n i n t h e t i m e d o m a i n a s a p p l i e d t o a c o a x i a l c a b l e J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 0 7 4 9 1 3 1 1 1 3 S AVA G E E G I L B E R T J R A D A S KY W T h e e a r l y-t i m e E 1 h i g h-a l t i t u d e e l e c t r o m a g n e t i c p u l s e HEMP a n d i t s i m p a c t o n t h e U S p o w e r g r i d M e t a-R-3 2 0 R M e t a t e c h C o r p o r a t i o n G o l e t a C A U S A 2 0 1 0 1 4 X I E H Y WAN G J G F AN R Y e t a l A p p l i c a t i o n o f a S P I C E m o d e l f o r m u l t i c o n d u c t o r t r a n s m i s s i o n l i n e s i n e l e c t r o m a g n e t i c t o p o l o g y C P I E R S P r o c e e d i n g s C a m b r i d g e 2 0 0 8 1 5 X I E H Y WAN G J G F AN R Y e t a l S t u d y o f l o s s e f f e c t 6-201020 谢海燕:系统级HEMP耦合分析方法研究进展第2期o f t r a n s m i s s i o n l i n e s a n d v a l i d i t y o f a S P I C E m o d e l i n e l e c t r o m a g n e t i c t o p o l o g y J P I E R 2 0 0 9 9 0 8 9 1 0 3 1 6 X I E H Y WAN G J G F AN R Y e t a l A h y b r i d F D T D-S P I C E m e t h o d f o r t r a n s m i s s i o n l i n e s e x c i t e d b y a n o n u n i f o r m i n c i d e n t w a v e J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 0 9 5 1 3 8 1 1 8 1 7 1 7 X I E H Y WAN G J G F AN R Y e t a l S P I C E m o d e l s t o a n a l y z e r a d i a t e d a n d c o n d u c t e d s u s c e p t i b i l i t i e s o f s h i e l d e d c o a x i a l c a b l e s J I E E E T r a n s E l e c t r o m a g n C o m p a t 2 0 1 0 5 2 1 2 1 5 2 2 2 1 8 X I E H Y W A N G J G F A N R Y e t a l S P I C E m o d e l s f o r r a d i a t e d a n d c o n d u c t e d s u s c e p t i b i l i t y a n a l y s e s o f m u l t i c o n d u c t o r s h i e l d e d c a b l e s J P I E R 2 0 1 0 1 0 3 2 4 1 2 5 7 1 9 X I E H Y WAN G J G F AN R Y e t a l S P I C E m o d e l s f o r p r e d i c t i o n o f d i s t u r b a n c e s