有限单元法的一般原理和表达格式03.ppt

1、绕节点平均法 0 1 2 3 4 A B E C F D 边界线 内节点：表征性较好 边界点：宜用外差值公式 每一单元均是常应力，相邻单元应力阶梯状变化 x y 0 1 2 3 2、两单元平均法 0 0 1 3 2 A B C D E F （1）以相邻单元应力的平均值作为其边中点的应力 （2）用外插值法求边界单元中点的应力 § 3-3 广义坐标有限单元法的一般格式 § 3-4 有限单元法解的性质和收敛性 § 3-5 平面问题高次单元 回顾 三结点三角形单元 位移模式是线性的 单元的应力和应变是常量 计算精度较低 计算精度 有限元分析时，当单元数目确定后，精度主要取决于单元形式。 提高精度也就是减少分析误差，误差的主要来源之一是单元的位移函数与实际位移的差异。 增加单元节点，采用高阶次的位移模式，提高计算精度 1 2 3 1 2 3 4 5 6 1 2 3 4 5 6 7 8 9 10 （a）线性单元 （b）二次单元 （c）三次单元 一、六结点三角形单元 x y 0 m j i 1 2 3 1、位移模式 位移模式的项数与自由度数相当；阶次对称。 (3-37) 常应变三角形单元是在直角坐标下分析的。 若仍在直角坐标系中分析三角形高阶单元，形函数的建立及单元分析得到的公式将比较复杂。 面积坐标 2、面积坐标 面积坐标系就是利用三角形单元面积比的关系，来确定三角形单元中任一点在单元中的位置，称作面积坐标 直角坐标 面积坐标 （1）面积坐标的概念 y x 0 m j p H i h (a) 三角形三结点单元 (3-38) (b) (c) （2）面积坐标和直角坐标的关系 (i)用直角坐标表示面积坐标 （ii）用面积坐标表示直角坐标 (3-39) 3、六结点三角形单元的位移插值函数 式中： (3-41） （3-42） 建立形函数的方法 间接法：如建立常应变三角形单元形函数是采用了广义坐标法。 直接法：利用形函数应满足的条件建立形函数。 边线方程法 （1）求出除该点之外的其余各结点的边线或连线方程； (2) 边线方程连乘得到有待定常数的形函数表达式，并检查是否与位移模式同阶次； (3) 根据形函数在本身结点为1的特性求得选定常数。 如求Ni 边线方程 将结点i的面积坐标（1，0，0）及（Ni)i=1代入上式得 同理可求得 请同学们推导N1,N2,N3 4、六结点三角形单元的刚度矩阵 （1）应变矩阵 （1，2，3） 应变矩阵的元素是面积坐标的一次式，也是直角坐标的一次式。 应变是线性变化的。 2、应力矩阵 (3-26) 例：如图所示为一平面应力状态的直角三角形单元，试求（a）形状函数矩阵Ｎ;（b）应变矩阵B;（c）应力矩阵S（d）单元刚度矩阵ke。 y x j b a m(0,0) i(a,0) （0,b) 解 （a） （b）应变矩阵B （C）应力矩阵S 式中 S=DB 由于 所以 （d）单元刚度矩阵 单元刚度矩阵的性质 1、单元刚度矩阵的物理意义 上标1表示x方向自由度，2表示y方向自由度，后一上标代表单位位移的方向，前一上标代表单位位移引起节点力方向 元素值为单位位移所引起的节点力 J节点产生单位水平位移时在i节点引起的水平节点力分量 J节点产生单位水平位移时在i节点引起的竖直节点力分量 2、单元刚度矩阵为对称矩阵 由功的互等定理中的反力互等可知 3、单元刚度矩阵与单元位置无关（但与方位有关) 单元刚度矩阵与单元位置（刚体平移)无关 4、奇异性 物理意义：在无约束的条件下，单元可以做刚体运动，其位移是不定的。 等效节点荷载 目的：非节点荷载移置为等效节点荷载 移置原则：静力等效 静力等效就是原荷载与等效节点荷载在任何虚位移上所作的虚功相等 y 0 x m i j P 集中力的移置 假想单元发生了虚位移 M点虚位移 单元节点虚位移 原荷载作的虚功 等效节点荷载作的虚功 按照静力等效原则有 (3-27) 面力的移置 体力的移置 单位体积的体力 单位面积的面力 线性位移模式下的载荷移置 （1）y方向的重力W (2)Ij边承受x方向的均布面力q (3)Jm边承受x方向的线性分布 y 0 x i m j l q j l m q i y 0 x s i ① ② ③ y 0 x i ① i ③ 总体刚度方程 方法一：节点平衡法 i i ② 根据外荷载与节点力的平衡条件，得： 节点力是节点位移和刚度系数的乘积 同理可对其它节点进行平衡 总体平衡方程 优点：力学概念清晰，以每一个节点为单位进行平衡 缺点：编写程序比较麻烦 8 7 5 4 1 2 3 5 6 7 8 9 ① ② ③ ④ ⑤ ⑥ ⑦ ⑧ 方法二：直接刚度法 （1）统一使用整体节点编号 总体码 局部码 反映了局部量与整体量的对应关系 （2）依照结构总体的节点自由度数2n