纯位移平面弹性问题二阶收敛Locking-free有限元.pdfVIP

个反例，说明离散的Korn不等式对GR有限元空间不成立，从而C-R元不能逼近纯应 题提出了一种新的方法，他们用协调线性三角形元来逼近一个位移分量，而用GR元 来逼近另一位移分量。一方面利用非协调单元的自由度使得在保持约束的同时还具 有最优逼近性质；另一方面利用协调元能够去掉单元间非平凡的刚性位移。用这种 思想，明平兵在[644中研究了非协调旋转Q1元。设计一种有限元格式，用双线性协 调元逼近位移的一个分量，用非协调旋转Q1元逼近位移的另外一个分量，从而证明 free的。DongwooSheen等人在【25】中研究了～组C-R形三角形元和类旋转口1四边形元： 过修改变分形式取得了类．Wilson元的最优误差估计；毛士鹏博士在f43】中对非协调 旋转Q1元进行了双参数离散，得到了一个QB 在f411【631中通过构造一个分片驴投影算子的方法，提出了一个八节点和一个十节点的 一阶收敛的，二阶收敛的单元构造都过份复杂，如何构造二阶收敛的、且具有较少自 中都有重要价值。本文通过在不完全三次多项式空间中将算子div邵艮制在B或Q1空间 中的方法，提出了四个二阶收敛的非协调三角形单元和矩形单元，本文分别称之为14． 明它们都有非常好的数值表现。 。各向异性有限元是有限元方法中的一个热点，传统的有限元逼近依赖于正则性条 件【16l和非退化条件[31，即存在常数G满足 hK／pg≤C，VK． 其中hⅣ是单元K的直径，船是K的最大内切球的直径。但是近来的一些研究l瓴柏，鹌，591 表明该条件对某些单元的收敛性并不是必要的，即这些单元可以用在窄形的剖 理。2004年陈绍春教授、石东洋教授在[481中提出了一个更易于判定的各向异性插值定 理，并做出了大量非常有价值的工作[47,49,54，55，硼。 众所周知，wilson元对二阶问题是一个非协调元，它比双线性元有更好的性 4 第1章前言 形剖分收敛。为了将该单元扩展到任意四边形剖分，人们提出了若干改进方法， 出了一个三维类：-Wilson元，并通过给出一个三线性各向异性插值算子，在不需要正则 性条件下，证明了该单元对于扁六面体剖分是收敛的。 本文分为八章，主要工作及章节安排如下： 第1章，介绍平丽弹性问题的有限元方法的发展和现状。 第2章，介绍一些弹性理论，Locking现象，弹性微分方程组的定性理论，其中的 结论在后文中要用到。 第3章，介绍Sobolev空间理论的一些结果以及有限元方法中的一些结论，这些结 论在后文中用到。 细地给出了单元的性质以及能量模和工2一模的误差估计。 敛的。给出了详细构造方法和误差估计。 第6章，给出了三维类一Wilson元的三线性各向异性插值算子和各向异性分析结果。 第7章，提出了一种求解齐次椭圆边值问题方程组的数值算法。 第8章，给出两个数值算例对本文构造的四个单元以及C-R元和w—Q元进行了数值 实验和分析，最后给出了两种Locking现象的数值结果。 5 第2章平面弹性问题概述 平面弹性问题研究弹性体在体力和面力作用下的应变与应力，对于一些特殊的弹 性体，如壳、薄板和细杆等，我们往往把它们简化为平面模型甚至一维的模型，其中 就包括平面弹性问题，这种方法在理论上和实际上都是很有意义的。平面弹性问题模 型不再是单一的方程，而是一个偏微分方程组。本文主要内容是用有限元方法求解纯 位移平面弹性问题，因此首先在本章引入本文所研究的平面弹性问题以及线性协调元 的Locking现象，同时引入本文所采用的一些算子、符号和记号。 符号说明：文中的系数e一般取不同的常数，在不同的地方可能取不同的值，本文 不作区分。 2．1-平面弹性问题 础=盼一㈢ 一。踟1。讹 ．．． 踟1．锄 m删。丽+面， ∞协2一面+万， 掣aa方=Ovt Orl)，一orloth)， 酮=豆1(grad痧+(删弘(拳害哿，争雩》)， 2．1平