有限元求解二维静电场问题

简介

笔者开源了正在编写的SIMag理论文档以及部分MATLAB验证程序，完整的求解器不开源。这篇博文摘取了SIMag二维静电场求解的理论文档，并设计了一个算例，编制了MATLAB程序验证算法。

笔者自知水平有限，欢迎各位朋友批评斧正。

理论

几何模型

利用Gmsh建立三轴电缆二维模型，其中Circle 2、Circle 3、Circle 4为电缆。

网格剖分

利用Gmsh将上述的二维几何模型进行剖分：

可以看到有四组线单元组，分别对应三个三轴电缆以及外侧的壳，我们对四组线单元组施加第一类边界条件，其中组2施加电位为1V，其它边界的电位为零

计算结果

编制MATLAB程序求解该二维静电场问题，计算得到得到电位云图为：

源代码

源代码见：

https://gitee.com/YoHen-Fu/simag-doc/tree/master/example

IP矩阵：斗沙片刻/斗沙片刻H

CSDN主页、知乎主页、gitee主页、github主页、bilibili
微信公众号(斗沙片刻H):