
Citation: | YANG Shuqi, ZHANG Xu, PENG Wenyang, SHU Junxiang, QIN Shuang, ZHONG Bin. Impact Initiation Characteristics of TATB Based Insensitive Explosives Mixed with HMX by Electromagnetic Velocity Gauges[J]. Chinese Journal of High Pressure Physics, 2020, 34(3): 033403. doi: 10.11858/gywlxb.20190852 |
球体的冲击破碎是颗粒物质力学的基本问题。球体结构材料广泛存在于日常生活和工程应用中,例如煤粒破碎[1-2]、落石灾害[3]、化工产品[4]等。研究球体材料在冲击载荷下的破碎,对于提高颗粒物质的加工技术、节约能源、促进经济发展、有效监控及预防自然灾害等具有重要价值。
早期国内外研究人员对单颗粒的冲击破碎开展了大量的实验研究,研究的重点是颗粒的临界破碎速度、颗粒的破碎模式、碎片形貌以及碎块的统计特征等。Andrews等[5-6]对钠钙玻璃球和陶瓷球进行了研究,在临界速度下,球体从刚性靶板回弹;超过临界速度后,钠钙玻璃球碎裂为大量尺寸不均的细小碎片,而陶瓷球则碎裂为若干大小均匀形似“月牙状”的碎块,并将钠钙玻璃球体破碎归因于球体压缩端的剪切应力,而将陶瓷球体碎裂归因于球体中部的冲击拉伸应力。Salman等[7-9]对不同类型材料和大小的脆性球体进行试验,结果表明大部分脆性材料的破坏模式是相似的,即:在中低速下,球体的局部裂纹扩展导致局部破碎,在撞击端形成压缩平台;在高速下,球体碎裂为大量难以识别的粉状碎片。Chau[10]、Wu[11]等根据冲击动能的转换分析了脆性球体破碎后碎块的统计特征。
在数值计算方面,易洪昇等[12]使用ABAQUS软件模拟了冲击过程中的球体内部应力情况,由于颗粒破碎涉及大变形,采用有限元方法会出现网格畸变问题,因此不适合模拟大变形和损伤现象。离散元方法(Discrete element method, DEM)能够满足岩体工程中破裂和裂纹发展等问题的研究,适合于从本质上研究固体介质的力学特性[13]。Shen等[14]利用离散元软件分析了球形岩石在不同应变率情况下碎块特征、碎块数量、损伤率等的变化情况。Carmona等[15]采用离散元方法研究了球体碎裂过程中裂纹的扩展过程,发现碎裂的大块碎片质量符合两参数的Weibull分布。
以石英玻璃球作为典型的脆性球体,开展了较广泛冲击速度下石英玻璃球撞击刚性壁的破坏现象研究,利用离散元数值软件PFC3D分析了石英玻璃球在不同撞击速度下的撞击碎裂过程。
实验采用高速气枪作为石英玻璃球的加载发射装置,撞击区域使用碳化钨材料作为近似刚性靶板,碳化钨的弹性模量为720 GPa,密度为15.6 g/cm3,泊松比为0.21。碳化钨周围附有厚壁钢块,如图1所示。采用Kirana超高速摄像机记录球体的破碎过程,拍摄幅频为2×105帧每秒,高速相机采用断路方式触发,触发线设置在气枪的出口位置。在球体撞击区域的外围设计了一个有机玻璃防护装置用于收集碎片。
高速气枪出口距离靶板250 mm,通过圆柱形长杆调节高速气枪和靶板的垂直度,从而实现玻璃球的正撞破坏。通过水平放置和位于球体上方放置的高速相机,实验观察石英玻璃球的飞行姿态,结果表明玻璃球体未发生明显偏转,撞击靶板位置为初始预设位置。根据超高速相机拍摄可获得玻璃球体的飞行时间,从而确定玻璃球体飞行的平均速度,近似等效为球体撞击速度。
实验试件采用昊天玻璃制品有限公司生产的超纯石英玻璃球,其SiO2的质量分数大于99.99%。球体直径为(10.00±0.02)mm,质量为(1.150±0.015)g,具有良好的透光性,且表面无明显缺陷。
为了研究玻璃球的临界冲击破碎速度,在低速下开展了大量的实验。实验结果表明:当石英玻璃球以11.87 m/s撞击刚性靶板时,碰撞为近似弹性碰撞,球体以略低于原速(回弹速度10.25 m/s)从刚性壁回弹,除却撞击点应力集中破坏外,球体内部没有产生宏观可见裂纹。石英玻璃球撞击前和回弹后如图2所示。
图3是冲击速度为78 m/s的冲击破碎过程。从图3可以看出,球体的破碎可以分为3个阶段。在初始阶段,撞击端首先发生破坏,形成经过子午面的宏观裂纹,同时形成细条状碎片,这些碎片以撞击点为中心向四周扩散,而在远离撞击端部分还保持相对完整;随着撞击的持续进行,在撞击端的碎片继续向四周运动,在远离撞击端的部分产生“龟甲状”的裂纹,并最终碎裂为若干“圆锥状”的碎块。
当冲击速度为35 m/s时,球体从刚性壁回弹,在撞击端产生了一些细小的碎片,如图4(a)所示。图5(a)表明,实验回收后的碎片存在3个破坏区域,分别为压缩破碎区、表面剥落区以及剪切破坏区。其中压缩破碎区是一个圆锥状的区域,这块区域的石英玻璃已不具备透光性,说明这一区域的石英玻璃发生了剧烈的压缩破碎,但是这些碎块还有部分镶嵌在球体上并没有从球体上脱落。在压缩破碎区的周围是一个环形的表面剥落区,受到表面波作用一些片状的碎片从球体上剥落,产生层片状的碎片。在球体内部,石英玻璃的折光性发生了改变,说明在剪切力的作用下石英玻璃内部产生剪切裂纹带,剪切裂纹带由压缩破碎区开始起裂,沿着子午面由撞击端向试件内部扩展。
当冲击速度为50 m/s时,从球体的透光性可以看出,球体完全碎裂为若干碎块,如图4(b)所示。在撞击端同样会出现类似于35 m/s时的压缩破碎现象,并且涉及的范围更大。由于冲击动能的增加,剪切破碎区完全贯穿了球体的子午面,最终导致石英玻璃球碎裂成若干类似于“月牙状”的碎块,如图5(b)所示。
当冲击速度为78 m/s时,与上述几种破碎情况完全不同,石英玻璃球不再从刚性壁回弹,而是发生了“坍塌式”的破碎,如图4(c)所示。在撞击端,产生大量细小的碎块,在球体的另一端表面有许多类似于“经线形”和“纬线形”的宏观裂纹,最终产生大量类似于“龟甲状”的碎片,内部也产生大量细小碎片,如图5(c)所示。
当冲击速度为135 m/s时,发生“蘑菇云”似的破碎现象,如图4(d)所示。产生的碎片更加细小,在远离撞击端的地方发生明显的“层裂”现象。球体在高速撞击过程中碎裂为极细小的碎片,透明球体不易观察裂纹发展,将石英玻璃球表面涂黑,有利于观测球体表面裂纹扩展。由图5(d)可知,随着速度的增加,球体内部破碎更为剧烈,产生大量粉末状细小碎片。
综上所述,动态冲击下石英玻璃球的破坏形式主要有4种特征:(1)局部破碎,在球体上形成3个破碎区域;(2)中心劈裂,沿着冲击轴线碎裂为类似于“月牙状”的碎片;(3)坍塌破坏,石英玻璃球整体破碎,碎裂为若干碎块;(4)层裂破坏,在远离撞击端的地方产生“蘑菇云”似的“层裂”破坏。
选择目前离散元中使用广泛的平行黏结模型(Linear parallel bond mode)。平行黏结模型本质上是用一层线性的黏结体将相互重叠的颗粒绑定起来,这些黏结的颗粒能很好地模拟岩石、玻璃等脆性材料的压缩破坏过程[16]。
为了使模型具有较高的密实度,将模型中的孔隙率设定为0.2,最大颗粒与最小颗粒的粒径比为3∶2。颗粒间的摩擦系数对材料的宏观参数影响不大[17],故取经验值0.577。
根据文献[18-19]的标定方法,通过单轴压缩标定了模型的弹性模量、泊松比和抗压强度,通过单轴拉伸标定抗拉强度,由三点弯标定弯曲强度,带缺口的三点弯标定断裂韧性和颗粒半径。最终得到的主要微观参数如表1所示。标定的宏观参数如表2所示,得到石英玻璃的宏观参数与制造商提供的材料参数基本一致,因此可以用于模拟石英玻璃的冲击破碎。
Effective modulus of linear contact/GPa | Normal to shear stiffness ratio of linear contact | Porosity | Minimum radius of particles/mm | Size ratio of maximum and minimum particles | Tensile strength of contact/MPa | Shear strength of contact/MPa |
55 | 2.9 | 0.2 | 0.1 | 1.5 | 300 | 600 |
Method | Equivalent density/(kg·m–3) | Elastic modulus/GPa | Poisson’s ratio | Compressive strength/MPa | Tensile strength/MPa | Bending strength/MPa | Fracture toughness/(N·m–3/2) |
Manufacture provide | 2.203 | 77.8 | 0.170 | 860 | 50 | 67.0 | 0.78 |
DEM numerical simulation | 2.203 | 78.0 | 0.172 | 798 | 50 | 67.4 | 0.85 |
根据上述标定的离散元微观参数建立与实验同尺寸的几何模型。采用刚性墙面模拟实验中的刚性壁,在模拟初始时刻,给所有颗粒施加垂直撞击刚性壁的速度,使其匀速撞击刚性壁。
图6是球体以78 m/s撞击刚性壁的过程。模拟中不同颜色代表的是碰撞过程中产生的不同碎块。对比图3和图6可知,模拟与实验吻合较好,模拟中玻璃球的破坏过程存在明显的时序性。在初始时刻,在撞击端发生局部的破碎,产生细小的碎片,同时裂纹向远离撞击端的地方传播;随着撞击的进行,撞击端产生的碎片向外运动,远离撞击端的地方碎裂为较大的锥形碎块,当撞击端处的细小碎片运动到更远处时,这些较大的锥形碎块以一定的速度继续撞击刚性壁。
球体以78 m/s撞击刚性壁的过程中,球体速度(左端粒子平均速度)、球体内部裂纹数量(粒子之间的力键)以及球体撞击力的时程曲线如图7所示。石英玻璃球撞击破碎过程大致可分为3个阶段:第1阶段,可近似为弹性压缩过程,球体所受到的冲击载荷也近似线性增加,大约经历
图8为在不同速度下石英玻璃球速度随时间变化的曲线。图8(a)显示了撞击速度低于临界破碎速度时球体速度的变化过程。可以看出,球体撞击靶板后匀速下降,最终几乎以原速从靶板回弹。图8(b)显示了撞击速度高于临界破碎速度时球体速度的变化过程。可以看出,球体左端粒子在撞击刚性壁后平均速度急剧下降并保留部分残余速度,撞击速度越大,速度下降越快,残余速度也越大。球体左端粒子以残余速度自由飞行,石英玻璃球撞击端(右端粒子)完全破碎,左端粒子在向右运动过程中没有任何阻力,将球体匀速运动部分对时间积分,发现球体的运动距离均约为球体直径的1/2。
图9给出了0~300 m/s撞击速度下石英玻璃球撞击力随撞击速度的变化规律。从图9可以看出,撞击速度为20 m/s左右时为撞击力增长规律的临界转折点,此前球体未产生明显破碎,此后球体破碎使得撞击力上涨趋势大幅衰减。
在临界转折速度前,石英玻璃球撞击刚性靶板过程可近似认为是弹性碰撞,通过Hertz接触理论可得弹性碰撞接触力[20]
Pm=(53πρ)35(34KI)−25V65R2 |
(1) |
式中:
图10给出了临界转折速度附近石英玻璃球撞击力随撞击速度的变化规律。从图10可以看出:当速度低于临界转折速度时,球体不破碎,内部几乎不产生裂纹,球体整体近似弹性碰撞,Hertz接触理论可以较好地预测石英玻璃球的撞击力;而速度高于临界转折速度时,球体发生破碎,不再具有球体整体性,破碎过程将卸载撞击力,并消耗部分能量,飞散的碎片也将带走大量动能;随着冲击速度的增加,石英玻璃球的极限撞击将等效于流体撞击。
冲击速度在临界破坏速度下为弹性碰撞,球体以略低于原速从刚性壁回弹,除却撞击点应力集中破坏外,球体内部没有产生宏观可见裂纹。低速碰撞时,球体呈现“压缩破碎区-表面剥落区-剪切破坏区”的破坏结构,剪切裂纹沿着子午面由撞击端向试件内部扩展,并还不足以形成碎片;中速碰撞时,由于剪切破坏区充分扩展,致使球体碎裂为若干“月牙状”的碎块;高速碰撞时,石英玻璃球发生坍塌式破碎,表面产生大量类似于“龟甲状”的碎片,内部也产生大量细小碎片;当速度更高时,在远离撞击端有明显的层裂现象产生,球体破碎产生大量粉末状碎片。
离散元软件PFC3D再现了石英玻璃球撞击刚性壁的过程,球体在高速碰撞下破碎可以分为弹性压缩、整体破碎和二次撞击3个阶段。石英玻璃球撞击力随撞击速度呈两段式增长,在球体碎裂前,撞击过程近似为弹性碰撞过程,Hertz接触理论可以较好地描述其撞击力,在球体碎裂后,破碎过程将卸载撞击力,并消耗部分能量,飞散的碎片也将带走大量动能,撞击力远低于Hertz接触理论值,并且随着撞击速度增大,偏差越大。
[1] |
张琪敏, 张旭, 赵康, 等. TATB基钝感炸药JB-9014的冲击起爆反应增长规律 [J]. 爆炸与冲击, 2019, 39(4): 041405.
ZHANG Q M, ZHANG X, ZHAO K, et al. Law of reaction growth of shock initiation on the TATB based insensitive explosive JB-9014 [J]. Explosion and Shock Waves, 2019, 39(4): 041405.
|
[2] |
HILL L G, GUSTAVSEN R L, ALCON R R, et al. Shock initiation of new and aged PBX 9501 measured with embedded electromagnetic particle velocity gauges: LA-13634-MS [R]. New Mexico, US: Los Alamos National Laboratory, 1999.
|
[3] |
ZHANG X, WANG Y F, HUANG W B, et al. Reaction buildup of PBX explosives JOB-9003 under different initiation pressures [J]. Journal of Energetic Materials, 2017, 35(2): 197–212. doi: 10.1080/07370652.2016.1250841
|
[4] |
GUSTAVSEN R L, GEHR R J, BUCHOLTZ S M, et al. Shock initiation of the tri-amino-tri-nitro-benzene explosive PBX9502 cooled to -55 ℃ [J]. Journal of Applied Physics, 2012, 112(7): 074909. doi: 10.1063/1.4757599
|
[5] |
张涛, 赵继波, 伍星, 等. 未反应JBO-9021炸药冲击雨贡纽曲线的研究 [J]. 高压物理学报, 2016, 30(6): 457–462. doi: 10.11858/gywlxb.2016.06.004
ZHANG T, ZHAO J B, WU X, et al. Hugoniot curve of unreacted JBO-9021 explosive [J]. Chinese Journal of High Pressure Physics, 2016, 30(6): 457–462. doi: 10.11858/gywlxb.2016.06.004
|
[6] |
GUSTAVSEN R L, SHEFFIELD S A, ALCON R R. Measurements of shock initiation in the tri-amino-tri-nitro-benzene based explosive PBX 9502: wave forms from embedded gauges and comparison of four different material lots [J]. Journal of Applied Physics, 2006, 99(11): 114907. doi: 10.1063/1.2195191
|
[7] |
张涛, 谷岩, 赵继波, 等. JBO-9021炸药的化学反应区宽度 [J]. 爆炸与冲击, 2017, 37(3): 415–421. doi: 10.11883/1001-1455(2017)03-0415-07
ZHANG T, GU Y, ZHAO J B, et al. Chemical reaction zone length of JBO-9021 [J]. Explosion and Shock Waves, 2017, 37(3): 415–421. doi: 10.11883/1001-1455(2017)03-0415-07
|
[8] |
SHEFFIELD S A, BLOOMQUIST D D, TARVER C M. Subnanosecond measurements of detonation fronts in solid high explosives [J]. The Journal of Chemical Physics, 1984, 80(8): 3831–3844. doi: 10.1063/1.447164
|
[9] |
裴红波, 黄文斌, 覃锦程, 等. 基于多普勒测速技术的JB-9014炸药反应区结构研究 [J]. 爆炸与冲击, 2018, 38(3): 485–490.
PEI H B, HUANG W B, QIN J C, et al. Reaction zone structure of JB-9014 explosive measured by PDV [J]. Explosion and Shock Waves, 2018, 38(3): 485–490.
|
[10] |
赵同虎, 张新彦, 李斌, 等. 用光电法研究钝感炸药JB-9014反应区结构 [J]. 高压物理学报, 2002, 16(2): 111–119. doi: 10.3969/j.issn.1000-5773.2002.02.005
ZHAO T H, ZHANG X Y, LI B, et al. Detonation reaction zones tructure of JB-9014 [J]. Chinese Journal of High Pressure Physics, 2002, 16(2): 111–119. doi: 10.3969/j.issn.1000-5773.2002.02.005
|
[11] |
LOBOIKO B G, LUBYATINSKY S N. Reaction zones of detonating solid explosives [J]. Combustion, Explosion, and Shock Waves, 2000, 36(6): 716–733. doi: 10.1023/A:1002898505288
|
[1] | HUANG Min, ZHU Benhao, XIAO Gesheng, QIAO Li. Simulation on Deformation Damage and Strain Rate Effect of Nb3Sn Composite Superconductors under Cycling Load at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2024, 38(2): 024201. doi: 10.11858/gywlxb.20230755 |
[2] | YANG Dong, JIANG Ziwei, ZHENG Zhijun. Dynamic Behavior and Constitutive Relationship of Titanium Alloy Ti6Al4V under High Temperature and High Strain Rate[J]. Chinese Journal of High Pressure Physics, 2024, 38(1): 014101. doi: 10.11858/gywlxb.20230743 |
[3] | YAN Jitao, CHEN Wenfei, ZHANG Hao, YOU Jian, SHI Chaoming, JIANG Haocheng, ZHU Jue. Dynamic Mechanical Behavior of G550 Cold-Formed Steel under High Temperature and High Strain Rate[J]. Chinese Journal of High Pressure Physics, 2023, 37(2): 024101. doi: 10.11858/gywlxb.20220705 |
[4] | QIAO Jingyan, LI Jinzhu, ZHANG Xihuang, YAO Zhiyan, SHEN Shiliang. Dynamic and Static Tensile Mechanical Properties of Glass Fiber Reinforced Plastics[J]. Chinese Journal of High Pressure Physics, 2023, 37(3): 034102. doi: 10.11858/gywlxb.20230618 |
[5] | WEI Jiawei, SHI Xiaopeng, FENG Zhenyu. Strain Rate Dependent Constitutive Model of Rubber[J]. Chinese Journal of High Pressure Physics, 2022, 36(2): 024205. doi: 10.11858/gywlxb.20210815 |
[6] | YU Wenfeng, LI Jinzhu, YAO Zhiyan, HUANG Fenglei. Mechanical Behaviors and Constitutive Model of Polymide under Quasi-Static and Dynamic Compressive Loading[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 044101. doi: 10.11858/gywlxb.20210922 |
[7] | TANG Ruitao, XU Liuyun, WEN Heming, WANG Zihao. A Macroscopic Dynamic Constitutive Model for Ceramic Materials[J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 044201. doi: 10.11858/gywlxb.20190863 |
[8] | WU Qijian, ZHI Xudong. Strain Rate Effect of GFRP-Reinforced Circular Steel Tube under Low-Velocity Impact[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 044203. doi: 10.11858/gywlxb.20180653 |
[9] | LI Xueyan, LI Zhibin, ZHANG Duo. Constitutive Model of Aluminum Foams Considering Temperature Effect under Quasi-Static Compression[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 044103. doi: 10.11858/gywlxb.20170642 |
[10] | CHEN Qing-Shan, MIAO Ying-Gang, GUO Ya-Zhou, LI Yu-Long. Comparative Analysis of 3 Constitutive Models for 93 Tungsten Alloy[J]. Chinese Journal of High Pressure Physics, 2017, 31(6): 753-760. doi: 10.11858/gywlxb.2017.06.010 |
[11] | GAO Guang-Fa. Effect of Strain-Rate Hardening on Dynamic Compressive Strength of Plain Concrete[J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 261-270. doi: 10.11858/gywlxb.2017.03.007 |
[12] | GAO Guang-Fa. Hardening Effect of the Strain Rate on the Dynamic Tensile Strength of the Plain Concrete[J]. Chinese Journal of High Pressure Physics, 2017, 31(5): 593-602. doi: 10.11858/gywlxb.2017.05.013 |
[13] | LI Ding-Yuan, ZHU Zhi-Wu, LU Ye-Sen. Mechanical Properties and Constitutive Relation for 42CrMo Steel under Impact Load[J]. Chinese Journal of High Pressure Physics, 2017, 31(6): 761-768. doi: 10.11858/gywlxb.2017.06.011 |
[14] | WANG Li-Wen, PANG Bao-Jun, CHEN Yong, ZHANG Kai. Study on Dynamic Mechanical Behavior and Constitutive Model of Reactive Powder Concrete after Exposure in High Temperature[J]. Chinese Journal of High Pressure Physics, 2012, 26(4): 361-368. doi: 10.11858/gywlxb.2012.04.001 |
[15] | LI Jun-Ling, LU Fang-Yun, FU Hua, ZHAO Yu-Gang, TAN Duo-Wang. Research on the Dynamic Behavior of a PBX Explosive[J]. Chinese Journal of High Pressure Physics, 2011, 25(2): 159-164 . doi: 10.11858/gywlxb.2011.02.012 |
[16] | PANG Bao-Jun, YANG Zhen-Qi, WANG Li-Wen, CHI Run-Qiang. Dynamic Compression Properties and Constitutive Model with Strain Rate Effect of Rubber Material[J]. Chinese Journal of High Pressure Physics, 2011, 25(5): 407-415 . doi: 10.11858/gywlxb.2011.05.005 |
[17] | XIAO Da-Wu, LI Ying-Lei, HU Shi-Sheng. Constitutive Model of Pure Zirconium under High Temperature and High Strain Rate[J]. Chinese Journal of High Pressure Physics, 2009, 23(1): 46-50 . doi: 10.11858/gywlxb.2009.01.008 |
[18] | LIU Yao-Dong, YU Ji-Lin, ZHENG Zhi-Jun. Effect of Inertia on the Dynamic Behavior of Cellular Metal[J]. Chinese Journal of High Pressure Physics, 2008, 22(2): 118-124 . doi: 10.11858/gywlxb.2008.02.002 |
[19] | ZHANG Peng, ZHAO Feng, BAI Shu-Lin, FU Bin, ZHANG Guo-Hua. Experimental Study on the Microstructures and Dynamic Behavior of a PBX Substitute Material[J]. Chinese Journal of High Pressure Physics, 2007, 21(1): 20-28 . doi: 10.11858/gywlxb.2007.01.004 |
[20] | WANG Zheng, NI Yu-Shan, CAO Ju-Zhen, ZHANG Wen. Building of a Constitutive Model for Concrete under Dynamic Impact[J]. Chinese Journal of High Pressure Physics, 2006, 20(4): 337-344 . doi: 10.11858/gywlxb.2006.04.001 |
Effective modulus of linear contact/GPa | Normal to shear stiffness ratio of linear contact | Porosity | Minimum radius of particles/mm | Size ratio of maximum and minimum particles | Tensile strength of contact/MPa | Shear strength of contact/MPa |
55 | 2.9 | 0.2 | 0.1 | 1.5 | 300 | 600 |
Method | Equivalent density/(kg·m–3) | Elastic modulus/GPa | Poisson’s ratio | Compressive strength/MPa | Tensile strength/MPa | Bending strength/MPa | Fracture toughness/(N·m–3/2) |
Manufacture provide | 2.203 | 77.8 | 0.170 | 860 | 50 | 67.0 | 0.78 |
DEM numerical simulation | 2.203 | 78.0 | 0.172 | 798 | 50 | 67.4 | 0.85 |
Effective modulus of linear contact/GPa | Normal to shear stiffness ratio of linear contact | Porosity | Minimum radius of particles/mm | Size ratio of maximum and minimum particles | Tensile strength of contact/MPa | Shear strength of contact/MPa |
55 | 2.9 | 0.2 | 0.1 | 1.5 | 300 | 600 |
Method | Equivalent density/(kg·m–3) | Elastic modulus/GPa | Poisson’s ratio | Compressive strength/MPa | Tensile strength/MPa | Bending strength/MPa | Fracture toughness/(N·m–3/2) |
Manufacture provide | 2.203 | 77.8 | 0.170 | 860 | 50 | 67.0 | 0.78 |
DEM numerical simulation | 2.203 | 78.0 | 0.172 | 798 | 50 | 67.4 | 0.85 |