FELAC数据结构介绍.ppt

FELAC数据结构介绍 有限元语言编译程序计算中，需要三种类型的输入数据：有限元网格相关数据、材料参数数据、时间及计算控制数据。 有限元网格相关数据包含了节点坐标、单元编码、规格数、初值、边值等信息，如果以 GID作为前处理，则这些数据都将存储在 GID格式的数据文件 DAT中，通过系统接口函数 gidpre能够自动将 DAT数据读入程序中，这些数据就是 coordinates、id、ubf、unod、element。 材料参数数据是单元计算中与物理方程相关的系数，这些数据存储在表格文件MAT中，本系统提供了MAT文件自动生成功能，用户只需要在PRE文件中给出这些材料参数就行。 控制参数time0文件 TIME0文件内容如下： 3.0 tmax 2.0 dt 1 nstep 1000 itnmax 2.0e-8 tolerance 2.0e-3 dampalfa 2.0e-4 dampbeta 共有 7行，第一行给出计算最大时刻，系统通过该时刻来判断计算是否结束；第二行给出时间增量，隐含的计算开始时刻为 0；第三行给出结果存储计算步间隔，对于时间相关的问题，计算将可能产生大量的结果数据，为了减少存储量和结果可视化的数据量，可以设定每隔多少计算步保存结果，中间计算步的结果将跳过不保存；第四行给出非线性计算允许最大迭代步，超过这个设定的迭代次数，系统将结束迭代，给出结果；第五行给出非线性迭代收敛性容差；第六行和第七行给出固体力学动力计算中采用瑞利阻尼矩阵时的参数（质量矩阵前的系数和刚度矩阵前的系数）。 计算结果图形显示pos文件 有限元计算将把各种物理场的计算结果以文件的形式保存在计算机的磁盘，文件将保留物理场的各个分量，为了告诉可视化程序绘制哪些物理场和分量的图形，用户需要通过一个以 POS为扩展名的文件（以下简称 POS文件给出）。 POS文件要由两段信息组成：第一段信息给出显示多少个计算步（预留参数，无用），第二段信息给出显示（对应于某个物理场）的结果指针名及其显示内容。 第一段内容只有一行，为整型数。 第二段内容紧接第一段，可以有多行，每行的编写规则为： 指针名自由度名 …自由度名整型数 …整型数 这里指针名代表用于存储计算结果的指针数组名，如 unod，自由度名代表结果中所有自由度的名称，后面的整型数代表显示这些自由度的序号。如： unod u v w 1 2 3 就代表把 unod指针中的 u、v、w都显示出来，而 unod u v w 1 2 则表示只把 u、v显示出来。 有限元计算程序根据 POS文件设置，将结果适时写出到磁盘文件中，对于瞬态问题计算，写结果与计算是同步完成的。 Thank you