线性聚能装药金属射流侵彻岩石全过程的数值模拟.docxVIP

线性聚能装药金属射流侵彻岩石全过程的数值模拟 矿岩切割爆破过程中的力学性质研究 聚能服装是一种使用包装的“孔”，以改善炸药爆炸对局部造成的破坏。20世纪末, 国内一些科技工作者开始将聚能装药应用于矿岩切割爆破, 并进行了理论和实验研究, 说明聚能药包爆炸切割岩石的可能性和实用性, 然而对聚能药包爆炸切割岩体的物理过程和力学条件还缺乏必要的研究。本文运用改进的MOCL程序, 数值模拟线性聚能装药金属射流的形成及射流侵彻岩石的全过程, 旨在为线性聚能装药切割岩石机理分析与应用研究提供参考。 1 计算方法 1.1 爆进行空气状态方程的一般规定 聚能射流的形成过程可用二维不定常弹塑性流体力学方程组描述。在柱对称情况下, 使用Euler坐标的基本方程如下: 质量守恒方程 ?ρ?t+uz?ρ?z+ur?ρ?r+ρ(?uz?z+?rαurrα?r)=0 动量守恒方程 ρ(?uz?t+uz?uz?z+ur?uz?r)=-?Ρ?z+?Szz?z+?(rαSrz)rα?rρ(?ur?t+uz?ur?z+ur?ur?r)=-?Ρ?r+?Srz?z+?(rαSrr)rα?r-αSθθr 能量守恒方程 ρ(?e?t+uz?e?z+ur?e?r)=-Ρ(?(rαur)rα?r+?uz?z)+Szz?uz?z+Srr?ur?r+Srz(?uz?r+?ur?z)+αurSθθr 式中,t为时间;r为径向欧拉坐标;z为轴向欧拉坐标;ur为径向质点速度;uz为轴向质点速度;Srr为径向偏应力;Sθθ为环向偏应力;Szz为轴向偏应力;Srz为径向或轴向剪应力;ρ为密度;e为比内能。 屈服条件:F=Y0/√3(S2rr+S2zz+S2rz+SrrSzz), 当F1时物质为弹性,F≤1时物质进入塑性状态。 熔化条件:当e≥em时, 物质熔化, 物质强度Y=0;当eem时,Y=Y0(1-e/em) 状态方程:在聚能射流的数值模拟中, 涉及到炸药、金属和空气3种介质。 对于空气, 近似采用理想气体状态方程,P= (kα-1)ρe,kα为空气的等熵指数。 对于炸药, 通常采用JWL状态方程, 其一般形式为 ΡS=Aexp(-R1V)+Bexp(-R2V)+CV-(w+1)ES=AR1exp(-R1V)+BR2exp(-R2V)+CwVw 式中,PS为压力;ES为爆轰产物的内能;V为爆轰产物的相对体积 (即爆轰产物体积与初始体积之比) ;下标S表示等熵过程;A、B、C、D、R1、R2、w为待定常数。 对金属采用Mie-Grüneisen状态方程, 其一般形式如下, Ρ=ΡΗ(1-Γμ2)Γρ(e-e0) 式中, ΡΗ={k1μ+k2μ2+k3μ3? ;μ=ρρ0-1;Γ为Grüneisen系数;k1、k2、k3为材料特性常数;ρ为当地密度;ρ0为初始密度;e0为单位质量的初始比内能;e为当地单位质量的比内能。 对于岩石, 压力P和体积应变εPv的关系为P=keαεPV,k,α为常数。 1.2 计算值的方法 1.2.1 轴向网向网内网格的网格分布 采用柱坐标系, 将计算空间域分成固定的以 (i,j) 编号的Euler矩形网格, 采用非等间距的网格分布, 各网格在轴向 (z方向) 和径向 (r方向) 的长度分别为Δzi、Δrj, 如图1所示。网格中心为 (zi,rj) 。密度为ρ,z向和r向的速度分量为uz和ur, 压力为P, 比内能为e, 总能为E, 定义在网格中心, 是网格量。人工粘性则定义在网格边界中点, 是格变量。 1.2.2 网格速度、能量和质量控制方法 计算采用“和”分裂差分格式, 对时间采用向前差分, 对空间采用中心差分。压力效应步 (Lagrange步) , 考虑压力梯度和偏应力的影响, 得到网格的中间速度和内能。输运步 (Euler步) , 质量、动量和能量在网格上的重新分配, 得到新的网格速度和网格内各介质的质量及内能。于是通过压力效应步 (Lagrange步) 和输运步 (Euler步) , 可得该问题的一阶精度显式格式。 1.2.3 网格线界面处理 由于计算域中含有炸药、空气、金属3种介质, 所以在计算域中, 同一网格内可能同时存在两种或3种介质。对于含两种物质的混合格, 运用网格线示踪点法界面处理技术, 即在整个计算过程中, 对每种物质用一串首尾相接的线段来逼近物质表面的曲线, 且每条线段的端点 (即示踪点) 均定义在网格线上, 根据网格中示踪线的位置来计算网格间材料的输运量。这样, 显示材料的分界面由Lagrange型示踪点跟踪物质变为沿Euler网格线的示踪点跟踪物质。 1.3 轴、固壁、自由面或远场的连续边界 在数值计算中, 只能取有限的计算域来模拟具体的物理问题。在计算域的四周和计算域的内部都可能存在有边界, 边界可以是对称轴、固壁、自由面或远