有限元的自動生成及其应用演示.pptVIP

;一、通用有限元软件的局限性 ;二、自己编程的种种困难;三、开发FEPG的动机与目的;四、由有限元公式自动产生有限元程序;五、软件实现思想;数据转换程序;有限元程序流程框图;有限元程序流程框图;有限元程序流程框图;有限元程序流程框图;自动生成系统由近百个元件程序组成。 元件程序是以文件名为输入输出参数的可独立运行的程序。 由脚本（批命令）文件把元件程序组成一个“程序系统”。 大大降低程序的复杂性，大大提高程序的可读性和再用性。 实现软件自由化（即自由软件）的一种软件设计方法。 ;提供一种用以描述有限元问题的语言，以下称为有限元（描述）语言。有限元语言是一种有层次的脚本语言。有限元语言应包括微分方程表达式，形函数表达式，算法表达式等有限元公式的书写形式。有限元语言应与有限元的专业语言尽可能一致，以便广大科技工作者易于理解与接受。有限元语言要适用于各种领域和各种有限元方法与有限元问题。关键问题是如何描述有限元问题，采用什么样的方式准确表达有限元问题及其求解算法。;用FEPG生成求解二维线弹性平面应变问题的有限元程序;问题的方程描述;虚功方程;需填写的FEPG公式文件;Strn.gcn文件;Strn.gio文件;Disp.vde文件;Disp.vde文件;Disp.vde文件;Strs.vde文件;Strs.vde文件;三维正交各向异性材料计算;三维弹性本构关系;其中;其中;需填写的FEPG公式文件;GCN文件;Gio文件;Disp.vde文件;Disp.vde文件;Disp.vde文件;平面二维水沙运动问题;问题描述;该问题的偏微分方程;该问题的偏微分方程;算子分裂法求解该问题;需填写的FEPG公式文件;GCN文件;if exist stop del stop :1 bft solvsin a copy unod unoda if exist end del end :2 solvsin b if not exist end goto 2 solvsin c solvsin d solvsin e call post if not exist stop goto 1;GIO文件;FLD.VDE文件;FLD.VDE文件;算子分裂法求解该问题;FLC.VDE文件;func $cv fhun1=+{un_j/x_j}*hun+{hun/x_j}*un_j $cv fhun1=+fhun1*dt+hun1 $cv fhvn1=+{un_j/x_j}*hvn+{hvn/x_j}*un_j $cv fhvn1=+fhvn1*dt+hvn1 fhu=+[hu]+[hu_j]*{un/x_j}*dt+[hu_j/x_j]*un*dt +[hu]*{un_j/x_j}*dt+[hu/x_j]*un_j*dt ? fhv=+[hv]+[hu_j]*{vn/x_j}*dt+[hu_j/x_j]*vn*dt +[hv]*{un_j/x_j}*dt+[hv/x_j]*un_j*dt ? stif dist=+[fhu_i;fhu_i] ? load=+[fhu_i]*fhun1_i ? end;算子分裂法求解该问题;SDD.VDE文件;func $c6 vol = 1.0d0 gradhsx = +[hs/x] ? gradhsy = +[hs/y] ? stif $c6 u2v2=dsqrt(u*u+v*v) $c6 ao=rfai*omega $c6 sstar=ek*(u2v2**3.0d0/(eg*hn1*omega))**em dist=+[gradhs_i;gradhs_i]*vol*edt+[hs;hs]*rfai*omeg*rho/hn1/rous load=+[hs]*sstar*rfai*omeg*rho/rous end ;算子分裂法求解该问题;SDC.VDE文件;func fhs=+[hs]+[hs/x_j]*un1_j*dt +[hs]*{un1_j/x_j}*dt ? stif dist=+[fhs;fhs] ? load=+[fhs]*(hs1-rho*(hn-hn1)) ? end ;算子分裂法求解该问题;ZBH.VDE文件;stif $c6 u2v2=dsqrt(u*u+v*v) $c6 ao=rfai*omeg $cv sstar=ek*(u2v2**3.0d0/(eg*h*omeg))**em dist=+[zb;zb]*0.0d0 ? load=+[zb]*ao*(hs/h-sstar) ? end ;结果图;线性正弦稳态涡流问题（TEAM Workshop 问题7） ;问题描述;求解方程及边界条件;需填写的FEPG公式文件;GCN文件;GIO文件;方程的afi.vde