一种吸纳不同步长接口延时的并行计算分网方法
【技术领域】
[0001] 本发明涉及一种电磁暂态仿真的仿真方法，具体涉及一种吸纳不同步长接口延时 的并行计算分网方法。
【背景技术】
[0002] 不同步长仿真技术算法的并行化程度不高，系统交互信息量较大，实现实时仿真 有难度。而并行接口算法能把两侧网络解耦，提升不同步长仿真的并行度，降低系统交互信 息，是实现实时仿真的可能解决方案。
[0003] 经过调研，目前可采用的并行不同步长仿真接口包括如下方法：
[0004] ITM法是一种直接的方法，它将电路按照联络线分开，一侧用电压源表示，另一侧 用电流源表示。根据大仿真步长侧网络使用的电压源或电流源，ITM方法又可分为电压主导 型和电流主导型。ITM最明显的弱点是有一个仿真时步的延迟，接口误差较大，同时，子网电 路参数对接口稳定性影响较大。
[0005] TFA法是假设两侧电路都可以表示为简单一阶线性电路，通过历史值实时计算一 阶电路的等值参数。本质上它是预测一阶电路的状态变量的值，与其他预测法类似，TFA方 法在解决非线性和高频信号时有局限性。另外，TFA方法存在不稳定的缺点。
[0006] PCD法是一种非常复杂的网络解耦方法，它依托了松弛迭代技术，把原始网络模拟 成一个多端口的网络矩阵，并通过松弛迭代确定矩阵的参数。此方法的并行程度较低，同时 在处理扰动时准确性较差。
[0007] 上述方法均不具有稳定性，且接口时延会影响仿真的稳定性，不具有通用性。
【发明内容】
[0008] 为解决上述现有技术中的不足，本发明的目的是提供一种吸纳不同步长接口延时 的并行计算分网方法。
[0009] 本发明的目的是采用下述技术方案实现的：
[0010] 本发明提供一种吸纳不同步长接口延时的并行计算分网方法，其改进之处在于， 所述方法包括下述步骤：
[0011]步骤1:确定并入子网的戴维南等值电路的传输线端口方程；
[0012] 步骤2:传输线解耦的不同步长仿真；
[0013] 步骤3:对子网I和子网II采用预置的仿真步长离散化；
[0014]步骤4:对离散化的仿真步长进行仿真。
[0015] 进一步地，所述步骤1包括：
[0016] 电压和电流的时域分析函数描述如下：
(1)
[0017]
[0018]式中：v+表不传输线上的入射波电压，τ表不传输线上的反射波电压;v+和Γ波形 相位一致，增益为f ;ν?ΡΓ波形相位一致，增益为# ; χ表示传输线上任意一点；t表示 当前时刻;V表示行波速度；Γ、Γ分别表示传输线上的入射波电流和反射波电流、Co、Lo分别 表示线路电容和电感；
[0019]忽略研究传输线内的波过程，研究在端口边界的传输线，得到：
[0020]
2)
[0021 ]式中：￥_表示的是从对端端口而来的反射波；I是流出端口的电流，V是端口处的电 压;式(2)表示传输线端口处的戴维南等值电路；
[0022]所述传输线端口方程如下所示：
[0023]
3)
[0024] 把式(2)表示的戴维南等值电路并入子网，求出端口电压VA和VB后，就能通过端口 之间的联络方程求出下一仿真步长的输入戴维南等值电路;V/T、V B_分别表示从对端端口 A 和B而来的反射波;Va+、Vb+分别表不从对端端口 A和B而来的入射波；iA、iB分别表不端口的A 和B的电流;η表示端口包含的节点个数。
[0025] 进一步地，所述步骤2包括:通过传输线解耦以后被分为子网I和II，两个子网通过 传输线L连接在一起;传输线端口在子网I中等效为等值戴维南电路，其电压源为2U f，其阻 抗为A，在子网II中等效为等值戴维南电路，其电压源为2US，其阻抗为Zu
[0026] 假设子网I是快速网络，子网II是慢速网络，通过传输线解耦的延迟时间为慢速子 网的仿真步长;子网I和子网II按照并入接口的戴维南等值电路，建立状态空间表示的状态 方程式如下：
(4)
[0027] (5)
[0028]
[0029]
[0030] (6)
[0031] 式中，Uf+表示子网I传输线的反射电压波，US+表示子网II传输线的反射电压波， 表示子网I传输线的入射电压波，U，表示子网II传输线的入射电压波;乂％表示子网II中的 状态变量，包括电感电流，电容电压和控制器中的积分变量;六3表示子网II中的网络状态矩 阵，Bs和Ds都是系数矩阵；Us是子网I的输入变量，其中U sint表示子网I内部的注入源，其中 Bsint表不Bs矩阵中与Usint相对应的系数矩阵;Cs表不子网II的电容对角矩阵;X ' f表不子网I 中的状态变量，包括电感电流，电容电压和控制器中的积分变量;Af表示子网I中的网络状 态矩阵，Bf和Df均为系数矩阵;Uf是子网I的输入变量，其中U fint表示子网I内部的注入源， Bf int表;^Bf矩阵中与Uf int相对应的系数矩阵，Df int表;^Df矩阵中与Uf int相对应的系数矩阵， Cf表示子网I的电容对角矩阵;Xf为子网I中的状态变量，Xs为子网I中的状态变量。
[0032] 进一步地，所述步骤3包括:按照后退欧拉法进行离散化，并考虑接口方程(7)，得 状态方程式如下(8):
[0033] (7)
[0034] (8)
[0035] 式中，Xf为子网I中的状态变量，t//为子网I中的传输线接口的戴维南等值电压源， Ufint为快速子网的自身电源，Af，Bfint，BfT和Dfint为子网I自身的状态空间表不的参数矩阵， If为传输线端口处的注入电流，Uf+为快速子网I传输线的反射电压波;xs为子网I中的状态 变量，U，为子网I中的传输线接口的戴维南等值电压源，U sint为慢速子网的自身电源，As， Bsint，BsT和Dsint为子网I自身的状态空间表示的参数矩阵，Is为传输线端口处的注入电流，U s +为快速子网I传输线的反射电压波;mh为子网II的仿真步长，h为子网I的仿真步长;DsT和DfT 为传输线两端电压源的关系，为1或者-1。
[0036] 进一步地，所述步骤4包括下述步骤：
[003 7 ] 步骤1 >，按照式(8)把子网II从t=mk-m时刻积分到t = mk时刻，计算出Xs (mk)和传 递给对端接口的电压；
[0038] 步骤2>，按照式(7)把子网I连续积分m步，从t = mk-m时刻积分至ljt = mk时刻，计算 出mk-m到mk之间所有mk-i时刻的值；
[0039] 步骤3>，在mk时刻完成子网I和子网II的信息交换，子网II接收算出的电压值发给 子网I，完成从mk-m到mk的数值积分，并返回步骤1;
[0040] 其中：t表示当前时刻;m表示子网II的步长大小，k表示子网II的走了多少步数，mk 表示子网Π 走了 k个m步、i分别表示中间任意需要计算的时刻;mk+i表示中间插值的时刻； mk和mk+m是跨度为m的两个离散时刻。
[0041]与最接近的现有技术相比，本发明提供的技术方案具有的优异效果是：
[0042]本发明通过实际存在的传输线或者由联络电容和联络电感等效成的传输线把系 统分开的方法。由于传输线模型的本质是隐式梯形积分法的近似，故本发明拥有很高的稳 定性。其次它的接口时延不会影响仿真的稳定性，具有较强通用性。
[0043]本发明能实现完全并行仿真，是实时仿真中理想的仿真方法。应用基于传输线的 不同步长仿真算法需要找到网络中适合分网的传输线，并改造为一个大步长延时的传输线 模型，通过模型把网络解耦。能够实现不同步长电磁暂态并行分网时接口没有延时导致的 计算误差，减少或避免因接口延时造成数值震荡和仿真误差，因为接口延时被模型特性所 消化。
【附图说明】
[0044]图1是本发明提供的传输线TLM模型图；
[0045]图2是本发明提供的传输线TLM等效戴维南电路图；
[0046] 图3是本发明提供的传输线解耦的不同步长仿真示例网络图；
[0047] 图4是本发明提供的不同步长仿真的运算时标图；
[0048] 图5是本发明提供的基于传输线解耦不同步长仿真的流程图。
【具体实施方式】
[0049]下面结合附图对本发明的【具体实施方式】作进一步的详细说明。
[0050]以下描述和附图充分地示出本发明的具体实施方案，以使本领域的技术人员能够 实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例 仅代表可能的变化。除非明确要求，否则单独的组件和功能是可选的，并且操作的顺序可以 变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发 明的实施方案的范围包括权利要求书的整个范围，以及权利要求书的所有可获得的等同 物。在本文中，本发明的这些实施方案可以被单独地或总地用术语"发明"来表示，这仅仅是 为了方便，并且如果事实上公开了超过一个的发明，不是要自动地限制该应用的范围为任 何单个发明或发明构思。
[0051] 本发明提出一种吸纳不同步长接口延时的并行计算分网方法，可以在不同步长电 磁暂态并行计算分网时吸纳接口延时，减少或避免因接口延时造成数值震荡和仿真误差， 包括下述步骤：
[0052] 步骤1:确定并入子网的戴维南等值电路的传输线端口方程；
[0053]发明利用分布参数线路的相关理论，可以知道图1上的电流和电压分布满足了入 射和反射波的叠加。电压和电流的时域分析函数可描述如下：
[0054]
(1)
[0055] 从上式可知，V+表示传输线上的入射波，T表示传输线上的反射波。V+和Γ波形相 位一致，增益为#。Γ和Γ波形相位一致，增益为g。
[0056] 忽略研究传输线内的波过程，研究在端口边界的传输线，可以得到：
[0057]
(2)
[0058] F表示的是从对端端口而来的反射波。I是流出端口的电流，V是端口处的电压;式 (2) 表示传输线端口处的戴维南等值电路，两端口传输线的端口可以等效为图2的等效戴维 南电路。
[0059] 在每一步仿真迭代的过程中，Vi首先影响A网络，然后在边界条件的制约下，输出 一个V+A脉冲给传输线，这一脉冲需要经过一个仿真步长移动到传输线的末端B处，并作为下 一个仿真时步的B网端口的输入脉冲VI』端口的过程与A端口类似。这个过程可以总结为式 (3) 的传输线端口方程。
[0060]
(3)
[0