1、第六章 用有限单元法解平面问题第五节第五节 单元的结点力列阵与劲度矩阵单元的结点力列阵与劲度矩阵第四节第四节 单元的应变列阵和应力列阵单元的应变列阵和应力列阵 第三节第三节 单元的位移模式与解答的收敛性单元的位移模式与解答的收敛性第二节第二节 有限单元法的概念有限单元法的概念第一节第一节 基本量及基本方程的矩阵表示基本量及基本方程的矩阵表示概述概述第六节第六节 荷载向结点移置荷载向结点移置 单元的结点荷载列阵单元的结点荷载列阵第六章 用有限单元法解平面问题例题例题第十一节第十一节 应用变分原理导出有限单元法的基本方程应用变分原理导出有限单元法的基本方程第十节第十节 计算实例计算实例第九节第九节
2、 计算成果的整理计算成果的整理第八节第八节 解题的具体步骤解题的具体步骤 单元的划分单元的划分第七节第七节 结构的整体分析结点平衡方程组结构的整体分析结点平衡方程组第六章 用有限单元法解平面问题 对工程问题，力学研究涉及到对工程问题，力学研究涉及到工程简化工程简化、物理模型物理模型和和力学分析力学分析，而解决力学问题的三大支柱为，而解决力学问题的三大支柱为实验手段实验手段，理论理论分析分析和和计算手段计算手段。用计算手段解决力学问题（用计算手段解决力学问题（计算力学）计算力学）是是计算机科学、计算数学和力学学科计算机科学、计算数学和力学学科交叉、交叉、相互渗透的产物。相互渗透的产物。一般认为计
3、算力学始于一般认为计算力学始于有限元方法的出现有限元方法的出现。数值计算方法是计算力学的核心内容数值计算方法是计算力学的核心内容，它是解决工程实，它是解决工程实际力学问题的有效手段，已被学术界和工程界广泛认可作际力学问题的有效手段，已被学术界和工程界广泛认可作为一种力学状态的分析工具。近几十年来数值方法发展迅为一种力学状态的分析工具。近几十年来数值方法发展迅速，相继出现了：速，相继出现了：第六章 用有限单元法解平面问题变分法（变分法（Variational Method）有限差分法（有限差分法（Finite Difference Method,FDM）有限元法（有限元法（Finite Elem
4、ent Method,FEM）边界元法（边界元法（Boundary Element Method,BEM）无限元法（无限元法（Infinite Element Method,IEM）刚体弹簧模型或刚性有限元法（刚体弹簧模型或刚性有限元法（Rigid-Spring Model RBSM）界面应力元模型（界面应力元模型（Interface Stress Element Model,ISEM）离散元法（离散元法（Distinct Element Method,DEM）关键块理论（关键块理论（Key Block Theory,KBT）非连续变形分析非连续变形分析(Discontinuous Defor
5、mation Analysis,DDA)无单元法（无单元法（Meshless Element Free Method）流形方法（流形方法（Manifold Method,MM）广义有限元法（广义有限元法（Generalized Finite Element Method,GFEM）混合数值方法（混合数值方法（Mixed Numerical Method）现有的数值分析方法现有的数值分析方法第六章 用有限单元法解平面问题第六章第六章 用有限单元法解平面问题用有限单元法解平面问题1.有限元法有限元法（Finite Element Method-FEM）FEM2.FEM的特点的特点 概述概述（1 1
6、）具有）具有通用性和灵活性通用性和灵活性。首先将连续体变换为离散化结构，然后再利用首先将连续体变换为离散化结构，然后再利用分片插值技术分片插值技术与与虚功原理虚功原理或变分方法进行求解。或变分方法进行求解。简称简称FEM，是弹性力学的一种是弹性力学的一种近似解法。近似解法。第六章 用有限单元法解平面问题简史3.FEM简史简史 （2 2）对同一类问题，可以编制出）对同一类问题，可以编制出通用程序通用程序，应用计算机进行计算。应用计算机进行计算。（3 3）只要适当加密网格，就可以达到工程）只要适当加密网格，就可以达到工程要求的精度。要求的精度。1943 1943年柯朗年柯朗(德国著名数学家德国著名
7、数学家)第一次提出第一次提出了了FEMFEM的概念。的概念。FEM FEM是上世纪中期才出现，并得到迅速发展是上世纪中期才出现，并得到迅速发展和广泛应用的一种数值解法。和广泛应用的一种数值解法。第六章 用有限单元法解平面问题 1970 1970年后，年后，FEMFEM被引入我国，并很快地得到应用被引入我国，并很快地得到应用和发展。和发展。简史 1956 1956年，特纳等人提出了年，特纳等人提出了FEMFEM。2020世纪世纪5050年代，平面问题的年代，平面问题的FEMFEM建立，并应用建立，并应用于工程问题。于工程问题。1960 1960年提出了年提出了FEMFEM的名称。的名称。20 2
8、0世纪世纪6060年代后，年代后，FEMFEM应用于各种力学问题和应用于各种力学问题和非线性问题，并得到迅速发展。非线性问题，并得到迅速发展。第六章 用有限单元法解平面问题导出方法5.5.本章介绍平面问题的本章介绍平面问题的FEMFEM4.FEMFEM的主要导出方法的主要导出方法 应用应用静力方法静力方法或变分方法导出。或变分方法导出。仅叙述按位移求解的方法。仅叙述按位移求解的方法。且一般都以平面应力问题来表示。且一般都以平面应力问题来表示。第六章 用有限单元法解平面问题6-1 基本量和基本方程的基本量和基本方程的 矩阵表示矩阵表示 本章无特别指明，均表示为本章无特别指明，均表示为平面应力平面
9、应力问题问题的公式。的公式。采用采用矩阵表示矩阵表示,可使公式统一、简洁，可使公式统一、简洁，且便于编制程序。且便于编制程序。第六章 用有限单元法解平面问题。Tyxff)(f。Tyxvyxu),(,),(d。Txyyx)(。Txyyx)(。Tjjiivuvu)(。TjyjxiyixFFFF)(F基本物理量基本物理量：。Tyxff)(f体力体力:基本物理量位移函数位移函数:应变应变:应力应力:结点位移列阵结点位移列阵:结点力列阵结点力列阵:面力面力:第六章 用有限单元法解平面问题 物理方程物理方程:)(bD)(2100010112cED FEM中应用的方程：中应用的方程：)()(ayvxuyvx
10、uT几何方程几何方程:应用的方程其中其中D D为弹性矩阵，对于平面应力问题是为弹性矩阵，对于平面应力问题是:第六章 用有限单元法解平面问题 -结点虚位移结点虚位移;-对应的虚应变。对应的虚应变。ATTdxdytF*)()(*应用的方程图6-1yxoij*,iiyvF*,iixuF*,jjyvF*,jjxuF虚功方程虚功方程:其中其中:在在FEMFEM中，用结点的平衡方程代替平衡中，用结点的平衡方程代替平衡微分方程，后者不再列出。微分方程，后者不再列出。第六章 用有限单元法解平面问题 3 3.整体分析整体分析。6-2 6-2 有限单元法的概念有限单元法的概念 FEMFEM的概念，可以简述为：的概
11、念，可以简述为：采用有限自由度采用有限自由度的离散单元组合体模型去描述实际具有无限自由的离散单元组合体模型去描述实际具有无限自由度的考察体，是一种在力学模型上进行近似的数度的考察体，是一种在力学模型上进行近似的数值计算方法。值计算方法。其理论基础是其理论基础是分片插值技术分片插值技术与与虚功原理虚功原理或变分原理。或变分原理。FEM的概念1.1.将连续体变换为离散化结构；将连续体变换为离散化结构；2.2.单元分析；单元分析；FEMFEM的分析过程：的分析过程：第六章 用有限单元法解平面问题(a)桁架(b)深梁（连续体）结构力学研究的对象结构力学研究的对象是是离散化结构离散化结构。如桁架，。如桁
12、架，各单元（杆件）之间除结点铰结外，没有其他联各单元（杆件）之间除结点铰结外，没有其他联系（图（系（图（a a）。）。结构离散化1.结构离散化结构离散化将连续体变换为离散化结构将连续体变换为离散化结构第六章 用有限单元法解平面问题 将连续体变换为离散化结构将连续体变换为离散化结构:即将连续体划分即将连续体划分为有限多个、有限大小的单元，并使这些单元仅为有限多个、有限大小的单元，并使这些单元仅在一些结点处用绞连结起来，构成所谓在一些结点处用绞连结起来，构成所谓离散化离散化结构结构。结构离散化弹力研究的对象弹力研究的对象，是，是连续体连续体（图（图（b b）)。(a)桁 架(b)深 梁（连 续 体
13、）第六章 用有限单元法解平面问题 图图（c）与图与图(a）相比，两者都是离散化结相比，两者都是离散化结构；区别是，桁架的单元是杆件，而图构；区别是，桁架的单元是杆件，而图（c）的的单元是三角形块体（注意：单元是三角形块体（注意：三角形单元内部仍三角形单元内部仍是连续体是连续体）。）。(c)深 梁（离 散 化 结 构）结构离散化例如例如：将深梁划分为许多三角形单元，这些单将深梁划分为许多三角形单元，这些单元仅在角点用元仅在角点用铰铰连接起来。连接起来。(a)桁 架(b)深 梁（连 续 体）第六章 用有限单元法解平面问题2.2.单元分析单元分析 求解方法 每个三角形单元仍然假定为每个三角形单元仍然
14、假定为连续的、均匀的、连续的、均匀的、各向同性的完全弹性体各向同性的完全弹性体。因单元内。因单元内部仍是连续体，部仍是连续体，应按弹性力学方法进行分析。应按弹性力学方法进行分析。取各结点位移取各结点位移 为基本未为基本未知量知量。然后对每个单元。然后对每个单元,分别求出各物理量分别求出各物理量,并均并均用用 来表示。来表示。),2,1()(ivuTiii),2,1(ii第六章 用有限单元法解平面问题（1）应用插值公式应用插值公式,由单元结点位移由单元结点位移 ，求单元的位移函数，求单元的位移函数Tmjie)(。Tyxvyxu),(),(d求解方法这个插值公式称为单元的这个插值公式称为单元的位移
15、模式位移模式，为：，为：。ed 单元分析的主要内容：单元分析的主要内容：第六章 用有限单元法解平面问题（4 4）应用虚功方程，由单元的应力）应用虚功方程，由单元的应力 ，求出求出单元的结点力单元的结点力，表示为，表示为（3 3）应用物理方程，由单元的应变）应用物理方程，由单元的应变 ，求出求出单元的应力单元的应力，表示为，表示为（2 2）应用几何方程，由单元的位移函数）应用几何方程，由单元的位移函数d d，求出求出单元的应变单元的应变，表示为，表示为。eS。eB求解方法。emjiekFFFF(第六章 用有限单元法解平面问题单元对结点的单元对结点的作用力，与作用力，与 数数值相同值相同,方向相反
16、，方向相反，作用于结点。作用于结点。-结点对单元的作用力，作用结点对单元的作用力，作用 于单元，称为于单元，称为结点力结点力，以，以正标向为正正标向为正。TiyixFF(iFTiyixFF(iF求解方法iFimjxyoiixFiyFjxFjyFmxFmyFiyFixFivmvjviumuju第六章 用有限单元法解平面问题（5 5）将每一单元中的各种外荷载，按）将每一单元中的各种外荷载，按虚功虚功 等效原则等效原则移置到结点上，化为移置到结点上，化为结点荷结点荷 载载，表示为，表示为 .(eLmLjLieLFFFF求解方法第六章 用有限单元法解平面问题 为已知值为已知值,是用结点位移表示的值。是
17、用结点位移表示的值。通过求解联立方程，得出各结点位移值，从而求通过求解联立方程，得出各结点位移值，从而求出各单元的应变和应力。出各单元的应变和应力。各单元移置到各单元移置到i i 结点上的结点荷载结点上的结点荷载 其中其中 表示对围绕表示对围绕i i 结点的单元求和；结点的单元求和；iF求解方法LiF3.3.整体分析整体分析,iF,FLi),2,1(,ieLieiFFe各单元对各单元对i i 结点的结点力结点的结点力作用于结点作用于结点i i上的力有：上的力有：第六章 用有限单元法解平面问题求解方法 3.3.整体分析整体分析 2.2.对单元进行分析对单元进行分析 1.1.将连续体变换为离散化结
18、构将连续体变换为离散化结构归纳起来，归纳起来，FEMFEM分析的主要步骤分析的主要步骤：（1 1）单元的位移模式）单元的位移模式（2 2）单元的应变列阵）单元的应变列阵（4 4）单元的结点力列阵单元的结点力列阵（5 5）单元的等效结点荷载列阵）单元的等效结点荷载列阵建立结点平衡方程组，求解各结点的位移。建立结点平衡方程组，求解各结点的位移。（3 3）单元的应力列阵）单元的应力列阵第六章 用有限单元法解平面问题思考题 1.1.桁架的单元为杆件，而平面体的单元为三角形桁架的单元为杆件，而平面体的单元为三角形块体，在三角形内仍是作为连续体来分析的。块体，在三角形内仍是作为连续体来分析的。前者可用结构
19、力学方法求解，后者只能用弹性前者可用结构力学方法求解，后者只能用弹性力学方法求解，为什么？力学方法求解，为什么？2.2.在平面问题中，是否也可以考虑其它的单在平面问题中，是否也可以考虑其它的单 元形状，如四边形单元？元形状，如四边形单元？第六章 用有限单元法解平面问题应用插值公式，可由应用插值公式，可由 求出位移求出位移 。首先必须解决：首先必须解决：由由单元的结点位移单元的结点位移来求出单元的位移函数来求出单元的位移函数 FEMFEM是取结点位移是取结点位移 为基本未知数。为基本未知数。这个插值公式表示了单元中位移的分布形式，这个插值公式表示了单元中位移的分布形式，因此称为因此称为位移模式位
20、移模式。Tmjie(i。Tyxvyxu),(),(de6-3 单元的位移模式与单元的位移模式与 解答的收敛性解答的收敛性 位移模式d第六章 用有限单元法解平面问题 插值公式（插值公式（a a）在结点在结点 应等于结应等于结点位移值点位移值 。由此可求出。由此可求出 泰勒级数展开式中，低次幂项是最重要的。所泰勒级数展开式中，低次幂项是最重要的。所以以三角形单元的位移模式三角形单元的位移模式，可取为：，可取为：。yxvyxu654321,),(,mjiyxii),(,mjivuii。61三角形单元（a a）第六章 用有限单元法解平面问题将式（将式（a a）按未知数按未知数 归纳为归纳为:其中其中
21、包含包含 。及,iiiivuyx,iivu。mmjjiimmjjiivNvNvNvuNuNuNu,三角形单元61或用矩阵表示为或用矩阵表示为:（b b）第六章 用有限单元法解平面问题N 称为形（态）函数矩阵。称为形（态）函数矩阵。eNdmmjjiimjimjivuvuvuNNNNNNvu000000三角形单元（c c）第六章 用有限单元法解平面问题 A A为三角形为三角形 的面积（图示坐标系中，的面积（图示坐标系中，按逆时针编号），有：按逆时针编号），有：其中其中:),(,2)(mjiAycxbaNiiii11,(,)11jjjjiiimmmmxyyxabci j mxyyx ijmmji,。
22、mmjjiiyxyxyxA1112三角形单元iijjN第六章 用有限单元法解平面问题 三结点三角形单元的位移模式，略去了三结点三角形单元的位移模式，略去了2 2次以次以上的项，因而其上的项，因而其误差量级是误差量级是 且其中只包含且其中只包含了了 的的1 1次项，所以在单元中次项，所以在单元中 的分布如图的分布如图（a a）所示，所示，的分布如图（的分布如图（b b）、（）、（c c）所示。所示。jimjjmmii);(2xo yx,iNvu和三角形单元(a)(b)(c)ivmvjviumuju1第六章 用有限单元法解平面问题 所以当单元趋于很小时，即所以当单元趋于很小时，即 时，为了使时，为
23、了使FEMFEM之解逼近于真解。则为了之解逼近于真解。则为了保保证证FEMFEM收敛性收敛性,位移模式应满足下列条件：位移模式应满足下列条件：FEMFEM中以后的一系列工作，都是以位移中以后的一系列工作，都是以位移模式为基础的。模式为基础的。0,yx收敛性条件第六章 用有限单元法解平面问题 因为当单元因为当单元 时，单元中的位移和时，单元中的位移和应变都趋近于基本量应变都趋近于基本量刚体位移和常量刚体位移和常量位移。位移。（1 1）位移模式必须能反映单元的刚体位移。）位移模式必须能反映单元的刚体位移。0收敛性条件（2 2）位移模式必须能反映单元的常量应变。）位移模式必须能反映单元的常量应变。单
24、元的位移由形变位移和刚体位移组成单元的位移由形变位移和刚体位移组成 单元的应变由变量应变和常量应变组成单元的应变由变量应变和常量应变组成 第六章 用有限单元法解平面问题。xxyvyyxu22,22353564353521,00 xvvyuu收敛性条件可见刚体位移项在式（可见刚体位移项在式（a a）中均已反映。中均已反映。与刚体位移相比，与刚体位移相比，将式（将式（a a）写成写成。yxvyxu654321,第六章 用有限单元法解平面问题（3 3）位移模式应尽可能反映位移的连续性。位移模式应尽可能反映位移的连续性。即应尽可能反映原连续体的位移连续即应尽可能反映原连续体的位移连续性。在三角形单元内
25、部，位移为连续（性。在三角形单元内部，位移为连续（位移位移模式是坐标的单值连续函数模式是坐标的单值连续函数）；在两单元边）；在两单元边界界ijij 上，上，之间均为线性变化，也为之间均为线性变化，也为连续。连续。对式（对式（a a）求应变，得：求应变，得：,5362xyyxji 和收敛性条件可见常量应变也已反映。可见常量应变也已反映。yxvyxu654321,第六章 用有限单元法解平面问题 （1）和（）和（2）是必要条件，而）是必要条件，而加上（加上（3）就为充分条件。）就为充分条件。收敛性条件 为了保证为了保证FEM的收敛性：的收敛性：第六章 用有限单元法解平面问题思考题1.1.应用泰勒级数
26、公式来选取位移模式，为什么必应用泰勒级数公式来选取位移模式，为什么必须从低次项开始选取？须从低次项开始选取？2.2.试考虑：将结构力学解法引入到求解连续体的试考虑：将结构力学解法引入到求解连续体的问题时，位移模式的建立是一个关键性工作，问题时，位移模式的建立是一个关键性工作，它使得单元它使得单元(连续体连续体)内部的分析工作都有可能内部的分析工作都有可能进行了。进行了。第六章 用有限单元法解平面问题6-4 6-4 单元的应变列阵和应力列阵单元的应变列阵和应力列阵 。mmjjiimmjjiivNvNvNvuNuNuNu,),(2/)(mjiAycxbaNiiii。位移函数其中，单元中的位移函数单
27、元中的位移函数用位移模式表示为第六章 用有限单元法解平面问题应用应用几何方程几何方程，求出，求出单元的应变列阵：单元的应变列阵：()00010002TiiijmjijmjiijjmmmmuvvuxyxyuvbbbucccvAcbcbcbuveB。）（a应变第六章 用有限单元法解平面问题)(),(bmjiBBBB)(),(0021cmjibccbAiiii。iB)(,deeSDB D应变S称为应力转换矩阵应力转换矩阵，写成分块形式为再应用物理方程，求出单元的应力列阵：B 称为应变矩阵应变矩阵，用分块矩阵表示，第六章 用有限单元法解平面问题 对于线性位移模式，求导后得到的应变和应对于线性位移模式，
28、求导后得到的应变和应力，均成为常量，因此，称为力，均成为常量，因此，称为常应变（应力）单常应变（应力）单元元。应变和应力的误差量级是。应变和应力的误差量级是 其精度比位其精度比位移低一阶，且相邻单元的应力是跳跃式的。移低一阶，且相邻单元的应力是跳跃式的。)(),(emjiSSSS)(),(2121)1(22fmjibccbcbAEiiiiii。iiDBS),(xo 应力第六章 用有限单元法解平面问题思考题1.1.如果在位移模式中取到泰勒级数中的二如果在位移模式中取到泰勒级数中的二次幂项，略去次幂项，略去 高阶小量，试考虑位移、高阶小量，试考虑位移、应变和应力的误差量级。应变和应力的误差量级。3
29、x第六章 用有限单元法解平面问题6-5 6-5 单元的结点力列阵与劲度矩阵单元的结点力列阵与劲度矩阵 现在来考现在来考虑其中一个单虑其中一个单元：元：模型oyxjmiiixFiyFjxFjyFmxFmyFiyFixF)(在在FEMFEM中，首先将中，首先将连续体变换为离散化连续体变换为离散化结构的模型。结构的模型。第六章 用有限单元法解平面问题（2 2）单元与周围的单元在边界上已没有联）单元与周围的单元在边界上已没有联 系，只在结点系，只在结点 互相联系。互相联系。mji,（1 1）将作用于）将作用于单元上的各种外荷载单元上的各种外荷载，按静，按静 力等效原则移置到结点上去，力等效原则移置到结
30、点上去，化为等化为等 效结点荷载。效结点荷载。故单元内已没有外荷载。故单元内已没有外荷载。第六章 用有限单元法解平面问题假想将单元与结点假想将单元与结点i i 切开，则：切开，则：),(,)(mjiFFTiyixiF),(,)(mjiFFTiyixiF其数值与其数值与 相同，而方向相反。相同，而方向相反。iF结点力以沿正坐标向为正。以沿正坐标向为正。对单元而言，这是作对单元而言，这是作 用于单元上的用于单元上的“外力外力”。结点作用于单元上的力结点作用于单元上的力，称为结点力结点力，单元作用于结点的力，单元作用于结点的力，为：为：第六章 用有限单元法解平面问题;)(TmjieFFFF ijm。
31、Txyyx)(按虚功方程，在虚位移上，外力的虚外力的虚功等于应力的虚功功等于应力的虚功。结点力而其内部有应力作用，考察已与结点切开后的单元 ，则此单元上作用有外力结点力，应用虚功方程，求单元的结点力：第六章 用有限单元法解平面问题 假设发生一组假设发生一组结点虚位移结点虚位移 则单元内则单元内任一点（任一点（x x,y y）的虚位移为的虚位移为 单元单元内内任一点（任一点（x x,y y）的虚应变为的虚应变为 代入代入虚虚功方程：在单元中，功方程：在单元中，外力（结点力外力（结点力 ）在虚在虚位移（结点虚位移位移（结点虚位移 ）上的虚功，等于应）上的虚功，等于应力力 在虚应变在虚应变 上的虚功
32、，上的虚功，即：即：,)(e*,)(e*Nd,)(e*B eF)()(*e)(*虚功方程。ATTedxdytF*e*)()(）（a第六章 用有限单元法解平面问题其中其中 与与 无关，故式无关，故式(a)a)成为成为式式(b b)是是由应力求结点力的一般公式由应力求结点力的一般公式。因为因为 是独立的任意的虚位移，虚是独立的任意的虚位移，虚功方程对任意的功方程对任意的 均应满足，可得出均应满足，可得出,)()()(TTeTeTBB*e)(*yx,)()(。ATeTedxdytBFT*e*TAdxdyteFB.e)(*e)(*代入代入(b)第六章 用有限单元法解平面问题式（式（c c）是是由结点位
33、移求结点力的一般公式，由结点位移求结点力的一般公式，称为单元的劲度矩阵称为单元的劲度矩阵K。元素)66(其中：其中：再将应力公式代入上式，得再将应力公式代入上式，得单元劲度矩阵（c）eeATetdxdykDBBFtdxdyATDBBk（d）第六章 用有限单元法解平面问题对于三角形单元，对于三角形单元，B B 矩阵内均为常数，矩阵内均为常数，有有,tADBBkT）（e 代入代入 B B，D D，得出得出 k k 如书中（如书中（6-376-37）及）及（6-386-38）所示。）所示。第六章 用有限单元法解平面问题（1 1）是是6 66 6的的方阵方阵，中每一个元素都表示中每一个元素都表示单元单
34、元各结点沿坐标方向发生单位位移时所引起的结点力各结点沿坐标方向发生单位位移时所引起的结点力。（2 2）由反力互等定理，）由反力互等定理，所以所以 是是对称矩对称矩阵阵，以对角线为对称轴。，以对角线为对称轴。k,srrsTkkkk单元劲度矩阵单元劲度矩阵k k的性质的性质（3 3）中每一行（或列）的元素之和为零。中每一行（或列）的元素之和为零。（4 4）由（）由（3 3）可导出行列式）可导出行列式 。0k（5 5）的元素与的元素与 单元的形状和方位等单元的形状和方位等有关，但有关，但与单元的大小和刚体的平动以及作与单元的大小和刚体的平动以及作 度转动无关度转动无关。,tEnkk第六章 用有限单元
35、法解平面问题 （书中P.117页），以直角三角形单元为例，计算了应力转换矩阵S和单元劲度矩阵 。从例题中可以看出，将单元边界上的应力向将单元边界上的应力向结点移置，化为作用于结点上的力，正好就是结点移置，化为作用于结点上的力，正好就是结点力。结点力。在在FEM中，单元边界之间的联系和相中，单元边界之间的联系和相互作用力，都向结点简化，归结成为结点的铰互作用力，都向结点简化，归结成为结点的铰结和结点力。结和结点力。思考题例题k试求出书中例题的位移模式。第六章 用有限单元法解平面问题。TLmyLmxLjyLjxLiyLixTeFFFFFF)()(LmLjLiLFFFF6 66 6荷载向结点移置荷载
36、向结点移置 单元的结点荷载列阵单元的结点荷载列阵 在FEM中，须将作用于单元中的外荷载向结点移置，化为等效结点荷载等效结点荷载，第六章 用有限单元法解平面问题（2）变形体静力等效原则变形体静力等效原则在任意的虚位移上，使原荷载与移置荷载的虚功相等。1 1、等效原则等效原则（1）刚体静力等效原则刚体静力等效原则使原荷载与移置荷载的主矢量以及对同一点的主矩也相同。移置原则 刚体静力等效原则只从运动效应来考虑，得出移置荷载不是唯一的解；变形体的静力等效原则考虑了变形效应，在一定的位移模式下，其结果是唯一的，且也满足了前者条件的。所以在FEM中，采用变形体的静力等效原则变形体的静力等效原则。第六章 用
37、有限单元法解平面问题 假设发生一组结点虚位移 ，则点的虚位移为 。使移置荷载的虚功等于原荷载的虚功：2 2、集中力的移置公式集中力的移置公式 原荷载 作用于单元中任一点 为单位厚度上的作用力；移置荷载 作用于结点 ,)(TPyPxffPf,)(TLmLjLieLFFFF。mji,e)(*Nd。ttTeTeTePT*P*L*fNfdF)()()(e)(*),(yx集中力第六章 用有限单元法解平面问题 对于任意的虚位移 ，虚功方程都必须满足，得：。tePTLfNFe)(*)(a面力?eLF第六章 用有限单元法解平面问题3、单元边界单元边界 上面力上面力 的移置公式的移置公式 应用式 ，将 代之为
38、并在边界 上积分，得:。tePTLfNFSftPf,dstfs。SedstfNFTL)(a)(b)(a面力?eLF?eLF第六章 用有限单元法解平面问题 应用式 ，将 代之为 并对单元域A 积分，得。AedxdytfNFTL,dxdytf)(c)(atPf4 4、单元内体力、单元内体力 的移置公式的移置公式 f体力。tePTLfNF)(a?eLF第六章 用有限单元法解平面问题思考题1.试导出书中例题的荷载移置公式。当位移模式为当位移模式为线性函数线性函数时，由时，由虚功方虚功方程程得出的移置荷载，与按得出的移置荷载，与按刚体静力等效原刚体静力等效原则则得出的结点荷载相同。得出的结点荷载相同。第
39、六章 用有限单元法解平面问题 在单元分析中，从单元的结点位移求位移分布求应变求应力求结点力，为单元的内力分析；外荷载移置到结点荷载，为单元的外力分析。6 67 7结构的整体分析结构的整体分析 结点平衡方程组结点平衡方程组 iFLiF 假设将结点i与周围的单元切开，则围绕i结点的每个单元对i 结点有结点力(）的作用,也有外荷载移置的结点荷载（）的作用。下面考虑整体分析整体分析。第六章 用有限单元法解平面问题对某一个单元 ，其中 是对围绕i 结点的单元求和。i 结点的平衡条件结点的平衡条件为 )(a结点平衡条件),2,1(,nieLieiFFeijm,mjinninikF)(a代入式 ，可表示为
40、),2,1()(,nieLiemjinnin。Fk)(b 是单元结点的局部编号；是整体结点的整体编号。mji,ni,2,1第六章 用有限单元法解平面问题),2,1()(,nieLiemjinnin。Fk)(b)(b将式 按整体结点编号排列，得整个结构的平衡方程组。考虑结构的约束条件后，从式考虑结构的约束条件后，从式 求出求出 ，就可以求出各单元的位移和应力。，就可以求出各单元的位移和应力。整体结点位移列阵，整体结点位移列阵，整体结点荷载列阵，整体结点荷载列阵，整体劲度矩阵。整体劲度矩阵。,LFK Tn21)(TLnLLL)(21FFFFK)(c)(c结点平衡方程组第六章 用有限单元法解平面问题
41、例2例1列出图示结构i 结点的平衡条件。（见书中P.121）psjmi第六章 用有限单元法解平面问题有限单元法的具体计算步骤：有限单元法的具体计算步骤：6 68 8解题的具体步骤解题的具体步骤 单元的划分单元的划分 1、划分单元网格，对单元和结点编号。2、选定直角坐标系，按程序要求填写和输入有关信息。单元内的ijm的局部编号应按书中规定的右手规则编号。否则会使三角形的面积出现负号等问题。第六章 用有限单元法解平面问题 3、使用已编好的程序进行上机计算。事先须将有限单元法的公式，计算方法和步骤都编入程序。4、对成果进行整理、分析。对第1和第4步的工作，也尽可能让计算机执行，以减少人工的工作量。如
42、自动划分网格，整理成果等。第六章 用有限单元法解平面问题 关于单元的划分，注意几点：单元的划分，注意几点：（8）结构具有凹槽或孔洞等应力集中处等。（1）单元大小问题；（2）单元在不同部位的合理布置问题；（3）三角形三个内角最好较接近；（4）利用对称性和反对称性；（5）厚度突变之处和材料不同之处；（6）载荷作用（集中力或突变分布载荷）处；（7）水利闸坝工程问题；第六章 用有限单元法解平面问题三峡升船机上闸首三峡升船机上闸首锦屏拱坝网格图锦屏拱坝网格图 龙滩地下洞室网格龙滩地下洞室网格 400米高尾矿坝有限米高尾矿坝有限元网格元网格第六章 用有限单元法解平面问题反应堆有限元模型反应堆有限元模型第六
43、章 用有限单元法解平面问题第六章 用有限单元法解平面问题第六章 用有限单元法解平面问题 在有限单元法中，位移的精度较高，其误差量级是，即与单元尺度的二次幂成正比。应力的误差量级是，即与单元的大小成正比。6 69 9计算成果的整理计算成果的整理 )(2xo)(xo 第六章 用有限单元法解平面问题 三结点三角形单元的应力的成果，不但应力的精度较低，而且还产生了所谓应应力的波动性力的波动性。对于结点位移的成果，可以直接采用。第六章 用有限单元法解平面问题 应力的波动性在三结点三角形单元中较为显应力的波动性在三结点三角形单元中较为显著著。由于计算出的应力的精度较低。假由于计算出的应力的精度较低。假设设
44、单元的应力成果为单元的应力成果为 ，其中，其中 为真解，为真解，为误差。则由于在结点都列出了平衡方程并令为误差。则由于在结点都列出了平衡方程并令其满足，从而使相邻的其满足，从而使相邻的单元的应力趋近于单元的应力趋近于 。这就产生了应力的波动性。这就产生了应力的波动性。原因是，第六章 用有限单元法解平面问题 为了提高应力的精度，解决应力波动性问题，可以采用两种应力成果的整理方法：一般地讲，两相邻单元平均法的精度较好，因为它涉及的区域范围较小。（1）两相邻单元平均法。（2）绕结点平均法。在受面力边界线附近，求得的应力误差较在受面力边界线附近，求得的应力误差较大。可采用向外插值的方法（例抛物线插大。
45、可采用向外插值的方法（例抛物线插值）来解决。值）来解决。第六章 用有限单元法解平面问题 为了提高应力的精度，可以采用两种方法。是加密网格，减少单元的尺寸，以提高应力的精度。是可以采用较多结点的单元，并使 位移模式中包含一些高幂次的项，从而提 高位移和应力的精度。二一第六章 用有限单元法解平面问题锦屏拱坝锦屏拱坝第六章 用有限单元法解平面问题第六章 用有限单元法解平面问题 书中应用三结点三角形单元，计算了下列例题：6 61010计算实例计算实例 1.楔形体受自重及齐顶水压力。2.简支梁受均布荷载。3.圆孔附近的应力集中。第六章 用有限单元法解平面问题 在整理应力成果时，读者应注意，应用三角形单元
46、时，（1）采用两单元平均法和绕结点平均法的 应力成果比较接近，但前者的精度略 好于后者。（2）边界面的应力，宜采用向外插值的方 法求出。第六章 用有限单元法解平面问题 在FEMFEM中，将连续体变换为离散化结构之后，有两种导出FEM公式的主要方法：6 61111应用变分原理导出应用变分原理导出 有限单元法基本方程有限单元法基本方程 第六章 用有限单元法解平面问题（2）建立单元的位移模式，求出单元中的 位移分布，;eB）（b;eS）（ci;eNd）（a1.按静力方法导出按静力方法导出FEMFEM公式公式（1）取结点位移为基本未知数；（3）由几何方程求出单元的应变，（4）由物理方程求出单元的应力，
47、按结构力学方法导出FEM公式第六章 用有限单元法解平面问题;AsedxdytdsttfNfNfNFTTPeL）（e。LFk ）（f;eekF）（deLF（5）由虚功方程求出单元的结点力，（6）由虚功方程求出单元的结点荷载 ，（7）建立结点平衡方程组，按结构力学方法导出FEM公式第六章 用有限单元法解平面问题（1）变分原理中的极小势能原理是。minVUEP）（g)().(,21hdxdytdsttVdxdytUATsTTAfdfdfdPT2.按变分方法导出按变分方法导出FEMFEM公式公式 保留上述（1）-（4）步骤，然后应用极小势能原理导出FEM基本方程。按变分法导出FEM公式对于平面问题，第
48、六章 用有限单元法解平面问题对于连续体,变分的宗量是位移函数 变分方程 可表示为总势能 对 的导数等于0，即PEvu,.,vu）（g,0uEP。0vEP）（i第六章 用有限单元法解平面问题变分宗量由 变换成（2）将经典变分原理应用到离散化结构，则).,2,1(nii,eePPEE,eeUU。eeVV）（jvu,总势能、形变势能和外力势能，可以用单元的势能之和来表示第六章 用有限单元法解平面问题其中 为三角形单元的面积。应用前面记号,)()(21)(2121 eAeTTeeAeTeeAeeeeedxdytdxdytdxdytUUDBBDB TeA。eeTeUk)(21）（k内力势能为第六章 用有
49、限单元法解平面问题其中 为三角形单元的受面力边界。引用前面记号,)()()(eATSPTPTTeeATSPTPTeeeedxdytdsttdxdytdsttVVfdfdfNfdfdfds。eeLTeVF)(）（leeLTeTePVUE)()(21Fk）（m外力势能为 总势能为第六章 用有限单元法解平面问题故总势能极小值条件 变换为（3）对于离散化结构，泛函数 的宗量变 换为 PE)(i)(),2,1(0nniEiP。.)(,2)(ccaababaaaTT),2,1(nii则式(n)成为引用矩阵运算公式，第六章 用有限单元法解平面问题)(),2,1(0)(oniEieTeP。,001,)(iee
50、LmLjLimjieeLeePEFFFFFFFk 其中第六章 用有限单元法解平面问题 代入式(o)，得出与结构力学方法导出的相同方程，)(),2,1(pnieLiei。FF 从物理意义上讲，将连续体的经典变分原 理(g)或(i)应用到离散化结构，成为式(p)。第六章 用有限单元法解平面问题 比较物理意义：凡是与微分方程对应的变分原理存在的任何问题，均可应用变分法导出FEM。式(p)表示总势能在所有结点处的极值条件。式(g)表示总势能的整体极值条件;第六章 用有限单元法解平面问题例题1例题2例题3例题4例题第六章 用有限单元法解平面问题 例题1 平面问题中采用的四结点矩阵单元，如图所示。该单元的
