12参双参数直角梯形板元的构造与收敛性分析.pdfVIP

第 37卷第 4期 郑 州 大 学 学 报 (理 学 版) Vo1．37NO．4 2005年 12月 J．ofZhengzhouUniv．(Nat．Sci．Ed．) Dec．2OO5 12参双参数直角梯形板元的构造及收敛性分析 孙会霞 ， 尹 丽 (1．河南工业大学理学院 郑州 450052；2．郑州轻工业学院数理系 郑州 450052) 摘要：考虑特殊的四边形单元 一 直角梯形单元，构造出一类 12自由度直角梯形板元；为了降低 自由度，节省计算 量，利用双参数法将 自由度进行离散 ，得到一类双参数 12参直角梯形板元．证明了该单元对薄板弯曲问题 的收敛性． 关键词 ：双参数法 ；直角梯形 ；收敛性 中图分类号：0242．21 文章编号 ：1671—6841(2OO5)04—0022—06 0 引言 1991年，陈绍春教授和石钟慈院士提出了构造非协调元的双参数法口]，利用这种方法 已经构造 出不少 有效的双参数矩形板元和三角形板元l2]，但对于任意四边形板元 尚不多见．这是因为 ，一方面要在参考元 R上构造形函数 ，而不可能在任意四边形上构造形 函数；另一方面 ，为了保证收敛性 ，单元 K上 的形函数 (而 非R上的)跨越 K 的边界应具有某种连续性．但是，由于从参考元 R 到一般单元K 上的变换是非线性 的， 形函数在参考元 R上的连续性在一般单元K 上得不到保证．因此 ，要构造参考元 R上适定的形函数并且使 其在K上具有保证收敛所需要的连续性是很困难的．本文从特殊的四边形单元——直角梯形单元开始，构 造出一类 12自由度双参数直角梯形板元．为了降低 自由度 ，节省计算量，利用双参数法将 自由度进行离散； 并以薄板弯曲问题为例 ，讨论 了它 的收敛性． 1 12参直角梯形板元的构造 设 K是xoy平面的一般直角梯形单元，其顶点为 a (0，O)，a (n，O)，a。(n—e，6)，a (O，6)，四边 中点为 1(n +n)，n。一 1(n +n。)，n一 1(n。+n )，n一 1( n 一 n+n)，z一 是梯形斜边的长度．设 R一 [O，1]×Eo，1]是参考单元，其顶点为 丘 (0， O)，a2(1，0)，3(1，1)， (0，1)，四边 中点为 a5(1 o)，6(1，1) ， ， 丘(告，1)，丘8(o，告)，如 图 1所示．不失一般性 ，总假设原点 a (0，O) 是直角梯形的直角顶点 ，这是因为，对任意的 直角梯形，总可以通过平移或旋转使其直角 顶点位于原点位置．设从 R 到K 的双线性 变换 F 为 z—n ——e刀，y=br] (1) 图l 直角梯形单元和参考元示意图 收稿 日期 ：2005—03—12 基金项目：国家 自然科学基金资助项 目，编号河南省 自然科学基金资助项 目，编号 411011300 作者简介 孙会霞 (1963一)，女，副教授，博士 ，主要从事有限元研究；E-mail：shx@haut．edu．cn． 从 而 第4期 孙会霞等 ：12参双参数直角梯形板元的构造及收敛性分析 23 定义 上的有限元空间{ ，户，}如下 ：取 户()一span{P ，P ，…，P }，其中， Pl一 (t— )(1一 r／)(1— 25)(1— 27)， P2一一 (1— 7／)(1— 27／)(1— 25) ， Ps一 (1— 25)(1— 27) ， P 一一 (1一 )(1— 27)(1—