循环荷载作用下砂土液化特性的非圆颗粒数值模拟 Numerical simulation of liquefaction of sands with non circular particles under cyclic loading.pdfVIP

水 利 学 报 2008年9月 SHUlUXUEBAO 第39卷第9期 文章编号：0559．9350(2008)09—1074．09 循环荷载作用下砂土液化特性的非圆颗粒数值模拟 史旦达1，周 健2”，刘文白1，贾敏才2’3 (1．上海海事大学海洋环境与工程学院，上海200135；2．同济大学地下建筑与工程系，上海200092； 3．同济大学岩土及地下工程教育部重点实验室，上海200092) 摘要：针对纯圆颗粒模拟的缺陷，利用颗粒流计算程序的团颗粒方法开发了形状近似实际角粒砂的角粒团颗粒，研 究了颗粒形状变化对数值试样液化特性的影响及其细观机理。利用自行开发的随机数制样方法制备出多种形状 颗粒随机混合的仿真试样，分析了振动液化过程中细观组构的演化规律，并与砂土可视化模型试验结果进行了对 比。结果表明，颗粒棱角度越大，试样抗液化能力越强，其在细观机理上与初始平均接触数大小有关；随机混合仿 真试样可以从定性角度反映实际砂土振动液化过程中的细观组构演化规律，为今后进一步的定量研究奠定了基 础。 关键词：颗粒形状；角粒团颗粒；平均接触数；颗粒棱角度；细观组构 中图分类号：TU441 文献标识码：A 1研究背景 实际砂土颗粒一般都具有不同程度的角粒特征，而颗粒形状对砂土力学性质的影响也一直是土力 等托1针对4种不同形状的砂土分别进行了三轴剪切试验，研究了颗粒形状对抗剪强度的影响。Vaid 等¨1针对相同级配的角粒砂和圆粒砂，研究了颗粒形状对砂土液化特性的影响。Ishibashi等¨1利用循环 扭剪试验，研究了颗粒圆度指标对砂土抗液化强度的影响。然而，试验研究一方面依赖于高精度试验设 备和细观量测技术的研发，另一方面由于实际砂土颗粒形状各异且外形复杂，各个形状指标相互影响， 很难对单一形状指标展开定量研究。以颗粒离散元为代表的数值分析方法可以在保证初始制样“同一 性”的前提下，根据不同的研究目的，快速生成单一形状或多种形状混合的仿真试样，为定量研究颗粒形 状对砂土力学性质的影响提供了条件。例如，Williams等b1针对分别由圆盘颗粒、椭圆颗粒和多边形颗 deformationofdisks)方法生成二维圆盘组合颗粒，研究了颗粒形状对剪切带形成的影响。 analysis 单元是二维圆盘和三维圆球单元，采用纯圆颗粒模拟存在颗粒过度旋转、组构信息单一等突出问题]， 因此，PFC2阴。程序中提供了“CLUMP”方法(译为“团颗粒”方法)可供用户对非圆组合颗粒进行二次开 发。贾敏才№’利用“团颗粒”方法构造三种不同形状的非圆颗粒，研究了振冲激振作用下颗粒排列定向 的演化规律。Sukulnaran等阳1将“团颗粒”方法与数字图像技术结合起来，将实际砂颗粒数字图像导人颗 粒流程序，为定量模拟砂颗粒提供了启示。关于“团颗粒”方法的具体原理可详见文献[10]。 收稿13期：2008—05—08 面上项目(2008C23023) 作者简介：史旦达(1979一)。男，浙江舟山人，博士，讲师，主要从事细观岩土力学、海洋岩土力学、港121软基处理等方面的研究。 E·mail：shidanda@163．com 一1074— 万方数据 本文利用“团颗粒”方法开发形状近似实际角粒砂的“角粒团颗粒”，研究颗粒形状变化对数值试样 液化特性的影响，初步探讨其细观力学机理，利用自行开发的随机数制样方法制备多种形状颗粒随机混 合的仿真试样，分析振动液化过程中细观组构的演化规律，并与砂土可视化模型试验结果进行对比。 2角粒团颗粒开发及数值试样制备 2．1 角粒团颗粒开发PFC20通过圆盘颗粒的团聚构造“团颗粒”，“团颗粒”在振动过程中不发生破碎。 本文生成“角粒团颗粒”时，圆盘颗粒采用同一粒径，且彼此不叠合。角粒团颗粒根据采用的圆盘颗粒数 量进行命名，如由3个圆盘构成称为“三圆团颗粒”，用c3表示；由4个圆盘构成称为“四圆团颗粒”，用 C4表示；依次类推。生成的角粒团颗粒如图1所示，其中c1表示纯圆颗粒。 Cl C3 C4 C5 C6