孔隙岩石损伤破坏行为数值模拟研究.pdfVIP

北京力学会第18届学术年会论文集：计算力学 孔隙岩石损伤破坏行为的数值模拟研究 张玉军 刘坚志 彭瑞东 张文博 (中国矿业大学(北京)力学与建筑工程学院，100083) 摘要：运用数值模拟的方法来研究孔隙岩石的损伤破坏行为是一种比较新颖的方法。本文 利用MATLAB程序生成了具有某种统计分布特征的随机数序列，并将其作为狭长椭圆孔隙 的长轴，短轴、倾角及位置坐标，建立了孔隙岩石的几何模型。通过选择合适的材料本构 模型，借助单元生死技术，模拟了孔隙岩石的损伤破坏过程。研究指出孔隙岩石的损伤破 坏行为不仅与孔隙率．孔隙位置等几何参数有关，也和基体的弹性模量、屈服强度等材料 属性有关。随着孔隙率增大，岩石细观单元损伤现象更加明显，岩石强度急剧下降，且损 伤阈值也有所下降。 关键词：岩石，孔隙，损伤，数值模拟，单元生死 一、引言 随着我国大规模经济建设的不断发展，采矿、交通运输、水利水电、城市建设中的岩 石工程规模和难度日益增大。研究揭示岩石损伤破坏行为的规律特点已成为岩石力学研 究的难点和热点问题。Chung和Waas／lJ用有限元对聚碳酸醋脂蜂窝材料在面内受压缩的 情况进行了数值模拟。鞠杨等12谴过砂岩及泡沫混凝土的CT扫描实验研究了岩石孔隙 的几何特征与分布规律。但目前有限元模拟研究还主要局限于对含有等轴圆孔的高孔隙 率岩石材料的研究【3】。较之圆形孔隙，这些具有狭长孔隙且孔隙率较低的岩石材料给数 值模拟带来一些新的困难，因此有必要进一步对此展开研究。本文将借助ANSYS有限 元软件，通过单元生死技术模拟孔隙岩石的损伤破坏行为，研究孔隙结构对岩石损伤破 坏行为的影响。 孔隙岩石损伤破坏模拟 2．1孔隙岩石的二维建模 根据以前的实验结果，本文构建了二维岩石细观单元，单元边长为100I．tm，其中含 有随机分布的狭长孔隙，以长短轴比为lO的椭圆孔表示，孔隙位置采用均匀分布，孔 径服从指数分布，并考虑了不同的孔隙率。 2．2孔隙岩石的损伤破坏过程模拟 图1所示为孔隙率为0．130％的岩石细观单元在采用单元生死技术后的损伤破坏过 程。图示为第～主席力分布云图，岩石基体的虑力绝大部分是一致的，只是孔边附近出 现应力集中，但由于损伤开裂，应力值并不会很大。且由于孔隙之间的相互作用，微裂 隙的演化形式比较复杂。当孔隙率较低时，微裂纹只是在孔隙的两端发展。当孔隙率较 高时，可能会出现微裂纹的汇聚贯通，产生新的微裂隙。 Ⅵ．127 m京∞≠鲁第l8％学术年§论i集+计算力学 器■圉豳 幽I扎隙率为0I300／,*J岩6细观■i在小蚓载荷F∞单“±她模拟 2j孔障岩石的损伤演化规律 为了定量研究孔隙岩石损伤破坏过程q一的演化规律，定义损伤变量为： D：生：—量!L (11 ^0(1一≠) 式中A。为细观单元的原始面积，％为细观单元的边长100}Irn，p为孔隙率，Ao为细观单 兀中损伤区域的面积，Am为损伤失效的有 限元单元的面张。图2所示为箨种孔隙牢 的岩石细观单兀在小同载荷步下计算得到 的损伤变量演化曲线。 由图2可以看出．刘于不肯孔隙的岩 石，不会出现损伤面秘积累，而是突然全 部破坏，对r古青孔隙的岩石，随着应变 的增加．扎隙岩石的损伤面#济不断增大， 而且损伤面般的增长速率越来越大。返就 表明随着变形的发胜，孔睬岩石的损伤程 度由慢到快不断加剧。 图2不同孔隙率下岩i细观单元的损伤演化曲线 三、结论 (1)选择合适的村料牟构模掣，借助单元生死技术，可以模拟材料的损伤破坏行 为，其发展过程实际上就是孔隙结构附近较高水甲麻』而导致的裂纹萌生、扩展乃拿贯 通．从而培终导致多孔材料的失效破坏。 (2)岩石细现雌元的力学响麻反映了岩石中孔隙的影响，研究表明，随着孔隙率 增大．损伤现象更加明