有限元方法学时.ppt

电磁场中的有限元法 南京航空航天大学 何小祥 Outline 有限元发展历史 有限元的优点 二维节点有限元（有限元处理过程） 二维高阶节点有限元 三维节点有限元方法 三维高阶节点有限元方法 节点有限元的缺点 矢量有限元方法 二维、三维一阶矢量基函数 高阶矢量基函数 有限元中的现代技术 时域有限元方法 有限元方程组求解-CG 1、有限元发展历史 2、有限元的优点 非结构网格-可以模拟任意形状的边界 泛函的思想-利于复杂介质特性模拟 4、 二维高阶节点有限元 5、三维节点有限元方法 6、三维高阶节点有限元方法 略 7、节点有限元的缺点 11、高阶矢量基函数 12、有限元中的现代技术 PML Method DDM Method FEM-BI FEM-PO，IPO etc. 13、时域有限元方法 14 FEM求解器-CG Comments on CG 迭代方法-理论上N步收敛 单调收敛 与矩阵性态有关-条件数-谱半径-最大特征值与最小特征值比 预处理技术-SSOR-ICCG 改变内积的定义 (一)研究背景 时域数值方法 时域有限差分法(FDTD) 传输线矩阵法(TLM) 时域积分方程法(TDIE) 时域有限元法(TDFEM) 多分辨率时域技术(MRTD) 其他 时域伪谱方法(PSTD) TDFEM 研究现状 J. M. Jin 洪 伟 金亚秋 南京理工大学 电子科技大学 南京航空航天大学 (二)理论分析 对于一般的时变场，在电流源激励下，Maxwell两个旋度方程为： 由上述两式可以推出关于H的矢量波动方程： 该矢量波动方程的等价泛函为： 使用非结构网格离散计算区域，使用矢量基函数Ni(ρ)对各个单元内的磁场H进行插值： 在完纯导体表面满足第二类边界条件： 应用里兹方法进行处理，最终可以得到常微分方程： 其中， 采用恒稳的Newmark方法离散时间步，得到如下差分方程： 其中，△t为时间步长，β为Newmark离散参数，un+1为n+1时间步的待求向量，un，un-1是n, n-1时间步的已求向量 (三)数据分析 其中中心时间t0=25.9ns，脉冲宽度τ=4ns 观察点在 x y o 1m(40) 1m(40) 图1 空谐振器中的磁场响应曲线 图2 部分介质加载谐振器中磁场响应曲线 x y o 0.55m 0.45m εr=9.0 图3 部分介质加载谐振器中响应频谱曲线 0.14038GHz 0.20142GHz 图4 全部介质加载(εr=2.25)时磁场响应曲线 图5 全部介质加载时磁场响应频谱曲线 0.14038GHz 0.22583GHz 0.28076GHz 0.31128GHz 0.36011GHz 0.125 0.36011 0.36056 TM23 0.125 0.31128 0.31623 TM13 0.735 0.28076 0.28284 TM22 0.993 0.22583 0.22361 TM21 0.735 0.14038 0.14142 TM11 误差(%) 仿真值(GHz) 精确值(GHz) 模式 图6 CG收敛曲线对比 等同于最小余量法了 一般频域FEM的系数矩阵不具备自伴性 位了求解（11.53）泛函 得到如下线性方程组 但是一、基函数p未知，二、我们为了求原来的方程组，必须求解新的方程组，无意义 谐振腔尺寸为12mm*28mm时，应用FEM进行分析，划分网格数为5*4，三角单元总数为40，场节点总数为99。 0.15% 0.25% 0.15% 0.11% 误 差 0.02207 0.024 0.028 0.056 解析解 (m) 0.02204 0.024 0.028 0.056 (m) 波型 特征值 计算误差 标量有限元法的缺点： 计算结果中会出现不易辨别的非物理解或伪解； 在介质分界面或导体表面不易强加边界条件； 由于场存在奇异性，处理导体及介质边缘和棱边有困难。 9 二维矢量有限元方法 10 三维矢量有限元方法 矢量有限元性态差 略 * * * * 有限元处理过程 边值问题 等价泛函 目标建模，网格剖分 参数提取 矩阵填充 局部矩阵计算 FEM方程组求解 后处理技术 差值基函数只需要在单元内满足精度要求   U代表位函数或者某个场分量 1 2 一般计算区域比较复杂，积分不容易实现，采用整体基函数很难满足齐次边界条件，采用了广义变分原理 3、二维节点有限元（有限元处理过程） x y o 3 Patran Ansys 所需参数-坐标，单元编号，节点或棱边编号 4 将波导横截面用M个三角形面元离散， 则各面元