[发明专利]基于EB连续‑不连续伽辽金混合的时域有限元方法

权利要求书

1.一种基于EB连续-不连续伽辽金混合的时域有限元方法，其特征在于包含如下步骤：

第一步，建立求解模型，使用四面体网格对模型进行离散，得到模型的结构信息，包括四面体的节点信息以及单元信息；

第二步，执行程序前处理，将模型通过八叉树进行分组，确定每个子区域内单元信息，循环所有子区域，对连续区域和不连续区域的棱边以及面进行整合编码；把模型任意分成N个子区域，每个子区域由几个相连的单元组成，N个子区域内采用连续有限元方法，子区域之间的数值信息传递由区域外表面上强加的不连续伽辽金的数值通量实现；

第三步，从一阶麦克斯韦旋度方程出发，对等式两边采用伽辽金法测试，对电场和磁通量用基函数展开，引入电场和磁通量的连续性条件，化简之后的方程使用蛙跳格式展开，得到最终的迭代公式；通过连续-不连续伽辽金时域有限元的迭代公式进行时间迭代，迭代结束得到空间中的电场和磁通量值；

第四步，数据后处理，根据计算出的场值提取相关的物理参数。

2.根据权利要求1所述的基于EB连续-不连续伽辽金混合的时域有限元方法，其特征在于：所述步骤三中，根据一阶麦克斯韦旋度方程：

ϵ∂E∂t=1μ▿×B---(1)]]>

∂B∂t=-▿×E---(2)]]>

上式中ε、μ分别表示离散单元的介电常数和磁导率，E和B分别表示电场强度和磁感应强度；

计算子区域应用连续伽辽金技术，并将E和B用基函数展开，最终得到

[Te]∂e∂t=[Pe]b---(3)]]>

[Tb]∂b∂t=[Pb]e---(4)]]>

计算子区域之间的外表面上应用不连续伽辽金技术，在交界面处强加迎风通量连续性条件，并将E和B用基函数展开，最终得到

[Te]∂e∂t=[Pe]b+[Se]b+[Sse]b++[Se1]e+[Sse1]e+---(5)]]>

[Tb]∂b∂t=[Pb]e+[Sb]e+[Ssb]e++[Sb1]b+[Ssb1]b+---(6)]]>

(3)、(4)和(5)、(6)耦合在一起矩阵求解，得到的最终矩阵形式与(5)、(6)式一样；其中e和b分别代表本单元中的电场值和磁通量值，e⁺和b⁺分别代表相邻单元的电场值和磁通量值；[T_e]、[T_b]、[P_e]、[P_b]、[S_e]、[S_b]、[Ss_e]、[Ss_b]、[Se1]、[Sb1]、[Sse1]、[Ssb1]分别为形成的矩阵。