二维双线性浸入式有限元研究-力学专业论文.docxVIP

Classified Index: O242.21 U.D.C: 52 Dissertation for the Master Degree in Engineering RESEARCH OF BILINEAR IMMERSED FINITE ELEMENTS FOR INTERFACE PROBLEMS Candidate： Zhang Xiaoshi Supervisor： Associate Prof. Zheng Liang Academic Degree Applied for： Master of Engineering Speciality： Mechanics Affiliation： Shenzhen Graduate School Date of Defence： June, 2013 Degree-Conferring-Institution： Harbin Institute of Technology 摘 要 在求解二维浸入式有限求解电磁场问题时，对于仿真区域通常采用三角 形单元或四边形单元，为了保证计算的精确度，往往要将仿真区域划分的非 常细，如果对整个区域都进行划分，至少要有几万甚至到几百万的网格 单元， 再加上多次循环，这样的计算将很费时并且需要巨大的内存空间， 受到计算 机内存容量的限制将无法进行计算。 对于二维双线性浸入式有限元求解泊松 方程问题，本文在 Dirichlet 边界条件和 Neumann 边界条件的基础上构造了一 个新的边界条件，即 周期性边界条件并且解决了周期性边界条件问题。对于 有周期性的特点的电场，可以找到区域的对称线，根据其周期性变化，可以 利用周期性边界条件，取一个周期计算即可，这样不仅加快了计算速度，也 减少了内存空间的需求。 对于电磁场接触面问题中，导体是电磁场中的一个主要的物体，在等离 子的环境中可以被充电。对于一个接触面问题，为了获得更好的收敛效果通 常是使用有限元法和有限差分法并采用贴体网格。双线性浸入式有限元在直 角坐标系下不仅可以获得同样的收敛效果，而且 无需重复的对网格进行划分 。 二维双线性浸入式有限元是 可以有效解决接触面问题。本文建立 了一个适用 于物体是导体的双线性浸入式有限元方法，并且通过数值仿真模拟对新算法 进行了分析。这种新方法在仿真区域可以获得精确的结果。 本文最后构建了二维双线性轴对称浸入式有限元 方法。二维双线性轴对 称浸入式有限元可用于求解三维轴对称接触面问题，并且 可以把把三维问题 简化为二维问题求解。应用里兹法证明了 建立的模型，然后通过 Fortran 编程 使得二维双线性浸入式有限元算法在计算机上实现，最后 进行数值仿真求解。 通过计算解析解与理论解求解误差的 L2 ，H1 范数，并通过做出误差的线性回 归图，可以看出 L2 ， H1 范数分别以 O(h2 ) 和 O(h) 速率收敛，表明仿真结果是 正确的。 关键词：双线性；浸入式有限元；周期性边界条件；导体；二维轴对称 Abstract For the applications with moving interfaces, the mesh has to be refined again and again according to the varying interface, which prevents the conventional finite difference and finite element method from working efficiently. For large scale problems, the computational expense for ref ining the mesh again and again is prohibitive. Therefore, many methods have been developed to get rid of this limitation so th at interface problems can be solved with a mesh independent of interface. A periodic boundary condition is developed based on Dirichlet boundary condition and Neumann boundary condition. The periodic boundary condition can solve the interface field boundary value problem on a Cartesian mesh. It also requires