固体力学中的变分差分方法.pdfVIP

第 40 卷 　第 2 期 中山大学学报 ( 自然科学版) Vol40 　No2 　 2001 年 　3 月 ACTA 　SCIENTIARUM 　NATURALIUM 　UNIVERSITATIS 　SUNYATSENI Mar 　2001 　 固体力学中的变分差分方法 富明慧 , 刘祚秋 ( 中山大学应用力学与工程系 ,广东 广州 510275) 摘 　要 : 以弹性力学平面问题为例阐述变分差分方程的建立方法 、求解过程 , 介绍变分差分法在固体力 学各领域中的应用. 实例分析表明: 变分差分方法对复杂区域具有较强的适应性和较高的求解精度. 关键词 : 变分差分 ; 复杂区域 ; 迭代格式 O2422 　　 A 　　 (2001) 中图分类号 : 文献标识码 : 文章编号 : 　　变分差分法 (又称能量差分法) 是利用变分 能, U 为应变能, V 为外力势能, Cij kl , ( i , j , k , l = 原理建立差分方程的一种数值方法[ 1～3 ] , 它由于 1, 2) 为弹性张量,f i 和p i 分别为体力和面力分量. 具有简便易行 、所需存贮量小且精度高等优点 , 本文采用求和约定. 因而在 固体力学 的许多领域获得 了广泛 的应 首先考虑区域 A 为矩形的情况. 利用均匀 用[4～8 ] . 变分差分法与有限元法各有优劣 , 可互 网格将区域 A 离散, 用 ( xI , yJ) 表示网格上第 为补充. 对于条件数较大 (离散后差分方程极度 I 行 、第 J 列结点的坐标. 则 ( ) ( ) 病态) 的问题 , 直接解法的适应性较差 , 很难保 xI = x 0 + I - 1 h1 , 0 ≤I ≤N 1 , 证数值解的精度 ; 然而高效迭代法在一定程度上 ( ) ( ) yJ = y 0 + J - 1 h2 , 0 ≤J ≤N 2 可以弥补这一不足. 而就构造高效迭代法而言 , h1 , h2 分别是 x 和y 方向的步长. 变分差分法较有限元法更有优势. 固体力学中的 许多问题由于几何或者物理上的因素 , 在离散后 差分方程具有较大的条件数 (如裂纹尖端 、板壳 的三维应力分析以及具有强各向异性的复合材料 ) 结构等 . 对于此类 问题