单元应变矩阵应力矩阵.ppt

参照式（7-2）和（7-3）P6，此基本单元位移函数可写为： （7-40） ? ? 1 2 3 4 ?= -1 ?= 1 ?= -1 ?= 1 图7-12 基本单元 为了得出实际单元的位移模式和局部坐标与整体坐标之间的变换关系，引入一个四节点的正方形单元，称基本单元（图7-12） 。 其中，形函数应为： 引入新变量 ??i、??i （i=1,2,3,4) 基本单元的形函数被写成： （7-41） 现在，把基本单元的位移模式（7-40）和形函数式（7-41）移用于图（7-11）所示的实际单元，则实际单元的位移模式取为： （7-40） 在结点处 ： 在其它结点处： 且，式中的形函数Ni仍由式（7-41）确定。而把式（7-41）中的?、?理解为图7-11所示实际单元的局部坐标，?i、?i便是实际单元中节点i的局部坐标。 （7-42） 利用形函数的上述性质，可以将任意四边形的整体坐标写成： 任意四边形单元中结点的整体坐标， 如果它已知，那么(7-42)表示了局部坐标与整体坐标的变换 另一方面，式（7-42）表明了实际单元中局部坐标（?、?）与整体坐标（x、y）的一一对应关系，是一个坐标变换式。 实际单元是任意四边形四节点单元，基本单元是正方形单元，可以认为：实际单元是对基本单元通过变换得来的。由于实际单元的位移模式中采用了基本单元等同的形函数，这个实际单元就称为等参数单元。 类似于本章3.2节进行的四结点矩形单元的特性分析，可以建立等参单元的应变矩阵、应力矩阵、刚度矩阵、节点力向量等的计算公式。与前面不同之处在于，在等参数单元法中，要将对整体坐标x、y的导数计算和积分计算转换为对局部坐标?、?的导数计算和积分计算。 例：实际单元的结点整体坐标如图(a)中括号内数字所示，基本单元的结点局部坐标如图(b)中括号内数字所示。 图(a) 实际单元 1(0,0) (1) 试验证基本单元上的结点局部坐标与实际单元上对应点的整体坐标的对应关系。 (2) 求基本单元的局部坐标原点( ),在实际单元上的整体坐标(x,y)是多少？ 4(0,1) 3(1,2) 2(2,0) x y 0(3/4,3/4) 0 1(－1,－1) 2(1,－1) 3(1,1) 4(－1,1) 图(b) 基本单元 解 （1）以3结点为例，根据形函数的性质： 在3结点处 ： 应用式（7-42） 说明了由基本单元上结点的局部坐标可映射出实际单元上对应的结点整体坐标 （2）由于 ，由（7-40）P33 得： 代入（7-42） 由上例可知：利用（7-42）在基本单元上任意一点 ，都可以在实际单元上找到一个对应点的坐标(x,y), 这样就把实际单元与基本单元紧密地联系起来。反之，则比较困难，这是因为形函数是一个二次函数。为了避开这个困难，一般都假定基本单元上已知点去求实际单元上的对应点。 * 第7章 平面问题高阶单元 7.1 位移模式阶次的选择 在前面两章中讨论了平面问题三结点三角形单元，其位移模式的最高阶是坐标x、y的一次项。这种位移模式导致单元常应变、常应力特性，单元应变矩阵、应力矩阵、刚度矩阵均为常数矩阵，因此计算非常简单。但这种单元难以反映应力梯度的迅速变化。要想提高计算精度，必须细分网格，增加单元数和点数，因而加大输入数据的工作量。 提高计算精度的另一条有效途径是采用高阶单元。由于高阶单元的应变、应力不再是常数，因此采用少量单元就 可能达到较高的精度。图7-1悬臂梁分别采用高、低阶单元计算就是一个典型的例子。 h 4h P A B 悬臂深梁 解析解：?A=1.0 ?B=1.0 常应变单元：?A=0.866 ?B=0.619 高阶单元：?A=0.99 ?B=0.99 ? ? ? ? ? ? 图7-1 选择位移模式时，第2章提到要考虑解的收敛性，即要考虑到位移模式的完备性和协调性。实际操作中，一般应考虑位移模式的对称性。这是因为，有限元位移模式的选择实际是以帕斯卡（Pascal）三解形基础上的（如图7-2所示），由低价至高阶，顺序选取，组成多项式。多项式中的项数等于单元节点自由度数。如三节点三角形单元，位移 模式取完全一次式，共3项。 六节点三角形单元，位移模 式取完全二次式共6项。如 果某一阶次不能全取，则 应按对称性原则适当选取。