电阻层析成像有限元节点编号方法研究.pdfVIP

卜………～…………………………一………一E虢矿善撬尝………………………………………一……… 电阻层析成像有限元节点编号方法研究 韦冲，胡继宝，陶娇娇，黄露 (徐州工程学院信电工程学院江苏徐州221111) 【摘要】：有限元网格节点编号规则对整体刚度矩阵的带宽有着非常重要的影响，而正问题计算 效率在很大程度上取决于整体刚度矩阵的带宽值。为了减少电阻层析成像中有限元运算时间，节省计 算机的内存空间，本文在仔细研究了目前相关学者常用的一些编号方法基础上，提出了一种启发式全 新“s”形编号规则，在一定程度上降低了整体刚度矩阵带宽值，减少了数据存储量，大大提高了正问题 的计算效率。 【关键词】：有限元；节点编号；带宽；电阻层析成像 0引言 矩阵K的半带宽定义为： Resistance 电阻层析成像(ElectricalTomography，∥(K)=max(沙，(K)，1≤f≤，2) f31 简称ERT)技术是电学层析成像技术存在的四种基本 半带宽B(K)与矩阵带宽1(K)的关系为： 形式之一，是过程层析成像技术中的一个分支，它与 y(K)=2∥(K)一1 (4) 电阻抗层析成像(Eleetrical 对于带状线性方程组AM嘏=B，其计算时间受 ImpedanceTomography，简 称EIT)技术最重要的一点区别就是在重构场域内部MB(K)2影响【蚓。 的阻抗分布时，电阻层析成像技术只保留了实部信 目前，国内外学者对于该问题的研究普遍采用的 息。因此从一定程度上讲，电阻层析成像技术是电阻 一个方法就是从图形的几何中心节点出发，逆时针方 抗层析成像技术的一种简化形式。该技术中，由于采 向按照层数由内而外进行编号。根据这种编号方法编 用的是安全电流激励，属非侵入式检测技术，在研究 写出的程序有一定的规律性，编程思路也较为清晰， 人体生理功能与疾病诊断方面具有很重要的临床价 但该方法在降低带宽方面的效果并没有从理论上证 值【Ⅲ。 明是最优的。本文提供两种编号思路：①将平面网格 为了提高正问题的计算效率，目前国内外学者在 图形通过平面直角坐标系分成为四块区域，从横轴的 有限元编号优化方面进行了大量的研究。本文首先通 正向开始，逆时针依次为每个象限的节点进行编号。 过直角坐标系对平面区域划分为四个象限，在传统的 ②基于传统编号方案及编号思路①，对所有节点划分 电阻层析成像拓扑结构的节点编号方法上进行改动， 不同的所属带，然后按照“S”形对所有有限元节点按 提供了一种新的节点编号方法，该方法虽然在降低整 照带的顺序依次进行编号。 体刚度矩阵带宽方面的效果不及传统方法，但它提供 2象限划分式编号方法 了一种节点编号方法的思路，并在此基础上研究出一 根据特定网格图形的特点，在传统方法的启发下 种启发式全新编号方案，有效降低了整体刚度矩阵带 进行改动，其八层均匀网格图形的节点编号具体过程 宽。 如下： 1带宽的定义与优化方法 步骤l：首先将圆形网格图形按照平面坐标系划 设有限元网络系统的节点总数为n，有限元方程 分为四个象限。其中横轴和纵轴的正半轴分属于第一 组的整体刚度矩阵K是一个对称方阵，定义fi(K)为与第二象限，订方向上的节点属于第三象限，3州2方 K中的第i行最左边非零元素的下标，那么有： 向上的节点属于第四象限。 Z(K)=min{，：％≠o) (1) 步骤2：借鉴传统的逆时针由内而外按层编号方 y，(必)=f—Z(K) (2) 法，本方法亦采用逆时针顺序编号。但基于象限的编