《2022年有限元学习心得 .pdf》由会员分享，可在线阅读，更多相关《2022年有限元学习心得 .pdf（11页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、实用总结1/11 有限元学习心得吴清鸽车辆工程短短八周的有限元课已经结束。关于有限元，我一直停留在一个很模糊的概念。我知道这是一个各个领域都必须涉及的点，只要有关于分析的，几乎都要涉及有限元。总体来说，这是一门非常重要又有点难度的课程。有限元方法()或有限元分析()，是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将它应用于工程技术中，可成为工程设计和分析的可靠工具。本课程教学基本内容有固体力学和结构力学简介。有限元法基础。桁架、梁、刚架、二维固体、板和壳、三维固体的有限元法。建模技术。热传导问题的有限
2、元分析。软件的使用.通过有限元分析课程学习使我了解和掌握了一些有限元知识：简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程，即表示位移应变关系的几何方程，表示应力应变关系的本构方程和表示内力外力关系的平衡方程。了解利用能量法形成有限元离散系统方程的基本原理，即哈密尔顿原理。掌握有限元分析的基本方法及步骤，包括域的离散、位移插值、构造形函数、单元有限元方程的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、板和壳结构、三维固体的有限元法分析技术，包括他们具体的形函数构造，应变矩阵，局部坐标系和整体坐标系中的单元矩
3、阵。各种结构的实例研究。了解并掌握建立高质量建模所涉及的各种关键技术。包括单元类型的选择，单元畸形的限制，不同阶数单元混用时网格的协调性问题，对称性的应用（平面对称、轴对称、旋转对称、重复对称），由多点约束方程形成刚域及应用（模拟偏移、不同自由度单元的连接、网格协调性的施加）等，以及多点约束方程的求解。以有限元通用软件为例了解一般商业有限元软件的组成及结构。掌握软件的基本使用。利用软件上机实践完成两个上机练习：刚架结构有限元分析和三维固体有限元分析。课程的具体学习内容：内容：1.三节点三角形单元：单元分析、总刚度矩阵组装、引入约束条件修正总刚度矩阵、载荷移置、方程求解。2.四边形单元分析、四节
4、点四面体单元分析、八节点六面体单元分析。3.其他常用单元形函数、自由度。名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 11 页 -实用总结2/11、三节点三角形单元1.1.单元分析分析步骤单元分析的任务是建立单元平衡方程，形成单元刚度矩阵。不失一般性，从图三角形离散结构中任取一个单元，设单元编号为，单元节点按右手法则顺序编号为 ,在定义的坐标系中，节点坐标分别为()，()，()，节点位移和节点力表示如图所示。取结点位移作基本未知量。由结点位移求结点力：其中，转换矩阵称为单元刚度矩阵。单元分析的主要目的就是要求出单元刚度矩阵。单元分析的步骤可表示如下：位移模式和形函数对于平面问题
5、，单元任意一点的位移可用位移分量,描述，他们是坐标,的函数。假定三节点单元的位移函数为,的线性函数，六个节点位移只能确定六个多项式的系数，所以平面问题的结点三角形单元的位移函数如下：所选用的这个位移函数，将单元内部任一点的位移定为座标的线性函数，位移模式很简单。位移函数写成矩阵形式为：将水平位移分量和结点坐写成矩阵：代入位移函数第一式：令则有为三角形单元的伴随矩阵为名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 11 页 -实用总结3/11 令则有同样，将垂直位移分量与结点坐标代入位移插值公式：最终确定六个待定系数：令（下标，轮换）称为形态矩阵，称为位移的形态函数位移函数的收敛性
6、选择单元位移函数时，应当保证有限元法解答的收敛性，即当网格逐渐加密时，有限元法的解答应当收敛于问题的正确解答。因此，选用的位移模式应当满足下列两方面的条件：()必须能反映单元的刚体位移和常量应变。个参数到反映了三个刚体位移和三个常量应变。()必须保证相邻单元在公共边界处的位移连续性。(线性函数的特性)应变矩阵和应力矩阵利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。用结点位移表示单元的应变的表达式为：矩阵称为几何矩阵名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 11 页 -实用总结4/11 由物理方程，可以得到单元的应力表达式：为应力矩阵单元刚度矩阵讨论单元内
7、部的应力与单元的结点力的关系，导出用结点位移表示结点力的表达式。由应力推算结点力，需要利用平衡方程。用虚功方程表示出平衡方程。考虑上图三角形单元的实际受力，任意虚设位移，节点位移结点力和内部应力为：与内部应变为：令实际受力状态在虚设位移上作虚功，外力虚功为微小矩形的内力虚功为根据虚功原理，得这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。虚应变可以由结点虚位移求出：代入虚功方程接上式，将应力用结点位移表示出有令则建立了单元的结点力与结点位移之间的关系，称为单元刚度矩阵。它是*矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性
8、常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。总刚度矩阵组装ivimUjUiUmvjvmj*ivi*mU*jU*iU*mv*jvmjys*xy*y*x?geexytxs(a)?|?2?|(b)D?D|?名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 11 页 -实用总结5/11 整体刚度矩阵是单元刚度矩阵的集成。、刚度集成法的物理概念：刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点力。、刚度矩阵的集成规则：先对每个单元求出单元刚度矩阵，然后将其中的每个子块送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵 的子块，从而得出结构刚度矩阵。关键是
9、如何找出中的子块在 中的对应位置。这需要了解单元中的结点编码与结构中的结点编码之间的对应关系。结构中的结点编码称为结点的总码，各个单元的三个结点又按逆时针方向编为,称为结点的局部码。单元刚度矩阵中的子块是按结点的局部码排列的，而结构刚度矩阵中的子块是按结点的总码排列的。因此，在单元刚度矩阵中，把结点的局部码换成总码，并把其中的子块按照总码次序重新排列。以单元为例，局部码对应于总码，因此中的子块按照总码重新排列后，得出扩大矩阵为上图所示：用同样的方法可得出其他单元的扩大矩阵将各单元的扩大矩阵迭加，即得出结构刚度矩阵集成规则包含搬家和迭加两个环节：、将单元刚度矩阵中的子块搬家，得出单元的扩大刚度矩
10、阵。、将各单元的扩大刚度矩阵迭加，得出结构刚度矩阵。引入约束条件修正总刚度矩阵整体刚度矩阵 求出后，结构的结点力 可表示为：在无支杆的结点处，结点力就等于已知的结点载荷。在有支杆的结点处，则求结点力时，还应把未知的支杆反力考虑在内。如果用 表示结点载荷和支杆反力组成的向量，则结点的平衡方程为根据支承条件对平衡方程加以处理。先考虑结点有水平支杆的情况。与结点水平方向对应的平衡方程是第个方程，根据支承情况，上式应换成，即在中，第行的对角线元素应改为，该行全部非对角线元素应改为。在 中，第个元素应改为。此外，为了保持矩阵 的对称性，则第列全部非对角线元素也改为。同理，如果结点有竖向支杆，则平衡方程的
11、第个方程应改为，为此，)1(jjK)2(mmK)2(miK)1(jmK)1(jiK)2(jmK)4(iiK)4(m iK)4(jiK2m1j2m4i126543216543局部码总码321i,j,m431j,m,i432m,j,i321ijm1j431jmi432mji4i)1(m mK)2(jjK)3(iiK)1(miK)3(imK)1(iiK)3(mmK)4(jjK)2(jiK)3(ijK)3(mjK)4(jmK)2(iiK)3(jjK)4(m mK名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 11 页 -实用总结在矩阵 中，第行的对角线元素改为，该行全部非对角线元素改为，
12、同时，第列全部非对角线元素也改为。在 中，第个元素改为。载荷移置将载荷移置到结点上，必须遵循静力等效的原则。静力等效是指原载荷与结点载荷在任意虚位移上做的虚功相等。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。单元的虚位移可以用结点的虚位移表示为令结点载荷为集中力的移置如图所示，在单元内任意一点作用集中力由虚功相等可得由于虚位移是任意的，则体力的移置令单元所受的均匀分布体力为由虚功相等可得：、四边形、四节点四面体、八节点六面体单元分析四边形单元分析四边单元可取矩形作为研究对象，矩形单元也是一种常用的单元，它采用了比常应变三角形单元次数更高的位移模式，因而可以更好地反映弹性体中的
13、位移状态和应力状态。矩形单元如图所示，其边长分别为和，两边分别平行于、轴。若取该矩形的四个角点为节点，因每个节点位移有两个分量，所以矩形单元共有个自由度。采用三角形单元中的方法，同样可以完成对这种单元的力学特性分析。然而，如果我们引入一个局部坐标系、，那么就可以推出比较简洁的结果。在局部坐标系中，节点的坐标是(,)，其值分别为。取位移模式：用与三角形单元相似的方法建立形函数，则位移模式可写成：式中：由几何方程可以求得单元的应变可推出：式中由虎克定律我们可以得出用节点位移表示的单元应力，即u2(U2)u3(U3)v3(V3)v2(V2)v1(V1)v4(V4)u1(U1)u4(U4)xoyo12
14、432a2b名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 11 页 -实用总结7/11 式中：对于平面应力问题有：若将单元刚度矩阵写成分块形式：则其中的子矩阵可按下式进行计算：四边形单元的节点位移与单元节点力之间的关系仍为：其中载荷列阵 与上节中的()式相同，仍可按上式计算等效节点力。但是，需要注意的是，矩形单元有四个节点（，），所以具有个元素，即：四节点四面体单元分析单元划分及位移模式如图所示的四面体单元，单元结点的编码为,。每个结点的位移具有三个分量,。这样单元结点的位移列阵可表示成：单元的位移模式采用线性多项式：式中，为待定系数，由单元结点的位移和坐标决定。将四个结点的坐
15、标（,）、（,）、（,）、（,）和结点位移（,）、（,）、（,）、（,）代入（）式可得个联立方程，解方程组便可求出。将这十二个系数回代到式，则得到由结点位移和形函数表示的单元内任一点的位移表达式：式中名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 11 页 -实用总结8/11，为四面体单元的形函数位移模式可以用矩阵形式表示：式中，为三阶单位阵，为形函数矩阵。上式即为单元结点位移和单元任意点位移之间的关系。单元应变和应力知道单元内任意一点位移后，可利用几何方程确定单元内该点的应变。将位移矩阵式代入空间问题几何方程得：其中上式表明几何矩阵 中的元素都是常量，因此单元中的应变也是常量。
16、也就是说，采用线性位移模式的四面体单元是常应变单元。将上式代入物理方程，就得到单元的应力列阵：式中：为四面体单元的应力矩阵，其分块形式为：其中单元刚度矩阵对于四面体单元，利用虚功原理，采用类似平面问题的处理方法可以得到其单刚矩阵。其中：为单元刚度矩阵名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 11 页 -实用总结9/11 写成分块形式为：式中子矩阵 由下式计算：可以看出，单元刚度矩阵是由单元结点的坐标和单元材料的弹性常数所决定的，是一个常数矩阵。八节点六面体单元分析形函数与坐标变换形函数坐标变换位移插值函数与几何矩阵简记为：单元刚度矩阵与等效节点载荷向量单元刚度矩阵可以表示为
17、：将上式中的替换为则有：写成高斯积分形式为：名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 11 页 -实用总结10/11、其他常用单元形函数、自由度3.1 轴对称单元轴对称结构体可以看成由任意一个纵向剖面绕着纵轴旋转一周而形成。此旋转轴即为对称轴，纵向剖面称为子午面，如图表示一圆柱体的子午面被分割为若干个三角形单元，再经过绕对称轴旋转，圆柱体被离散成若干个三棱圆环单元，各单元之间用圆环形的铰链相连接。对于轴对称问题，采用圆柱坐标较为方便。以弹性体的对称轴为轴，其约束及外载荷也都对称于轴，因此弹性体内各点的各项应力分量、应变分量和位移分量都与环向坐标无关，只是径向坐标和轴向坐标的
18、函数。也就是说，在任何一个过轴的子午面上的位移、应变和应力的分布规律都相同。因此轴对称问题可把三维问题简化为以（，）为自变量的二维问题。zjabriijmmdc由于轴对称性，弹性体内各点只可能存在径向位移和轴向位移。此时，位移、只是、的函数，而环向位移。即：轴对称问题的物理方程可写为：由于轴对称性，我们只需分析任意一个子午面上的位移、应力和应变情况。其有限元分析计算步骤和平面问题相似。首先进行结构区域的有限元剖分。采用的单元是三角形、矩形或任意四边形环绕对称轴旋转一周而得到的整圆环，通常采用的单元是三角形截面的整圆环。在单元类型确定之后，单元剖分可以在子午面内进行，如图表示的子午面被分割为若干
19、个三角形，绕对称轴旋转后即形成若干个三棱圆环单元。相邻的单元由圆环形的铰链相连接。单元的棱边都是圆，故称为结圆。每个结圆与平面的交点称为结点。名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 11 页 -实用总结11/11 这样，各单元在子午面平面上形成三角形网格，就如同平面问题中在平面上的网格一样。采用位移法有限元分析，其基本未知量为结点位移。单元的结点位移列阵如下：对于每一个环形单元，需要假定其位移模式。仿照平面三角形单元，取线性位移模式：类 似 于 平 面 三 角 形 单的 结 点 坐 标 及 结 点 位 移代入式中,可以解出六个待定系数。再将这些待定系数回代到式中，就可以
20、得到由结点位移和形函数所表示的单元内任一点的位移表达式：其中形函数：形函数矩阵的表达式为：有上面分析可知，轴对称单元自由度有六个。以上就是关于课程汇总报告的全部内容。通过这八周的学习，我已经对有限元的基础有了一个大致掌握，关于用有限元进行具体分析也掌握了一些最基本的方法。其中应用到很多矩阵变换之中的知识，我会加强这方面知识的巩固。在以后的研究方向重，我也会对有限元分析的方法勤加练习。这门课对我以后的课题方向和分析方法有着举足轻重的作用。感谢雷老师严谨认真的教学，把理论课学习与上机练习紧密结合起来，是我们更加容易掌握要点，更加容易记住方法。在此表示衷心的感谢。名师资料总结-精品资料欢迎下载-名师精心整理-第 11 页，共 11 页 -