考虑岩石应变率效应的DDA程序.pdf

２０１５年１０月　 　 　 　 　 　 　 　 　 　 　 考虑岩石应变率效应的ＤＤＡ程序　 余德运ꎬ等　 　 　 　 　 　 　 　 　 　 　 ４５ ｄｏｉ:１０.３９６９/ｊ.ｉｓｓｎ.１００１￣８３５２.２０１５.０５.０１１ 考虑岩石应变率效应的ＤＤＡ程序❋ ① ① ② 余德运 　 刘殿书 　 何成龙 ①中国矿业大学(北京)力学与建筑工程学院(北京ꎬ１０００８３) ②北京理工大学爆炸科学与技术国家重点实验室(北京ꎬ１０００８１) [摘　 要]　 针对原经典ＤＤＡ程序在研究岩石动力学问题时不考虑岩石的应变率效应ꎬ在ＤＤＡ程序中引入与岩石 应变率相关的动态强度表达式ꎬ将块体间接触弹簧破坏准则发展为动态破坏准则ꎬ进而完善ＤＤＡ程序对岩石动力 学破坏问题的处理能力ꎮ 用原ＤＤＡ程序和改进ＤＤＡ程序对动载作用下岩石单轴拉、压破坏实验进行模拟ꎬ两者的 模拟结果对比分析表明:在不同的加载速率作用下ꎬ原 ＤＤＡ程序得到的失效应力和破坏形式基本没有变化ꎬ这与 实验结果不相符ꎮ 但是ꎬ在不同的加载速率作用下ꎬ改进的ＤＤＡ程序得到破坏形式能与实验结果相符ꎬ且能体现 动载作用下岩石动态强度随加载速率的提高而增大的强度特性ꎮ [关键词]　 非连续变形分析ꎻ应变率ꎻ数值模拟 [分类号]　 ＴＤ２３５ꎻＯ３４２ 引言 １　 原ＤＤＡ 程序中的岩石断裂处理方法及存在问 ＤＤＡ(ｄｉｓｃｏｎｔｉｎｕｏｕｓ ｄｅｆｏｒｍａｔｉｏｎ ａｎａｌｙｓｉｓ)方法是 题 [１] １９８４年由石根华 提出的一种对诸如岩石等不连 １.１　 原ＤＤＡ 程序中的岩石断裂处理方法 续块体系统的力学响应问题进行数值计算的方法ꎮ ＤＤＡ程序研究岩石力学问题时ꎬ将岩石节理分 该方法以块体的位移为未知数ꎬ基于最小势能原理 为真实节理和虚拟节理ꎬ如图１所示ꎮ 建立平衡方程ꎬ计算效果高ꎬ并且能充分考虑研究对 象的不连续特性ꎬ可模拟不连续块体的非连续变形 和大位移运动情况ꎬ已成为研究岩石这种典型的不 连续介质的力学响应问题的重要方法ꎮ 岩石在动载与准静载条件下的力学响应的区别 在于:动载荷下的力学问题需要考虑两种动态效应ꎬ [２] 即惯性效应和与材料本构相关的应变率效应 ꎮ 　 　 原经典ＤＤＡ程序研究岩石动态力学问题时ꎬ岩石的 图１　 ＤＤＡ程序中的岩石节理划分 强度固定不变ꎬ仅通过添加动态系数来考虑惯性效 Ｆｉｇ.１　 Ｄｉｖｉｓｉｏｎ ｏｆ ｒｏｃｋｊｏｉｎｔ ｉｎ ＤＤＡ ｐｒｏｇｒａｍ [３] 应 ꎬ不考虑应变率效应ꎬ因此精确度不高ꎮ 一直 　 　 直实节理即为岩体中客观存在的节理ꎬ并假设 以来ꎬ从事 ＤＤＡ 程序学习和研究的专家、学者用 其节理面的抗拉强度为零ꎮ 虚拟节理是为了划分块 ＤＤＡ程序模拟研究岩石动力学问题时ꎬ都没有考虑 体单元而虚构的节理ꎬ其客观上是不存在的ꎬ被虚拟 [４￣９] 岩石在动载条件下的应变率效应 ꎮ 节理分开的块体实际上是连续的ꎬ并假设虚拟节理 为了完善ＤＤＡ程序对岩石动力学破坏问题的 面强度参数等于岩石的强度参数ꎬ从而实现块体之 处理能力ꎬ以期能够将ＤＤＡ程序更好地应用到岩石 爆破等动态破坏问题的数值模拟当中ꎬ本文在原经 间的黏结ꎬ反映其连续特性ꎮ 如果接触弹簧力超过 典ＤＤＡ程序的基