论文部分内容阅读
在求解二维浸入式有限求解电磁场问题时,对于仿真区域通常采用三角形单元或四边形单元,为了保证计算的精确度,往往要将仿真区域划分的非常细,如果对整个区域都进行划分,至少要有几万甚至到几百万的网格单元,再加上多次循环,这样的计算将很费时并且需要巨大的内存空间,受到计算机内存容量的限制将无法进行计算。对于二维双线性浸入式有限元求解泊松方程问题,本文在Dirichlet边界条件和Neumann边界条件的基础上构造了一个新的边界条件,即周期性边界条件并且解决了周期性边界条件问题。对于有周期性的特点的电场,可以找到区域的对称线,根据其周期性变化,可以利用周期性边界条件,取一个周期计算即可,这样不仅加快了计算速度,也减少了内存空间的需求。 对于电磁场接触面问题中,导体是电磁场中的一个主要的物体,在等离子的环境中可以被充电。对于一个接触面问题,为了获得更好的收敛效果通常是使用有限元法和有限差分法并采用贴体网格。双线性浸入式有限元在直角坐标系下不仅可以获得同样的收敛效果,而且无需重复的对网格进行划分。二维双线性浸入式有限元是可以有效解决接触面问题。本文建立了一个适用于物体是导体的双线性浸入式有限元方法,并且通过数值仿真模拟对新算法进行了分析。这种新方法在仿真区域可以获得精确的结果。 本文最后构建了二维双线性轴对称浸入式有限元方法。二维双线性轴对称浸入式有限元可用于求解三维轴对称接触面问题,并且可以把把三维问题简化为二维问题求解。应用里兹法证明了建立的模型,然后通过Fortran编程使得二维双线性浸入式有限元算法在计算机上实现,最后进行数值仿真求解。通过计算解析解与理论解求解误差的L2,H1范数,并通过做出误差的线性回归图,可以看出L2,H1范数分别以2O(h)和O(h)速率收敛,表明仿真结果是正确的。