杆系结构的有限元法分析

有限元单元法

FiniteElementAnalysis第2章杆系结构的有限元法分析精选ppt2.1概述有限单元法的基本思想是从整体到局部，再回到整体，即对我们分析的整体对象，根据其结构特点，对其进行离散化，得到有限个独立的单元，然后对每个单元进行单元分析，最后根据单元分析的结果对结构物进行整体分析，求得结构物的某些参数。在所有结构中，杆系结构是最简单的一类结构，也是我们在工程上最常见的一类结构。如平面桁架、平面刚架、连续梁、空间刚架、空间桁架等都属于此类结构，以此类结构为基础介绍有限单元法的分析过程。首先了解一下有限单元法分析问题的基本步骤。精选ppt第一步：对结构物进行离散化，划分为有限个单元

弯曲杆件系统

截面连续变化的杆件系统

以直代曲

若干微小的等截面杆单元

精选ppt第二步：对各结点和单元进行编码单元划分示意图

精选ppt第三步：建立整体坐标系和各单元的局部坐标系整体坐标系和各单元的局部坐标系精选ppt第四步：对已知参数进行准备和整理对于各单元，需要准备的数据包括：单元截面积:单元长度:单元弹性模量:单元剪切模量:单元惯性矩:等。

精选ppt第五步：对结点位移进行编码结点位移进行编码前处理法后处理法精选ppt第六步：进行单元分析精选ppt

我们进行单元分析的最终目的是要对结构进行整体分析，因此必须由单元特性矩阵构成整体特性矩阵。注意的是，如果局部坐标系与整体坐标系不一致，则需进行坐标变换，将局部坐标系下的单元特性转换为整体坐标系下的单元特性。第七步：进行整体分析，形成整体刚度矩阵精选ppt第八步：引入边界条件边界条件的引入可以使问题具有解的唯一性，否则，我们的问题就是不适定的。精选ppt第九步，求解方程组，计算结构的整体结点位移列阵，并进一步计算各单元的应力分量及主应力、主向。

第十步，求单元内力，对计算成果进行整理、分析，用表格、图线示出所需的位移及应力。精选ppt2.2局部坐标系中杆单元分析

所谓杆件是指从构造上来说其长度远大于其截面尺寸的一维构件。在结构力学上我们通常将承受轴力或扭矩的杆件称为杆，而将承受横向力和弯矩的杆件称为梁。在有限单元法中这两种情况的单元分别称为杆单元和梁单元。但由于在实际工程结构中，同一构件上，上述几种受力状态往往同时存在，因此为方便起见，本书都称之为杆单元。并且，本书所讨论的杆单元均是指等截面直杆单元，对于变截面杆和弯曲杆件，我们在进行单元划分时可以将其分为若干等截面杆单元。因此本书的分析方法仍然对其适应。精选ppt2.2.1拉压杆单元拉压杆单元示意图精选ppt①用结点位移表示单元上任意截面的位移u

其中、为待定系数。

由位移的边界条件：用矩阵表示为：精选ppt精选ppt②进行应力、应变分析

根据应变的定义，有：由虎克定律，其应力为：精选ppt③求单元刚度矩阵

利用虚位移原理求单元刚度矩阵：精选ppt精选ppt2.2.2扭转杆单元扭转杆单元示意图精选ppt精选ppt精选ppt精选ppt2.2.3只计弯曲的杆单元只计弯曲的杆单元示意图精选ppt精选ppt精选ppt精选ppt精选ppt其中的单元刚度矩阵求得为：

精选ppt2.2.4平面一般杆件单元一般杆单元示意图精选ppt精选ppt精选ppt2.2.5空间杆件单元空间杆单元示意图精选ppt精选ppt2.2.6单元刚度矩阵的性质精选ppt精选ppt2.3杆系结构的整体分析2.3.1平面问题坐标变换矩阵

（a）

（b）平面问题两种坐标系下杆端力转换关系示意图

角度转动的正负由右手定则确定，本书中以顺时针方向转动为正精选ppt精选ppt正交矩阵，其逆矩阵等于其转置矩阵精选ppt同样的推导，可以得到两种坐标系下的杆端位移之间的转换关系为：整体坐标系下的单元刚度方程

精选ppt2.3.2空间问题坐标变换矩阵空间问题两种坐标系下杆端力转换关系示意图精选ppt精选ppt精选ppt关系矩阵精选ppt整体坐标系下单元杆端力矩阵与局部坐标系下单元杆端力矩阵具有如下的关系表达式：空间坐标系的单元转换矩阵正交矩阵精选ppt2.3.3杆系结构的整体分析对杆系结构进行单元分析，仅仅是有限元分析中的第一步。我们的目的是要对整个结构进行分析，研究结构的整体性能。因此，在对结构的各单元分析完成后，必须将单元分析的结果进行整合，对结构进行整体分析。整体分析的过程实际上是如何将单元分析的结果进行有效组合，建立整体刚度方程并求解结点位移的过程。根据对结点位移的编码方式，可以采用“先处理法”和“后处理法”来建立整体刚度方程。精选ppt2.3.3.1后处理法

所谓后处理法，就是由单元刚度矩阵形成整体刚度矩阵，建立刚度方程后再引入支承条件，进而求解结点位移的方法。运用这种方法时，假设所有结点位移均为未知量，按照顺序统一进行编码。后处理法位移编码示意图精选ppt精选ppt奇异矩阵，不能求逆矩阵，即可得到无穷多个解

必须引入边界条件

精选ppt边界条件：

精选ppt精选ppt2.3.3.2先处理法

先处理法位移编码示意图利用先处理法对单元结点位移编码时，仅对独立的位移分量按自然数顺序编号，若某些位移分量由于连接条件的限制彼此相等，则编为同一位移号，在支座处，由于刚性约束而使位移分量为零时，则对应的编号为0。精选ppt2.3.3.3杆系结构整体刚度矩阵

刚度集成法首先求出各单元的贡献矩阵然后将它们叠加起来形成整体刚度矩阵边定位，边累加单元的定位数组精选ppt精选ppt精选ppt精选ppt精选ppt精选ppt单元划分、建立局部坐标系和整体坐标系。并对数据进行整理，对单元和结点编号。精选ppt求局部坐标系中的单元刚度矩阵

由于单元①、②的尺寸完全一样，因此其单元刚度矩阵

精选ppt求整体坐标系中的单元刚度矩阵

精选ppt对于单元②，由于其局部坐标系与整体坐标系一致，因此两种坐标系下的单元刚度矩阵相同精选ppt形成整体刚度矩阵精选ppt2.4等效结点荷载和边界条件的处理2.4.1非结点荷载的处理根据有限元方法的离散思想，我们需要将作用于单元上的外荷载（包括集中荷载、面分布荷载、体分布荷载、力偶等）按照虚功等效的原则移植到结点上，成为等效结点荷载。这里的虚功等效，是指原力系与等效结点荷载在任何可能的微小位移（虚位移）上所做的虚功相等。

等效结点荷载

精选ppt精选ppt2.4.2边界条件的处理2.4.2.1铰结点

（a）

（b）

铰结点的处理示意图

精选ppt2.4.2.2弹性支承点

弹性支承点的处理示意图精选ppt2.5杆系结构分析算例精选ppt整体刚度矩阵精选ppt求总结点荷载_局部坐标系单元①精选ppt求总结点荷载_整体坐标系精选ppt刚架等效结点荷载矩阵

直接作用在结点上的荷载总的结点荷载矩阵精选ppt整体刚度方程

精选ppt各单元在整体坐标系下的杆端