ansys子结构分析实例解析

ANSYS中的超单元

从8.0版开始，ANSYS中增加了超单元功能，本文通过一些实际例子，探讨了ANSYS 中超单元的具体使用。

1 使用超单元进行静力分析

根据ANSYS帮助文件，使用超单元的过程可以划分为三个阶段(称为Pass)：

(1) 生成超单元模型(Generation Pass)

(2) 使用超单元数据(Use Pass)

(3) 扩展模型(Expansion Pass)

以下摘自htbbzzg邹老师博客，请勿乱传!

下面以一个例子加以说明：

一块板，尺寸为20×40×2，材料为钢，一端固支，另一端承受法向载荷。

首先生成原始模型se_all.db，即按照整个结构进行分析，以便后面与超单元结果进行比较：

首先生成两个矩形，尺寸各为20×2。然后定义单元类型shell63；

定义实常数1为： 2 (板厚度)。

材料性能：弹性模量E=201000；波松比μ=0.3；密度ρ=7.8e-9；

单位为mm-s-N-MPa。

采用边长1划分单元；一端设置位移约束all，另一端所有(21个)节点各承受Z向力5。计算模型如下图：

静力分析的计算结果如下：

为了后面比较的方便，分别给出两个area上的结果：

超单元部分，按照上述步骤操作如下：

(1) 生成超单元

选择后半段作为超单元，前半段作为非超单元(主单元)。

按照ANSYS使用超单元的要求，超单元与非超单元部分的界面节点必须一致(重合)，且最好分别的节点编号也相同，否则需要分别对各节点对建立耦合方程，操作比较麻烦。

实际上，利用ANSYS中提供的mesh200单元，对超单元和非超单元的界面实体，按照同一顺序，先于所有其它实体划分单元，很容易满足界面节点编号相同的要求。对于多级超单元的情况，则还要结合其它操作(如偏移节点号等)以满足这一要求。

对于本例，采用另一办法，即先建立整个模型，然后再划分超单元和非超单元。即：将上述模型分别存为se_1.db (超单元部分)和se_main.db (非超单元部分)两个文件，然后分别处理。

对于se_1.db模型，按照超单元方式进行处理。由于模型及边界条件已建立，只需删除前半段上的划分，结果就是超单元所需的模型。

然后直接进入创建超单元矩阵的操作，首先说明一下创建超单元矩阵的一般步骤：

A进入求解模块：

命令：/Solu

GUI：Main menu -> Solution

B设置分析类型为“子结构或部件模态综合“

命令：ANTYPE

GUI ：Main menu -> Solution -> Analysis Type -> New Analysis

选择Substructuring/CMS (子结构或部件模态综合)

C设置子结构选项

命令：SEOPT

GUI：Main menu -> Solution -> Analysis Type -> Analysis Options

设置内容有：

方程求解方法(波前法–默认，或稀疏矩阵法–推荐)；

要求生成的矩阵，可以是：刚度矩阵；刚度和质量矩阵；刚度、质量和阻尼矩阵(稀疏矩阵法不能生成阻尼矩阵)；

要求输出到output窗口的矩阵(矩阵和载荷矢量或仅输出载荷矢量，默认不输出) ;

扩展超单元的方法(后续操作时，对本超单元扩展时使用的方法)，默认是回代法(backsubstitution method)，需要存储主自由度的三角矩阵(波前法为sebname.tri，稀疏矩阵法为sename.LNxx)；完全求解方法(Full resolve method不存储任何三角矩阵)。

随着问题规模加大，三角矩阵可能很大；完全求解方法不需要三角矩阵，但花费时间较长，且不能重启动子结构分析(对于不同的载荷矢量)。

质量矩阵公式：是否使用集中质量矩阵，推荐用一致质量矩阵。

D选择主自由度

命令：M

GUI：Main menu -> Solution -> Master DOFs -> User Selected -> Define

主自由度是超单元和其它单元之间的界面，必须确保将超单元与其它单元的界面上的所有节点自由度定义为主自由度。假如模型中除了超单元外没有其它单元，也必须定义主自由度。

在将超单元用于动力分析时，主自由度的动力特性应能代表超单元的动力特性。对于在使用超单元时将被定义为约束或施加载荷的节点，应该定义为主自由度(如果它们也位于当前的超单元中)。

对于大挠度问题或使用SETRAN命令的情况，主节点的6个自由度都应作为主自由度。

本例中选择超单元外端面上所有节点的所有自由度作为主自由度。一般情况均推荐选择相应主节点的All自由度。

E在超单元上施加载荷

作为超单元的结构可以施加任意载荷，但要记住：

程序会生成一个载荷矢量，其中包含了施加的所有载荷。对于每个载荷步，在超单元文件中写一个等效的载荷矢量，最多可以写31个载荷矢量。

在超单元结构中可以使用非零位移约束，并作为载荷矢量的一部分。求解后对超单元进行结果扩展时，对于有非零位移的载荷步，其数据库必须匹配，否则需再次施加非零位移约束后再进行结果扩展。

可以在使用超单元时再定义载荷和约束，此时应注意将相应节点选为主自由度节点。

同样也可以在使用超单元时再定义惯性力(加速度或角速度)，但在超单元部分要生成质量矩阵。如果在使用超单元时需要对其进行旋转，推荐这样做，因为超单元的载荷矢量是随着超单元一起转动的。

注意：如果生成质量矩阵，推荐在使用超单元时再施加约束，但要注意将约束节点选为