传感器有限元计算中的自动分割技术.pdf

电子测量与仪器学报 2000增刊 传感器有限元计算中的自动分割技术 陈书旺 庄庆德 曲国明 (河北科技大学电子信息工程系，石家庄050054) 【摘要1本文阐述了一种用有限元法来计算传感器特性的元素自动分割方法.它能按用户指定的元素密度分 布函数自动地生成节点与元素，井且对于不适当的分创有纠错功能。 关位词:传感器 有限元 分]PPJ 挤协: 1引言 传感器广泛应用于工Ati踌斗学研究中，它的信息捕获功能，在自动检测、自动控制、遥 感遥测等领域中起着决定性作用。传感器的种类萦多，既有物质微观结构又有材料宏观特性，技 术更新层出不穷，结构设计千差万别。正是由于传感器的这些特点，使得在传感器的设计及制造 中有相当大的困难。 为解决这一问题，我们采用了有限元数值计算方法来计算传感器的特性。其中，如何对传感 器整体进行自动分割，是一个重要的课题。我们所研究的分割方法可以自动地按用户指定的密度 分布进行.从而为传感器的计算及最佳设计提供了一种有效的方法。 2计算方法 有限元法的基本思想是将连续的求解区域离散为一组有限个且按一定方式相互联结在一起 的单元的组合体。由于单元能按不同的联结方式进行组合。且单元本身又可以有不同形状，因此 可以模型化几何形状复杂的求解域。有限单元法作为数值分析方法的另一个重要特点是利用在每 一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。因此，一个问题的有限 元分析中，未知场函数及其导数在各个节点上的数值就成为新的未知量，从而使一个连续的无限 自由度问题变成离散的有限自由度问题。一经求解出这些未知量.就可以通过插值函数计算出各 个单元内场函数的近似值，从而得到整个求解域上的近似解。随着单元数目的增加，解的近似程 度将不断改进。如果单元是满足收敛条件的，近似解最后将收效于精确解。 (1)编程思路 由于传感器的形状、大小、功能各异，对它们的分析大多只能依据工程人员的经验进行，运 算量相当大，且重复劳动率高。采用有限元法对传感器进行元素自动分割，按设计者指定的元素 密度分布函数自动生成节点与元素，将传感器分解为许多相关联的小单元，进而对传感器的特性 做出描述，帮助我们来正确评价传感器的综合性能。对传感器的设计也有一定的指导作用。 (2)网格生成技术 本程序的一个主要功能就是可方便地将有限元分析的计算模型进行网格划分，得到离散模型 因此，网格生成的方便程度是有限元分析程序使用性能的重要指标之一。本例中所 的有关数据。 采用的网格生成方法为几何分区生成法。传感器的外形轮廓可按几何特征划分成若干区域，如三 和丽为河北省自然科学墓金资助项目600249号. 陈书旺:讲师;庄庆德:教授;曲国明:副教授。 459 -~~~~~~~~~~..闻 电子测量与仪器学报 __ — 一~一一一一一一一一一2000Itf+1 角形、矩形、扇形等。我们把这些区域规格化为多边形，每个区域当其端点坐标及反映形状的参 数 (如:半径、圆心等)确定时，便被唯一确定，即可对其进行网格生成。 本研究中的网格元素均采用三角形。三角形单元对复杂边界有较强的适应能力，因此很容易 将一个二维域离散成有限个三角形单元，在边界上以若千段直线近似原来的曲线边界. 为防止在分割中出现导致计算发散的钝角三角形，本程序可以自动识别三角形形状。必要时， 可将原设定的区域分割成几个子区域，再对各子区域分别加以处理，以保证分割正确地实行。 (3)均匀分割的程序分析 首先，以平面四边形区域、均匀元素尺寸为例，对程序加以介绍。如图1所示。四边形区域 为圆环形，内、外弧所对应的公共圆心坐 标为 (x,y)，内、外弧的半径分别为r1, r2，弧的起始、终止角度为al,a2，内弧 边上的分割数为k，在直线边上的分割数 为k1。这些参数均由用户预先输入设定。 计算步骤为: ①依据环形区域中内环半径r1及外 环半径r2、圆心坐标((x,Y)扇形起止角度 (a