结构力学：有限单元法-5.pptVIP

四.采用面积坐标时的单元分析 j i k y x P 1 .面积坐标 三角形单元中任一点P 可用直角坐标 (x , y)表示。 连P i、 P j、 P k，则可得三个小三角形。它们和 大三角形?123的面积比，记作 ?Pij ?Pik ?Pjk L i= ?P jk/ ?ijk L j= ?P ik/ ?ijk L k= ?P ij/ ?ijk 面积坐标 由于 Li+ Lj + Lk = 1，只有两个是独立的。 三角形中任一点P 的位置可用面积坐标Li、 Lj 确定。 当P 点在i结点时Lj = Lk= 0， Li= 1。余类推。 可见面积坐标具有“形函数”的性质。 Li+ Lj + Lk=1 四.采用面积坐标时的单元分析 j i k y x P 2 .位移模式 ?Pij ?Pik ?Pjk 由于面积坐标有形函数性质，因此根据试凑法可得到形函数矩阵。 形函数 Ni=Li 面积坐标 如果结点 i 位移为ui、vi，（i=i，j，k）则单元位移模式（位移场）为 u=? Niui ； v=? Nivi 面积坐标和直角坐标关系: j i k y x P ?Pij ?Pik ?Pjk j i k y x P ?Pij ?Pik ?Pjk Li+ Lj + Lk=1 x=Lixi+ Ljxj + Lkxk y=Liyi+ Ljyj + Lkyk 后面的分析过程与结果与 前面广义坐标法一致. 3.1 常应变三角形单元 3 平面问题的有限元分析 一.离散化 单元结点位移向量 3.2 矩形双线性单元 水坝 1 2 3 4 (x1 , y1) a a b b (x2 , y2) (x3 , y3) (x4 , y4) 单元结点力向量 同理,有 二.单元分析 若用广义坐标法,则与三角形单元类似的可得到 下面用试凑法确定形函数矩阵 1.单元位移 设单元内位移为 1 2 3 4 (x1 , y1) a a b b (x2 , y2) (x3 , y3) (x4 , y4) 4 1 2 3 (-1,-1) 1 1 1 1 (1,-1) (1,1) (-1 ,1) 令 由形函数性质 可设 正则(自然) 坐标系 或 思考题:这种单元是收敛的单元吗?为什么? 2.单元应力与应变 同理,有 4 1 2 3 (-1,-1) 1 1 1 1 (1,-1) (1,1) (-1 ,1) 或 思考题:这种单元是收敛的单元吗?为什么? 2.单元应力与应变 2.单元应力与应变 应变矩阵 应变矩阵 应力矩阵 对于平面应力问题 应力矩阵 对于平面应力问题 从应变矩阵和应力矩阵可见: 单元内应力和应变沿x方向是线性变化的, 沿y向也是线性变化的. 这是由所设位移模式所决定的. 3.单元刚度矩阵和单元等效结点荷载 用虚位移原理或势能原理可推得 对于平面应力问题