大作业报告参考2有限元学习心得.doc

大作业报告参考2有限元学习心得.doc

《大作业报告参考2有限元学习心得.doc》由会员分享，可在线阅读，更多相关《大作业报告参考2有限元学习心得.doc（14页珍藏版）》请在冰豆网上搜索。

有限元学习心得

吴清鸽车辆工程50110802411

短短八周的有限元课已经结束。

关于有限元，我一直停留在一个很模糊的概念。

我知道这是一个各个领域都必须涉及的点，只要有关于CAE分析的，几乎都要涉及有限元。

总体来说，这是一门非常重要又有点难度的课程。

有限元方法（finiteelementmethod）或有限元分析（finiteelementanalysis），是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。

将它用于在科学研究中，可成为探究物质客观规律的先进手段。

将它应用于工程技术中，可成为工程设计和分析的可靠工具。

本课程教学基本内容有固体力学和结构力学简介；有限元法基础；桁架、梁、刚架、二维固体、板和壳、三维固体的有限元法；建模技术；热传导问题的有限元分析；PATRAN软件的使用.

通过有限元分析课程学习使我了解和掌握了一些有限元知识：

1．简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程，即表示位移-应变关系的几何方程，表示应力-应变关系的本构方程和表示内力-外力关系的平衡方程。

2．了解利用能量法形成有限元离散系统方程的基本原理，即哈密尔顿原理。

掌握有限元分

析的基本方法及步骤，包括域的离散、位移插值、构造形函数、单元有限元方程的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。

3．具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、板和壳结构、三维固体的有限元法分析技术，包括他们具体的形函数构造，应变矩阵，局部坐标系和整体坐标系中的单元矩阵。

各种结构的实例研究。

4．了解并掌握建立高质量建模所涉及的各种关键技术。

包括单元类型的选择，单元畸形的限制，不同阶数单元混用时网格的协调性问题，对称性的应用（平面对称、轴对称、旋转对称、重复对称），由多点约束方程形成刚域及应用（模拟偏移、不同自由度单元的连接、网格协调性的施加）等，以及多点约束方程的求解。

以PATRAN有限元通用软件为例了解一般商业有限元软件的组成及结构。

掌握PATRAN软件的基本使用。

利用PATRAN软件上机实践完成两个上机练习：

刚架结构有限元分析和三维固体有限元分析。

课程的具体学习内容：

内容：

1.三节点三角形单元：

单元分析、总刚度矩阵组装、引入约束条件修正总刚度矩阵、载荷移置、方程求解；

2.四边形单元分析、四节点四面体单元分析、八节点六面体单元分析；

3.其他常用单元形函数、自由度。

1、三节点三角形单元

1.1.单元分析

1.1.1分析步骤

单元分析的任务是建立单元平衡方程，形成单元刚度矩阵。

不失一般性，从图1-1三角形离散结构中任取一个单元，设单元编号为e，单元节点按右手法则顺序编号为i,j,m,在定义的坐标系xOy中，节点坐标分别为（xi+yi），（xj+yj），（xm+ym），节点位移和节点力表示如图1-1所示。

取结点位移作基本未知量。

由结点位移求结点力：

其中，转换矩阵称为单元刚度矩阵。

单元分析的主要目的就是要求出单元刚度矩阵。

单元分析的步骤可表示如下：

1.1.2位移模式和形函数

对于平面问题，单元任意一点的位移可用位移分量u,v描述，他们是坐标x,y的函数。

假定三节点单元的位移函数为x,y的线性函数，六个节点位移只能确定六个多项式的系数，所以平面问题的3结点三角形单元的位移函数如下：

所选用的这个位移函数，将单元内部任一点的位移定为座标的线性函数，位移模式很简单。

位移函数写成矩阵形式为：

将水平位移分量和结点坐写成矩阵：

代入位移函数第一式：

令 则有

A为三角形单元

[T]的伴随矩阵为

令 则有

同样，将垂直位移分量与结点坐标代入位移插值公式：

最终确定六个待定系数：

令 （下标i，j，m轮换）

[N]称为形态矩阵，

Ni称为位移的形态函数

1.1.3位移函数的收敛性

选择单元位移函数时，应当保证有限元法解答的收敛性，即当网格逐渐加密时，有限元法的解答应当收敛于问题的正确解答。

因此，选用的位移模式应当满足下列两方面的条件：

（1）必须能反映单元的刚体位移和常量应变。

6个参数到反映了三个刚体位移和三个常量应变。

（2）必须保证相邻单元在公共边界处的位移连续性。

（线性函数的特性）

1.1.4应变矩阵和应力矩阵

利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。

用结点位移表示单元的应变的表达式为：

[B]矩阵称为几何矩阵

由物理方程，可以得到单元的应力表达式：

为应力矩阵

1.1.5单元刚度矩阵

讨论单元内部的应力与单元的结点力的关系，导出用结点位移表示结点力的表达式。

由应力推算结点力，需要利用平衡方程。

用虚功方程表示出平衡方程。

考虑上图三角形单元的实际受力，任意虚设位移，节点位移结点力和内部应力为：

与内部应变为：

令实际受力状态在虚设位移上作虚功，外力虚功为

微小矩形的内力虚功为

根据虚功原理，得

这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。

虚应变可以由结点虚位移求出：

代入虚功方程

接上式，将应力用结点位移表示出

有

令

则

建立了单元的结点力与结点位移之间的关系，称为单元刚度矩阵。

它是6*6矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。

1.2总刚度矩阵组装

整体刚度矩阵是单元刚度矩阵的集成。

1、刚度集成法的物理概念：

刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点力。

2、刚度矩阵的集成规则：

先对每个单元求出单元刚度矩阵，然后将其中的每个子块送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵[K]的子块，从而得出结构刚度矩阵[K]。

关键是如何找出中的子块在[K]中的对应位置。

这需要了解单元中的结点编码与结构中的结点编码之间的对应关系。

结构中的结点编码称为结点的总码，各个单元的三个结点又按逆时针方向编为i,j,m,称为结点的局部码。

单元刚度矩阵中的子块是按结点的局部码排列的，而结构刚度矩阵中的子块是按结点的总码排列的。

因此，在单元刚度矩阵中，把结点的局部码换成总码，并把其中的子块按照总码次序重新排列。

以单元②为例，局部码i,j,m对应于总码5,2,4，因此中的子块按照总码重新排列后，得出扩大矩阵为上图所示：

用同样的方法可得出其他单元的扩大矩阵

将各单元的扩大矩阵迭加，即得出结构刚度矩阵[K]：

集成规则包含搬家和迭加两个环节：

1、将单元刚度矩阵中的子块搬家，得出单元的扩大刚度矩阵。

2、将各单元的扩大刚度矩阵迭加，得出结构刚度矩阵[K]。

1.3引入约束条件修正总刚度矩阵

整体刚度矩阵[K]求出后，结构的结点力{F}可表示为：

在无支杆的结点处，结点力就等于已知的结点载荷。

在有支杆的结点处，则求结点力时，还应把未知的支杆反力考虑在内。

如果用{P}表示结点载荷和支杆反力组成的向量，则结点的平衡方程为

根据支承条件对平衡方程加以处理。

先考虑结点n有水平支杆的情况。

与结点n水平方向对应的平衡方程是第2n-1个方程，

根据支承情况，上式应换成，即在[K]中，第2n-1行的对角线元素应改为1，该行全部非对角线元素应改为0。

在{P}中，第2n-1个元素应改为0。

此外，为了保持矩阵[K]的对称性，则第2n-1列全部非对角线元素也改为0。

同理，如果结点n有竖向支杆，则平衡方程的第2n个方程应改为，为此，在矩阵[K]中，第2n行的对角线元素改为1，该行全部非对角线元素改为0，同时，第2n列全部非对角线元素也改为0。

在{P}中，第2n个元素改为0。

1.4载荷移置

将载荷移置到结点上，必须遵循静力等效的原则。

静力等效是指原载荷与结点载荷在任意虚位移上做的虚功相等。

在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。

单元的虚位移可以用结点的虚位移表示为

令结点载荷为

集中力的移置

如图所示，在单元内任意一点作用集中力

由虚功相等可得

由于虚位移是任意的，则

体力的移置

令单元所受的均匀分布体力为

由虚功相等可得：

2、四边形、四节点四面体、八节点六面体单元分析

2.1四边形单元分析

四边单元可取矩形作为研究对象，矩形单元也是一种常用的单元，它采用了比常应变三角形单元次数更高的位移模式，因而可以更好地反映弹性体中的位移状态和应力状态。

矩形单元1234如图所示，其边长分别为2a和2b，两边分别平行于x、y轴。

若取该矩形的四个角点为节点，因每个节点位移有两个分量，所以矩形单元共有8个自由度。

采用三角形单元中的方法，同样可以完成对这种单元的力学特性分析。

然而，如果我们引入一个局部坐标系x、h，那么就可以推出比较简洁的结果。

在局部坐标系中，节点i的坐标是（xi,hi），其值分别为±1。

取位移模式：

用与三角形单元相似的方法建立形函数，则位移模式可写成：

式中：

由几何方程可以求得单元的应变

可推出：

式中

由虎克定律我们可以得出用节点位移表示的单元应力，即

式中：

对于平面应力问题

有：

若将单元刚度矩阵写成分块形式：

则其中的子矩阵可按下式进行计算：

四边形单元的节点位移与单元节点力之间的关系仍为：

其中载荷列阵{R}e与上节中的（c）式相同，仍可按上式计算等效节点力。

但是，需要注意的是，矩形单元有四个节点（1，2，3，4），所以{R}e具有8个元素，即：

2.2四节点四面体单元分析

2.2.1单元划分及位移模式

如图1所示的四面体单元，单元结点的编码为i,j,m,n。

每个结点的位移具有三个分量u,v,w。

这样单元结点的位移列阵可表示成：

单元的位移模式采用线性多项式：

式中，为待定系数，由单元结点的位移和坐标决定。

将四个结点的坐标（xi,yi,zi）、（xj,yj,zj）、（xm,ym,zm）、（xn,yn,zn）和结点位移（ui,vi,wi）、（uj,vj,wj）、（um,vm,wm）、（un,vn,wn）代入

（2）式可得12个联立方程，解方程组便可求出。

将这十二个系数回代到式，则得到由结点位移和形函数表示的单元内任一点的位移表达式：

式中

Ni，Nj，Nm，Nn为四面体单元的形函数

位移模式可以用矩阵形式表示：

式中，[I]为三阶单位阵，[N]为形函数矩阵。

上式即为单元结点位移和单元任意点位移之间的关系。

2.2.2单元应变和应力

知道单元内任意一点位移后，可利用几何方程确定单元内该点的应变。

将位移矩阵式代入空间问题几何方程得：

其中

上式表明几何矩阵[B]中的元素都是常量，因此单元中的应变也是常量。

也就是说，采用线性位移模式的四面体单元是常应变单