并行有限元系统介绍.ppt

用FEPG生成求解二维线弹性平面应变问题的有限元程序 不连续岩体拱坝系统的动力过程数值模拟 北京大学地球物理系 北京大学力学与工程系 中国水利水电科学研究院 地下厂房洞室群的三维围岩稳定分析 —— FEPG系统在水电工程中的应用 五、PFEPG实例 More than 97,600 hexahedral element partitioned to 8 sub-domains 1.并行FEPG 弹性问题在并行机计算结果 (32 子域近280,000 自由度) * 并行有限元程序生成系统 （PFEPG 5.2） 注：本材料取自飞箭公司有关资料 1983年开始研制 1993年研制成功 1994年获得中科院科技进步二等奖 1995年国家科技进步二等奖 1996年在数学家杨乐的建议下开始商品化 1999年成立北京飞箭软件公司 至今已推出FEPG5.0版本，国内用户单位100多家 软件发明人梁国平研究员 FEPG的开发历史 一、FEPG简介 一、通用有限元软件的局限性 通用软件并不通用，只能求解很少的一部分有限元问题。 需要有一支训练有素、组织严密、长期稳定的有限元软件队伍，不断维护和改进其通用有限元程序。 只能求解今天已有的微分方程，不能求解明天的有限元问题, 只能采用今天已出现的算法，不可能采用明天才有的算法。 不公开源代码，对用户来说程序是不可改动的，极大地限制了用户的创造性。 二、自己编程的种种困难 有限元程序复杂，代码庞大。 编程者缺乏软件专业训练，代码可读性差，再用性差，程序难以维护，程序缺乏严密组织与发展。 大多数工程师、教授与研究人员，他们主要是研究有限元方法或把有限元方法应用于工程设计与科学研究，而不是致力于有限元软件的开发，因此他们不可能花很多时间编写程序。 三、开发FEPG的动机与目的 许多专家希望有自己的有限元程序，但不愿花许多时间编程。 FEPG应该不受领域的限制，各种专业的有限元问题都可采用本系统求解。 本系统不可能象通用有限元软件那样，通过不断增加单元模块和算法模块来逐步扩大其求解范围。 不管什么样的有限元问题，不管采用什么样的算法，都能采用本系统求解。 本系统是一种公开源代码的自由软件，任何非本系统开发者都可维护和修改。 程序生成——实例之 问题的方程描述 平衡方程 几何方程 本构方程 虚功方程 虚功方程弱形式 将本构关系代入上式即可得以位移为基本未知量的虚功方程 需填写的FEPG公式文件 GIO文件 VDE文件 GCN文件 Strn.gcn文件 defi a ell b str a startsin a solvsin a stress b Strn.gio文件 disp strs #elemtype t3 2dxy Disp.vde文件 disp u v coor x y func ex ey exy shap %1 %2 gaus %3 mate pe pv fu fv 1.0e10;0.3;0.0;-100.0; $c6 fact = pe/(1.+pv)/(1.-2.*pv) vect u u v vect ed ex ey vect fs fu fv matr de 2 2 (1.-pv) pv pv (1.-pv) Disp.vde文件 func ex=+[u/x] ey=+[v/y] exy=+[u/y]+[v/x] stif dist=+[ed_i;ed_j]*de_i_j*fact+[exy;exy]*(0.5-pv)*fact load=+[u_i]*fs_i end Disp.vde文件 func ex=+[u/x] ey=+[v/y] exy=+[u/y]+[v/x] stif dist=+[ed_i;ed_j]*de_i_j*fact+[exy;exy]*(0.5-pv)*fact load=+[u_i]*fs_i end Strs.vde文件 disp sa,sb,sc, coef u v coor x y shap %1 %2 gaus %3 mass %1 mate pe pv fu fv 1.0e10;0.3;0.0;-100.0; $c6 fact = pe/(1.+pv)/(1.-2.*pv) vect u u v vect ed ex ey vect fs fsa fsb vect ss sa sb sc vect ff fsa fsb fsc matr de 2 2 (1.-pv) pv pv (1.-pv) Strs.vde文件 stif $cv ex=+{u/x} $cv ey=+{v/y} $cv fs_i=+de_i_j*ed_j*fact $cv func=+{u/y}+{