计算力学课程设计.docxVIP

前言

有限单元法是在当今工程分析中获得最广泛应用的数值计算方法，由于其通用性和有效性，受到工程技术界的高度重视。伴随着计算机科学和技术的快速开展，现已成为计算机辅助设计和计算机辅助制造的重要组成局部。

由于有限元法是通过计算机实现的，因此有限元程序的编制及相关软件的研发就变得尤为重要。从二十世纪五十年代以来，有限元软件的开展按目的和用途可分为专用软件和大型通用商业软件，而且软件往往集成了网格自动划分，结果分析和显示等前后处理功能，而且随着时间的开展，大型通用商业软件的功能由线性扩展到非线性，由结构扩展到非结构等等，这一系列强大功能的实现与运用都要求我们对有限元法的根底理论知识有较为清楚的认识以及对程序编写的根本能力有较好掌握。

一、有限元分析程序的理论根底

作为弹性力学微分方程的等效积分形式，虚位移原理和虚应

原理分别是平衡方程与力的边界条件和几何方程与位移边界条件的等效积分形式。将物理方程引入虚位移原理和虚应力原理可以分别导出最小位能原理和最小余能原理，它们本质上和等效积分的伽辽金“弱”形式相一致，这是建立弹性力学有限原方程一般表达格式的理论根底。

对于通过弹性力学变分原理建立的弹性力学问题有限元方法，其未知场变量是位移，以节点位移为根本未知量，并以最小位能原理为根底建立的有限单元为位移元。弹性力学平面问题有限元分析表达格式的建立步骤一般为：

对所给的模型进行单元离散，一般采用三角形单元或四边形单元，不过由于四边形单元具有更高精度，应用更为普遍，一般而言一次或二次单元已足以满足精度要求。离散后对单元进行编号，并且给定单元各节点的整体编码以及局部坐标系下按逆时针局部编码。

在局部坐标系下对各单元构造形函得到由单元各节点位移表示的单元位移形式，进而得到单元刚度矩阵，利用等参元性质和雅阁比矩阵进行组集，建立整体刚度矩阵。

建立单元等效结点载荷列阵并组集成结构节点等效载荷列阵

得到单元各节点位移与结构节点等效载荷列阵的线代方程组，求解得到位移，进而可得应力应变。

根据以上一般步骤，编制相应的计算机程序并采用数值积分方法处理有关数学计算就可以得到完整的弹性力学问题有限元分析程序。

二、平面问题有限元分析教学程序

本程序可对二维弹性力学问题行有限元分析计算。计算采用的单元形式为四边形四节点单元或者四边形八节点单元，对于对称和非对称矩阵采用变带宽存储方法，最终输出结果为单元各节点的位移。

C---------------------------------------------------------------------

C|

CFEA2DP----AFINITEELEMENTANAYLYSISPROGRAMFOR|

C2DELASTICPROBLEMS|

CTANGENTMATRIXISSTOREDWITHVARIOUSBANDMETHOH|

CTHISPROGRAMISUSEDTODEMONSTRATETHEUSAGEOFVRIOUSBAND|

CSTRORAGESCHEMOFSYMMETRICANDUNSYMMETRICTANGENTSMATRIX|

C|

CXIAOBAOHE|

CATCHONGQINGUNIVERSITY,CHINA.(12/6/2012)|

C|

C---------------------------------------------------------------------

C

PROGRAMFEA2DP

C

CA(1)-A(N1-1):X(NDM,NUMAP);A(N1)-A(N2-1):F(NDF,NUMAP)

CA(N2)-A(N3-1):B(