应用FEPG求解物理问题初步.ppt

应用FEPG求解物理问题初步 本讲的内容及目标 上讲回顾 FEPG提供了一套有限元语言 结合生成系统，FEPG成为了一个编制有限元程序的平台，用户可以得到全部用于计算的Fortran源代码 FEPG提供了多个学科领域的公式库 有限元分析的基本流程 FEPG的五个元件程序 START BFT ---- Boundary Force Time E ---- Equation, Element SOLV ---- SOLVE U START元件程序的功能 Start 给出每个节点的各个自由度与将来要形成的代数方程组的变量（即方程号）的对应关系。也就是说，哪个节点的那个自由度对应于方程组的哪个变量。 给出解的初值 BFT元件程序的功能 BFT 给出每一时刻解的边值，即指定位移和指定载荷 步长累加，更新计算时间 生成保存计算结果的命令文件post.bat E 元件程序的功能 E 计算单元刚度矩阵、单元质量矩阵、单元载荷向量。并把他们由节点各自由度表示转换成由代数方程组的变量表示 处理边界约束条件，形成代数方程组的右端项 U元件程序的功能 U 把求解器求出的位移转化为各节点自由度的位移 其它后处理计算 五个元件程序 应用有限元方法求解问题 计算程序的内容取决于待解问题及所采用的有限元方法 有限元问题可以是静态的，也可以是动态的 有限元问题可以是线性的，也可以是非线性的 有限元问题求解流程图:线性稳态问题 有限元问题求解流程图:非线性稳态问题 有限元问题求解流程图:线性瞬态问题 有限元问题求解流程图:非线性瞬态问题 单元计算是有限元计算的基本任务 单元计算或单元刚度矩阵的生成是有限元离散的直接体现 单元计算的代码一般不多，但是保证其正确、有效是不容易的 PDE类型文件 描述待解物理问题的方程表达式 这个表达式基于虚位移原理 共包括五类：PDE VDE CDE FDE FBC 要描述什么？ 方程表达式 待解未知量 坐标系及坐标变量 方程中的参数、系数 采用什么类型的单元 矩阵部分 向量部分 方程所能展现的所有内容 未知量名称：u 坐标名称：x, y 几个系数：ek eq ec 形函数：q 4 数值积分方法：q 刚度矩阵: ?? 质量矩阵:?? 右端向量:?? PDE文件编写 PDE文件编写 举例说明：瞬态热传导 PDE文件编写 PDE文件书写 PDE文件的基本段落 DEFI : definition, 定义 ，声明 FUNC: function， 定义函数 STIF: stiff matrix ， 刚度矩阵 LOAD : load ，载荷向量 PDE文件编写 PDE文件－－DEFI信息段 PDE文件－－DEFI信息段 PDE文件－－DEFI信息段 PDE文件－－FUNC信息段 PDE文件－－FUNC信息段 PDE文件－－FUNC信息段 PDE文件－－FUNC信息段 PDE文件－－STIF信息段 PDE文件－－STIF信息段 PDE文件－－STIF信息段 PDE文件－－LOAD信息段 给出虚功方程弱解表达式中对应的荷载项表达式，用以计算单元荷载向量。 LOAD = [未知函数或其导数或自定义函 数]*系数表达式 … … [未知函数或其导数或自定义函数]*系数表达式 PDE文件－－LOAD信息段 PDE文件－－END信息段 PDE文件以空行后跟一个 END 表示文件的结束。 PDE文件－－附加fortran源程序 PDE文件允许用户在DEFI、FUNC、STIF这三个信息段中使用$c6、$c0、$cv、$cc插入fortran源程序和FEPG语句。 此外，PDE中用到的fortran子程序（subroutine）和函数（function），系统还提供了两种附加方法： 在PDE文件的结束字END之后空一行填写关键字FORT，然后附加这些fortran子程序和函数。 把这些fortran源程序放在另外一个文件中，单独编译，最后与PDE文件产生的单元子程序一起连接到单元计算程序中（即E程序），方法见GIO文件书写。 FBC文件编写 PDE文件是用来书写虚功方程弱形式的体积分项，FBC文件是用来书写虚功方程弱形式的边界积分项（相应于第二、第三类边界条件）。给出计算边界条件的单元刚度矩阵、单元质量矩阵、单元阻尼矩阵、单元荷载向量的计算表达式。 边界单元刚度矩阵 STIF DIST = 单元刚度矩阵（分布矩阵）计算表达式 FBC文件编写 边界单元质量矩阵（以下三选一） MASS DIST = 单元阻尼矩阵（分布矩阵）计算表达式 或者