Chip Collaborative Protection Design Method Based on Piecewise Linear Model for Transmission Line Pulse Transient Interference Signal

FU Lu; YAN Zhaowen; LIU Yuzhu; SU Lixuan

doi:10.11999/JEIT220975

Volume 45 Issue 9

Sep. 2023

Turn off MathJax

Article Contents

Article Navigation > Journal of Electronics & Information Technology > 2023 > 45(9): 3263-3271

Qiu Changxing, Liang Changhong, Han Yi. GENERALIZED MEASUREMENT EQUATION FOR MULTI PORT LOSSLESS NETWORK[J]. Journal of Electronics & Information Technology, 1990, 12(4): 352-360.

Citation:

FU Lu, YAN Zhaowen, LIU Yuzhu, SU Lixuan. Chip Collaborative Protection Design Method Based on Piecewise Linear Model for Transmission Line Pulse Transient Interference Signal[J]. Journal of Electronics & Information Technology, 2023, 45(9): 3263-3271. doi: 10.11999/JEIT220975

Qiu Changxing, Liang Changhong, Han Yi. GENERALIZED MEASUREMENT EQUATION FOR MULTI PORT LOSSLESS NETWORK[J]. Journal of Electronics & Information Technology, 1990, 12(4): 352-360.

Citation:

PDF( 2718 KB)

Chip Collaborative Protection Design Method Based on Piecewise Linear Model for Transmission Line Pulse Transient Interference Signal

doi: 10.11999/JEIT220975

School of Electronics and Information Engineering, Beihang University, Beijing 100191, China

Funds: The Time Domain Testing System of Electromagnetic Compatibility Conductivity Sensitivity (61427803), The Developing a Behavioral Co-Analysis Methodology for Studying Electromagnetic Susceptibility of 3D IC (TSV) and PCB (61271044)

Received Date: 2022-07-21
Rev Recd Date: 2022-11-18

Available Online: 2022-11-30

Publish Date: 2023-09-27

Abstract

Abstract

With the development trend of miniaturization, high density and high speed of electronic equipment, integrated circuit, as the basic core unit of electronic equipment, is also developing in this direction, which brings more and more serious problems of electromagnetic compatibility. Among them, the problem of electrostatic discharge has attracted more and more attention of designers, producers and users. In this paper, the chip is tested by Transmission Line Pulse (TLP) method, and the volt ampere characteristic data of the device in response to electrostatic discharge interference are obtained. Based on the TLP test data, the piecewise linear modeling method is applied to build the model of the chip to deal with the electrostatic discharge interference. Based on the equivalent circuit of the diode and the volt ampere characteristic data in its data book, the Transient Voltage Suppression (TVS) diode model is constructed and verified by TLP test. At the same time, combined with the above two models, this paper carries out the research on the collaborative protection design method of chip electrostatic discharge interference, and obtains the collaborative protection design process and examples of the chip. This method realizes the collaborative protection design of the chip by simulation, which can save the design cost and time.
- Integrated circuit,
- Electromagnetic compatibility,
- Transmission Line Pulse (TLP),
- Electrostatic discharge,
- Protective design

FullText(HTML)

1. 引言

微波管广泛应用于卫星通信、雷达系统、电子对抗以及科学研究领域^[1-3]，其中永磁聚焦系统常作为外部磁场用于电子束聚焦^[4]。目前，在大量的科学及工程应用中，有限元方法针对复杂结构仍然是主流的数值分析工具^[5,6]。在永磁聚焦系统仿真设计领域，国外Maxwell^[7]是最流行的3维有限元仿真设计商业软件，国内电子科技大学开发的微波管模拟器套装^[8]里面的永磁聚焦模拟器^[4]则是具有国家自主知识产权的3维磁场仿真设计软件。它们都拥有永磁聚焦系统仿真能力，但是当求解大规模和多尺度问题时，其有限元矩阵求解通常会花费大量的时间和内存，有时甚至由于缺少有效的预处理导致矩阵无法求解。

非重叠区域分解法采用“分而治之”的思想，将原求解域划分成若干个互不重叠的子区域进行求解，其具有天然的并行性，非常适合仿真大型复杂结构^[9]。本文采用基于磁标势有限元的非重叠区域分解方法来进行永磁聚焦系统的仿真计算。归根结底，基于磁标势的永磁聚焦系统仿真属于求解泊松方程边值问题。在这类边值问题的基于有限元的区域分解方法中，目前的处理方式主要有两种，一种是拉格朗日乘子类型的区域分解法^[10,11]；一种是基于内罚方式的区域分解法^[12,13]。前者会有两种不同类型的推导过程：一是保留拉格朗日乘子的方式，这种方式会额外增加未知量个数，并且会生成一个对称不定的鞍点矩阵系统，不利于方程的求解；二是在推导过程中消去拉格朗日乘子，这种方式会生成一个对称正定的系数矩阵，但是这种方式有可能会因为巨大的计算代价而不能显式地计算系数矩阵^[14]。基于内罚的区域分解方法则不需要引入诸如拉格朗日乘子类型的辅助变量，只需要将传输条件引入基于内罚方式的有限元弱形式推导过程中，目前采用的主要是Robin传输条件，由Lions^[15]首次提出，但目前基于Robin传输条件的内罚区域分解法需要考虑法向偏导项的计算，并且最终形成的有限元矩阵是非对称的。

本文提出的区域分解方法同样是基于内罚方式的，但是引入了一种新型传输条件，其来源于接触热阻的定义。相比于现有的方法，该区域分解方法的主要优势包括：(1)不需要引入多余的未知量，使得有限元矩阵维数更少；(2)有限元矩阵集成过程更加简单，只需要考虑区域交界面上的物理量，而且不需要进行法向偏导项的计算，更重要的是最终产生的有限元矩阵满足对称正定性，矩阵性质更好，非常适合采用共轭梯度法进行求解。通过对多个永磁结构的仿真计算可以发现，相比于商业软件Maxwell，本文所提出的区域分解方法在保证求解精度的同时，可以更加高效地实现对微波管永磁聚焦系统的仿真。

2. 基于有限元的区域分解方法

2.1 基于内罚的区域分解有限元弱形式

永磁磁场的磁标势有限元分析的边值问题为

$\left. \begin{aligned} & { - \nabla \cdot \mu \nabla \varphi = \nabla \cdot \mu {{{H}}_{\rm{c}}},\;\;\;{\rm{in}}\;\varOmega } \\ & {\varphi \; = \;{\rm{0,}}\quad\quad\quad\quad\quad\quad\quad{\rm{on}}\;{\varGamma _{\rm{v}}}} \\ & { - \mu \nabla \varphi \cdot {{n}} = \mu {{{H}}_{\rm{c}}} \cdot {{n}}{\rm{,}}\;\;\;{\rm{on}}\;{\varGamma _{\rm{m}}}} \end{aligned} \right\}$

(1)

其中， $\mu$ 为磁导率， $\varphi$ 为标量磁位， ${{{H}}_{\rm{c}}}$ 为永磁材料的矫顽场强， ${{n}}$ 为永磁边界的外法向矢量， $\varOmega$ 为求解域， ${\varGamma _{\rm{v}}}$ 为真空的边界， ${\varGamma _{\rm{m}}}$ 为永磁材料的边界。

为了便于推导区域分解有限元弱形式，将求解域分成2个子区域，如图1所示，其中 ${\varOmega _{\rm{1}}}$ , ${\varOmega _{\rm{2}}}$ 代表2个子区域， ${\varGamma _{{\rm{v1}}}}$ ， ${\varGamma _{{\rm{v2}}}}$ 为2个子区域的真空边界， ${\varGamma _{{\rm{m1}}}}$ , ${\varGamma _{{\rm{m2}}}}$ 为2个子区域的永磁边界， $\varGamma$ 为2个子区域的交界面， ${{{n}}_1}$ , ${{{n}}_2}$ 为交界面上的外法向矢量。

图 1 单个区域分成2个子区域示意图

下载: 全尺寸图片幻灯片

扩展单个区域的边值问题式(1)到2个子区域，可以得到

$- \nabla \cdot {\mu _{\rm{1}}}\nabla {\varphi _{\rm{1}}} = \nabla \cdot {\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}},\;\;\;{\rm{in}}\;{\varOmega _{\rm{1}}}$

(2)

$- \nabla \cdot {\mu _{\rm{2}}}\nabla {\varphi _{\rm{2}}} = \nabla \cdot {\mu _{\rm{2}}}{{{H}}_{{\rm{c}}{\rm{2}}}},\;\;\;{\rm{in}}\;{\varOmega _{\rm{2}}}$

(3)

$- {\mu _{\rm{1}}}\frac{{\partial {\varphi _{\rm{1}}}}}{{\partial {n_{\rm{1}}}}} = {\mu _{\rm{2}}}\frac{{\partial {\varphi _{\rm{2}}}}}{{\partial {n_{\rm{2}}}}},\;\;\;{\rm{on}}\;\varGamma$

(4)

$\;{\varphi _{\rm{1}}} = {\varphi _{\rm{2}}},\;\;{\rm{on}}\;\varGamma$

(5)

$\;{\varphi _{\rm{1}}}\; = \;{\rm{0}},\;\;{\rm{on}}\;{\varGamma _{{\rm{v1}}}}$

(6)

$\;{\varphi _{\rm{2}}}\; = \;{\rm{0}},\;\;{\rm{on}}\;{\varGamma _{{\rm{v2}}}}$

(7)

$- {\mu _{\rm{1}}}\nabla {\varphi _{\rm{1}}} \cdot {{{n}}_{\rm{1}}} = {\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}} \cdot {{{n}}_{\rm{1}}},\;\;\;{\rm{on}}\;{\varGamma _{\rm{m1}}}$

(8)

$- {\mu _{\rm{2}}}\nabla {\varphi _{\rm{2}}} \cdot {{{n}}_{\rm{2}}} = {\mu _{\rm{2}}}{{{H}}_{{\rm{c}}{\rm{2}}}} \cdot {{{n}}_{\rm{2}}},\;\;\;{\rm{on}}\;{\varGamma _{\rm{m2}}}$

(9)

其中，式(4)和式(5)用来保证区域交界面上物理量的连续性，但是其收敛性很差，常用的Robin传输条件也是通过两式的线性变换得到的。本文抛弃了之前的传输条件构造方式，而是从接触热阻的定义^[16]出发，开创性地提出了一种新型传输条件，其具体表达式为

$- {\mu _{\rm{1}}}\frac{{\partial {\varphi _{\rm{1}}}}}{{\partial {n_{\rm{1}}}}} = \gamma \left( {{\varphi _{\rm{1}}} - {\varphi _{\rm{2}}}} \right),\;\;{\rm{on}}\;\varGamma$

(10)

$- {\mu _{\rm{2}}}\frac{{\partial {\varphi _{\rm{2}}}}}{{\partial {n_{\rm{2}}}}} = \gamma \left( {{\varphi _{\rm{2}}} - {\varphi _{\rm{1}}}} \right),\;\;{\rm{on}}\;\varGamma$

(11)

其中， $\gamma$ 表示区域交界面上物理量连续程度的物理量，理论上当 $\gamma$ 无穷大时，区域交界面上的物理量就会完全连续^[17]，在实际应用过程中，只需要取一个比较大的值，10⁶量级基本可以满足需求。

为了推导2个区域的磁标势有限元弱形式，用式(10)和式(11)代替式(4)和式(5)，可以得到残差表达式

$R_{\rm{1}}^{{\varOmega _{\rm{1}}}} = \nabla \cdot {\mu _{\rm{1}}}\nabla {\varphi _{\rm{1}}}{\rm{ + }}\nabla \cdot {\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}},\;\;\;{\rm{in}}\;{\varOmega _{\rm{1}}}$

(12)

$R_{\rm{2}}^{{\varOmega _{\rm{2}}}} = \nabla \cdot {\mu _{\rm{2}}}\nabla {\varphi _2}{\rm{ + }}\nabla \cdot {\mu _{\rm{2}}}{{{H}}_{{\rm{c}}{\rm{2}}}},\;\;\;{\rm{in}}\;{\varOmega _{\rm{2}}}$

(13)

$R_{\rm{3}}^\varGamma = {\mu _{\rm{1}}}\frac{{\partial {\varphi _{\rm{1}}}}}{{\partial {n_{\rm{1}}}}} + \gamma \left( {{\varphi _{\rm{1}}} - {\varphi _{\rm{2}}}} \right),\;\;\;{\rm{on}}\;\varGamma$

(14)

$R_{\rm{4}}^\varGamma = {\mu _{\rm{2}}}\frac{{\partial {\varphi _{\rm{2}}}}}{{\partial {n_{\rm{2}}}}} + \gamma \left( {{\varphi _{\rm{2}}} - {\varphi _{\rm{1}}}} \right),\;\;\;{\rm{on}}\;\varGamma$

(15)

$R_{\rm{5}}^{{\varGamma _{{\rm{v1}}}}} = {\varphi _{\rm{1}}},\;\;{\rm{on}}\;{\varGamma _{{\rm{v1}}}}$

(16)

$R_{\rm{6}}^{{\varGamma _{{\rm{v2}}}}} = {\varphi _{\rm{2}}},\;\;{\rm{on}}\;{\varGamma _{{\rm{v2}}}}$

(17)

$R_{\rm{7}}^{{\varGamma _{{\rm{m1}}}}} = {\mu _{\rm{1}}}\nabla {\varphi _{\rm{1}}} \cdot {{{n}}_{\rm{1}}} + {\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}} \cdot {{{n}}_{\rm{1}}},\;\;\;{\rm{on}}\;{\varGamma _{\rm{m1}}}\!\!\!$

(18)

$R_{\rm{8}}^{{\varGamma _{{\rm{m2}}}}} = {\mu _{\rm{2}}}\nabla {\varphi _{\rm{2}}} \cdot {{{n}}_{\rm{2}}} + {\mu _{\rm{2}}}{{{H}}_{{\rm{c}}{\rm{2}}}} \cdot {{{n}}_{\rm{2}}},\;\;\;{\rm{on}}\;{\varGamma _{\rm{m2}}}\!\!\!\!\!\!$

(19)

首先定义体积分和面积分如式(20)

$\left. \begin{aligned} & {{{\left( {u,v} \right)}_\varOmega } = \int_\varOmega {\left( {u \cdot v} \right)} {\rm{d}}v} \\ & {{{\left\langle {u,v} \right\rangle }_\varGamma } = \int_\varGamma {\left( {u \cdot v} \right)} {\rm{d}}s} \end{aligned} \right\}$

(20)

其中， $u$ 是权函数， $v$ 为残差项。由式(12)—式(19)可以得到如式(21)所示的线性组合方程

$\begin{aligned} {\left( {{w_{\rm{1}}},R_{\rm{1}}^{{\varOmega _{\rm{1}}}}} \right)_{{\varOmega _{\rm{1}}}}}\,& {\rm{ + }}{\left( {{w_{\rm{2}}},R_{\rm{2}}^{{\varOmega _{\rm{2}}}}} \right)_{{\varOmega _{\rm{2}}}}}{\rm{ + }}{c_{\rm{1}}}{\left\langle {{w_{\rm{1}}},R_{\rm{3}}^\varGamma } \right\rangle _\varGamma }\\ & {\rm{ + }}{c_{\rm{2}}}{\left\langle {{w_{\rm{2}}},R_{\rm{4}}^\varGamma } \right\rangle _\varGamma } + {c_{\rm{3}}}{\left\langle {{w_{\rm{1}}},R_{\rm{5}}^{{\varGamma _{{\rm{v1}}}}}} \right\rangle _{{\varGamma _{{\rm{v1}}}}}} \\ & {\rm{ + }}{c_{\rm{4}}}{\left\langle {{w_{\rm{2}}},R_{\rm{6}}^{{\varGamma _{{\rm{v2}}}}}} \right\rangle _{{\varGamma _{{\rm{v2}}}}}} + {c_{\rm{5}}}{\left\langle {{w_{\rm{1}}},R_{\rm{7}}^{{\varGamma _{{\rm{m1}}}}}} \right\rangle _{{\varGamma _{{\rm{m1}}}}}} \\ & + {c_{\rm{6}}}{\left\langle {{w_{\rm{2}}},R_{\rm{8}}^{{\varGamma _{{\rm{m2}}}}}} \right\rangle _{{\varGamma _{{\rm{m2}}}}}} = {\rm{0}} \\[-15pt] \end{aligned}$

(21)

其中， ${c_1}$ , ${c_2}$ , ${c_3}$ , ${c_4}$ , ${c_5}$ 和 ${c_6}$ 为待定系数。

由格林公式，式(21)中的前两项可以写为

$\begin{split} {\left( {{w_i},R_i^{{\varOmega _i}}} \right)_{{\varOmega _i}}} =\,& - {\left( {\nabla {w_i},{\mu _i}\nabla {\varphi _i}} \right)_{{\varOmega _i}}} \\ & + {\left\langle {{w_i},{\mu _i}\frac{{\partial {\varphi _i}}}{{\partial {n_i}}}} \right\rangle _{\varGamma + {\varGamma _{{\rm{v}}i}} + {\varGamma _{{\rm{m}}i}}}} \\ \,& - {\left( {\nabla {w_i},{\mu _i}{{{H}}_{{\rm{ci}}}}} \right)_{{\varOmega _i}}} \\ & + {\left\langle {{w_i},{\mu _i}{{{H}}_{{\rm{ci}}}} \cdot {{{n}}_i}} \right\rangle _{{\varGamma _{{\rm{m}}i}}}},\;i = {\rm{1}},{\rm{2}} \end{split}$

(22)

对于式(22)中第1类边界条件项，由于基函数 ${w_i}$ 具有任意性，令 ${w_i} = 0$ ，有

${\left\langle {{w_i},{\mu _i}\frac{{\partial {\varphi _i}}}{{\partial {n_i}}}} \right\rangle _{{\varGamma _{{\rm{v}}i}}}} = {\rm{0}}$

(23)

也就是说，第1类边界条件项在有限元弱形式推导过程中可以不考虑，但需要采用强强加的方式添加到最终的有限元矩阵方程里。接下来，对于区域交界面上的项有

$\left. \begin{aligned} & {{\left\langle {{w_{\rm{1}}},R_{\rm{3}}^\varGamma } \right\rangle }_\varGamma } = {{\left\langle {{w_{\rm{1}}},{\mu _{\rm{1}}}\frac{{\partial {\varphi _{\rm{1}}}}}{{\partial {n_{\rm{1}}}}}} \right\rangle }_\varGamma } \\ &\qquad\qquad\quad + {{\left\langle {{w_{\rm{1}}},\gamma \left( {{\varphi _{\rm{1}}} - {\varphi _{\rm{2}}}} \right)} \right\rangle }_\varGamma } \\ & {{\left\langle {{w_{\rm{2}}},R_{\rm{4}}^\varGamma } \right\rangle }_\varGamma } = {{\left\langle {{w_{\rm{2}}},{\mu _{\rm{2}}}\frac{{\partial {\varphi _{\rm{2}}}}}{{\partial {n_{\rm{2}}}}}} \right\rangle }_\varGamma } \\ &\qquad\qquad\quad+ {{\left\langle {{w_{\rm{2}}},\gamma \left( {{\varphi _{\rm{2}}} - {\varphi _{\rm{1}}}} \right)} \right\rangle }_\varGamma } \end{aligned} \right\}$

(24)

对于永磁的边界项有

$\left. \begin{aligned} & {{\left\langle {{w_{\rm{1}}},R_{\rm{7}}^{{\varGamma _{{\rm{m1}}}}}} \right\rangle }_{{\varGamma _{{\rm{m1}}}}}} = {{\left\langle {{w_{\rm{1}}},{\mu _{\rm{1}}}\frac{{\partial {\varphi _{\rm{1}}}}}{{\partial {n_{\rm{1}}}}}} \right\rangle }_{{\varGamma _{{\rm{m1}}}}}} \\ &\qquad\qquad\qquad\quad + {{\left\langle {{w_{\rm{1}}},{\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}} \cdot {{{n}}_{\rm{1}}}} \right\rangle }_{{\varGamma _{{\rm{m1}}}}}} \\ & {{\left\langle {{w_{\rm{2}}},R_{\rm{8}}^{{\varGamma _{{\rm{m2}}}}}} \right\rangle }_{{\varGamma _{{\rm{m2}}}}}} = {{\left\langle {{w_{\rm{2}}},{\mu _{\rm{2}}}\frac{{\partial {\varphi _{\rm{2}}}}}{{\partial {n_{\rm{2}}}}}} \right\rangle }_{{\varGamma _{{\rm{m2}}}}}}\\ &\qquad\qquad\qquad\quad + {{\left\langle {{w_{\rm{2}}},{\mu _{\rm{2}}}{{{H}}_{{\rm{c2}}}} \cdot {{{n}}_{\rm{2}}}} \right\rangle }_{{\varGamma _{{\rm{m2}}}}}} \end{aligned} \right\}$

(25)

令 ${c_1} = - 1$ , ${c_2} = - 1$ , ${c_5} = - 1$ , ${c_6} = - 1$ ，可以得到如式(26)所示的有限元弱形式表达式

$\begin{split} & {\left( {\nabla {w_{\rm{1}}},{\mu _{\rm{1}}}\nabla {\varphi _{\rm{1}}}} \right)_{{\varOmega _{\rm{1}}}}} + {\left( {\nabla {w_{\rm{2}}},{\mu _{\rm{2}}}\nabla {\varphi _{\rm{2}}}} \right)_{{\varOmega _{\rm{2}}}}} + {\left\langle {{w_{\rm{1}}},\gamma {\varphi _{\rm{1}}}} \right\rangle _\varGamma } \\ & \qquad + {\left\langle {{w_{\rm{2}}},\gamma {\varphi _{\rm{2}}}} \right\rangle _\varGamma } - {\left\langle {{w_{\rm{1}}},\gamma {\varphi _{\rm{2}}}} \right\rangle _\varGamma } - {\left\langle {{w_{\rm{2}}},\gamma {\varphi _{\rm{1}}}} \right\rangle _\varGamma } \\ &\quad = - {\left( {\nabla {w_{\rm{1}}},{\mu _{\rm{1}}}{{{H}}_{{\rm{c}}{\rm{1}}}}} \right)_{{\varOmega _{\rm{1}}}}} - {\left( {\nabla {w_{\rm{2}}},{\mu _{\rm{2}}}{{{H}}_{{\rm{c}}{\rm{2}}}}} \right)_{{\varOmega _{\rm{2}}}}} \\[-10pt] \end{split}$

(26)

不难发现式(26)可以扩展到任意多个子区域的情形。

2.2 有限元离散

由于四面体单元在处理复杂边界时具有良好的适应性，同时为了使用较少的网格和自由度得到较高的计算精度，本文采用了基于四面体单元的2阶叠层标量基函数进行有限元离散。首先定义体积坐标，四面体内的体积坐标满足式(27)

${L_1} + {L_2} + {L_3} + {L_4} = 1$

(27)

2阶叠层标量基函数则从体积坐标出发，由1阶基函数构造出2阶基函数，包括了4个顶点基函数和6个边基函数^[18]，构造形式为

$\begin{split} W_{10}^2 =\,& \left\{ {{L_1},} \right.{L_2},{L_3},{L_4},{L_1}{L_2},{L_1}{L_3},{L_1}{L_4},\\ &{L_2}{L_3}, {L_2}{L_4},{L_3}{L_4}\} \end{split}$

(28)

用2阶基函数去离散式(26)，便可以得到式(29)所示的矩阵方程

$\left[ {\begin{array}{*{20}{c}} {{{{A}}_1}}&{{{{C}}_{12}}}\\ {{{{C}}_{21}}}&{{{{A}}_2}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{{{x}}_1}}\\ {{{{x}}_2}} \end{array}} \right]{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{{{y}}_1}}\\ {{{{y}}_2}} \end{array}} \right]$

(29)

由于式(29)的系数矩阵是对称正定的，本文采用了包括块雅可比和多波前块不完全楚列斯基分解^[19]的两层预处理的共轭梯度算法来进行矩阵方程的求解，相比于传统有限元法，可以大幅提高求解效率和减少内存消耗。经过有限元分析得到标量磁位 $\varphi$ 的值，就可以根据式(30)得到磁感应强度 ${{B}}$ 的值

${{B}} = \mu \left( { - \nabla \varphi + {{{H}}_{\rm{c}}}} \right)$

(30)

2.3 区域划分

本文使用METIS软件包^[20]进行区域的划分，经过大量的对比分析，区域划分过程需要考虑以下几点：

(1) 区域划分的个数对并行效率的影响很大，随着子区域数目的增加，并行计算效率会逐步增加，虽然理论上子区域数目可以随意取值，但是实际上随着区域数目的进一步增加，线程之间的资源竞争会更加激烈，并且线程切换花销也随之增大，会使得并行效率降低。

(2) 考虑到求解过程中每个子区域矩阵都需要进行预处理，为了避免线程等待，划分区域时应尽量使得每个子区域大小相当。

(3) 划分区域时应尽量使得区域交界面数量少，可以加快矩阵求解过程收敛速度，从而提高计算效率。

3. 仿真实例

本节通过仿真多个微波管永磁聚焦系统，并与商业软件 Maxwell对比，来验证所提出的基于有限元的区域分解方法的准确性和高效性。区域分解法中的因子 $\gamma$ 取值为10⁸，多波前块不完全楚列斯基分解残差为10^–4，预处理共轭梯度法收敛残差和Maxwell一样为10^–6。所有的仿真计算都是在一台小型工作站(Windows 10, Intel Xeon 5122 3.60 GHz 3.59 GHz 双处理器，16 threads, 128 GB RAM)上完成的。

3.1 单周期结构仿真

如图2，选取了一个典型的单周期结构^[4]，采用的区域划分方式为沿着Z轴方向并尽量使得每个子区域的大小相当。首先与商业软件Maxwell进行精度上的比较，图3绘制了本文提出的区域分解法和Maxwell软件轴切面上的磁感应强度云图分布，可以看到其磁感应强度云图分布趋势一致。由于磁钢、极靴与真空交界处磁场变化比较剧烈，此处的磁感应强度由于网格因素会产生奇异值，所以将显示范围固定为0～1 T。

图 2 单周期结构计算模型及区域分解示意图

下载: 全尺寸图片幻灯片

图 3 区域分解法与Maxwell软件轴切面磁感应强度云图对比

下载: 全尺寸图片幻灯片

如图4所示，将轴线上的磁场与Maxwell进行对比，可以看到其吻合情况很好，在两个峰值点处的相对误差分别为0.12%和0.06%。此外还与Maxwell进行了计算性能对比，如表1所示，随着子区域个数的增加，计算时间和内存相比于Maxwell的优势越来越大，当划分为20个子区域时，时间加速比达到了11.4倍，而内存只有Maxwell的53.5%. 这里需要注意的是：在线程数更多的计算机上，区域数的增加会带来更加优越的计算性能，这里划分到20个区域已经足以说明所提出的区域分解方法具有非常好的计算优势。

图 4 单周期结构轴线磁场

${B_z}$ 分布

下载: 全尺寸图片幻灯片

表 1 单周期结构区域分解法与Maxwell软件性能对比

求解方法	子区域数	网格数	计算时间(s)	峰值内存(MB)
Maxwell		2428919	435	10342
区域分解法	8	2718101	291	8129
	12	2718101	135	7111
	16	2718101	132	6461
	20	2718101	38	5532

下载: 导出CSV

| 显示表格

3.2 Wiggler结构仿真

本实例考虑两周期Wiggler结构^[4]的仿真计算，磁钢材料为SmCo28，为了展示区域分解法的计算效率并兼顾区域划分的方便快捷，将计算模型固定划分为20个子区域，其计算模型和区域分解示意图如图5所示。首先进行磁感应强度云图的对比，如图6所示，可以看到两者轴切面上的云图分布趋势相同。为了进一步证明所提出的区域分解方法的准确性，选取了轴线上的磁感应强度与Maxwell进行对比，由图7可以看到其吻合程度非常好，并且其峰值处的相对误差最大不超过0.14%。

图 5 Wiggler计算模型和区域划分示意图

下载: 全尺寸图片幻灯片

图 6 区域分解法与Maxwell轴切面磁感应强度云图分布对比

下载: 全尺寸图片幻灯片

图 7 Wiggler结构轴线B_y分布和峰值相对误差曲线

下载: 全尺寸图片幻灯片

此外，如表2所示，与Maxwell对比了3组网格数目相当情况下的计算时间和峰值内存，可以看到3组不同实例下的时间加速比分别为3.7, 3.2和4.2，但是其峰值内存分别只有Maxwell的77%, 82%和73%，充分证明了所提出的区域分解法的高效性。

表 2 Wiggler结构区域分解法与Maxwell性能对比

实例	网格数	求解方法	计算时间(s)	峰值内存(MB)
实例1	5987880	Maxwell	2171	30515
	6476933	区域分解法	589	23407
实例2	7784252	Maxwell	2717	35942
	8200780	区域分解法	858	29387
实例3	9014971	Maxwell	4766	45875
	9158627	区域分解法	1129	33466

下载: 导出CSV

| 显示表格

4. 结束语

本文针对微波管中的永磁聚焦系统仿真，提出了一种先进的基于有限元的区域分解求解技术，并对其理论进行了详细的描述，还给出了实际应用中区域划分的相关原则和技巧。通过对多个永磁结构的建模与仿真计算，并与商业软件Maxwell进行详细的对比，验证了本文提出的区域分解方法的准确性和高效性。本文给出的针对永磁聚焦系统仿真的区域分解求解技术后续有望集成到微波管模拟器套装^[8]中，为管型设计师提供更好的仿真设计平台。

References(28)

References

[1]	付蕊. 基于场路协同仿真的电子设备电磁敏感性研究[D]. [硕士论文], 西安电子科技大学, 2020. FU Rui. Research on electromagnetic susceptibility of electronic equipment based on field-circuit cooperated simulation[D]. [Master dissertation], Xidian University, 2020.
[2]	JEON C, LIM H, JIN J, et al. A new method of component-level ESD test to assess system-level ESD risk for ICs[C]. 2019 Joint International Symposium on Electromagnetic Compatibility, Sapporo and Asia-Pacific International Symposium on Electromagnetic Compatibility (EMC Sapporo/APEMC), Sapporo, Japan, 2019: 681–684.
[3]	BOUANGEUNE D, HONG W K, CHOI S S, et al. ESD robustness of low-voltage/high-speed TVS devices with epitaxial grown films[C]. The 1st IEEE Global Conference on Consumer Electronics 2012, Tokyo, Japan, 2012: 189–192.
[4]	KHURANA N, MALONEY T, and YEH W. ESD on CHMOS devices - equivalent circuits, physical models and failure mechanisms[C]. The 23rd International Reliability Physics Symposium, Orlando, USA, 1985: 212–223.
[5]	MAAS J S, PRATT D J, and BOXLEITNER W. Furniture ESD-the forgotten parameter in ESD testing[C]. IEEE 1991 International Symposium on Electromagnetic Compatibility, Cherry Hill, USA, 1991: 248–252.
[6]	TS EN 60749–28-2008. Semiconductor devices - Mechanical and climatic test methods - Part 28: ElectroStatic Discharge (ESD) sensitivity testing - Direct contact charged device model (DC - CDM)[S]. 2018.
[7]	ESD. ANSI/ESD STM5.5. 1-2016 ESD association standard test method for elcetrostatic discharge sensitivity testing-Transmission Line Pulse (TLP) - device level[S]. Rome: EOS/ESD Association, 2016.
[8]	WANG A Z H. On-Chip ESD Protection for Integrated Circuits[M]. Boston: Kluwer Academic Publishers, 2002.
[9]	MUHONEN K, GRUND E, and ASHTON R. High-Speed TLP and ESD Characterization of ICs[C]. 2021 IEEE BiCMOS and Compound Semiconductor Integrated Circuits and Technology Symposium (BCICTS), Monterey, USA, 2021: 1–6.
[10]	BSI. IEC 62615: 2010 Electrostatic discharge sensitivity testing - Transmission Line Pulse (TLP)-Component level[S]. London: BSI, 2010.
[11]	DI BICCARI L, CERATI L, ZULLINO L, et al. HV ESD diodes investigation under vf-TLP stresses: TCAD approach[C]. 2016 38th Electrical Overstress/Electrostatic Discharge Symposium (EOS/ESD), Garden Grove, USA, 2016: 1–10.
[12]	CHEN Qi, LI Cheng, LU Fei, et al. TLP measurement and analysis of graphene NEMS switches for on-chip ESD protection[C]. 2017 IEEE 12th International Conference on Nano/Micro Engineered and Molecular Systems (NEMS), Los Angeles, USA, 2017: 370–374.
[13]	HUANG Shaochang, LI C H, CHEN Lifan, et al. The correlations between ESD and TLP in large array devices[C]. 2020 IEEE International Conference on Consumer Electronics - Taiwan (ICCE-Taiwan), Taoyuan, China, 2020: 1–2.
[14]	MA Rui, WANG Li, ZHANG Chen, et al. TLP and HBM ESD test correlation for power ICs[C]. 2013 IEEE International Conference of Electron Devices and Solid-state Circuits, Hong Kong, China, 2013: 1–2.
[15]	XI Yunfeng, MALOBABIC S, VASHCHENKO V, et al. Correlation between TLP, HMM, and system-level ESD pulses for Cu metallization[J]. IEEE Transactions on Device and Materials Reliability, 2014, 14(1): 446–450. doi: 10.1109/TDMR.2013.2292039
[16]	SANTORO R P. Piecewise-linear modeling of I-V characteristics with SPICE[J]. IEEE Transactions on Education, 1995, 38(2): 107–117. doi: 10.1109/13.387211
[17]	AUSTERMANN C, SCHEIER S, FREI S, et al. Modellierung von ESD schutzelementen mit snapback-verhalten[J]. Tagungsband 14. ESD-Forum, 2015: 115–120.
[18]	王志明. 应用于CAN总线静电保护的低电容TVS二极管设计[D]. [硕士论文], 电子科技大学, 2021. WANG Zhiming. Design of low capacitance TVS diodes for CAN bus electrostatic protection[D]. [Master dissertation], University of Electronic Science and Technology of China, 2021.
[19]	KASHYAP A S, SANDVIK P, MCMAHON J, et al. Silicon carbide transient voltage suppressor for next generation lightning protection[C]. 2014 IEEE Workshop on Wide Bandgap Power Devices and Applications, Knoxville, USA, 2014: 147–150.
[20]	黄薇, 罗启元, 蒋同全, 等. 瞬态电压抑制二极管的建模及仿真[J]. 太赫兹科学与电子信息学报, 2014, 12(6): 927–931,936. doi: 10.11805/TKYDA201406.0927 HUANG Wei, LUO Qiyuan, JIANG Tongquan, et al. Modeling and simulation of the transient voltage suppression diode[J]. Journal of Terahertz Science and Electronic Information Technology, 2014, 12(6): 927–931,936. doi: 10.11805/TKYDA201406.0927
[21]	周月宾, 饶宏, 许树楷, 等. 一种二极管箝位型MMC的高效等值建模方法[J]. 中国电机工程学报, 2016, 36(7): 1925–1932. doi: 10.13334/j.0258-8013.pcsee.2016.07.020 ZHOU Yuebin, RAO Hong, XU Shukai, et al. An equivalent efficient modeling approach for diode clamp sub-module based MMC[J]. Proceedings of the CSEE, 2016, 36(7): 1925–1932. doi: 10.13334/j.0258-8013.pcsee.2016.07.020
[22]	SHEN Yushu, KER M D, and JIANG H C. Transient Voltage Suppressor (TVS) on signal integrity of microelectronics system with CMOS ICs under system-level ESD test[C]. 2020 IEEE International Symposium on the Physical and Failure Analysis of Integrated Circuits (IPFA), Singapore, 2020: 1–4.
[23]	SHEN Yushu and KER M D. The impact of holding voltage of Transient Voltage Suppressor (TVS) on signal integrity of microelectronics system with CMOS ICs under system-level ESD and EFT/burst tests[J]. IEEE Transactions on Electron Devices, 2021, 68(5): 2152–2159. doi: 10.1109/TED.2021.3063208
[24]	MAGHLAKELIDZE G, MARATHE S, HUANG Wei, et al. Effect of RF signals on TVS diode trigger voltage for ESD protection[C]. 2020 IEEE International Symposium on Electromagnetic Compatibility & Signal/Power Integrity (EMCSI), Reno, USA, 2020: 194–199.
[25]	ZENG Xun, ZHOU Liang, and ZHANG Chengrui. Investigation of protection effects using transient voltage suppressor diodes-based circuits under high power microwave pulses[C]. 2021 IEEE International Joint EMC/SI/PI and EMC Europe Symposium, Raleigh, USA, 2021: 524–527.
[26]	张杨. 基于ESD的TVS冲击电压特性机理研究[J]. 电子质量, 2022(6): 188–193. doi: 10.3969/j.issn.1003-0107.2022.06.046 ZHANG Yang. Research on impulse voltage characteristic of TVS based on ESD[J]. Electronics Quality, 2022(6): 188–193. doi: 10.3969/j.issn.1003-0107.2022.06.046
[27]	陈强, 徐可, 陈真真, 等. 基于TLP建模的系统级静电放电效应仿真[J]. 强激光与粒子束, 2019, 31(10): 103208. doi: 10.11884/HPLPB201931.190113 CHEN Qiang, XU Ke, CHEN Zhenzhen, et al. System-level electrostatic discharge simulation based on transmission line pulse modeling[J]. High Power Laser and Particle Beams, 2019, 31(10): 103208. doi: 10.11884/HPLPB201931.190113
[28]	李景钦. 二极管电磁兼容模型建立方法研究[D]. [硕士论文], 电子科技大学, 2021. LI Jingqin. Research on establishing method of diode electromagnetic compatibility model[D]. [Master dissertation], University of Electronic Science and Technology of China, 2021.