资源描述:
《多场耦合问题的程序生成.ppt》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、多场耦合问题的程序生成GCN文件GIO文件学习目的1.掌握GCN文件的填写格式;2.掌握用GCN功能生成程序的一般过程;3.了解非线性问题线性化的过程。GCN文件的结构aNFE文件名耦合场名&bNFE文件名耦合场名cNFE文件名耦合场名............空一行命令流库中NFE文件命名基本原则NFE算法文件是FEPG系统用来生成单元计算程序(即E元件程序)和后处理计算程序(即U元件程序)描叙算法,系统提供大量供选择的算法。算法是与方程的类型有关的(方程类型可参看关于《PDE文件填写》的pow
2、erpoint文档),库中NFE文件命名基本原则:ell——椭圆par——抛物wave——波动n*——非线性ELL.NFE求解线性椭圆型方程的算法程序。PAR.NFE求解线性抛物型方程的算法程序,时间离散采用Crank-Nicolson格式。PARB.NFE求解线性抛物型方程的算法程序,时间离散采用向后差分格式。WAVE.NFE求解线性双曲型方程的算法程序,时间离散采用Wave速度法格式NEWMARK.NFE求解线性双曲型方程的算法程序,时间离散采用Newmark格式。NELL.NFE求解非线性
3、椭圆型方程的算法程序。NPAR.NFE求解非线性抛物型方程的算法程序,时间离散采用Crank-Nicolson格式。NPARB.NFE求解非线性抛物型方程的算法程序,时间离散采用向后差分格式。NWAVE.NFE求解非线性双曲型方程的算法程序,时间离散采用Wave速度法格式。NNW.NFE求解非线性双曲型方程的算法程序,时间离散采用Newmark格式。STR.NFE已知位移场求应力场的算法程序,采用最小二乘法。也是求解已知场量的梯度的算法程序。命令流命令流段给出耦合场的计算流程。该段由两种命令组成
4、,一种是DOS命令;另一种是本系统的专有命令。较常用的专有命令可参看《FEPG中级教程》关于GCN文件填写的部分。线性稳态问题命令流段startcsolvcSTART元件程序E元件程序SOLV求解器U元件程序命令流段线性稳态问题流程线性动态问题命令流段startcifexiststopdelstop:1bftsolvcifnotexiststopgoto1命令流段线性动态问题流程START元件程序BFT元件程序E元件程序SOLV求解器U元件程序时间计算是否结束否是非线性稳态问题命令流段start
5、cifexistenddelend:1solvcifnotexistendgoto1命令流段非线性稳态问题流程START元件程序E元件程序SOLV求解器U元件程序迭代是否结束否非线性动态问题命令流段startcifexiststopdelstop:1bftifexistenddelend:2solvcifnotexistendgoto2ifnotexiststopgoto1命令流段非线性动态问题流程START元件程序BFT元件程序E元件程序SOLV求解器U元件程序迭代是否结束否时间计算是否结束是
6、否GIO文件的结构a场的PDE文件名b场的PDE文件名c场的PDE文件名............#elemtype单元类型坐标系单元类型包括一维有l2,l3,l2g2,l3g2,l3g3二维有t3,t6,q4,q8,q9,q4g2,q8g3,q9g3三维有w4,w10,c8,c20,c27,c8g2,c20g3,c27g3坐标系包括一维有1dx,1dr,1ds二维有2dxy,2dro,2drz,2drs三维有3dxyz,3droz,3drso如果GIO文件没有第二段对一维问题体单元类型取l2,坐
7、标系取1dx;对二维问题,体单元类型取t3,坐标系取2dxy;对三维,体单元类型取w4,坐标系取3dxyz。关于单元类型和坐标系的说明如何生成有限元程序1.GIOGCN文件名三个例子1.poisson方程的求解2.热力耦合问题3.非线性问题对于如下的piosson方程边界条件如下准备poisson.vde文件dispu,coorx,y,zshap%1%2gaus%3load(-6.0)vectxxyzstifdist=+[u/x_i;u/x_i]end准备ell.gcn文件defiaell&St
8、artsinaSolvsina准备ell.gio文件Poisson#elemtypec83dxyz运行gioell生成全部程序2.按照某个特定问题,建好几何模型并施加边界条件完成poisson方程的定解问题描述:平板长1米,宽0.5米;左端温度为0℃,右端温度为100℃,下端完全固定。求在此条件下的板的温度分布、变形和应力。板的线膨胀系数1.0e-5/℃,弹性模量为1000Mpa,泊松比为0.3,热传导系数为10W/m/℃。不计板的体力和内热源。热固耦合问题问题分析:1.该问题包括了三个物理场,