0.100E+05 0 0.00000 0.00000 0.00000 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ TIME HISTORY $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*DATABASE_BINARY_D3PLOT 0.1000E+03

*DATABASE_BINARY_D3THDT 0.1000E+02

*DATABASE_EXTENT_BINARY 0 0 3 1 0 0 0 0 0 0 4 0 0 0 $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ SECTION DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*SECTION_SOLID 2 1

*SECTION_SOLID_ALE 1 5 $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ MATERIAL DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*MAT_RIGID

2 7.80 2.10 0.280000 0.0 0.0 0.0 1.00 7.00 7.00

*MAT_RIGID

3 7.80 2.10 0.280000 0.0 0.0 0.0 1.00 6.00 7.00

*MAT_RIGID

4 7.80 2.10 0.280000 0.0 0.0 0.0 1.00 6.00 7.00

*MAT_NULL

1 1.0000000 0.0000000 1.00000-8 0.0000000 0.0000000 0.0000000 0.0000000 *EOS_LINEAR_POLYNOMIAL

1 1.00000-6 1.92100-3 0.0000000 0.0000000 0.4000000 0.4000000 0.0000000 0.0000000 0.0000000 $ $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ PARTS DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ *PART

Part 1 for Mat 4 and Elem Type 1 1 1 1 1 0 0 0 $ *PART

Part 2 for Mat 2 and Elem Type 1 2 2 2 0 0 0 0 $ *PART

Part 3 for Mat 3 and Elem Type 1 3 2 3 0 0 0 0 *PART

Part 3 for Mat 3 and Elem Type 1 4 2 4 0 0 0 0 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ RIGID BOUNDRIES $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*CONTACT_SURFACE_TO_SURFACE 1 2 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 *CONTACT_SURFACE_TO_SURFACE 1 3 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 *CONTACT_SURFACE_TO_SURFACE 1 4 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 $

*DEFINE_CURVE

1 0 1.000 1.000 0.000 0.000

0.000000000000E+00 1.000000000000E-04 1.000000000000E+05 1.000000000000E-04 *BOUNDARY_PRESCRIBED_MOTION_RIGID 3 2 0 1 -1.00 0 0.000 5.00e3 *BOUNDARY_PRESCRIBED_MOTION_RIGID 4 2 0 1 -1.00 0 5.00e3 0.000 $

通常我们处理液面晃动采用上述方法，LS-DYNA求解器在下一个版本LS970中在MESHFREE功能上增加了更为稳定和高效的EFG技术，在解决此类问题上将给予我们更大的灵活性和更多的选择。

24复杂几何模型的系列网格划分技术

众所周知，对于有限元分析来说，网格划分是其中最关键的一个步骤，网格划分的好坏直接影响到解算的精度和速度。在ANSYS中，大家知道，网格划分有三个步骤：定义单元属性（包括实常数）、在几何模型上定义网格属性、划分网格。在这里，我们仅对网格划分这个步骤所涉及到的一些问题，尤其是与复杂模型相关的一些问题作简要阐述。 一、 自由网格划分

自由网格划分是自动化程度最高的网格划分技术之一，它在面上（平面、曲面）可以自动生成三角形或四边形网格，在体上自动生成四面体网格。通常情况下，可利用ANSYS的智能尺寸控制技术（SMARTSIZE命令）来自动控制网格的大小

和疏密分布，也可进行人工设置网格的大小（AESIZE、LESIZE、KESIZE、ESIZE等系列命令）并控制疏密分布以及选择分网算法等（MOPT命令）。对于复杂几何模型而言，这种分网方法省时省力，但缺点是单元数量通常会很大，计算效率降低。同时，由于这种方法对于三维复杂模型只能生成四面体单元，为了获得较好的计算精度，建议采用二次四面体单元（92号单元）。如果选用的是六面体单元，则此方法自动将六面体单元退化为阶次一致的四面体单元，因此，最好不要选用线性的六面体单元（没有中间节点，比如45号单元），因为该单元退化后为线性的四面体单元，具有过刚的刚度，计算精度较差；如果选用二次的六面体单元（比如95号单元），由于其是退化形式，节点数与其六面体原型单元一致，只是有多个节点在同一位置而已，因此，可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元，减少每个单元的节点数量，提高求解效率。在有些情况下，必须要用六面体单元的退化形式来进行自由网格划分，比如，在进行混合网格划分（后面详述）时，只有用六面体单元才能形成金字塔过渡单元。对于计算流体力学和考虑集肤效应的电磁场分析而言，自由网格划分中的层网格功能（由LESIZE命令的LAYER1和LAYER2域控制）是非常有用的。

二、 映射网格划分

映射网格划分是对规整模型的一种规整网格划分方法，其原始概念是：对于面，只能是四边形面，网格划分数需在对边上保持一致，形成的单元全部为四边形；对于体，只能是六面体，对应线和面的网格划分数保持一致；形成的单元全部为六面体。在ANSYS中，这些条件有了很大的放宽，包括：

1 面可以是三角形、四边形、或其它任意多边形。对于四边以上的多边形，必须用LCCAT命令将某些边联成一条边，以使得对于网格划分而言，仍然是三角形或四边形；或者用AMAP命令定义3到4个顶点（程序自动将两个顶点之间的所有线段联成一条）来进行映射划分。

2 面上对边的网格划分数可以不同，但有一些限制条件。 3 面上可以形成全三角形的映射网格。

4 体可以是四面体、五面体、六面体或其它任意多面体。对于六面以上的多面体，必须用ACCAT命令将某些面联成一个面，以使得对于网格划分而言，仍然是四、五或六面体。

5 体上对应线和面的网格划分数可以不同，但有一些限制条件。 对于三维复杂几何模型而言，通常的做法是利用ANSYS布尔运算功能，将其切割成一系列四、五或六面体，然后对这些切割好的体进行映射网格划分。当然，这种纯粹的映射划分方式比较烦琐，需要的时间和精力较多。 面的三角形映射网格划分往往可以为体的自由网格划分服务，以使体的自由网格划分满足一些特定的要求，比如：体的某个狭长面的短边方向上要求一定要