Numerical simulation on the cracking process of columnar jointed rock masses in contour blasting
-
摘要: 以白鹤滩水电站柱状节理玄武岩为研究对象,根据现场的柱状节理分布统计资料生成数值模型,模拟边坡开挖轮廓爆破下柱状节理岩体的开裂过程。结果表明:轮廓爆破时,柱状节理面对岩体裂纹的发展起着导向和控制作用;爆炸荷载作用下炮孔周围的节理面首先张开,然后裂缝沿着炮孔连心线附近的节理面扩展,最后形成贯通的由大量柱状节理面连接而成的裂缝;预裂爆破条件下,预裂缝宽度先增大后减少,最后达到稳定的开度,而且预裂爆破在柱状节理岩体中产生的裂纹范围比光面爆破产生的大。Abstract: With Baihetan columnar joints basalt as the research object, numerical model is generated according to the local distribution statistics. Based on numerical analysis, the cracking process of columnar jointed rock masses with the slope excavated by contour blasting is simulated. The results show that the columnar joints control and lead the generation of cracking in contour blasting. The joints around the hole are firstly opened by the explosion load. Then they extend along the joints near the defiled line of holes. Finally the cracks formed by a large number of columnar joints are connected. The process of the width of pre-splitting crack is that it reduces to stabile opening after increases under the condition of pre-splitting blasting. And the range of cracks in pre-splitting blasting is larger than these in smooth blasting.
-
Key words:
- mechanics of explosion /
- cracking process /
- contour blasting /
- rock /
- columnar joints
-
智能引信通常用于判断战斗部最佳起爆位置、最优起爆时机以及确保能够安全引爆战斗部,而钽电容在智能引信中起着贮能、点火、解耦、延时等作用。因此,钽电容的质量对智能引信的起爆控制系统起着决定性的作用。当炮弹发射时,会产生突发性高机械冲击载荷,钽电容经受高冲击载荷作用时会发生特性变化。钽电容在高冲击载荷下的参数变化特性是决定智能引信实现其功能的关键因素之一。
针对钽电容在环境载荷下的参数变化特性已做了许多研究。徐建军研究了钽电容在高温下的漏电流变化现象[1], Alexander研究了钽电容在机械振动、反偏电压作用下的电参数变化特性[2-3], 程融发现钽电容在冲击应力下易出现瞬时短路的现象[4], 何荣华等对高过载下军用电容的失效进行了探讨,并分析了电容失效对引信的影响[5]。
虽然关于钽电容的参数变化问题已做了大量研究,但有关参数变化的影响因素主要集中在振动、温度等环境载荷方面,而针对冲击载荷对钽电容影响的研究比较少。即使有关于钽电容经受冲击载荷作用的研究也主要集中在实验现象上,而缺少对引起参数变化的内部机理的分析。本研究首先通过实验手段分析了钽电容在冲击载荷作用下电容量、漏电流的变化特性,然后对电参数变化特性分析其内在机理。
1. 实验设计
1.1 钽电容结构
固态钽电容主要包括阳极Ta、电介质Ta2O5、阴极MnO2。阳极Ta是由钽粉压制烧结而成,因此具有多孔性,有助于增大内部接触面积,提高电容容值。Ta2O5介质通过阳极氧化,形成于阳极Ta的外表面。将钽棒浸入锰盐溶液中并热分解,形成阴极MnO2。钽电容的结构如图 1所示。
1.2 钽电容实验测试
钽电容用于智能引信中,需要经受瞬间强冲击载荷的作用,在冲击载荷作用下钽电容性能的变化为其电容量和漏电电流的变化。马歇特锤试验装置能够产生瞬态冲击载荷,载荷范围在0~50 000 g之间,能够在试验中模拟钽电容的实际使用工况。因此,使用马歇特锤击设备测试钽电容的电容量和漏电流。实验过程中,选取4种不同类型的钽电容进行测试,将被测钽电容焊接在电路板上,电路板采用螺纹连接与击锤固定在一起。击锤安装于固定在偏心轮的锤杆上。实验时使用偏心轮将击锤升至一定高度,用卡齿固定。卡齿放开后,依靠偏心轮下方的重锤作为动力,击锤击打在铁锭上产生冲击载荷。使用压电式加速度计实时记录冲击过载值。使用漏电流测试仪(Keithley 6487 Picoammeter)监测钽电容漏电流随冲击作用的变化,动态信号分析仪记录皮安表输出。使用阻抗分析仪(Agilent 4294A)记录电容量的变化。传感器型号为PCB 350B21,灵敏度为0.05 mV/g,量程为±100 000 g。在电容充电5 min后进行实验,实验过程中,采用步进冲击方法,逐步增大冲击g值,观察钽电容电参数的变化。实验系统如图 2所示。钽电容参数及测量范围如表 1所示。
表 1 钽电容参数及测量范围Table 1. Tantalum capacitors and sensor电容 型号 数量 电容量 漏电流 测量范围/μF 测量精度/μF 误差/% 测量范围/mA 测量精度/mA 误差/% 16 V22 μF
16 V47 μF
25 V22 μF
25 V47 μF16K226
16K476
25K226
25K4765
5
5
50~47 ±0.01 99.98~100.02 0~2 ±0.01 99.5~100.5 2. 实验结果
图 3为钽电容的额定电容随冲击载荷变化的情况。从实验结果发现,电容量随冲击载荷的变化而发生变化。对于16 V22 μF钽电容,在施加1 000 g冲击后,电容量有0.1 μF的升高;持续增加冲击载荷,在3 000 g时电容量增加0.15 μF;当冲击值达到10 000 g时,电容量升高0.4 μF。对于16 V47 μF钽电容,在施加1 000 g冲击后,电容量有0.15 μF的升高;在3 000 g时电容量增加0.2 μF;当冲击值达到10 000 g时,电容量升高0.45 μF。对于25 V22 μF钽电容,与冲击载荷对应的电容量变化为0.12、0.17和0.45 μF;对于25 V47 μF,电容量对应的变化为0.13、0.19和0.50 μF。冲击完成后,电容量恢复至原来的额定值。
图 4为漏电流随冲击载荷的变化情况。电容接通电源后,涌入的电流为充电电流,充电电流开始很大,随着时间延长(≥5 min)而下降稳定在远小于额定电流的特定值,此时的电流为漏电流。实验过程中,冲击值从1 000 g开始逐渐增加。从实验结果发现,漏电流变化趋势与冲击载荷变化相一致。对于16 V22 μF电容施加1 000 g的冲击后,漏电流增加一个数量级,当冲击值增加到3 000 g时,漏电流约为0.15 mA,随着冲击值的增大,漏电流逐渐增加。当冲击值达到10 000 g时,漏电流达到0.25 mA。对于16 V47 μF电容,施加1 000 g的冲击后,漏电流达到0.28 mA,当冲击值增加到3 000 g时,漏电流约为0.4 mA,当冲击值达到10 000 g时,电容漏电流达到0.6 mA。对25 V22 μF电容施加1 000 g的冲击后,漏电流约为0.13 mA,当冲击值增加到3 000 g时,漏电流约为0.44 mA,当冲击值达到10 000 g时,漏电流达到0.9 mA。对于25 V47 μF电容,施加1 000 g的冲击后,漏电流达到0.4 mA,当冲击值增加到3 000 g时,漏电流约为0.6 mA,当冲击值达到10 000 g时,电容漏电流达到1.7 mA。
如图 4所示,当冲击载荷结束后,4组钽电容的漏电流经过一段时间后恢复到原始水平。即使漏电流达到很高的量级,在钽电容未被击穿时,冲击过后,漏电流依然能恢复到原来的水平。这种特性称为钽电容的自愈特性。
智能引信能够引爆战斗部,要求钽电容中存储的能量必须大于电雷管的安全起爆能量。由于钽电容在冲击载荷作用下漏电流增大,则电容中储存的能量以漏电流的形式出现损耗。当钽电容剩余的能量大于电雷管需要的发火能量时,才能可靠引爆电雷管,可靠实现智能引信的功能。
CU2/2−[U(t)Ileatb+wother]≥I2fRdtb (1) 电容充满电后存储的总能量为CU2/2, 电雷管爆炸所需要的能量为If2Rdtb,If=0.45 A为最小发火电流,Rd为电雷管阻抗,tb为点火时间。U(t)Ileatb+wother为漏电流和其他元件损耗的能量,U(t)为电压变化值。为使电雷管能够引爆,4种电容允许的最大漏电流为1.1、2.4、3.5和7.3 mA。
针对钽电容允许的最大漏电流,继续加大冲击应力,得出与钽电容对应的临界失效冲击载荷,如图 6所示。
由图 6可知,钽电容漏电流在冲击载荷作用下呈指数型升高。钽电容临界失效载荷及钽电容漏电流在冲击载荷下的变化模型如表 2所示。
表 2 漏电流-冲击模型Table 2. Numerical model of leakage current-shock load电容 漏电流-冲击模型 最大漏电流/mA 临界冲击载荷/(104g) 16 V22 μF I=A1+A2exp(α g) A1=1.09×10-4,A2=8.17×10-6,α=3 429 1.1 1.65 16 V47 μF I=B1+B2exp(β g) B1=2.14×10-4,B2=5.28×10-5,β=4 506 2.4 1.72 25 V22 μF I=C1+C2exp(γ g) C1=1.17×10-4,C2=8.98×10-5,γ=4 427 3.5 1.74 25 V47 μF I=D1+D2exp(λ g)D1=2.14×10-4,D2=5.28×10-5,λ=4 506 7.3 1.85 3. 电参数变化机理分析
由实验结果可知,冲击会引起钽电容电参数的变化:电容量的升高及漏电流随冲击载荷的升高呈指数型升高。分析钽电容电参数变化,可得冲击载荷对钽电容的作用机理为:弹性变形范围内,冲击引发的钽电容变形导致电容的变化;冲击引起介质层微裂缝以及冲击引起陷阱能级的升高导致漏电流增大。
3.1 冲击引起的钽电容变形
由于钽电容是由钽颗粒压制烧结而成,因此钽电容可以看做多个钽颗粒电容的并联。图 7为单个钽颗粒模型。钽电容在受到冲击载荷作用时,赋能颗粒发生变形,在弹性变形范围内,由球形颗粒变形为椭球形颗粒,然后恢复原状。
图 7(a)所示的球形电容的电容量为:
Cs=4πεRr/(R−r) (2) 式中:ε为介电常数, F/m;R为Ta2O5外径, mm;r为Ta半径, mm。图 7(b)所示椭球型电容的电容量为[7]:
Ce=QU=16πεr√R/w/(arctanR−w2r√R/w−arctanr−rw/R2r√R/w) (3) 式中:Q为电荷量, C;U为内外极板之间的电压, V;w为变形量, mm。
在冲击过程中,电容的变化量为:
ΔC=Ce−Cs=4πε{4r√R/w/(arctanR−w2r√R/w−arctanr−rw/R2r√R/w)−Rr/(R−r)} (4) 冲击作用下球形颗粒的变形见图 8。根据赫兹接触理论得,在力F(t)作用下,在弹性变形范围内,两球心O1、O2相互接近的距离为:
δ3=916R1+R2R1R2(1−μ21E21+1−μ22E22)F2(t) (5) 式中:R1、R2为两球的半径,μ1、μ2、E1和E2分别为泊松比及弹性模量,GPa。
则电容球形颗粒在受到冲击作用时,球心相互接近的距离为:
δ3=9F2(t)(1−μ2)/(4RE2) (6) 则单个球形颗粒变形量为:
δ=123√9F2(t)(1−μ2)/(4RE2) (7) 在冲击作用下,电容值的变化量为:
{ΔC=Cs−Ce=4πε[Rr/(R−r)−4r√R/w/(arctanR−w2r√R/w−arctanr−rw/R2r√R/w)]w=123√9F2(t)4(1−μ2RE2) (8) 图 9为钽颗粒电容随冲击载荷的变化。
假设钽电容由半径均匀一致的球体钽颗粒电容组成,则每个钽电容所包含的球形电容的数量N为:
N=V/∑Vi (9) 式中:V为未封装的钽电容体积,Vi为Ta、Ta2O5、MnO2的体积。不同类型钽电容包含球形钽颗粒电容的数量及电容变化值如表 3所示, 图 10为不同类型电容与冲击载荷的关系。
表 3 冲击载荷下钽电容容量的改变量Table 3. Capacitance variation for different capacitors电容 体积/mm3 数量/105 电容变化/μF 1 000 g 3 000 g 10 000 g 16 V22 μF 2.5×2.0×1.5 1.07 0.06 0.12 0.24 16 V47 μF 4.0×3.0×0.8 1.37 0.08 0.16 0.36 25 V22 μF 3.5×2.8×1.5 2.10 0.12 0.24 0.46 25 V47 μF 4.0×3.0×1.7 2.92 0.18 0.32 0.76 由表 3可知,计算值与实验值存在偏差。分析可知造成数值偏差的原因为:由于加工钽电容时,在外层添加石墨和银层,而在计算时,忽略石墨和银层,造成数值出现偏差;钽电容变形时,不仅在纵向出现变形,横向同样存在微小形变,计算时未考虑横向变形因素。
3.2 冲击引起的漏电流变化
冲击引起的漏电流变化是由两种机制决定的:介质层在冲击下产生微裂缝;冲击下电子陷阱的产生,Ta—O化合键的断裂以及氧气空穴的产生。陷阱浓度的增加使参与导电的电子增多,漏电流从而增大。
实验中发现,在钽电容未被击穿时,冲击过后,漏电流依然能恢复到原来的水平。这种现象可由钽电容的自愈机理解释。在冲击载荷作用下,钽电容介质层产生“薄弱”区,使该区域漏电流增大。持续增大的电流使局域温度升高,当温度达到400 ℃时,导电性良好的MnO2转变为导电性能差的Mn2O3,覆盖到缺陷区域上,使“薄弱”区域处的漏电流恢复至原来的量级。
4. 结论
利用马歇特落锤实验对固体钽电容进行高过载测试,从测试结果发现,冲击载荷对固体钽电容电参数的影响主要表现为电容量的升高以及漏电流的增大。电容量随着冲击载荷的增加而增加,钽电容漏电流随冲击载荷的升高呈指数型升高,冲击过后,电容量恢复额定容值,经过一段时间,钽电容漏电流同样恢复至原来的量级。同时,分析钽电容电参数变化原因,得出引起钽电容电参数变化的机制为:冲击引起的钽电容的弹性变形使电容量增大,冲击引起的介质层Ta2O5中的微裂缝以及冲击引发的介质层陷阱浓度的变化。钽电容球形模型在冲击载荷作用下发生弹性变形,使电容量发生改变。在冲击载荷的作用下,介质层缺陷处容易引发微裂缝,形成导电微沟道,使漏电流增大。另外,由于冲击载荷的作用使介质层中化合键产生断裂导致电子传导机制中参与导电的电子增多,同时由于氧空穴在Ta2O5和MnO2界面处的积累,使界面产生极化,界面势垒降低进而导致注入的电子数量增多,电流密度增大,漏电流升高。
-
表 1 岩石力学参数
Table 1. Mechanical parameters of the rock
ρ/(kg·m-3) E/GPa μ pR/MPa α/(°) σs/MPa 2 800 60 0.2 0.9 45 5 表 2 节理面力学参数
Table 2. Mechanical parameters of the joint mass
Kn/GPa Kt/GPa αi/(°) pJ/MPa σt/MPa 52.9 12.8 28 0.75 0.78 -
[1] 张倬元, 王士天, 王兰生.工程地质分析原理[M].北京: 地质出版社, 1994. [2] Mallet R. On the origin and mechanism of production of the prismatic(or columnar)structure of basalt[J]. Philosophical Magazine, 1875, 4(50): 122-135, 201-226. doi: 10.1080/14786447508641280 [3] 刘海宁, 王俊梅, 王思敬.白鹤滩柱状节理岩体真三轴模型试验研究[J].岩土力学, 2010, 31(1): 163-171. http://www.cnki.com.cn/article/cjfdtotal-ytlx2010s1028.htmLiu Hai-ning, Wang Jun-mei, Wang Si-jing. Experimental research of columnar jointed basalt with true apparatus at Baihetan Hydropower Station[J]. Rock and Soil Mechanics, 2010, 31(1): 163-171. http://www.cnki.com.cn/article/cjfdtotal-ytlx2010s1028.htm [4] 石安池, 唐鸣发, 周其健.金沙江白鹤滩水电站柱状节理玄武岩岩体变形特性研究[J].岩石力学与工程学报, 2008, 27(10): 2079-2086. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200810017.htmShi An-chi, Tang Ming-fa, Zhou Qi-jian. Research of deformation characteristics of columnar jointed basalt at Baihetan Hydropower Station[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(10): 2079-2086. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200810017.htm [5] 徐卫亚, 郑文棠, 石安池.水利工程中的柱状节理岩体分类及质量评价[J].水利学报, 2011, 42(3): 262-279. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=slxb201103002Xu Wei-ya, Zheng Wen-tang, Shi An-chi. Classification and quality assessment of irregular columnar jointed basaltic rock mass for hydraulic engineering[J]. Journal of Hydraulic Engineering, 2011, 42(3): 262-279. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=slxb201103002 [6] 徐卫亚, 郑文棠, 宁宇, 等.柱状节理坝基岩体三维各向异性数值分析[J].岩土力学, 2010, 31(3): 949-955. http://d.wanfangdata.com.cn/Periodical/ytlx201003048Xu Wei-ya, Zheng Wen-tang, Ning Yu, et al. 3D anisotropic numerical analysis of rock mass with columnar for dam foundation[J]. Rock and Soil Mechanics, 2010, 31(3): 949-955. http://d.wanfangdata.com.cn/Periodical/ytlx201003048 [7] 朱珍德, 秦天昊, 王士宏.基于Cosserat理论的柱状节理岩体各向异性本构模型研究[J].岩石力学与工程学报, 2010, 29(2): 4068-4076. http://www.cqvip.com/QK/96026X/2010A02/35616589.htmlZhu Zhen-de, Qin Tian-hao, Wang Shi-hong, et al. Study of anisotropic constitutive model of columnar jointed rock masses based on cosserat theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(2): 4068-4076. http://www.cqvip.com/QK/96026X/2010A02/35616589.html [8] Wang Z L, Konietky H. Modeling of blast-induced fractures in jointed rock masses[J]. Engineering Fracture Mechanics, 2009, 76: 1945-1955. doi: 10.1016/j.engfracmech.2009.05.004 [9] 刘红岩, 秦四清, 杨军.爆炸荷载下岩石破坏的数值流形方法模拟[J].爆炸与冲击, 2007, 27(1): 50-56. doi: 10.11883/1001-1455(2007)01-0050-07Liu Hong-yan, Qin Si-qing, Yang Jun. Simulation of rock failure by numerical manifold method under blasting load[J]. Explosion and Sock Waves, 2007, 27(1): 50-56. doi: 10.11883/1001-1455(2007)01-0050-07 [10] 成艳.长白山第四纪玄武岩隧洞光面爆破研究[D].吉林: 吉林大学建设工程学院, 2012. [11] 齐景嶽, 刘正雄.隧道现代爆破技术[M].北京: 中国铁道出版社, 1995. [12] Ding D, Zhu C. Estimaing amount of explosive for fracture plane control blasting with notched boreholes[J]. Transactions of NFsoc, 1993, 3(2): 18-21. http://www.cqvip.com/QK/85276X/199302/4001541744.html [13] 邴凤山.中国水力发电年鉴2010[M].北京: 中国电力出版社, 2011. [14] 朱道建, 杨林德, 蔡永昌.柱状节理岩体压缩破坏过程模拟及机制分析[J].岩石力学与工程学报, 2009, 28(4): 716-724. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200904012.htmZhu Dao-jian, Yang Lin-de, Cai Yong-chang. Simulation of compressive failure of columnar jointed rock mass and its failure mechanism analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(4): 716-724. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200904012.htm [15] 金李, 卢文波, 陈明, 等.节理岩体的爆破松动机理[J].爆炸与冲击, 2009, 29(5): 274-280. doi: 10.11883/1001-1455(2009)05-0474-07Jin Li, Lu Wen-bo, Chen Ming, et al. Mechanism of jointed rock loosing under blasting load[J]. Explosion and Shock Waves, 2009, 29(5): 274-280. doi: 10.11883/1001-1455(2009)05-0474-07 期刊类型引用(8)
1. 叶紫阳,吴红波,杨仕春,黄菓树,李天浩,孙翼,马成帅,任梦雨. 延期状态下电子雷管电子控制模块的高过载加载实验. 高压物理学报. 2025(01): 100-110 . 百度学术
2. 王家乐,李洪伟,王小兵,梁昊,周恩,苏洪,赵金耀. 冲击载荷作用下钽电容的电压瞬变特性及微观机理. 爆炸与冲击. 2024(04): 17-25 . 本站查看
3. 李洪伟,王家乐,郭子如,鹿华雨,柳祥威. 冲击载荷下液态铝电解电容失效规律研究. 安徽理工大学学报(自然科学版). 2024(02): 36-41 . 百度学术
4. 刘波,杨荷,赵慧,吴学星,李华梅,程祥利. 单轴静压条件下高压多层陶瓷电容的容值变化. 兵工学报. 2023(06): 1858-1866 . 百度学术
5. 贾丰州,牛少华,孙远程,陶逢刚,刘天国. 冲击载荷作用下的固体钽电容力-电响应特性. 探测与控制学报. 2022(05): 20-25 . 百度学术
6. 李波,张宏翔,李世中,张亚. 极端温度条件典型引信电参数变化规律和失效机理. 装备环境工程. 2022(11): 34-40 . 百度学术
7. 伍权,郑跃,舒启佳,田东斌. 钽电容器高频振动结构可靠性分析与设计. 电子元件与材料. 2021(04): 339-343+349 . 百度学术
8. 张美云,石庚辰,刘强,张力,邱强强. 侵彻硬目标引信侵彻信号获取与处理技术. 制导与引信. 2020(03): 1-10+24 . 百度学术
其他类型引用(2)
-