粘弹性围岩力学参数反分析的一种数值法岩石力学与工程学报.PDF

第16 卷　第 5 期 ( ) 岩石力学与工程学报 16 5 : 478～ 482 1997 年 10 月 Ch inese J ou rna l of R ock M echan ics and E ng ineering O c t. , 1997 粘弹性围岩力学参数反分析的一种数值法 1 2 1 1 朱浮声 　薛　琳 　李　宏 　李国辉 ( 1 ) (2 ) 东北大学土木系　沈阳　110006 　 青岛建筑工程学院　青岛　266033 摘要　依据蠕变柔量和广义蠕变柔量概念, 提出了由实测位移反演原岩应力场与粘弹性力学 参数的边界元算法, 实现了粘弹性围岩力学参数反演的半解析化。 关键词　蠕变柔量, 广义蠕变柔量, 位移反分析, 边界元法 1　前言 为了研究软弱地层中岩土工程围岩的时变特性, 可以采用粘弹性力学模型。在正确进 行岩土工程分析和设计前, 必须首先确定合理的围岩力学参数和原岩应力场。以现场量测 信息为基础发展起来的反分析方法是提供这些实用参数的重要手段。 由于岩土工程问题的复杂性, 反分析方法使用的主要工具是工程数值方法。自从 [ 1 ] Saku rai 提出弹性岩体等效弹性模量和原岩应力场的有限单元位移反分析算法以来 , 位 移反分析在弹塑性、粘弹性和粘弹塑性岩体力学参数反演计算中取得了显著进展[ 2～ 4 ] , 并 对反演建模、模型识别、模型可信度与检验理论等进行了广泛研究[ 5～ 7 ] 。对于粘弹性反分 [ 3 ] 析问题, 王芝银等提出了一种有限元算法 , 李云鹏等则给出了基于“综合模量”概念的边 [ 8 ] 界元算法 。 作者发现, 在上述边界元位移反分析法中, 由于“综合模量”反映的是各粘弹性参数的 综合特征, 为了得到给定力学模型的粘弹性参数, 必须采用参数回归或参数优化等分离方 ( ) 法进行优化计算。这样不仅增加了计算工作量 机时 和计算机源程序的复杂性, 而且, 产 生了新的计算误差。事实上, 如果将粘弹性边界元法的基本解表示为蠕变柔量和广义蠕变 柔量的函数, 由此构成的位移反分析边界元法不仅可大大提高计算精度, 减少编程困难, 并且可以直接用于应力反分析法中。 2　粘弹性边值问题的蠕变柔量边界元法 考虑平面应变问题。本文作者提出的粘弹性边值问题的蠕变柔量边界元的实质是, 将 粘弹性平面问题的基本解表示为蠕变柔量和广义蠕变柔量的解析函数, 对于几种常见简化 粘弹性力学模型,