一种提高数值解精度的方法探讨.pdfVIP

第21 卷 第1 期 岩石力学与工程学报 21(1) ：136～139 2002 年1 月 Chinese Journal of Rock Mechanics and Engineering Jan .，2002 一种提高数值解精度的方法探讨 1，2 2 2 1 闵 弘 章 光 王水林 王元汉 (1 华中科技大学土木工程学院 武汉 430074) (2 中国科学院武汉岩土力学研究所 武汉 430071) 的数学覆盖和物理覆盖。 1 前 言 数学覆盖由覆盖整个材料体的许多有限覆盖组 成。常规的网格和域，诸如规则的格子、有限元的 数值流形方法是石根华博士提出的一种涵盖有 网格或级数的收敛域，均能转换为数学覆盖。 限元与离散元方法的全新的数值方法[1~8] 。该方 如果裂缝或块体的边界把一个数学覆盖分成两 法采用有限覆盖体系(数学覆盖和物理覆盖)去覆盖 个或更多的完全不连续的区域，这些区域定义为物 材料区域。在对材料区域完成覆盖之后，流形单元 理覆盖。因此物理覆盖是不连续缝对数学覆盖的再 上的位移函数可根据各物理覆盖上的位移函数加权 剖分。 平均得到，由此可以求出系统内的各项势能，再对 2.1 数值流形方法的有限元覆盖 系统的势能求极值即可得到系统的控制方程。 2.1.1 由有限元网格构成的有限覆盖 本文将数值流形方法应用于p 型自适应分析， 有限元网格可定义流形方法的有限覆盖。如考 以探讨一种更有效地提高数值解精度的新途径。 虑任一节点，含有这一节点的所有单元形成一个数 学覆盖。 2 数值流形方法简介 图1 采用具有4 个节点的4 个四边形单元的有 限元网格覆盖材料区域，区域内有一裂隙尖灭其中。 数值流形方法是利用现代数学“流形”(包括拓 节点5 的数学覆盖V5 由4 个四边形组成，即2365 ， 朴流形和微分流形)分析中的技术建立起来的。数值 1254，5689 和5874 。节点1 的数学覆盖V1 仅由四 流形和传统的微分流形有以下几点区别：(1) 微分 边形 1254 组成。图1 的数学覆盖见表1。 流形的总体函数是高度可微的，且完全可被定义而 与覆盖无关；(2) 数值流形的总体函数是在覆盖基 7 8 9 7 础上定义的，分段可微，在接触面上几乎都是不连 1 81