CN114002052A - 一种全尾砂充填体抗剪切强度的计算方法 - Google Patents

一种全尾砂充填体抗剪切强度的计算方法 Download PDF

Info

Publication number
CN114002052A
CN114002052A CN202111298345.6A CN202111298345A CN114002052A CN 114002052 A CN114002052 A CN 114002052A CN 202111298345 A CN202111298345 A CN 202111298345A CN 114002052 A CN114002052 A CN 114002052A
Authority
CN
China
Prior art keywords
shear
shear strength
shearing
area
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111298345.6A
Other languages
English (en)
Inventor
孙伟
高通
卢开放
赵建光
刘增
李兆宇
姜明归
樊锴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202111298345.6A priority Critical patent/CN114002052A/zh
Publication of CN114002052A publication Critical patent/CN114002052A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02—Details
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02—Details
    • G01N3/04—Chucks
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种全尾砂充填体抗剪切强度的计算方法，属于矿业工程技术领域。本发明设置应变控制式直剪仪的剪切速率V，确定剪切盒半径r和测力环系数C，制备全尾砂充填体试件，对全尾砂充填体试件进行直剪试验，记录剪切时间T和测力环读数M，计算出剪切位移L；将动态变化的有效剪切面积SE带入抗剪强度τ的表达式中，得出考虑剪切面积动态变化的抗剪强度表达式τE，通过剪切时间T和测力环读数M计算出试件考虑剪切面积动态变化的充填体抗剪强度τE；绘制抗剪强度τE‑法向应力σ的关系曲线，通过摩尔‑库伦准则计算出黏聚力c和内摩擦角

Description

一种全尾砂充填体抗剪切强度的计算方法
技术领域
本发明涉及一种全尾砂充填体抗剪切强度的计算方法，属于矿业工程技术领域。
背景技术
近年来，充填采矿法因具有控制地表沉降、维护采场稳定性、实现废弃尾砂二次利用等优势，在金属矿山绿色开采中得到了广泛的应用。充填采矿其关键技术在于充填体的强度，大量的现场勘察表明，充填体破环失稳多是由于其抗剪强度不足导致的剪切破坏。
目前充填体抗剪强度的获取多采用传统的直剪实验，即通过电动(手动)直剪仪得到试件破坏时的抗剪强度，利用摩尔-库伦准则
Figure BDA0003337261280000011
绘制抗剪强度τ与法向应力σ的关系曲线，进一步计算得到抗剪强度指标黏聚力c及内摩擦角
Figure BDA0003337261280000012
传统直剪试验因未考虑随着剪切的进行，充填体剪切面积动态变化不断变小的问题，致使所获取的抗剪强度及其指标黏聚力c、内摩擦角
Figure BDA0003337261280000013
与真实值存在一定的误差。
发明内容
本发明针对充填体抗剪强度测试时其指标黏聚力c、内摩擦角
Figure BDA0003337261280000014
与真实值存在较大误差的问题，提出了一种全尾砂充填体抗剪切强度的计算方法，即通过考虑直剪实验中剪切面积的动态变化，采用积分法计算剪切面积动态变化的充填体抗剪强度及黏聚力c、内摩擦角
Figure BDA0003337261280000015
值。
一种全尾砂充填体抗剪切强度的计算方法，具体步骤如下：
(1)设置应变控制式直剪仪的剪切速率V，确定剪切盒半径r和测力环系数C；根据《JJG144-2007标准测力仪检定规程》对测力环系数进行标定；
(2)制备全尾砂充填体试件：
1)对全尾砂充填体试件所需材料进行取样，并对尾砂物理性质和化学成分进行分析；
2)根据剪切盒尺寸，制备试件，养护3d后分别在100kPa、200kPa、300kPa、400kPa法向应力下进行直剪实验，为降低实验误差，每组进行3个平行试验，最后对剪切应力选取平均值；
(3)对全尾砂充填体试件进行直剪试验，记录剪切时间T和测力环读数M，计算出剪切位移L；
(4)将动态变化的有效剪切面积SE带入抗剪强度τ的表达式中，得出考虑剪切面积动态变化的抗剪强度表达式τE，通过剪切时间T和测力环读数M计算出试件考虑剪切面积动态变化的充填体抗剪强度τE；
(5)绘制抗剪强度τE-法向应力σ的关系曲线，通过摩尔-库伦准则计算出黏聚力c和内摩擦角
Figure BDA0003337261280000021
值。
所述步骤(1)剪切速率V的设置以使全尾砂充填体试件在5～8min内剪断；
根据直剪仪型号确定剪切盒半径r；
所述步骤(3)剪切位移L为
L＝V·T
式中：V为剪切速率(mm/min)，T为剪切时间(min)。
所述步骤(4)抗剪强度τ的表达式为
Figure BDA0003337261280000022
式中：τ为抗剪强度，C为测力环系数，M为测力环读数，SE为有效剪切面积；
动态变化的有效剪切面积SE通过积分法确定；
如图2所示，O1、O2距离为剪切位移L，剪切盒半径为r，以O1为原点建立直角坐标系，通过积分法确定动态变化的有效剪切面积SE：
Figure BDA0003337261280000023
化简得：
Figure BDA0003337261280000024
式中：SE为动态变化的有效剪切面积，r为剪切盒半径，L为剪切位移；
考虑剪切面积动态变化的充填体抗剪强度τE为、
Figure BDA0003337261280000025
式中：V为剪切速率，T为剪切时间，C为测力环系数，M为测力环读数，SE为有效剪切面积。
所述摩尔-库伦准则为
Figure BDA0003337261280000026
式中：τE为考虑剪切面积动态变化的充填体抗剪强度，σ为法向应力，
Figure BDA0003337261280000027
为内摩擦角，c为黏聚力。
本发明的有益效果是：
(1)本发明针对传统直剪试验中未考虑剪切面积的动态变化对充填体抗剪强度计算结果的影响，通过考虑直剪实验中剪切面积的动态变化，采用积分法计算剪切面积动态变化的充填体抗剪强度及黏聚力c、内摩擦角
Figure BDA0003337261280000031
值，本发明方法与传统抗剪强度的计算方法得出的充填体抗剪强度误差较大，高达8.2～12.36％；
(2)本发明方法具有计算简单，待定参数少，结果可靠等优点，可为全尾砂充填体抗剪强度及其指标c、
Figure BDA0003337261280000032
值的计算提供更准确的结果。
附图说明
图1为直剪实验中剪切面积动态变化示意图；
图2为积分法计算动态变化的剪切面积示意图；
图3为100kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图；
图4为200kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图；
图5为300kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图；
图6为400kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图；
图7为实施例2与传统方法的抗剪强度-法向应力关系对比曲线图。
具体实施方式
下面结合具体实施方式对本发明作进一步详细说明，但本发明的保护范围并不限于所述内容。
实施例1：一种全尾砂充填体抗剪切强度的计算方法，具体步骤如下：
(1)设置应变控制式直剪仪的剪切速率V，确定剪切盒半径r和测力环系数C；剪切速率V的设置以使全尾砂充填体试件在5～8min内剪断；根据直剪仪型号确定剪切盒半径r；根据《JJG144-2007标准测力仪检定规程》对测力环系数进行标定；
(2)制备全尾砂充填体试件：
1)对全尾砂充填体试件所需材料进行取样，并对尾砂物理性质和化学成分进行分析；
2)根据剪切盒尺寸，制备试件，养护3d后分别在100kPa、200kPa、300kPa、400kPa法向应力下进行直剪实验，为降低实验误差，每组进行3个平行试验，对剪切应力选取平均值；
(3)对全尾砂充填体试件进行直剪试验，记录剪切时间T和测力环读数M，计算出剪切位移L：
L＝V·T
式中：V为剪切速率(mm/min)，T为剪切时间(min)；
(4)将动态变化的有效剪切面积SE带入抗剪强度τ的表达式中，得出考虑剪切面积动态变化的抗剪强度表达式τE，通过剪切时间T和测力环读数M计算出试件考虑剪切面积动态变化的充填体抗剪强度τE；
其中抗剪强度τ的表达式为
Figure BDA0003337261280000041
式中：τ为抗剪强度，C为测力环系数，M为测力环读数，SE为有效剪切面积；
动态变化的有效剪切面积SE通过积分法确定；
直剪实验中剪切面积动态变化示意图见图1，积分法计算动态变化的剪切面积示意图如图2所示，O1、O2距离为剪切位移L，剪切盒半径为r，以O1为原点建立直角坐标系，通过积分法确定动态变化的有效剪切面积SE：
Figure BDA0003337261280000042
化简得：
Figure BDA0003337261280000043
式中：SE为动态变化的有效剪切面积，r为剪切盒半径，L为剪切位移；
考虑剪切面积动态变化的充填体抗剪强度τE为、
Figure BDA0003337261280000044
式中：V为剪切速率，T为剪切时间，C为测力环系数，M为测力环读数，SE为有效剪切面积；
(5)绘制抗剪强度τE-法向应力σ的关系曲线，通过摩尔-库伦准则计算出黏聚力c和内摩擦角
Figure BDA0003337261280000045
值；摩尔-库伦准则为
Figure BDA0003337261280000046
式中：τE为考虑剪切面积动态变化的充填体抗剪强度，σ为法向应力，
Figure BDA0003337261280000047
为内摩擦角，c为黏聚力。
实施例2：一种全尾砂充填体抗剪切强度的计算方法，具体步骤如下：
(1)设置应变控制式直剪仪的剪切速率V为0.8mm/min以使全尾砂充填体试件在5～8min内剪断，确定剪切盒半径r为30.9mm，根据《JJG144-2007标准测力仪检定规程》对测力环系数C进行标定，标定结果为C为1.667kPa/0.01mm，经换算取C为0.5001kN/mm，检定结果如表1所示；
表1测力环检定结果(室温：22℃±2°)
Figure BDA0003337261280000051
(2)制备全尾砂充填体试件：
1)对全尾砂充填体试件所需材料进行取样，并对尾砂物理性质和化学成分进行分析；试验尾砂样品取自云南某矿，物理性质及化学成分如表2和表3所示；
表2尾砂物理性质参数
Figure BDA0003337261280000052
表3尾砂化学组成成分(％)
Figure BDA0003337261280000053
2)根据剪切盒尺寸，选取73％质量浓度制备4组试件，其中试件直径R为61.8mm，高为45mm，养护3d后分别在100kPa、200kPa、300kPa、400kPa法向应力下进行直剪实验，为降低实验误差，每组进行3个平行试验，最终对剪切应力选取平均值；
(3)对全尾砂充填体试件进行直剪试验，记录剪切时间T和测力环读数M，计算出剪切位移L：
L＝V·T
式中：V为剪切速率(mm/min)，T为剪切时间(min)；
(4)将动态变化的有效剪切面积SE带入抗剪强度τ的表达式中，得出考虑剪切面积动态变化的抗剪强度表达式τE，通过剪切时间T和测力环读数M计算出试件考虑剪切面积动态变化的充填体抗剪强度τE；
其中抗剪强度τ的表达式为
Figure BDA0003337261280000061
式中：τ为抗剪强度，C为测力环系数，M为测力环读数，SE为有效剪切面积；
动态变化的有效剪切面积SE通过积分法确定；
直剪实验中剪切面积动态变化示意图见图1，积分法计算动态变化的剪切面积示意图如图2所示，O1、O2距离为剪切位移L，剪切盒半径为r，以O1为原点建立直角坐标系，通过积分法确定动态变化的有效剪切面积SE：
Figure BDA0003337261280000062
化简得：
Figure BDA0003337261280000063
式中：SE为动态变化的有效剪切面积，r为剪切盒半径，L为剪切位移；
考虑剪切面积动态变化的充填体抗剪强度τE为、
Figure BDA0003337261280000064
式中：V为剪切速率，T为剪切时间，C为测力环系数，M为测力环读数，SE为有效剪切面积；
100kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图见图3，200kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图见图4，300kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图见图5，400kPa法向应力下本方法与传统方法的剪切应力-剪切位移关系对比曲线图见图6，本方法充分考虑直剪实验过程中有效剪切面积动态变化，不断变小的问题，更符合实际情况，由图3～6可知，不同法向应力下，使用本方法对充填体抗剪强度的计算结果均高于传统方法，且随着法向应力的增大，误差也逐步增大，达到8.2～12.6％；
(5)绘制抗剪强度τE-法向应力σ的关系曲线(见图7)，通过摩尔-库伦准则计算出黏聚力c和内摩擦角
Figure BDA0003337261280000065
值(见表4)；摩尔-库伦准则为
Figure BDA0003337261280000066
式中：τE为考虑剪切面积动态变化的充填体抗剪强度，σ为法向应力，
Figure BDA0003337261280000067
为内摩擦角，c为黏聚力；
表4黏聚力c和内摩擦角
Figure BDA0003337261280000071
值
Figure BDA0003337261280000072
从表1可知，本方法充分考虑直剪实验过程中有效剪切面积动态变化，不断变小的问题，更符合实际情况，与传统抗剪强度的计算方法得出的充填体抗剪强度指标误差较大，其中黏聚力c误差达到4.3％，内摩擦角
Figure BDA0003337261280000073
误差达到12.1％。
以上结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (7)

1.一种全尾砂充填体抗剪切强度的计算方法，其特征在于，具体步骤如下：
(1)设置应变控制式直剪仪的剪切速率V，确定剪切盒半径r和测力环系数C；
(2)制备全尾砂充填体试件；
(3)对全尾砂充填体试件进行直剪试验，记录剪切时间T和测力环读数M，计算出剪切位移L；
(4)将动态变化的有效剪切面积SE带入抗剪强度τ的表达式中，得出考虑剪切面积动态变化的抗剪强度表达式τE，通过剪切时间T和测力环读数M计算出试件考虑剪切面积动态变化的充填体抗剪强度τE；
(5)绘制抗剪强度τE-法向应力σ的关系曲线，通过摩尔-库伦准则计算出黏聚力c和内摩擦角
Figure FDA0003337261270000014
值。
2.根据权利要求1所述全尾砂充填体抗剪切强度的计算方法，其特征在于：步骤(1)剪切速率V的设置以使全尾砂充填体试件在5～8min内剪断。
3.根据权利要求1所述全尾砂充填体抗剪切强度的计算方法，其特征在于：步骤(3)剪切位移L为
L＝V·T
式中：V为剪切速率(mm/min)，T为剪切时间(min)。
4.根据权利要求3所述全尾砂充填体抗剪切强度的计算方法，其特征在于：步骤(4)抗剪强度τ的表达式为
Figure FDA0003337261270000011
式中：τ为抗剪强度，C为测力环系数，M为测力环读数，SE为有效剪切面积。
5.根据权利要求4所述全尾砂充填体抗剪切强度的计算方法，其特征在于：动态变化的有效剪切面积SE通过积分法确定；
Figure FDA0003337261270000012
化简得：
Figure FDA0003337261270000013
式中：SE为动态变化的有效剪切面积，r为剪切盒半径，L为剪切位移。
6.根据权利要求5所述全尾砂充填体抗剪切强度的计算方法，其特征在于：考虑剪切面积动态变化的充填体抗剪强度τE为
Figure FDA0003337261270000021
式中：V为剪切速率，T为剪切时间，C为测力环系数，M为测力环读数，SE为有效剪切面积。
7.根据权利要求6所述全尾砂充填体抗剪切强度的计算方法，其特征在于：摩尔-库伦准则为
Figure FDA0003337261270000022
式中：τE为考虑剪切面积动态变化的充填体抗剪强度，σ为法向应力，
Figure FDA0003337261270000023
为内摩擦角，c为黏聚力。
CN202111298345.6A 2021-11-04 2021-11-04 一种全尾砂充填体抗剪切强度的计算方法 Pending CN114002052A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111298345.6A CN114002052A (zh) 2021-11-04 2021-11-04 一种全尾砂充填体抗剪切强度的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111298345.6A CN114002052A (zh) 2021-11-04 2021-11-04 一种全尾砂充填体抗剪切强度的计算方法

Publications (1)

Publication Number Publication Date
CN114002052A true CN114002052A (zh) 2022-02-01

Family

ID=79927632

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111298345.6A Pending CN114002052A (zh) 2021-11-04 2021-11-04 一种全尾砂充填体抗剪切强度的计算方法

Country Status (1)

Country Link
CN (1) CN114002052A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593460A (zh) * 2018-04-28 2018-09-28 长沙理工大学 基于直剪试验确定土体抗剪强度的动态精确计算方法
CN109163985A (zh) * 2018-08-15 2019-01-08 三峡大学 一种测试岩样抗剪强度和长期蠕变变形的原位测试装置及方法
CN110411861A (zh) * 2019-08-27 2019-11-05 山东农业大学 一种应变控制式直剪试验数据优化方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593460A (zh) * 2018-04-28 2018-09-28 长沙理工大学 基于直剪试验确定土体抗剪强度的动态精确计算方法
CN109163985A (zh) * 2018-08-15 2019-01-08 三峡大学 一种测试岩样抗剪强度和长期蠕变变形的原位测试装置及方法
CN110411861A (zh) * 2019-08-27 2019-11-05 山东农业大学 一种应变控制式直剪试验数据优化方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张亮亮等: "直剪实验剪切面积的修正及误差分析", 《力学与实践》 *

Similar Documents

Publication Publication Date Title
CN110296892B (zh) 基于能量分析的岩石损伤演化过程中特征应力确定方法
CN108287112B (zh) 一种基于三轴压缩试验测定岩石损伤参数的方法
CN103940666A (zh) 一种模拟断续裂隙岩石力学特性的细观参数确定方法
CN110726608B (zh) 一种基于应力-应变曲线能量演化的页岩脆性评价方法
CN107907409A (zh) 一种确定岩石起裂应力的方法、设备及存储设备
CN114002052A (zh) 一种全尾砂充填体抗剪切强度的计算方法
CN111325482A (zh) 一种大跨度地下洞库的围岩快速分级方法
CN109612899A (zh) 一种压力修正式气体渗透率计算方法
CN108593460B (zh) 基于直剪试验确定土体抗剪强度的动态精确计算方法
CN105160165B (zh) 一种应变软化型边坡动态稳定性评价方法
CN109387433A (zh) 基于稳态蠕变速率倒数确定岩石单轴长期强度的方法
CN111398116B (zh) 一种利用特定方向压汞法表征页岩各向异性的方法
CN108225915A (zh) 一种含裂纹类缺陷结构的动态断裂评定方法
CN110348167B (zh) 一种测定金属细节疲劳额定强度截止值DFRcutoff的概率模型方法
CN103984853A (zh) 空气预热器漏风率对锅炉排烟温度的修正方法
CN116106091B (zh) 深部含瓦斯煤体失稳灾变时的顶板弹性能贡献率估算方法
CN111595677A (zh) 软岩长期强度值确定方法以及软岩承载力确定方法
CN114894671B (zh) 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法
CN215066336U (zh) 基于管带机的烧结矿FeO含量在线检测装置
CN111855338B (zh) 一种高相似度相似材料制作方法
CN115795773A (zh) 一种煤体灾变时顶板弹性能贡献率影响因素的分析方法
Wang et al. Numerical analysis on rules of dynamic gas emission during boring process
CN112380484A (zh) 一种测定及计算岩体原位模量的方法
CN107463722B (zh) 一种薄层电阻率测井响应校正的新方法
CN201485417U (zh) 一种用于石油沥青生产的降蜡设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20220201

RJ01 Rejection of invention patent application after publication