带有沙漏控制的相对自由度壳元-工程力学-清华大学.PDFVIP

第19 卷第3 期 工 程 力 学 Vol.19 No.3 2002 年 6 月 ENGINEERING MECHANICS June 2002 文章编号 1000-4750(2002)03-122-06 带有沙漏控制的相对自由度壳元 1 1 2 1 陈丽华 程建钢 黄文彬 姚振汉 (1. 清华大学工程力学系, 北京 100084; 2. 中国农业大学(东校区)应用力学教研室, 北京 100083) 摘 要 针对一点积分的八节点相对自由度壳单元存在的沙漏现象 提出采用拟应变法解决该问题的方法 并对 锁死问题进行研究 给出了带有沙漏控制的八节点相对自由度壳元内的坐标 位移插值公式 推导了拟应变的表 达式 通过Hu— Washizu 变分原理 建立了有限元求解方程 利用Wilson 非协调位移模式 单元的计算精度得 到了明显改善 算例表明 基于八节点相对自由度壳单元 本文给出的沙漏控制算法能够有效的解决线性静力问 题 并且具有较高的计算精度 关键词 相对自由度壳单元 拟应变 沙漏 锁死 中图分类号 0347.3 文献标识码 A 1 引言 量 具有显著的简便性和通用性 针对这一特点 本文对带有沙漏控制的一点积分的相对自由度壳 在有限元分析中 高阶插值单元通常比线性插 元进行了研究 值单元更为精确 但是从计算量的角度出发 考虑 到变量的个数和问题的规模 后者在实际应用中更 2 变分原理 为常见 线性静力系统的Hu— Washizu 变分原理[7] 在形成单元刚度矩阵和求解内力时 采用完全 T 积分方案能够保证结果的稳定性和收敛性 但是该 ä P(u,å ,ó ) V ä å ô d V (1) 方案的计算量相当大 而且对于不可压缩材料 可 V ä [ó T ( s u å )] d V ä W ext 0 能导致体积锁死 对于薄的受弯结构 也可能出现 ä 表示变分 u,å 和ó 分别是独立的位移场 拟应 剪切锁死 为了克服锁死 减少计算量 通常采用 变场和拟应力场 ô 是由本构关系ô Cå 得到的 缩减积分方案 但它可能导致结构的不稳定 或出 ext Cauchy 应力 s u 是位移梯度的对称部分 W 为 现沙漏现象(hourglass) 沙漏现象是指由于一点积 作用节点上外力所做的功 设拟应力场ó 满足正交 分而引起的伪零能模式的存在 导致单元产生 沙 条件[10] 漏 ( 古代一种计时器)形状的变形 从而使计算结 T ó ( u å ) 0 (2)