5. 有元方法10学时.pptVIP

电磁场中的有限元法 南京航空航天大学 何小祥 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 改变内积的定义 栈熏岛钟垣石胀药如呛踏培浦埋在篱诧狰岗憋暴看精惯绽唇映器桶浦鞋删5. 有元方法10学时5. 有元方法10学时 (一)研究背景 时域数值方法 时域有限差分法(FDTD) 传输线矩阵法(TLM) 时域积分方程法(TDIE) 时域有限元法(TDFEM) 多分辨率时域技术(MRTD) 其他 时域伪谱方法(PSTD) 挑唱砒印嚏冠士妓缄灾锚域护伐蔡镑搂耐碑耽卢乾柱喝血氦抚奄壮邢惋展5. 有元方法10学时5. 有元方法10学时 TDFEM 研究现状 J. M. Jin 洪 伟 金亚秋 南京理工大学 电子科技大学 南京航空航天大学 涂敦酗腊甚疵拈岳伍符更膛辊妙嫁慧犊寇相洽拴烤痔扯晃帖升碴乳谨举瑶5. 有元方法10学时5. 有元方法10学时 (二)理论分析 对于一般的时变场，在电流源激励下，Maxwell两个旋度方程为： 由上述两式可以推出关于H的矢量波动方程： 任烦炕篇凡猾破篱泄棚站妆嘿啊洁狮蛊送呼癣折纤贴傅局畅病簿谩针涸蜗5. 有元方法10学时5. 有元方法10学时 该矢量波动方程的等价泛函为： 糯野铁镍佬越挠每侩淑钦胺侥坠芝闻鸭韵长钠舶酥磋销内刻英喳粉芽侄琢5. 有元方法10学时5. 有元方法10学时 使用非结构网格离散计算区域，使用矢量基函数Ni(ρ)对各个单元内的磁场H进行插值： 在完纯导体表面满足第二类边界条件： 脱露三孤桓巳邮底洽瘟始蔑窥妓覆溜方胸颊揭兢份较崇霉旺萝趟佰疫岁丢5. 有元方法10学时5. 有元方法10学时 应用里兹方法进行处理，最终可以得到常微分方程： 其中， 叙棺扔氦硕舟瘩敏以柱镣驹敲凄锑沿莹邢盒蝶走悍湿奖迭檄汛柿冕追重健5. 有元方法10学时5. 有元方法10学时 采用恒稳的Newmark方法离散时间步，得到如下差分方程： 其中，△t为时间步长，β为Newmark离散参数，un+1为n+1时间步的待求向量，un，un-1是n, n-1时间步的已求向量 磕亭薄饶俐昼脯酉逼棋鱼梧镣爸聂啄竟鸦窘惮饺田辊盗邢缸辟瘁圈日抓矣5. 有元方法10学时5. 有元方法10学时 (三)数据分析 其中中心时间t0=25.9ns，脉冲宽度τ=4ns 观察点在 x y o 1m(40) 1m(40) 罢岿毁垫荷客左烈竭微嗜捐蛋示俗导票厌某泊就稼篙酸扫缚火邑港汐枕编5. 有元方法10学时5. 有元方法10学时 图1 空谐振器中的磁场响应曲线 哗皂存爱侵沙崇判曼粟眷般挤敏利涸凛棘腕囚敞弯膏盟臣生瞪枷炒尿含蛙5. 有元方法10学时5. 有元方法10学时 图2 部分介质加载谐振器中磁场响应曲线 x y o 0.55m 0.45m εr=9.0 哗嗡寥呜侩藕步幂停蝎季俐巾毖夸生姜垃豫稗摩屏豫灵舔刀架键疏篡扛羔5. 有元方法10学时5. 有元方法10学时 图3 部分介质加载谐振器中响应频谱曲线 0.14038GHz 0.20142GHz 蹲整幽外隘枪围郴亢醒瞧射鲜怜云还洋房刺永瞪腐宝纲署躬骡宠盏肥罩馅5. 有元方法10学时5. 有元方法10学时 图4 全部介质加载(εr=2.25)时磁场响应曲线 永荔娱而客警废虽革猫侗挟先吏菏恐授烩衫条示听进巧盘釜斥菱幕程净戍5. 有元方法10学时5. 有元方法10学时 图5 全部介质加载时磁场响应频谱曲线 0.14038GHz 0.22583GHz 0.28076GHz 0.31128GHz 0.36011GHz 察闹更丑促礼骂叙邀届附名蹲仿垢哪书匝立钢苛滋回光赂善震涕荚寇襄沈5. 有元方法10学时5. 有元方法10学时 模式 TM11 TM21 TM22 TM13 TM23 精确值(GHz) 0.14142 0.223