平面四边形4结点等参有限单元法(共38页).doc

精选优质文档-倾情为你奉上有限元程序设计平面四边形4结点等参有限单元法程序设计1、程序功能及特点a.该程序采用四边形4节点等参单元，能解决弹性力学的平面应力应变问题。b.前处理采用网格自动划分技术，自动生成单元及结点信息。b.能计算受集中力、自重体力、分布面力和静水压力的作用。c.计算结点的位移和单元中心点的应力分量及其主应力。d.后处理采取整体应力磨平求得各个结点的应力分量。e.算例计算结果与ANSYS计算结果比较，并给出误差分析。f.程序采用Visual Fortran 5.0编制而成。2、程序流程及图框图2- 程序流程图图2-子程序框图其中，各子程序的主要功能为：INPUT输入原始数据HUAFEN自动网格划分，形成COOR(2,NP),X,Y的坐标值与单元信息CBAND形成主元素序号指示矩阵MA（*）SKO形成整体刚度矩阵KCONCR计算集中力引起的等效结点荷载RBODYR