第2章 单元分析.pptVIP

有限元思路框图;思考题;（1）剖分结构时应对单元、节点分别用连续正整数编号。;（2）从结构中取出单元，进行单元分析;第二章 单元分析 ——平面问题常应变单元;单 元 分 析 的 内 容;（2-1）;图2-1;5、单元内任意点的应力列阵???;7、物理方程矩阵式;（2-9）;单 元 分 析 的 内 容;2.2 位移函数和形函数; 不同类型结构会有不同的位移函数。这里，仍以平面问题三角形单元（图2-2）为例，说明设定位移函数的有关问题。; 一个三角形单元有3个节点（以 i、j、m为 序），共有6个节点位移分量。其单元位移或单元节点位移列阵为：;本问题选位移函数为：;3、选取位移函数应考虑的问题 ;（4）位移函数中必须包含单元的刚体位移。 ;例：平面应力矩形板被划分为若干三角形单元。;①、②、③、④单元的位移函数都是;4、形函数;（2-13） ;（2-14） ;式中，A为三角形单元的面积，有:;将式（2-14）代入式（2-12）的第一式，整理后得;（2-12）;（2-16）;形 函 数;（2-16）;[N]为形函数矩阵，进一步写成分块形式：;x;性质1 形函数Ni在节点i上的值等于1，在其它节点上的值等于0。对于本单元，有: ;性质2 在单元中任一点，所有形函数之和等于1。对于本单元，有;x;性质3 在三角形单元的边界ij上任一点（x，y），有: ;性质4 形函数在单元上的面积分和边界上的线积分公式为 ;2.3 单元应变矩阵和应力矩阵;对位移函数（式（2-16））;（2-25）;2、单元应力矩阵;（2-30）;（2-32）;几何关系位移函数;2.4 单元应变能和外力势能的矩阵表达;由于???和???T是常量，提到积分号外，上式可写成 ;2.6节中将明确[k]的力学意义是单元刚度矩阵。式（2-33）便是计算单元刚度矩阵的基本矩阵式。它适合于各种类型的单元。;（1） 体积力势能;注意到式（2-20）; 现在，只考虑弹性体边界上的表面力，它只在部分单元上形成表面力（右下图）。设边界单位长度上受到的表面力如式（2-1）。;（3） 集中力势能; 由单元的应变能U（2-34）和外力势能V（2-36）, 可得单元的总势能?;式（2-38）是从能量原理导出的单元平衡方程。这个方程表达了单元力与单元位移之间的关系。其中，?Fd?和单元节点力?F?具有相同的意义。 ;2.6 单元刚度矩阵; 对于平面应力三角形单元，应变矩阵[B]是常数矩阵，同时弹性矩阵[D]也是常数矩阵，于是式（2-33）可以化简为 ; 将式（2-9）和（2-26）代入上式，;式中子矩阵为2×2矩阵，有 ;（2-42）;图2-7;当坐标系旋转任意角度时，由于公式（2-40）、（2-41）（2-42）中已经包含坐标系的影响，它们仍然照用。;第一步：计算bi、ci和单元面积A ;第二步：求子矩阵 (由式（2-41），算得 );第三步：形成[k] 将[kii]等按式（2-40）组集成[k] 。;（2-43a） ;当a=b时，即等腰直角三角形单元，有 ;② 子程序框图;③ 子程序;3、单元刚度矩阵性质;（2-38）;（2）[k]的每一行或每一列元素之和为零;（3）[k]是对称矩阵 由k单元的表达式，可见，由此可知[k]具有对称性。 ;证明：;（4）单元刚度矩阵是奇异矩阵（ 即[k]的行列式为零） 单元刚度矩阵是在单元处于平衡状态的前提下得出的。单元作为分离体看待，作用在它上面的外力（单元力）必定是平衡力系。然而，研究单元平衡时没有引入约束。承受平衡力系作用的无约束单元，其变形是确定的，但位移不是确定的。所以出现性质（3）中的“平动问题”，即单元可以发生任意的刚体运动。从数学上讲，方程（2-38）的解不是唯一的或不能确定的。由此，单元刚度矩阵一定是奇异的。;2.7 等价节点力;（2-44）;3、等价节点力计算举例;注意到形函数的性质4：;（2-47）;（2）均布面力;注意到形函数性质4 ：; 式（2-48）表明：在ij边上受均布面力的平面问题三角形单元，其等价节点力等于将均布面力合力之半简单地简化到i、j节点上，方向与分布力方向相同。m节点上为零。;（3）线性分布面力;i;（2-49） ; 式（2-49）表明：ij边受线性分布面力： i点为[qsx, qsy]T，j点为0时，其等价节点力可将总载荷的2/3分配给i点，1/3分配给j点，m点为零得出。 ; 单元上的体积力和表面力向结点的移置都是符合直观