《第三讲平面问题有限元分析.ppt》由会员分享，可在线阅读，更多相关《第三讲平面问题有限元分析.ppt（72页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第三讲平面问题有限元分析现在学习的是第1页，共72页第三章第三章 平面问题的有限单元法平面问题的有限单元法3-23-2、有限单元法的计算步骤、有限单元法的计算步骤3-33-3、单元位移函数、单元位移函数3-43-4、单元载荷移置、单元载荷移置3-53-5、单元应力矩阵、单元应力矩阵3-63-6、单元刚度矩阵、单元刚度矩阵3-73-7、单元刚度矩阵的物理意义及其性质、单元刚度矩阵的物理意义及其性质3-83-8、整体分析、整体分析3-93-9、整体刚度矩阵的形成、整体刚度矩阵的形成3-103-10、支承条件的处理、支承条件的处理3-113-11、整体刚度矩阵的特点、整体刚度矩阵的特点现在学习的是第
2、2页，共72页2-2 有限单元法的计算步骤有限单元法的计算步骤弹性力学平面问题的有限单元法包括三个主要步骤：弹性力学平面问题的有限单元法包括三个主要步骤：1、离散化、离散化 2、单元分析、单元分析 3、单元综合、单元综合现在学习的是第3页，共72页2-2 有限单元法的计算步骤有限单元法的计算步骤 1、离散化、离散化 有限单元法的基础是用所谓有限个单元的集合体来有限单元法的基础是用所谓有限个单元的集合体来代替原来的连续体，因而必须将连续体简化为由有限代替原来的连续体，因而必须将连续体简化为由有限个单元组成的离散体。对于平面问题，最简单，因而个单元组成的离散体。对于平面问题，最简单，因而最常用的单
3、元是三角形单元。这些单元在结点处用铰最常用的单元是三角形单元。这些单元在结点处用铰相连，荷载也移置到结点上，成为结点荷载。在结点相连，荷载也移置到结点上，成为结点荷载。在结点位移或其某一分量可以不计之处，就在结点上安置一位移或其某一分量可以不计之处，就在结点上安置一个铰支座或相应的连杆支座。个铰支座或相应的连杆支座。现在学习的是第4页，共72页5.5.相邻单元的尺寸尽可能接近。相邻单元的尺寸尽可能接近。6.6.结点所连接的单元个数尽可能一致。结点所连接的单元个数尽可能一致。宜宜不宜不宜结点的选择和单元划分原则结点的选择和单元划分原则1.1.集中力作用点、分布力突变点、支承点应选作结点。集中力作
4、用点、分布力突变点、支承点应选作结点。2.2.不同厚度、不同材料的部分不应划在同一个单元。不同厚度、不同材料的部分不应划在同一个单元。3.3.应力变化大处单元应密集一些。结点的多少与疏密要考虑计算应力变化大处单元应密集一些。结点的多少与疏密要考虑计算 机的容量和计算精度。机的容量和计算精度。4.4.单元边界的边长之比应尽可能靠近单元边界的边长之比应尽可能靠近1 1。宜宜不宜不宜有限元分析应注意的问题有限元分析应注意的问题现在学习的是第5页，共72页7 7、充分利用结构的对称性、充分利用结构的对称性PPPPP现在学习的是第6页，共72页2-2 有限单元法的计算步骤有限单元法的计算步骤 2、单元分
5、析、单元分析 对三角形单元，建立结点位移与结点力之间的转换关系。对三角形单元，建立结点位移与结点力之间的转换关系。结点位移结点位移 结点力结点力 现在学习的是第7页，共72页2-2 有限单元法的计算步骤有限单元法的计算步骤 2、单元分析、单元分析-单元刚度矩阵单元刚度矩阵 取结点位移作基本未知量。由结点位移求结点力：取结点位移作基本未知量。由结点位移求结点力：其中，转换矩阵称为单元刚度矩阵。单元分析的主要目的就其中，转换矩阵称为单元刚度矩阵。单元分析的主要目的就是要求出单元刚度矩阵。是要求出单元刚度矩阵。单元分析的步骤可表示如下：单元分析的步骤可表示如下：现在学习的是第8页，共72页2-2 有
6、限单元法的计算步骤有限单元法的计算步骤 3、单元综合、单元综合 将离散化了的各个单元合成整体结构，利用结点平衡方程求将离散化了的各个单元合成整体结构，利用结点平衡方程求出结点位移。出结点位移。在位移法中，主要的任务是求出基本未知量在位移法中，主要的任务是求出基本未知量-结点位移。为此需要结点位移。为此需要建立结点的平衡方程。建立结点的平衡方程。现在学习的是第9页，共72页2-2 有限单元法的计算步骤有限单元法的计算步骤 3、单元综合、单元综合 i点总的结点力应为：点总的结点力应为：根据结点的平衡条件，得根据结点的平衡条件，得 单元单元e e的结点力，可按式的结点力，可按式(2-2)(2-2)用
7、结点位移表示，代入得到用用结点位移表示，代入得到用结点位移表示的平衡方程。结点位移表示的平衡方程。每个可动结点有两个未知位移，有两个平衡方程，所以方程每个可动结点有两个未知位移，有两个平衡方程，所以方程总数与未知位移总数相等，可以求出所有的结点位移。总数与未知位移总数相等，可以求出所有的结点位移。单元综合的目的就是要求出结点位移。结点位移求出后，可进一单元综合的目的就是要求出结点位移。结点位移求出后，可进一步求出各单元的应力。步求出各单元的应力。现在学习的是第10页，共72页2-3 单元位移函数单元位移函数 如果弹性体的位移分量是座标的已知函数，则可用几何方程求应变如果弹性体的位移分量是座标的
8、已知函数，则可用几何方程求应变分量，再从物理方程求应力分量。但对一个连续体，内部各点的位移变化分量，再从物理方程求应力分量。但对一个连续体，内部各点的位移变化情况很难用一个简单函数来描绘。情况很难用一个简单函数来描绘。有限单元法的基本原理是分块近似，即将弹性体划分成若干细小网格，有限单元法的基本原理是分块近似，即将弹性体划分成若干细小网格，在每一个单元范围内，内部各点的位移变化情况可近似地用简单函数来描在每一个单元范围内，内部各点的位移变化情况可近似地用简单函数来描绘。对每个单元，可以假定一个简单函数，用它近似表示该单元的位移。绘。对每个单元，可以假定一个简单函数，用它近似表示该单元的位移。这
9、个函数称为这个函数称为位移函数，或称为位移模式、位移模型、位移场位移函数，或称为位移模式、位移模型、位移场。对于平面问题，单元位移函数可以用多项式表示，对于平面问题，单元位移函数可以用多项式表示，多项式中包含的项数越多，就越接近实际的位移分布，越精确。但选取多多项式中包含的项数越多，就越接近实际的位移分布，越精确。但选取多少项数，要受单元型式的限制。少项数，要受单元型式的限制。现在学习的是第11页，共72页2-2 单元位移函数单元位移函数 三结点三角形单元三结点三角形单元六六个个节节点点位位移移只只能能确确定定六六个个多多项项式式的的系系数数，所所以以平平面面问问题题的的3 3结结点点三角形单
10、元的位移函数如下，三角形单元的位移函数如下，所所选选用用的的这这个个位位移移函函数数，将将单单元元内内部部任任一一点点的的位位移移定定为为座座标标的的线线性性函函数，位移模式很简单。数，位移模式很简单。位移函数写成矩阵形式为：位移函数写成矩阵形式为：现在学习的是第12页，共72页将水平位移分量和结点坐标代入位移函数第一式，写成矩阵形式，现在学习的是第13页，共72页令 则有 A为三角形单元的面积。现在学习的是第14页，共72页T的伴随矩阵为，令 现在学习的是第15页，共72页则同样，将垂直位移分量与结点坐标代入位移插值公式现在学习的是第16页，共72页2-3 单元位移函数单元位移函数最终确定六
11、个待定系数最终确定六个待定系数现在学习的是第17页，共72页2-3 单元位移函数单元位移函数令令 （下标（下标i i，j j，m m轮换）轮换）简写为简写为II是单位矩阵，是单位矩阵，NN称为形态矩阵，称为形态矩阵，N Ni i称为位移的形态函数称为位移的形态函数现在学习的是第18页，共72页2-3 单元位移函数单元位移函数 选选择择单单元元位位移移函函数数时时，应应当当保保证证有有限限元元法法解解答答的的收收敛敛性性，即即当当网网格格逐逐渐渐加加密密时时，有有限限元元法法的的解解答答应应当当收收敛敛于于问问题题的的正正确确解解答答。因因此此，选用的位移模式应当满足下列两方面的条件：选用的位移
12、模式应当满足下列两方面的条件：(1)(1)必须能反映单元的刚体位移和常量应变。必须能反映单元的刚体位移和常量应变。6 6个参数个参数 到到 反映了三个刚体位移和三个常量应变。反映了三个刚体位移和三个常量应变。(2)(2)必须保证相邻单元在公共边界处的位移连续性。必须保证相邻单元在公共边界处的位移连续性。(线性函数的特性线性函数的特性)现在学习的是第19页，共72页 形态函数Ni具有以下性质：1)在单元结点上形态函数的值为1或为0。2）在单元中的任意一点上，三个形态函数之和等于1。用 来计算三角形面积时，要注意单元结点的排列顺序，当三个结点i，j，m取逆时针顺序时，当三个结点i，j，m取顺时针顺
13、序时，现在学习的是第20页，共72页图 4-3现在学习的是第21页，共72页2-3 单元位移函数单元位移函数 作业：图示等腰三角形单元，求其形态矩阵作业：图示等腰三角形单元，求其形态矩阵N,N,应变矩阵应变矩阵,应力矩应力矩阵阵,单元刚度矩阵。单元刚度矩阵。现在学习的是第22页，共72页2-3 单元位移函数单元位移函数 由三角形的面积由三角形的面积现在学习的是第23页，共72页2-4 单元载荷移置单元载荷移置 连续弹性体离散为单元组合体时，为简化受力情况，需连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向结点把弹性体承受的任意分布的载荷都向结点移置移置(分解分解
14、)，而成为，而成为结点载荷。如果弹性体受承受的载荷全都是集中力，则将所有集中结点载荷。如果弹性体受承受的载荷全都是集中力，则将所有集中力的作用点取为结点，就不存在移置的问题，集中力就是结点载荷。力的作用点取为结点，就不存在移置的问题，集中力就是结点载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在结点上。但实际问题往往受有分布的面力和体力，都不可能只作用在结点上。因此，必须进行载荷移置。如果集中力的作用点未被取为结点，该因此，必须进行载荷移置。如果集中力的作用点未被取为结点，该集中力也要向结点移置。集中力也要向结点移置。将载荷移置到结点上，必须遵循将载荷移置到结点上，必须遵循静力等效的
15、原则静力等效的原则。静力等效。静力等效是是指原载荷与结点载荷在任意虚位移上做的虚功相等指原载荷与结点载荷在任意虚位移上做的虚功相等。在一定的位。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。移模式下，移置结果是唯一的，且总能符合静力等效原则。现在学习的是第24页，共72页单元的虚位移可以用结点的虚位移 表示为 表示为，（3-15）令结点载荷为 令结点载荷为现在学习的是第25页，共72页集中力的移置如图所示，在单元内任意一点作用集中力由虚功相等可得，由于虚位移是任意的，则现在学习的是第26页，共72页例题1：在均质、等厚的三角形单元ijm的任意一点p（xp，yp）上作用有集中载荷
16、。现在学习的是第27页，共72页体力的移置令单元所受的均匀分布体力为 由虚功相等可得，现在学习的是第28页，共72页分布面力的移置设在单元的边上分布有面力，同样可以得到结点载荷，现在学习的是第29页，共72页例题：设有均质、等厚的三角形单元ijm，受到沿y方向的重力载荷qy的作用。求均布体力移置到各结点的载荷。现在学习的是第30页，共72页 同理，现在学习的是第31页，共72页例题：在均质、等厚的三角形单元ijm的ij边上作用有沿x方向按三角形分布的载荷，求移置后的结点载荷。取局部坐标s，在i点s=0，在j点s=L，L为ij边的长度。在ij边上，以局部坐标表示的插值函数为，现在学习的是第32页
17、，共72页载荷为 现在学习的是第33页，共72页 例：例：总载荷的总载荷的2/32/3移置到结点移置到结点i i，1/31/3移置到结点移置到结点j j，与原载荷同向，与原载荷同向现在学习的是第34页，共72页现在学习的是第35页，共72页2-5 单元应力矩阵单元应力矩阵 本节利用几何方程、物理方程，实现用结点位移表示单元的应变本节利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。和单元的应力。用结点位移表示单元的应变的表达式为用结点位移表示单元的应变的表达式为 ，BB矩阵称为几何矩阵。矩阵称为几何矩阵。现在学习的是第36页，共72页2-5 单元应力矩阵单元应力矩阵 由物理方程
18、，可以得到单元的应力表达式由物理方程，可以得到单元的应力表达式 为应力矩阵为应力矩阵现在学习的是第37页，共72页2-6 单元刚度矩阵单元刚度矩阵 讨论单元内部的应力与单元的结点力的关系，导出用结点位移表讨论单元内部的应力与单元的结点力的关系，导出用结点位移表示结点力的表达式。示结点力的表达式。由应力推算结点力，需要利用平衡方程。第一章中已经用虚功方由应力推算结点力，需要利用平衡方程。第一章中已经用虚功方程表示出平衡方程。程表示出平衡方程。现在学习的是第38页，共72页2-6 单元刚度矩阵单元刚度矩阵 考虑上图三角形单元的实际受力，结点力和内部应力为：考虑上图三角形单元的实际受力，结点力和内部
19、应力为：任意虚设位移，结点位移与内部应变为任意虚设位移，结点位移与内部应变为现在学习的是第39页，共72页2-6 单元刚度矩阵单元刚度矩阵 令实际受力状态在虚设位移上作虚功，外力虚功为令实际受力状态在虚设位移上作虚功，外力虚功为现在学习的是第40页，共72页2-6 单元刚度矩阵单元刚度矩阵 计算内力虚功时，从弹性体中截取微小矩形，边长为计算内力虚功时，从弹性体中截取微小矩形，边长为dxdx和和dydy，厚度，厚度为为t t，图示微小矩形的实际应力和虚设变形。，图示微小矩形的实际应力和虚设变形。现在学习的是第41页，共72页2-6 单元刚度矩阵单元刚度矩阵 微小矩形的内力虚功为微小矩形的内力虚功
20、为 整个弹性体的内力虚功为整个弹性体的内力虚功为现在学习的是第42页，共72页2-6 单元刚度矩阵单元刚度矩阵 根据虚功原理，得根据虚功原理，得 这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。虚应变可以由结点虚位移求出：虚应变可以由结点虚位移求出：代入虚功方程代入虚功方程现在学习的是第43页，共72页2-6 单元刚度矩阵单元刚度矩阵 接上式，将应力用结点位移表示出接上式，将应力用结点位移表示出 有有 令令 则则 建立了单元的结点力与结点位移之间的关系，建立了单元的结点力与结点位移之间的关系，称为单称为单元刚度矩阵。它
21、是元刚度矩阵。它是6*66*6矩阵，其元素表示该单元的各结点沿坐标方矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。行移动而改变。现在学习的是第44页，共72页2-6 单元刚度矩阵单元刚度矩阵 由于由于DD中元素是常量，而在线性位移模式下，中元素是常量，而在线性位移模式下，BB中的元素也是中的元素也是常量，且常量，且 因此因此 可以进一步得出平面应力问题和平面应
22、变问题中的单元刚度矩阵。可以进一步得出平面应力问题和平面应变问题中的单元刚度矩阵。现在学习的是第45页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 已经求出了下列关系已经求出了下列关系现在学习的是第46页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 结点力和结点位移的关系：结点力和结点位移的关系：(以简单平面桁架为例以简单平面桁架为例)平面问题中，离散化的单元组合体极为相似，单元组合体在结点平面问题中，离散化的单元组合体极为相似，单元组合体在结点载荷的作用下，结点对单元、单元对结点都有作用力与反作用力存在，载荷的作用下，结点对
23、单元、单元对结点都有作用力与反作用力存在，大小相等方向相反，统称为结点力。大小相等方向相反，统称为结点力。结点力和结点位移的关系前面已经求出：结点力和结点位移的关系前面已经求出：现在学习的是第47页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义：单元刚度矩阵的物理意义：将将 写成分块矩阵写成分块矩阵 写成普通方程写成普通方程 其中其中 表示结点表示结点S(S=i,j,m)S(S=i,j,m)产生单位位移时，在结点产生单位位移时，在结点r(r=i,j,m)r(r=i,j,m)上所需要施加的结点力的大小。上所需要施加的结点力的大小。现在学习的
24、是第48页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义：单元刚度矩阵的物理意义：将结点力列矩阵将结点力列矩阵 与结点位移列矩阵与结点位移列矩阵 均展开成均展开成(6*1)(6*1)阶列阶列矩阵，单元刚度矩阵相应地展开成矩阵，单元刚度矩阵相应地展开成(6*6)(6*6)阶方阵：阶方阵：元素元素K K的脚码，标有的脚码，标有“-”的表示水平方向，没有标的表示水平方向，没有标“-”的表示的表示垂直方向。垂直方向。现在学习的是第49页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义：单元刚度
25、矩阵的物理意义：单元刚度矩阵的每一个元素都有明显的物理意义。单元刚度矩阵的每一个元素都有明显的物理意义。表示结点表示结点S(S=i,j,m)S(S=i,j,m)在水平方向、垂直方向产生单在水平方向、垂直方向产生单位位移时，在结点位位移时，在结点r(r=i,j,m)r(r=i,j,m)上分别所要施加的水平结点力和垂直结上分别所要施加的水平结点力和垂直结点力的大小。例如点力的大小。例如 表示结点表示结点j j在垂直方向产生单位位移时，在在垂直方向产生单位位移时，在结点结点i i所需要施加的水平结点力的大小。所需要施加的水平结点力的大小。现在学习的是第50页，共72页2-7 单元刚度矩阵的物理意义及
26、其性质单元刚度矩阵的物理意义及其性质 单元刚度矩阵的性质：单元刚度矩阵的性质：1)1)对称性：对称性：是对称矩阵是对称矩阵 2)2)奇异性：奇异性：是奇异矩阵，是奇异矩阵，单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。零。由对称性可知，各行元素的总和也为零。现在学习的是第51页，共72页2-7 单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义及其性质 单元刚度矩阵的性质：单元
27、刚度矩阵的性质：例题：求下图所示单元的刚度矩阵，设例题：求下图所示单元的刚度矩阵，设1 1、求、求BB2 2、求、求 DD3 3、求、求 SS4 4、求、求 现在学习的是第52页，共72页2-8 整体分析整体分析 将各单元组合成结构，进行整体分析。将各单元组合成结构，进行整体分析。整体分析分整体分析分4 4个步骤个步骤1 1、建立整体刚度矩阵；、建立整体刚度矩阵；2 2、根据支承条件修改整体、根据支承条件修改整体刚度矩阵；刚度矩阵；3 3、解方程组，求出结点位、解方程组，求出结点位移；移；(消去法与叠加法消去法与叠加法)4 4、根据结点位移求出应力。、根据结点位移求出应力。现在学习的是第53页
28、，共72页2-8 整体分析整体分析 图示结构的网格共有四个图示结构的网格共有四个单元和六个结点。在结点单元和六个结点。在结点1 1、4 4、6 6共有四个支杆支承。结构的共有四个支杆支承。结构的载荷已经转移为结点载荷。载荷已经转移为结点载荷。整体分析的四个步骤：整体分析的四个步骤：1 1、建立整体刚度矩阵；、建立整体刚度矩阵；2 2、根据支承条件修改整体刚度、根据支承条件修改整体刚度矩阵；矩阵；3 3、解方程组，求结点位移；、解方程组，求结点位移；4 4、根据结点位移求出应力。、根据结点位移求出应力。对单元的分析得出单元刚度矩阵，下面，将各单元组合成结构，进行对单元的分析得出单元刚度矩阵，下面
29、，将各单元组合成结构，进行整体分析。整体分析。现在学习的是第54页，共72页2-8 整体分析整体分析 1 1、建立整体刚度矩阵、建立整体刚度矩阵(也叫作结构刚度矩阵也叫作结构刚度矩阵)上图中的结构有六个结点，共有上图中的结构有六个结点，共有1212个结点位移分量和个结点位移分量和1212个结点个结点力分量。由结构的结点位移向量求结构的结点力向量时，转换关系力分量。由结构的结点位移向量求结构的结点力向量时，转换关系为：为：分块形式为：分块形式为：其中子向量其中子向量 和和 都是二阶向量，子矩阵都是二阶向量，子矩阵 是二行是二行二列矩阵。整体刚度矩阵二列矩阵。整体刚度矩阵KK是是12*1212*1
30、2阶矩阵。阶矩阵。现在学习的是第55页，共72页2-8 整体分析整体分析 2 2、根据支承条件修改整体刚度矩阵。、根据支承条件修改整体刚度矩阵。建立整体刚度矩阵时，每个结点的位移当作未知量看待，没有考虑具体的支建立整体刚度矩阵时，每个结点的位移当作未知量看待，没有考虑具体的支承情况，因此进行整体分析时还要针对支承条件加以处理。承情况，因此进行整体分析时还要针对支承条件加以处理。在上图的结构中，支承条件共有四个，即在结点在上图的结构中，支承条件共有四个，即在结点1 1、4 4、6 6的四个支杆处相的四个支杆处相应位移已知为零：应位移已知为零：建立结点平衡方程时，应根据上述边界条件进行处理。建立结
31、点平衡方程时，应根据上述边界条件进行处理。3 3、解方程组，求出结点位移。、解方程组，求出结点位移。通常采用消元法和迭代法两种方法。通常采用消元法和迭代法两种方法。4 4、根据结点位移求出应力。、根据结点位移求出应力。现在学习的是第56页，共72页2-9 整体刚度矩阵的形式整体刚度矩阵的形式 整体刚度矩阵整体刚度矩阵 是单元刚度矩阵是单元刚度矩阵 的集成。的集成。1 1、刚度集成法的物理概念：、刚度集成法的物理概念：刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点力。力。由由2-82-8节的例题可见，与结点节的例题可见，与结点2 2和
32、和3 3相关的单元有单元相关的单元有单元和和，当，当结点结点3 3发生单位位移时，相关单元发生单位位移时，相关单元和和同时在结点同时在结点2 2引起结点力，引起结点力，将相关单元在结点将相关单元在结点2 2的结点力相加，就得出结构在结点的结点力相加，就得出结构在结点2 2的结点力。的结点力。由此看出，结构的刚度系数是相关单元的刚度系数的集成，结构刚度由此看出，结构的刚度系数是相关单元的刚度系数的集成，结构刚度矩阵中的子块是相关单元的对应子块的集成。矩阵中的子块是相关单元的对应子块的集成。现在学习的是第57页，共72页2-9 整体刚度矩阵的形式整体刚度矩阵的形式 2 2、刚度矩阵的集成规则：、刚
33、度矩阵的集成规则：先对每个单元求出单元刚度矩阵先对每个单元求出单元刚度矩阵 ，然后将其中的每个子块，然后将其中的每个子块 送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵KK的子块，从而得出结构刚度矩阵的子块，从而得出结构刚度矩阵KK。关键是如何找出关键是如何找出 中的子块在中的子块在KK中的对应位置。这需要了解单元中的对应位置。这需要了解单元中的结点编码与结构中的结点编码之间的对应关系。中的结点编码与结构中的结点编码之间的对应关系。现在学习的是第58页，共72页2-9 整体刚度矩阵的形式整体刚度矩阵的形式 2 2
34、、刚度矩阵的集成规则：、刚度矩阵的集成规则：结构中的结点编码称为结构中的结点编码称为结点的总码，各个单元的三结点的总码，各个单元的三个结点又按逆时针方向编为个结点又按逆时针方向编为i,j,m,i,j,m,称为结点的局部码。称为结点的局部码。单元刚度矩阵中的子块是单元刚度矩阵中的子块是按结点的局部码排列的，而结按结点的局部码排列的，而结构刚度矩阵中的子块是按结点构刚度矩阵中的子块是按结点的总码排列的。因此，在单元的总码排列的。因此，在单元刚度矩阵中，把结点的局部码刚度矩阵中，把结点的局部码换成总码，并把其中的子块按换成总码，并把其中的子块按照总码次序重新排列。照总码次序重新排列。现在学习的是第5
35、9页，共72页2-9 整体刚度矩阵的形式整体刚度矩阵的形式 以单元以单元为例，局部码为例，局部码i,j,mi,j,m对应于总码对应于总码5,2,45,2,4，因此，因此 中的子块按照总码重新排列后，得出扩大矩阵中的子块按照总码重新排列后，得出扩大矩阵 为：为：现在学习的是第60页，共72页2-9 整体刚度矩阵的形式整体刚度矩阵的形式 用同样的方法可得出其他单元的扩大矩阵用同样的方法可得出其他单元的扩大矩阵 将各单元的扩大矩阵迭加，即得出结构刚度矩阵将各单元的扩大矩阵迭加，即得出结构刚度矩阵KK：集成规则包含搬家和迭加两个环节：集成规则包含搬家和迭加两个环节：1 1、将单元刚度矩阵、将单元刚度矩
36、阵 中的子块搬家，得出单元的扩大刚度矩阵中的子块搬家，得出单元的扩大刚度矩阵 。2 2、将各单元的扩大刚度矩阵、将各单元的扩大刚度矩阵 迭加，得出结构刚度矩阵迭加，得出结构刚度矩阵KK。(例题略例题略)现在学习的是第61页，共72页现在学习的是第62页，共72页2-10 支承条件的处理支承条件的处理 整体刚度矩阵整体刚度矩阵KK求出后，结构的结点力求出后，结构的结点力FF可表示为可表示为 在无支杆的结点处，结点力就等于已知的结点载荷。在有支杆的结点在无支杆的结点处，结点力就等于已知的结点载荷。在有支杆的结点处，则求结点力时，还应把未知的支杆反力考虑在内。如果用处，则求结点力时，还应把未知的支杆
37、反力考虑在内。如果用PP表示结表示结点载荷和支杆反力组成的向量，则结点的平衡方程为点载荷和支杆反力组成的向量，则结点的平衡方程为 根据支承条件对平衡方程加以处理。先考虑结点根据支承条件对平衡方程加以处理。先考虑结点n n有水平支杆的情有水平支杆的情况。与结点况。与结点n n水平方向对应的平衡方程是第水平方向对应的平衡方程是第2n-12n-1个方程，个方程，根据支承情况，上式应换成根据支承情况，上式应换成 ，即在，即在KK中，第中，第2n-2n-1 1行的对角线元素行的对角线元素 应改为应改为1 1，该行全部非对角线元素应改为，该行全部非对角线元素应改为0 0。在在PP中，第中，第2n-12n-
38、1个元素个元素 应改为应改为0 0。此外，为了保持矩阵此外，为了保持矩阵KK的对称性，则第的对称性，则第2n-12n-1列全部非对角线元列全部非对角线元素也改为素也改为0 0。现在学习的是第63页，共72页2-10 支承条件的处理支承条件的处理 同理，如果结点同理，如果结点n n有竖向支杆，则平衡方程的第有竖向支杆，则平衡方程的第2n2n个方程应改为个方程应改为 ，为此，在矩阵，为此，在矩阵KK中，第中，第2n2n行的对角线元素改为行的对角线元素改为1 1，该行全部非对角，该行全部非对角线元素改为线元素改为0 0，同时，第，同时，第2n2n列全部非对角线元素也改为列全部非对角线元素也改为0 0
39、。在。在PP中，中，第第2n2n个元素改为个元素改为0 0。现在学习的是第64页，共72页2-10 支承条件的处理支承条件的处理 2-8 2-8节中的结构，结点节中的结构，结点1 1有水平支杆，结点有水平支杆，结点2 2有两个支杆，结点有两个支杆，结点3 3有有竖向支杆。对支承条件处理后，矩阵修改为：竖向支杆。对支承条件处理后，矩阵修改为：现在学习的是第65页，共72页2-10 支承条件的处理支承条件的处理 最后考虑支点最后考虑支点n n的水平位移的水平位移 为已知非零值为已知非零值 的情况，的情况，这时的支承条件为这时的支承条件为 对平衡方程的第对平衡方程的第2n-12n-1个方程作如下修改
40、：个方程作如下修改：对角线系数对角线系数 乘以一个大数乘以一个大数A(A(例如例如A=A=101010 10)；右边自由项换成右边自由项换成 ；其余各项保持不变，即有其余各项保持不变，即有 实际上，上式中除包含大数实际上，上式中除包含大数A A的两项外，其他各项相对地都比的两项外，其他各项相对地都比较小，可以忽略不计，因此与支承条件等价。较小，可以忽略不计，因此与支承条件等价。对于支点的竖向位移对于支点的竖向位移 的支承条件，也可用同样的方的支承条件，也可用同样的方法进行修改。法进行修改。现在学习的是第66页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 在有限元法中，整体刚度矩阵的
41、阶数通常是很高的，在解算在有限元法中，整体刚度矩阵的阶数通常是很高的，在解算时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体刚度矩时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体刚度矩阵的特性达到节省存贮容量的途径。阵的特性达到节省存贮容量的途径。1 1、对称性。、对称性。只存贮矩阵的上三角部分，节省近一半的存贮容量。只存贮矩阵的上三角部分，节省近一半的存贮容量。2 2、稀疏性。、稀疏性。矩阵的绝大多数元素都是零，非零元素只占一小部分。矩阵的绝大多数元素都是零，非零元素只占一小部分。现在学习的是第67页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 2 2、稀疏性。、稀疏性。矩阵的绝
42、大多数元素都是零，非零元素只占一小部分。矩阵的绝大多数元素都是零，非零元素只占一小部分。结点结点5 5只与周围的六个结点只与周围的六个结点(2(2、3 3、4 4、6 6、8 8、9)9)用三角形单元相连，用三角形单元相连，它们是它们是5 5的相关结点。只有当这七的相关结点。只有当这七个相关结点产生位移时，才使该结个相关结点产生位移时，才使该结点产生结点力，其余结点发生位移点产生结点力，其余结点发生位移时并不在该结点处引起结点力。因时并不在该结点处引起结点力。因此，在矩阵此，在矩阵KK中，第中，第5 5行的非零子行的非零子块只有七个块只有七个(即与相关结点对应即与相关结点对应的七个子块的七个子
43、块)。现在学习的是第68页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 2 2、稀疏性。、稀疏性。一般，一个结点的相关结点一般，一个结点的相关结点不会超过九个，如果网格中有不会超过九个，如果网格中有200200个结点，则一行中非零子块的个结点，则一行中非零子块的个数与该行的子块总数相比不大个数与该行的子块总数相比不大于于9/2009/200，即在，即在5%5%以下，如果网以下，如果网格的结点个数越多，则刚度矩格的结点个数越多，则刚度矩阵的稀疏性就越突出。阵的稀疏性就越突出。利用矩阵利用矩阵KK的稀疏性，可设法的稀疏性，可设法只存贮非零元素，从而可大量地节只存贮非零元素，从而可大量地
44、节省存贮容量。省存贮容量。现在学习的是第69页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 3 3、带形分布规律。、带形分布规律。上图中，矩阵上图中，矩阵KK的非零元素分布在以对角线为中心的带的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中形区域内，称为带形矩阵。在半个带形区域中(包括对角线元素在包括对角线元素在内内)，每行具有的元素个数叫做半带宽，用，每行具有的元素个数叫做半带宽，用d d表示。半带宽的一般计表示。半带宽的一般计算公式是：算公式是：半带宽半带宽 d=(d=(相邻结点码的最大差值相邻结点码的最大差值 +1)*2+1)*2 上图中相邻结点码的
45、最大差值为上图中相邻结点码的最大差值为4 4，故，故d=(4+1)*2=10d=(4+1)*2=10 利用带形矩阵的特点并利用对称性，可只存贮上半带的元利用带形矩阵的特点并利用对称性，可只存贮上半带的元素，叫半带存贮。素，叫半带存贮。现在学习的是第70页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 图图(a)(a)中的矩阵中的矩阵KK为为n n行行n n列矩阵，半带宽为列矩阵，半带宽为d d。半带存贮时从。半带存贮时从KK中取出上半带元素，按图中取出上半带元素，按图(b)(b)中的矩阵中的矩阵 的排列方式进行存贮，即将的排列方式进行存贮，即将上半部斜带换成竖带。存贮量上半部斜带换成
46、竖带。存贮量n*dn*d，存贮量与，存贮量与KK中元素总数之比为中元素总数之比为d/nd/n，d d值越小，则存贮量约省。值越小，则存贮量约省。矩阵矩阵K K 矩阵矩阵 对角线对角线 第第1 1列列 r r行行 r r行行 r r列列 4545度斜线度斜线r r行行s s列列 r r行行s-r+1s-r+1列元素列元素 元素元素现在学习的是第71页，共72页2-11 整体刚度矩阵的特点整体刚度矩阵的特点 同一网格中，如果采用不同的结点编码，则相应的半带宽同一网格中，如果采用不同的结点编码，则相应的半带宽d d也可能也可能不同。如图，是同一网格的三种结点编码，相邻结点码的最大差值分别为不同。如图，是同一网格的三种结点编码，相邻结点码的最大差值分别为4 4、6 6、8 8，半带宽分别为，半带宽分别为1010、1414、1818。因此，应当采用合理的结点编。因此，应当采用合理的结点编码方式，以便得到最小的半带宽，从而节省存贮容量。码方式，以便得到最小的半带宽，从而节省存贮容量。现在学习的是第72页，共72页