多场耦合问题的程序生成GCNGIO多场耦合问题的程序生成课件学习目的1.掌握GCN文件的填写格式；2.掌握用GCN功能生成程序的一般过程；3.了解非线性问题线性化的过程。多场耦合问题的程序生成课件GCN文件的结构a NFE文件名 耦合场名 & b NFE文件名 耦合场名c NFE文件名 耦合场名. . . . . . . . . . . .空一行命令流多场耦合问题的程序生成课件库中NFE文件命名基本原则NFE算法文件是FEPG系统用来生成单元计算程序(即E元件程序)和后处理计算程序(即U元件程序)描叙算法，系统提供大量供选择的算法。算法是与方程的类型有关的（方程类型可参看关于PDE文件填写的powerpoint文档），库中NFE文件命名基本原则：ell椭圆par抛物wave波动n*非线性多场耦合问题的程序生成课件nELL.NFE 求解线性椭圆型方程的算法程序。求解线性椭圆型方程的算法程序。nPAR.NFE 求解线性抛物型方程的算法程序，时间求解线性抛物型方程的算法程序，时间 离散采用离散采用Crank-Nicolson格式。格式。nPARB.NFE 求解线性抛物型方程的算法程序，时间求解线性抛物型方程的算法程序，时间离散采用向后差分格式。离散采用向后差分格式。nWAVE.NFE 求解线性双曲型方程的算法程序，时间求解线性双曲型方程的算法程序，时间离散采用离散采用Wave速度法格式速度法格式 nNEWMARK.NFE 求解线性双曲型方程的算法程序，求解线性双曲型方程的算法程序，时间离散采用时间离散采用Newmark格式。格式。nNELL.NFE 求解非线性椭圆型方程的算法程序。求解非线性椭圆型方程的算法程序。nNPAR.NFE 求解非线性抛物型方程的算法程序求解非线性抛物型方程的算法程序,时间时间离散采用离散采用Crank-Nicolson格式。格式。 多场耦合问题的程序生成课件nNPARB.NFE 求解非线性抛物型方程的算法程序，求解非线性抛物型方程的算法程序，时间离散采用向后差分格式。时间离散采用向后差分格式。nNWAVE.NFE 求解非线性双曲型方程的算法程序，求解非线性双曲型方程的算法程序，时间离散采用时间离散采用Wave速度法格式。速度法格式。nNNW.NFE 求解非线性双曲型方程的算法程序，时求解非线性双曲型方程的算法程序，时间离散采用间离散采用Newmark格式。格式。nSTR.NFE 已知位移场求应力场的算法程序，采用已知位移场求应力场的算法程序，采用最小二乘法。也是求解已知场量的梯度的算法程序。最小二乘法。也是求解已知场量的梯度的算法程序。多场耦合问题的程序生成课件命令流命令流段给出耦合场的计算流程。该段由两种命令组成，一种是 DOS命令；另一种是本系统的专有命令。较常用的专有命令可参看FEPG中级教程关于GCN文件填写的部分。多场耦合问题的程序生成课件线性稳态问题命令流段startcsolvcSTART元件程序E元件程序SOLV求解器U元件程序命令流段线性稳态问题流程多场耦合问题的程序生成课件线性动态问题命令流段startcif exist stop del stop:1bftsolvcif not exist 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 goto 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 文件没有第二段对一维问题体单元类型取 l2，坐标系取 1dx；对二维问题，体单元类型取 t3，坐标系取 2dxy；对三维，体单元类型取 w4，坐标系取 3dxyz。关于单元类型和坐标系的说明多场耦合问题的程序生成课件如何生成有限元程序1. GIO GCN文件名多场耦合问题的程序生成课件三个例子1.poisson方程的求解2.热力耦合问题3.非线性问题多场耦合问题的程序生成课件对于如下的对于如下的piosson方程方程边界条件如下多场耦合问题的程序生成课件准备准备poisson.vde文件文件disp u, coor x, y, zshap %1 %2gaus %3load (-6.0)vect x x y zstifdist=+u/x_i;u/x_iend多场耦合问题的程序生成课件准备准备ell.gcn文件文件defia ell &Startsin aSolvsin a 多场耦合问题的程序生成课件准备准备ell.gio文件文件Poisson#elemtype c83dxyz多场耦合问题的程序生成课件1.运行运行gio ell生成全部程序生成全部程序2.2.按按照照某某个个特特定定问问题题，建建好好几几何何模模型型并并施施加加边边界界条条件件完成完成poisson方程的定解方程的定解多场耦合问题的程序生成课件问题描述：平板长1米，宽0.5米；左端温度为0，右端温度为100，下端完全固定。求在此条件下的板的温度分布、变形和应力。板的线膨胀系数 1.0e-5/,弹性模量为1000Mpa,泊松比为0.3，热传导系数为10W/m/。不计板的体力和内热源。热固耦合问题多场耦合问题的程序生成课件问题分析：1.该问题包括了三个物理场，温度场、位移场、应力场。需要准备描述三个物理场场方程的VDE或PDE文件。2.准备描述场之间耦合关系的GCN文件。温度场不考虑固体变形和应力分布，但由于材料的热胀冷缩会引起材料变形和热应力。3.准备GCN文件和PDE/VDE文件联接关系的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温度场ETB.VDE文件多场耦合问题的程序生成课件二维直角坐标下的平面热应力描述方程平衡方程几何方程本构方程多场耦合问题的程序生成课件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.-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 变形场ETA.VDE文件多场耦合问题的程序生成课件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 x,ycoef u,v,tnshap %1 %2gaus %3mass %1 1.0mate 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应力场ETC.VDE文件多场耦合问题的程序生成课件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_j$cv fm_i=+fm_i-fe_i$cv fsxy=fact*exy*(1.-pv)/2dist =+sx;sx*0.0 load=+fs_i*ff_i end 应力场VDE文件多场耦合问题的程序生成课件defia ell b &b ellc str a bstartsin bstartsin asolvsin bsolvsin astress c ET.GCN文件多场耦合问题的程序生成课件etaetbetc#elemtype q42dxy ET.GIO文件多场耦合问题的程序生成课件非线性薛定鄂能谱方程（0,0）非线性问题nVa=2.0nf=1.0(32,32)多场耦合问题的程序生成课件多场耦合问题的程序生成课件n使用Newton迭代法线性化后的虚功方程为多场耦合问题的程序生成课件准备准备ell.vde文件文件disp u, coor x, y, shap %1 %2gaus %3coef unmate va f 2.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)end多场耦合问题的程序生成课件Le.gcndefia nell &STARTsin aIF EXIST END DEL END:2SOLVsin aIF NOT EXIST END GOTO 2多场耦合问题的程序生成课件Le.gioell#elemtype q42dxy多场耦合问题的程序生成课件多场耦合问题的程序生成课件程序生成练习（0,0）按照下面的方程用GCN功能生成计算程序，并按给定的参数给出终了时刻为10的压力分布多场耦合问题的程序生成课件程序生成练习（0,0）将下列方程线性化后给出PDE或VDE文件，再用GCN功能计算在下列边值条件和给定区域下的解多场耦合问题的程序生成课件