杆系结构有限元法解析课件

3杆系结构有限元法3杆系结构有限元法1

杆系结构定义：当结构长度尺寸比两个截面方向的尺寸大得多时，这类结构称为杆件。工程中常见得轴、支柱、螺栓、加强肋以及各类型钢等都属于杆件。

是在节点处通过铆接、焊接或用其他方法把若干个杆件连接起来组成一个能共同承担外部载荷的结构。石油工程中的井架、管汇结构等。

杆系结构定义：当结构长度尺寸比两个截面方向的尺寸2杆件结构可分为桁架和刚架两种有限元法对杆系结构离散，通常采用自然离散的形式，也就是把等截面的杆件作为单元。当单元的两端为铰接，杆件内力只有轴力存在，“杆单元”－“桁架”和其他结构采用铰连接的杆称为桁杆。桁杆的连接处可以自由转动，因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的几何因素主要是截面面积，与截面形状无关。当单元两端可以承受弯矩和剪力作用时称为“梁单元”－“刚架”和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动，因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关，而且与截面形状和方位有很大关系。建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。杆件结构可分为桁架和刚架两种3挖掘机桥梁挖掘机桥梁4鸟巢空间立体网架奥运鸟巢的有限元模型鸟巢空间立体网架奥运鸟巢的有限元模型5工程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹簧。弹簧系统力F与弹簧伸长量（位移）之间关系由胡克定律有式中k为弹簧的刚度，是弹簧的固有参数。它对应于力－位移图中F-关系直线的斜率。当k和力F已知时，可由下式求出弹簧伸长量弹簧力－位移间关系(4—1)3-1引言工程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹6

当处理比较复杂的铰支杆系统时，要确定系统在力F的作用下，节点B、C、D和E处的变形。以便计算各杆件的内应力及各杆所受的轴向力，可假设整个杆件系统也具有像式(4-1)中k值一样的刚度，这样在力F的作用下各点的位移就可以用类似式(4-1)的公式计算了，不过．这时的系统刚度应采用一个矩阵来表示，即，同理，各点的位移也应采用一个矩阵来表示，即，再加上矩阵，就构成了称为对应于施加在系统上各节点力的刚度矩阵。F当处理比较复杂的铰支杆系统时，要确定系统在力F的作用下，节7问题：1、复杂结构其刚度矩阵是多少阶的？2、如何求出？3、为什么着重讨论系统的刚度矩阵？系统的整体刚度矩阵－求出所受外力作用下各杆件节点处的位移－计算各杆件的受力和应力问题：8ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹簧的刚度矩阵是2x2阶的。为求出它们，将图2-4所示弹簧系统看作两个简单的系统，然后合成。一、单个弹簧的刚度矩阵3-2弹簧系统的刚度矩阵ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹9由力的平衡有ku1F1aF2aAA‘(a)u2=0ku1=0F1bF2bu2BB‘(b)ku1F1u2F2AA‘BB‘1)只有节点1可以变形,点2固定2)只有节点2可以变形,点1固定3)根据线弹性系统的叠加原理，叠加1)2)两种情况，就得到与原始问题一样的结构，如图（c），叠加结果为：(c)作用于节点1上的合力作用于节点2上的合力刚度矩阵对成、奇异矩阵（2－5）（2－6）由力的平衡有ku1F1aF2aAA‘(a)u2=0ku10二、组合弹簧的刚度矩阵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无位移，有二、组合弹簧的刚度矩阵kakbu1,F1u2,F2u3,F311组合弹簧的刚度矩阵4)合成。对整个系统来说有3个节点，每个节点只有一个方向的位移。因此方程式应用如下形式：利用线弹性系统的叠加原理，找出3×3阶刚度矩阵各元素的表达式节点1处的合力节点2处的合力节点3处的合力对成、奇异矩阵（2－8）组合弹簧的刚度矩阵4)合成。对整个系统来说有3个节点，每个12用同样的方法可以求解具有更多个弹簧的串连系统，推导过程乏味。知道单个弹簧的刚度矩阵－－直接叠加出多个串联系统的总刚度矩阵。用同样的方法可以求解具有更多个弹簧的串连13知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵对整个系统来说有3个节点，将上述方程扩大成3阶方程：整个系统有3个节点（位移），将上述方程扩大成3阶方程，按矩阵相加原理将两式叠加，（2－9）矩阵扩大办法单元数量增多时，相应扩大后的矩阵就相当大，扩大后的非零元素在矩阵的什么位置，概念上就不很清楚了。知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵对整个系统来14按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办法来叠加。以上面两个弹簧系统为例，系统共三个节点，每个节点有一个自由度，因此，该系统总刚度矩阵应该是3×3阶的矩阵。第1个单元的节点号为1和2，则单元刚度矩阵中的元素在总刚度矩阵中应在位置第1行、第2行的第1列，第2列第2个单元的节点号为2和3，则单元刚度矩阵叠加到总刚度矩阵的第2行、第3行的第2列、第3列元素上按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办15三、方程求解（约束条件的引入）由式（2－6）和式（2－8）可知，刚度矩阵是一个奇异阵，即它的行列式的值为零，矩阵的逆不存在。

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

物理概念解释：对整个系统的位移u1、u2和u3，没有加以限制，从而在任何外力的作用下系统会发生刚体运动。u1＝u2＝

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

pa＝ka×（u2-u1）

有限元方法求解弹簧系统受力问题的基本步骤：①形成每个单元的刚度矩阵②各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵③引入约束条件④以节点位移为未知量求解线性代数方程组⑤用每个单元的力－位移关系求得单元力。弹簧1-2受力pa＝ka×（弹簧1－2长度的变化量）

173－3杆件系统的有限元法一、杆单元分析1、局部坐标系下杆单元刚度矩阵上面求解弹簧系统的有限元方法可以直接用力求解受轴向力的杆件系统。均质等截面铰支杆，刚度值可由材料力学中力与变形的关系中获得均质等截面铰支杆的力－位移方程可写为3－3杆件系统的有限元法一、杆单元分析均质等截面铰支杆的182、坐标变换、整体坐标系下的单元刚度矩阵为建立整个结构的刚度矩阵，需要在一个共同的统一坐标系（即总体坐标系）中建立平衡方程。由于刚架各单元的空间位置不同，各个单元的局部坐标系一般也不相同。实际杆件系统都是互相成一定角度排列的杆件连接在一起的每个杆件的单元坐标系统所有杆件的都适用的整体坐标系统2、坐标变换、整体坐标系下的单元刚度矩阵实际杆件系统都是互相19

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

（2－16）将式（2－13）代入式（2－12）有有（2－14）用左乘上式两边（2－15）再将式（2－14）代入23二、整体（平面桁架）分析求解整体坐标系下结构受力与位移方程组可得到各节点的位移。从而可求出每根杆的受力。i，j——整体坐标系中任一杆单元的两个节点号。（2－17）（2－18）二、整体（平面桁架）分析可得到各节点的位移。从而可求出每根杆24有限元方法求解平面桁架系统受力问题的基本步骤：①形成局部坐标系下每个单元的刚度矩阵②由坐标变化矩阵，获得整体坐标系下单元刚度矩阵③各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵④引入约束条件⑤以节点位移为未知量求解线性代数方程组⑥用每个单元的力－位移关系求得单元力。有限元方法求解平面桁架系统受力问题的基本步骤：25例题例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，节点2处受力Fx2、Fy2，所有杆件材料相同，弹性模量为E，截面积均为A，求各杆受力。（单元3和1之间夹角为45°）yxFx21231①②③Fy2单元①＝0°＝1＝1单元②＝90°＝0＝1单元③＝135°将它们代人（2－17），得到单元①例题例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，26单元②单元③单元②单元③27这里e1表示单元①（具有节点1，2）….，将实际值代人，总刚度矩阵为整个系统有6个自由度，整体刚度矩阵是6×6阶的。将上述单元刚阵按节点号叠加到6×6阶矩阵中，就得到整体刚度矩阵。这里e1表示单元①（具有节点1，2）….，将实际值代人28将u2，v2代人原方程，则其它力可表示为只有u2，v2需要求解，因此上述方程可简化为最后各杆所受到的力，由式（2-18）可求出将u2，v2代人原方程，则其它力可表示为只有u2，v2需要求29例题（ANSYS3-1）例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，节点2处受力Fx2、Fy2，所有杆件材料相同，弹性模量为E，截面积均为A，求各杆受力。（单元3和1之间夹角为45°）yxFx21231①②③Fy2单元类型；Link8单元实常数：A＝1材料属性：EX=3.5E10PRXY=0.1667例题（ANSYS3-1）例3-1：平面三杆桁架如下图所示，30ANSYS算例3-2临时墩拼装场地临时墩ANSYS算例3-2临时墩拼装场地临时墩31杆系结构有限元法解析课件32Beam188单元

图5平面模型顶推140m轴力分布图Beam188单元

图5平面模型顶推140m轴力分布图333-4平面刚架有限元法各杆间是刚性固结的，当刚架结构受外力作用时，杆件内不仅有轴力，还有剪力和弯矩，称这类杆件为梁。3-4平面刚架有限元法各杆间是刚性固结的，当刚架结构受外力341、直接刚度法推导梁单元有限元格式L12平面刚架结构——梁单元1、直接刚度法推导梁单元有限元格式L12平面刚架结构——梁35材料力学或结构力学：梁所受弯矩与变形之间的关系，列方程由静力平衡，列方程以上两式写成矩阵形式即为局部坐标系下梁单元刚度矩阵材料力学或结构力学：梁所受弯矩与变形之间的关系，列方程由静力36扩展后的局部坐标系下的梁单元6×6刚度矩阵局部坐标到整体坐标的变换矩阵其中，单元刚度矩阵在整体坐标系下的形式为扩展后的局部坐标系下的梁单元6×6刚度矩阵局部坐标到整体坐标371、力学条件建立单元受力和位移之间的关系式2、局部坐标系下的单元刚度矩阵3、整体坐标系下的单元刚度矩阵4、单元刚度叠加，构成总体刚度矩阵5、引入边界条件，求线性方程6、得到系统各节点处的位移7、进而得到每根梁所受力和力矩1、力学条件建立单元受力和位移之间的关系式382、位移函数－虚功原理推导梁单元有限元计算格式第一步：写出单元的位移、节点力向量局部坐标系下，节点1的位移向量和力向量对节点2也类似，从而梁1-2的节点位移和节点力向量为这些向量每个包含4项，因此单元刚度矩阵应该是4×4阶的。2、位移函数－对节点2也类似，从而梁1-2的节点位移39第二步：选择适当的位移函数选择一个简单函数，用节点上的位移来表示单元上各点的位移。这一位移函数一般情况下可选择多项式。多项式的系数个数应与单元自由度数目相同，使各点的位移可以用节点处的位移所唯一确定。整个单元具有四个自由度，而且只与x坐标有关，由于(3-7)(3-8)对于梁单元，设位移函数写成第二步：选择适当的位移函数整个单元具有四个自由度，而且只与40写成矩阵形式，有x=L处对梁单元，x=0处，第三步：求单元中任一点的位移与节点位移的关系将梁单元两个节点处，对应的位移代入到(3-7)、(3-8)式中，(3-10)求解方程式（3－10）得从而单元上任一点的位移可用节点位移表示：写成矩阵形式，有x=L处对梁单元，x=0处，第三步：求单元中41第四步：求单元应变-单元位移-节点位移间的关系单元内任一点的应变可以通过对该点处的位移微分得到。将代入，写成矩阵形式几何矩阵第四步：求单元应变-单元位移-节点位移间的关系将42第五步：求应力-应变--节点位移间的关系(弹性力学物理方程)弹性矩阵对于梁的弯曲问题，由材料力学知识可知，应力-应变相当于内力矩与曲率关系，近似表达式为：第五步：求应力-应变--节点位移间的关系(弹性力学物理方程43虚位移原理虚位移原理是能量原理在力学特性分析中的一种具体形式。本章将利用这种原理建立单元特性方程。1.虚功与虚应变能：虚位移原理虚位移原理是能量原理在力学特性分析中的一种具体形式44虚位移原理虚位移原理是能量原理在力学特性分析中的一种具体形式。本章将利用这种原理建立单元特性方程。虚位移是指在约束条件允许的范围内弹性体可能发生的任意微小位移。外力要做的虚功，大小为：同时，弹性体内将产生虚应变={δε}，应力在虚应变上所做的虚功是储存在弹性体内的虚应变能，若用δ表示单位体积虚应变能，δW

={δf}T{R}δ={δε

}T{σ}虚位移原理虚位移原理是能量原理在力学特性分析中的一种具体形式45如果在虚位移发生之前弹性体是平衡的，那么虚位移发生时，外力在虚位移上所做的虚功就等于弹性体的虚应变能，即

δW＝δU对于虚位移原理，在虚位移发生过程中，原有的外力、应力、温度及速度应保持不变。外力的形式有集中力、体力和表明力。如果在虚位移发生之前弹性体是平衡对于虚位移原理46第六步：求节点力与节点位移关系虚功原理：系统保持平衡的充要条件是外力在虚位移上所做的功等于内力在相应虚位移上所做的功。系统各节点虚位移向量节点外力在虚位移上所做的虚功任一点处虚位移引起的虚应变为，该处应力为内应力所做的功（单位体积上的应变能）为：虚位移引起的虚应变同样成立，虚功原理，整个体积上功的平衡有：其中，节点虚位移和节点位移都与积分无关第六步：求节点力与节点位移关系系统各节点虚位移向量节点外力47其中，单元刚度矩阵的表达式：具体到梁单元，积分区域是一维的，且从x＝0～L单元刚度矩阵的表达式：其中，单元刚度矩阵的表达式：具体到梁单元，积分区域是一维的，48第七步：求节点位移与应力关系第七步：求节点位移与应力关系491、单元节点位移和节点力列阵的坐标变换任意平面向量V（如平面梁单元的节点位移或节点力列向量），在总体坐标系xoy中的分量为Vx和Vy，ox轴沿逆时针方向旋转α角与轴同向，如图iθiuijθjuj1、单元节点位移和节点力列阵的坐标变换iθiuijθjuj50则有写成矩阵形式为令为平面梁单元的节点力和节点位移列阵的坐标变换矩阵。则有写成矩阵形式为令为平面梁单元的节点力和节点位移列阵的坐标51设总体坐标系中第e个单元的节点力和节点位移列阵分别为局部坐标系为，单元刚度矩阵为2、单元刚度矩阵的坐标变换由坐标转换可知由2、单元刚度矩阵的坐标变换由坐标转换可知由52得令则上式即为单元刚度矩阵的坐标转换式，即为总体坐标系中第e个单元的刚度矩阵，它是一个对称矩阵。得令则上式即为单元刚度矩阵的坐标转换式，即为总体53ANSYS算例3-3ANSYS算例3-3543杆系结构有限元法3杆系结构有限元法55

杆系结构定义：当结构长度尺寸比两个截面方向的尺寸大得多时，这类结构称为杆件。工程中常见得轴、支柱、螺栓、加强肋以及各类型钢等都属于杆件。

是在节点处通过铆接、焊接或用其他方法把若干个杆件连接起来组成一个能共同承担外部载荷的结构。石油工程中的井架、管汇结构等。

杆系结构定义：当结构长度尺寸比两个截面方向的尺寸56杆件结构可分为桁架和刚架两种有限元法对杆系结构离散，通常采用自然离散的形式，也就是把等截面的杆件作为单元。当单元的两端为铰接，杆件内力只有轴力存在，“杆单元”－“桁架”和其他结构采用铰连接的杆称为桁杆。桁杆的连接处可以自由转动，因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的几何因素主要是截面面积，与截面形状无关。当单元两端可以承受弯矩和剪力作用时称为“梁单元”－“刚架”和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动，因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关，而且与截面形状和方位有很大关系。建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。杆件结构可分为桁架和刚架两种57挖掘机桥梁挖掘机桥梁58鸟巢空间立体网架奥运鸟巢的有限元模型鸟巢空间立体网架奥运鸟巢的有限元模型59工程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹簧。弹簧系统力F与弹簧伸长量（位移）之间关系由胡克定律有式中k为弹簧的刚度，是弹簧的固有参数。它对应于力－位移图中F-关系直线的斜率。当k和力F已知时，可由下式求出弹簧伸长量弹簧力－位移间关系(4—1)3-1引言工程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹60

当处理比较复杂的铰支杆系统时，要确定系统在力F的作用下，节点B、C、D和E处的变形。以便计算各杆件的内应力及各杆所受的轴向力，可假设整个杆件系统也具有像式(4-1)中k值一样的刚度，这样在力F的作用下各点的位移就可以用类似式(4-1)的公式计算了，不过．这时的系统刚度应采用一个矩阵来表示，即，同理，各点的位移也应采用一个矩阵来表示，即，再加上矩阵，就构成了称为对应于施加在系统上各节点力的刚度矩阵。F当处理比较复杂的铰支杆系统时，要确定系统在力F的作用下，节61问题：1、复杂结构其刚度矩阵是多少阶的？2、如何求出？3、为什么着重讨论系统的刚度矩阵？系统的整体刚度矩阵－求出所受外力作用下各杆件节点处的位移－计算各杆件的受力和应力问题：62ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹簧的刚度矩阵是2x2阶的。为求出它们，将图2-4所示弹簧系统看作两个简单的系统，然后合成。一、单个弹簧的刚度矩阵3-2弹簧系统的刚度矩阵ku1,F1u2,F2弹簧的作用力向量为位移向量为从而这个弹63由力的平衡有ku1F1aF2aAA‘(a)u2=0ku1=0F1bF2bu2BB‘(b)ku1F1u2F2AA‘BB‘1)只有节点1可以变形,点2固定2)只有节点2可以变形,点1固定3)根据线弹性系统的叠加原理，叠加1)2)两种情况，就得到与原始问题一样的结构，如图（c），叠加结果为：(c)作用于节点1上的合力作用于节点2上的合力刚度矩阵对成、奇异矩阵（2－5）（2－6）由力的平衡有ku1F1aF2aAA‘(a)u2=0ku64二、组合弹簧的刚度矩阵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无位移，有二、组合弹簧的刚度矩阵kakbu1,F1u2,F2u3,F365组合弹簧的刚度矩阵4)合成。对整个系统来说有3个节点，每个节点只有一个方向的位移。因此方程式应用如下形式：利用线弹性系统的叠加原理，找出3×3阶刚度矩阵各元素的表达式节点1处的合力节点2处的合力节点3处的合力对成、奇异矩阵（2－8）组合弹簧的刚度矩阵4)合成。对整个系统来说有3个节点，每个66用同样的方法可以求解具有更多个弹簧的串连系统，推导过程乏味。知道单个弹簧的刚度矩阵－－直接叠加出多个串联系统的总刚度矩阵。用同样的方法可以求解具有更多个弹簧的串连67知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵对整个系统来说有3个节点，将上述方程扩大成3阶方程：整个系统有3个节点（位移），将上述方程扩大成3阶方程，按矩阵相加原理将两式叠加，（2－9）矩阵扩大办法单元数量增多时，相应扩大后的矩阵就相当大，扩大后的非零元素在矩阵的什么位置，概念上就不很清楚了。知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵对整个系统来68按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办法来叠加。以上面两个弹簧系统为例，系统共三个节点，每个节点有一个自由度，因此，该系统总刚度矩阵应该是3×3阶的矩阵。第1个单元的节点号为1和2，则单元刚度矩阵中的元素在总刚度矩阵中应在位置第1行、第2行的第1列，第2列第2个单元的节点号为2和3，则单元刚度矩阵叠加到总刚度矩阵的第2行、第3行的第2列、第3列元素上按节点号将相应单元的刚度矩阵中元素kij写到总刚度矩阵中的办69三、方程求解（约束条件的引入）由式（2－6）和式（2－8）可知，刚度矩阵是一个奇异阵，即它的行列式的值为零，矩阵的逆不存在。

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

物理概念解释：对整个系统的位移u1、u2和u3，没有加以限制，从而在任何外力的作用下系统会发生刚体运动。u1＝u2＝

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

pa＝ka×（u2-u1）

有限元方法求解弹簧系统受力问题的基本步骤：①形成每个单元的刚度矩阵②各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵③引入约束条件④以节点位移为未知量求解线性代数方程组⑤用每个单元的力－位移关系求得单元力。弹簧1-2受力pa＝ka×（弹簧1－2长度的变化量）

713－3杆件系统的有限元法一、杆单元分析1、局部坐标系下杆单元刚度矩阵上面求解弹簧系统的有限元方法可以直接用力求解受轴向力的杆件系统。均质等截面铰支杆，刚度值可由材料力学中力与变形的关系中获得均质等截面铰支杆的力－位移方程可写为3－3杆件系统的有限元法一、杆单元分析均质等截面铰支杆的722、坐标变换、整体坐标系下的单元刚度矩阵为建立整个结构的刚度矩阵，需要在一个共同的统一坐标系（即总体坐标系）中建立平衡方程。由于刚架各单元的空间位置不同，各个单元的局部坐标系一般也不相同。实际杆件系统都是互相成一定角度排列的杆件连接在一起的每个杆件的单元坐标系统所有杆件的都适用的整体坐标系统2、坐标变换、整体坐标系下的单元刚度矩阵实际杆件系统都是互相73

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

（2－16）将式（2－13）代入式（2－12）有有（2－14）用左乘上式两边（2－15）再将式（2－14）代入77二、整体（平面桁架）分析求解整体坐标系下结构受力与位移方程组可得到各节点的位移。从而可求出每根杆的受力。i，j——整体坐标系中任一杆单元的两个节点号。（2－17）（2－18）二、整体（平面桁架）分析可得到各节点的位移。从而可求出每根杆78有限元方法求解平面桁架系统受力问题的基本步骤：①形成局部坐标系下每个单元的刚度矩阵②由坐标变化矩阵，获得整体坐标系下单元刚度矩阵③各个单元的刚度矩阵按节点号叠加成整体系统的刚度矩阵④引入约束条件⑤以节点位移为未知量求解线性代数方程组⑥用每个单元的力－位移关系求得单元力。有限元方法求解平面桁架系统受力问题的基本步骤：79例题例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，节点2处受力Fx2、Fy2，所有杆件材料相同，弹性模量为E，截面积均为A，求各杆受力。（单元3和1之间夹角为45°）yxFx21231①②③Fy2单元①＝0°＝1＝1单元②＝90°＝0＝1单元③＝135°将它们代人（2－17），得到单元①例题例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，80单元②单元③单元②单元③81这里e1表示单元①（具有节点1，2）….，将实际值代人，总刚度矩阵为整个系统有6个自由度，整体刚度矩阵是6×6阶的。将上述单元刚阵按节点号叠加到6×6阶矩阵中，就得到整体刚度矩阵。这里e1表示单元①（具有节点1，2）….，将实际值代人82将u2，v2代人原方程，则其它力可表示为只有u2，v2需要求解，因此上述方程可简化为最后各杆所受到的力，由式（2-18）可求出将u2，v2代人原方程，则其它力可表示为只有u2，v2需要求83例题（ANSYS3-1）例3-1：平面三杆桁架如下图所示，节点1、节点3处固定，节点2处受力Fx2、Fy2，所有杆件材料相同，弹性模量为E，截面积均为A，求各杆受力。（单元3和1之间夹角为45°）yxFx21231①②③Fy2单元类型；Link8单元实常数：A＝1材料属性：EX=3.5E10PRXY=0.1667例题（ANSYS3-1）例3-1：平面三杆桁架如下图所示，84ANSYS算例3-2临时墩拼装场地临时墩ANSYS算例3-2临时墩拼装场地临时墩85杆系结构有限元法解析课件86Beam188单元

图5平面模型顶推140m轴力分布图Beam188单元

图5平面模型顶推140m轴力分布图873-4平面刚架有限元法各杆间是刚性固结的，当刚架结构受外力作用时，杆件内不仅有轴力，还有剪力和弯矩，称这类杆件为梁。3-4平面刚架有限元法各杆间是刚性固结的，当刚架结构受外力881、直接刚度法推导梁单元有限元格式L12平面刚架结构——梁单元1、直接刚度法推导梁单元有限元格式L12平面刚架结构——梁89材料力学或结构力学：梁所受弯矩与变形之间的关系，列方程由静力平衡，列方程以上两式写成矩阵形式即为局部坐标系下梁单元刚度矩阵材料力学或结构力学：梁所受弯矩与变形之间的关系，列方程由静力90扩展后的局部坐标系下的梁单元6×6刚度矩阵局部坐标到整体坐标的变换矩阵其中，单元刚度矩阵在整体坐标系下的形式为扩展后的局部坐标系下的梁单元6×6刚度矩阵局部坐标到整体坐标911、力学条件建立单元受力和位移之间的关系式2、局部坐标系下的单元刚度矩阵3、整体坐标系下的单元刚度矩阵4、单元刚度叠加，构成总体刚度矩阵5、引入边界条件，求线性方程6、得到系统各节点处的位移7、进而得到每根梁所受力和力矩1、力学条件建立单元受力和位移之间的关系式922、位移函数－虚功原理推导梁单元有限元计算格式第一步：写出单元的位移、节点力向量局部坐标系下，节点1的位移向量和力向量对节点2也类似，从而梁1-2的节点位移和节点力向量为这些向量每个包含4项，因此单元刚度矩阵应该是4×4阶的。2、位移函数－对节点2也类似，从而梁1-2的节点位移93第二步：选择适当的位移函数选择一个简单函数，用节点上的位移来表示单元上各点的位移。这一位移函数一般情况下可选择多项式。多项式的系数个数应与单元自由度数目相同，使各点的位移可以用节点处的位移所唯一确定。整个单元具有四个自由度，而且只与x坐标有关，由于(3-7)(3-8)对于梁单元，设位移函数写成第二步：选择适当的位移函数整个单元具有四个自由度，而且只与94写成矩阵形式，有x=L处对梁单元，x=0处，第三步：求单元中任一点的位移与节点位移的关系将梁单元两个节点处，对应的位移代入到(3-7)、(3-8)式中，(3-10)求解方程式（3