多场耦合问题程序生成.pptVIP

ELL.NFE 求解线性椭圆型方程的算法程序。 PAR.NFE 求解线性抛物型方程的算法程序，时间 离散采用Crank-Nicolson格式。 PARB.NFE 求解线性抛物型方程的算法程序，时间离散采用向后差分格式。 WAVE.NFE 求解线性双曲型方程的算法程序，时间离散采用Wave速度法格式 NEWMARK.NFE 求解线性双曲型方程的算法程序，时间离散采用Newmark格式。 NELL.NFE 求解非线性椭圆型方程的算法程序。 NPAR.NFE 求解非线性抛物型方程的算法程序,时间离散采用Crank-Nicolson格式。 NPARB.NFE 求解非线性抛物型方程的算法程序，时间离散采用向后差分格式。 NWAVE.NFE 求解非线性双曲型方程的算法程序，时间离散采用Wave速度法格式。 NNW.NFE 求解非线性双曲型方程的算法程序，时间离散采用Newmark格式。 STR.NFE 已知位移场求应力场的算法程序，采用最小二乘法。也是求解已知场量的梯度的算法程序。 Le.gcn defi a nell STARTsin a IF EXIST END DEL END :2 SOLVsin a IF NOT EXIST END GOTO 2 Le.gio ell #elemtype q4 2dxy 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)/2 dist =+[sx;sx]*0.0 load=+[fs_i]*ff_i end 应力场VDE文件 defi a ell b b ell c str a b startsin b startsin a solvsin b solvsin a stress c ET.GCN文件 eta etb etc #elemtype q4 2dxy ET.GIO文件 非线性薛定鄂 能谱方程 （0,0） 非线性问题 Va=2.0 f=1.0 (32,32) 使用Newton迭代法线性化后的虚功方程为 准备ell.vde文件 disp u, coor x, y, shap %1 %2 gaus %3 coef un mate va f 2.0;1.0; vect x x y Stif $c6 eg=va*(sin(x)**2+sin(y)**2)+3.0*un**2 dist=+[u/x_i;u/x_i]+[u;u]*eg load=+[u]*(f+2.0*un**3) end * * 多场耦合问题的程序生成 GCN文件 GIO文件 学习目的 1.掌握GCN文件的填写格式； 2.掌握用GCN功能生成程序的一般过程； 3.了解非线性问题线性化的过程。 GCN文件的结构 a NFE文件名 耦合场名 b NFE文件名 耦合场名 c NFE文件名 耦合场名 . . . . . . . . . . . . 空一行 命令流 库中NFE文件命名基本原则 NFE算法文件是FEPG系统用来生成单元计算程序(即E元件程序)和后处理计算程序(即U元件程序)描叙算法，系统提供大量供选择的算法。算法是与方程的类型有关的（方程类型可参看关于《PDE文件填写》的powerpoint文档），库中NFE文件命名基本原则： ell——椭圆 par——抛物 wave——波动 n*——非线性 命令流 命令流段给出耦合场的计算流程。该段由两种命令组成，一种是 DOS命令；另一种是本系统的专有命令。较常用的专有命令可参看《FEPG中级教程》关于GCN文件填写的部分。 线性稳态问题命令流段 startc solvc START元件程序 E元件程序 SOLV求解器 U元件程序 命令流段 线性稳态问题流程 线性动态问题命令流段 startc if exist stop del stop :1 bft solvc if not exist stop goto 1 命令流段 线性动态问题流程 START元件程序 BFT元件程序 E元件程序 SOLV求解器 U元件程序 时间计算是否结束 否 是 非线性稳态问题命令流段 startc if exist end del end :1 solvc if not exist end goto 1 命令流段 非线性稳态问题