弹性力学轴对称问题的有限元线法Ξ.PDF

第 13 卷　第 3 期 应 用 力 学 学 报 . 13　 . 3 V o l N o 1996 年 9 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Sep. 1996 弹性力学轴对称问题的有限元线法 牛忠荣　于红光　李景高 (合肥工业大学　合肥　230009) 摘　　要 给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了 24 条结线的等参数单元, 沿结线方向的两点边值问题采用插值矩阵法解之。算例表明, 本法具有良好的收敛性和较高的计算精度。 关键词: 有限元线法; 轴对称; 线弹性; 插值矩阵法 1　引　言 众所周知, 有限元法以分域多项式插值对结构采用完全离散化, 使之几乎能应付各种复杂 边界问题, 且具有通用性。使力学及场域问题的计算得到飞跃发展。有限元法放弃利用所求问 题的解析信息, 招致计算量大, 精度较差。利用解析良好的函数配合结构的部分离散已建立若 干半解析方法, 如有限条法、边界元法等, 对某些问题, 其效率优于有限元法, 但适用范围无法 与有限元法比拟。近年来, 一些学者[ 1 ] [ 2 ] 基于常微分方程求解器( ) , 提出了有限元 OD E So lver [ 3 ] 线法解二维偏微分方程。该法精度较高, 适用于复杂区域, 现已发开了若干线弹性单元 , 显出 广阔的应用前景。 本文给出解弹性力学空间轴对称问题的有限线元。采用笔者建立插值矩阵法的OD E 求 解器 IVMM [ 4 ] [ 5 ]。 2　轴对称问题有限元线法 ( ) 空间轴对称结构用柱坐标 r, z 描述, 该结构的轴对称线元为一关于对称轴z 的旋转环形 ( ) ( ) 体, 单元子午面为 rz 平面, 如图 1 示, 由P 条结线L 1 , …,L P 和两端边S 1 , S 2 构成, 与图 2 母 单元的映射关系 P P r = N ( ) ( ) = N ( ) ( ) ( ) i ri ; 　Z i Z i 1 i= 1 i= 1 本文得到机械部科学基金资助 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 3 期