从能量改进到精度提高：采用线性位移的杂交宏元方法.pdfVIP

四川大学硕士学位论文 1．引言 早在70年代初，研究者们就致力于设计构造简单、粗网格精度高的有限元方 法来处理固体和结构力学问题。对给定的网格，人们可以构造不同有限元方法， 而这些方法往往具有不同的数值性能。近年来这种所谓高性能有限元的研究再次 人…提出了丰富单元内部非协调位移的非协调方法，其提出的Q6元大大改进了 通过分片检验。 13]导出了著名的假定 应力杂交有限元(P—S)。由于恰当的应力模式选择，该方法具有良好的的数值性能。 与假定应力杂交法最接近的一类方法是Simo和砌蹦【4】提出、后又被其他一些研 究者发展的所谓增强应变方法，其变分基础为胡海昌一鹫津久一郎原理，包含了 位移、应力、增强应变三类变量。2002年周和谢【91对具有高数值性能的丰富应力／ 应变杂交有限元方法进行了统一的理论分析。2001年，基于组合杂交变分原理 调条件导出： (1．1) r哺可以=0，Vf，V曰 t 力模式，得到了一个精度与CH(0．1)相当的高性能4节点杂交应力元ECQ4。 的便利，又能达到比相应位移格式更高的数值精度。 然而，对于三角形网格，通过添加bubble位移的增强应变方法不能改进标准的 四川大学硕士学位论文 三节点线性元即常应变三角形元。2004年周和谢一般性的讨论了对传统位移格式 的改进。按照力学几何观点，利用组合杂交变分原理的一个奇妙表达，阐明了其 增强低阶元的粗网格精度和稳定性的内在机制，并将该机制成功的应用到Allman 三角形元的数值精度改进中。 本文深入探讨了组合杂交方法的零能量误差机制。并将其应用到常应变三角形 元的改进中。我们提出了一类宏元，由相邻的若干个三角形构成，在宏元上采用 分片常应力模式和连续线性位移。因为应力参数可以在单元水平消去，所以计算 量并为增加。数值试验表明，通过调整组合参数，改进的版本可以达到比常应变 三角形元更高的精度。 力四边形宏元方法[7]进行了理论分析。该宏元采用连续分片线性位移插值逼近和 分片独立设计的5参数自平衡应力模式。我们证明了有限元解的存在唯一性，导 出了相应的误差估计。我们的分析也表明，在使用自平衡应力模式时，单元上采 用连续线性位移与采用连续等参双线性位移是等价的。 组合杂交格式和能量精度 考虑线弹性问题： Q (2．1) 一diva=7，盯=D6(ff)in (2．2) o‘iIn=T，iIro=0， OilF=铀=‘Ur0， 应变，D弹性模量，于为给定的体力，f为给定的表面力，rl为F作甩在Q ’ 的部分边界． 维逼近空间，使得 研c址#peHl㈣2 l：ilro=O}， 四川大学硕士学位论文 叫(足)：=span{bubblefunctions} rE矗 KEn 其中Sobolev空间 ‘bubble’函数表示单元K上的非协调内部位移， divr=(OrlI／Oxl+Orl2／Ox2，aq2／Oxl+Or22／叙2)2． 组合杂交变分原理可表述为[11,12,14】： (2·3) n擗(盯6，孬6)_叫inf～su—pn兽(f，哥)， 其中组合杂交能量泛函 K K 。 r。n8K 雠 r．D-lr]as， +口扣．s(i)一i1 i，∈U?表示单元非协调内部位移，f∈矿6为分片独立应力