基于代表性体积单元模型的矿石组分界面破碎特征

蔡改贫 郝书灏 余成 宣律伟

蔡改贫, 郝书灏, 余成, 宣律伟. 基于代表性体积单元模型的矿石组分界面破碎特征[J]. 高压物理学报, 2022, 36(3): 035302. doi: 10.11858/gywlxb.20210896
引用本文: 蔡改贫, 郝书灏, 余成, 宣律伟. 基于代表性体积单元模型的矿石组分界面破碎特征[J]. 高压物理学报, 2022, 36(3): 035302. doi: 10.11858/gywlxb.20210896
CAI Gaipin, HAO Shuhao, YU Cheng, XUAN Lyuwei. Fracture Characteristics of Ore Components Interface Based on Representative Volume Unit Model[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 035302. doi: 10.11858/gywlxb.20210896
Citation: CAI Gaipin, HAO Shuhao, YU Cheng, XUAN Lyuwei. Fracture Characteristics of Ore Components Interface Based on Representative Volume Unit Model[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 035302. doi: 10.11858/gywlxb.20210896

基于代表性体积单元模型的矿石组分界面破碎特征

doi: 10.11858/gywlxb.20210896
基金项目: 国家自然科学基金(51464017);江西省重点研发计划(20181ACE50034)
详细信息
    作者简介:

    蔡改贫(1964-),男,博士,教授,主要从事散体-矿冶智能装备、物料破碎动力学关键技术和装备开发等研究. E-mail:664431653@qq.com

  • 中图分类号: O347; TD912

Fracture Characteristics of Ore Components Interface Based on Representative Volume Unit Model

  • 摘要: 矿物组分的界面破碎特征一般指黏结界面在外部载荷作用下产生的应力、应变等,其对研究矿物组分解离、提高矿石破碎效率具有重要意义。针对矿石内部组分矿物聚集、有用矿物非均匀分布的特点,开展了岩石内部岩相分析实验和矿物界面原位加载实验,在此基础上,利用非线性、多尺度建模平台DIGIMAT构建符合组分矿物微观结构的代表性体积单元(RVE)模型,并通过DIGIMAT-ABAQUS耦合开展矿石RVE模型原位破碎模拟。结果表明:(1)黑钨矿石中有用矿物以颗粒状分布在矿石内部,主要分布在石英矿物内及其与硅质岩矿物黏结界面处;(2)不同组分黏结界面的力学性质存在差异,且与组成矿物的物理属性、形态特征等相关,石英-硅质岩界面的最小破碎应力(界面发生破坏时的应力)范围为1.1785~1.4826 GPa,石英-钨界面的最小破碎应力范围为1.3355~1.5420 GPa;(3)加载速率为0.010或0.005 kN/s时,矿石破碎峰值应力无明显变化,但其对矿石内部形变的影响较大,且加载速率为0.010 kN/s时,强化阶段应力易突然降低并不断波动;(4)原位载荷产生的破坏主要发生在载荷作用区域边界,且两组界面中石英矿物破碎力学性能参数大于钨矿物和硅质岩矿物,即界面组成矿物中石英矿物优先形成破坏。

     

  • 全浸水带间隙发射作为一种新的水下发射方式[1],通过间隙燃气在膛内贴壁运动,卷吸回流后逐渐汇聚成弹前气幕,排出弹前水柱,将射弹在膛口的发射环境由水介质转化为气体介质。当气体射流流出枪口后,膛口处射流迅速膨胀成球形气体空腔,射弹穿过气体空腔与水介质接触,产生强烈的冲击载荷。水下高速射弹的弹体主要由硬铝合金尾杆和钨合金头部组成,两者镶嵌连接,连接强度有限,入水瞬间弹头会承受巨大冲击,因此射弹的入水冲击载荷成为水下射弹设计中的一个重要问题。

    Karman[2]最先开始对入水冲击现象进行研究,采用动量定理并引入附加质量的概念,推导出入水冲击载荷的计算公式。Wagner[3]将Karman的方法理论化,提出了近似平板理论及自相似解法,得出了冲击压力在结构沾湿面的分布情况,使理论分析更加符合实际情况,为后来学者的理论研究奠定了基础。在国内,秦洪德等[4]、王永虎等[5]对入水冲击问题的现状和进展进行了详细的分析。宋保维等[6]基于不可压缩的非定常势流理论,建立了空投水雷入水冲击计算的数学模型。卢炽华等[7]利用不同浸深的附加质量,对刚性细长体斜姿态落水冲击进行建模,得出其入水角很小,会使弹体处在最危险的状态。王永虎等[8-9]先后对刚性尖拱体垂直姿态高速入水和斜入水的冲击理论进行了建模和仿真。魏卓慧等[10]建立了刚性截锥弹体垂直入水冲击载荷的数学模型,并对其进行了数值计算。陈诚等[11]对超空泡航行器倾斜入水冲击载荷进行了试验研究,得出了峰值时刻的阻力系数。朱珠等[12]利用商业软件FLUENT,建立了柱体回转体高速入水冲击的数值模拟模型,得到了速度对入水冲击载荷的影响规律。然而以上研究主要针对由空中入水的冲击载荷分析,对于水下入水冲击问题研究较少。本工作在此前提条件下,计算分析全水下发射高速射弹入水的冲击载荷,这对于水下发射武器研究具有一定的现实意义。

    本研究拟建立锥形弹体水平及斜入水的冲击载荷理论模型,模型中考虑弹体重力、弹体浮力、附加质量、弹头锥角及入水攻角的影响,对不同头部结构参数的锥形弹体以不同入水速度入水的冲击载荷进行计算,分析入水速度、弹头锥角和入水攻角对冲击载荷的影响。研究结果对于弹体入水冲击载荷的预测及全水下发射方式发射的射弹头部结构设计具有参考价值。

    假设射弹为刚体,不考虑射弹入水时空泡的影响,根据动量定理,射弹高速入水冲击时的动量方程为

    Mv0=(M+m)v+Mgt+Fbt+Fdt
    (1)

    式中:M为射弹的质量,m为射弹的附加质量,Fd为射弹入水时受到的阻力,Fb为射弹所受的浮力。

    对式(1)等号两边进行微分,得到射弹入水冲击时的动力学方程

    (M+m)d2hdt2+dmdtdhdt+12ρCdA(dhdt)2=0
    (2)

    式中:h为射弹侵入的距离,A为阻力面积,Cd为阻力系数,ρ为水的密度。

    射弹沿 x 轴水平入水时,由于入水冲击过程的瞬时性,入水初期其运动方向基本保持不变,流体动力主要作用于射弹轴线方向,在射弹轴线方向上重力和浮力对射弹影响很小,基本可忽略不计。

    射弹入水冲击过程中,射弹浸没在水中的体积会排挤液面流体产生隆起现象,如图1所示,有效液面决定了自由水平液面的抬高程度,这取决于射弹锥头的外形和入水角等初始状态。沾湿因子定义为有效液面与实际液面的比值。

    图  1  射弹水平入水冲击示意图
    Figure  1.  Schematic diagram of horizontal impact of the projectile

    利用轴长体假设和Tayler关于在不同浸深时附加质量的表达式[13],参考垂直入水相关文献[10, 14],求出锥形射弹入水的附加质量为

    m=43ρ(Cwh)3tan3β
    (3)
    A=π(Cwh)2tan2β
    (4)
    dmdt=dmdhdhdt=4ρtan3β(Cwh)2Cwdhdt
    (5)

    式中:Cw为沾湿因子,β为射弹头部的半锥角,如图1所示。

    将式(3)~式(5)代入式(2),得到总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2=0
    (6)

    采用MATLAB软件,利用龙格-库塔方法进行求解,可以得到弹体入水时的冲击载荷。

    图2所示,采用全新的发射方式时,燃气排出并在膛口形成球形气体空腔。以地面为坐标系,枪口斜向上(即y轴正方向)、斜向下(即y轴负方向)发射时,根据圆切线定理,穿过球形气腔仍可看作垂直于液面入水,然而射弹入水在有攻角的情况下,轴向上会受到重力和浮力的作用分力影响,攻角正负值相反时,重力与浮力在射弹轴向上的作用分力方向也相反。

    图  2  射弹攻角示意图
    Figure  2.  Angle of attack of the projectile

    其他条件与水平入水时保持不变,攻角为正时,根据动量方程,得到动力学方程

    (M+m)d2hdt2+dmdtdhdt+12ρCdA(dhdt)2+MgsinαFbsinα=0
    (7)

    式中:α为射弹攻角,射弹斜向上发射时为正值,斜向下时为负值。

    射弹所受浮力为

    Fb=13πρg(Cwh)3tan2β
    (8)

    当攻角为正值时,得到的总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2+Mgsinα13πρg(Cwh)3tan2βsinα=0
    (9)

    同理,攻角为负值时的总方程为

    [M+43ρ(Cwh)3tan3β]d2hdt2+4ρtan3β(Cwh)2Cw(dhdt)2+12ρCdπ(Cwh)2tan2β(dhdt)2Mgsinα+13πρg(Cwh)3tan2βsinα=0
    (10)

    对式(9)、式(10)进行求解,可以得到不同攻角下的入水冲击载荷。

    由于水下射弹的质量较轻,入水速度较大,水下超空泡射弹入水平均速度约为600 m/s,因此入水速度对射弹入水冲击载荷的影响很大。本计算中,设射弹质量为0.14 kg,射弹头部半锥角为6°,计算得到不同速度时射弹的入水冲击载荷曲线,如图3所示,其中用射弹加速度反映入水冲击载荷。可以看出:射弹的入水冲击载荷先增大后减小,载荷变化主要发生在射弹入水前1 ms,最后渐渐趋于稳定;入水速度越大,冲击载荷峰值越大,入水后达到峰值的时间越短。此外,计算了不同速度下的入水冲击载荷峰值,如图4所示。入水速度在400~700 m/s范围时,射弹入水冲击载荷峰值一般为103g量级(g为重力加速度)。从冲击载荷峰值与速度的关系可以得出入水冲击载荷的峰值与速度基本呈线性关系。

    图  3  不同速度时的入水冲击载荷
    Figure  3.  Driving impact loads at different velocities
    图  4  不同速度时的冲击载荷峰值
    Figure  4.  Peak impact loads at different velocities

    设射弹质量为0.14 kg,通过计算获得了不同锥角的锥形射弹以600 m/s入水时的冲击载荷曲线,如图5所示。可以看出:射弹的入水冲击载荷先增大后减小,最后趋于稳定;半锥角β越大,冲击载荷峰值越大,并且入水后达到峰值的时间越短。改变半锥角,计算出不同半锥角情况下射弹入水冲击载荷峰值,如图6所示。从冲击载荷峰值与半锥角的关系可以得出入水冲击载荷峰值与半锥角基本呈线性关系。

    图  5  不同锥角时的入水冲击载荷
    Figure  5.  Driving impact loads at different cone angles
    图  6  不同锥角时的冲击载荷峰值
    Figure  6.  Peak impact loads at different cone angles

    当射弹质量为0.14 kg,入水速度为600 m/s,攻角α分别取45°、−45 °和0°(即水平入水)时,计算得到的射弹入水冲击载荷曲线如图7所示。可见,3条曲线基本重叠,差值在1g量级,相比于速度和半锥角对冲击载荷的影响,攻角对冲击载荷的影响几乎可以忽略不计。

    图  7  不同攻角时的入水冲击载荷
    Figure  7.  Inlet impact loads at different angles of attack

    此外,射弹在高速入水状态下,从液面进入水中的实际深度与理想深度存在一定的偏差,为了表达该误差,引入沾湿因子,沾湿因子的取值对冲击载荷的结果也存在影响。沾湿因子不同时射弹的入水冲击载荷如图8所示。从图8可以看出:沾湿因子越大,入水冲击载荷峰值越大;但沾湿因子对入水冲击载荷的影响较小,当沾湿因子变化值为0.1时,入水冲击载荷的变化在10%以内。

    图  8  沾湿因子不同时的入水冲击载荷
    Figure  8.  Inlet impact loads with different wet factors

    (1)射弹头部锥角相同时,入水冲击载荷峰值与速度呈正线性相关,入水速度越大,冲击载荷达到峰值的时间越短;射弹入水速度相同时,入水冲击载荷的大小与锥角呈正线性相关,锥角越大,冲击载荷达到峰值的时间越短。

    (2)射弹锥角和入水速度相同、入水攻角不同时,入水冲击载荷曲线与水平入水曲线基本重合,说明射弹重力和浮力在轴向上对入水冲击载荷的影响相对入水阻力几乎可以忽略不计。当射弹带有攻角入水后,重力和浮力更多的是对射弹产生径向力矩影响。

    (3)沾湿因子越大,入水冲击载荷峰值越大,但对入水冲击载荷影响较小,当沾湿因子变化为0.1时,入水冲击载荷的变化在10%以内。沾湿因子作为一个变量,其大小反过来也取决于入水冲击载荷,入水冲击载荷越大,沾湿因子也越大。

    (4)理论模型借鉴了高速弹体垂直入水的理论模型,数值模拟计算结果与已报道的高速射弹垂直入水冲击载荷数值计算和仿真结果高度一致,验证了本数学建模和数值计算的准确性。

    (5)探讨了全水下发射高速射弹入水瞬间的冲击载荷,高速射弹入水后形成超空泡,冲击载荷迅速减小。本工作对射弹未形成超空泡的情况进行了模型推导和数值模拟计算,对形成超空泡之前的理论研究具有重要意义。

  • 图  处理前黑钨矿石试样

    Figure  1.  Wolframite ore sample before treatment

    图  黑钨矿石内部钨颗粒岩相

    Figure  2.  Tungsten granular facies inside wolframite

    图  石英和硅质岩的断裂界面

    Figure  3.  Quartz and siliceous rock interface fracture map

    图  原位加载实验装置

    Figure  4.  In-situ loading experimental device

    图  自制原位加载压头示意图(单位:mm)

    Figure  5.  Sketch map of self-made in-situ loading indenter (Unit: mm)

    图  原位加载实验示意图

    Figure  6.  Schematic diagram of in-situ loading experiment

    图  不同加载速率下石英-硅质岩组分界面的应力-应变曲线

    Figure  7.  Stress-strain curves of quartz-siliceous rock component interface under different loading rates

    图  不同加载速率下石英-钨组分界面的应力-应变曲线

    Figure  8.  Stress-strain curves of quartz-tungsten component interface under different loading rates

    图  三维矿石RVE模型和网格模型

    Figure  9.  3D ore RVE model and mesh model

    图  10  石英-硅质岩界面原位加载模拟结点的应力、应变分布

    Figure  10.  Stress and strain distributions of in-situ loading simulation node of the quartz-siliceous rock interface

    图  11  石英-钨界面原位加载模拟结点的应力、应变分布

    Figure  11.  Stress and strain distributions of in-situ loading simulation node of the quartz-tungsten interface

    表  1  钨颗粒的特征尺寸

    Table  1.   Characteristic sizes of different tungsten particles

    W particle No.Length/mmWidth/mmPerimeter/mmArea/mm2
    11.6091.5096.9901.2295
    22.9662.14410.980 3.6125
    32.6452.23313.930 2.8330
    41.8491.1867.6342.3356
    50.9220.1372.2300.0703
    61.1980.7754.7200.4160
    71.0470.6223.4420.6056
    81.0400.6743.4360.6587
    下载: 导出CSV

    表  2  矿物界面的最小破碎应力

    Table  2.   Minimum fracture stress of component mineral interface

    Interface typesMinimum fracture stress/GPa
    12345678
    Quartz-siliceous rock interface1.40341.17851.33551.27321.28451.31711.47841.4826
    Quartz-tungsten interface1.54201.43731.44581.46711.36381.33551.43171.4571
    下载: 导出CSV

    表  3  不同黑钨矿石矿物组分的材料属性

    Table  3.   Mineral material properties of different wolframite components

    MineralDensity/(g·cm−3)Poisson’s ratioElastic modulus/GPaConstitutive modelStructural symmetry
    Quartz2.650.134ElastoplasticIsotropic
    Siliceous rock2.720.235ElastoplasticIsotropic
    Tungsten7.250.2035 ElastoplasticIsotropic
    下载: 导出CSV

    表  4  内含物的形态参数

    Table  4.   Inclusion morphological parameters

    MineralShapeVolume fractionAspect ratioDiameter/cm
    Siliceous rockCylinder0.502
    TungstenEllipsoid0.0510.2
    下载: 导出CSV

    表  5  石英-硅质岩界面结点的应力和应变

    Table  5.   Quartz-siliceous rock interface junction stress and strain

    Area of loadingNode No.Stress/GPaStrain/10−4 Area of loadingNode No.Stress/GPaStrain/10−4
    Boundary115981.398791.32 Internal607521.546680.12
    116051.405421.28607531.613140.14
    116101.408671.30*61394 1.589470.26
    116121.155261.06*61396 1.217560.09
    116131.697451.53*61398 1.365770.09
    607461.457941.48*61400 1.621670.30
    607471.639801.53External*61390 0.249880.42
    607481.655261.61*61404 0.230540.38
    607491.612491.58115970.420940.63
    607541.537911.42116010.243540.63
    *61392 1.330851.35116040.321360.47
    *61402 1.194141.29607450.379220.62
    Internal115991.387960.47607500.526460.87
    116081.616470.54607550.293990.55
    116091.300320.25607590.367650.69
    607511.612760.26615200.412910.72
    下载: 导出CSV

    表  6  石英-钨界面结点的应力和应变

    Table  6.   Quartz-tungsten interface junction stress and strain

    Area of loadingNode No.Stress/GPaStrain/10−4 Area of loadingNode No.Stress/GPaStrain/10−4
    Boundary363671.548420.21 Internal608301.961830.54
    363691.583280.24*613491.426061.89
    363701.851870.28*613511.614531.84
    608131.728781.59*613531.267901.46
    608141.983272.18External608100.316630.54
    608151.517711.55608120.370340.64
    608281.916921.96608260.410690.71
    608312.006892.53608290.504120.85
    608362.001152.31613590.586311.02
    *613470.969731.13613610.524180.85
    *613551.482551.90*613330.575040.72
    Internal363661.594340.09*613430.513550.57
    608161.722960.46*613450.677650.95
    608172.143090.82*613571.135981.54
    下载: 导出CSV
  • [1] 张成良, 刘磊, 王超. 高等岩石力学及工程应用 [M]. 长沙: 中南大学出版社, 2016.

    ZHANG C L, LIU L, WANG C. Advanced rock mechanics and engineering application [M]. Changsha: Central South University Press, 2016.
    [2] GAUDIN A M. Principles of mineral dressing [M]. New York: McGraw-Hill, 1939: 70−89.
    [3] 王雅蓉, 周乐光. 相界特征对矿物单体解离度的影响 [J]. 东北大学学报(自然科学版), 1996, 17(3): 88–92.

    WANG Y R, ZHOU L G. Effects of interfacial characteristics on mineral liberation [J]. Journal of Northeastern University (Natural Science), 1996, 17(3): 88–92.
    [4] 刘建远, 应平. 颗粒床压载粉碎对某硫化铜矿石矿物解离的影响 [J]. 有色金属(选矿部分), 2018(5): 81–87. doi: 10.3969/j.issn.1671-9492.2018.05.016

    LIU J Y, YING P. Impact of comminution by particle-bed compressive stressing on mineral liberation of a copper sulfide ore [J]. Nonferrous Metals (Mineral Processing Section), 2018(5): 81–87. doi: 10.3969/j.issn.1671-9492.2018.05.016
    [5] SUN W J, WANG L B, WANG Y Q. Mechanical properties of rock materials with related to mineralogical characteristics and grain size through experimental investigation: a comprehensive review [J]. Frontiers of Structural and Civil Engineering, 2017, 11(3): 322–328. doi: 10.1007/s11709-017-0387-9
    [6] 杨小彬, 周杰, 宋义敏, 等. 循环加载岩石界面滑移位移演化特征试验研究 [J]. 煤炭学报, 2019, 44(10): 3041–3048. doi: 10.13225/j.cnki.jccs.2018.1579

    YANG X B, ZHOU J, SONG Y M, et al. Evolution characteristics of sliding displacement of rock interface under cyclic loading [J]. Journal of China Coal Society, 2019, 44(10): 3041–3048. doi: 10.13225/j.cnki.jccs.2018.1579
    [7] ALNEASAN M, BEHNIA M, BAGHERPOUR R. Analytical investigations of interface crack growth between two dissimilar rock layers under compression and tension [J]. Engineering Geology, 2019, 259: 105188. doi: 10.1016/j.enggeo.2019.105188
    [8] 蔡改贫, 赵小涛. 基于细观力学的矿石颗粒破碎特性研究 [J]. 应用力学学报, 2020, 37(4): 1792–1797. doi: 10.11776/cjam.37.04.B097

    CAI G P, ZHAO X T. Study on the fracture characteristics of ore particles based on micromechanics [J]. Chinese Journal of Applied Mechanics, 2020, 37(4): 1792–1797. doi: 10.11776/cjam.37.04.B097
    [9] 王光勇, 余锐, 马东方, 等. 饱水细砂岩动态抗拉与抗压强度试验对比研究 [J]. 高压物理学报, 2020, 34(4): 044101. doi: 10.11858/gywlxb.20190857

    WANG G Y, YU R, MA D F, et al. Comparative study on dynamic tensile and compressive strength of the saturated fine sandstone [J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 044101. doi: 10.11858/gywlxb.20190857
    [10] 纪杰杰, 李洪涛, 吴发名, 等. 冲击荷载作用下岩石破碎分形特征 [J]. 振动与冲击, 2020, 39(13): 176–183, 214. doi: 10.13465/j.cnki.jvs.2020.13.026

    JI J J, LI H T, WU F M, et al. Fractal characteristics of rock fragmentation under impact load [J]. Journal of Vibration and Shock, 2020, 39(13): 176–183, 214. doi: 10.13465/j.cnki.jvs.2020.13.026
    [11] 蔡改贫, 宣律伟, 张雪涛, 等. 多尺度内聚颗粒模型破碎过程研究 [J]. 岩土力学, 2020, 41(6): 1809–1817. doi: 10.16285/j.rsm.2019.1218

    CAI G P, XUAN L W, ZHANG X T, et al. Investigation into the crushing process in multi-scale cohesive particle model [J]. Rock and Soil Mechanics, 2020, 41(6): 1809–1817. doi: 10.16285/j.rsm.2019.1218
  • 加载中
图(11) / 表(6)
计量
  • 文章访问数:  1158
  • HTML全文浏览量:  439
  • PDF下载量:  27
出版历程
  • 收稿日期:  2021-11-03
  • 修回日期:  2021-11-26
  • 录用日期:  2021-12-03
  • 刊出日期:  2022-05-30

目录

/

返回文章
返回