·1432· 北京科技大学学报 第36卷 12000a 5000[间 4500 10000 4000 一剪切k居=0.5 3500 一拉伸法，=0.5 8000 =一一剪切kk=1D 一一=拉伸k=1.0 ==剪切kk=2.0 3000 =一拉伸：k传=2.0 …剪切k/k=10.0 ·拉伸kk-10.0 6000 2500 2000 4000 1500 1000 2000 500 0051.0 1.52.02.53.03540455.0 0005101520253035404530 轴向应变/% 轴向应变% 图8颗粒接触刚度变化时平行黏结破坏数一轴向应变曲线.(）剪切破坏：(b)拉伸破坏 Fig.8 Broken-bond number-axial strain curve with the change in particle stiffness ratio:(a)shear failure:(b)tensile failure (a) d 图9颗粒接触刚度变化时剪切过程中Y-Z平面内平行黏结变化和接触网络.(a)￡，=1%：(b)￡，=2%：(c)E。=3.5%:(d)￡，= 4.0%;(e)e,=4.5%:(0e,=5.0% Fig.9 Parallel bond evolution and contact networks in the y-Z plane during shearing with the change in particle stiffness ratio:(a)s=1%:(b) 6.=2%:(c)6,=3.5%;(d)E.=4.0%:(e)6a=4.5%:(f06.=5.0% 抗会增大整个砂岩模型的峰前刚度Ea:当 的时间越早 Kp/K。=0.25、0.5、1和2时，计算得到的初始刚 图11描述了加载过程中体应变的变化规律 度Emia=1.598、1.855、2.37和3.40GPa.可见，初 可以看出模型均会出现剪缩的现象，随着Kp,K。 始刚度的值是随着Kp,b/K。增大而增大，与图10的 的增大，剪缩最大值对应的应变越小.说明平行黏 结果相吻合.其结果说明颗粒间黏结物质的刚度越 结的刚度越大，外载作用下的脆性性质越强，越容易 大，则砂岩脆性越强，在相同外荷载作用下出现破坏 出现破坏，与图10的应变变化一致.北 京 科 技 大 学 学 报 第 36 卷 图 8 颗粒接触刚度变化时平行黏结破坏数--轴向应变曲线． ( a) 剪切破坏; ( b) 拉伸破坏 Fig． 8 Broken-bond number--axial strain curve with the change in particle stiffness ratio: ( a) shear failure; ( b) tensile failure 图 9 颗粒接触刚度变化时剪切过程中 Y--Z 平面内平行黏结变化和接触网络． ( a) εa = 1% ; ( b) εa = 2% ; ( c) εa = 3. 5% ; ( d) εa = 4. 0% ; ( e) εa = 4. 5% ; ( f) εa = 5. 0% Fig． 9 Parallel bond evolution and contact networks in the Y--Z plane during shearing with the change in particle stiffness ratio: ( a) εa = 1% ; ( b) εa = 2% ; ( c) εa = 3. 5% ; ( d) εa = 4. 0% ; ( e) εa = 4. 5% ; ( f) εa = 5. 0% 抗会增大整个砂岩模型的峰前刚度 Einitiall ． 当 Kep，pb /Kc = 0. 25、0. 5、1 和 2 时，计算得到的初始刚 度 Einitiall = 1. 598、1. 855、2. 37 和 3. 40 GPa． 可见，初 始刚度的值是随着 Kep，pb /Kc 增大而增大，与图10 的 结果相吻合． 其结果说明颗粒间黏结物质的刚度越 大，则砂岩脆性越强，在相同外荷载作用下出现破坏 的时间越早． 图 11 描述了加载过程中体应变的变化规律． 可以看出模型均会出现剪缩的现象，随着 Kep，pb /Kc 的增大，剪缩最大值对应的应变越小． 说明平行黏 结的刚度越大，外载作用下的脆性性质越强，越容易 出现破坏，与图 10 的应变变化一致． · 2341 ·