有限元第三章 单元类型及单元刚度矩阵.ppt

三、二维单元及其单元刚度阵 ●形状函数 利用形状函数的性质，可得 验证 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 ●位移函数 其中 为几何矩阵 ●单元应变 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 其中 ●几何矩阵 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 ●几何矩阵 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 平面应力问题的应力为 单元刚度矩阵 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 其中 对于平面应变问题 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 2.矩形单元 ●一次矩形单元 三、二维单元及其单元刚度阵 2.矩形单元 ●二次矩形单元 x y 1 2 3 o 4 8 7 6 5 令 八个节点新坐标 三、二维单元及其单元刚度阵 2.矩形单元 ●二次矩形单元 x y 1 2 3 o 4 8 7 6 5 ●形状函数 利用形状函数的性质，可得 三、二维单元及其单元刚度阵 2.矩形单元 ●二次矩形单元 1.三角形单元 ●一次三角形单元 三、二维单元及其单元刚度阵 ●面积坐标 如图所示，在三角形单元A1A2A3中，有任意一点P(x,y)连接PA1、PA2、PA3,得到三个小三角形：ΔPA2A3、ΔPA3A1、ΔPA1A2,记面积比为： 称λ1、λ2、λ3为P点的面积坐标 x y A1 A2 A3 P(x,y) o 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标 由于λ1+λ2+λ3=1，因此该三个坐标不独立。其负号的规定为：分子分母对应的三角形顶点编号转向相同时为正，反之为负，由于三角形A1A2A3的顶点编号一般规定为逆时针，因此，子三角形顶点编号为逆时针时面积坐标为正，反之为负。 x y A1 A2 A3 P(x,y) o 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标 三角形中一些特殊点的面积坐标 顶点 边中点 形心 x y A1 A2 A3 G A6 A5 A4 o 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 三角形三个顶点在直角坐标系中的坐标为（xi,yi）,则ΔA1A2A3的面积为 类似地三个小三角形的面积依次为 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 其中 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 同理 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 其中 令 i=1,2,3;i、j、k按1,2,3轮转 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 则 于是有 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 对于单元的三个角点，应有 A1点 A2点 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 A3点 重写上述三式 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●面积坐标与直角坐标的关系 所以有 于是有 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●形状函数 根据形状函数的定义 i=1时，对于A1、A2、A3点 i=1,2,3 所以 三、二维单元及其单元刚度阵 1.三角形单元 ●一次三角形单元 ●形状函数 同理，i=2、3时，对于A1、A2、A3点 所