桥梁软件应用有限元简介课件

第二章结构分析的有限元法

第二章结构分析的有限元法12.1有限元法发展简况

利用定义在三角形区域上的分片连续函数和最小位能原理St.Venant扭转问题的近似解有限元法的研究现代有限元法飞机结构分析1943Courant应用数学家、物理学家、工程师1960Tumer、Clough第一次成功尝试第一次用三角形单元平面应力问题解答提出了有限单元法的名称各种非线性问题多物理场耦合问题多尺度问题商品化有限元软件20世纪70年代国外几何非线性：因几何变形引起结构刚度改变材料非线性：弹性（超弹和多线性弹性）、粘弹性、非弹性状态非线性：接触问题2.1有限元法发展简况利用定义在三角形区域上的分片连续函22.1有限元法发展简况

固体力学流体力学传热学电磁学学科应用力学计算结构优化计算功能计算技术纯粹数值技术前、后处理技术的高度智能化和与CAD的集成化2.1有限元法发展简况固体力学流体力学学科应用力学计算结32.2有限元法的基本思路及其求解步骤

经典的解析法从连续体的微分方程入手，寻求满足微分方程和定解条件的适合全域的解析解，一旦得到解析解，就可知道域内任意点的解大多数问题，特别是实际问题很难甚至无法用解析法得到问题的解析解在整个求解域上满足控制方程在边界上满足边界条件的场函数寻找很困难有限元法单元节点有限元模型2.2有限元法的基本思路及其求解步骤经典的解析法42.2有限元法的基本思路及其求解步骤

有限元法基本思路抛弃寻找一个满足整个求解域的场函数的思路把求解域划分成有限个四边形单元对每一个单元通过插值的方法，用其节点上的位移建立该单元的位移函数123每个单元都有与其对应的位移函数表达式用全部单元域之和代替整个求解域，用全部单元的位移函数之和代替满足整个求解域的位移函数4对单元进行力学特性分析，建立单元节点力与单元节点位移的关系，并将结构的外载荷等效移植到节点上，再在节点上建立力的平衡方程，求解后得到节点上的位移，继而得到各个单元的应力5以节点位移为未知量，通过求解力的平衡方程获得节点位移，然后按单元计算应力2.2有限元法的基本思路及其求解步骤有限元法基本思52.2有限元法的基本思路及其求解步骤

有限元法求解步骤1离散化将结构（求解域）划分为有限个单元，让全部单元的集合与原结构近似等价划分单元时，二者在几何形体上越逼近越好，特别是在位移和应力急剧变化的地方2选择单元位移函数在有限元法中，需要用单元节点位移通过插值方法建立单元位移函数（单元位移模式），即用单元节点位移来描述单元位移。单元位移函数的合理与否，直接关系到有限元分析的计算精度、效率和收敛性。通常取为多项式形式2.2有限元法的基本思路及其求解步骤有限元法求解步62.2有限元法的基本思路及其求解步骤

3单元特性分析（1）依照应变与位移之间的几何关系，根据所选择的单元位移函数，建立单元应变与单元节点位移之间的关系式。据此式，在求出节点位移后，可以求得单元应变。（2）依照物理关系（胡克定律），建立单元应力与单元节点位移之间的关系式。据此式，在求出节点位移后，可以求得单元应力。（3）根据虚位移原理或最小势能原理，建立单元刚度方程，即单元节点力与单元节点位移之间的关系式。此步骤核心是计算单元刚度矩阵。4外载荷处理将外载荷（体力、面力等）等效移植到节点上。2.2有限元法的基本思路及其求解步骤3单元特性分析72.2有限元法的基本思路及其求解步骤

5建立节点上的力平衡方程按照有限元法的统一格式，形成如下形式的以节点位移为未知量的代数方程组由各个单元的刚度矩阵组装成的总体刚度矩阵待求的节点位移列阵按节点编号顺序形成的节点载荷列阵6处理边界条件、解算节点位移(2.1)按照实际位移边界条件，对式(2.1)进行整理，解之，可得单元节点位移。有了节点位移，即可根据单元特性分析中建立的关系式，求应力、应变、内力等。后处理：对所选应力、应变等，以彩色云图或图表的形式显示计算结果。2.2有限元法的基本思路及其求解步骤5建立节点上的82.3有限元程序的结构简介

对一个题目或一个实际工程问题进行有限元分析，大体上分3个主要步骤有限元建模有限元求解计算结果分析与整理前处理求解器后处理程序结构2.3有限元程序的结构简介对一个题目或一个实际工程92.3有限元程序的结构简介

前处理几何模型的建立定义约束条件网格剖分确定材料参数和载荷有限元建模形成有限元分析所需用的有限元计算数据可视化有限元模型前处理中，可以用图形显示所建立的几何模型、单元网格、约束条件等2.3有限元程序的结构简介前处理几何模型的建立定义10求解器2.3有限元程序的结构简介

有限元程序的核心部分主要完成有限元模型的力学计算，即根据前处理形成的有限元计算数据，完成以下工作：计算单元刚度矩阵计算节点载荷组装总体刚度矩阵将载荷等效简化到节点上形成总体有限元平衡方程求解节点位移计算应力、应变、内力等求解器2.3有限元程序的结构简介有限元程序的核心部112.3有限元程序的结构简介

后处理根据计算者的要求对计算结果进行检查、分析、整理、打印输出等进行数据检索响应量合成绘制变形图、应力图、应变图、曲线图等可视化的方式分析、观察计算结果计算者进行有限元分析的工作量主要体现在前处理和后处理方面2.3有限元程序的结构简介后处理根据计算者的要求对122.4算例

2.4.1平面三角形常应变单元任意区域三角形单元网格剖分示意图典型三角形单元yxo单元内任意点(x,y)的位移坐标x和y的函数建立单元位移函数通过插值方法建立，即用单元的节点位移来表示单元内任意点的位移1.单元位移函数2.4算例2.4.1平面三角形常应变单元任意区域三角132.4算例

典型三角形单元yxo单元位移函数选用坐标x和y的一次多项式待定系数未知量求解(2)(3)得到2.4算例典型三角形单元yxo单元位移函数选用坐标142.4算例

2.4算例152.4算例

代入将求得的得到用单元的节点位移表示的单元位移函数式中是单元形状函数，简称形函数是常数，取决于单元的三个节点坐标返回P242.4算例代入将求得的得到用单元的节点位移表示的单162.4算例

三角形单元的面积A2.4算例三角形单元的面积A17单元位移函数表达式2.4算例

写成矩阵形式简写为其中，表示单元内任意点处位移的单元位移函数列阵为形函数矩阵返回P28单元位移函数表达式2.4算例写成矩阵形式简写为其中182.4算例

形函数的性质在节点上形函数的值是式（6）表示形函数Ni在其自身节点上的值等于1，在其他节点上的值等于0，即12.4算例形函数的性质在节点上形函数的值是式（6）19单元中任意一点上的各个形函数之和等于1，即22.4算例

由单元中任意一点上的各个形函数之和等于1，即22.4算例202.4算例

小结（1）本节的三角形单元，形函数是线性的，为x、y的一次函数；（2）在单元内部和各条单元边上，位移也是线性的，可由两个节点的位移唯一确定；（3）相邻单元的公共节点的节点位移是相等的，因此，能保证相邻单元在公共边界上以及单元内部的位移连续性。2.4算例小结（1）本节的三角形单元，形函数是线性212.4算例

单元位移函数确定后，根据几何方程求得单元内任意点处的应变，即单元应变2.单元应变和单元应力2.4算例单元位移函数确定后，根据几何方程求得单元222.4算例

由形函数对坐标变量求偏导式（8）代入式（7）中，得到返回P282.4算例由形函数对坐标变量求偏导式（8）代入式（23由式（9）得到2.4算例

单元应变矩阵（几何矩阵）分块矩阵参数由单元的节点坐标确定，因此，它们取决于单元形状，当单元的节点坐标确定后，它们都是常量，所以，3节点三角形单元的应变矩阵[B]是常数矩阵由式（9）得到2.4算例单元应变矩阵（几何矩阵）分242.4算例

根据物理方程其中为平面应力（平面应变）问题的弹性矩阵平面应力问题平面应变问题应力矩阵[S]也是常数矩阵单元应力矩阵返回P282.4算例根据物理方程其中为平面应力（平面应变）问25（1）3节点平面三角形单元，应变矩阵和应力矩阵都为常数矩阵；（2）3节点平面三角形单元，各点的应变和应力都是相同的，且是常数，所以3节点三角形单元是常应变单元，也是常应力单元；（3）采用3节点三角形单元时，在应力变化剧烈或应力梯度较大的部位，单元划分应适当加密。2.4算例

小结（1）3节点平面三角形单元，应变矩阵和应力矩阵都为常数矩阵；263.单元刚度矩阵2.4算例

分析结构在载荷作用下产生变形和应力，于是在各单元之间就产生相互作用。实际上，各单元之间的相互作用是通过相邻边界上（即，单元的边，实际是面）的分布力而产生的。按照有限元方法，结构离散化为一个个单元后，单元之间的相互作用就由单元的节点力来实现，即用单元节点力等效代替相邻边界上的相互作用力，这样，节点力就与单元应力相关，而单元应力与节点位移相关，因此，单元节点力与单元节点位移相关。建立单元节点力与单元节点位移之间的关系表示单元节点力单元节点位移3.单元刚度矩阵2.4算例分析结构在载荷作用下产27由（9）式2.4算例

由（12）式把一个单元作为分析对象时，可以把节点力看作外力。单元节点力和单元节点位移之间的关系可由虚位移原理导出在外力作用下，处于平衡状态的变形体，当发生约束允许的任意微小的虚位移时，外力在虚位移上所做的虚功等于整个体积内的应力在虚应变上所做的虚功。推导1令单元的节点虚位移为由P18由(9)式由（9）式2.4算例由（12）式把一个单元作为分析282.4算例

2节点力在虚位移上所做的虚功为3单元应力在虚应变上所做的虚功为单元体积将代入(19)式由于节点虚位移是任意的4建立单元的虚功方程为2.4算例2节点力在虚位移上所做的虚功为3单元应力29由单元的虚功方程2.4算例

任意性相等节点力和节点位移之间的关系式就建立起来了令则(21)式变为单元刚度方程注意：这里，节点力不是结构上的外载荷，而是按虚位移原理把单元边界上的分布力近似等效到单元节点上的一种节点力。节点力在实际结构中是不存在的总结：式（22）、（23）是由三角形常应变单元推导得到的，但是，这两式及其推导过程所基于的原理和方法具有普遍性。原则上说（22）式是位移有限元分析中普遍适用的单元刚度矩阵表达式，对于不同单元，只是其中的具体计算细节不同。由单元的虚功方程2.4算例任意性相等节点力和节点位30三角形常应变单元的刚度矩阵分析2.4算例

一般情况，单元应变矩阵[B]是坐标的函数矩阵。这里，三角形常应变单元，[B]是常数矩阵。如果材料是线性的、匀质的，矩阵[D]也是常数矩阵。单元厚度t是常量,则dV=tdxdy,因此，三角形常应变单元的刚度矩阵可以写成：三角形常应变单元的刚度矩阵分析2.4算例一般情况，31将单元刚度矩阵写成分块形式：2.4算例

将单元刚度矩阵写成分块形式：2.4算例322.4算例

（1）对称性——单元刚度矩阵是对称矩阵；（2）奇异性——单元刚度矩阵是奇异矩阵，它不存在逆矩阵；（3）主元恒正——单元刚度矩阵对角元素的数值恒大于0，可由（25）式看出；（4）单元刚度矩阵的元素具有明确的物理意义，如图所示单元，其刚度矩阵第一列的6个元素ki1(i=1,2,3,4,5,6)的物理意义是，当单元的第1个节点位移(节点i的ui)为1，而其他节点位移全为0时，需要在6个节点位移方向上施加的节点力的大小。单元刚度矩阵的特性2.4算例（1）对称性——单元刚度矩阵是对称矩阵；332.4算例

4.单元等效节点载荷应用虚位移原理进行载荷等效移植：移植后的节点载荷和移植前的载荷在约束允许的任意虚位移上所做的功相等。单元的节点虚位移为单元的虚位移为划分单元时，一般都将作用有集中力的地方划分为节点，集中力即可直接施加到节点上。以下说明体积力和表面力向节点移植：2.4算例4.单元等效节点载荷应用虚位移原理进行342.4算例

(1)体积力等效移植

令单位体积的力为：单位体积力在x轴和y轴方向的分量令单元体积力等效的移植到单元节点上的等效节点载荷为：由虚位移原理，和在虚位移上所做的虚功相等，即单元面积单元厚度(16)式代入(28)式：2.4算例(1)体积力等效移植令单位体积的力为：352.4算例

2.4算例36特殊地，体积力是重力，且重力方向为负y方向，单元的单位体积力是2.4算例

特殊地，体积力是重力，且重力方向为负y方向，单元的单位体积力37(2)表面力等效移植

2.4算例

工程问题中，表面力一般都垂直于其作用面，所以，在有限元法中，要求定义的表面力垂直于其作用面，这样，可以将表面力分解到沿坐标轴方向。令{q}为表面力矢量，则它可以表示为：表面力在x轴和y轴方向的分量令表面力等效的移植到单元节点上的等效节点载荷为：由虚位移原理，和在虚位移上所做的虚功相等，即表面力作用的边界单元的虚位移为代入(33)式得到：(2)表面力等效移植2.4算例工程问题中，表面力382.4算例

2.4算例392.4算例

沿单元边界均匀分布的表面力mi边上有垂直于边界的均匀分布的表面力q,将分解为x和y方向的均布力qx和qy，这样，mi边上的表面力可以表示为：边的边长为2.4算例沿单元边界均匀分布的表面力mi边上有垂直402.4算例

5.总体平衡方程的建立厚度为t的正方形板，左边固定，在其右上角分别作用有x、y方向的集中力F1x，F1y，建立其总体有限元平衡方程。(1)划分单元，所有节点总体编号121234(2)对号入座12324123单元号局部编码整体编码i14j32m23122.4算例5.总体平衡方程的建立厚度为t的正方形41单元刚度矩阵，其分块形式为2.4算例

整体刚度矩阵，其分块形式为单元号局部编码整体编码i14j32m2312单元刚度矩阵，其分块形式为2.4算例整体刚度矩阵，422.4算例

(3)总体平衡方程121234单元号局部编码整体编码i14j32m2312……总体平衡方程整体刚度矩阵整体节点位移整体节点力2.4算例(3)总体平衡方程121234单元号局部432.4算例

1212342.4算例121234446.位移边界条件的处理及总体平衡方程的求解2.4算例

(1)零位移边界条件的处理划行划列法6.位移边界条件的处理及总体平衡方程的求解2.4算例45对角线元素改1法2.4算例

对角线元素改1法2.4算例462.4算例

2.4算例47(2)非零位移边界条件的处理——乘大数法2.4算例

如果与第i个方程对应的节点位移不为0，而是某一个数值可将kii乘以一个大数M，同时将载荷列阵中对应的元素改为(2)非零位移边界条件的处理——乘大数法2.4算例482.4算例

4.习题作业如图所示的平面应力直角三角形单元直角边长分别为a,b,厚度为t，弹性模量为E，泊松比为求该单元的刚度矩阵。2.4算例4.习题作业如图所示的平面应力直角三角49第二章结构分析的有限元法

第二章结构分析的有限元法502.1有限元法发展简况

利用定义在三角形区域上的分片连续函数和最小位能原理St.Venant扭转问题的近似解有限元法的研究现代有限元法飞机结构分析1943Courant应用数学家、物理学家、工程师1960Tumer、Clough第一次成功尝试第一次用三角形单元平面应力问题解答提出了有限单元法的名称各种非线性问题多物理场耦合问题多尺度问题商品化有限元软件20世纪70年代国外几何非线性：因几何变形引起结构刚度改变材料非线性：弹性（超弹和多线性弹性）、粘弹性、非弹性状态非线性：接触问题2.1有限元法发展简况利用定义在三角形区域上的分片连续函512.1有限元法发展简况

固体力学流体力学传热学电磁学学科应用力学计算结构优化计算功能计算技术纯粹数值技术前、后处理技术的高度智能化和与CAD的集成化2.1有限元法发展简况固体力学流体力学学科应用力学计算结522.2有限元法的基本思路及其求解步骤

经典的解析法从连续体的微分方程入手，寻求满足微分方程和定解条件的适合全域的解析解，一旦得到解析解，就可知道域内任意点的解大多数问题，特别是实际问题很难甚至无法用解析法得到问题的解析解在整个求解域上满足控制方程在边界上满足边界条件的场函数寻找很困难有限元法单元节点有限元模型2.2有限元法的基本思路及其求解步骤经典的解析法532.2有限元法的基本思路及其求解步骤

有限元法基本思路抛弃寻找一个满足整个求解域的场函数的思路把求解域划分成有限个四边形单元对每一个单元通过插值的方法，用其节点上的位移建立该单元的位移函数123每个单元都有与其对应的位移函数表达式用全部单元域之和代替整个求解域，用全部单元的位移函数之和代替满足整个求解域的位移函数4对单元进行力学特性分析，建立单元节点力与单元节点位移的关系，并将结构的外载荷等效移植到节点上，再在节点上建立力的平衡方程，求解后得到节点上的位移，继而得到各个单元的应力5以节点位移为未知量，通过求解力的平衡方程获得节点位移，然后按单元计算应力2.2有限元法的基本思路及其求解步骤有限元法基本思542.2有限元法的基本思路及其求解步骤

有限元法求解步骤1离散化将结构（求解域）划分为有限个单元，让全部单元的集合与原结构近似等价划分单元时，二者在几何形体上越逼近越好，特别是在位移和应力急剧变化的地方2选择单元位移函数在有限元法中，需要用单元节点位移通过插值方法建立单元位移函数（单元位移模式），即用单元节点位移来描述单元位移。单元位移函数的合理与否，直接关系到有限元分析的计算精度、效率和收敛性。通常取为多项式形式2.2有限元法的基本思路及其求解步骤有限元法求解步552.2有限元法的基本思路及其求解步骤

3单元特性分析（1）依照应变与位移之间的几何关系，根据所选择的单元位移函数，建立单元应变与单元节点位移之间的关系式。据此式，在求出节点位移后，可以求得单元应变。（2）依照物理关系（胡克定律），建立单元应力与单元节点位移之间的关系式。据此式，在求出节点位移后，可以求得单元应力。（3）根据虚位移原理或最小势能原理，建立单元刚度方程，即单元节点力与单元节点位移之间的关系式。此步骤核心是计算单元刚度矩阵。4外载荷处理将外载荷（体力、面力等）等效移植到节点上。2.2有限元法的基本思路及其求解步骤3单元特性分析562.2有限元法的基本思路及其求解步骤

5建立节点上的力平衡方程按照有限元法的统一格式，形成如下形式的以节点位移为未知量的代数方程组由各个单元的刚度矩阵组装成的总体刚度矩阵待求的节点位移列阵按节点编号顺序形成的节点载荷列阵6处理边界条件、解算节点位移(2.1)按照实际位移边界条件，对式(2.1)进行整理，解之，可得单元节点位移。有了节点位移，即可根据单元特性分析中建立的关系式，求应力、应变、内力等。后处理：对所选应力、应变等，以彩色云图或图表的形式显示计算结果。2.2有限元法的基本思路及其求解步骤5建立节点上的572.3有限元程序的结构简介

对一个题目或一个实际工程问题进行有限元分析，大体上分3个主要步骤有限元建模有限元求解计算结果分析与整理前处理求解器后处理程序结构2.3有限元程序的结构简介对一个题目或一个实际工程582.3有限元程序的结构简介

前处理几何模型的建立定义约束条件网格剖分确定材料参数和载荷有限元建模形成有限元分析所需用的有限元计算数据可视化有限元模型前处理中，可以用图形显示所建立的几何模型、单元网格、约束条件等2.3有限元程序的结构简介前处理几何模型的建立定义59求解器2.3有限元程序的结构简介

有限元程序的核心部分主要完成有限元模型的力学计算，即根据前处理形成的有限元计算数据，完成以下工作：计算单元刚度矩阵计算节点载荷组装总体刚度矩阵将载荷等效简化到节点上形成总体有限元平衡方程求解节点位移计算应力、应变、内力等求解器2.3有限元程序的结构简介有限元程序的核心部602.3有限元程序的结构简介

后处理根据计算者的要求对计算结果进行检查、分析、整理、打印输出等进行数据检索响应量合成绘制变形图、应力图、应变图、曲线图等可视化的方式分析、观察计算结果计算者进行有限元分析的工作量主要体现在前处理和后处理方面2.3有限元程序的结构简介后处理根据计算者的要求对612.4算例

2.4.1平面三角形常应变单元任意区域三角形单元网格剖分示意图典型三角形单元yxo单元内任意点(x,y)的位移坐标x和y的函数建立单元位移函数通过插值方法建立，即用单元的节点位移来表示单元内任意点的位移1.单元位移函数2.4算例2.4.1平面三角形常应变单元任意区域三角622.4算例

典型三角形单元yxo单元位移函数选用坐标x和y的一次多项式待定系数未知量求解(2)(3)得到2.4算例典型三角形单元yxo单元位移函数选用坐标632.4算例

2.4算例642.4算例

代入将求得的得到用单元的节点位移表示的单元位移函数式中是单元形状函数，简称形函数是常数，取决于单元的三个节点坐标返回P242.4算例代入将求得的得到用单元的节点位移表示的单652.4算例

三角形单元的面积A2.4算例三角形单元的面积A66单元位移函数表达式2.4算例

写成矩阵形式简写为其中，表示单元内任意点处位移的单元位移函数列阵为形函数矩阵返回P28单元位移函数表达式2.4算例写成矩阵形式简写为其中672.4算例

形函数的性质在节点上形函数的值是式（6）表示形函数Ni在其自身节点上的值等于1，在其他节点上的值等于0，即12.4算例形函数的性质在节点上形函数的值是式（6）68单元中任意一点上的各个形函数之和等于1，即22.4算例

由单元中任意一点上的各个形函数之和等于1，即22.4算例692.4算例

小结（1）本节的三角形单元，形函数是线性的，为x、y的一次函数；（2）在单元内部和各条单元边上，位移也是线性的，可由两个节点的位移唯一确定；（3）相邻单元的公共节点的节点位移是相等的，因此，能保证相邻单元在公共边界上以及单元内部的位移连续性。2.4算例小结（1）本节的三角形单元，形函数是线性702.4算例

单元位移函数确定后，根据几何方程求得单元内任意点处的应变，即单元应变2.单元应变和单元应力2.4算例单元位移函数确定后，根据几何方程求得单元712.4算例

由形函数对坐标变量求偏导式（8）代入式（7）中，得到返回P282.4算例由形函数对坐标变量求偏导式（8）代入式（72由式（9）得到2.4算例

单元应变矩阵（几何矩阵）分块矩阵参数由单元的节点坐标确定，因此，它们取决于单元形状，当单元的节点坐标确定后，它们都是常量，所以，3节点三角形单元的应变矩阵[B]是常数矩阵由式（9）得到2.4算例单元应变矩阵（几何矩阵）分732.4算例

根据物理方程其中为平面应力（平面应变）问题的弹性矩阵平面应力问题平面应变问题应力矩阵[S]也是常数矩阵单元应力矩阵返回P282.4算例根据物理方程其中为平面应力（平面应变）问74（1）3节点平面三角形单元，应变矩阵和应力矩阵都为常数矩阵；（2）3节点平面三角形单元，各点的应变和应力都是相同的，且是常数，所以3节点三角形单元是常应变单元，也是常应力单元；（3）采用3节点三角形单元时，在应力变化剧烈或应力梯度较大的部位，单元划分应适当加密。2.4算例

小结（1）3节点平面三角形单元，应变矩阵和应力矩阵都为常数矩阵；753.单元刚度矩阵2.4算例

分析结构在载荷作用下产生变形和应力，于是在各单元之间就产生相互作用。实际上，各单元之间的相互作用是通过相邻边界上（即，单元的边，实际是面）的分布力而产生的。按照有限元方法，结构离散化为一个个单元后，单元之间的相互作用就由单元的节点力来实现，即用单元节点力等效代替相邻边界上的相互作用力，这样，节点力就与单元应力相关，而单元应力与节点位移相关，因此，单元节点力与单元节点位移相关。建立单元节点力与单元节点位移之间的关系表示单元节点力单元节点位移3.单元刚度矩阵2.4算例分析结构在载荷作用下产76由（9）式2.4算例

由（12）式把一个单元作为分析对象时，可以把节点力看作外力。单元节点力和单元节点位移之间的关系可由虚位移原理导出在外力作用下，处于平衡状态的变形体，当发生约束允许的任意微小的虚位移时，外力在虚位移上所做的虚功等于整个体积内的应力在虚应变上所做的虚功。推导1令单元的节点虚位移为由P18由(9)式由（9）式2.4算例由（12）式把一个单元作为分析772.4算例

2节点力在虚位移上所做的虚功为3单元应力在虚应变上所做的虚功为单元体积将代入(19)式由于节点虚位移是任意的4建立单元的虚功方程为2.4算例2节点力在虚位移上所做的虚功为3单元应力78由单元的虚功方程2.4算例

任意性相等节点力和节点位移之间的关系式就建立起来了令则(21)式变为单元刚度方程注意：这里，节点力不是结构上的外载荷，而是按虚位移原理把单元边界上的分布力近似等效到单元节点上的一种节点力。节点力在实际结构中是不存在的总结：式（22）、（23）是由三角形常应变单元推导得到的，但是，这两式及其推导过程所基于的原理和方法具有普遍性。原则上说（22）式是位移有限元分析中普遍适用的单元刚度矩阵表达式，对于不同单元，只是其中的具体计算细节不同。由单元的虚功方程2.4算例任意性相等节点力和节点位79三角形常应变单元的刚度矩阵分析2.4算例

一般情况，单元应变矩阵[B]是坐标的函数矩阵。这里，三角形常应变单元，[B]是常数矩阵。如果材料是线性的、匀质的，矩阵[D]也是常数矩阵。单元厚度t是常量,则dV=tdxdy,因此，三角形常应变单元的刚度矩阵可以写成：三角形常应变单元的刚度矩阵分析2.4算例一般情况，80将单元刚度矩阵写成分块形式：2.4算例

将单元刚度矩阵写成分块形式：2.4算例812.4算例

（1）对称性——单元刚度矩阵是对称矩阵；（2）奇异性——单元刚度矩阵是奇异矩阵，它不存在逆矩阵；（3）主元恒正——单元刚度矩阵对角元素的数值恒大于0，可由（25）式看出；（4）单元刚度矩阵的元素具有明确的物理意义，如图所示单元，其刚度矩阵第一列的6个元素ki1(i=1,2,3,4,5,6)的物理意义是，当单元的第1个节点位移(节点i的ui)为1，而其他节点位移全为0时，需要在6个节点位移方向上施加的节点力的大小。单元刚度矩阵的特性2.4算例（1）对称性——单元刚度矩阵是对称矩阵；822.4算例

4.单元等效节点载荷应用虚位移原理进行载荷等效移植：移植后的节点载荷和移植前的载荷在约束允许的任意虚位移上所做的功相等。单元的节点虚位移为单元的虚位移为划分单元时，一般都将作用有集中力的地方划分为节点，集中力即可直接施加到节点上。以下说明体积力和表面力向节点移植：2.4算例4.单元等效节点载荷应用虚位移原理进行832.4算例

(1)体积力等效移植

令单位体积的力为：单位体积力在x轴和y轴方向的分量令单元体积力等效的移植到单元节点上的等效节点载荷为：由虚位移原理，和在虚位移上所做的虚功相