数值流形方法对象设计.pdfVIP

第 18 卷第 1 期 计 算 力 学 学 报 V o l. 18 N o. 1 2001 年 2 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 2001 文章编号:(2001) 数值流形方法的对象设计 1 2 1 张湘伟 , 　蔡永昌 , 　骆少明 ( 1. 汕头大学工学院, 广东 汕头 515063; 　2 重庆大学 建筑工程学院, 重庆 400044) 摘　要: 研究了数值流形方法的对象设计方法和组织方式, 为流形方法的理论研究和向 三维问题的扩展打下良好的基础。研究发现数值流形方法具有编程和前后处理简单的 特点, 且仅用三角形流形单元和一阶近似的覆盖位移函数就可达到有限元多结点等参 元的求解精度, 具有深远的工程意义。计算结果表明, 数值解与理论解吻合。 关键词: 数值流形方法; 对象设计; 有限元法 中图分类号: T P 311　　文献标识码:A 1　引　　言 [1 ] 在有限覆盖系统基础上新发展起来的数值流形方法 , 在岩土工程界引起了广泛的关注。它以拓扑流形 为基础, 联合物理网格和数学网格, 通过在分析域内各物理覆盖上建立一般覆盖函数和加权求和形成总体位 移函数, 可以将解析法、有限元法和适用于块体分析的D D A 法包含并融合于统一的形式中, 是一种新的、较 通用的数值方法。 在石根华 和国内外一些学者大量研究工作[1, 2, 3, 4 ] 的基础上, 目前二维流形方法的理论已经基本形成。但 是作为一种新近提出的数值方法, 数值流形方法仍有大量的问题需要去研究和解决, 流形方法的理论也需要 得到进一步的完善和发展。而这些研究的进行必胯依靠具体的程序或软件来实现, 因此程序设计是流形方法 的一个很重要的环节。关于流形方法的程序设计, 国内外还没有专门的文献介绍。由于流形方法与有限元法有 着极大的相似性, 本文直接借鉴目前研究较多的面向对象有限元法程序设计理论, 首次将面向对象的思想引 入流形方法的程序设计, 利用面向对象程序设计方法的数据抽象、信息隐蔽机理和面向对象语言的继承性、封 装性等特性, 为流形方法的理论研究和向三维问题的扩展打下良好的基础。 本文的计算程序在物理覆盖上使用了常量和完全一阶近似的覆盖位移函数, 分析结果表明在数值流形 方法三角形流形单元上采用完全一阶近似的覆盖位移函数可以达到有限元四边形八结点等参元的求解精 度。这就意味着流形方法能够充分利用简单三角形单元的基本理论和其对任意复杂边界网格划分的适应性, 通过增加位移覆盖函数的阶数而不是单元的结点数来提高数值解的精度, 极大地简化了有限元的分析程序, 其工程实际意义是深远的。 2　 有限覆盖函数和位移函数 2 1　 矩阵形式的有限覆盖位移函数 流形方法具有两套网格, 即数学网格和物理网格。在将有限元网格转换为 流形方法的有限覆盖之后, 裂隙等可作为新的物理网格输入。物理网格是材料 体的边界和裂缝, 它与数学覆盖的交构成物理覆盖, 几个物理覆盖的交称为流 形单元。对无裂缝的连续体数值分析, 流形方法可以采用重合的数学、物理覆 盖。 收稿日期: ( ) 基金项目: 国家自然科学基金资助项目 ( ) 作者简介: 张湘伟 1950 , 男, 教授, 博士导师 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 1 期