Please wait a minute...
浙江大学学报(工学版)
J4  2010, Vol. 44 Issue (3): 463-467    DOI: 10.3785/j.issn.1008973X.2010.03.009
自动化技术、计算机技术     
基于无网格有限元方法的心磁场计算仿真研究
李重视1, 闫丹丹1, 朱善安1, Bin He2
1. 浙江大学 电气工程学院,浙江 杭州 310027； 2. 明尼苏达大学 生物医学工程系,  明尼苏达 明尼阿波利斯 55455
LI Zhongshi1, YAN Dandan1, ZHU Shanan1, HE Bin 2
 全文: PDF 
摘要：

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

The magnetic field outside produced by the current dipole within a volume conductor of an arbitrary geometry was calculated by using a meshless finite element method (FEM) and BioSavart law. In order to take the advantage of BioSavart law, the current density within the whole medium has to be calculated first. According to the boundary conditions, the meshless FEM was used to calculate the potential at each node. The current density was achieved by the gradient of potential. The components of the induced field were realized by the integral of current density. The algorithm was evaluated on a singlelayer homogeneous sphere model with a satisfactory result achieved. This algorithm was also used on the hearttorso model from magnetic resonance imaging of a healthy person for the latent application on the analysis of magnetocardiogram (MCG) and magnetoencephalogram(MEG).
Key words:  meshless FEM    magnetocardiography    forward problem    body volume    Bio-Savart law
出版日期: 2010-04-01
:  TP 391.9  
基金资助:

国家自然科学基金资助项目(50577055)
通讯作者: 朱善安,男,教授,博导     E-mail: zsa@zju.edu.cn
作者简介: 李重视（1978—）,男,浙江临安人,博士生,从事生物电仿真计算研究
服务　　
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章　　
李重视
闫丹丹
朱善安
Bin He

引用本文:

李重视, 闫丹丹, 朱善安, Bin He. 基于无网格有限元方法的心磁场计算仿真研究[J]. J4, 2010, 44(3): 463-467.

LI Chong-Shi, YAN Dan-Dan, SHU Shan-An, Bin He. . J4, 2010, 44(3): 463-467.

链接本文:

http://www.zjujournals.com/xueshu/eng/CN/10.3785/j.issn.1008973X.2010.03.009    或    http://www.zjujournals.com/xueshu/eng/CN/Y2010/V44/I3/463

［1］ GREENSITE F. Heart surface electrocardiographic inverse solutions modeling and imaging of bioelectric activity: principles and applications［M］. New York: Kluwer/Plenum. 2004.
［2］ HE B. Electrocardiographic tomographic imaging modeling and imaging of bioelectric activity: principles and applications［M］. New York: Kluwer/ Plenum, 2004.
［3］ HE B, MUSHA T, OKAMOTO Y, et al. Electric dipole tracing in the brain by means of the boundary element method and its solution accuracy ［J］. IEEE Trans Biomed Eng, 1987, 34: 406414.
［4］ AOKI M, OKAMOTO Y, MUSHA T, et al. Threedimensional simulation of the ventricular depolarization and repolarization processes and body surface potentials: normal heart and bundle branch block［J］. IEEE Trans Biomed Eng, 1987, 34(6): 454.
［5］ ZHANG Yingchun, ZHU Shanan, HE B. A highorder finite element algorithm for solving the threedimensional EEG forward problem ［J］. Phys Med Biol, 2004, 49: 29752987.
［6］ HAUEISEN J, RAMON C, EISELT M, et al. Influence of tissue resistivities on neuromagnetic fields and eletric potentials studied with a finite element model of the head ［J］. IEEE Trans Biomed Eng, 1997, 44(8): 727735.
［7］ ROSENFELD M, TANAMI R, ABBOUD S. Numerical solution of the potential due to dipole sources in volume conductors with arbitrary geometry and conductivity ［J］. IEEE Trans Biomed Eng, 1996, 43(7): 679689.
［8］ LUCY L B. A numerical approach to the testing of the fission hypothesis ［J］. The Astron J, 1977, 8(12): 10131024.
［9］ SWEGLE J W, HICKS D L, ATTAWAY S W. Smoothed particle hydrodynamics stability analysis ［J］. J Comput Phys, 1995, 116: 123134.
［10］ CHEN J K, BERAUN J E, JIH C J. An improvement for tensile instability in smoothed particle hydrodynamics ［J］. Comput Mech, 1999, 23: 279287.
［11］ BELYTSCHKO T, GU L, LU Y Y. Fracture and crack growth by elementfree Galerkin methods ［J］. Model Simul Mater Sci Eng, 1994, 2: 307318.
［12］ DURATE C A, ODEN J T. Hp clouds:a meshless method to solve boundaryvalue problems\
[R\]. Austin: University of Texas at Austin, 1995.
［13］ BELYTSCHKO T, KRONGAUZ Y, ORGAN B, et al. Meshless methods: an overview and recent developments ［J］, Comput Methods in App Mech and Eng, 1996, 139(14): 347.
［14］ ZHANG G M, BATRA R C. Modified smoothed particle hydrodynamics method and its application to transient problems ［J］. Computational Mechanics, 2004, 34: 137146.
［15］ ZHANG X, LIU X H, SONG K Z, et al. Leastsquare collocation meshless method［J］. Int J Numer Methods Eng, 2001, 51(9): 10891100.
［16］ von ELLENRIEDER N, MURAVCHIK C H, NEHORAI A, A meshless method for solving the EEG forward problem［J］. IEEE Trans on Biomed Eng, 2005, 52(2): 249257.
［17］ WANG Y, RUDY Y. Application of the method of fundamental solutions to potentialbased inverse electrocardiography［J］ Ann Biomed Eng, 2006,34(8): 12721288.
［18］ LI Z S, ZHU S A,HE B, Solving the ECG forward problem by means of a meshless finite element method［J］.Phys Med Biol, 2007, 52: N287N296.
［19］ CUFFIN BN, COHEN D, Magnetic fileds of a dipole in spherical volume conductor shapes［J］, IEEE Trans Biomed Eng, 1977, 24(4): 373381.
[1] 张俊红, 张玉声, 王健, 徐喆轩, 胡欢, 赵永欢. 考虑热机耦合的排气歧管多目标优化设计[J]. 浙江大学学报(工学版), 2017, 51(6): 1153-1162.
[2] 任迪, 万健, 殷昱煜, 周丽, 高敏. 基于贝叶斯分类的Web服务质量预测方法研究[J]. 浙江大学学报(工学版), 2017, 51(6): 1242-1251.
[3] 宋瑞祥, 张庆国, 于海敬, 徐丽, 施悯悯. 遥感数据的城市不透水面估算及增温效应[J]. 浙江大学学报(工学版), 2017, 51(5): 1051-1056.
[4] 李静, 王哲. 似平面应力条件下混凝土的变形特性[J]. 浙江大学学报(工学版), 2017, 51(4): 745-751.
[5] 蒋鑫龙, 陈益强, 刘军发, 忽丽莎, 沈建飞. 面向自闭症患者社交距离认知的可穿戴系统[J]. 浙江大学学报(工学版), 2017, 51(4): 637-647.
[6] 李明, 刘扬, 唐雪松. 疲劳裂纹的跨尺度分析[J]. 浙江大学学报(工学版), 2017, 51(3): 524-531.
[7] 张捷, 肖新标, 王瑞乾, 金学松. 高速列车铝型材声振特性测试及等效建模[J]. 浙江大学学报(工学版), 2017, 51(3): 545-553.
[8] 李晓东, 祝跃飞, 刘胜利, 肖睿卿. 基于权限的Android应用程序安全审计方法[J]. 浙江大学学报(工学版), 2017, 51(3): 590-597.
[9] 苏星, 王慧泉, 金仲和. 实时高可靠综合电子系统的逻辑架构设计[J]. 浙江大学学报(工学版), 2017, 51(3): 628-636.
[10] 潜龙昊, 胡士强, 杨永胜. 多节双八面体变几何桁架臂逆运动学解析算法[J]. 浙江大学学报(工学版), 2017, 51(1): 75-81.
[11] 张伟, 胡友德, 郑立荣. 基于频率可调驻波振荡器的芯片时钟系统设计[J]. 浙江大学学报(工学版), 2017, 51(1): 168-176.
[12] 谢罗峰, 徐慧宁, 黄沁元, 赵越, 殷国富. 应用双树复小波包和NCA-LSSVM检测磁瓦内部缺陷[J]. 浙江大学学报(工学版), 2017, 51(1): 184-191.
[13] 刘海宾, 王勇, 马鹏磊, 谢玉东. 基于平行式振荡翼系统参数耦合分析[J]. 浙江大学学报(工学版), 2017, 51(1): 153-159.
[14] 厉文榜,方梦祥,岑建孟,肖平,时正海,山石泉,王勤辉,骆仲泱. K-Fe复合催化剂对煤半焦气化速率与产物的影响[J]. 浙江大学学报(工学版), 2016, 50(9): 1746-1751.
[15] 江衍铭,郝偌楠,蔡文柄. 台湾洪水预报进展及模型实务应用[J]. 浙江大学学报(工学版), 2016, 50(9): 1784-1790.