杆件结构的有限元法

1、杆件结构的有限元法第1页，共28页，2022年，5月20日，3点46分，星期四第一篇 有限元法第2页，共28页，2022年，5月20日，3点46分，星期四第一篇 有限元法第二章 杆件结构的有限元法第3页，共28页，2022年，5月20日，3点46分，星期四 当结构长度尺寸比两个截面方向的尺寸大得多时，这类结构称为杆件。工程中常见得轴、支柱、螺栓、加强肋以及各类型钢等都属于杆件。第4页，共28页，2022年，5月20日，3点46分，星期四杆件结构可分为珩杆和梁两种。和其他结构采用铰连接的杆称为珩杆。珩杆的连接处可以自由转动，因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的几何因素主要是

2、截面面积，与截面形状无关。和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动，因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关，而且与截面形状和方位有很大关系。建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。第5页，共28页，2022年，5月20日，3点46分，星期四由杆件组成的机构体系称为杆系，如起重机、桥梁等。由珩杆组成的杆系称为珩架，由梁组成的杆系称为刚架。第6页，共28页，2022年，5月20日，3点46分，星期四奥运会场馆鸟巢空间立体网架第7页，共28页，2022年，5月20日，3点46分，星期四工

3、程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹簧。弹簧系统力F与弹簧伸长量 （位移）之间关系由胡克定律有式中k为弹簧的刚度，是弹簧的固有参数。它对应于力位移图中F- 关系直线的斜率。当k和力F已知时，可由下式求出弹簧伸长量弹簧力位移间关系(41)2-1 引 言第8页，共28页，2022年，5月20日，3点46分，星期四 当处理比较复杂的铰支杆系统时，要确定系统在力P的作用下，节点B、C、D和E处的变形。以便计算各杆件的内应力及各杆所受的轴向力，可假设整个杆件系统也具有像式(41)中k值一样的刚度，这样在力P的作用下各点的位移就可以用类似式(41)的公式计算了，不过这时的系统刚度应采

4、用一个矩阵来表示，即 ，同理，各点的位移也应采用一个矩阵来表示，即 ，再加上矩阵 ，就构成了称为对应于施加存系统上各节点力的刚度矩阵。第9页，共28页，2022年，5月20日，3点46分，星期四 问题： 1、复杂结构其刚度矩阵是多少阶的？ 2、如何求出？ 3、为什么着重讨论系统的刚度矩阵？ 系统的整体刚度矩阵求出所受外力作用下各杆件节点处的位移计算各杆件的受力和应力第10页，共28页，2022年，5月20日，3点46分，星期四ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹簧的刚度矩阵是2x 2阶的。为求出它们，将图24所示弹簧系统看作两个简单的系统，然后合成。一、单个弹簧的刚度矩

5、阵2-2 弹簧系统的刚度矩阵第11页，共28页，2022年，5月20日，3点46分，星期四由力的平衡有ku1F1aF2aA A(a)u2=0ku1=0F1bF2bu2B B(b)ku1F1u2F2A AB B1)只有节点1可以变形,点2固定2)只有节点2可以变形,点1固定3)根据线弹性系统的叠加原理，叠加1) 2)两种情况，就得到与原始问题一样的结构，如图（c），叠加结果为：(c)作用于节点1上的合力作用于节点2上的合力刚度矩阵对成、奇异矩阵（25）（26）第12页，共28页，2022年，5月20日，3点46分，星期四二、组合弹簧的刚度矩阵kakbu1,F1u2,F2u3,F31233u1,F

6、1akaF2aF3akbu20u30F1bkakbu2,F2bF3bu10u30F1ckakbF2cu3,F3cu10u20(a)(b)(c)1) 只允许节点1有位移u1，力F1a与位移u1之间的关系由于u1 u20，没有力作用于节点3，因此，考虑弹簧1-2，由静力平衡条件有2) 只允许节点2有位移u2，这时由于位移的连续性，每个弹簧在节点2要求有相同的位移，即，弹簧1-2的伸长量与弹簧2-3的缩短量相等。对弹簧1-2 有拉力kau2，对弹簧2-3 有压力kbu2分别对两弹簧求静力平衡，有3) 只允许节点3有位移u3，类似于情况1），有由于节点1、2无位移，有第13页，共28页，2022年，5

7、月20日，3点46分，星期四组合弹簧的刚度矩阵4) 合成。对整个系统来说有3个节点，每个节点只有一个方向的位移。因此方程式应用如下形式：利用线弹性系统的叠加原理，找出33阶刚度矩阵各元素的表达式节点1处的合力节点2处的合力节点3处的合力对成、奇异矩阵（28）第14页，共28页，2022年，5月20日，3点46分，星期四 用同样的方法可以求解具有更多个弹簧的串连系统，推导过程乏味。 知道单个弹簧的刚度矩阵直接叠加出多个串联系统的总刚度矩阵。第15页，共28页，2022年，5月20日，3点46分，星期四知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵对整个系统来说有3个节点，将上述方程扩大成3阶方

8、程：整个系统有3个节点（位移），将上述方程扩大成3阶方程，按矩阵相加原理将两式叠加，（29）矩阵扩大办法单元数量增多时，相应扩大后的矩阵就相当大，扩大后的非零元素在矩阵的什么位置，概念上就不很清楚了。第16页，共28页，2022年，5月20日，3点46分，星期四按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办法来叠加。以上面两个弹簧系统为例，系统共三个节点，每个节点有一个自由度，因此，该系统总刚度矩阵应该是33阶的矩阵。第1个单元的节点号为1和2，则单元刚度矩阵中的元素在总刚度矩阵中应在位置第1行、第2行的第1列，第2列第2个单元的节点号为2和3，则单元刚度矩阵叠加到总刚度矩阵的第

9、2行、第3行的第2列、第3列元素上第17页，共28页，2022年，5月20日，3点46分，星期四三、方程求解（约束条件的引入）由式（26）和式（28）可知，刚度矩阵是一个奇异阵，即它的行列式的值为零，矩阵的逆不存在。对应线性代数方程组式（27）和式（29）无定解。物理概念解释：对整个系统的位移u1、 u2和 u3，没有加以限制，从而在任何外力的作用下系统会发生刚体运动。u1u2 u3u,且u没有定值，所以方程无定解。为使方程组有定解，只需给系统加上一定的约束（称为约束条件或边界条件）例如：两弹簧系统，节点1固定不动，有u10，则式（29）成为从而可得到定解。通过解上述方程可得到各个节点的位移，

10、利用已求得的位移就可计算出每个弹簧所受力的大小。第18页，共28页，2022年，5月20日，3点46分，星期四弹簧12受力 paka（弹簧12长度的变化量） paka（u2-u1） 有限元方法求解弹簧系统受力问题的基本步骤：形成每个单元的刚度矩阵各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵引入约束条件以节点位移为未知量求解线性代数方程组用每个单元的力位移关系求得单元力。第19页，共28页，2022年，5月20日，3点46分，星期四23 杆件系统的有限元法 一、铰支杆系统的有限元计算格式上面求解弹簧系统的有限元方法可以直接用力求解受轴向力的杆件系统。均质等截面铰支杆，刚度值可由材料力学中力

11、与变形的关系中获得均质等截面铰支杆的力位移方程可写为第20页，共28页，2022年，5月20日，3点46分，星期四坐标变换为建立整个结构的刚度矩阵，需要在一个共同的统一坐标系（即总体坐标系）中建立平衡方程。由于刚架各单元的空间位置不同，各个单元的局部坐标系一般也不相同。实际杆件系统都是互相成一定角度排列的杆件连接在一起的每个杆件的单元坐标系统所有杆件的都适用的整体坐标系统第21页，共28页，2022年，5月20日，3点46分，星期四 12对应局部坐标，x，y对应整体坐标系统对应局部坐标系的位移和作用力，对应整体坐标系的位移和作用力。注意：(1)图中 角是从整体坐标系x轴正向起算逆时针转到杆件方

12、向。（2）铰支连接的杆中能承受轴向力 和产生轴向位移 ,因此局部坐标系下 ， 。第22页，共28页，2022年，5月20日，3点46分，星期四方便矩阵运算，将力和位移的矩阵用四阶方程表示：将上式从局部坐标系转换到整体坐标系，表示为：类似地可写出节点2处的表达式。第23页，共28页，2022年，5月20日，3点46分，星期四令 ， ，则节点力的变换关系为（213）或称为变换矩阵。与力的坐标变换式类似，斜杆在两节点的位移有同样的坐标变换式（214）利用式（213）和式（214）可以把局部坐标系下方程（212）表示成整体坐标系下的方程。整体坐标系下单元的刚度矩阵。第24页，共28页，2022年，5月20日，3点46分，星期四用 左乘上式两边（215）再将式（214）代入式（215），有单元刚度矩阵 在整体坐标系下的表达式可以用局部坐标系下的表达式求出， （216）将式（213）代入式（212）有有第25页，共28页，2022年，5月20日，3点46分，星期四求解整体坐标系下结构受力与位移方程组可得到各节点的位移。从而可求出每根杆的受力。i，j整体坐标系中任一杆单元的两个节点号。第26页，共28页，2022年，5月20日，3点46分，星期四二、刚阵存储与节点排列n根杆