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

1、多场耦合问题的程序生成,GCN文件,GIO文件,学习目的 1.掌握GCN文件的填写格式； 2.掌握用GCN功能生成程序的一般过程； 3.了解非线性问题线性化的过程。,GCN文件的结构,a NFE文件名 耦合场名 u/x_i end,准备ell.gcn文件,defi a ell 0.1;0.0; load q vect x x y stif dist=+u/x_i;u/x_i*ek end,温度场ETB.VDE文件,二维直角坐标下的平面热应力描述方程,平衡方程,几何方程,本构方程,disp u,v coor x,y func ex ey exy shap %1 %2 gaus %3 coef t

2、n mate 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*alfa vect em ex ey vect u u v vect fd fx fy vect fe ft*tn ft*tn matr de 2 2 1 pv pv 1,变形场ETA.VDE文件,FUNC ex=+u/x ey=+v/y exy=+u/y+v/x stif dist =+em_i;em_j*de_i_j*fact+exy;exy*fact*(1.-pv)/2 load=+u_i

3、*fd_i+em_i*fe_i end,变形场VDE文件,应力场描述方程,已知变形和温度分布后，代入几何方程求应变，再通过本构方程求出热应力,几何方程,本构方程,disp sx,sy,sxy coor x,y coef u,v,tn shap %1 %2 gaus %3 mass %1 1.0 mate 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 sxy vect ff fsx fsy fsxy vect em ex ey vect fm fsx fsy vect

4、fe ftt ftt matr de 2 2 1 pv pv 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)/2 dist =+sx;sx*0.0 load=+fs_i*ff_i end,应力场VDE文件,defi a ell b 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,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,程序生成练习,（0,0）,按照下面的方程用GCN功能生成计算程序，并按给定的参数