第二章__杆系结构的有限元法解决方案.pptx

有限元单元法 Finite Element Analysis 第2章　杆系结构的有限元法分析 2.1　概　　述 有限单元法的基本思想是从整体到局部，再回到整体，即对我们分析的整体对象，根据其结构特点，对其进行离散化，得到有限个独立的单元，然后对每个单元进行单元分析，最后根据单元分析的结果对结构物进行整体分析，求得结构物的某些参数。 在所有结构中，杆系结构是最简单的一类结构，也是我们在工程上最常见的一类结构。如平面桁架、平面刚架、连续梁、空间刚架、空间桁架等都属于此类结构，以此类结构为基础介绍有限单元法的分析过程。 首先了解一下有限单元法分析问题的基本步骤。 第一步：对结构物进行离散化，划分为有限个单元 弯曲杆件系统 截面连续变化的杆件系统 以直代曲 若干微小的等截面杆单元 第二步：对各结点和单元进行编码 单元划分示意图 第三步：建立整体坐标系和各单元的局部坐标系 整体坐标系和各单元的局部坐标系 第四步：对已知参数进行准备和整理 对于各单元，需要准备的数据包括： 单元截面积: 单元长度: 单元弹性模量: 单元剪切模量: 单元惯性矩: 等。 第五步：对结点位移进行编码 结点位移进行编码 前处理法 后处理法 第六步：进行单元分析 我们进行单元分析的最终目的是要对结构进行整体分析，因此必须由单元特性矩阵构成整体特性矩阵。注意的是，如果局部坐标系与整体坐标系不一致，则需进行坐标变换，将局部坐标系下的单元特性转换为整体坐标系下的单元特性。 第七步：进行整体分析，形成整体刚度矩阵 第八步：引入边界条件 边界条件的引入可以使问题具有解的唯一性，否则，我们的问题就是不适定的。 第九步，求解方程组，计算结构的整体结点位移列阵 ，并进一步计算各单元的应力分量及主应力、主向。第十步，求单元内力，对计算成果进行整理、分析，用表格、图线示出所需的位移及应力。 2.2　局部坐标系中杆单元分析 所谓杆件是指从构造上来说其长度远大于其截面尺寸的一维构件。在结构力学上我们通常将承受轴力或扭矩的杆件称为杆，而将承受横向力和弯矩的杆件称为梁。在有限单元法中这两种情况的单元分别称为杆单元和梁单元。但由于在实际工程结构中，同一构件上，上述几种受力状态往往同时存在，因此为方便起见，本书都称之为杆单元。并且，本书所讨论的杆单元均是指等截面直杆单元，对于变截面杆和弯曲杆件，我们在进行单元划分时可以将其分为若干等截面杆单元。因此本书的分析方法仍然对其适应。 2.2.1　拉压杆单元 拉压杆单元示意图 ① 用结点位移表示单元上任意截面的位移u 其中 、 为待定系数。 　　　　　 由位移的边界条件： 用矩阵表示为： ② 进行应力、应变分析 根据应变的定义，有： 由虎克定律，其应力为： ③ 求单元刚度矩阵 利用虚位移原理求单元刚度矩阵： 2.2.2　扭转杆单元 扭转杆单元示意图 2.2.3　只计弯曲的杆单元 只计弯曲的杆单元示意图 其中的单元刚度矩阵求得为： 　　　　　 2.2.4　平面一般杆件单元 一般杆单元示意图 2.2.5　空间杆件单元 空间杆单元示意图 2.2.6　单元刚度矩阵的性质 2.3　杆系结构的整体分析 2.3.1　平面问题坐标变换矩阵 （a） （b） 平面问题两种坐标系下杆端力转换关系示意图 角度转动的正负由右手定则确定，本书中以顺时针方向转动为正 正交矩阵，其逆矩阵等于其转置矩阵 同样的推导，可以得到两种坐标系下的杆端位移之间的转换关系为： 整体坐标系下的单元刚度方程 2.3.2　空间问题坐标变换矩阵 空间问题两种坐标系下杆端力转换关系示意图 关系矩阵 整体坐标系下单元杆端力矩阵与局部坐标系下单元杆端力矩阵具有如下的关系表达式： 空间坐标系的单元转换矩阵 正交矩阵 2.3.3　杆系结构的整体分析 对杆系结构进行单元分析，仅仅是有限元分析中的第一步。我们的目的是要对整个结构进行分析，研究结构的整体性能。因此，在对结构的各单元分析完成后，必须将单元分析的结果进行整合，对结构进行整体分析。整体分析的过程实际上是如何将单元分析的结果进行有效组合，建立整体刚度方程并求解结点位移的过程。根据对结点位移的编码方式，可以采用“先处理法”和“后处理法”来建立整体刚度方程。 2.3.3.1　后处理法 所谓后处理法，就是由单元刚度矩阵形成整体刚度矩阵，建立刚度方程后再引入支承条件，进而求解结点位移的方法。运用这种方法时，假设所有结点位移均为未知量，按照顺序统一进行编码。 后处理法位移编码