弹性力学平面问题的有限元法优秀PPT.ppt

《弹性力学平面问题的有限元法优秀PPT.ppt》由会员分享，可在线阅读，更多相关《弹性力学平面问题的有限元法优秀PPT.ppt（91页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、弹性力学平面问题的有限元法1你现在浏览的是第一页，共91页9.1弹性力学平面问题的基本方程弹性力学平面问题的基本方程弹性力学弹性力学:是研究弹性体（变形体）在约束和外载荷作用下是研究弹性体（变形体）在约束和外载荷作用下应力和变形分布规律的一门学科。应力和变形分布规律的一门学科。方法方法:在弹性力学中针对微小的单元体（在弹性力学中针对微小的单元体（dxdydz）建立基本方）建立基本方程，把复杂形状弹性体的受力和变形分析问题归结为偏微程，把复杂形状弹性体的受力和变形分析问题归结为偏微分方程组的边值问题。分方程组的边值问题。弹性力学的基本方程弹性力学的基本方程:平衡方程、几何方程、物理方程。平衡方程

2、、几何方程、物理方程。有限元方法所处理的对象有限元方法所处理的对象:任意变形体任意变形体变形体变形体:即物体内任意两点之间可发生相对移动。即物体内任意两点之间可发生相对移动。2你现在浏览的是第二页，共91页1）连续，）连续，2）均匀，）均匀，3）各向同性，）各向同性，4）完全弹性，）完全弹性，5）小变形。）小变形。5)小变形假定小变形假定:物体变形远小于物体的几何尺寸，在建立方程时，可以物体变形远小于物体的几何尺寸，在建立方程时，可以略去高阶小量略去高阶小量(二阶以上二阶以上)。1)物体内的物质连续性假定物体内的物质连续性假定:物质无空隙，可用连续函数来描述。物质无空隙，可用连续函数来描述。2

3、)物体内的物质均匀性假定物体内的物质均匀性假定:物体内各个位置的物质具有相同特性。物体内各个位置的物质具有相同特性。3)物体内的物质物体内的物质(力学力学)特性各向同性假定特性各向同性假定:物体内同一位置的物质在各个物体内同一位置的物质在各个方向上具有相同特性。方向上具有相同特性。4)线性弹性假定线性弹性假定:物体的变形与外力作用的关系是线性的，外力去除物体的变形与外力作用的关系是线性的，外力去除后，物体可恢复原状。后，物体可恢复原状。弹性力学的基本假定如下：弹性力学的基本假定如下：3你现在浏览的是第三页，共91页9.1.1基本变量基本变量弹性力学中的基本变量为弹性力学中的基本变量为体力体力、

4、面力面力、应力应力、位移位移、应变应变1）体力：体力：是分布在物体体积内部的力，例如重力和惯性力。是分布在物体体积内部的力，例如重力和惯性力。2）面力：面力：是作用在物体表面上的力，例如两物体间接触力、是作用在物体表面上的力，例如两物体间接触力、流体压力。流体压力。3）应力：应力：物体受到约束和外力作用，其内部将产生内力。物体内物体受到约束和外力作用，其内部将产生内力。物体内某一点的内力就是应力。某一点的内力就是应力。应力应力S在其作用截面上的法向分量称为在其作用截面上的法向分量称为正应力正应力，用，用表示；在作用截面上的表示；在作用截面上的切向分量称为切向分量称为剪应力剪应力，用，用表示。表

5、示。外力外力4你现在浏览的是第四页，共91页外力外力内力内力内力5你现在浏览的是第五页，共91页将每个面上的应力分解为一个正应力和将每个面上的应力分解为一个正应力和两个剪应力，分别与三个坐标轴平行。两个剪应力，分别与三个坐标轴平行。研究受外力作用的物体中某点的应力状态研究受外力作用的物体中某点的应力状态单元体应力分量问题：问题：1）、下标表示？）、下标表示？2）、剪应力互等关系？）、剪应力互等关系？6你现在浏览的是第六页，共91页剪应力互等：剪应力互等：作用在两个互相垂直的面上并且垂直于该作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的（大小、方向）。两面交线的剪应力是互等的（大小

6、、方向）。物体内任意一点的物体内任意一点的应应力状力状态态可以用六个独立的可以用六个独立的应应力分量力分量来表示。来表示。7你现在浏览的是第七页，共91页4）位移位移:包括刚体位移、相对位移。包括刚体位移、相对位移。由于物体受力后发生了变形，物体内个点间的相对移动。用位移在由于物体受力后发生了变形，物体内个点间的相对移动。用位移在x，y，z坐标轴上的投影坐标轴上的投影u、v、w表示。表示。5）应变应变:物体的形状改变可以归结为长度和角度的改变物体的形状改变可以归结为长度和角度的改变。各线段的单位长度的伸缩，称为各线段的单位长度的伸缩，称为正应变正应变，用，用表示。表示。两个垂直线段之间的直角的

7、改变，用弧度表示，称为两个垂直线段之间的直角的改变，用弧度表示，称为剪应变剪应变，用用表示。表示。与与应应力定力定义类义类似，物体内任意一点的似，物体内任意一点的变变形，可以用形，可以用六个六个应变应变分量表示。分量表示。8你现在浏览的是第八页，共91页9.1.3平衡方程（平衡方程（应力应力体力之间关系）体力之间关系）弹性力学中，在物体中取出一个微小单元体建立平衡方程。平衡弹性力学中，在物体中取出一个微小单元体建立平衡方程。平衡方程代表了方程代表了力的平衡力的平衡关系，建立了关系，建立了应力分量应力分量和和体力分量体力分量之间的关之间的关系。对于系。对于平面问题平面问题，在物体内的任意一点有，

8、在物体内的任意一点有，9你现在浏览的是第九页，共91页三维应力情况下的平衡微分方程三维应力情况下的平衡微分方程若：若：则平面问题则平面问题10你现在浏览的是第十页，共91页9.1.4几何方程（应变几何方程（应变位移关系）位移关系）对于平面问题，总的变形可分解为长度变化和角度变化：对于平面问题，总的变形可分解为长度变化和角度变化：由几何方程可以得到由几何方程可以得到位移和变形位移和变形之间的关系。之间的关系。定义定义x方向的相对伸长量方向的相对伸长量定义定义y方向的相对伸长量方向的相对伸长量11你现在浏览的是第十一页，共91页定义夹角的变化定义夹角的变化则定义夹角的总变化为则定义夹角的总变化为则

9、则平面问题平面问题的几何变形方程为的几何变形方程为:12你现在浏览的是第十二页，共91页3D问题的几何变形方程为问题的几何变形方程为:13你现在浏览的是第十三页，共91页变形协调方程（变形连续方程、相容方程）变形协调方程（变形连续方程、相容方程）描述六描述六个应变个应变分量之分量之间的关间的关系。系。14你现在浏览的是第十四页，共91页9.1.5物理方程（应力与应变关系，本构方程）物理方程（应力与应变关系，本构方程）对弹性体，应力对弹性体，应力-应变应变线性关系线性关系广义虎克定理：广义虎克定理：15你现在浏览的是第十五页，共91页虚位移原理虚位移原理若弹性体在已知的面力和体力的作用下处于平衡

10、状态，那么使若弹性体在已知的面力和体力的作用下处于平衡状态，那么使弹性体产生虚位移时，所有作用在弹性体上的外力在虚位移上弹性体产生虚位移时，所有作用在弹性体上的外力在虚位移上所做的功就等于弹性体所具有的虚应变能。所做的功就等于弹性体所具有的虚应变能。同样当虚位移发生时，在弹性体单位体积内应力在相应的虚应变上所作的功同样当虚位移发生时，在弹性体单位体积内应力在相应的虚应变上所作的功为为16你现在浏览的是第十六页，共91页9.1.2平面应力和平面应变问题平面应力和平面应变问题弹性体在满足一定条件时，其弹性体在满足一定条件时，其变形和应力的分布规律可以用在某一平面内变形和应力的分布规律可以用在某一平

11、面内的变形和应力的分布规律来代替的变形和应力的分布规律来代替，这类问题称为，这类问题称为平面问题平面问题。平面应变问题平面应变问题平面问题平面问题平面应力问题平面应力问题1.平面应力问题平面应力问题设设有有很很薄薄的的等等厚厚薄薄板板（某某一一方方向向的的尺尺寸寸较较另另外外两两个个方方向向尺尺寸寸小小很很多多），只只在在板板边边上上受受到到平平行行于于板板面面并并且且不不沿沿厚厚度度变变化化的的面面力力，体力也平行于板面且不沿厚度变化。体力也平行于板面且不沿厚度变化。即即：当结构满足以下两个条件时，则认为是平面应力问题。当结构满足以下两个条件时，则认为是平面应力问题。(1)(1)几何条件：厚

12、度尺寸远远小于截面尺寸，即结构形状呈薄板形。几何条件：厚度尺寸远远小于截面尺寸，即结构形状呈薄板形。(2)(2)载荷条件：载荷平行于板平面且沿厚度方向均匀分布，板平载荷条件：载荷平行于板平面且沿厚度方向均匀分布，板平面不受任何外力作用。面不受任何外力作用。17你现在浏览的是第十七页，共91页设板的厚度为设板的厚度为t，在板面上：，在板面上：由于平板很薄，外力不沿厚度变化，因此在整块板上有，由于平板很薄，外力不沿厚度变化，因此在整块板上有，工程中的许多结构都可作为平面应力问题来处理，如工程中的许多结构都可作为平面应力问题来处理，如链传链传动中的链片动中的链片、发动机中的连杆发动机中的连杆、内燃机

13、的飞轮内燃机的飞轮、轧机的机轧机的机架架和和齿宽较小齿宽较小的直齿圆柱齿轮的直齿圆柱齿轮等。等。18你现在浏览的是第十八页，共91页19你现在浏览的是第十九页，共91页2）平面应变问题）平面应变问题设有很长的柱形体（长度远大于它的横向尺寸），支承情况不沿长度变化，在柱面上设有很长的柱形体（长度远大于它的横向尺寸），支承情况不沿长度变化，在柱面上受到平行于横截面（受到平行于横截面（z轴）而且不沿长度变化的面力，体力也如此分布。轴）而且不沿长度变化的面力，体力也如此分布。凡满足以下两个条件的结构可视为平面应凡满足以下两个条件的结构可视为平面应变问题。变问题。(1)(1)几何条件几何条件:沿厚度方向

14、的截面沿厚度方向的截面形状和大小相同且厚度尺寸远远大于截形状和大小相同且厚度尺寸远远大于截面尺寸，即结构呈等截面的细长形。面尺寸，即结构呈等截面的细长形。(2)(2)载荷条件载荷条件:载荷垂直于厚度方载荷垂直于厚度方向向(平行横截面平行横截面)且沿厚度均匀分布，两个且沿厚度均匀分布，两个端面不受力。端面不受力。20你现在浏览的是第二十页，共91页以柱体的任一横截面为以柱体的任一横截面为XY平面，任一纵线为平面，任一纵线为Z轴。假定该柱体为无限轴。假定该柱体为无限长，则任一截面都可以看作对称面。由对称性，长，则任一截面都可以看作对称面。由对称性，工程中工程中滚针轴承的滚针滚针轴承的滚针、轧钢机的

15、轧辊轧钢机的轧辊、水坝水坝、受内压管道受内压管道、齿宽较大的直齿轮齿宽较大的直齿轮等都可按平面应变等都可按平面应变问题来处理。问题来处理。位移分量都不沿位移分量都不沿z方向变化。方向变化。21你现在浏览的是第二十一页，共91页1 1、平面应力问题中（、平面应力问题中（Z Z轴垂直于该平面），诸应力分量中为零的是（轴垂直于该平面），诸应力分量中为零的是（）。）。2、在平面应力问题中，沿板厚方向（、在平面应力问题中，沿板厚方向（）。）。A应变为零，但应力不为零应变为零，但应力不为零B应力为零，但应变不为零应力为零，但应变不为零C应力、应变都为零应力、应变都为零D应变、应力都不为零应变、应力都不为零

16、3 3、从作图的结构体中取出单元体进行应力状态分析，正确的是、从作图的结构体中取出单元体进行应力状态分析，正确的是()()A.x=y=0,xy0B.xy=yz=0,x=y0C.yz=xz=0,z=0D.x=y0,xy=0Ax,y，zBxy,xz,yzCx,y,xyDz,yz,xz22你现在浏览的是第二十二页，共91页弹性力学平面问题的物理方程由广义虎克定律得到。弹性力学平面问题的物理方程由广义虎克定律得到。1 1）平面应力问题的物理方程）平面应力问题的物理方程23你现在浏览的是第二十三页，共91页令：令：D=C-1D称为称为弹性矩阵弹性矩阵。对称矩阵，与材料性能参数对称矩阵，与材料性能参数E、

17、有关。有关。由应变求应力的由应变求应力的弹性方程弹性方程。24你现在浏览的是第二十四页，共91页2 2）平面应变问题的物理方程）平面应变问题的物理方程25你现在浏览的是第二十五页，共91页在平面应力问题的物理方程中，将在平面应力问题的物理方程中，将E替换为替换为，替换为替换为可以得到平面应变问题的物理方程；可以得到平面应变问题的物理方程；在平面应变问题的物理方程中，将在平面应变问题的物理方程中，将E替换为替换为、替换替换，可以得到平面应力问题的物理方程。可以得到平面应力问题的物理方程。26你现在浏览的是第二十六页，共91页平面问题的基本解法平面问题的基本解法平面问题的未知变量平面问题的未知变量

18、平衡方程平衡方程几何方程几何方程物理方程物理方程平面应力平面应力平面应变平面应变27你现在浏览的是第二十七页，共91页一、应力法：应力作为基本未知量一、应力法：应力作为基本未知量位移位移应变应变应力应力几何方程几何方程物理方程物理方程应力应力应变应变位移位移物理方程物理方程几何方程几何方程求解求解二、位移法：位移分量作为基本未知量二、位移法：位移分量作为基本未知量求解求解弹性力学问题的基本解法弹性力学问题的基本解法外力外力平衡微分方程平衡微分方程力力平衡微分方程平衡微分方程28你现在浏览的是第二十八页，共91页9.2单元位移函数单元位移函数9.2.1平面问题的平面问题的3节点三角形单元节点三角

19、形单元3结结点三角形点三角形单单元元节节点点i、j、m的坐的坐标标分分别为别为节节点位移分别为点位移分别为逆时针方向编码为正逆时针方向编码为正uvP29你现在浏览的是第二十九页，共91页节点节点位移位移单元内位单元内位移分量移分量应变应变应力应力节点力节点力几何方程几何方程物理方程物理方程平衡方程平衡方程K?怎样描述位移的变化规律？怎样描述位移的变化规律？位移模式位移模式30你现在浏览的是第三十页，共91页uvP9.2.1位移模式位移模式由节点位移求内部任一点位移由节点位移求内部任一点位移弹性体内实际的位移分布可以用单元内的位移分布函数来分弹性体内实际的位移分布可以用单元内的位移分布函数来分块

20、近似地表示。块近似地表示。在单元内的位移变化可以假定一个函数来表示，这个函数称在单元内的位移变化可以假定一个函数来表示，这个函数称为为单元位移函数单元位移函数、或、或单元位移模式单元位移模式。（单元位移和节点位移之间的关系）（单元位移和节点位移之间的关系）假设单元内任一点假设单元内任一点P(x,y)的位移的位移u,v为坐标的某种函数为坐标的某种函数u(x,y)、v(x,y)将将3个结点上的坐标和位移分量代入公个结点上的坐标和位移分量代入公式，就可以将六个待定系数用节点坐标式，就可以将六个待定系数用节点坐标和位移分量表示出来。和位移分量表示出来。31你现在浏览的是第三十一页，共91页将水平位移分

21、量和节点坐标公式中的第一式，将水平位移分量和节点坐标公式中的第一式，（下标（下标i，j，m轮换）轮换）其中其中:32你现在浏览的是第三十二页，共91页同理同理:令令（下（下标标i，j，m轮换轮换）Ni称为形态函数称为形态函数N称为形态矩阵称为形态矩阵33你现在浏览的是第三十三页，共91页单元内的位移函数用节点位移表示，可以简写成，单元内的位移函数用节点位移表示，可以简写成，单元的结点位移记为单元的结点位移记为单元内的位移记为单元内的位移记为34你现在浏览的是第三十四页，共91页选择单元位移函数应满足以下条件：选择单元位移函数应满足以下条件：1）位移模式必须反映单元的刚体位移：常数项）位移模式必

22、须反映单元的刚体位移：常数项2）位移模式必须反映单元的常量应变：线性）位移模式必须反映单元的常量应变：线性3）相邻单元在公共边界上的位移连续，单元之间不能重叠，也）相邻单元在公共边界上的位移连续，单元之间不能重叠，也不能脱离。即位移函数在单元之间连续，称为协调性条件。不能脱离。即位移函数在单元之间连续，称为协调性条件。满足满足1）、）、2）条原则的称为完备性单元。）条原则的称为完备性单元。同时满足三条原则的称为完备协调单元。同时满足三条原则的称为完备协调单元。35你现在浏览的是第三十五页，共91页形态函数形态函数N(x,y)具有以下性质：具有以下性质：1）形态函数在单元节点上形态函数的值，具有

23、）形态函数在单元节点上形态函数的值，具有“本点为本点为1，它点，它点为为0”。2）形态函数在单元中任一点，三个形态函数之和为）形态函数在单元中任一点，三个形态函数之和为1。36你现在浏览的是第三十六页，共91页3）三角形单元任意一边上的形函数，仅与该边的两）三角形单元任意一边上的形函数，仅与该边的两端节点坐标有关，而与其他节点坐标无关。端节点坐标有关，而与其他节点坐标无关。例如：在例如：在ij边上（即边上（即i，j节点），有节点），有37你现在浏览的是第三十七页，共91页例题例题1：如图所示等腰三角形单元，求其形态矩阵：如图所示等腰三角形单元，求其形态矩阵N。38你现在浏览的是第三十八页，共9

24、1页三角形三角形积为积为形态函数为形态函数为形形态态矩矩阵为阵为39你现在浏览的是第三十九页，共91页例例2 2、证明：对平面三角形单元形函数存在下列关系、证明：对平面三角形单元形函数存在下列关系40你现在浏览的是第四十页，共91页节点节点位移位移单元内位单元内位移分量移分量应变应变应力应力节点力节点力几何方程几何方程物理方程物理方程平衡方程平衡方程K位移模式位移模式N41你现在浏览的是第四十一页，共91页9.3单元刚度矩阵单元刚度矩阵用节点位移表示单元内部各点位移：用节点位移表示单元内部各点位移：1由节点位移求应变由节点位移求应变42你现在浏览的是第四十二页，共91页记为记为，B矩矩阵阵称称

25、为为几何矩几何矩阵阵。B矩矩阵阵可以表示可以表示为为分分块块矩矩阵阵的形式的形式由几何方程可以得到单元的应变表达式由几何方程可以得到单元的应变表达式43你现在浏览的是第四十三页，共91页关于几何矩阵关于几何矩阵1)几何矩阵中每一元素几何矩阵中每一元素bi的物理意义：单元某一节点有单位位移的物理意义：单元某一节点有单位位移且其它节点位移皆为且其它节点位移皆为0时，所引起单元内部的时，所引起单元内部的应变分布应变分布。2)几何矩阵中各元素为常数（单元位移确定后，节点坐标值为几何矩阵中各元素为常数（单元位移确定后，节点坐标值为定值），因而定值），因而应变分量为常量应变分量为常量，平面三角形单元为常应

26、变单平面三角形单元为常应变单元元。（这是由于采用线性位移函数的结果）。（这是由于采用线性位移函数的结果）。1、推导杆单元的形态函数、形态矩阵、几何矩阵。、推导杆单元的形态函数、形态矩阵、几何矩阵。练习：练习：EAl12x44你现在浏览的是第四十四页，共91页45你现在浏览的是第四十五页，共91页2、已知三角形单元的形态矩阵为、已知三角形单元的形态矩阵为根据形态矩阵求三角形单元的几何矩阵根据形态矩阵求三角形单元的几何矩阵46你现在浏览的是第四十六页，共91页2由应变求应力由应变求应力由物理方程，可以得到单元的应力表达式由物理方程，可以得到单元的应力表达式D称为称为弹性矩阵弹性矩阵，对于平面应力问

27、题，对于平面应力问题对于平面应变问题对于平面应变问题47你现在浏览的是第四十七页，共91页3由节点位移求应力由节点位移求应力定定义义为为应应力矩力矩阵阵。将应力矩阵分块表示为将应力矩阵分块表示为关于应力矩阵关于应力矩阵1)应力矩阵中每一元素的物理意义：单元某一节点有单位位移且其它应力矩阵中每一元素的物理意义：单元某一节点有单位位移且其它节点位移皆为节点位移皆为0时，所引起单元内部的时，所引起单元内部的应力分布应力分布。2)应力矩阵中各元素为常数，因而应力矩阵中各元素为常数，因而应力分量为常量应力分量为常量，平面三角形单元为平面三角形单元为常应力单元常应力单元。48你现在浏览的是第四十八页，共9

28、1页4由应力求节点力由应力求节点力根据虚位移原理根据虚位移原理49你现在浏览的是第四十九页，共91页5由节点位移求节点力由节点位移求节点力单元刚度矩阵单元刚度矩阵50你现在浏览的是第五十页，共91页例题例题2：如图所示等腰三角形单元，求其刚度矩阵：如图所示等腰三角形单元，求其刚度矩阵K，设，设=0。解：解：1）求）求B51你现在浏览的是第五十一页，共91页2）求）求D3）求）求S52你现在浏览的是第五十二页，共91页4）求）求K53你现在浏览的是第五十三页，共91页9.4单元刚度矩阵的性质与物理意义单元刚度矩阵的性质与物理意义（一）单元刚度矩阵的物理意义（一）单元刚度矩阵的物理意义假假设单设单

29、元的元的节节点位移如下：点位移如下：由由，得到，得到节节点力如下：点力如下：Kix,ix表示表示i节点在水平方向产生单位位移节点在水平方向产生单位位移时，在时，在节节点点i的的水平方向水平方向上需要施加的上需要施加的节节点点力。力。Kiy,ix表示表示i节节点在水平方向产生单位位移时，点在水平方向产生单位位移时，在在节节点点i的的垂直方向垂直方向上需要施加的上需要施加的节节点力。点力。因此单元刚度矩阵中每个元素都可以理解为刚度系数，即在结点产生单因此单元刚度矩阵中每个元素都可以理解为刚度系数，即在结点产生单位位移时需要施加的力。位位移时需要施加的力。54你现在浏览的是第五十四页，共91页（二）

30、单元刚度矩阵的性质（二）单元刚度矩阵的性质1)对称性对称性2)奇异性奇异性3)分块性分块性55你现在浏览的是第五十五页，共91页9.5整体刚度矩阵的形成整体刚度矩阵的形成基本方法是基本方法是刚度集成法刚度集成法，由单元刚度矩阵中的元素累加得到整体，由单元刚度矩阵中的元素累加得到整体刚度矩阵中的元素，即整体刚度矩阵是单元刚度矩阵的集成。刚度矩阵中的元素，即整体刚度矩阵是单元刚度矩阵的集成。如何得到整体刚度矩阵？如何得到整体刚度矩阵？刚刚度度集集成成法法即即结结构构中中的的节节点点力力是是相相关关单单元元节节点点力力的的叠叠加加，整整体体刚刚度度矩矩阵阵的的元元素素是是相相关关单单元元的的单单元元

31、刚刚度度矩矩阵阵元元素素的的集集成成。节节点点3在在整整体体刚刚度度矩矩阵阵的的对对应应系系数数，应应该该是是单单元元（1）、（3）、（4）中中对应系数的集成。对应系数的集成。56你现在浏览的是第五十六页，共91页iijijjmjmmim单元单元(1):i，j，m单元的局部码与整体的总码的对应关系：单元的局部码与整体的总码的对应关系：1，2，3单元单元(2):i，j，m2，4，5单元单元(3):i，j，m2，5，3单元单元(4):i，j，m3，5，6ijmijm12312345645657你现在浏览的是第五十七页，共91页123123456456123123232452455325356356

32、58你现在浏览的是第五十八页，共91页kii1kij1kij1kji1kjj1+kii2+kii3kjm1+kim3kij2kim2+kij3kmi1kmj1+kmi3kmm1+kmm3+kii4kij1+kmj3+kij4kij1+kim4kji2kjj2kjm2kmi2+kji3kij1+kjm3+kji4kmj2kmm2+kjj3+kjj4kij1+kjm4kij1+kmi4kij1+kmj4kij1+kmm4123456123456整体刚度矩阵如下所示：整体刚度矩阵如下所示：总码总码59你现在浏览的是第五十九页，共91页刚度矩阵集成的规则刚度矩阵集成的规则确定整体结构节点数确定整体结构

33、节点数n，整体刚度矩阵为，整体刚度矩阵为2nx2n阶矩阵，整体阶矩阵，整体刚度矩阵分块矩阵为刚度矩阵分块矩阵为nxn阶矩阵。阶矩阵。确定局部码与总码的对应关系，确定局部码与总码的对应关系，i、j、m总码位置。总码位置。将各单元刚度矩阵中的每个子矩阵将各单元刚度矩阵中的每个子矩阵Kije放到整体刚度矩放到整体刚度矩阵的对应位置。阵的对应位置。出现在同一位置上的分块矩阵，迭加成整体刚度矩阵中出现在同一位置上的分块矩阵，迭加成整体刚度矩阵中的一个子矩阵。的一个子矩阵。60你现在浏览的是第六十页，共91页9.6整体刚度矩阵的特点与存储方法整体刚度矩阵的特点与存储方法1、整体刚度矩阵的特点：、整体刚度矩

34、阵的特点：对称性对称性稀疏性稀疏性非零系数带形分布非零系数带形分布对称性：对称性：利用对称性，只保存整体矩阵上三角部分的系数即可。利用对称性，只保存整体矩阵上三角部分的系数即可。61你现在浏览的是第六十一页，共91页稀疏性稀疏性:在整体刚度矩阵中，存在大量的零元素，非零元素的个数只占较小的部在整体刚度矩阵中，存在大量的零元素，非零元素的个数只占较小的部分分。这是由结构体存在大量互不相关的节点造成的。这是由结构体存在大量互不相关的节点造成的。62你现在浏览的是第六十二页，共91页整体刚度矩阵的非零元素分布在以对角线为中心的带形区域内，这种矩阵称为整体刚度矩阵的非零元素分布在以对角线为中心的带形区

35、域内，这种矩阵称为带形矩阵带形矩阵。在包括对角线元素的半个带形区域内，每行具有的元素个数叫做在包括对角线元素的半个带形区域内，每行具有的元素个数叫做半带宽半带宽。63你现在浏览的是第六十三页，共91页2D连续体问题总体刚度矩阵的半带宽连续体问题总体刚度矩阵的半带宽:64你现在浏览的是第六十四页，共91页9.7约束条件的处理约束条件的处理位移边界条件在大多数情形下有两类位移边界条件在大多数情形下有两类:第一类第一类:零位移边界零位移边界第二类第二类:给定具体数值的给定具体数值的位移边界位移边界65你现在浏览的是第六十五页，共91页第一类第一类:零位移边界零位移边界方法方法:对角元素置对角元素置“

36、1”法法不变不变不变不变66你现在浏览的是第六十六页，共91页第二类第二类:给定具体数值的位移边界给定具体数值的位移边界方法方法:对角元素乘大数法对角元素乘大数法67你现在浏览的是第六十七页，共91页9.8单元等效节点载荷列阵单元等效节点载荷列阵1、体积分布力、体积分布力68你现在浏览的是第六十八页，共91页2、分布面力、分布面力xqmjiys69你现在浏览的是第六十九页，共91页xqmjiys70你现在浏览的是第七十页，共91页xqmjiys71你现在浏览的是第七十一页，共91页3、集中力、集中力xqmjiyslilj72你现在浏览的是第七十二页，共91页节点载荷列阵的集成节点载荷列阵的集成

37、单元单元(1)123456单元单元(2)123456单元单元(3)123456单元单元(4)123456P73你现在浏览的是第七十三页，共91页结构的平衡方程结构的平衡方程74你现在浏览的是第七十四页，共91页9.9应用应用iijijjmjmmimP1KN/m1、如图所示结构为一平面应力问题离散、如图所示结构为一平面应力问题离散化后的结构图，用有限元法计算节点位化后的结构图，用有限元法计算节点位移、单元应变、单元应力。（移、单元应变、单元应力。（=0,t=1）解：解：载荷载荷P作用在节点作用在节点1上，边界条件：上，边界条件：u1=u2=u4=v4=v5=v6=01)单元分析：单元分析：75你

38、现在浏览的是第七十五页，共91页由于由于=076你现在浏览的是第七十六页，共91页2)整体分析：整体分析：对称77你现在浏览的是第七十七页，共91页对称78你现在浏览的是第七十八页，共91页3)单元应变计算：单元应变计算：79你现在浏览的是第七十九页，共91页4)单元应力计算：单元应力计算：80你现在浏览的是第八十页，共91页5)单元节点力计算：单元节点力计算：81你现在浏览的是第八十一页，共91页2、等腰直角三角形单元、等腰直角三角形单元ABC，AC=BC=10mm，厚度为，厚度为5mm，已知材，已知材料弹性模量料弹性模量E2x105N/mm2，=0。如。如A点沿点沿x方向位移为方向位移为u

39、A=2x10-5mm，B点沿点沿y方向位移为方向位移为vB=10-5mm，而，而C点的位移为点的位移为0。试计算。试计算此单元此单元3个节点的节点力。个节点的节点力。ijm82你现在浏览的是第八十二页，共91页3、如图所示一厚度为、如图所示一厚度为t的等厚度矩形薄板。该矩形薄板一端的等厚度矩形薄板。该矩形薄板一端固定，一端承受均匀拉力固定，一端承受均匀拉力q(吨吨/米米)，长为，长为2米，宽为米，宽为1米，材料米，材料弹性常数为弹性常数为E和和(=1/3)。在不计自重的情况下，试用有限元法。在不计自重的情况下，试用有限元法求其应力分量。求其应力分量。83你现在浏览的是第八十三页，共91页84你

40、现在浏览的是第八十四页，共91页85你现在浏览的是第八十五页，共91页4、如图所示二杆平面桁架，杆长为、如图所示二杆平面桁架，杆长为L，弹性模量为，弹性模量为E，杆截面积，杆截面积为为A，试求（，试求（1）整体刚度矩阵；（）整体刚度矩阵；（2）在）在1、2节点处引入支承条节点处引入支承条件，写出总体平衡方程。件，写出总体平衡方程。123Pxy45086你现在浏览的是第八十六页，共91页5、三角形单元的面积为、三角形单元的面积为1，厚度为，厚度为1，已知三角形单元的形态矩，已知三角形单元的形态矩阵为阵为利用单元的形态矩阵求三角形单元的刚度矩阵。利用单元的形态矩阵求三角形单元的刚度矩阵。87你现在

41、浏览的是第八十七页，共91页1、在一平面桁架中，已知节点、在一平面桁架中，已知节点3处铅直方向位移为零。若用划行划列法引入支承条件，处铅直方向位移为零。若用划行划列法引入支承条件，则应划去总体刚度矩阵中的（则应划去总体刚度矩阵中的（）第第3行和第行和第3列列第第6行和第行和第6列列第第3行和第行和第6列列第第6行和第行和第3列列2、对于每个节点具有三个位移分量的杆单元，两节点局部码为、对于每个节点具有三个位移分量的杆单元，两节点局部码为1，2，总码为，总码为4和和1。其。其单元刚度矩阵中的元素单元刚度矩阵中的元素k32应放入总体刚度矩阵应放入总体刚度矩阵K的（的（）第第3行第行第2列上列上第第

42、4行第行第1列上列上第第9行第行第6列上列上第第12行第行第11列上列上3、在一平面刚架中共有、在一平面刚架中共有9个杆单元，个杆单元，12个节点，则其总体刚度矩阵个节点，则其总体刚度矩阵K是（是（）9阶方阵阶方阵12阶方阵阶方阵36阶方阵阶方阵912阶矩阵阶矩阵4、若把平面应力问题的弹性矩阵改为平面应变问题的弹性矩阵只需将（）、若把平面应力问题的弹性矩阵改为平面应变问题的弹性矩阵只需将（）E换成换成E/(1-2),换成换成/(1-2)E换成换成E/(1-2),换成换成/(1-)E换成换成E/(1-),换成换成/(1-2)E换成换成E/(1-),换成换成/(1-)88你现在浏览的是第八十八页，

43、共91页5、刚架杆单元与平面三角形单元（、刚架杆单元与平面三角形单元（）单元刚度矩阵阶数不同单元刚度矩阵阶数不同局部坐标系的维数不同局部坐标系的维数不同无任何不同无任何不同节点载荷和位移分量数不同节点载荷和位移分量数不同6、图图示平面示平面结结构的构的总总体体刚刚度矩度矩阵阵和和竖带竖带矩矩阵阵K*的元素的元素总总数数分分别别是（是（）。）。400和和200400和和160484和和200484和和1607、材料性质均匀的三节点三角形单元，其内部各点（、材料性质均匀的三节点三角形单元，其内部各点（）应力和应变均不随位置变化应力和应变均不随位置变化应力和应变均随位置变化应力和应变均随位置变化应力

44、不随位置变化，应变随位置变化应力不随位置变化，应变随位置变化应力随位置变化，应变不随位置变化应力随位置变化，应变不随位置变化89你现在浏览的是第八十九页，共91页8、描述位移与应变关系的方程称（、描述位移与应变关系的方程称（）弹性方程弹性方程几何方程几何方程平衡方程平衡方程虚功方程虚功方程9、在以平面刚架中，支承节点、在以平面刚架中，支承节点4的水平方向位移为已知，若用置大数法引的水平方向位移为已知，若用置大数法引入支承条件，则应将总体刚度矩阵中的（入支承条件，则应将总体刚度矩阵中的（）第第4行和第行和第4列上的元素换为大数列上的元素换为大数A第第4行和第行和第4列上的所有元素换为大数列上的所

45、有元素换为大数A第第10行、第行、第10列上的元素换为大数列上的元素换为大数A第第10行、第行、第10列上的所有元素换为大数列上的所有元素换为大数A10、图示的四根杆组成的平面刚架结构，用杆单元进行有限元分析，单元和节点的划、图示的四根杆组成的平面刚架结构，用杆单元进行有限元分析，单元和节点的划分如图示，则总体刚度矩阵的大小为（分如图示，则总体刚度矩阵的大小为（）8x8阶矩阵阶矩阵10 x10阶矩阵阶矩阵12x12阶矩阵阶矩阵16x16阶矩阵阶矩阵90你现在浏览的是第九十页，共91页11、在弹性力学平面刚架问题中，已知相邻节点总码的最大差值为在弹性力学平面刚架问题中，已知相邻节点总码的最大差值为5，则半，则半宽值为宽值为()1018151212、图示平面应力问题的结构中、图示平面应力问题的结构中,单元刚度矩阵（）单元刚度矩阵（）KI=KIII,KII=KIV,但KIKIIKI=KII,KIII=KIV,但KIKIIIKIKIIKIIIKIVKI=KII=KIII=KIV91你现在浏览的是第九十一页，共91页