一、平面三角形3结点有限元程序.docVIP

一、平面三角形３结点有限元程序 1、程序名：FEM3.FOR，FEM3.EXE 2、程序功能 本程序能计算弹性力学的平面应力问题和平面应变问题；能考虑自重和结点集中力两种荷载的作用，在计算自重时轴取垂直向上为正； 能处理非零已知位移，如支座沉降的作用。主要输出的内容包括：结点位移、单元应力、主应力、第一主应力与轴的夹角以及约束结点的支座反力。 程序采用Visual Fortran 5.0编制而成，输入数据全部采用自由格式。 3、程序流程及框图 图１-1 程序流程图 图１-2 程序图 INPUT——输入结点坐标、单元信息和材料参数； MR——形成结点自由度序号矩阵； FORMMA——形成指标矩阵MA（N）并调用其他功能子程序，相当于主控程序； DIV——取出单元的3个结点号码和该单元的材料号并计算单元的等； MGK——形成整体劲度矩阵并按一维存放在SK（NH）中； LOAD——形成整体结点荷载列阵F； OUTPUT——输出结点位移或结点荷载； TREAT——由于有非零已知位移，对K和F进行处理； DECOMP——整体劲度矩阵的分解运算； FOBA——前代、回代求出未知结点位移； ERFAC——计算约束结点的支座反力； KRS——计算单元劲度矩阵中的子块Krs。 4、输入数据及