岩石爆破破碎机理的分形研究
- 格式:pdf
- 大小:222.58 KB
- 文档页数:6
第16卷第5期1997年10月焦作工学院学报JOURNAL OF JIAOZUO INSTITUTE OF TECHNOLOGYVol116№15Oct. 1997
收稿日期:1997-03-25
岩石爆破破碎机理的分形研究员小有1 杨小林1 李国富2(1 焦作工学院建筑工程系 焦作454000 2 焦作矿务局科研所)
摘要 应用分形几何理论,研究了爆炸载荷作用下岩石块度的分布特征,发现在双对数坐标下,爆破块体粒径的质量累积百分数与粒径间呈直线关系,表明其块度分布具有分形结构。根据该直线的斜率b,由公式D=3-b可求得相应的分维。分维可作为描述块度分布特征的一个有序参量,其值反映了爆破块度分布的离散程度和块度的均匀程度。关键词 分形维数 爆破块度分布 块度预报预测中图法分类号 TD23511
第一作者简介 员小有,男,1965年生,工学硕士,工程师,主要从事岩土爆破工程的教学和科研工作,曾发表论文多篇。FractalStudiesonBlastingBrokingMechanismofRocks/YunXiaoyouetal1//Dept1ofConstructionEngineering,JiaozuoInstituteofTechnology,Jiaozuo454000Abstract Thefragmentdistributionofrock,whichisgroundbyblasting,arestudiedbyfractalgeome2try1Itisfoundthatthereisalinearrelationshipbetweentheaccumulativemasspercentageandthefrag2mentsizeinthelog-logplot,anditshowsthatthefragmentdistributionofrockhasfractalstructure1Theslope“b”intheplotcanbeusedtocalculatetherelativefractaldimension“D”bytheequationD
=3-b1Dmaybetakenasanorderedparameter,whichdescribesthefragmentdistributioncharacter2isticsofrock1ThevalueofDreflectsthescattertrendofdistributionaswellastheuniformityoffrag2ments1Keywords fractaldimension;fragmentdistributionofblasting;fragmentforecastingandcalculating
岩石爆破中块度分布是评定爆破质量的重要指标,它不仅影响矿山生产过程的铲装、运输和粗碎等工序的效率和生产成本,而且还影响定向爆破筑坝形成坝体的渗透性、沉陷性和稳定性。关于爆破块度的研究已经做了很多工作,采取了块度分布函数、爆堆摄影、大块统计、模糊数学、灰色理论、爆前节理裂隙分析及过程数值模拟等研究手段,现已能够确定现场爆破块度分布和评价爆破效果。但以上这些方法实施起来都比较麻烦,而且在爆破块度预报方面往往达不到要求。因此关于爆破块度的研究还需要不断开拓新的领域。绝大多数材料在形成、演化过程中由于受各种随机因素的影响,其内部结构存在着大量不同尺度、不同类型、并且呈无序分部的微观缺陷。在冲击荷载情况下,这些微观缺陷的周围将产生应力集中,使其生长、扩展甚至贯通,同时微观缺陷的无序性被强烈地放大。岩石类脆性材料爆破块度形成,显然是在爆炸载荷作用下产生大量断裂面相互贯通的结果。实验表明,岩体裂纹的分布及其演化在一定标度范围内具有很好的统计自相似性,
而这种自相似行为导致了材料宏观破碎后所形成碎块的块度分布具有分形特征,所以可用块体的分形构造来描述爆破岩块的形成过程。因此近年来岩石材料的原始损伤及其在爆炸载荷作用下的演化规律引起了爆破理论工作者的普遍关注。
1 块度呈分形分布的微观机理裂纹在几何上的非规则性是脆性材料的共同特征,但其分布可以看成是在不同层次、不同尺度并且按某种规则构造的裂纹嵌套,即第k+1步的小裂纹是在第k步的大裂纹群基础上生成的子群,从而构成一个统计自相似的分布系统,因而可以用线性分形来模拟这种非规则性现象。岩石材料在爆炸载荷作用下,裂纹将在极短的时间内进行触发-扩展-
触发链式生长模型的演化,从而使嵌套尺寸、层次急剧扩大和增加,而整体裂纹分布的分维D保持不变。文献〔3〕表明,在一定范围内,裂纹分布的统计自相似特征,在尺频关系上可用分形的基本关系--幂定律很好地描述。即根据分形的定义,为
n=ns
a
a0
篇。
-D(1)
式中,n———裂纹尺寸大于、等于a的裂纹数;ns———试样中裂纹总数;a———裂纹核尺寸或晶粒尺寸;D———裂纹分布的分维。以上是基于在断裂力学中常将裂纹简化为二维平面上一条直线来处理的考虑。脆性材料在爆炸冲击载荷作用下,裂纹将扩展、贯通,且伴有新裂纹产生,但由于其演化的统计自相似性,仍有类似于(1)式的关系。为了分析的方便,可把爆炸加载过程人为地按时间划分为若干个等长的阶段,假定在爆炸加载过程中,材料中的低层次大裂纹比高层次小裂纹优先贯通断裂,那么在第k阶段有nk个大于、等于尺寸ymin的贯通裂纹断裂,则
nk=nsk
ymin
a0
tc
-D(2)
式中,ymin———此阶段材料的最小贯通裂纹尺寸;nsk———此阶段试样中裂纹总数。由此产生的碎块总数Nk为
Nk=2nk-(nk-1)=nk+1≈nskymina0im-D (当nk
很大时)(3)
ymin可看作此阶段产生的最小颗粒尺寸x
min,
因此
nk=Cx-Dmin(4)
进一步可推出在爆炸载荷作用到此阶段时的块度分布为N=Cx-D (x≥x
min)(5)
式中,N———块度尺寸≥x的碎块数;C———比例常数。
・73・第5期员小有等:岩石爆破破碎机理的分形研究 此阶段块度小于x的碎块数占总碎块数的比率Y
N(x)为
YN(x)=1-Cx-DNk
(6)
从上面的推导可以看出,脆性材料在爆炸载荷作用的每个阶段后,其块度分布均呈分形分布,并且只与材料内部的原始裂纹有关。即在爆炸载荷作用过程的不同阶段,产生碎块块度分布的分维值是不变的,都等于材料原始裂纹分布的分维。
2 块度分布分维的推导根据公式(6)有dYN(x)∝x-D-1dx(7) 尺度在x~x+dx之间的碎块数目dN表示为dN=NkdYN(x)(8) 若假定岩石的密度为ρ,则尺寸在x~x+dx之间的碎块质量dW可表示为dW=ρKVx3NkdYN(x)(9)式中,KV———体积形状因子。另外dW还可写成下述形式dW=WdYW(x)(10)式中,W———爆破碎块的总质量;Y
W(x)———块度小于x的碎块质量占总质量的比
率。由式(9)和(10)得YW(x)=ρNk(KV/W)x3・dYN(x)(11) 将式(7)代入式(11)可得dYW(x)∝x2-Ddx(12) 对上式积分得YW(x)∝x3-D
(13)
所以,如果爆破块度的分布满足
YW(x)∝xb(14)
即在双对数坐标下Y
W(x)与x间存在直线段,这就表明爆破块度的分布为分形分布,
根据其斜率b即可求得块度分布的分维值为D=3-b(15)
3 块度分布分维的计算我们使用不同强度的砼试件,采用不同的不偶合装药系数进行爆破,对所形成碎块的块度分布进行统计分析后发现,同种强度试件不同装药条件下的爆破块度,在一定尺寸范围内,其块度分布呈现出分形分布特征。部分试验数据及计算结果如图1、2、3和表1所示,其中分维值D及相关系数r由最小二乘法求得。
・83・ 焦 作 工 学 院 学 报1997年图1试件强度S=19170MPa图2试件强度S=40180MPa图3
试件强度S=6213MPa
表1 不同不偶合系数下爆破块度分布和分形维数
Kd抗压强度MPa筛 下 累 计 百 分 数<20<40<60<80<100<120bDrα=σx
11001917401862138184717771421716510154101442811418.0317159391633512434177721427015070108100100100018693019131018754211307210869211246019730199301958015832017213017312
113319174018621391627106311118.7312104517832183191509141551662912112183781286919591173100100100110070018910018724119930211090211276019850196701972016569018200116737
1160191740186213915061299165161491510916125221582813938134461893116946191711257118074141100100100018541018972019099211459211028210901019850197801968015588017908015153
2100191740186213812291266182151801918717171221903119731187431904117842197741406615063144100100100019004019325019228210994210675210772019630198601996015668015519015888
216719174018621371241011711108161431713816178341542716829191481433714049133701117416968127100100100019556019424019657210444210576210343019930197301979017612016575016921
4 结果讨论及分维意义411 G-G-S分布均匀性系数的分维性质从表中可以看出,各相关系数之值均在1~0199之间,这种强相关性说明爆破块度分布的分形特征是客观存在的。事实上,对于大多数脆性材料,实验表明,在爆炸载荷作用下形成的块度分布,多符合R-R(Rosin-Rammler)和G-G-S(Gate-Gaudin-Scuh2mann)分布。一般认为R-R分布倾向于粗粒端,G-G-S分布倾向于细粒端。若将二
・93・第5期员小有等:岩石爆破破碎机理的分形研究