收稿日期:2002-09-04;修改稿收到日期:2 00 3-04-16.
基金项目:国家自然科学基金重点项目(50135030);面上项
目(10072026)和江苏省自然科学基金
(
BK
2002090)资助项目 .
作者简介:郑世杰 (1968-),男, 副 教 授.
第21卷第3期
2004 年 6 月
计算力学学报
Chinese Journal of Computational M echanics
Vol
.21,
No
.3
Ju n e
2004
文章编号:1007-4708(2004)03-0308-06
用于超弹性分析的高效杂交应变
——
EAS
固体壳单元的研究
郑世杰
1
, 佘锦炎
2
(1.南京航空航天大学 智能材料与结构航空科技重点实验室, 江苏 南京 210016;2.香港大学 机械工程系)
摘 要:提出了一个用于橡胶材料分析的弱变分方程,基于一个可压缩
Neo
-
H ookean
模型,将
EA S
模式和杂交
元法有机地结合起来,推导了一个高效、稳定的杂交应变—— EA S 固体壳单元,通过巧妙地选择应变插值函数
和本文中提出的一个精化措施,克服了橡胶材料所表现出的大应变超弹性本构为杂交元正交化法的实施所带来
的困难,保证了整个单元列式都仅采用低阶高斯积分,显著提高了计算效率,确保橡胶材料不可压缩性计算的顺
利进行,并克服了固体壳元的厚度自锁问题。
关键词:超弹性;杂交固体壳;厚度自锁;变分原理;稳定
中图分类号:
O
316 文献标识码:
A
1 引 言
在现代工业中,橡胶材料获得了广泛的应用,
但由于橡胶材料所表现出的不可压缩性和大应变
超弹性,关于它的本构模型和数值分析的研究较之
其它材料向结构分析人员提出了更严峻的挑
战
[1-5 ]
,伴随大应变而来的是壳体厚度方向上的显
著变化和有限转动。基于平面应力假设的退化壳
元,忽略了厚度方向的变形,显然难以满足橡胶材
料的大变形超弹性分析的要求
[2]
。因此研究者们转
而发展放弃平面应力假设的、考虑厚度方向变形的
壳单元
[2, 3]
,但列式过程曲折、复杂。最近新兴的固
体壳元
[6-8]
自然地考虑了壳体厚度方向的变形,用
于壳体的大转动分析时可避免退化壳元的繁杂的
矢量转动计算
[9]
,并且可采用大的迭代步,且用于
同时需要固体和壳单元进行分析的结构时,不需要
为处理平动和转动自由度的连接而使用过渡性单
元或定义额外的代数约束方程。毋庸置疑,使用固
体壳元进行复杂板、壳结构的非线性分析是结构分
析人员的一种理想选择,但随之而来的是固体壳元
的厚度自锁问题(即发生弯曲变形时沿厚度方向线
性分布的位移插值函数和沿厚度方向线性分布的
应力的矛盾)。由于橡胶材料的本构方程是应变的
函数,通过修改本构方程克服厚度自锁是不可能
的,只有修改位移插值函数或应变模式。但由于超
弹性问题迭代求解的巨大计算量,现有的超弹性壳
元研究基本上集中于低阶单元
[4 ,5]
,目前还未见到
用于橡胶材料超弹性大应变分析的高阶固体壳元
和杂交元。
正交化方法
[10]
是近年来多变量有限元研究领
域取得的一个主要成果,它能明显地提高低、高阶
单元的计算效率。但橡胶材料的本构方程大多是右
C auchy-G reen 张量的不变量或主伸长量的隐函数
形式
[1]
,即在不同的高斯点取值不同,为正交化方
法的实现,也为多变量有限元赖以生存的基础-变
分泛函的建立带来了巨大的困难,这也是目前难以
见到用于橡胶材料超弹性大应变分析的杂交单元
的根本原因。本文建立了一个新的修正的广义变分
原理,引入 E A S 模式克服固体壳元的厚度自锁问
题,通过精心选择的插值模式保证了正交化方法的
完全实现,并在高阶稳定阵的计算中采用了一个精
化措施,保证了整个单元列式都仅采用低阶高斯积
分,从而既克服了橡胶材料所固有的不可压缩性自
锁问题,又显著提高了计算效率。
2 变分原理
∏
e
HW
=〈W (ε)- σ
T
(ε-ε
-
)〉- (U )
T
F (1)
其中 〈
°
〉=
∫
V
e
(
°
)dv
e
,W(ε) 为应变能函数,ε为当
前构型的独立假定的 G reen 应变矢量,ε
-
为与位移