pFEPG学习提纲.docVIP

pFEPG学习提纲 有限元方法基础部分 目前，应用于求解偏微分方程的数值方法主要有三种，即有限元方法，有限体积方法和差分法。经过几十年的发展，有限元方法的基础理论和方法已经比较成熟，已成为当今工程技术领域中应用最为广泛，成效最为显著的数值分析方法。 在工程或物理问题的数学模型（基本变量、基本方程、求解域和边界条件等）确定以后，有限元方法作为对其进行分析的数值计算方法的要点可归纳如下 将一个表示结构或连续体的求解区域离散为若干个子域（单元），并通过他们边界上的节点相互联结成为组合体。 用每个单元内所假设的近似函数来分片地表示全求解区域内待求的未知变量。而每个单元内的近似函数由未知场函数（或其导数）在单元个各界点上的数值和其对应的插值函数连表达。这些插值函数通常采用分片多项式 通过和原问题数学模型等效的变分原理或虚功原理，建立求解未知量的代数方程组或常微分方程组。 pFEPG系统基本思想和架构部分 pFEPG系统的目标是用于输入微分方程表达式和算法之后，就可以得到所有有限元计算的程序源代码。 pFEPG系统采用一种语言三种技术开发而成。这一种语言就是有限元语言，三种技术分别是自动生成技术，元件化技术和公式库技术。 pFEPG系统目前有两个系列，即网络版和并行版。 pFEPG系统的有限元语言部分 pFEPG系统提供的有限元语言包括两个部分，即微分方程表达式的描述语言和算法描述语言。基本功能如下： 描述偏微分方程表达式 用一个文件描述微分方程的虚功方程 描述有限元算法 如何把一个非线性的偏微分线性化 把依赖时间问题对时间离散 多物理场耦合的描述 计算流程和迭代的控制 应用有限元语言描述微分方程表达式将会形成一个PDE类文件 pFEPG的微分方程表达式考虑的是偏微分方程基于虚功原理的弱解形式。 pFEPG采用有限元语言编写文件来描述这个虚功方程。基本的文件类型有PDE, VDE, CDE, FDE, FBC, GES. 单元子程序自动生成系统根据这些文件自动生成计算单元刚度矩阵，单元阻尼矩阵和单元荷载向量等的单元子程序。 PDE文件用来描述虚功方程中的体积分部分的表达式。VDE, CDE, FDE是PDE文件另外的简化形式。 FBC文件用来描述虚功方程中的边界积分部分的表达式。其语法与PDE文件基本相同，但坐标分量少一维。 PDE文件一般包含六个段落，分别是DEFI, MASS, DAMP, STIF, LOAD, END. MASS, DAMP段落不是必需的， 所描述的内容会写在DEFI段里，用MASS和DAMP行表示。 应用有限元语言描述有限元算法将形成两类文件 单物理场的算法描述文件：NFE文件 对于线性问题，算法文件的作用是组集刚度矩阵和载荷向量等；对方程组的解进行简单的运算或处理。线性问题的NFE文件应该包含下面的信息 - 给出如何由单元的刚度矩阵，质量矩阵等计算线性代数方程组矩阵 - 给出如何计算方程组的右端项 - 给出后处理的一些简单计算 - 读写文件 - 给出VECTOR和VAR定义语句 - 嵌入FORTRAN源程序行 非线性问题的算法文件的作用是：读入上一个迭代时刻未知函数的值，并将其传入单元子程序中；组集刚度矩阵和载荷向量等；对方程组的解进行处理，为下一个迭代步准备数据。非线性问题的算法文件应该包含的信息： - 给出单元计算所需的前一迭代时刻的未知函数名及其先后次序 - 给出如何由单元的刚度矩阵，质量矩阵等计算线性代数方程组矩阵 - 给出如何计算方程组的右端项 - 需要从文件读入前一时刻的未知函数的值 - 属于不同微分方程组的未知函数可能对应同一个节点有不同数目的自由度，因此要分别说明 - 给出存放计算结果的节点位移名 - 给出后处理的一些简单计算 - 读写文件 - 给出VECTOR和VAR定义语句 - 嵌入FORTRAN源程序行 多物理场耦合的算法描述文件：GIO, GCN文件 PDE文件和NFE文件的COEF功能的联合使用能够实现多物理场耦合计算。但缺点是需要用户自己编写复杂的NFE文件 通过GCN文件实现多物理场耦合。其优点是用户可以直接使用NFE算法库 GIO文件的结构和语法 待解微分方程声明 单元和积分方法的选取 坐标系类型的声明 GCN文件的结构和语法 第一段信息表明微分方程和算法的对应关系以及物理场之间的耦合关系 第二段信息给出计算流程。同时也是给出了元件程序的组装方式 pFEPG系统的元件程序部分 系统采用元件化程序设计方法，把有限元计算程序分解为若干个元件程序。每个元件程序都是一个完整的FORTRAN程序，单独进行编译、连接与运行，它们之间的通信完全通过磁盘文件，由元件程序组成的程序系统是通过由命令流自动生成的批命令方式实现的，即把所要执行的元件程序按照一定的顺序以操作系统所能接受的命令流方式写