平面问题高阶单元.pptVIP

现在，我们从任意四边形单元着手，介绍等参数单元的概念。1234xy图7-10任意四边形单元前面讲到的四节点矩形单元虽然比较简单，但难以应用于斜线边界。图7-10所示四节点任意四边形单元容易适应这种边界，但要在整体坐标系内，写出它的统一的形函数又是相当复杂和困难的。第30页，共54页，星期日，2025年，2月5日但是若能找到它与一个规则正方形的关系，就能写出它的统一的位移模式，这可以通过坐标变换来解决。在图7-10所示四边形单元上，用等分四边的两族直线分割该四边形，以两族曲线的中心（?=0、?=0）为原点，沿?、?增大的方向作?轴和?轴，并令四边的?=±1、?=±1，就得出一组新坐标系??（图7-11）。1234xy图11实际单元?=-1?=1?=1?=-1??这里，?、?是一种局部（单元）坐标，它只应用于单元范围内。而x，y是整体（结构）坐标，它适用于所有的单元。图中的任意四边形单元是研究对象，称为实际单元。第31页，共54页，星期日，2025年，2月5日参照式（7-2）和（7-3）P6，此基本单元位移函数可写为：（7-40）??1234?=-1?=1?=-1?=1图7-12基本单元为了得出实际单元的位移模式和局部坐标与整体坐标之间的变换关系，引入一个四节点的正方形单元，称基本单元（图7-12）。第32页，共54页，星期日，2025年，2月5日其中，形函数应为：引入新变量??i、??i（i=1,2,3,4)基本单元的形函数被写成：第33页，共54页，星期日，2025年，2月5日（7-41）现在，把基本单元的位移模式（7-40）和形函数式（7-41）移用于图（7-11）所示的实际单元，则实际单元的位移模式取为：（7-40）在结点处：在其它结点处：第34页，共54页，星期日，2025年，2月5日且，式中的形函数Ni仍由式（7-41）确定。而把式（7-41）中的?、?理解为图7-11所示实际单元的局部坐标，?i、?i便是实际单元中节点i的局部坐标。（7-42）利用形函数的上述性质，可以将任意四边形的整体坐标写成：任意四边形单元中结点的整体坐标，如果它已知，那么(7-42)表示了局部坐标与整体坐标的变换第35页，共54页，星期日，2025年，2月5日另一方面，式（7-42）表明了实际单元中局部坐标（?、?）与整体坐标（x、y）的一一对应关系，是一个坐标变换式。实际单元是任意四边形四节点单元，基本单元是正方形单元，可以认为：实际单元是对基本单元通过变换得来的。由于实际单元的位移模式中采用了基本单元等同的形函数，这个实际单元就称为等参数单元。类似于本章3.2节进行的四结点矩形单元的特性分析，可以建立等参单元的应变矩阵、应力矩阵、刚度矩阵、节点力向量等的计算公式。与前面不同之处在于，在等参数单元法中，要将对整体坐标x、y的导数计算和积分计算转换为对局部坐标?、?的导数计算和积分计算。第36页，共54页，星期日，2025年，2月5日例：实际单元的结点整体坐标如图(a)中括号内数字所示，基本单元的结点局部坐标如图(b)中括号内数字所示。图(a)实际单元1(0,0)(1)试验证基本单元上的结点局部坐标与实际单元上对应点的整体坐标的对应关系。(2)求基本单元的局部坐标原点(),在实际单元上的整体坐标(x,y)是多少？4(0,1)3(1,2)2(2,0)xy0(3/4,3/4)01(－1,－1)2(1,－1)3(1,1)4(－1,1)图(b)基本单元第37页，共54页，星期日，2025年，2月5日解（1）以3结点为例，根据形函数的性质：在3结点处：应用式（7-42）说明了由基本单元上结点的局部坐标可映射出实际单元上对应的结点整体坐标（2）由于，由（7-40）P33得：第38页，共54页，星期日，2025年，2月5日平面问题高阶单元第1页，共54页，星期日，2025年，2月5日可能达到较高的精度。图7-1悬臂梁分别采用高、低阶单元计算就是一个典型的例子。h4hPAB悬臂深梁解析解：?A=1.0?B=1.0常应变单元：?A=0.866