有限元三维静电场算法验证

简介

上次设计了一个小例子，做了有限元二维静电场算法的验证。这次设计了一个简单的三维算例，验证静电场三维算法。之所以选择简单的算例，主要是为了方便调试，验证算法及程序正确性后，复杂算例与简单算例在底层公式算法上没有本质区别。

理论

注：上述体单元插值均基于通用的三维四面体网格。

几何模型

三维算例如下图所示，其中内导线电位为 1V，外导线电位为 0V，导线之间为一层介质。

网格剖分

依然使用开软程序 Gmsh 进行三维四面体网格剖分。

计算结果

z=0 切面电位云图：

y=0，z=0 沿线电位：


算法结果已经与成熟商软进行了对比，相关对比结果后续补充。

源代码

程序源码见T2.m：
https://gitee.com/YoHen-Fu/simag-doc/tree/master/example