基于无网格有限元方法的心磁场计算仿真研究

摘要

利用无网格有限元方法求解毕奥一萨伐定律来计算电流偶极子源在任意形状均质容积导体外产生的磁场分布,通过毕奥-萨伐定律计算容积导体外产生的磁场分布须先得到整个媒质中的电流密度分布,依据相应的边界条件,采用无网格有限元方法求解媒质容积导体有限元模型中各节点的电势;通过对电势求梯度得到媒质内部的电流密度分布;媒质内部感应磁场分布的各个分量可通过对电流密度进行数值积分计算来实现. 为了验证提出的无网格有限元方法的有效性,在单层均质球模型上进行了仿真研究,并将该方法应用于从核磁共振成像(MRI)得到的真实形状的心脏-躯干心脏模型上,取得了令人满意的计算结果,表明了算法的有效性及其在心磁图(MEG)和磁场分布(MCG)上的潜在应用前景.