11_平面问题有限元分析.pptVIP

11.7.2 整体刚度矩阵 假设整体结构被划分为ne个单元和n个节点，n1个自由度，在整体坐标系下，对于每个单元均有： i j m 图12 x y s l 在ij边上有： 将?qs?和上式代入式(1-46），有 (1-49） 式(1-49）表明：ij边受线性分布面力： i点为[qsx, qsy]T，j点为0 时，其等价节点力可将总载荷的2/3分配给i点，1/3分配给j点，m点为零得出。 x y i j m qsi qs 体积力和表面力向节点的移置符合静力等效原理的前提条件是：线性位移模式。 § 11.7 系统分析 11.7.1 坐标系 研究各离散单元集合成整体结构，集合整体结构的平衡和变形协调，建立整体结构平衡方程。 单元分析时采用的坐标系成为局部坐标或单元坐标（单元刚度矩阵的通用性）。而结构系统分析时，必须在统一的坐标系内进行（各力学量才能叠加），称为“结构坐标”或“整体坐标”，如图13所示。 单元坐标系下，单元位移、单元力、单元刚度矩阵表示为： 整体坐标系下，单元位移、单元力、单元刚度矩阵表示为： X Y X Y P ○ ○ ○ ○ ○ P 图（13） (a) 平面桁架(杆件单元) 悬臂深梁 (平面三角形单元) x y x y x y x y 如何从单元坐标转化为结构坐标将在第13章中讨论。 将上述这些方程集合起来（在整体坐标下叠加），便可得到整个结构的平衡方程。为此，需要将[k]、{δ}、{F}维度膨胀，分别扩大为n1×n1、n1×1和n1×1的矩阵才能相加。膨胀后，原有节点号对应位置的元素不变，而其它元素均为零。 单元应变能写成 (1-34） 2、 单元外力势能 单元受到的外力一般包括体积力、表面力和集中力。自重属于体积力范畴。表面力指作用在单元表面的分布载荷，如风力、压力，以及相邻单元互相作用的内力等。 （1） 体积力势能 单位体积中的体积力如式(1-35）所示。 单元上体积力具有的势能Vv为 (1-35） i j m x y · qVx qVy i j m x y · u v 注意到式(1-20） 有 (1-20） （2） 表面力势能 面积力虽然包括单元之间公共边上互相作用的分布力，但它们属于结构内力，成对出现，集合时互相抵消，在结构整体分析时可以不加考虑，因此单元分析时也就不予考虑。 现在，只考虑弹性体边界上的表面力，它只在部分单元上形成表面力（右下图）。设边界面上单位面积受到的表面力如下式： l—单元边界长度 h—单元厚度 A—表面力作用面积 ① ② ③ ④ ?qs? ?qs?沿厚度均匀分布，则单元表面力的势能Vs为 （3） 集中力势能 当结构受到集中力时，通常在划分单元网格时就把集中力的作用点设置为节点。于是单元集中力?Pc?的势能Vc为 p ① ② ③ ④ （4）总势能 把(1-35）式中原括符内的部分用列阵?Fd?代替， 综合以上诸式，单元外力的总势能V为 (1-35） ?Fd?具有和???相同的行、列数。则 (1-36） 由单元的应变能U(1-34）和外力势能V(1-36），可得单元的总势能? (1-37） 将式(1-37）代入， 根据弹性力学最小势能原理：结构处于稳定平衡的必要和充分条件是总势能有极小值。 3、单元平衡方程 于是有， (1-34） (1-36） 式(1-38）是从能量原理导出的单元平衡方程。这个方程表达了单元力与单元位移之间的关系。其中，?Fd?和单元节点力?F?具有相同的意义。 (1-38） 即得单元平衡方程 § 11.5 单元刚度矩阵 平衡方程(1-38）中的矩阵[k]是单元力和单元位移关系间的系数矩阵，代表了单元的刚度特性，称为单元刚度矩阵。单元刚度矩阵的维数为nj ×nj， nj 是单元位移总数。其一般计算公式为： 1、一般计算公式 它与单元应变矩阵[B]和弹性矩阵[D]有关。 (1-33） 对于平面应力三角形单元，应变矩阵[B]是常数矩阵，同时弹性矩阵[D]也是常数矩阵，于是式(1-33）可以化简为 式中A表示三角形单元的面积。h是单元厚度。 2、平面问题三角形单元刚度矩阵 （1）平面应力三角形单元 (1-39） 将式(1-9）和(1-26）代入上式， 即得平面应力三角形单元刚度矩阵。写成分块形式，有 (1-40） (1-9） (1-26） 式(1-40)中子矩阵[krs]为2×2矩阵，有 (1-41） （2）平面应变