有限元教学程序及使用说明.docVIP

有限元教学程序及使用说明 附录 有限元教学程序及使用说明 ?A-1平面三角形,结点有限元程序 1、程序名:FEM3.FOR，FEM3.EXE 2、程序功能 本程序能计算弹性力学的平面应力问题和平面应变问题;能考虑自重和结点集中力两种 荷载的作用，在计算自重时轴取垂直向上为正; 能处理非零已知位移，如支座沉降的作y 用。主要输出的内容包括:结点位移、单元应力、主应力、第一主应力与轴的夹角以及约x束结点的支座反力。 程序采用Visual Fortran 编制而成，输入数据全部采用自由格式。 3、程序流程及框图 启 动 输入结点坐标、单元信息等 形成整体刚度矩阵K 形成整体结点荷载R 求解方程组 Ka,R 计算各单元应力及主应力 FFFFFFFFFFFFFFFFFFFE NFEFEN 输出结点位移和单元应力 li力 结 束 图A-1 程序流程图 277 图A-2 程序框图 其中，各子程序的功能如下: INPUT——输入结点坐标、单元信息和材料参数; MR——形成结点自由度序号矩阵; FORMMA——形成指标矩阵MA(N)并调用其他功能子程序，相当于主控程序; DIV——取出单元的3个结点号码和该单元的材料号并计算单元的，等; bciiMGK——形成整体劲度矩阵并按一维存放在SK(NH)中; LOAD——形成整体结点荷载列阵F; OUTPUT——输出结点