高级技术分析指南chap4自适应网格划分

1、第四章4.1 何为？ANSYS 程序提供了近似的技术自动估计特定分析类型中因为网格划分带来的误差。（误差估计在“ANSYS 基本过程手册”5.3.5 分析计算误差一节中。）通过这种误差估计，程序可以确定网格是否足够细。如果不够的话，程序将自动细化网格以减少误差。这一自动估计网格划分误差并细化网格的过程就叫做，然后通过一系列的求解过程使得误差低于用户指定的数值（或直到用户指定的最大求解次数）。4.2的先决条件中包含一个预先写好的宏，ADAPT.MAC，完成ANSYS的功能。用户的模型在使用这个宏之前必须满足一些特定的条件（在一些情况下，不满足要求的模型也可以用修正的过程完成这些要求包括：，下面还

2、要）。标准的 ADAPT 过程只适用于单次求解的线性静力结构分析和线性稳态热分析。模型最好应该使用一种材料类型，因为误差计算是根据平均结点应力进行的，在不同材料过渡位置往往不能进行计算。而且单元的能量误差是受材料弹性模量影响的。因此，在两个相邻单元应力连续的情况下，其能量误差也可能由于材料特性不同而不一样。在模型中同样应该避免壳厚突变，这也可能造成在应力平均是发生问题。模型必须使用支持误差计算的单元类型。（见表 41）模型必须是可以划分网格的：即模型中不能有引起网格划分出错的部分。表 4-1可用单元4-1类型单元说明二维结构实体PLANE2 PLANE25 PLANE42 PLANE82 PL

3、ANE832-D 6-Node Triangular Solid Axisymmetric Harmonic Solid2-D 4-Node Isoparametric Solid 2-D 8-Node SolidAxisymmetric Harmonic 8-Node Solid三维结构实体SOLID45 SOLID64 SOLID73 SOLID92 SOLID953-D 8-Node Isoparametric Solid 3-D Anisotropic Solid3-D 8-Node Solid with Roional DOF 3-D 10-Node Tetrahedral Soli

4、d3-D 20-Node Isoparametric Solid三维结构壳体S43S63S93Plastic quadrilateral S Elastic Quadrilateral S 8-Node Isoparametric S二维热实体PLANE35 PLANE752-D 6-Node Triangular Solid Axisymmetric Harmonic Solid4.3 如何使用进行：基本过程的基本过程包括如下步骤：1. 象其他线性静力分析或稳态热分析一样，先进入前处理器（/PREP7 或MaenuPreprosor）。然后指定单元类型，实参和材料特性，要满足上面提到的条件。

5、2. 用实体建模过程建立模型，用可以划分网格的面或体建模。用户不需指定单元大小也不用划分网格，ADAPT 宏会自动划分网格。（如果要同时划分面和体网格，生成 ADAPTMSH.MAC 用户子程序见后。）3. 在 PREP7 中或在 SOLUTION(/SOLU 或 MaenuSolution)中指定分析类型，分析选项，载荷和载荷步选项。在一个载荷步中仅施加实体模型载荷和惯性载荷（加速度，角加速度和角速度）。（通过 ADAPTBC.MAC 用户子程序可以施加有限单元载荷，固连和约束方程。通过 ADAPTSOL.MAC 用户子程序可以加入多个载荷步。这些子程序在本章稍后还要。）如果在 PREP7

6、中，退出前处理器FINISH。（可以在 SOLUTION 或在初始状态下调用 ADAPT 宏）。用下列方法激活自适应求解。d: ADAPTGUI: MaenuSolutionAdaptive Mesh注意，可以在热或结构分析中使用 ADAPT 宏，但不能在一次自适应分析中同时进行这两种不同类型的计算。在的迭代过程中，单元的大小将作调整（在 FACMN 和 FACMX 指定的范围内）以减小或增加单元能量误差，直到误差满足指定的数值（或指定的最大求解次数）为止。6. 当自适应网格计算收敛时，程序自动将单元形状检查打开SHPP,ON。然后返回 SOLUTION 或初始状态，这取决于调用 ADAPT

7、的状态。接下来可以进入 T1 用标准操作进行后处理。4.4 修改基本过程4.4.1 选择自适应性如果用户清楚某个部分网格划分的误差相对影响较小时（如应力水平较低且变化较小），可以将这些区域从中排除，以加快分析速度。同样，用户也许想将接近应力奇异点的部分（如集中载荷）排除掉。选择逻辑操作可以解决这类问题。4-2PLANE55 PLANE77 PLANE782-D 4-Node Isoparametric Solid 2-D 8-Node SolidAxisymmetric Harmonic 8-Node Solid三维热实体SOLID70 SOLID87 SOLID903-D 8-Node Is

8、oparametric Solid 3-D 10-Node Tetrahedral Solid3-D 20-Node Isoparametric Solid三维热壳体S57Plastic Quadrilateral S图 4-1选择自适应能改进有应力集中的模型如果用户选择了一个关键点集，ADAPT 宏仍将包含进所有的这些关键点（即 ADAPT 宏仍将对选择的和未选择的关键点都作网格改动），除非将 ADAPT 命令（MaenuSolutionAdaptive Mesh）中的 KYKPS 设为 1。如果用户选择了一个面或体集，ADAPT 宏将只在选择的区域调整网格大小。此时必须在调用 ADAPT

9、宏之前在 PREP7 中对整个模型进行网格划分。4.4.2 用用户子程序定制 ADAPT 宏标准的 ADAPT 宏并不能满足特定的分析需要。例如，用户可能想同时对面和体进行网格划分，这在标准宏当中是不可以的。对于这种或其他一些类似情况，可以对 ADAPT 宏进行修改使之适用于特定的分析。ANSYS 程序这种方式完成，本身就使得用户可以对其进行相应的修改以适应不同的要求。方便的是，用户不用总是通过修改 ADAPT 代码的方式来定制宏。宏的三个部分可以用用户子程序的方法来修改，这个方法将 ADAPT 宏和用户文件分开，用户可以生成子程序由 ADAPT 宏来调用。这三个部分是：网格划分命令序列，边界

10、条件命令序列和求解命令序列。相应的用户子程序名为 ADAPTMSH.MAC,ADAPTBC.MAC 和 ADAPTSOL.MAC。下面看一下这三个子程序的功能：4.4.2.1 生成用户网格划分子程序（ADAPTMSH.MAC）缺省情况下，如果模型中包含一个或多时，ADAPT 宏将只对体划分网格而不对面进行划分。如果当前选择集中没有体，宏才对面进行划分。如果要同时对面和体进行划分的话，就要生成用户子程序 ADAPTMSH.MAC 来提供相应的操作。在重新划分网格之前要清楚所有这些实体中划分过网格的实体。子程序大致如下：C* Subroutine ADATMSH.MAC-Your name-Job

11、 Name-Dreated4-3TYPE,1 ACLEAR,3,5,2 VCLEAR,ALL AMESH,3,5,2 TYPE,2 VMESH,ALL!指定划分网格的单元类型属性!清除本程序中要重新划分网格的面和体的网格!对面3和5 划分网格（ADAPT 不对其他面划分网格）!改变体单元类型!对所有体划分网格请查阅 TYPE,ACLEAR,VCLEAR,AMESH 和 VMESH 命令得到更详细的解释。建议将 C*行包含在文件中以区别不同的宏。这一行将在任务的输出中出现，就可以确认 ADAPT 宏正确地调用了用户子程序。4.4.2.2 生成用户边界条件子程序（ADAPTBC.MAC）ADAPT

12、 宏在每次求解循环时都要清除并重新划分网格，因此模型的结点和单元也要不断的改变。这就使得所有的有限单元载荷，度固连，约束方程等施加在结点和单元的边界条件都不能使用。如果要包含这些有限单元约束时，就要使用用户子程序，ADAPTBC.MAC。在这个子程序中，可以选择结点然后定义有限单元载荷，度固连和约束方程。下面是一个 ADAPTBC.MAC 的例子：C* Subroutine NSEL,S,LOC,X,0 D,ALL,UX,0 NSEL,S,LOC,Y,0 D,ALL,UX,0 NSEL,ALLADAPTBC.MAC-Your name-Hob Name-D!选择 X 坐标为 0 的结点!令选择

13、的所有结点 UX 为 0!选择 Y 坐标为 0 的结点!令选择的所有结点 UY 为 0!选择所有结点reated4.4.2.3 生成用户求解子程序（ADAPTSOL.MAC） ADAPT 宏中的求解命令集很简单：/SOLU SOLVE FINISH这个缺省令集只能求解单个载荷步。用户可以将其他令集加入用户子程序 ADAPTSOL.MAC 中。4.4.2.4 关于用户子程序的一些其他说明用户可以象生成其他文件一样生成这些子程序。就是说可以用 APDL 语言的*CREATE 命令（UtilityMacroCreate Macro）和 APDL 的*END 命令，也可以用外部的文件编辑器。当 ADA

14、PT 宏调用这些子程序时，先搜索 ANSYS 根目录，再搜索用户根目录，最后是当前目录。因此，要确保其他目录中不包括与所用文件同名的文件。注释行（C*）会在输出文件中出现，可以通过它来检查是否使用了正确的文件。另外，通过在运行 ADAPT 宏之前用/PSEARCH,OFF(UtilityMacroMacro Search Path)指定只搜索 ANSYS 根目录和用户当前目录可以从某种程度上减少文件混用的可能。不管这些子程序在什么位置，它们都能被找到，除非将 ADAPT 命令的 KYMAC 选项设为 1。4.4.3 定制 ADAPT 宏(UADAPT.MAC)有些情况下用户需要修改 ADAPT

15、 宏但不能通过单独的用户子程序的方式，4-4那么就需要直接修改 ADAPT 宏的主体。但是，由于某些原因，不直接对ADAPT 宏进行修改。（例如，别的用户和你同时使用一个，在调用 ADAPT 宏时会发现宏被修改了！）因此，在 ANSYS 安装中支持一个宏的拷贝文件 UADAPT.MAC，便于用户修改。如果对 UADAPT.MAC 文件进行了修改，建议对修改后的文件取一个新的文件名。然后在调用时输入这个文件名。要知道的是，如果新文件名是一个“unknownd”，ANSYS 将搜索目录，然后是登录的目录，最后是工作目录，直到找到这个宏为止。如果修改的宏只能为一个用户使用，那么置应在用户登录目录的层

16、次之下（不能等于或高于这个目录层次）。这样，的位的低层次的文件可以通过*USE 命令（Utility来调用。macroExecute Data Block）4.5的几点说明下面的建议可能有助于的使用：不需指定初始网格大小，但指定大小可能有利于自适应收敛。如果用户指定了关键点网格大小，ADAPT 宏在第一次循环时使用这个值，然后在随后的循环中进行调整。用下列命令指定单元大小：d: KESIZEGUI: MaenuPreprosor-Meshing-Size Cntrls-Keypos-All KPsMaenuPreprosor-Meshing-Size Cntrls-Keypos-Picked

17、KPs如果定义了线分段数或大小比例，ADAPT 宏将在每次循环中都使用这个数值而不作改变。如果没有定义的网格份数，在初始网格划分时将使用缺省的网格大小参见 SMRTSIZE 和 DESIZE 命令。用下列方法指定线分段数或大小比例：d: LESIZEGUI: MaenuPreprosor-Meshing-Size Cntrls-Lines-All LinesMaenuPreprosor-Meshing-Size Cntrls-Lines-Picked Lines网格划分适用于 2-D 实体和 3-D 壳单元。但面的划分效果不明显。划分效果要好得多。网格划分适用于 3-D 实体。对体进行划分比总

18、体上说，在自适应网格中有中间结点的单元比线性单元要好。不要用集中载荷或尖角等引起奇异性的结构，因为此时 ADAPT 在这些奇异点处能量值将不收敛。如果模型中有集中载荷时，将其用施加在一个小面上的压力等效。（或通过选择将奇异部分排除在之外。）在许多情况下，用一系列相对小的区域替代少数几个大的区域将得到更好的网格划分。如果最大响应位置已知或事先可以推测，就在附近放置一个关键点。如果是在交互方式下运行 ADAPT，而 ANSYS 在没有提示出错信息时突然退出，可以在 Jobmame.T 文件中查看部分以确定出错原因。同样，在批处理方式下运行 ADAPT 时，可以看 Jobname.T 确定出错原因。如果模型中有些区域有过度的时，在网格划分中就会出错。在这种情况4-5下 ， 用KESIZE命 令 (MainPreprosor-Meshing-SizeCntrls-Keypos-Picked KPs)中 SIZE 域指定区域附近关键点的最大单元长度。同时，ADAPT 命令中的 FACMX 将设为 1，大小增加。过度部分单元应当结果文件（Jobname.RST 或 Jobnam