岩石材料剪切破坏问题的剪切应变率梯度解析解.docxVIP

岩石材料剪切破坏问题的剪切应变率梯度解析解 1 非局部塑性剪应变率 包括文献[1、2、3、4、5、6、7、8、9、10、11、12、13、14、15、16和17]。文献求解了直杆单向拉伸问题，通过“非率型方式”求得了塑性应变，之后，对塑性应变进行求导，得到塑性应变率。文献在研究岩石力学问题时也采用这一方式，即并非完全采用“率型方式”。 对于同一问题而言，笔者认为，完全采用“率型方式”可以使得理论分析一脉相承。和“非率型方式”相比，除了侧重于研究（塑性）应变率、速度及应力率以外，采用“率型方式”研究同一问题在理论分析方面可能要简洁一些。 大多数关于梯度塑性理论的文献都普遍采用了“非率型方式”，限于篇幅无法列出全部文献。鲜见的文献直接采用“率型方式”建立梯度塑性理论。受此启发，笔者尝试根据非局部模型导出了一维二阶剪切应变率梯度格式，并对单轴压缩岩石试样剪切破坏的轴向响应进行了理论分析。 2 非局部塑性剪应变率 连续介质应变梯度模型可由连续介质非局部模型通过泰勒级数展开而获得。由于微结构之间存在相互作用，且其影响范围广，因此，材料上某一点的应力不仅取决于该点的状态，而且，还取决于与该点临近的诸多点的状态。根据Kroner,Eringen和Edelen提出的早期非局部模型：非局部塑性剪应变率可以表达为其本身所对应的参数的加权平均，该计算在y周围长度为0L的邻域内进行： 式中ξ为从点y到其邻域内任一点的距离；g(ξ)代表某一权函数，它随ξ至y距离的增加而单调下降，微结构相互作用的广泛程度由权函数中所包含的内部长度参数决定。将γp(y+ξ)在y点采用泰勒级数展开，并代入式（1），考虑到各向同性假设，可以将各奇次项消去，仅保留前二项，可以得到： 取非局部权函数为 式中l称之为材料内部长度参数。 将式（3）代入式（2），经过积分运算可得 可见，非局部塑性剪应变率，不仅和局部塑性剪应变率有关，而且，还和其二阶应变率梯度有关。将上式中替代传统本构模型的，就可在传统塑性理论框架之内考虑应变梯度率效应。 3 剪应力的本构模型 当剪切应力超过强度后，假设两端受压应力作用的岩样在内部产生与压应力σ之间的夹角为α，宽度为w的稳定平直剪切带。试件高度为L。假设剪切带内部的材料被剪切破坏，并沿剪切带反对称线发生相对滑动；带外为弹性体。剪切滑动受控于与剪切带方向平行的剪应力τ。假设剪切带边缘所受到的外部弹性体的剪应力τ为均布的。 假设准脆性材料剪应力与剪应变的本构关系为双线性，见图1。剪切弹性模量和剪切降模量分别为G和λ，λ0；τc为剪应力峰值强度。剪切带内部的岩体在剪应力τ（垂直于OY轴）的作用下，将发生剪切屈服，将其简化为沿剪切带法向的一维问题。根据图1，剪切应力率与塑性剪切应变率的关系为 式中。将式（5）中用式（4）中替代，可以解出式（5）非齐次常系数微分方程的通解，考虑到是坐标y的偶函数，故需将解中的非对称项消去，再利用弹性区与剪切塑性区交界处的边界条件：在y=±w/2时，，可解得 仅关心塑性剪切应变率达到最大值时的剪切带宽度，即，可以得到剪切带宽度w的最小非零解及塑性剪切应变率分别为 设剪切带上下盘相对错动速度为u，则有 4 iii类变形行为 在应变软化阶段，岩石试样端部的位移量主要是由岩样整体的压缩变形量和剪切带外部岩体沿剪切带发生滑移而引起的，岩石试件的轴向劈裂不会引起较大的岩石试样端部的位移量。因此，岩石试样端部的速度等于将岩样视为一体的压缩速度和剪切带外部岩体沿剪切带发生滑移而引起的岩石试样端部的速度之和。根据线弹性虎克定律，可表示为 根据式（10）和及c=Gλ(G+λ)，可以表示为 由于弹性模量与剪切模量有如下关系： 因而，岩石试样端部的速度为 式中v为泊松比。 在应变软化阶段，有 因此，当K0，则岩石试样端部的速度同号，均为负。为负意味着岩样发生突然回跳，发生II类变形行为。文献式（12）与K0等价，其推导过程在此给出。当K0，则岩石试样端部的速度异号。为正，发生I类变形行为。 为研究试件高度及泊松比对岩样轴向响应的影响，参数取值如下：l=0.002 m,G=2 GPa，λ=.02 GPa，α=30o。计算结果见图2。试样高度越大，岩石试件越容易发生脆性破坏；反之，岩石试件越容易发生韧性破坏。泊松比越小，岩石试件越容易发生脆性破坏；反之，亦然。 为研究剪切降模量及岩石材料内部长度对岩样轴向响应的影响，参数取值如下：L=0.15 m,G=2 GPa,v=.025，α=30o。计算结果见图3。岩石材料内部长度越小，岩石试件越容易发生脆性破坏；反之，亦然。岩石材料剪切降模量越大，岩石试件越容易发生脆性破坏。反之，亦然。 为研究剪切带倾角及剪切带宽度对岩样轴向响应的影响，参数取值如下：λ=0.2 GPa,L=0.15 m,G=2 GPa,v=.0