生物组织中任意倾斜多柱状电极的电场计算.docxVIP

生物组织中任意倾斜多柱状电极的电场计算 针灸和针灸广泛应用于医学。这是因为针灸和针灸可以临床上的伤口小，使用灵活，可以治疗浅表、深部肿瘤和内脏器官病变。它是一种有效的补救措施，逐渐成为癌症和其他疾病综合治疗的有效补充手段。由于单针电极的时间分布是不均匀的[10.12]，治疗组织非常大，需要扩大治疗范围，通常需要增加针灸和脊柱针的数量。同时，为了杀死被破坏的组织的薄细胞，在杀死组织中需要强烈的横向分布，以最小的损伤。为了达到不同的治疗效果和目的，对多针倾斜电极的治疗要求非常严格。获得的文献主要是单元网格阻抗法和有限元法。在现有的限元法中，只使用了两个针电极的计算示例。多针电极的计算示例是二维平面模型，垂直针灸和莫哈布西翁系统[9.13]。近年来，人们对针灸和莫哈布西翁系统空间分布的定量研究较少。在文献中，主要有单元网格阻抗法和有限元法。在有限元法中，只使用了两个针电极。多针电极的计算示例是二维平面模型和垂直针灸电极[9.13]。近年来，人们主要使用积分法原理对多针电极的垂直分布进行了计算，但目前还不够。此外，还有许多计算针灸中产生的电流的计算方法，这些算法主要基于简化模型而使用电路理论。这是接近计算方法或对称电离作用的问题。然而，如果这种设计的电离作用是不规则的，或多针倾斜电极的电离作用是困难的，它是无法计算复杂电离作用的。现在，医学上有许多关于应用的报道，但许多关于临床实践报告的定性描述。在临床上，为了改变电压和电极距离来改变治疗效果，人们需要对多针倾斜电极的场分布进行定量分析和计算。这对医学的继续应用非常重要。此外，对于多针倾斜电极的动态分布特征进行了进一步的探索。 在数值计算中所使用的有限元或有限差分法, 其缺点是方程涉及的未知量太多, 使计算的速度和精度受到限制.若用于计算倾斜多针状电极系的场分布, 把针状电极当柱状处理, 则网格需划分的非常小;如果把针状电极当作线状处理, 则会出现奇点;倾斜针状电极的存在还会出现不规则立体网格的划分问题, 利用有限元或有限差分法计算倾斜多针状电极系的场, 在边界的处理上存在困难, 因此目前尚无倾斜多针状电极系所产生场的算例.积分方程法所涉及的未知量只分布在不同介质的边界上, 可使计算量减小, 从而较大幅度地提高计算速度[23~32].同时积分方程法将电流分布按线状处理, 可以很好地解决用有限元素或有限差分法计算倾斜多针状电极系边界处理所遇到的困难.但是用积分方程法计算静电类问题时, 有时离散方程会存在病态问题.我们利用算例对边界元积分方程系数矩阵进行了考查, 所选算例属于良态方程.本文实现了利用积分方程法计算倾斜多针状电极系的电场分布, 并获得了较快的计算速度. 1 计算电流分布的性能 1.1 含多个含电极的x-u-t-ln 点电流源I在生物导电组织中任意一点的电势可以写为 其中r为场点到源点的距离, σ是生物组织的电导率. 图1为任意多针状电极的不规则空间分布图, 若有L1, L2, …, Lm, …, Ln个针状电极 (编号) , 并用L1, L2, …, Lm, …, Ln代表针状电极的长度.设由单位长度针状电极流入生物组织的电流为I1, I2, …, Im, …In, 且I1, I2, …, Im, …, In在同一个电极上的分布是不均匀的、随空间坐标而变化, 则有 其中U=[U1, U2, L, Um, L, Un]T, U1, U2, …, Um, …, Un分别为针状电极上的电势, 上标T表示矩阵的转置, r1, r2, …, rm, …, rn分别为场点到源点的距离. 1.2 转化单元i,j 设场点和源点都在针状电极上, 将L1, L2, …Lm, …, Ln个电极分别分为M1, M2, …, Mm, …, Mn个网格, 采用块中心网格系统, 总网格数为M=M1+M2+L+Mm+L+Mn, 则第j个网格在第i个网格所产生的电势为 其中ir是网格编号为i的中心点的矢径, r′为第j个网格上任意一微分元dl的矢径.若i和j分别为多针状电极系中任意两针状电极Lm1和Lm2上的任意两个网格 (i和j也可以在同一根针状电极上, 如图2) , 设任意两针状电极Lm1和Lm2下端点的坐标为 (x01, y01, z01) 和 (x02, y02, z02) , 第i个网格中心点的坐标为 (x1, y1, z1) , 针状电极Lm1和Lm2的方向余弦为 (cosα1, cosβ1, cosγ1) 和 (cosα2, cosβ2, cosγ2) , α11β, γ1和α2, β2, γ2分别为针状电极Lm1和Lm2与坐标轴x, y, z的夹角.在第j个网格中认为I (j) 为常数, 若I (j) 系数用a (i, j) 表示, 则有 其中 在这里il和lj表示由针状电极Lm1和Lm2的