三角形单元有限元程序设计.doc

三角形单元有限元程序设计 %******************************************************************* % NNODE 单元节点数 NPION 总结点数 % NELEM 单元数 NVFIX 受约束边界点数 % FIXED 约束信息数组 NFORCE 节点力数 % FORCE 节点力数组 COORD 结构节点坐标数组 % LNODS 单元定义数组 % YOUNG 弹性模量 POISS 泊松比 % THICK 厚度 B 单元应变矩阵(3*6) % D 单元弹性矩阵(3*3) S 单元应力矩阵(3*6) % A 单元面积 ESTIF 单元刚度矩阵 % ASTIF 总体刚度矩阵 ASLOD 总体荷载向量 % ASDISP 节点位移向量 ELEDISP 单元节点位移向量 % STRESS 单元应力 FP1 数据文件指针 %******************************************************************* format short e %设定输出类型 clear %清除内存变量 FP1=