有限元程序设计(共18页).doc

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