数值模拟中岩体力学特性的工程与实验关系.docVIP

《数伉模拟屮岩体力学特性的工程与实验关系》 前言 随着更便宜、高效的软件和硬件的应用，数值模拟正越来越多的应用于岩土工程设计。然而， 数值模拟中关键的一步在于确定岩体的力学特性参数，即更精确的岩石硬度和强度特 性。 这篇文章给出了一系列用来岩体数伉模拟的硬度与强度特性。文章屮给出的实验模拟特性都 选自人量通用的模拟文献。具体来说,将从实验到现场使用挽算参数值的文献标记出来搜索, 其结果令人吃惊：在数以千记的数值模拟文献，只有几百篇提到了实验岩体特性，而且在这 些文章当巾，只有40篇似乎使用丫某种程度上的简化。那些使用换算参数的文章被做成曲 线图构成了本文的主要内容。岩石硬度特性与强度特性被分别分析，它们各自换算系数呈现 出有趣的差异。 方法 在审阅文献过程中，已经广泛的研究过历史实例和数值模拟反演分析例子。 每一篇审阅过的文章都依据实验室测定岩石特性和数值模拟岩体输入特性以及其他相 关的定量数据被收入数据库。大多数的文章都只提供丫不完整数据，它们不是忽略丫关键 参数就是，有些文章给岀了实验岩体特性，只有很少数文章解释了实验岩体特性根据岩体等 级适应岩体的过程。人们只能认为这些特性的选取与该模型的来源或模型建立者有关，例 如钢材从没有尺寸效应的环境中选取。儿乎没有岩石力学专家招绝承认即使是强度最 大的岩石在其岩体特性上也需要调整，因此本文给出的图表和数据都集屮在那些使 用换算系数的文章上。一些具有最有效的相关数值的文章也在本文的结尾处给出。 结果 图1为实验测出的杨氏模量与模型使用的杨氏模量的曲线。曲线中每个案例都根椐 其来源编号。如果将这些数据直线拟合可得到一个显著的趋势：模型使用的硬度平 均是实验测出硬度的0.469倍(图2)，类似的数据曲线也可用换算系数绘出，如图 3所示。在这里，低硬度岩石使用的换算系数有明显增高的趋势，在一些低硬度岩 石中甚至可见到很高的换算系数。 图4给出丫实验测出的单轴压缩强度与模型使用的单轴压缩强度的曲线。曲线中每 个案例都根据其来源编号。如果将这些数据直线拟合可得到一个显著的趋势：模型 使用的强度平均是实验测出强度的0.284倍(图5)，类似的数据曲线也可用换算系 数绘出，如图6所示。在这里，软弱岩石使用的换算系数也有明显增高的趋势。 图7给岀了抗拉强度的趋势，表明实验值大至是巾两个因素而减小的，图8则表明泊松比没 有明显的下降趋势。 4换算方法 许多 提出了一些实验室和现场之间的岩石特性关系，有些在关系屮还包括了岩 体分级。由Bieniawski提出的广泛应用的推导变形模量的公式1当RMR大于50时有 18.2%的误差。然而，当RMR小于或等于50时，Bieniawski公式因会导致变形模量 小于或等于0而变得不再适用。Serafim和Pereira [39]使用Bieniawski的岩体分级 系统(RMR)得到了用于完整RMR变化的另一表达式，公式(2)， Em=2RMR-湖(GPa)⑴ R淑-10 Em = 10^^ (GPa) (2) 阁9给出了根据硬度数据得出的所有绘阁表达式。双轴压缩也被用来比较这些数据。这就要求 RMR与实验弹性模量相关。一个简单的线性关系已经在所有典型特性中使用。(RMR = 0-100, E = 0-120 GPa.） Nicholson和Bieniawski[40],给出了一个换筇因素的经验公式，公式（3）。这个因数被计算出 來用以得到用来计算RMR和实验杨氏模量的岩体变形模量。 TOC \o 1-5 \h \z P D WD HYPERLINK \l bookmark10 \o Current Document \h RF=^il = o^Q28RMR2 + 0.9exp( ) (3) ￡int 22.82 Mitri etal. [33]使用下jfli的公式（4）来得到岩体的变形模量和汁算出胡克-布朗参数来表现工程情 况下的RMR RFnn0.5x 1 - COS(7TXint RF nn 0.5x 1 - COS(7TX int RMR 100 被绘出（图9）的公式3和公式4采用以上相似的数据处理方法。当RMR低于3?时，公式3可以 用以硬度的大规模换算。公式4可能更适合于在低RMR和低硬度范围吋的实际换算。尽管比 较的公式曲线和数据是由简单的RMR和实验硬度线性关系组成的，但是仍然可以淸楚的看出 所有公式在低RMR范围时硬度都减少了很多。 Matsui [9]提出了一个基于最小误差函数的直接方法，公式5。这个函数表示了n向位移W，（根 拿 据实际巷道围岩测量）和n向位移+ （由有限元分析获得）的最小二乘差换算。因为数值模 * 型中输岀的~值是依据在有限元分析屮假设的弹性模M参数值确定的，误差￡也对应函数 巾的参数（即￡ = 因此，矢量弹性模量的减