岩样变形和破坏过程的二元介质模拟.pdfVIP

水 利 水 运 工 程 学 报 2004年3月 图1 3种基本元件 图2 理想材料的应力应变关系 Fig．1 Three basic components Fig．2 Stress—strain curves of the ideal material 图3 大理岩的应力应变曲线 图4 砂岩的应力应变曲线 Fig．3 Stress·strain curves of marble Fig．4 Stress—strain curves of sandstone 1 破裂强度随机变化的二元介质模型 假设材料试样由大量微元组合而成．每个微元破裂以前均为弹脆性体，破裂以后转化为弹塑性体，且所 有微元破裂前的弹性模量E。和破裂后的弹性模量E 均为常量，破裂前、后的泊松比也相同，只是破裂强度 q以某一规律分布．现只考虑常围压下轴向应力和轴向应变之间的关系．令：平均应力or =or 一or 和平均应 变 = 。一 ( 为加围压or 时的轴向应变)．若用or 和or ， 以分别代表未破裂和已破裂微元的偏应 力和偏应变，则有： ord=(1一A)ordl+Aor (1) d=(1一A) dl+A (2) 其中：A为已破损微元所占的体积率(体积破损率)．令 ，=CoC (c为局部应变系数)．由(2)式可得 s =s ———■一 OCd (3)j 再分别代入ord。=E。8d。和or =E 后，由(1)式可得平均应力 ord=(1一A)cEl d+(1一c+cA)E2 d (4) 现研究破损率A的变化规律．设q 为材料基准破裂强度，即or =q 时，微元发生破裂的比例为 (1一 ) 0·632，并假定破损率A服从下列两种规律： (1)Weibull规律 A=…xp[一( ] (5) 当ordl=0时，A=0；ordI：q 时，A=0．632；ordl=∞时，A=1． (2)线性规律 A=·一 [(·+ )一 ] (6) 第1期 沈珠江，等：岩样变形和破坏过程的二元介质模拟 3 当odl≤[(1+1／m)一2．718／m]q 时，A：0；o-dl=q 时，A=0．632； dl≥(1+1／m)q 时，A=1．(5)式中 Weibull曲线的最大斜率为m(见图5)，恰为(6)式直线的斜率．由于微元的局部应力引起微元破裂，所以 (5)和(6)式中应当用 。而不用平均应力 ．将 。：cE。 代人(5)和(6)式后，可得 A= 一pf-(c ] (7) A= 一 [( + 1)一c ] (8) 其中： =q,／E。；c应当小于1，且应满足初始条件A：0时，c=1．假设 c：1一mlA (9) 式中的系数m。也应小于或等于1．将(9)式代A(7)和(8)式，可得到破损率A与平均应变 之间的关系．再将 其代A(4)式，即可得出平均应力与平均应变之间的关系．但是(4)式中的 该如何取