杆系结构的有限元法分析.pptx

有限元单元法FiniteElementAnalysis,土木工程与建筑学院:陈月顺CollegeofCivilEngineeringandArchitecture,第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后处理法,所谓后处理法，就是由单元刚度矩阵形成整体刚度矩阵，建立刚度方程后再引入支承条件，进而求解结点位移的方法。运用这种方法时，假设所有结点位移均为未知量，按照顺序统一进行编码。,后处理法位移编码示意图,奇异矩阵，不能求逆矩阵，即可得到无穷多个解,必须引入边界条件,边界条件：,2.3.3.2先处理法,先处理法位移编码示意图,利用先处理法对单元结点位移编码时，仅对独立的位移分量按自然数顺序编号，若某些位移分量由于连接条件的限制彼此相等，则编为同一位移号，在支座处，由于刚性约束而使位移分量为零时，则对应的编号为0。,2.3.3.3杆系结构整体刚度矩阵,刚度集成法,首先求出各单元的贡献矩阵,然后将它们叠加起来形成整体刚度矩阵,边定位，边累加,单元的定位数组,单元划分、建立局部坐标系和整体坐标系。并对数据进行整理，对单元和结点编号。,求局部坐标系中的单元刚度矩阵,由于单元、的尺寸完全一样，因此其单元刚度矩阵,求整体坐标系中的单元刚度矩阵,对于单元，由于其局部坐标系与整体坐标系一致，因此两种坐标系下的单元刚度矩阵相同,形成整体刚度矩阵,2.4等效结点荷载和边界条件的处理,2.4.1非结点荷载的处理,根据有限元方法的离散思想，我们需要将作用于单元上的外荷载（包括集中荷载、面分布荷载、体分布荷载、力偶等）按照虚功等效的原则移植到结点上，成为等效结点荷载。这里的虚功等效，是指原力系与等效结点荷载在任何可能的微小位移（虚位移）上所做的虚功相等。,等效结点荷载,2.4.2边界条件的处理,2.4.2.1铰结点,（a）（b）铰结点的处理示意图,2.4.2.2弹性支承点,弹性支承点的处理示意图,2.5杆系结构分析算例,整体刚度矩阵,求总结点荷载_局部坐标系,单元,求总结点荷载_整体坐标系,刚架等效结点荷载矩阵,直接作用在结点上的荷载,总的结点荷载矩阵,整体刚度方程,各单元在整体坐标系下的杆端位移,计算杆端内力,根据结点平衡条件，计算支座反力,2.6程序设计方法,