对一维有限元问题.docxVIP

4. 在总体坐标系中有一个任意三角形，设三个节点的坐标分别为： 节点1：（x1, y1） 节点2：（x2, y2） 节点3：（x3, y3） （1）试推导其总体坐标系中的形函数。 （2）证明（1）给出的形函数与用面积坐标表达的形函数等价。 解答： （1）对于任意三角形单元，如右图所示，节 点为1、2、3.节点位移为（u1，v1）、（u2，v2）、 （u3，v3），单元节点位移为： u 设单元位移场函数为： u= v= 其中： ~ 为待定常数，在节点处可以解出： 其中，=det DA= + + = 这里其实可以通过计算是三角形面积的两倍。 ， ， ， ， ， ， 把这些结果带入单元内假定的位移场中，得到 类似可以得到 形函数即为： （2）用面积坐标来表达形函数： ox o x y 图７－9 １ ２ ３ M(x,y) 　　　　　 三角形１－２－３的面积记作　　, 令　 则Ｍ的位置确定后，L1、L2、L3的值即完全确定。反过来，L1、L2、L3给定后，Ｍ点在单元内的位置也完全确定。（L1、L2、L3）称为Ｍ点的面积坐标。 显然，节点１的面积坐标为　　（１，０，０） 节点２的面积坐标为　　（０，１，０） 节点３的面积坐标为　　（０，０，１） 点Ｍ的面积坐标（L1、L2、L3）和总体坐标（x、y）之间存在着线性关系 其中 用面积坐标作为三角形单元的自然坐标时，只要利用单元形状函数的性质就可以直接写出相应的插值函数，即针对某一节点的形状函数，有在该节点的函数值为1，而在其他的节点处为0这一性质，用数学表达式表示为： 对图示的三节点单元，形状插值函数由一个线性函数构成。由上式，对于每个角点，可用通过其他两个角点的直线方程的线性函数来构成。例如对于节点1，可用2-3边的边方程来构成它的插值函数，即，其他两个角点也类似，该单元的形状函数为 即线性函数单元的三个形状插值函数就是三角形单元的三个面积坐标。