一、划分方法
有限元积分法FIT(Finite Integration Technique,FIT)与传输线矩阵法(Transmission Line Matrix,TlM),这两种方法都可以划分六面体网格:
有限积分技术 (Finite Integration Technique, FIT)
基本概念:
- FIT是一种数值方法,用于求解Maxwell方程组。它通过在离散化网格上应用有限积分的方式来模拟电磁场。
- 离散化:FIT将整个仿真区域划分为若干小的单元(如立方体、长方体),并在这些单元上定义电场和磁场。
- Maxwell方程:FIT直接离散Maxwell方程,保持了其原有的形式。这使得FIT在处理复杂几何结构和材料时具有高精度。
优点:
- 精确性高:由于直接离散Maxwell方程,FIT在模拟复杂几何和材料时具有较高的精确性。
- 适用范围广:可以处理不同形状的物体和各种材料,包括各向异性和非均匀材料。
- 稳健性好:在数值计算中表现出很好的稳健性,不易出现数值不稳定问题。
缺点:
- 计算量大:由于需要对整个仿真区域进行离散化,计算量和存储需求较大,特别是在高频和大规模仿真中。
- 复杂度高:对用户来说,设置和使用可能较为复杂,需要一定的学习和操作经验。
应用场景:
- 高频电磁仿真,如天线设计、微波电路和器件仿真。
- 电磁兼容性(EMC)分析。
- 电磁散射和雷达截面(RCS)计算。
传输线矩阵法 (Transmission Line Matrix, TLM)
基本概念:
- TLM是一种基于传输线网络模拟电磁场传播的数值方法。它将仿真区域划分为网格,每个网格单元用等效的传输线表示。
- 传输线模型:在TLM中,电磁波的传播被模拟为在传输线网络中的电压和电流脉冲传播。
- 散射和传输:TLM通过计算每个网格节点的电压和电流散射,模拟电磁波在整个网络中的传播和相互作用。
优点:
- 直观性强:TLM基于传输线模型,容易理解和实现,特别适用于模拟电磁波在简单介质中的传播。
- 灵活性好:可以很容易地处理多种边界条件和材料属性,特别是对于非均匀介质和非线性材料。
- 适用性强:在时域分析中表现良好,适用于瞬态响应和宽带分析。
缺点:
- 计算精度:在处理复杂几何和高频问题时,计算精度可能不如FIT高。
- 网格要求:TLM方法对网格的要求较高,网格不均匀性可能影响计算精度和效率。
- 内存需求:对大规模三维问题,内存需求和计算量较大。
应用场景:
- 瞬态电磁仿真,如电磁脉冲(EMP)和雷电效应分析。
- 时域电磁兼容性(EMC)分析。
- 宽带天线和电磁波传播仿真。
有限元方法(Finite Element Method,FEM),常用于划分四面体网格和平面网格:
有限元方法(Finite Element Method,FEM)
-
离散化:
- 将连续的计算域划分为有限数量的简单几何形状(如三角形、四边形、六面体等),称为单元(element)。
- 在每个单元上定义基函数,通过这些基函数的线性组合近似求解域内的变量。
-
基函数:
- 基函数通常是低阶多项式,如线性、二次或三次多项式。
- 基函数的选择影响解的精度和计算复杂度。
-
组装:
- 将所有单元的局部方程组装成整个计算域的全局方程。
- 这一步通常涉及矩阵的拼接和求解,形成一个大型稀疏矩阵。
-
求解:
- 使用数值线性代数方法(如迭代法或直接法)求解全局方程,得到变量的近似值。
-
后处理:
- 对解进行后处理,如可视化、计算派生量(应力、流量等),并进行结果分析和验证。
二、网格设置的条件
(1)网格设置的充分条件(强条件)(无条件记住!!!)
a)满足每个波长20-35条网格线
b)结构细节一定要分辨
c)一个网格种满足三种情况:一种介质、两种介质、三种介质(空气-介质-空气)
(2)网格设置的充分条件(弱条件)(无条件记住!!!)
a)满足每个波长20-35条网格线
b)场强强处结构细节一定要分辨
c)场强强处无短路
e.g:
对于同轴线:内外导体之间至少要有一个网络,部分填充网络由PBA级数处理
对于微带线:纵向上取决于微带线的厚度,横向上微带宽度至少需要两个网格
对于差分线:两条走线之间至少要有一天网格线
对于离散端口:至少要有一个网格
场强强处无短路:
"cc"问题(即一个网格被PEC填充了):
- 网络粗糙造成PEC填充,产生短路
- 必须避免此情况
三、六面体网格的设置步骤
- 选择合适的模板(e.g 平面天线与波导天线,平面天线的网格加密要更密一些)
- 六面体网格手动设置三部曲
- 对于微带结构,设置35/35/50;对于波导结构,设置25/25/30
- 最小网格(min,mesh step)的最大化。此数越小,仿真时间越长。此数小一倍,则在所有其他条件不变的情况下,仿真时间增加一倍;相反地,此数大一倍,仿真时间小一倍。所以此数越大越好!若次数与关键结构中最小间距可以比拟,则设置完毕。若此数远小于结构的最小间距,则必须采用限制最小网格步长(Smallest mesh step)的选项。最小网格步长的设置原则:等于仿真关键结构中最小尺寸乘以0.5至0.8.如最精细结构为0.25,则需设置最小网格步长应在0.17左右。
- 设置最小网格。仅改变了次数为0.17,就使得总网格数下降了5倍, 加之最小网络增大了1.7倍,右边设置的仿真速度将比左边快8.5倍。
- 在上述设置完成后,倘若仍然有大量的细小网络,这些网络可能很接近Min,mesh step, 则可以进入网格属性中的special菜单,打开Fixpoints页面,将默认的1000改小。此设置慎用。
四、时域自适应网格加密
- 仿真频率范围设置。定义两个变量:设计频率上下限为f_high和f_low。在Freq Range Settings对话框中以Fmin=f_low/1.2 ,Fmax=f_high*1.2的方式设置仿真频率。即使得仿真频带略宽于设计指标所规定的频带。(由于边带不准,边界上收敛不佳)
- 将自适应网格加密的频率按照如下设置限制在指标频率的上下限内,以避免时域激励信号在信号边带附近收敛不佳的问题,从而影响到自适应网格终止条件的满足,使仿真时间加长。
五、六面体网格设置的流程(略)
六、六面体TLM网格以及具体模型
时域传输线矩阵算法(TLM):
- TlM使用立体的传输线网络来描述电磁场的空间离散模型,作为激励的电压脉冲沿着这些传输线传播,并在节点处发生反射。
- TLM vs FITD:
- 两者均是长方体六面体网络(HEX)
- FITD的HEX是像魔方一样,以每个小长方体为单位,电磁场在相邻长方体间传播
- TLM的HEX像钢筋混凝土中的钢筋一样,以每条拉通的线为单位
- TLM的HEX的长宽比可以很大,如100甚至上千,而FITD的最好在100以内
- TLM支持精简模型库,特别适用于EMC仿真;FITD适用于其他应用,比如天线、RCS等等
- 精简模型:
- 在EMC仿真中,很多模型的细节对仿真结果至关重要:
- 缝隙、搭接、线缆、通风孔、燕尾槽、导电橡胶、螺钉等等
- 可以使用非常精细的网格来划分这些细节,但会导致:
- 仿真时间长
- 占用很大内存
- 使用精简模型来置换这些细节可以极大地提高对这些精细结构的仿真效率:
- 用等效电路模型来描述孔缝、导线、屏蔽电缆和导电涂覆等
- 在EMC仿真中,很多模型的细节对仿真结果至关重要: