第四章平面问题的有限单元法课件.pptxVIP

第四章 平面问题的有限单元法 1 • 基本条件 • (1)等厚度的薄板； • (2) 体力作用于体内，平行于板的中面，沿板厚不变； • (3) 面力作用于板边，平行于板的中面，沿板厚不变； • (4) 约束作用于板边，平行于板的中面，沿板厚不变。 弹性力学的平面应力问题 2 弹性力学的平面应力问题 由于薄板很薄，应力是连续变化的，又无z向外力，可认为： 由于两板面上无面力和 约束作用： 简化为平面应力问题，仅剩： 其值与z无关 坐标系： 3 基本条件 (1)很长的常截面柱； (2) 体力作用于体内，平行于横截面，沿柱体长度方向不变； (3) 面力作用于柱面，平行于横截面，沿柱体长度方向不变； (4) 约束作用于柱面，平行于横截面，沿柱体长度方向不变。 弹性力学的平面应变问题 4 • 由于截面、外力、约 束沿z 向不变，外力、 约束平行xy面，柱体 非常长：故任何z面 (截面)均为对称面。 弹性力学的平面应变问题 简化为平面应变问题： 其值与z无关 坐标系： 5 以平面应力问题为例介绍 平面问题的有限单元法 6 x方向的位移u， 4 y方向的位移v 8 单元位移分量(4节点) 1 平面应力单元类型简介 3节点三角形单元 i 8节点4边形曲边单元 节点位移分量 每节点2个位移分量(自由度) 平面应力单元类型 3 6 2 8节点单元 7 单元e 5 k 三角形单元 4节点4边形单元 单元e j 单元e 四边形单元 k j l i 7 应力梯度变化比较大的地方，网格应密一些 有应力集中的地方，网格应密一些 单元边界长度不要相差过大 单元各边夹角不要太大 集中载荷处要设置节点 结构不同材料交界面处要设置节点并作为单元边界 结构厚度突变处要设置节点并作为单元边界 分布载荷突变处要设置节点 施加位移约束处要设置节点 注意单元间的连接 平面应力单元网格划分 8 平面应力单元网格划分 病态单元 a-边长差别太大 b-边长差别太大 c-边夹角太大 9 界面 材料B 这样不行 设置节点 设置节点 材料A b c 节点号 x y z 1 0 0 0 2 1 0 0 3 2 0 0 4 2 1 0 5 1 1 0 6 0 1 0 7 0 0.5 0 8 1 0.5 0 9 2 0.5 0 单元 号 节点 i 节点 j 节点k 节点 l 材料 编号 其它 常数 1 1 2 8 7 1 2 7 8 5 6 1 3 5 8 9 4 1 4 3 9 8 8 1 5 8 2 3 3 1 单元节点信息 单元拓扑信息 节点信息 10 以最经典的三角形单元为例 11 单元位移模式概念 单元内任一点的位移要用节点上的位移值近似表达出来，这 就需要假定一个近似函数来表示单元内的位移分布，所选择 的近似函数就称为单元位移函数或单元位移模式。 对于弹性力学平面问题，一般选择多项式 ( polynomial ) 来作为单元内的位移解或插值函数或位移模式。 单元位移函数(位移模式) 多项式的项数越多，结果就越精确，但取多少项由单元形式 2定。 一维单元 二维单元 节点上只有六个位移分量，所以 单元内部位移函数的待定参数不 能超过这个数目。可假设单元内 部位移为x、y的线性函数： 三角单元的位移函数 参数ai 由位移边界条件确定。 13 三角单元的位移函数 节点i 节点j 节点k 于是： 三角单元的位移函数 根据线性代数的知识，可知： 如果令 则： 15 三角单元的位移函数 T*为T的伴随矩阵 其中： 16 代入位移函数公式： 得到： 三角单元的位移函数 把求得的系数 17 Ni,Nj,Nk是坐标的函数，它们反映了单元的位移形态，故 称为三角单元的形态函数 (或形函数) 三角单元的位移函数 表达为矩阵形式： 这里： 18 如图所示三角单元IJK，P为三角单 元内任意一点，其坐标为(x,y) P点在三角单元各角点上产生的形 函数分别是Ni,Nj,Nk 三角单元的位移函数 形函数具有明确的几何意义： 同理： 求P点位移 P点的位移可由节点位移近似表示为 三角单元的位移函数 位移函数运用示例： 已知各节点位移为： 20 三角单元的位移函数 于是： 21 形函数的本质 计算点(x,y)的位移u(x,y)、v(x,y)可用单元内各节点的 位移值ui,vi的加权之和来近似表示，其中，各节点位移 加权系数为关于计算点(x,y)的函数，即为形函数 三角形单元形函数的性质 1、单元节点产生的形函数值为1或0 2、形函数之和等于1 三角单元的位移函数 22 应满足： 单元内位移模式必须是连续的，公共边上位移必 须协调 位移模式必须反映单元的刚体位移 位移模式必须反映单元的常