无限元理论.PDF

17. 无限元理论 马富银(译) 目录： 17.1、支持的拓扑结构 17.2 、自由度 17.3、无限元域的参数 17.4 、收敛条件 17.5、无限元的收敛判据 17.6、主轴正交性 17.7、几何匹配条件 17.8、模型离散 无限元可以接受无界声场域模型，确保一个自由度的条件。 17.1 支持的拓扑结构 有限元指定的单元维数比实际空间维数(或声学单元)少一维，如表17.1 所示。 表17.1：无限元单元 一维单元 二维单元 轴对称单元 三维单元 1 节点元 × 二节点线单元 √ 二节点线单元 √ 四节点四边形单元 √ 三节点线单元 √ 三节点线单元 √ 八节点四边形单元 √ 三节点三角形单元 √ 六节点三角形单元 √ 图17.1 基于线的用户自定义主节点(实心)后Actran 自动生成内部节点(空心) 的无限元(2 节 点或3 节点)示意图，在所有图中，p 的顺序沿径向增加 图17.2 基于四边形面的用户自定义主节点(实心)后Actran 自动生成内部节点(空心) 的无限 元(4 节点或8 节点)示意图 图17.3 基于三角形面的用户自定义主节点(实心)后Actran 自动生成内部节点(空心) 的无限 元(3 节点或6 节点)示意图 17.2 支持的拓扑结构 一个无限元节点(一个无限元节点对应无限元的一个单元节点)将增加一个自由度，即声 压，这里的无限元节点不包括有限元与无限元之间耦合面的节点。无限元节点总数等于有限 元—无限元耦合面节点数乘以径向插值序数减一。 备注：Actran 用一个潜在变量代替声压，这种选择的目的是需要处理对流波动方程(声 传播与流动参考第23 章) 。 17.3 无限元域的参数 Actran 通过无限元来模拟声传播区域的无边界传播条件，一个无限单元可以通过沿着一 个基面的节点(如图 17.1 到 17.3 所示)和这个基面的无限元参数进行定义。无限元域是一个 描述离散无限声传播区域边界的无限单元的集合，对于传统的声辐射问题，通常定义一个无 限元域，而对于传统的声传播问题，需要定义两个半无限元域(如图17.4)。 一个无限元域由下面的一些参数确定： (1)无限元域的材料参数； (2)插值阶次； (3)参考椭球坐标系(原点和三个主轴) ，支持不连续的多极扩展逼近； (4)可选择的定义一个统一的流动参数。 图17.4 声辐射中的单一无限元域(左)和声传播中的两个半球无限元域(右)示意图 用户需要通过这些参数来确定椭球坐标系，坐标系定义如下： (1)坐标原点( ， 和 三个坐标值) ； xc y c zc (2)坐标轴：由 ， ， 三个半轴长度和 ， ， 三个正交向量定义的三个正交方向。 a b c v v v 1 2 3 参考椭球表面由局部坐标系下的如下公式确定： 2 2 2 x y z 2  2  2 1 (17-1) a b c 根据 ， ， 之间的关系，参考椭球坐标系可以简化为椭球的形