有限元编程大作业报告.docx

本科生实验报告书四节点等参单元有限元分析的FORTRAN程序目录问题概述………………………………………………(1)四节点四边形等参单元介绍……………………………(1)单元应力磨平方法介绍…………………………………(4)程序流程设计………………………………………… (6)4.1 程序设计概述4.2 流程图程序结构及程序说明……………………………………(8)程序应用及算例分析……………………………………(9)6.1 算例概述6.2 算例ANSYS求解6.3 算例程序数值解6.4 算例分析7. 总结…………………………………………………(15)7. 附录……………………………………………………()问题概述等参单元是有限元方法中使用最广泛的单元类型。等参单元的位移模式和坐标变换均采用相同的形函数，这种坐标变换叫做等参变换。通过等参变换可以将自然（局部）坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中形状扭曲的单元，因而使得单元有较好的适应性。本问题首先对平面四节点四边形等参单元的形函数、应力矩阵和等效节点力矩阵、应力磨平公式等的推导和计算求解。并通过设计FORTRAN求解程序进行编程求解，最后给出算例（受集中荷载的悬臂梁）并进行求解，将解与ANSYS的解进行比较。在这个过程中，采用了高斯三点积分和高斯两点积分，这种积分方法的求解效率较高而且精度也较好。在问题的最后，尝试去分析引起