第四章子结构
什么是子结构？
子结构就是将一组单元用矩阵凝聚为一个单元的过程。这个单一的矩阵单元称为超单元。在ANSYS分析中，超单元可以象其他单元类型一样使用。唯一的区别就是必须先进行结构生成分析以生成超单元。子结构可以在ANSYS/Mutiphysics，ANSYS/Mechanical和ANSYS/Structural中使用。
使用子结构主要是为了节省机时，并且允许在比较有限的计算机设备资源的基础上求解超大规模的问题。原因之一如a)非线性分析和带有大量重复几何结构的分析。在非线性分析中，可以将模型线性部分作成子结构，这样这部分的单元矩阵就不用在非线性迭代过程中重复计算。在有重复几何结构的模型中（如有四条腿的桌子），可以对于重复的部分生成超单元，然后将它拷贝到不同的位置，这样做可以节省大量的机时。
子结构还用于模型有大转动的情况下。对于这些模型，ANSYS假定每个结构都是围绕其质心转动的。在三维情况下，子结构有三个转动自由度和三个平动自由度。在大转动模型中，用户在使用部分之前无须对子结构施加约束，因为每个子结构都是作为一个单元进行处理，是允许刚体位移的。
另外一个原因b)一个问题就波前大小和需用磁盘空间来说相对于一个计算机系统太庞大了。这样，用户可以通过子结构将问题分块进行分析，每一块对于计算机系统来说都是可以计算的。
如何使用子结构
子结构分析有以下三个步骤：
●生成部分
●使用部分
●扩展部分
生成部分就是将普通的有限元单元凝聚为一个超单元。凝聚是通过定义一组主自由度来实现的。主自由度用于定义超单元与模型中其他单元的边界，提取模型的动力学特性。图4-1是一个板状构件用接触单元分析的示意。由于接触单元需要迭代计算，将板状构件形成子结构将显著地节省机时。本例中，主自由度是板与接触单元相连的自由度。
图4-1 子结构使用示例
使用部分就是将超单元与模型整体相连进行分析的部分。整个模型可以是一个超单元，也可以象上例一样是超单元与非超单元相连的。使用部分的计算只
是超单元的凝聚（自由度计算仅限于主自由度）和非超单元的全部计算。
扩展部分就是从凝聚计算结果开始计算整个超单元中所有的自由度。如果在使用部分有多个超单元，那么每个超单元都需要有单独的扩展过程。
图4-2示出了整个子结构分析的数据流向和所用的文件。三个步骤的详细解释见以后的叙述。
图4-2 典型子结构分析中的数据流向
生成部分：生成超单元
本部分主要有两步：
1．建立模型。
2．施加边界条件，生成超单元矩阵。
第一步：建立模型
在这一步中，指定文件名和分析名称，用PREP7定义单元类型，单元实参，材料特性和模型几何结构。这些任务在ANSYS绝大多数分析中都是通用的，在ANSYS Basic Analysis Procedures Guide有所叙述。在生成部分，需要记住以下几点：
文件名——在子结构分析中很有用处。有效地使用文件名，在三部分分析中可以省略很多文件处理操作。
用以下方法指定文件名：
Command: /FILENAME
GUI: Utility Menu>File>Change Jobname
如：/FILENAME,GEN
将生成过程中所有文件名都定义为GEN。缺省的文件名是FILE（或file）或在进入ANSYS 后定义的任意文件名。
单元类型——ANSYS提供的绝大多数单元都可以用来生成超单元。唯一的限制是单元必须是线性的。如果生成超单元时有双线性单元的话，ANSYS将自动作为线性单元处理。
注意：在直接耦合中带载荷向量的耦合单元是不能做子结构分析的。可以用同种形状的单元来替代。细节参看ANSYS Coupled-Field Analysis Guide。
材料特性——定义所有必须的材料特性。例如，如果生成质量矩阵，就必须定义密度或其他形式的质量；如果要生成热传导矩阵，就要定义比热。同样，超单元是线性的，非线性材料将被忽略。
模型生成——在生成部分，主要生成模型的超单元部分。非超单元部分是在以后的使用部分生成的。但是，在建模的开始就需要对模型的两个部分有所规划，主要是确定超单元部分和非超单元部分如何连接。为了保证连接正确，应该保证接触部分结点号一致。（其他可以方便用户的方法在本章“使用部分”一节还有介绍。）
要生成整体模型应该这样做：将模型存储在数据库文件中，选择子结构部分进行生成计算。在以后的使用部分，RESUME（Utility Menu>File>Resume from）数据库文件，不选（unselect）子结构，用超单元矩阵代替。
第二步：施加边界条件，生成超单元矩阵。
生成部分的结果包含超单元矩阵。象其他分析一样，用户要定义分析类型和分析设置，施加边界条件，定义载荷步，开始计算。如何完成这些工作见下面的叙述：
1．进入求解器：
Command: /SOLU
GUI: Main Menu>Solution
2．定义分析类型和分析设置：
分析类型——选择生成超单元使用下列方法：
Command: ANTYPE
GUI: Main Menu>Solution>-Analysis Type-New Analysis
新的分析或重启动——如果是开始一个新的分析时，只要指定分析类型（如上所述）即可。如果是重启动计算，必须在ANTYPE命令中设定STATUS=REST （Main Menu>Solution>-Analysis Type-Restart）。如果要另外施加载荷时，可以用重启动。（重启动时，初始运算后的Jobname.EMAT，Jobname.ESA V和Jobname.DB文件要存在。）
超单元矩阵文件名——指定超单元矩阵文件名（Sename）。程序将自动添加后缀.SUB，因此完整的文件名是Sename.SUB。缺省是使用工作文件名[/FILENAME]定义超单元矩阵文件名，可以使用以下命令：
Command: SEOPT
GUI: Main Menu>Solution>Analysis Options
要生成的矩阵——可以指定仅生成刚度矩阵（或传导矩阵，电磁系数矩阵）；生成刚度和质量矩阵（或热传导矩阵等）；生成刚度，质量和阻尼矩阵。质量矩阵用于结构动力学分析和在使用部分有惯性载荷的情况下。在热分析中，只有瞬态热分析才用到热传导矩阵。对于其他分析和阻尼矩阵也大同小异。用SEOPT 命令或其GUI路径来定义。
输出矩阵——这个选项允许输出超单元矩阵。可以指定输出矩阵和载荷向量，也可以只输出载荷向量。缺省值是不输出任何矩阵。要输出矩阵，用SEOPT 命令或其相应的GUI路径。
质量矩阵形成——只在想生成质量矩阵时使用。用户可以选择缺省生成（取决于所用单元类型）或集中质量近似。对于绝大多数情况，推荐使用缺省生成的方式。但是，在极薄构件的分析中，如细长杆或极薄壳体，集中质量近似将得到
