有限元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. 附录… … … … … … … … … … … … … … … … … … …（ PAGE \* MERGEFORMAT - 19 - 问题概述 等参单元是有限元方法中使用最广泛的单元类型。等参单元的位移模式和坐标变换均采用相同的形函数，这种坐标变换叫做等参变换。通过等参变换可以将自然（局部）坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中形状扭曲的单元，因而使得单元有较好的适应性。 本问题首先对平面八结点四边形等参单元的形函数、应力矩阵和等效节点力矩阵、单元应力磨平等的推导和计算，并通过设计FORTRAN求解程序进行编程求解，最后给出算例（）并进行求解，将解与ANSYS解进行比较。在这个过程中，采用了高斯三点积分和高斯两点积分，这种积分方法的求解效率较高而且精度也较好。在问题的最后，尝试去分析引起数值解误差的原因，并分析八结点等参单元的若干特性。 八结点四边形等参单元介绍 边长为2的正方形单元（如下图所示），在其形心处安置一个局部坐标οξη。单元角结点i的坐标（ξi,ηi）分别为±1和0,因此单元四条边界的方程可用简单公式ξ±1=0和η±1=0逐一给出。 图1-1 等参单元示意图 形函数Ni： 位移函数： 坐标变换式： ， 单元应变矩阵 式中——单元节点的位移列阵； 的子块矩阵为 这里记号，。根据复合函数求导规则，有 从而有 N1,ξN1,η=1-η/41-ξ/4 N2,ξN2,η=η-1ξξ2-1/2 N3,ξN3,η=1-η/41+ξ/4 N4,ξN4,η=1-η2/2-η1+ξ N5,ξN5,η=1+η/41+ξ/4 N6,ξN6,η=-ξ1+η1-ξ2/2 N7,ξN7,η=1+η/41-ξ/4 N8,ξN8,η=η2-1/4ηξ-1 代入 Ni,xNi,y=J-1Ni,ξNi,η=1Jx,ξ-y,ξ-x,ηy,ηNi,ξNi,η 可得 N1,xN1,y=1x,ξy,η-y,ξx,ηx,ξ-y,ξ-x,ηy,η1-η/41-ξ/4=14(x,ξy,η-y,ξx,η)x,ξ1-η-y,ξ1-ξ-x,η1-η+y,η1-ξ N2,xN2,y=1x,ξy,η-y,ξx,ηx,ξ-y,ξ-x,ηy,ηη-1ξξ2-1/2=1(x,ξy,η-y,ξx,η)x,ξη-1ξ-y,ξξ2-1/2-x,ηη-1ξ+y,ηξ2-1/2 N3,xN3,y=1x,ξy,η-y,ξx,ηx,ξ-y,ξ-x,ηy,η1-η/41+ξ/4=14(x,ξy,η-y,ξx,η)x,ξ1-η-y,ξ1+ξ-x,η1-η+y,η1+ξ N4,xN4,y=1x,ξy,η-y,ξx,ηx,ξ-y,ξ-x