有限元编程综述资料.docVIP

有限元编程综述 指导老师： 日 期： 摘 要： 有限元法，主要用来解决复杂结构中力与位移的关系。其基本思想是将具有无限个自由度的连续的求解区域离散为具有有限个自由度、且按一定方式(节点)相互连接在一起的离散体(单元)，即将连续体假想划分为数目有限的离散单元，而单元之间只在数目有限的指定点处相互联结，用离散单元的集合体代替原来的连续体。本文主要介绍有限元分析的及编程的主要通过一个具体相关内容。 ：有限元流程 所谓有限元法（FEA）ANSYS、ABAQUS、CAD 的实体造型，优化设计等接口。现在可编写工程应用软件的计算机语言较多，其中C语言是一个较为优秀的语言，很容易满足现在有限元分析程序编程的要求。 C语言最初是为操作系统、编译器以及文字处理等编程而发明的。随着不断完善，它已应用到其它领域，包括工程应用软件的编程。近年来，C语言已经成为计算机领域最普及的一个编程语言，几乎世界上所有的计算机都装有C的编译器，从PC机到巨型机到超巨型的并行机，C与所有的硬件和操作系统联系在一起。用C 编写的程序，可移植性极好，几乎不用作多少修改，就可在任何一台装有ANSI、C编译器的计算机上运行。C既是高级语言，也是低级语言，也就是说，可用它作数值计算，也可用它对计算机存储进行操作。 编程思想 本程序采用C语言编程，编制平面四边形四节点等参元程序，用以求解平面结构问题。程序采用二维等带宽存储整体刚度矩阵，乘大数法引入约束，等带宽高斯消去法求解位移，然后求中间高斯点的应力，最后用绕节点平均法讲单元应力等效到节点上，再将结果写到tecplot文件中。 在有限元程序中，变量数据需赋值的可分为节点信息，单元信息，载荷信息等。对于一个节点来说，需以下信息：节点编号（整型），节点坐标（实型），节点已知位移（实型），节点载荷（实型），边界条件（实型）等。同样，对于一个单元来说，需以下信息：单元的节点联接信息（整型），材料信息（弹性模量，泊松比等）（实型）等。 在FORTRAN 程序中，以上这些变量混合在一起，很难辨认，使程序的可读性不好，如需要进行单元网络的自适应划分，节点及单元的修改将非常困难。在进行C语言编译过程中，采用结构struct 使每个节点信息存储在一个结构体数组中，提高程序的可读性，使数据结构更趋于合理。 平面四节点等参单元介绍 四节点等参单元实际单元与基本单元的映射关系如图 31所示 坐标的映射关系为： 其位移模式和坐标的映射有相同的插值函数，形函数为： 单元应变矩阵为： 上式一般简写为： 其中的子块矩阵为 由于是、的函数，在中的、要按照复合函数来求导，即 从而有： 因此，单元应力矩阵为： 单元刚度矩阵为： 其中积分采用三点高斯积分， 其中，和或是加权系数，和是单元内的坐标.。对于三点高斯积分，高斯积分点的位置： ，，。 结构刚度矩阵为： 结构结点荷载列阵为： 注意，对于上两式中的理解不是简单的叠加而是按照对应的自由度集成。 总刚平衡方程： 从式上式求出： 4 有限元分析流程 一个典型的有限元分析过程主要包括以下几个步骤： 读输入数据，定义节点及单元数组。 由边界条件计算方程个数，赋值荷载列阵。 读入在带状存储的总刚度矩阵中单元和载荷信息。 定义总刚度阵数组。 组装总刚度阵。 解方程组。 其流程图可见下图： 5计算流程图 6 程序变量及函数说明 6.1 控制信息 np：1，x方向；u2，y方向 n=nf*np ：结构的节点位移总数 ndf=nd*nf ：每个单元的节点自由度数 6.2 输入的原始数据 x(np):节点的x方向坐标 y(np):节点的y方向坐标 me（nd，ne）：单元节点的总体编号 nrr(nr1+nr2) ：约束为零的位移所对应的总体位移编号 p（n）：载荷向量 nu（nu1+nu2）：位移载荷 mat（6，nm）：材料参数 6.3程序中的其他标识符 LD（n）：存放结构刚度阵所以主对角线元素在A（nn）中的序号 IS（ndf）：单元节点位移和节点力在总体位移阵列和载荷阵列中对应的序号 EK（ndf，ndf）：总体坐标系下的单元刚度矩阵 A（nn）：架构刚度阵下三角变带宽一维压缩存储的数组 nn：数组A的元素个数 RSTG(3):高斯积分点的值 H(3):高斯积分点的加权系数 S（6，ne）：各单元的应力分量 XJAC（2,2）：雅阁比矩阵 RJAC（2,2）：雅阁比矩阵的逆 PN（2,4）：4个节点处形函数对局部坐标的导数 DNX（2,4）：4个节点处形函数对整体坐标的导数 FUN（4）：形函数的值 7 计算结果与Abaqus分析结果的比较 7.1、中间带圆孔平面应力板的分析 宽40m，长50m，圆孔位于板中心，半