冰碛土体强度值的合理选取.docxVIP

冰碛土体强度值的合理选取 1 利用计算机模拟试验求取动力参数的方法 在计算和分析岩石边坡的稳定性时，岩石强度参数是一个非常重要的问题，但其合理性一直是岩石工程领域的难题。在传统的试验方法中,人们曾试图从取样方法、试验方法以及测量精度等多方面进行改进,以期获得相对精确的取值,并且取得了大量成果。不过无法回避的事实是:(1)取样时,取不扰动样难。以目前的科学和试验水平,在取样过程中想避免应力释放和土体结构的扰动是不可能的。(2)试样不具备代表性。相对整个工程对象而言,试样的选取数量微不足道,以工程对象极小部分的强度特性表征整个工程对象的强度特性,明显不合理。扩大取样数量和取样地点,需耗费大量的人力、物力和财力,且不能从根本上解决这个问题。(3)试样尺寸太小。囿于试验仪器规格和试验条件,试样尺寸往往偏小,因此,必然要舍弃许多地质信息,试验结果自然无法反映岩、土体的宏观力学性质。 曾经有人采用现场调查和室内试验相结合的方法——地质强度指标(GSI)法来评估岩体力学强度,并被广泛接受。GSI法是一个半经验、半试验方法,它通过现场调查对工程岩体进行分级,依据广义HOEK-BROWN准则,采用室内试验得出的岩石强度参数求取岩体的强度参数。由于岩体分级主要依据经验判断,难以避免主观性和随意性,导致获取的岩体强度参数因人而异。 近十年来,随着计算机和计算技术的进步,计算机“模拟试验”逐渐成为热门研究课题。在三峡水利枢纽工程几个关键问题的应用基础研究中,董学晟等采用计算机模拟岩体力学试验,求取宏观的岩体力学参数,在国内较早地将计算机模拟试验方法推向工程应用;王学滨等采用直径为5 cm“花岗岩试件”模拟了抗压应力-应变全过程,采用5 cm×1 cm试件研究了变形局部化的剪切带图案和演化规律;郭志华等进行了水平层状2 m×2 m×2 m岩体“三轴压缩试验”,分析了岩层倾角对变形参数的影响。混凝土骨料形状、分布以及骨料剖分方法方面的研究也取得了可喜进展,他们采用的基本方法是:骨料形状和大小采用随机模拟方法,分布依据“互不靠”规则,能较好地反映搅拌后的混凝土结构特性。冰碛土沉积过程是极复杂的动力过程,至今无法用非线性方程来描述。面对这样一个纯天然沉积、粗细悬殊且具有砾石随机“聚团”结构特性的全级配冰碛土体,难以估计“砾石空间分布”的概率密度函数,很难实现冰碛土结构的随机模拟,由此导致三轴模拟试验成为“无米之炊”,无法进行。 元胞自动机(CA)是空间、时间和状态都离散的动力系统,无需构造复杂的数学模型,只需由一些极简单的局部相互作用规则来驱动整个系统的演化。元胞的状态变化只取决于自身及其邻居的状态,而对元胞自动机整个空间来说,元胞更新是彼此独立的,作用规则是通用的。因此,本研究首先采用元胞自动机模型来模拟冰碛土的沉积结构,为三轴数值模拟试验“制备”任意尺寸的试样;然后,借助FLAC3D实现冰碛土的三轴试验数值模拟,依据常规三轴试验的力学公式求取冰碛土的强度参数;最后,藉此提出笔者心中解决岩、土体材料强度问题的“虚拟岩土实验室”构想。 2 冰字土结构的元胞自动机模拟 通过试验方法获取工程对象的强度参数指标,首先要确保试样的代表性。现实困难在于实验室和现场试验方法都不能从数量、尺寸和原状性等方面满足这个要求,GSI方法则遭遇岩体结构如何精确量化的瓶颈。数值模拟方法由于具有可控性、无破坏性、安全性和可重复性等特点,可有效弥补这些缺点。 常规的有限元方法通常将岩土体当作各向同性的均质体处理。不过对于冰碛土,这样处理显然忽略了冰碛土中黏土和砾石这两种材料各自对冰碛土体强度指标作出的贡献,而且冰碛土自身组成和结构有其鲜明的特点,与混凝土有本质的区别。 冰碛土结构的元胞自动机模拟方法,是指在指定的二维空间内,以方形网格形式,根据砾石在冰碛土沉积过程的动力特性,以及工程现场勘察所描述和测定的砾石级配、分布以及结构特性,拟定砾石“沉积”的演化规则,通过元胞自动机演化,随机生成冰碛土的结构图。它能从总体上、近似现场观察结果中描述冰碛土结构特性,充分反映砾石分布的随机性,以及砾石沉积过程“聚团”的疑似成层性。虽不能具体表示出每块砾石的“可能形状”,但因“砾石土”的特定组成所决定,并不会因此而影响模拟的最终结果。依据这一方法,开发了能量化描述砾石沉积结构特性的元胞自动机程序GTSM。程序可在试件尺寸、单位长度、初始条件以及砾石含量(确定性的和不确定性的)等多方面为用户提供选择来“制备”三轴试验试样。图1为砾石含量在10%~60%区间内,随机选择砾石含量生成的试件尺寸为0.2 m×10 m×16 m时的2幅冰碛土结构图。从其外观上看,砾石“团聚”几何形态与现场观察结果颇为相似,表明采用该程序“制备”试样是十分合适的。 3 试验原理、方法和结果 3.1 本构模型的建立 三轴数值模拟