杆件结构的有限元法课件

有限元理论与应用第一篇有限元法第一篇有限元法第二章杆件结构的有限元法杆件结构可分为珩杆和梁两种。和其他结构采用铰连接的杆称为珩杆。珩杆的连接处可以自由转动，因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的几何因素主要是截面面积，与截面形状无关。和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动，因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关，而且与截面形状和方位有很大关系。建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。由杆件组成的机构体系称为杆系，如起重机、桥梁等。由珩杆组成的杆系称为珩架，由梁组成的杆系称为刚架。奥运会场馆鸟巢空间立体网架当处理比较复杂的铰支杆系统时，要确定系统在力P的作用下，节点B、C、D和E处的变形。以便计算各杆件的内应力及各杆所受的轴向力，可假设整个杆件系统也具有像式(4—1)中k值一样的刚度，这样在力P的作用下各点的位移就可以用类似式(4—1)的公式计算了，不过．这时的系统刚度应采用一个矩阵来表示，即，同理，各点的位移也应采用一个矩阵来表示，即，再加上矩阵，就构成了称为对应于施加存系统上各节点力的刚度矩阵。问题：1、复杂结构其刚度矩阵是多少阶的？2、如何求出？3、为什么着重讨论系统的刚度矩阵？系统的整体刚度矩阵－求出所受外力作用下各杆件节点处的位移－计算各杆件的受力和应力ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹簧的刚度矩阵是2x2阶的。为求出它们，将图2—4所示弹簧系统看作两个简单的系统，然后合成。一、单个弹簧的刚度矩阵2-2弹簧系统的刚度矩阵二、组合弹簧的刚度矩阵kakbu1,F1u2,F2u3,F31233u1,F1akaF2aF3akbu2＝0u3＝0F1bkakbu2,F2bF3bu1＝0u3＝0F1ckakbF2cu3,F3cu1＝0u2＝0(a)(b)(c)1)只允许节点1有位移u1，力F1a与位移u1之间的关系由于u1＝u2＝0，没有力作用于节点3，因此，考虑弹簧1-2，由静力平衡条件有2)只允许节点2有位移u2，这时由于位移的连续性，每个弹簧在节点2要求有相同的位移，即，弹簧1-2的伸长量与弹簧2-3的缩短量相等。对弹簧1-2有拉力kau2，对弹簧2-3有压力kbu2分别对两弹簧求静力平衡，有3)只允许节点3有位移u3，类似于情况1），有由于节点1、2无位移，有组合弹簧的刚度矩阵4)合成。对整个系统来说有3个节点，每个节点只有一个方向的位移。因此方程式应用如下形式：利用线弹性系统的叠加原理，找出3×3阶刚度矩阵各元素的表达式节点1处的合力节点2处的合力节点3处的合力对成、奇异矩阵（2－8）用同样的方法可以求解具有更多个弹簧的串连系统，推导过程乏味。知道单个弹簧的刚度矩阵－－直接叠加出多个串联系统的总刚度矩阵。按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办法来叠加。以上面两个弹簧系统为例，系统共三个节点，每个节点有一个自由度，因此，该系统总刚度矩阵应该是3×3阶的矩阵。第1个单元的节点号为1和2，则单元刚度矩阵中的元素在总刚度矩阵中应在位置第1行、第2行的第1列，第2列第2个单元的节点号为2和3，则单元刚度矩阵叠加到总刚度矩阵的第2行、第3行的第2列、第3列元素上三、方程求解（约束条件的引入）由式（2－6）和式（2－8）可知，刚度矩阵是一个奇异阵，即它的行列式的值为零，矩阵的逆不存在。

对应线性代数方程组式（2－7）和式（2－9）无定解。

物理概念解释：对整个系统的位移u1、u2和

u3，没有加以限制，从而在任何外力的作用下系统会发生刚体运动。u1＝u2＝

u3＝u,且u没有定值，所以方程无定解。为使方程组有定解，只需给系统加上一定的约束（称为约束条件或边界条件）例如：两弹簧系统，节点1固定不动，有u1＝0，则式（2－9）成为从而可得到定解。通过解上述方程可得到各个节点的位移，利用已求得的位移就可计算出每个弹簧所受力的大小。弹簧1－2受力pa＝ka×（弹簧1－2长度的变化量）

pa＝ka×（u2-u1）

有限元方法求解弹簧系统受力问题的基本步骤：①形成每个单元的刚度矩阵②各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵③引入约束条件④以节点位移为未知量求解线性代数方程组⑤用每个单元的力－位移关系求得单元力。坐标变换为建立整个结构的刚度矩阵，需要在一个共同的统一坐标系（即总体坐标系）中建立平衡方程。由于刚架各单元的空间位置不同，各个单元的局部坐标系一般也不相同。实际杆件系统都是互相成一定角度排列的杆件连接在一起的每个杆件的单元坐标系统所有杆件的都适用的整体坐标系统

12对应局部坐标，x，y对应整体坐标系统对应局部坐标系的位移和作用力，对应整体坐标系的位移和作用力。注意：(1)图中角是从整体坐标系x轴正向起算逆时针转到杆件方向。（2）铰支连接的杆中能承受轴向力和产生轴向位移,因此局部坐标系下，。方便矩阵运算，将力和位移的矩阵用四阶方程表示：将上式从局部坐标系转换到整体坐标系，表示为：类似地可写出节点2处的表达式。用左乘上式两边（2－15）再将式（2－14）代入式（2－15），有单元刚度矩阵