FEPG系统简介.docVIP

FEPG系统简介 系统主要思想和做法 计算机发展至今，在带动社会不断发展和变化的同时，其本身也在不断的发展与进步，从机器代码到高级语言（FORTRAN/C）再到面向对象的程序设计方法的发展，已经引领社会发生了巨大的变革，目前人们都在期待着更为方便的新一代语言的出现。针对有限元研究领域，FEPG为各专业的研究人员创造了一门独特的有限元语言，为各学科领域的数值计算提供了极大的方便，这种语言可以节省90%以上的软件编写代码量。例如针对一个三维稳态热传导问题，采用FEPG提供的有限元语言，用户只需填写以下VDE/PDE、GCN、GIO等几个简单的文件即可生成全部生成有限元计算源代码。 稳态热传导控制方程： 假设求解边界条件为： 按照有限元方法，转化为虚功方程的弱形式： （*） 对应的偏微分方程PDE文件书写如下： disp u u代表未知函数（温度）名 coor x y z 给出总体坐标系下的坐标变量X，Y，Z shap c 8 c表示六面体单元，8表示八节点单元， 也可相应替换为w 4 （四面体四节点单元）等 gaus 2 给出每个方向上积分点的个数 \ mate ek q ek,eq为材料参数名，后面可给定相应值 空一行 stif 给出刚度矩阵 dist=+[u/x;u/x]*ek +[u/y;u/y]*ek +[u/z;u/z]*ek 对应于方程（*）中的左端项 空一行 load=+[u]*q 给出载荷向量，对应于方程（*）中的右端项 空一行 end 结束符 计算流程由以下几行语言规定： Defi 定义物理场 a ell a场采用椭圆方程求解 空一行 STARTC a 初始化a场 SOLVC a 求解a场 连接GCN和PDE的GIO文件： PDE文件名 # elemtype c 8 单元类型 3dxyz 坐标系 以上为采用FEPG独特的有限元语言来求解一个三维稳态热传导问题所需要的文件，它的书写形式接近于通常习惯的写法，采用有限元语言描述有限元方程的代码量远远低于直接采用FRORTRAN/C编程，因此可以把我们从庞大复杂的编程劳动中解脱出来，将更多精力投入到寻找合理的物理模型与更准确的算法中。 为什么会有如此效果？原因就在于FEPG采用了元件化的程序设计方法和程序自动生成技术，有限元程序的元件化，使用户在保证程序完整性的同时可以根据需要更改核心的算法和公理论式，大大减少了工作量，又保证了程序的灵活性和适用性。 1、元件化思想 有限元程序总体可分为三个组成部分：前处理；有限元计算；后处理。 一般有限元程序框图 FEPG系统采用元件化程序设计方法，把庞大的有限元计算程序分解为若干个元件程序。每个元件程序都是一个完整的FORTRAN程序，可以单独进行编译、连接与运行，它们之间的通信完全通过磁盘文件。由元件程序组成的程序系统是通过命令流产生的批命令方式实现的，即把所要执行的元件程序按照一定的顺序以操作系统所能接受的命令流方式写在一个命令流文件中，由FEPG系统自动生成相应的批命令文件，然后运行这个批命令文件，由操作系统解释执行。 有限元计算程序主要由START、BFT、SOLV和E、U五个元件程序组成，这五个元件程序的功能如下： START程序：该程序给出每个节点的各个自由度与将来要形成的代数方程组的变量（即方程号）的对应关系（即哪个节点的哪个自由度将要对应于方程组的哪个变量）以及解的初值。 BFT元件程序：此元件程序的主要功能是给出每一时刻解的边值，即指定节点位移和载荷，以及对时间的更新和生成保存计算结果的批命令文件“post.bat”。 E元件程序：E元件程序用于计算单元刚度、质量