AG  >> Vol. 8 No. 4 (August 2018)

    层面接触刚度对叠合岩梁应力分布的影响规律
    Influence Law of Level Contact Stiffness on Stress Distribution of Laminated Rock Beam

  • 全文下载: PDF(1747KB) HTML   XML   PP.756-763   DOI: 10.12677/AG.2018.84081  
  • 下载量: 114  浏览量: 165   国家自然科学基金支持

作者:  

曹力,鲁海峰,姚多喜,徐泽栋,周乐:安徽理工大学地球与环境学院,安徽 淮南;煤炭工业协会矿井水害防治工程研究中心,安徽 淮南

关键词:
叠合岩梁剪切刚度法向刚度应力分布Laminated Rock Beam Shear Stiffness Normal Stiffness Stress Distribution

摘要:

为探讨层面接触刚度参数对叠合梁应力分布的影响规律,采用FLAC3D建立不同层间接触刚度参数下的双层叠合梁数值计算模型,分析了在均布荷载条件下剪切刚度Ks和法向刚度Kn对两端固定叠合岩梁应力分布的影响。研究结果表明:双层叠合固定岩梁在均布荷载作用下,应力分布与层面接触刚度关系明显。当剪切刚度Ks增大,岩梁水平方向最大拉应力和最大压应力相应减小;增大法向刚度Kn,水平方向拉应力逐渐增加,而压应力逐渐减小。同时当Kn和Ks增大到一定值时,岩梁应力分布受接触刚度改变而变化微小,趋于稳定。通过本文的研究,可对数值计算时合理选取接触刚度参数提供了一定的借鉴。

In order to investigate the influence of contact stiffness parameters on the stress distribution of laminated beams, a numerical model of double deck laminated beams under different contact stiffness parameters is established by FLAC3D, and the influence of shear stiffness Ks and normal stiffness Kn on the distribution of stress distribution on both ends of the laminated beam under uniform load conditions is analyzed. The results show that the stress distribution of the double deck fixed rock beam under uniform load is obviously related to the contact stiffness of the layer. When the shear stiffness Ks increases, the maximum tensile stress and the maximum compressive stress in the horizontal direction of the rock beam decrease, and the tensile stress in the horizontal direction increases gradually with the normal stiffness Kn, and the compressive stress gradually decreases. At the same time, when Kn and Ks increase to a certain value, the stress distribution of rock beam changes slightly with the change of contact stiffness, and tends to be stable. Through the research in this paper, it can provide some reference for selecting reasonable contact stiffness parameters in numerical calculation.

1. 引言

层状岩体在自然界中分布十分广泛,由于其存在着一组优势贯穿结构面,与完整岩石相比,其结构面上的抗剪强度较低,易发生倾倒、弯折以及离层等变形破坏 [1] [2] [3] [4] [5] 。求解上述层状岩体力学问题时,一般可将岩体简化为岩梁进行分析计算。目前,大部分学者针对均质梁分析较多,而对于非均质梁的研究主要集中在层间无滑动的组合梁方面。梁运培等 [6] 提出了采场上覆岩层移动的组合模型,对组合岩梁模型在岩层移动规律研究及地表沉降控制方面进行了深入分析。黄菊华等 [7] 讨论了不同叠合方式的叠合梁的应力分析问题,得出不同材料、不同叠合方式对应力的影响规律。徐铎等 [8] 对岩梁材料拉压弹性模量不尽相同的横向叠加组合岩梁分层截面应力大小分布情况作了分析。而分析层间有滑动的叠合梁方面,运用数值计算的方法研究时一般需要到层间接触刚度参数,即剪切刚度和法向刚度 [9] [10] 。而目前在鲜有文献讨论接触刚度参数对叠合梁应力分布的影响规律,鉴于此,本文以两端固定叠合岩梁为例,运用FLAC3D数值计算软件,建立叠合岩梁在均布外荷载作用下的不同接触刚度参数时的应力求解模型。研究成果可为争取选取接触刚度参数提供一定的借鉴与参考。

2. 层面的变形性质

在不同地质作用中,岩体中形成了不同的结构面,比如层理、岩层面、不整合面和软弱夹层。在荷载作用下,岩体中的应力分布受结构面及其性质的影响 [11] 。在此主要研究岩层之间的接触面即层面的变形性质对应力分布的影响。

2.2. 法向变形性质

当结构面法向力 σ n 为压应力时,结构面产生法向压缩变形。开始先为点或线接触,经挤压,局部破碎或劈裂,接触面增加。结构面压缩量呈指数曲线特征 [12] [13] ,其指数函数为:

u = u 0 ( 1 e σ n K n ) (1)

式(1)中u为结构面法向压缩量,cm;u0为法向压应力,Mpa;Kn为结构面法向刚度,Mpa/cm。

法向刚度Kn是反映结构面法向变形性质的重要参数。其定义为在法向应力作用下,结构面产生单位法向变形所需要的应力,数值上等于法向变形曲线上的一点的切线斜率,即:

K n = σ n Δ V j (2)

2.2. 剪切变形性质

在一定的法向压应力作用下,结构面在剪应力作用下产生沿层面方向的滑移形成层间剪切带。结构面在剪力作用下产生塑型变形和脆性破坏。可用应力—变形曲线关系表征其变形规律,用剪切刚度Ks表示其变形性质,剪切刚度Ks是反映结构面剪切变形性质的重要参数 [11] 。剪切刚度Ks是剪切变形曲线斜率:

K s = τ u s (3)

式(3)中Ks为结构面剪切刚度,MPa/cm。对于结构面的抗剪强度,可以用库伦—莫尔准则来表达:

τ n = σ v t g ϕ + C (4)

式(4)中 τ n 为平行结构面的剪应力,MPa; σ n 为垂直结构面的正应力; ϕ 为结构面内摩擦角,(˚);C为结构面的凝聚力,MPa。

3. 数值模型的建立

3.1. 模型概述

为寻找岩层层面接触刚度,即层面的剪切刚度和法向刚度对叠合岩层应力分布和位移的影响,采用弹性力学模型 [14] ,通过FLAC3D进行数值模拟。建立如图1所示模型,双层叠合固定梁由两层岩性不同的岩层上下自然叠合而成,叠合梁的两端为固定端,梁的跨度为l,上下层的高度分别为h1和h2,宽度b为1;两岩层的弹性模量为E1和E2,泊松比为v1,v2;梁表面受均布荷载q作用,上下岩层之间存在接触层面,接触面剪切刚度Ks/Mpa∙cm−1,法向刚度Kn/Mpa∙cm−1

3.2. 参数的选取

计算模型中上层岩梁的弹性模量E1,泊松比v1;下层岩梁的弹性模量E2,泊松比v2。叠合梁的跨度l,梁的上表面承受的均布荷载q,上层高度h1,下层高度h2,梁的宽度b,参数取值见表1。在上、下层之间设置接触面,接触面初始剪切刚度Ks = 5 Mpa/cm,初始法向刚度Kn = 10 Mpa/cm。

当保持接触层面法向刚度Kn初始值不变时,改变层面剪切刚度Ks (Ks = 5 + 5i, I = 0~9),当保持接触层面剪切刚度Ks初始值不变时,改变层面法向刚度Kn(Kn = 10 + 10i, I = 0~9)。

3.3. 监测点的布置

Flac3d模拟中把模型分为240个单元体,监测单元体应力和位移,选取A、B两个单元体作为监测点,监测其水平方向应力和垂直方向应力,监测A和B单元体垂向位移。A单元体坐标为(x, y, z),(x = 15, 16; y = 0, 1; z = 2, 3),B单元体坐标为(x, y, z),(x = 15, 16; y = 0, 1; z = 6, 7),A、B单元体在叠合岩梁模型中位置如图2

Figure 1. Analysis model of laminated rock beam

图1. 叠合岩梁分析模型

Table 1. Parameter value table for analysis model of laminated rock beam

表1. 叠合岩梁分析模型参数取值表

Figure 2. Location of monitoring points A and B

图2. 监测点A、B位置

4. 模拟结果分析

4.1. 应力分析

通过FLAC3D对初始条件下双层叠合岩梁的应力模拟,在接触面剪切刚度Ks = 5 Mpa/cm,法向刚度Kn = 10 Mpa/cm时,岩梁水平方向和垂直方向应力分布如图3(a)、图3(b)所示。在上覆均匀荷载q作用下,简支岩梁上、下层水平方向应力最大值集中在梁的两端,上下层的上缘和下缘分别处于受压和受拉状态,垂直方向应力由上到下分布,均匀减小,上层应力通过接触面传递到下层,通过接触面时由于接触面变形性质作用,传递到下层的应力减小。

1) 改变剪切刚度

当改变接触面剪切刚度Ks时,叠合岩梁应力分布发生变化,随着剪切刚度Ks的增大,岩梁上、下层两端应力向跨中集中,岩梁水平方向最大拉应力和最大压应力都减小,岩梁上、下层最大应力与剪切刚度Ks呈负相关;当剪切刚度Ks小于15 Mpa/cm时,岩梁垂直方向的最大应力随着Ks的增大而减小,当Ks大于15 Mpa/cm时,岩梁垂直方向的最大应力不变,如图4(a)、图4(b)所示。

叠合岩梁上层B单元体应力大于下层A单元体,岩梁上层跨中应力较集中,随着层面剪切刚度Ks的变化,单元体A和B的水平方向和垂直方向应力发生变化,如图5(a)、图5(b)所示。随着Ks的增大,单元体A的水平和垂直方向应力减小,当Ks大于15 Mpa/cm时,单元体A水平应力方向发生变化并随着Ks的继续增大而增大,垂直方向应力减小趋势变缓,当Ks大于30 Mpa/cm时垂直应力不变;单元体B的水平方向应力随着Ks的增大而增大,增速逐渐变缓,垂直方向应力增加,当Ks大于15 Mpa/cm时不变。

2) 改变法向刚度

当叠合岩梁层面剪切刚度Ks不变,改变法向刚度Kn时,叠合岩梁上、下层水平和垂直应力分布发

生变化,随着法向刚度Kn增大,叠合岩梁水平方向最大拉应力逐渐增加并趋于稳定达到14.15 Mpa,最大压应力逐渐减小到17.65 Mpa (图6(a)),其垂直方向最大拉应力和压应力均减小(图6(b)),当Kn大于30 Mpa/cm时,减速变缓。

随着层面法向刚度Kn的增大,单元体A和B的水平和垂直方向应力发生轻微变化,如图7(a)、图7(b)

(a)(b)

Figure 3. Stress distribution of laminated rock beams

图3. 叠合岩梁应力分布

(a) (b)

Figure 4. Maximum stress change of laminated rock beam

图4. 叠合岩梁最大应力变化

(a) (b)

Figure 5. Stress change of cell body A and B

图5. 单元体A、B应力变化

(a) (b)

Figure 6. Maximum stress change of laminated rock beam

图6. 叠合岩梁最大应力变化

(a) (b)

Figure 7. Stress change of cell body A and B

图7. 单元体A、B应力变化

所示,随着Kn的增大,当Kn小于30 Mpa/cm时,单元体A的水平应力增大,垂直应力减小,当Kn大于30 Mpa/cm时,单元体A应力基本不变。随着Kn的增大,单元体B的应力基本不变。

4.2. 位移分析

由于叠合岩梁两端为固定端,所以岩梁两端没有垂直位移,岩梁垂直位移由跨中向两端均匀减小。在接触面初始剪切刚度Ks = 5 Mpa/cm,初始法向刚度Kn = 10 Mpa/cm的条件下,岩梁跨中最大垂直位移14.743 mm,如图8。由于上下层接触面接触刚度的作用,距垮中等距时上层位移大于下层位移。

Figure 8. Vertical displacement of laminated rock beams

图8. 叠合岩梁垂直方向位移

(a) (b)

Figure 9. Vertical displacement of element A and B

图9. 单元体A、B垂直位移

A单元位于叠合岩梁下层,B单元位于岩梁上层,随着层面接触刚度的变化,在上覆均匀荷载作用下,岩梁垂直方向位移有较大变化,其水平方向位移忽略不计。

当保持法向刚度Kn不变,改变剪切刚度Ks,岩梁垂直位移变化如图9(a)所示。当剪切刚度Ks增大时,A单元和B单元的垂直位移都减小,A单元位移小于B单元。改变切向刚度Ks时,叠合岩梁上下层弯曲均减小。

当保持剪切刚度Ks不变,改变法向刚度Kn,岩梁垂直位移变化如图9(b)所示。当剪切刚度Kn增大时,A单元垂直位移增大,B单元垂直位移减小,A单元位移小于B单元。当Kn大于30 Mpa/cm时,A和B垂直位移增速和减速变缓。改变法向刚度Kn时,叠合岩梁上层弯曲减小,下层弯曲增加。

5. 结论

a) 双层叠合简支岩梁在上覆均匀荷载作用下,内部应力分布受中间层面接触刚度的影响,上层应力大于下层应力,上层应力通过接触面时减弱。层面变形性质包括层面的剪切刚度Ks和法向刚度Kn,二者的变化对于双层叠合岩梁内部应力分布有着很大的影响。

b) 当增加接触面剪切刚度Ks时,岩梁水平方向最大拉应力和最大压应力都减小,岩梁上、下层最大应力与剪切刚度Ks呈负相关;当剪切刚度Ks小于15 Mpa/cm时,岩梁垂直方向的最大应力随着Ks的增大而减小,反之则不变。随着剪切刚度Ks的增大,岩梁下层内部的水平和垂直方向应力减小,上层内部的应力增大。当Ks大于15 Mpa/cm时,岩梁内部应力变化趋势变缓直至趋于稳定值。

c) 当增加接触面法向刚度Kn时,叠合岩梁水平方向最大拉应力逐渐增加并趋于稳定达到14.15 Mpa,最大压应力逐渐减小到17.65 Mpa;其垂直方向最大拉应力和压应力均减小,当Kn大于30 Mpa/cm时,减速变缓。岩梁内部上、下层的水平和垂直方向应力发生轻微变化,当Kn小于30 Mpa/cm时,岩梁下层的水平应力增大,垂直应力减小,反之则基本不变。岩梁上层内部应力基本不变。

d) 岩梁垂直位移由跨中向两端均匀减小。由于上、下层接触面接触刚度的作用,距垮中等距时上层位移大于下层位移。当增加剪切刚度Ks时,岩梁上下层垂直位移都减小,下层位移小于上层。当增加法向刚度Kn,岩梁下层垂直位移增大,上层垂直位移减小,A单元位移小于B单元。当Kn大于30 Mpa/cm时,下层和上层垂直位移增速和减速变缓并趋于稳定。

基金项目

国家自然科学基金面上项目(51471008)。

文章引用:
曹力, 鲁海峰, 姚多喜, 徐泽栋, 周乐. 层面接触刚度对叠合岩梁应力分布的影响规律[J]. 地球科学前沿, 2018, 8(4): 756-763. https://doi.org/10.12677/AG.2018.84081

参考文献

[1] 贾蓬, 唐春安, 王述红. 巷道层状岩层顶板破坏机理[J]. 煤炭学报, 2006, 31(1): 11-15.
[2] 李德海. 近水平层状岩层移动规律的探究[J]. 采矿与安全工程学报, 1996(2): 39-42.
[3] 卜庆为, 高阳, 常瑞. 水平应力作用下层状顶板回采巷道围岩失稳机理及支护分析[J]. 煤炭技术, 2017, 36(11): 12-14.
[4] 王洪建, 刘大安, 黄志全, 等. 层状页岩岩石力学特性及其脆性评价[J]. 工程地质学报, 2017, 25(6): 1414-1423.
[5] 徐嘉谟, 马凤山, 曹家源, 等. 岩体在工程状态下间断面的力学效应[J]. 工程地质学报, 2018, 26(1): 10-23.
[6] 梁运培, 孙东玲. 岩层移动的组合岩梁理论及其应用研究[J]. 岩石力学与工程学报, 2002(5): 654-657.
[7] 黄菊华, 钱应平, 李厚民. 不同形式叠合梁的应力分析[J]. 湖北工业大学学报, 2008(3): 87-89.
[8] 徐铎, 赵建峰, 田慧强, 等. 横向叠加岩梁截面应力分析[J]. 煤炭工程, 2011(10): 113-116.
[9] 孟召平, 陆鹏庆, 贺小黑. 沉积结构面及其对岩体力学性质的影响[J]. 煤田地质与勘探, 2009, 37(1): 33-37.
[10] 孙宗颀. 不连续面应力–变形性质的研究[J]. 岩石力学与工程学报, 1987(4): 287-300.
[11] 刘佑荣, 唐辉明. 岩体力学[M]. 武汉: 中国地质大学出版社, 1999.
[12] 孙广忠. 岩体结构力学[M]. 北京: 科学出版社, 1988.
[13] 潭学术. 复合岩体力学理论及其应用[M]. 北京: 煤炭工业出版社, 1994.
[14] 苗丹, 刘一华. 双材料叠合简支梁受均布载荷作用时的解析解[J]. 应用力学学报, 2016, 33(4): 589-595, 735.