《有限元学习心得》由会员分享，可在线阅读，更多相关《有限元学习心得（18页珍藏版）》请在金锄头文库上搜索。
1、有限元学习心得吴清鸽车辆工程50110802411短短八周的有限元课已经结束。关于有限元，我一直停留在一个很模糊的概念。我知道这是一个各个领域都必须涉及的点，只要有关于分析的，几乎都要涉及有限元。总体来说，这是一门非常重要又有点难度的课程。有限兀方法(finiteelementmethod)或有限兀分析(finiteelementanalysis)，是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将它应用于工程技术中，可成为工程设计和分析的可靠工具。本课程教学基本内容有固体力学和结构力学简介；有限元法基
2、础；桁架、梁、刚架、二维固体、板和壳、三维固体的有限元法；建模技术；热传导问题的有限元分析；PATRAN软件的使用.通过有限元分析课程学习使我了解和掌握了一些有限元知识：1简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程，即表示位移-应变关系的几何方程，表示应力-应变关系的本构方程和表示内力-外力关系的平衡方程。2了解利用能量法形成有限元离散系统方程的基本原理，即哈密尔顿原理。掌握有限元分析的基本方法及步骤，包括域的离散、位移插值、构造形函数、单元有限元方程的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。3具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、
3、板和壳结构、三维固体的有限元法分析技术，包括他们具体的形函数构造，应变矩阵，局部坐标系和整体坐标系中的单元矩阵。各种结构的实例研究。4了解并掌握建立高质量建模所涉及的各种关键技术。包括单元类型的选择，单元畸形的限制，不同阶数单元混用时网格的协调性问题，对称性的应用(平面对称、轴对称、旋转对称、重复对称)，由多点约束方程形成刚域及应用(模拟偏移、不同自由度单元的连接、网格协调性的施加)等，以及多点约束方程的求解。以PATRAN有限元通用软件为例了解一般商业有限元软件的组成及结构。掌握PATRAN软件的基本使用。利用PATRAN软件上机实践完成两个上机练习：刚架结构有限元分析和三维固体有限元分析。
4、课程的具体学习内容：内容：1. 三节点三角形单元：单元分析、总刚度矩阵组装、引入约束条件修正总刚度矩阵、载荷移置、方程求解；2. 四边形单元分析、四节点四面体单元分析、八节点六面体单元分析；#3. 其他常用单元形函数、自由度。1、三节点三角形单元1.1.单元分析1.1.1分析步骤单元分析的任务是建立单元平衡方程，形成单元刚度矩阵。不失一般性，从图1-1三角形离散结构中任取一个单元，设单元编号为e,单元节点按右手法则顺序编号为i,j,m,在定义的坐标系xOy中，节点坐标分别为(xi+yi)，(xj+yj)，(xm+ym)，节点位移和节点力表示如图1-1所示。m取结点位移作基本未知量。由(a)Fe
5、Me位移求结点力:uiviiVufi=iujvVjjuPmUmV_m其中，转换矩阵称为单元刚度矩阵。单元分析的主要目的就是要求出单元刚度矩阵。单元分析的步骤可表示如下：1.1.2位移模式和形函数对于平面问题，单元任意一点的位移可用位移分量u,v描述，他们是坐标x,y的函数。假定三节点单元的位移函数为x,y的线性函数，六个节点位移只能确定六个多项式的系数，所以平面问题的3结点三角形单元的位移函数如下：u=a+ax+ay123v=a+ax+ay456所选用的这个位移函数，将单元内部任一点的位移定为座标的线性函数，位移模式很简单。位移函数写成矩阵形式为：uv1000xya5a6将水平位移分量和结点坐
6、写成矩阵代入位移函数第一式：u=aaxayi12i3iu=aaxayj12j3ju=aaxay2m则有yyiTjymm1ujxyaiiii=1xyVajjj2u1xyammm31xi1xj1xma1=bbbVv52Aijmjacccv6ijmmaaaau11ijmiVa=bbbVu22Aijmjacccu3ijmm最终确定六个待定系数：riraaaaav41ijmiVa=bbb=vf=0NiNi030Njuiviujvjumvm令T*=aiabibcicTaibajbamb则有Va1a1=aibajbambuiVujjjijm22Aijmjabccccacccummmijm3ijmm同样，将垂
7、直位移分量与结点坐标代入位移插值公式：令N=(abxcy)(下标ijm轮换)feje=ijmuiviuvjui2AiiiN称为形态矩阵，N.称为位移的形态函数i113位移函数的收敛性m选择单元位移函数时，应当保证有限元法解答的收敛性，即当网格逐渐加密时，有限元法的解答应当收敛于问题的正确解答。因此，选用的位移模式应当满足下列两方面的条件：(1) 必须能反映单元的刚体位移和常量应变。6个参数a1到6反映了三个刚体位移和三个常量应变。(2) 必须保证相邻单元在公共边界处的位移连续性。(线性函数的特性)114应变矩阵和应力矩阵利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。用结点位
8、移表示单元的应变的表达式为：uxvyuv.yx丄2Abi0ci0cibibj0cj0cjbjbm0cm0cmbmuiviuijrvjum.vm“二Bx(H)(Tx（丫）#根据虚功原理，得Ty*eTFe=J!*htdxdy这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。虚应变可以由结点虚位移求出：=（）代入虚功方程接上式，将应力用结点位移表示出汕=&O=B建立了单元的结点力与结点位移之间的关系，称为单元刚度矩阵。它是6*6矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而
9、改变。1.2总刚度矩阵组装-的集成。整体刚度矩阵K是单元刚度矩阵1、刚度集成法的物理概念：刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点力2、刚度矩阵的集成规则：先对每个单元求出单元刚度矩阵,然后将其中的每个子块到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵K的子块，从而得出结构刚度矩阵K。关键是如何找出中的子块在K中的对应位置。这需要了解单元中的结点编码与结构中的结点编码之间的对应关系。结构中的结点编码称为结点的总码，各个单元的三个结点又按逆时针方向编为i,j,m,称为结点的局部码。单元刚度矩阵中的子块是按结点的局部码排列的，而结构刚度矩阵中的子块是按结点的总码排列的。因此，在单元刚度矩阵中，把结点的局部码换成总码，并把其中的子块按照总码次序重新排列。局部码mKK(2)KJ)圧讪(2Kmm()%(：)K(2)K(2)K(:)j2i以单元为例，局部码悶严对应于总码5,2,4，因此中的子块按照总码重新排列后，得出扩大矩阵为上图所示：用同样的方法可得出其他单元的扩大矩阵、将各单元的扩大矩阵迭加，即得出结构刚度矩阵K：口+=，口集成规则包含搬家和迭加两个环节：1、将
