CN108009371B - 一种泥石流运动的底床阻力测算方法及其应用 - Google Patents
一种泥石流运动的底床阻力测算方法及其应用 Download PDFInfo
- Publication number
- CN108009371B CN108009371B CN201711334894.8A CN201711334894A CN108009371B CN 108009371 B CN108009371 B CN 108009371B CN 201711334894 A CN201711334894 A CN 201711334894A CN 108009371 B CN108009371 B CN 108009371B
- Authority
- CN
- China
- Prior art keywords
- debris flow
- box
- shear plate
- measured value
- shear
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Civil Engineering (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Measuring Volume Flow (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种泥石流运动的底床阻力测算方法及其应用。所述测算方法首先在水槽中部底床设置组合了两个压力传感器和一个拉力传感器的剪力盒,然后采集剪力盒的内置盒质量、剪力盒放置坡度、进行水槽实验前的拉力传感器和压力传感器测量值,接着进行水槽实验,并采集泥石流经过时两个压力传感器和一个拉力传感器的测量值,最后计算得到泥石流运动的底床阻力。本发明的泥石流运动底床阻力测算方法通过组合两个压力传感器和一个拉力传感器,能够直接测量并计算得到泥石流在运动过程中的底床阻力,为验证泥石流运动阻力模型以及认识泥石流运动、堆积机理提供可靠的数据,且本发明的剪力盒结构简单,能方便地用于泥石流动力学实验。
Description
技术领域
本发明涉及一种泥石流运动的底床阻力测算方法,及其在泥石流动力学实验中的应用。
背景技术
泥石流是山区常见的地质灾害,常常威胁山区人民的生命和财产安全。为了减轻泥石流灾害带来的损失,诸如拦挡坝等防治工程常被用来调节物质和能量,而这些工程措施的设计需要对泥石流的运动有较为深入的认识。
泥石流之所以会形成很大灾害的一个主要原因是这些颗粒物质的运动具有超强流动特性。这种超强运动性表现为很高的运动流速和很远的冲出距离。高速使得其具有很强的冲击力,对沿程的结构物具有极强的冲击破坏力;超远的冲出距离使得危害的范围大大增加。因此,对泥石流运动的阻力特性的研究是非常重要的一个研究课题。由于直接测量泥石流运动阻力存在困难,目前泥石流运动阻力的确定主要是基于一些理论推导,而缺乏针对泥石流运动阻力的直接测量计算,使得对泥石流运动的认识存在很大的局限性。
发明内容
本发明的目的就是针对现有技术的不足,提供一种在泥石流动力学实验中测算泥石流运动的底床阻力的方法,该方法能够通过泥石流水槽实验直接测量并计算得到泥石流运动的底床阻力,为验证泥石流运动阻力模型以及认识泥石流运动、堆积机理提供可靠的数据。
为实现上述目的,本发明的技术方案是:
泥石流运动的底床阻力的确定是研究泥石流运动及动力学特性的关键要素。为了能直接测量得到测算泥石流运动底床阻力所需参数,本发明设计了一个剪力盒。具体而言,本发明提出一种泥石流运动的底床阻力测算方法,具体步骤如下:
(一)设置水槽,将剪力盒设于水槽中部的底床中央位置。所述剪力盒包括一个外置盒和一个设于外置盒内的内置盒;内置盒底面与外置盒底面之间通过动滑轮连接,内置盒在外置盒内能自由滑动;内置盒顶面与外置盒顶面位于同一平面、且通过柔性防水布料相连接,内置盒顶面为剪力板;内置盒内设有支架(用于支撑能测量泥石流总压力的压力传感器A),支架上设有与剪力板无缝相连的压力传感器A;内置盒的上游端通过拉力传感器与外置盒连接,内置盒的下游端通过压力传感器B与外置盒连接。
(二)通过实际测量,确定内置盒的质量m、单位kg;通过现场调查,确定剪力盒放置坡度θ、单位°,由于野外泥石流沟道流通段的自然坡度为0-20°之间,因此水槽坡度为0-20°之间,剪力盒放置坡度θ一般为0°-20°;采集没有泥石流经过剪力板时拉力传感器的测量值f0拉、单位N;采集没有泥石流经过剪力板时压力传感器B的测量值f0压、单位N。
(三)配制相应容重和粘度的泥石流进行水槽实验,采集泥石流经过剪力板时拉力传感器的测量值f拉、单位N,泥石流经过剪力板时压力传感器B的测量值f压、单位N,及泥石流经过剪力板时压力传感器A的测量值σ1、单位N。
(四)通过以下公式确定泥石流运动的底床阻力f泥
式中,f泥—泥石流运动的底床阻力,单位N;
f压—泥石流经过剪力板时压力传感器B的测量值,单位N,由步骤(三)确定;
f拉—泥石流经过剪力板时拉力传感器的测量值,单位N,由步骤(三)确定;
m—内置盒(2)的质量,单位kg,由步骤(二)确定;
g—重力加速度,取值9.8m/s2;
θ—剪力盒放置坡度,单位°,由步骤(二)确定;
σ1—泥石流经过剪力板时压力传感器A的测量值,单位N,由步骤(三)确定;
f0压—没有泥石流经过剪力板时压力传感器B的测量值,单位N,由步骤(二)确定;
f0拉—没有泥石流经过剪力板时拉力传感器的测量值,单位N,由步骤(二)确定。
所述泥石流运动的底床阻力测算方法适用于泥石流动力学实验中测量泥石流运动的底床阻力。本发明的主要技术思想:当泥石流向下运动时,对内置盒做受力分析得
mgsinθ-(mgcosθ+σ1)μ+f泥-f压-f拉=0 公式1
公式1中,m为内置盒的质量,单位kg;g为重力加速度,取g=9.8m/s2;θ为剪力盒放置坡度,单位°;σ1为泥石流经过剪力板时压力传感器A的测量值,单位牛顿;μ为内置盒与外置盒之间的滑动摩擦系数;f泥为泥石流运动的底床阻力,单位牛顿;f压为泥石流经过剪力板时压力传感器B的测量值,单位牛顿;f拉为泥石流经过剪力板时拉力传感器的测量值,单位牛顿。
在测量泥石流底床运动阻力之前,需要用公式2对内置盒底面的滑动摩擦进行标定,如下式所示,
公式2中,μ为内置盒与外置盒之间的滑动摩擦系数;θ为剪力盒放置坡度,单位°;m为内置盒的质量,单位kg;g为重力加速度,取g=9.8m/s2;f0压为没有泥石流经过剪力板时压力传感器B的测量值,单位牛顿;f0拉为没有泥石流经过剪力板时拉力传感器的测量值,单位牛顿。
结合公式1和公式2可以计算得到泥石流底床的运动阻力f泥,如下式所示,
与现有技术相比,本发明的有益效果是:通过组合两个压力传感器和一个拉力传感器,能够直接测量并计算得到泥石流在运动过程中的底床阻力,且本发明的剪力盒结构简单,能方便地用于泥石流动力学实验。
附图说明
图1是实施例中泥石流动力学实验示意图。
图2是剪力盒的纵坡面结构示意图
图中标号如下:
1 外置盒 2 内置盒
3 动滑轮 4 柔性防水布料
5 剪力板 6 支架
7 压力传感器A 8 拉力传感器
9 压力传感器B 10 水槽
11 剪力盒 12 闸门
13 料斗
θ 剪力盒放置坡度
具体实施方式
下面对本发明的优选实施例作进一步的描述。
如图1、图2所示。采用本发明的测算方法在泥石流动力学实验中测算泥石流运动的底床阻力,具体步骤如下:
第一步,搭建泥石流动力学实验所需装置,如图1所示,其中水槽10长6m、宽0.4m、水槽10坡度为10°,将剪力盒11设于水槽10中部的底床、距闸门12位置3.0m处。所述剪力盒11包括一个外置盒1和一个设于外置盒1内的内置盒2;内置盒2底面与外置盒1底面之间通过动滑轮3连接,内置盒2在外置盒1内能自由滑动;内置盒2顶面与外置盒1顶面位于同一平面、且通过柔性防水布料4相连接,内置盒2顶面为面积S=0.01m2的剪力板5;内置盒2内设有支架6,支架6上设有与剪力板5相连的压力传感器A7;内置盒2的上游端通过拉力传感器8与外置盒1连接,内置盒2的下游端通过压力传感器B 9与外置盒1连接。
第二步,通过实际测量,确定内置盒2的质量m为0.1kg;通过现场调查,确定剪力盒11放置坡度θ为10°;采集没有泥石流经过剪力板5时拉力传感器8的测量值f0拉为0.020N;采集没有泥石流经过剪力板5时压力传感器B 9的测量值f0压为0.015N。
第三步,配制实验用泥石流进行水槽实验,泥石流流体密度为1700kg/m3,方量为0.19m3,置于料斗13中,提拉闸门12释放泥石流。采集泥石流经过剪力板5时拉力传感器8的测量值f拉为0.752N,泥石流经过剪力板5时压力传感器B 9的测量值f压为0.478N,及泥石流经过剪力板5时压力传感器A7的测量值σ1为13.328N。
第四步,通过以下公式确定泥石流运动的底床阻力f泥
计算得到泥石流运动的底床阻力f泥为1.360N。
将流体密度为1700kg/m3、方量为0.29m3的泥石流置于料斗13中,再次进行水槽实验,即再次进行上述方法的第三步和第四步;采集泥石流经过剪力板5时拉力传感器8的测量值f拉为0.784N,泥石流经过剪力板5时压力传感器B 9的测量值f压为0.481N,及泥石流经过剪力板5时压力传感器A7的测量值σ1为19.992N;计算得到泥石流运动的底床阻力f泥为1.535N;进一步计算得到泥石流运动底床的剪应力τ=153.5Pa,实测泥石流运动流速为5.2m/s,泥深为0.12m。
假定泥石流体为宾汉流体,即式中τ为泥石流底床剪应力,τb为泥石流体屈服应力,η为泥石流体粘滞系数,u为泥石流流速,h为泥石流流深。将两次实验测量数据代入宾汉模型分别得及联立两个方程得τb=23.5Pa,η=1.0Pa·s。
采集水槽实验所用泥石流样品进行室内流变实验,得到泥石流体屈服应力为22.7Pa,泥石流体粘滞系数为0.96Pa·s,这一结果与上述水槽实验得到的流变参数几乎一致。因此,可以验证本发明的泥石流运动的底床阻力测算方法能够用于泥石流动力学实验中流变参数的确定。
Claims (3)
1.一种泥石流运动的底床阻力测算方法,其特征在于:所述泥石流运动的底床阻力测算方法步骤如下:
(一)将剪力盒(11)设于水槽(10)中部的底床;所述剪力盒(11)包括一个外置盒(1)和一个设于外置盒(1)内的内置盒(2),内置盒(2)底面与外置盒(1)底面之间通过动滑轮(3)连接、内置盒(2)在外置盒(1)内能自由滑动,内置盒(2)顶面与外置盒(1)顶面位于同一平面、且通过柔性防水布料(4)相连接,内置盒(2)顶面为剪力板(5),内置盒(2)内设有支架(6),支架(6)上设有与剪力板(5)相连的压力传感器A(7),内置盒(2)的上游端通过拉力传感器(8)与外置盒(1)连接,内置盒(2)的下游端通过压力传感器B(9)与外置盒(1)连接;
(二)通过实际测量,确定内置盒(2)的质量m、单位kg;通过现场调查,确定剪力盒(11)放置坡度θ、单位°;采集没有泥石流经过剪力板(5)时拉力传感器(8)的测量值f0拉、单位N;采集没有泥石流经过剪力板(5)时压力传感器B(9)的测量值f0压、单位N;
(三)进行水槽实验,采集泥石流经过剪力板(5)时拉力传感器(8)的测量值f拉、单位N,泥石流经过剪力板(5)时压力传感器B(9)的测量值f压、单位N,及泥石流经过剪力板(5)时压力传感器A(7)的测量值σ1、单位N;
(四)通过以下公式确定泥石流运动的底床阻力f泥
式中,f泥—泥石流运动的底床阻力,单位N;
f压—泥石流经过剪力板(5)时压力传感器B(9)的测量值,单位N,由步骤(三)确定;
f拉—泥石流经过剪力板(5)时拉力传感器(8)的测量值,单位N,由步骤(三)确定;
m—内置盒(2)的质量,单位kg,由步骤(二)确定;
g—重力加速度,取值9.8m/s2;
θ—剪力盒(11)放置坡度,单位°,由步骤(二)确定;
σ1—泥石流经过剪力板(5)时压力传感器A(7)的测量值,单位N,由步骤(三)确定;
f0压—没有泥石流经过剪力板(5)时压力传感器B(9)的测量值,单位N,由步骤(二)确定;
f0拉—没有泥石流经过剪力板(5)时拉力传感器(8)的测量值,单位N,由步骤(二)确定。
2.根据权利要求1所述泥石流运动的底床阻力测算方法,其特征在于:剪力盒(11)放置坡度θ为0°-20°。
3.根据权利要求1所述泥石流运动的底床阻力测算方法,其特征在于:适用于泥石流动力学实验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711334894.8A CN108009371B (zh) | 2017-12-14 | 2017-12-14 | 一种泥石流运动的底床阻力测算方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711334894.8A CN108009371B (zh) | 2017-12-14 | 2017-12-14 | 一种泥石流运动的底床阻力测算方法及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108009371A CN108009371A (zh) | 2018-05-08 |
CN108009371B true CN108009371B (zh) | 2020-12-01 |
Family
ID=62059028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711334894.8A Active CN108009371B (zh) | 2017-12-14 | 2017-12-14 | 一种泥石流运动的底床阻力测算方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108009371B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627199B (zh) * | 2018-07-18 | 2020-05-05 | 蔡灏晖 | 一种可同时测量泥石流多项力学及运动学参数的装置 |
CN109211522B (zh) * | 2018-10-17 | 2023-10-17 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流泥舌沿程空间演化测量系统及测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699296A (zh) * | 2009-10-30 | 2010-04-28 | 中国科学院水利部成都山地灾害与环境研究所 | 一种泥石流流速测量系统与测量方法及应用 |
CN102147325A (zh) * | 2011-02-28 | 2011-08-10 | 浙江工业大学 | 非恒定泥石流实验装置 |
CN104952211A (zh) * | 2014-03-31 | 2015-09-30 | 中铁西北科学研究院有限公司深圳南方分院 | 一种泥石流灾害预警方法及无线预警系统 |
CN106529198A (zh) * | 2016-12-13 | 2017-03-22 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流全过程数值模拟及数值计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101564214B1 (ko) * | 2014-08-18 | 2015-10-29 | 한국지질자원연구원 | 간극수압-지표변위-전단면 완전연계 산사태 조기탐지 시험장치 |
-
2017
- 2017-12-14 CN CN201711334894.8A patent/CN108009371B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699296A (zh) * | 2009-10-30 | 2010-04-28 | 中国科学院水利部成都山地灾害与环境研究所 | 一种泥石流流速测量系统与测量方法及应用 |
CN102147325A (zh) * | 2011-02-28 | 2011-08-10 | 浙江工业大学 | 非恒定泥石流实验装置 |
CN104952211A (zh) * | 2014-03-31 | 2015-09-30 | 中铁西北科学研究院有限公司深圳南方分院 | 一种泥石流灾害预警方法及无线预警系统 |
CN106529198A (zh) * | 2016-12-13 | 2017-03-22 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流全过程数值模拟及数值计算方法 |
Non-Patent Citations (3)
Title |
---|
"尾矿库溃坝泥石流相似模拟试验台设计及验证";敬小非等;《中国安全生产科学技术》;20170731;第13卷(第7期);第24-29页 * |
"泥石流动力过程模拟及特征研究";樊赟赟;《中国博士学位论文全文数据库(电子期刊)》;20110815;A011-1 * |
"粘性泥石流阻力和运动方程验证分析";周必凡;《山地学报》;19990228;第17卷(第1期);第55-58页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108009371A (zh) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105510153B (zh) | 一种大型土体界面剪切试验模型及试验方法 | |
Rong et al. | Laboratory investigation of nonlinear flow characteristics in rough fractures during shear process | |
CN108009371B (zh) | 一种泥石流运动的底床阻力测算方法及其应用 | |
WO2016107360A1 (zh) | 一种可调多功能破碎岩体崩滑试验装置及其使用方法 | |
CN202339307U (zh) | 测定泥页岩吸水扩散系数装置 | |
CN103293179B (zh) | 基于悬吊法的混凝土早期热膨胀系数测试装置及测试方法 | |
Hotta | Basal interstitial water pressure in laboratory debris flows over a rigid bed in an open channel | |
CN106932325A (zh) | 动静荷载作用下岩石裂隙泥水流体渗流装置与试验方法 | |
Grabowski | Measuring the shear strength of cohesive sediment in the field | |
CN206990356U (zh) | 一种模拟坝肩软弱夹层的实验装置 | |
CN113466126B (zh) | 一种可进行饱和、固结并可考虑温度效应的多功能界面剪切装置 | |
Sadeghian et al. | Comparison of binomial and power equations in radial non-Darcy flows in coarse porous media | |
KR102471781B1 (ko) | 비점착성 혼합 토사의 침강속도 자동 측정 장치 | |
Van Tonder et al. | Seepage column hydraulic conductivity tests in the geotechnical centrifuge | |
CN104297085B (zh) | 一种数显式模具硬度测量装置 | |
Abbood | An experimental model for flow through porous media using water filter | |
Ter-Martirosyan et al. | Rheological models creation on the results triaxial tests of sands | |
CN204807421U (zh) | 颗粒柱顶加载底部质量测试系统 | |
CN114264792B (zh) | 一种突水溃砂水砂运移模拟试验装置及试验方法 | |
US8677815B2 (en) | Apparatus for a borehole influx fluid parameters measurement | |
CN103499456B (zh) | 防砂筛管过流性能检测装置及方法 | |
Shrivastava et al. | Shear behavior of rock under different normal stiffness | |
Yamazaki et al. | Geotechnical properties of deep-sea sediments and the in-situ measurement techniques | |
Gleason | Particle image velocimetry: A new technique to measure strain in loaded snow | |
Marks et al. | Design of protection structures: the role of the grainsize distribution |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |