有限元p8单元大作业报告教程.docxVIP

本科生实验报告书八节点等参单元有限元分析的FORTRAN程序目录问题概述……………………………………………（1）八节点四边形等参单元介绍…………………………（1）单元应力磨平方法介绍………………………………（7）程序流程设计…………………………………………（8）4.1 程序设计概述4.2 流程图程序变量说明及程序结构……………………………（10）5.1 子程序说明程序应用及算例分析…………………………………（11）6.1 算例概述6.2 算例ANSYS求解6.3 算例程序数值解6.4 算例分析7. 总结 …………………………………………………（18）7. 附录…………………………………………………（问题概述等参单元是有限元方法中使用最广泛的单元类型。等参单元的位移模式和坐标变换均采用相同的形函数，这种坐标变换叫做等参变换。通过等参变换可以将自然（局部）坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中形状扭曲的单元，因而使得单元有较好的适应性。本问题首先对平面八结点四边形等参单元的形函数、应力矩阵和等效节点力矩阵、单元应力磨平等的推导和计算，并通过设计FORTRAN求解程序进行编程求解，最后给出算例（）并进行求解，将解与ANSYS解进行比较。在这个过程中，采用了高斯三点积分和高斯两点积分，这种积分方法的求解效率较高而且精度也较好。在问题的最后，尝试去分析引起数值解误差的原因，并分析八结点等参单元的若干特性。八结点四边形等参单元介绍边长为2的正方形单元（如下图所示），在其形心处安置一个局部坐标。单元角结点i的坐标（,）分别为,因此单元四条边界的方程可用简单公式和逐一给出。 图1-1 等参单元示意图 形函数：位移函数：坐标变换式：，单元应变矩阵式中——单元节点的位移列阵；的子块矩阵为这里记号，。根据复合函数求导规则，有从而有代入可得迭入应变矩阵可得应变矩阵B因此，单元内的应力可以表示为其中单元刚度矩阵为其中积分采用三点高斯积分，（高斯积分点的总数），和是加权系数，和是单元内的坐标.。对于三点高斯积分，高斯积分点的位置：,，，。单元等效节点荷载结构刚度矩阵结构结点荷载列阵总刚平衡方程从上式求出将回代入和中，得到和。等效节点力集中力 将集中荷载作用点取为结点，集中荷载就可以直接转化为等效结点力。体积力 等效结点力按下式计算其中积分采用高斯2点积分，（高斯积分点的总数），和或是加权系数，和是单元内的坐标。对于2点高斯积分，高斯积分点的位置： ; ; ; 表面力 单元的ij边上两个结点的等效结点力按下式计算应力磨平方法介绍 对与平面八节点等参单元，进行应力磨平时，单元结点数可只取4个，即用二维双线性单元进行应力磨平。磨平式各应力分量分别进行，这时应力磨平插值函数Ni应采用双线性函数，即??　 　　 （3.1）　在此情况下，4个结点上的应力可由高斯点上的应力外推得到，即令在2×2高斯积分点上有。4个高斯积分点的座标(见图3-1)如下:　　图3-1 二维等参单元将高斯点坐标代入(3.1)式得到下面的等式：　　　　　　　　　　 (3.2）式中等号左边的应力列阵是有限元中已求出的4个高斯点相应的应力分量；是磨平后应力的结点值；转换矩阵由(3.1)式的第二式代入高斯点坐标后的插值函数值构成。由（3.2）式求逆可得（3.3）其中（3.4）各应力分量均可用（3.3）求解。这种改进结点应力的方法亦称之为应力插值外推。求得改进的应力结点值后，如需要求单元内部的应力值仍可按(3.1)第一式进行计算。　 采用单元局部应力磨平的方法，对于同一结点，由不同相邻单元求得的应力改进值通常是不相同的。可把相关单元求得的改进结点值再取平均作为最后的结点应力值。程序流程设计程序设计概述本程序的设计以四节点等参单元的FORTRAN程序为基础，在其架构之上修改而来。由于八节点等参单元在边线中点有节点，所以在计算改进节点值时要另外计算，然后把相关单元求得的改进结点值再取平均作为最后的结点应力值。流程图程序变量说明及程序结构子程序说明[P4INPUT] 读入单元数据、节点数据、节点约束条件及各类单元属性[DEA] 形成主对角元素地址[P8SSM] 形成总体刚度矩阵[P8STIFB] 计算单元刚度矩阵[P8BMATB] 计算当前单元的应变矩阵[P8MODB] 计算当前单元的弹性矩阵[P8DBE] 计算应力矩阵[P8XJACM] 计算当前单元的形函数和雅可比行列式[BOUNDARY] 边界条件处理[LDLT] LDLT分解结构刚度矩阵[P8LOAK] 计算单元荷载矩阵形成节点总荷载向量用置0置1法处理边界条件时产生的