[发明专利]一种联合充填采矿中充填体的充填高度分析方法有效
申请号: | 202111316462.0 | 申请日: | 2021-11-08 |
公开(公告)号: | CN113969802B | 公开(公告)日: | 2023-09-08 |
发明(设计)人: | 刘浪;侯东壮;张波;王美;郇超;朱梦博;阮仕山;方治余;杨潘;屈慧升;孙伟吉 | 申请(专利权)人: | 西安科技大学;西安弗尔绿创矿业科技有限责任公司 |
主分类号: | E21F15/00 | 分类号: | E21F15/00;G01N25/16;G06F17/11 |
代理公司: | 西安启诚专利知识产权代理事务所(普通合伙) 61240 | 代理人: | 李艳春 |
地址: | 710054 *** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 联合 充填 采矿 高度 分析 方法 | ||
1.一种联合充填采矿中充填体的充填高度分析方法,其特征在于,所述充填体包括强度充填体和膨胀充填体,所述膨胀充填体位于强度充填体的上方,所述充填高度分析方法包括以下步骤:
步骤一、计算强度充填体固结过程中的收缩率;
步骤二、计算膨胀充填体固结过程中的膨胀系数;
步骤三、建立充填体的平衡方程、变形协调方程和物理方程;
所述充填体的平衡方程为:
F1+ρ1gh1+ρ2gh2-F2=0
其中,F1为顶板给充填体提供的约束反力,F2为底板给充填体提供的约束反力,ρ1为膨胀充填体的密度,ρ2为强度充填体的密度,h1为膨胀充填体的充填高度,h2为强度充填体的充填高度,且h1+h2=H,H为充填体的总高度,g为重力加速度;
所述充填体的变形协调方程为:
Δl1=Δl2
其中,Δl1为膨胀充填体的变形量,Δl2为强度充填体的变形量;
所述充填体的物理方程为:
其中,FN1为膨胀充填体的内部轴力,且FN1=F1+ρ1gx,FN2为强度充填体的内部轴力,且FN2=F2-ρ2gy,E1为膨胀充填体的弹性模量,E2为强度充填体的弹性模量,Rc为强度充填体的收缩率,Re为膨胀充填体的膨胀系数,x为改性镁渣基充填体的竖向收缩量,y为原镁渣基充填体的竖向膨胀量;
步骤四、对充填体进行临界状态分析,得到膨胀充填体和强度充填体的充填高度;
所述临界状态包括第一临界状态和第二临界状态;
当充填体处于第一临界状态时,所述膨胀充填体与顶板处于刚刚接触状态,顶板提供的约束反力F1为0,所述强度充填体的充填体积达到最大值,此时FN1=ρ1gx,根据充填体的平衡方程、变形协调方程和物理方程,得到:
g(2E1ρ1-E2ρ1-E1ρ2)h12+2E1[E2(Re+Rc)-gH(ρ1-ρ2)]h1-E1H(2E2Rc+ρ2gH)=0
令参数a1=g(2E1ρ1-E2ρ1-E1ρ2),
b1=2E1[E2(Re+Rc)-gH(ρ1-ρ2)],
c1=E1H(2E2Rc+ρ2gH);
解得膨胀充填体的充填高度强度充填体的充填高度
当充填体处于第二临界状态时,所述强度充填体的底面即将被压缩碎,应力达到强度充填体的抗压强度σ2=[σ2],所述强度充填体的充填体积达到最小值,此时FN2=[σ2]·1,F2=[σ2],根据充填体的平衡方程、变形协调方程和物理方程,得到:
g(E2ρ1-2E2ρ2+E1ρ2)h12+2[E1E2(Re+Rc)+[σ2](E1-2E2)+gH(E2ρ2-E1ρ2)]h1-E1H(2E2Rc+2[σ2]-ρ2gH)=0
令参数a2=g(E2ρ1-2E2ρ2+E1ρ2),
b2=2[E1E2(Re+Rc)+[σ2](E1-2E2)+gH(E2ρ2-E1ρ2)],
c2=E1H(2E2Rc+2[σ2]-ρ2gH);
解得膨胀充填体的充填高度强度充填体的充填高度
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安科技大学;西安弗尔绿创矿业科技有限责任公司,未经西安科技大学;西安弗尔绿创矿业科技有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202111316462.0/1.html,转载请声明来源钻瓜专利网。