多场耦合问题的程序生成.ppt

《多场耦合问题的程序生成.ppt》由会员分享，可在线阅读，更多相关《多场耦合问题的程序生成.ppt（42页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、多场耦合问题的程序生成GCN文件GIO文件学习目的GCN文件的填写格式；GCN功能生成程序的一般过程；3.了解非线性问题线性化的过程。GCN文件的结构a NFE文件名 耦合场名&b NFE文件名 耦合场名c NFE文件名 耦合场名.空一行命令流库中NFE文件命名基本原则NFE算法文件是FEPG系统用来生成单元计算程序(即E元件程序)和后处理计算程序(即U元件程序)描叙算法，系统提供大量供选择的算法。算法是与方程的类型有关的（方程类型可参看关于PDE文件填写的powerpoint文档），库中NFE文件命名基本原则：ell椭圆par抛物wave波动n*非线性nELL.NFE 求解线性椭圆型方程的算

2、法程序。求解线性椭圆型方程的算法程序。nPAR.NFE 求解线性抛物型方程的算法程序，时间求解线性抛物型方程的算法程序，时间 离散采用离散采用Crank-Nicolson格式。格式。nPARB.NFE 求解线性抛物型方程的算法程序，时间求解线性抛物型方程的算法程序，时间离散采用向后差分格式。离散采用向后差分格式。nWAVE.NFE 求解线性双曲型方程的算法程序，时间求解线性双曲型方程的算法程序，时间离散采用离散采用Wave速度法格式速度法格式 nNEWMARK.NFE 求解线性双曲型方程的算法程序，求解线性双曲型方程的算法程序，时间离散采用时间离散采用Newmark格式。格式。nNELL.NF

3、E 求解非线性椭圆型方程的算法程序。求解非线性椭圆型方程的算法程序。nNPAR.NFE 求解非线性抛物型方程的算法程序求解非线性抛物型方程的算法程序,时间时间离散采用离散采用Crank-Nicolson格式。格式。nNPARB.NFE 求解非线性抛物型方程的算法程序，求解非线性抛物型方程的算法程序，时间离散采用向后差分格式。时间离散采用向后差分格式。nNWAVE.NFE 求解非线性双曲型方程的算法程序，求解非线性双曲型方程的算法程序，时间离散采用时间离散采用Wave速度法格式。速度法格式。nNNW.NFE 求解非线性双曲型方程的算法程序，时求解非线性双曲型方程的算法程序，时间离散采用间离散采用

4、Newmark格式。格式。nSTR.NFE 已知位移场求应力场的算法程序，采用已知位移场求应力场的算法程序，采用最小二乘法。也是求解已知场量的梯度的算法程序。最小二乘法。也是求解已知场量的梯度的算法程序。命令流命令流段给出耦合场的计算流程。该段由两种命令组成，一种是 DOS命令；另一种是本系统的专有命令。较常用的专有命令可参看FEPG中级教程关于GCN文件填写的部分。线性稳态问题命令流段startcsolvcSTART元件程序E元件程序SOLV求解器U元件程序命令流段线性稳态问题流程线性动态问题命令流段startcif exist stop del stop:1bftsolvcif not e

5、xist stop goto 1命令流段线性动态问题流程START元件程序BFT元件程序E元件程序SOLV求解器U元件程序时间计算是否结束否 是非线性稳态问题命令流段startcif exist end del end:1solvcif not exist end goto 1命令流段非线性稳态问题流程START元件程序E元件程序SOLV求解器U元件程序迭代是否结束否非线性动态问题命令流段startcif exist stop del stop:1bftif exist end del end:2solvcif not exist end goto 2if not exist stop got

6、o 1命令流段非线性动态问题流程START元件程序BFT元件程序E元件程序SOLV求解器U元件程序迭代是否结束否时间计算是否结束 是否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,3drson如果 GIO 文件没有

7、第二段对一维问题体单元类型取 l2，坐标系取 1dx；对二维问题，体单元类型取 t3，坐标系取 2dxy；对三维，体单元类型取 w4，坐标系取 3dxyz。关于单元类型和坐标系的说明如何生成有限元程序1.GIO GCN文件名三个例子1.poisson方程的求解2.热力耦合问题3.非线性问题对于如下的对于如下的piosson方程方程边界条件如下准备文件准备文件disp u,coor x,y,zshap%1%2gaus%3load (-6.0)vect x x y zstifdist=+u/x_i;u/x_iend准备文件准备文件defia ell&Startsin aSolvsin a 准备文件

8、准备文件Poisson#elemtype c83dxyz1.运行运行gio ell生成全部程序生成全部程序2.2.按按照照某某个个特特定定问问题题，建建好好几几何何模模型型并并施施加加边边界界条条件件完成完成poisson方程的定解方程的定解问题描述：e-5/,弹性模量为1000Mpa,泊松比为0.3，热传导系数为10W/m/。不计板的体力和内热源。热固耦合问题问题分析：1.该问题包括了三个物理场，温度场、位移场、应力场。需要准备描述三个物理场场方程的VDE或PDE文件。GCN文件。温度场不考虑固体变形和应力分布，但由于材料的热胀冷缩会引起材料变形和热应力。GCN文件和PDE/VDE文件联接关

9、系的GIO文件热固耦合问题热传导方程二维直角坐标下的平面热应力描述方程平衡方程几何方程本构方程disp ucoor x yshap%1%2gaus%3mate ek ec q 1.0e1;0.1;0.0;load qvect x x ystifdist=+u/x_i;u/x_i*ekend温度场文件二维直角坐标下的平面热应力描述方程平衡方程几何方程本构方程disp u,vcoor x,yfunc ex ey exyshap%1%2gaus%3coef tnmate pe pv alfa fx fy 1.0e9;0.3;1.0e-5;0.0;0.0;$c6 fact=pe/(1.+pv)/(1.

10、-pv)$c6 ft=(1.+pv)*fact*alfavect em ex eyvect u u vvect fd fx fyvect fe ft*tn ft*tnmatr de 2 21 pvpv 1 变形场文件FUNCex=+u/x ey=+v/y exy=+u/y+v/x stifdist=+em_i;em_j*de_i_j*fact+exy;exy*fact*(1.-pv)/2 load=+u_i*fd_i+em_i*fe_i end 变形场VDE文件应力场描述方程已知变形和温度分布后，代入几何方程求应变，再通过本构方程求出热应力几何方程本构方程disp sx,sy,sxycoor

11、x,ycoef u,v,tnshap%1%2gaus%3mate pe pv alfa fx fy 1.0e9;0.3;1.0e-5;0.0;0.0;$c6 fact=pe/(1.+pv)/(1.-pv)vect fs sx sy sxyvect ff fsx fsy fsxyvect em ex eyvect fm fsx fsyvect fe ftt fttmatr de 2 21 pvpv 1应力场文件stif$cv ftt=(1.+pv)*alfa*tn*fact$cv ex=+u/x$cv ey=+v/y$cv exy=+u/y+v/x$cv fm_i=+fact*de_i_j*em

12、_j$cv fm_i=+fm_i-fe_i$cv fsxy=fact*exy*(1.-pv)/2dist=+sx;sx load=+fs_i*ff_i end 应力场VDE文件defia ell b&b ellc str a bstartsin bstartsin asolvsin bsolvsin astress c 文件etaetbetc#elemtype q42dxy 文件非线性薛定鄂能谱方程（0,0）非线性问题nVa(32,32)n使用Newton迭代法线性化后的虚功方程为准备文件准备文件disp u,coor x,y,shap%1%2gaus%3coef unmate va f 2.

13、0;1.0;vect x x yStif$c6 eg=va*(sin(x)*2+sin(y)*2)+3.0*un*2dist=+u/x_i;u/x_i+u;u*egload=+u*(f+2.0*un*3)enddefia nell&STARTsin aIF EXIST END DEL END:2SOLVsin aIF NOT EXIST END GOTO 2ell#elemtype q42dxy程序生成练习（0,0）按照下面的方程用GCN功能生成计算程序，并按给定的参数给出终了时刻为10的压力分布程序生成练习（0,0）将下列方程线性化后给出PDE或VDE文件，再用GCN功能计算在下列边值条件和给定区域下的解