第十一讲有限元法程序实现(11).pptxVIP

1． 雅可比迭代法 ;D=[8,0,0;0,-11,0;0,0,12;]; L=[0,0,0;4,0,0;6,3,0;]; V=[0,-3,2;0,0,-1;0,0,0;]; B=-D^(-1)*(L+V); C=D^(-1)*F; u0=[0;0;0;]; u1=u0; u2=u0; k=0; e=1; while(e10^-6) u2=B*u1+C; e=norm(u1-u2); k=k+1; u1=u2; end;3． 高斯—赛德尔迭代 ;clear; clc; K=[8,-3,2; 4,-11,-1; 6,3,12;]; F=[20; 33; 36]; n=3; D=[8,0,0;0,-11,0;0,0,12;]; L=[0,0,0;4,0,0;6,3,0;]; V=[0,-3,2;0,0,-1;0,0,0;]; u0=[0;0;0;]; u1=u0; u2=u0; k=0; e=1; ;高斯赛德尔迭代收敛的充分条件：;坐标数据表格(Coor0) 坐标数据表格存放在由前处理程序生成的节点坐标文件 Coor0 中。 坐标数据表格内容包含： 节点总数（整型数） 节点坐标数（整型数） 所有节点坐标值（实型数） ;节点规格数表格(id0) 节点规格数表格存放在前处理程序生成的节点规格数文件 Id0 中。 节点规格数表格内容包含： 节点总数（整型数） 节点自由度数（整型数） 节点自由度规格数（整型数） 本系统用 1，0，-1 和大于 1 的整型数来表示节点规格数，它们的意义分别 是： 1：表示节点上有自由度，但无约束 0：表示节点上无自由度 -1：表示节点上有自由度，而且有约束 ;指定节点位移和节点荷载信息表格(disp0) 指定节点位移和荷载信息数据表格存放在前处理程序生成的指定节点位移和荷载信息数据文件 Disp0 中。 指定节点位移和节点荷载信息表格包含内容： 节点总数（整型数） 节点自由度数（整型数） 节点指定位移和荷载值（实型数） 说明： Id0 文件和 Disp0 文件都是对节点上的约束进行说明，二者相互对应： （1）节点规格数表格中规定为“-1”的节点自由度给出它的指定位移值。 （2）节点规格数表格中规定为“1”的节点自由度给出它的指定荷载值，如果没有荷载应输入零值。 ;初始值表格 FEPG 系统最多允许有三个初值表格文件：Disp1,Disp2,Disp3。 （1）对于不依赖于时间、非线性问题，要给出初始迭代值 Disp1 文件。 （2）对于依赖时间问题，要给出初始时刻值 Disp1,Disp2,Disp3 文件。 例如求解依赖于时间的温度场问题，要给出初始时刻温度值 Disp1 文件。 而对于求解某些用到特殊算法的依赖于时间的问题，例如用 NEWMARK 方法求解波动方程的问题则要给出 Disp1（初始时刻位移值），Disp2(初始 时刻速度值)，Disp3(初始时刻加速度值)文件。 ;Disp1,Disp2,Disp3 文件的输入内容如下： Disp1 文件： 节点总数（整型数） 节点自由度数（整型数） 节点初始时刻位移（实型数） Disp2 文件： 节点总数（整型数） 节点自由度数（整型数） 节点初始时刻速度（实型数） Disp3 文件： 节点总数（整型数） 节点自由度数（整型数） 节点初始时刻加速度（实型数） ;单元信息数据 (elem0) 单元信息数据是由多个数据表格组成，存放于前处理程序生成的 elem0 单 元信息文件中。 单元信息数据表格与材料有关，每种类型单元有两个数据表格： 单元 节点数据表格和材料数据表格。 单元节点数据表格的内容包含： 单元总数（整型数） 单元节点数（整型数） 单元节点号及材料编号（整型数） 材料数据表格的内容包含： 材料种类数（整型数） 材料参数个数（整型数） 材料参数值（实型数） ;FEPG6.0软件安装;FEPG软件数值模拟的一般步骤;静电场的例子;ell #elemtype q4 2dxy;defi a ell STARTsin a SOLVsin a;disp u coor x,y shap %1 %2 gaus %3 mate eu1 fz 3.5416d-11;0.0; stif dist=[u/x;u/x]*eu1+[u/y;u/y]*eu1 load=[u]*fz end; 渗流问题;dam;defi a ell STARTsin a SOLVsin a;disp u coor x,y mate pk eq 5.0d-7 0.0 shap %1 %2 gaus %3 stif dist