单向压缩条件下圆形颗粒体含孔洞试样的力学模拟 Numerical simulation of mechanical behavior of a specimen composed of circular granular material.pdfVIP

第3期 水利水运工程学报 No．3 2010年9月 HYDRO-SCⅡeNCEANDENGINEERING Sep．2010 单向压缩条件下圆形颗粒体含孔洞试样的力学模拟 伍小林，王学滨，潘一山 (辽宁工程技术大学力学与工程学院，辽宁阜新 123000) 摘要：通过编程，用连续介质快速拉格朗日分析(FLAC)实现了圆形颗粒体在单向压缩条件下的力学行为模 拟．颗粒体试样中间含一圆形孔洞，颗粒体被离散成正方形网格，若2个颗粒体发生接触，则需要在颗粒体之间 预置界面．研究了施加的荷载对试样应力线、应变线及受拉破坏区域等的影响．结果表明，试样内部的最大压应 力约为所施加压应力值的5倍，这有别于连续介质理论的结果．孔洞顶部及底部各有1个三角形受拉区，内部可 见明显的拉裂纹，裂纹的长度随试样端部所施加应力的增加而增大． 关 键词：颗粒体；单向压缩；应力线；应变线；拉应力区 中图分类号：TU45 文献标志码：A 颗粒体材料(沙粒、道渣、矿石等)是典型的非连续介质，模拟非连续介质力学行为常采用离散元法程 序．在常规的离散元法程序中，假设颗粒是刚性的，若颗粒发生接触则可以传递力，从而形成力线或力网．但 常规的离散元法程序不能获得连续介质模型中常见的一些力学量的分布及演变规律，如应力分布、应变分 布、破坏单元数及能量释放量等．连续介质快速拉格朗日分析(FLAC)，可以有效地处理连续介质的变形、破 坏等问题，已广泛应用于岩土力学．FLAC可以直接获得或通过编程获得应力分布、应变分布、能量释放总 量、能量释放率、释放的能量分布及声发射数¨。’等与岩土材料破坏及失稳密切相关的力学量信息． 本文采用针对连续介质而开发的FLAC模拟典型的非连续介质，即颗粒堆积体材料在单向压缩条件下 的力学行为．首先通过编程生成圆形颗粒堆积体，颗粒半径服从均匀分布；然后，将颗粒体离散为正方形单元 构成的集合，若颗粒之间发生接触，界面同样需要离散为正方形单元构成的集合．文中颗粒体试样由3部分 构成，即颗粒、界面和孔隙．颗粒体试样生成之后，即可按照FLAC的常规方法进行加载，研究多种力学量的 分布及演变规律，从而实现采用FLAC来模拟颗粒体材料的变形、破坏过程． 1 计算模型、本构关系及计算方案 1．1模型的生成 在某一矩形区域内生成圆形颗粒堆积体，颗粒的堆积过程由下而上进行．为确保任一颗粒i在重力作用 下能保持静止，颗粒i至少要与其下方(或两侧)的2个颗粒发生接触，且颗粒i的中心应在其下方(或两侧) 2个颗粒的中心之间．此外，矩形区域的左边界、右边界及下边界均可视为半径为无限大的颗粒．在颗粒堆积 过程中，需要记录当前表层的各个颗粒，当任一新的颗粒被堆积上后，需修改当前表层各颗粒的记录．在堆积 新颗粒时，优先堆积在表层各颗粒中距矩形区域下边界最近的颗粒附近． 1．2试样的切割和挖孔 设矩形区域为A，在A区域内任取1个边长为Im的正方形区域曰，选择中心位于正方形区域B之内的 收稿日期：2009—09—02 厅2009年度高等学校科研项目(2009A322) 作者简介：伍小林(1987一)，男，四川宜宾人，硕士研究生，主要从事非均质材料(岩石、混凝土及金属等)变形、破坏及稳 定性研究．通讯作者：王学滨(E—mail：wxbbb@263．net) 万方数据 第3期 伍小林，等：单向压缩条件下圆形颗粒体含孔洞试样的力学模拟 4l 颗粒集合体作为试样．因颗粒体的强度高，不易在切割机上被切成两半，故在切割试样的过程中，不允许1个 颗粒被切成两半．如果有些颗粒的中心位于区域B之内，但 仅与B内的其它颗粒有1个接触点，则需剔除这些颗粒．本 文计算方案中，区域曰内的颗粒数为2680，最大的颗粒半径 为1．250cm，最小的颗粒半径为0．625cm．编程中假设颗粒 半径的分布服从均匀分布，颗粒半径的统计结果见图1． 切割好试样后，以试样的中心为圆心挖1个半径为 颗粒半径／cm 0．2m的圆孔(见图2)．在图