Hybrid2Trefftz有限元法的研究进展.pdfVIP

第 28 卷 　第 1 期 力　学　进　展 Vol. 28 　No. 1 1998 年 2 月 25 日 ADVANCES IN M ECHAN ICS Feb. 25 , 　1998 HybridTrefftz 有限元法的研究进展 秦庆华 清华大学工程力学系 ,北京　100084 摘 　要 　系统概述了 HybridTrefftz 有限元法及其在弹性力学中的应用. 该单元模型由于在插值函 数上的灵活选择性使其比普通有限元能更有效地处理局部效应问题 ,如孔洞 ,集中荷载等. 通过适当 选择单元插值函数可构造出高精度的 p - 扩展元和一系列满足特殊条件的新单元 , 以在同等条件下 提高计算精度. 关键词 　HybridTrefftz 有限元 ,p - 扩展元 ,网线位移 ,局部效应 1 　引 　言 作为计算力学中的一种数值方法 , HybridTrefftz ( HT) 有限元法已开始引起人们的重视并 在工程中得到越来越广泛的应用. 该单元模型是基于在单元域内和边界上分别独立地假设两组 位移函数而建立的[ 1～ 3] ,这与唐立民等的“网线单元”有相似之处[ 4] . 其做法是 :首先寻求单元 ( 区域上的插值函数 ,这些函数应是全场方程的齐次解 ;其次是在单元边界上假设插值函数 这一 ) 函常形象地称为网线函数 ,使其满足单元间的位移协调条件. 最后用变分法或其它方法将二者 联系起来 ,使最终得到的刚度方程仍为以节点参数为未知数的线性方程. 该单元模型最早由Jirousek 和 Leon[ 1] 提出. 通过对文中的 24 自由度薄板单元的考察,他们 发现这种单元精度很高 ,且对网格畸变不敏感. 进一步研究还发现这种单元在处理孔洞、角点和 裂纹等带局部效应的问题时非常有效[ 3 , 5 , 6] . 对于带局部效应的特殊单元 ,该模型的做法是构造 一种特殊的域内插值函数 ,使它既满足全场方程 ,又能近似反映其局部奇异性. 这种函数可从许 ( [ 7] [ 3] [ 6] ) 多文献中找到 如各种角点 、圆孔 、椭园孔 . 于是 ,在处理这类由局部效应引起的应力集 中及其它奇异性时 ,不需另外细分网格 ,只需在插值函数上作调整 ,就能得到较理想的精度. 到 目前为止, 这一有限元模型已应用在平面弹性力学问题[ 2 , 8] 、薄板[ 3 , 9 , 10] 、中厚板[ 11～ 13] 、厚 板[ 14] 、薄壳[ 15] 、三维固体[ 16] 、轴对称问题[ 17] 、泊松方程[ 18] 、动态热传导[ 19] 、板弯曲动力问 题[ 20] 和非线性问题[ 21～ 23] . 2 　方法概述 与普通有限元法不同的是 , HT 有限元法是分别在单元域内和边界上独立地对位移插值. 在一给定单元 e 的区域 Ω 上 ,位移插值为 e 　收稿日期 : ,修回日期 : ·71 · © 1995-2007 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. ( Ω ) ( ) u = u+ N c = u + Nc 　 在