Matlab---四边形八节点等参元计算程序.docVIP

%主运行程序 %**************************************************************************************************** % 四边形八节点等参元计算程序 %全局变量定义 %**************************************************************************************************** global NNODE NPION NELEM NVFIX NFORCE COORD LNODS YOUNG POISS THICK % 单元节点数,总结点数,单元数,受约束自由度数,节点力数,结构节点坐标数组,单元定义数组,弹性模量,泊松比,厚度 global FORCE FIXED ELMENT % FORCE 节点力数组(n,3) n:受力节点数目,(n,1):作用点,(n,2):x方向,(n,3):y方向 % FIXED 约束信息数组(n,3) n:受约束节点数目, (n,1):约束点 (n,2)与(n,3)分别为约束点x方向和y方向的约束情况,受约束为1否则为0 global BMATX DMATX SMATX AREA % 单元应变矩阵(3*16),单元弹性矩阵(16*16),单元应力矩阵(3*16),单元面积 global ASTIF ASLOD ASDISP ESTIF COORT ELEDISP STRESS ASTRESS % 总体刚度矩阵,总体荷载向量,节点位移向量,单元刚度矩阵 global FP1 %数据文件指针 % 主程序 %**************************************************************************************************** format short e %设定输出类型 FP1=fopen(in.txt,rt); %打开数据文件 READING %% %调用子程读取数据 ASSEMBLE %% %调用子程生成单刚,组成总刚并加入约束信息 ASTIF; %约束处理后总刚 FORMLOAD %% %调用子程生成荷载向量 ASDISP=ASTIF\ASLOD %计算节点位移向量 WRITESTRESS %% %调用子程计算单元应力 fclose(FP1); %关闭文件 function READING() %读取数据文件中的数据 global NNODE NPION NELEM NVFIX NGAUS NNODE COORD COORT LNODS ASDISP YOUNG POISS THICK FIXED NFORCE FORCE FP1 NPION=fscanf(FP1,%d,1) NELEM=fscanf(FP1,%d,1) NFORCE=fscanf(FP1,%d,1) NVFIX=fscanf(FP1,%d,1), NGAUS=fscanf(FP1,%d,1) NNODE=fscanf(FP1,%d,1), YOUNG=fscanf(FP1,%e,1) POISS=fscanf(FP1,%f,1) THICK=fscanf(FP1,%f,1) LNODS=fscanf(FP1,%d,[9,NELEM]) COORT=fscanf(FP1,%d,[1,2*NNODE]) COORD=fscanf(FP1,%f,[3,NPION]) FORCE=fscanf(FP1,%f,[3,NFORCE]) FIXED=fscanf(FP1,%d,[3,NVFIX]) function ASSEMBLE() %生成单刚,组成总刚并加入约束信息 %所引用的全局变量 global NPION NELEM NVFIX LNODS ASTIF THICK ESTIF NVIFX global