(6432)--拉拔实验中土工格栅的变形与受力分析.pdfVIP

以土工格栅加筋砂垫层体的拉拔实验为研究对象，采用 PFC3D 程序建立其颗粒流数值模型。以平行黏结 的“ball”单元模拟土工格栅，通过分析拉拔过程中土工格栅颗粒的位移、速率以及黏结力揭示土工格栅的变形与 受力特性。拉拔过程中格栅的变形由拉拔端向末端发展，同时格栅网格产生变形导致网格中间变形大于两侧的变 形。格栅的变形可以分为两个阶段，分别为弹性变形阶段和弹性变形恢复阶段。由于土工格栅横肋的侧阻力及摩 擦力作用，纵肋的受力在纵/横肋节点处发生突变。 土工格栅；颗粒流；数值模拟；拉拔实验 U416.12 A DOI ：10.15913/j.cnki.kjycx.2019.13.004 粒长轴方向随机原则。④土工格栅模型。土工格栅采用具有 岩土工程数值计算主要包括连续与非连续介质力学方 平行黏结参数的“ball”单元模拟。 法，其中，连续介质力学方法主要有限差分法、有限元法、 加筋砂垫层体拉拔模型（砂垫层厚度为 10 cm，在土工 边界元法等，而非连续介质力学方法主要有颗粒离散元法、 格栅上下对称铺设）如图 1 （b ）所示。土工格栅数值模型 块体离散元法等。土体是一种非连续介质，颗粒流方法在模 如图 2 所示。 拟岩土工程时具有一定的优势，主要考虑了土体的碎散型， 将土体模拟为由众多颗粒组成的集合体，对离散颗粒及其接 触设置细观参数来反映实际土体宏观现象。 Itasca 公司开发的 PFC （Particle Flow Code ）程序作为 目前较为主流的颗粒流程序，其程序主要包含球体单元、墙 体单元及接触与黏结等基本元素。球体单元是程序内离散单 元的主要构成元素，墙体单元可作为边界约束等，设置合理 的接触与黏结类型及参数反映土体宏观问题。 （a ）加筋砂垫层模型颗粒直径分布图 加筋砂垫层是提高地基承载力[1] 、改善软土地基沉降[2] 的有效手段，本文采用 PFC3D 建立土工格栅加筋砂垫层体的 拉拔实验离散元模型，借助离散元模型在细观参数分析上的 优势，分析土工格栅在拉拔过程中的变形与受力特性。 2.1 拉拔实验模型建立 拉拔实验的颗粒流数值模型建立步骤如下：①实验箱模 型。采用“wall”单元建立拉拔实验模型箱，模型箱尺寸 为 300 mm×300 mm×300 mm。②黏土模型。黏土采用“ball” （b ）数值模型图 单元模拟，为了保证模型的真实性同时考虑模型计算效率， 图 1 拉拔实验数值模型图 采用分层填充，在土工格栅附近的土体采用较小直径球体， 其他区域适当放大颗粒直径，以降低模型颗粒总数，提高计 算效率，模型各区域颗粒尺寸如图 1 （a ）所示。③砂土模型。 采用非圆颗粒相比用圆形或球形颗粒模拟砂土更能真实反 映砂土的工程特性[3-4] 。因此，将模型内砂垫层部分的“ball” 置换成“clump” ，置换过程遵循体积相对、质量相等以及颗 图 2 土工格栅数值模型 科技与创新┃Science and Technology Innovation 2019 年 第 13 期 2.2 模型参数的选取 表 2 所示。墙体的参数主要包括刚度以及摩擦系数，采用光 通过反复调试，在保证数值模型有效性的同时提高模型 滑墙体，摩擦系数设为 0 ，为防止与墙体接触的颗粒溢出模 计算效率。筋材颗粒细观参数及模型土体参数分别如表 1 和 型，墙体刚度设为颗粒刚度的 10 倍。 表 1 筋材颗粒细观参数 半径/mm 颗粒比例/ （kg/m3 ）