有限元方法10学时课件.ppt

电磁场中的有限元法电磁场中的有限元法 南京航空航天大学南京航空航天大学 何小祥何小祥 Outline有限元发展历史有限元发展历史有限元的优点有限元的优点二维节点有限元（有限元处理过程）二维节点有限元（有限元处理过程）二维高阶节点有限元二维高阶节点有限元三维节点有限元方法三维节点有限元方法三维高阶节点有限元方法三维高阶节点有限元方法节点有限元的缺点节点有限元的缺点矢量有限元方法矢量有限元方法二维、三维一阶矢量基函数二维、三维一阶矢量基函数高阶矢量基函数高阶矢量基函数有限元中的现代技术有限元中的现代技术时域有限元方法时域有限元方法有限元方程组求解有限元方程组求解-CG1、有限元发展历史、有限元发展历史2、有限元的优点、有限元的优点非结构网格-可以模拟任意形状的边界泛函的思想-利于复杂介质特性模拟有限元处理过程有限元处理过程边值问题等价泛函目标建模，网格剖分参数提取矩阵填充局部矩阵计算FEM方程组求解后处理技术差值基函数只需差值基函数只需要在单元内满足要在单元内满足精度要求精度要求 U代表位函数或者某个场分量12一般计算区域比较复杂，积分不容易实现，采用整体基函数很难满足齐次边界条件，采用了广义变分原理3、二维节点有限元（有限元处理过程）、二维节点有限元（有限元处理过程）xyo3Patran Ansys所需参数所需参数-坐标，单元编号，节点或棱边编号坐标，单元编号，节点或棱边编号4将波导横截面用M个三角形面元离散，则各面元中的u 用线性插值函数表示为积分区域变为三角形单元内思想、思想、理解理解能量的概念5物理含义123参数提取优点优点利用横等式可以获得局部矩阵的解析解6单元一、节点2.4.1方程组矩阵填充方程组矩阵填充7（1）（2）（3）单元二、节点2.4.5单元三、节点3.5.2单元四、节点4.5.6单元节点与相邻单元节点必须吻合，过程中使用到各个单元中节点局部编号与整体编号的关系，列向量等可以如上进行填充最后形成的系数矩阵是一个对称的稀疏矩阵最后形成的系数矩阵是一个对称的稀疏矩阵可导出线性代数方程组由此方程解出本征值和对应的本征矢。984、二维高阶节点有限元二维高阶节点有限元当插值函数为高次多项式时，用自然坐标表示的单元一阶基函数二阶基函数三阶基函数 经过化简其中可以事先制定成表格，程序中调用波型 (m)0.0460.0230.0200.01829解析解 (m)0.0460.0230.0200.01834误差0.02%0.35%0.08%0.26%波导尺寸为10mm*23mm。分格数为4*3，应用2阶FEM分析得到特征值如下，三角单元总数为24，场节点总数为63 谐振腔尺寸为12mm*28mm时，应用FEM进行分析，划分网格数为5*4，三角单元总数为40，场节点总数为99。波型 (m)0.0560.0280.0240.02204解析解 (m)0.0560.0280.0240.02207误 差0.11%0.15%0.25%0.15%特征值特征值计算误差计算误差5、三维节点有限元方法、三维节点有限元方法6、三维高阶节点有限元方法、三维高阶节点有限元方法略7、节点有限元的缺点、节点有限元的缺点标量有限元法的缺标量有限元法的缺点：点：计算结果中会出计算结果中会出现不易辨别的非现不易辨别的非物理解或伪解；物理解或伪解；在介质分界面或在介质分界面或导体表面不易强导体表面不易强加边界条件；加边界条件；由于场存在奇异由于场存在奇异性，处理导体及性，处理导体及介质边缘和棱边介质边缘和棱边有困难。有困难。9 二维矢量有限元方法二维矢量有限元方法10 三维矢量有限元方法三维矢量有限元方法矢量有限元性态差矢量有限元性态差11、高阶矢量基函数、高阶矢量基函数略12、有限元中的现代技术、有限元中的现代技术PML MethodDDM MethodFEM-BIFEM-PO，IPO etc.13、时域有限元方法、时域有限元方法(一)研究背景时域数值方法时域有限差分法(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 观察点在 xyo1m(40)1m(40)图图1 空谐振器中的磁场响应曲线空谐振器中的磁场响应曲线图图2 部分介质加载谐振器中磁场响应曲线部分介质加载谐振器中磁场响应曲线 xyo0.55m0.45mr=9.0图图3 部分介质加载谐振器中响应频谱曲线部分介质加载谐振器中响应频谱曲线 0.14038GHz0.20142GHz图图4 全部介质加载全部介质加载(r=2.25)时磁场响应曲线时磁场响应曲线图图5 全部介质加载时磁场响应频谱曲线全部介质加载时磁场响应频谱曲线 0.14038GHz0.22583GHz0.28076GHz0.31128GHz0.36011GHz模式精确值(GHz)仿真值(GHz)误差(%)TM110.141420.140380.735TM210.223610.225830.993TM220.282840.280760.735TM130.316230.311280.125TM230.360560.360110.125图图6 CG收敛曲线对比收敛曲线对比14 FEM求解器求解器-CG等同于最小余量法了一般频域FEM的系数矩阵不具备自伴性位了求解（11.53）泛函得到如下线性方程组但是一、基函数p未知，二、我们为了求原来的方程组，必须求解新的方程组，无意义初始剩初始剩余向量余向量向量系向量系数数构造新构造新解解构造新的构造新的剩余向量剩余向量构造新的构造新的基函数向基函数向量量Comments on CG 迭代方法-理论上N步收敛单调收敛与矩阵性态有关-条件数-谱半径-最大特征值与最小特征值比预处理技术-SSOR-ICCG改变内积的定义CG和BICG方法的求解速度与矩阵的性态即条件数由很大关系。通过预处理技术可以有效提高矩阵的性态，从而提高计算速度