弹性近不可压缩问题的精化元法.pdfVIP

大连理工大学学报 Vol. 4. 7 ,No.4 第47垂第4 期 2001 年 7 月 Journal of Dalian University of Technology Jul. 2 0 0 7 文章编号， 1000-8608(2007)04-0479-04 弹性近不可压缩问题的精化元法 卓 陈万吉\冀 宾，赵 (大连理工大学工业装备结构骨折回辈重点实验室，辽宁大连 116021 ) 摘要对于弹性问题，材料近不可压缩可引起哺限元法的悻权问锁为解决此问题，以精化 元浩为基础，将单元应变正主舟解为常应变和高阶应变，其中常应变可以保证llJ:鼓，对于近平 可压缩问题只需起略高阶应变中的体应变，从而避免单元体租平可压缩向锁 按照上述方挂 静改了平面八节点等#元(IQ的，井用于平面应变厚壁简计算，通过与 AßAQUS 革统教件的 八节点等#元计算结果对比，表明 IQs 单兀用于可压和近平可压缩问题都有效. 关键词:不可压缩;精化元，厚壁筒;平面八节点等参元 中图分类号 0343 文献标iR码， A 工作，其中精化不协调元可以直接用于改进任何 。引言 协调和不协调元忡8] 用有限元法解弹性不可压缩问题，在材料不 按精化元法，先假定单元位移函数为 可压缩或近不可压缩情况下，ν → 0.5 ，→∞，体 u = u + p民(1) 积应变ι 必应当为零，但是常规有限元法不保证 式中:U 为任意已知的协调元或不协调的单元函 体积应变E 为零，这时将会出现体积闭锁现象，求 数;P 满足DTp=l， 设D 为微分算子向量，对于二 解所得位移偏小，甚至趋向于 0 咿 {if éf if\ T 目前，国内外已提出不少解决体积闭锁问题 维问题，D = [一 ÷ 一一|风为待定常 \ Jx2 ay Jx句1 的方法，例如减缩积分或减缩/选择积分法口、，]、 数. 罚函数法I气混合/杂交法。]等，前者应用较多. 精化不协调元法基于如下泛函· 但减缩积分会使单元刚度阵秩降低，常引起多余 零能模式为此，Belytschko 等提出 B-bar 方 1(;:p = 早(jl÷(DV 阳D叫w一 法[气可以避免多余零能模式，但这种方法常常不 能保证通过分片检验.以上解决体积闭锁的方法 j巾每一叫-w ω 都是以 ι=0 为约束条件.本文采用一种新的精