Dynamic response and parameter analysis of concrete-filled steel tubular structure under lateral impact loading
-
摘要: 采用刚塑性结构模态分析法和数值模拟方法,对侧向冲击载荷作用下的圆截面钢管混凝土结构进行了塑性动力分析。将钢管混凝土等效为刚塑性地基梁模型,给出了钢管混凝土构件跨中侧向变形的模态解析解,得到了冲击载荷作用下影响结构最终侧向变形的无量纲参数。利用ABAQUS/Explicit软件建立了钢管混凝土结构在侧向冲击作用下的动态响应数值模型并进行了计算分析,将理论预测值和数值模拟结果与试验结果进行了交叉对比。结合量纲分析和数值模型对影响构件最终变形的几何、物理参数及初始冲量进行了分析。结果表明:理论预测值和数值模拟结果与试验结果吻合较好,结构的塑性变形与理论假定的塑性铰分布一致。构件几何参数中,长径比和径厚比对其侧向最终变形有较大影响;冲击头相对宽度可改变构件的变形模态;相比于几何参数,钢管和混凝土芯层的物理参数对构件跨中挠度的影响较小;结构的侧向变形与初始冲量成二次幂相关。最后给出了理论分析参数的适用范围。刚塑性响应模态解可较好地预测钢管混凝土结构在侧向冲击载荷作用下的塑性变形行为。Abstract: By employing the mode approximation method for rigid-plastic structural dynamic behavior and numerical simulation, a dynamic response analysis was conducted on circular-section concrete-filled steel tubular (CFST) structures subjected to lateral impact loadings. The mechanical model of the CFST structure was equivalently represented as a rigid-plastic foundation beam model according to its plastic behavior. Under the linear velocity field assumption and the geometric similarity, the equivalently initial velocity for mode approximation of the structure was derived and compared with the existing experimental data. An analytical solution for the plastic lateral deformation at the mid-span of the CFST with two fixed ends by the rigid-plastic mode approximation method was provided, yielding non-dimensional geometric and physical parameters that influenced the ultimate lateral plastic deformation. A numerical model of the CFST structure under lateral impact was established using ABAQUS/Explicit. The theoretical and numerical predictions were both compared with existing experimental global deformations. Dimensional analysis and numerical modeling were combined to analyze the geometric and physical parameters, as well as the initial impact impulse, which influence the plastic deformation of the CFST structure. The results demonstrate a good agreement between the theoretical, numerical results, and experimental data, confirming that the plastic deformations of the structure align with the assumed distribution of plastic hinges. For geometric variables, the ratio of length to diameter and ratio of thickness to diameter exert a significant influence on the final lateral deformation. The relative width of the indenter can alter the deformation shape of the structure. The physical parameters of the steel tube and core concrete have less impact on the deflection at the mid-span compared with the geometric variables. The final lateral deformation of the CFST structure exhibits a quadratic correlation with the initial impact impulse. Finally, the applicable range of all the theoretical analysis variables is given according to the corresponding parameter analysis. The proposed mode solutions for rigid-plastic response provide a reliable prediction of the plastic deformation behavior of the CFST structures under lateral impact loadings.
-
泡沫金属材料的率相关性能研究主要采用分离式霍普金森压杆(SHPB)动态压缩实验技术。由于泡沫铝复杂的微结构以及低阻抗、低波速、大变形等特点, 使其动态性能测试难度较大, 实验结果的离散性很大。同时, 由于泡沫铝材料制备方法多样, 泡沫铝的力学性能差异很大。迄今为止, 对泡沫铝应变率效应的研究尚没有得出一致的结论[1-2]。
V.S.Deshpande等[3]进行了0.001~5 000 s-1应变率范围内Alulight泡沫和Duocel泡沫的动、静态压缩实验, 得到闭孔泡沫铝和开孔泡沫铝均对应变率不敏感的结论。郭伟国等[4]对相对密度8.4%和4.5%的Duocel泡沫在0.001~500 s-1应变率下的形变和力学性能进行了系统的实验研究, 同样得出该材料对应变率不敏感的结论。K.A.Dannemann等[5]在Alporas泡沫和6101铝合金开孔泡沫(Duocel)的动、静态压缩实验中则发现, 400~2 500 s-1的应变率范围内闭孔泡沫铝具有明显的应变率效应, 开孔泡沫的应变率效应不明显。T.Mukai等[6-7]同样也得出了闭孔泡沫铝对应变率敏感的结论。程和法等[8-9]分别对相对密度0.363~0.419的开孔泡沫Al-Mg合金以及相对密度0.341~0.419的通孔泡沫铝, 在0.001~1 600 s-1和0.001~2 000 s-1应变率范围内的动、静态压缩力学性能进行了实验研究, 指出两种开孔泡沫材料均具有明显的应变率效应。田杰等[10]开展了4种不同基体材料泡沫铝的准静态和动态压缩实验, 得到的压缩应力应变曲线显示, 无论基体材料对应变率是否敏感性, 泡沫材料都表现出明显的应变率效应。张健等[11]对孔隙率60%~79%的闭孔泡沫铝在0.001~2 500 s-1应变率范围内的实验研究结果表明, 闭孔泡沫铝的应变率敏感性受孔隙率大小影响, 随着孔隙率的降低, 闭孔泡沫铝表现出明显的应变率效应, 高孔隙率闭孔泡沫铝对应变率基本不敏感。
尽管有很多泡沫铝对应变率敏感的结论, 但对于敏感性产生原因的解释不尽相同。K.A.Dannemann等[5]将其归结为泡孔内气体的黏性流动。胡时胜等[12]认为泡孔内气体的黏性流动引起的应力增加可忽略不计, 不可能是泡沫材料应变率敏感的根本原因, 敏感性主要是由泡孔结构的变形特性产生, 基体材料的应变率效应及泡孔的形状大小并不能对泡沫材料应变率的敏感性起主导作用。张健等[11]认为基体材料的应变率敏感性是泡沫材料应变率敏感性的决定性因素, 但影响程度受孔隙率大小影响, 随着孔隙率的增大敏感性降低; 微惯性、波效应和孔内气体压力对泡沫铝的平台应力不产生明显影响。
总体来说, 倾向于泡沫铝对应变率敏感的研究通常将敏感性产生原因归结为基体材料率相关性能、微结构惯性以及孔穴内部气体压力变化。而根据L.J.Gibson等[13]的研究结论, 可以认为, 孔穴气体压力变化对泡沫铝材料力学性能的影响很小, 不是引起泡沫铝材料应变率敏感的主要原因。此外, 尽管相对密度、应变率等参数不是泡沫铝材料是否存在应变率效应的决定因素, 但相对密度的变化可能会影响泡沫铝对应变率的敏感程度; 应变率(或加载速率)的变化可能会改变泡沫铝的变形机制, 进而影响泡沫铝材料应变率敏感性的产生原因。为此, 本文中, 研究基体材料率相关性能、微结构惯性对泡沫铝率相关性的影响, 分析相对密度与泡沫铝应变率敏感程度的关系, 以及应变率变化对泡沫铝动态压缩力学性能的影响。
1. 模型
1.1 三维随机球形泡孔模型
泡沫铝材料通常为无序开、闭孔混合结构。本文中采用Matlab和Abaqus软件生成了三维随机分布球形泡孔模型模拟泡沫铝材料的微细观结构, 模型几何边界呈周期性排列, 并采用代表性体积单元边长与泡孔孔径的比描述三维随机模型中泡孔的规模大小, 即:
δ=l/(2r) (1) 式中:δ为模型的泡孔规模大小, l为模型的边长, r为球形泡孔的半径。
图 1为泡孔规模δ=4, 孔隙率分别为65%、50%和35%的泡沫材料三维球形泡孔模型。模型采用四面体单元离散, 单元数量分别为218 032、306 139和377 064。
1.2 基体材料本构模型
采用Cowper-Symonds模型描述铝合金基体材料的应变率效应[11]:
σdy=σsy(1+(˙εC)1/P) (2) 式中:
为材料的动态屈服应力,
为材料的静态屈服应力,
为应变率, C、P为Cowper-Symonds应变率参数。
基体材料的弹塑性本构采用理想弹塑性模型。计算过程中, 应变率参数C和P分别为6 500 s-1和4, 杨氏模量为69 GPa, 泊松比为0.33, 静态屈服应力为250 MPa, 密度为2.7 g/mm3。另外, 当应变率参数C和P取值为零时, 表示基体材料率无关。
1.3 泡孔规模影响分析
泡孔模型在理论上应包含尽可能多的泡孔结构, 但过多的微结构将会导致计算量过大而使模型不具有实际应用价值, 需要综合考虑模型泡孔规模与得到可信结论之间的关系。为此, 分别构建了泡孔规模δ=2、δ=3和δ=4, 孔隙率范围35%~65%的一系列三维随机泡孔模型, 并采用Abaqus有限元软件分析了泡沫杨氏模量随相对密度的分布规律, 见图 2。
可见, 泡孔规模δ=3时, 结果的离散性较小, 且模型没有明显尺度效应。所以采用泡孔规模δ=3的代表性体积单元模型用于后续分析。因此, 以下泡沫模型的动态力学性能模拟中, 采用的代表性体积单元模型尺寸为6 mm×6 mm×6 mm, 泡孔半径为1 mm。为便于施加压缩载荷, 在模型两端各添加一块刚性平板, 并对其中一块平板施加固支约束, 对另一块刚性平板施加恒定速率为v的位移载荷。刚性平板与泡孔模型接触, 泡孔模型自接触。泡沫材料的应变率可以表示为:
˙ε=vl (3) 式中:v为加载速率, l为代表性体积单元模型沿加载方向的尺寸。
以下计算采用Ansys/Ls-dyna有限元程序, 不考虑压缩过程中气体压力变化的影响。本文中应力、应变均指工程应力和工程应变。
2. 计算结果
2.1 应变率敏感因素
图 3为有限元计算的相对密度50%泡沫模型动态压缩应力应变曲线, 应变率为10~104 s-1。可见, 泡沫模型的压缩应力应变曲线显示出了典型的三阶段变形特征, 即线弹性段、塑性屈服平台段以及致密段; 且随着应变率的增加, 应力应变曲线的平台段逐渐抬高, 泡沫模型的屈服应力与应变率大小相关, 泡沫模型表现出明显的应变率敏感效应。
由此可以得到初步结论, 若基体材料对应变率敏感, 则泡沫铝材料也存在应变率敏感效应。至于泡沫铝应变率敏感效应是由基体材料的率相关性能引起的还是微结构惯性导致的, 还需要进一步分析。
当基体材料率无关时, 相对密度50%泡沫模型的动态压缩应力应变曲线见图 4。可见, 泡沫模型的压缩应力应变曲线同样表现出三阶段变形特征; 在中、低应变率(10~1 000 s-1)下, 随着应变率的增加, 应力应变曲线的平台段变化不大, 泡沫模型的压缩力学性能对应变率不敏感(见图 4(a)); 在高应变率(2 500~10 000 s-1)下, 随着应变率的增加, 应力应变曲线的平台段逐渐抬高, 泡沫模型表现出一定的应变率敏感特征(见图 4(b))。与图 3相比, 图 4(b)中平台应力随应变率的变化幅度小, 可见在高应变率下, 基体材料的率相关效应对泡沫模型的率相关性能仍具有重要影响。
结合前面的分析可以得到以下结论:(1)在中、低应变率下, 泡沫铝应变率的敏感效应主要由基体材料的率相关性能引起, 微结构惯性的影响很小。(2)在高应变率下, 微结构惯性对泡沫铝应变率效应的影响逐渐增大, 此时泡沫铝率相关性能受基体材料的率相关性能和微结构惯性联合作用。
2.2 相对密度的影响
图 5为10~10 000 s-1应变率范围动态压缩载荷条件下, 泡沫铝流动应力随相对密度的分布。其中, 流动应力为各模型应力应变曲线上与0.2压缩应变所对应的平台段应力[9]。
L.J.Gibson等[13]给出了静态压缩条件下泡沫材料的塑性坍塌强度与相对密度关系式:
σplσys=aρn (4) 式中:σpl为泡沫材料的塑性坍塌强度, σys为基体材料的屈服强度, a和n为材料常数, ρ为泡沫相对密度。
采用上式对图 5中计算结果进行拟合, 可以得到不同应变率下a和n的结果, 见表 1, 拟合曲线见图 5。由表 1中的拟合计算结果发现:(1)随着应变率的提高, 比例系数a逐渐增大, 泡沫模型表现出应变率强化效应。(2)在中、低应变率下, 随着应变率的提高, 指数n保持不变, 表明相对密度变化对泡沫模型的应变率敏感程度基本没有影响。(3)在高应变率下, 随着应变率的提高, 指数n逐渐减小, 表明泡沫模型的应变率敏感程度与相对密度相关。
表 1 拟合结果Table 1. Fitting results of parameters˙ε/s-1 a n 10 0.698 1.72 100 0.784 1.72 1 000 0.942 1.72 2 500 1.03 1.68 5 000 1.13 1.66 10 000 1.27 1.64 2.3 应变率的影响
由式(4)可以得到应变率由
增加到
时, 泡沫铝材料塑性坍塌强度的相对变化量为:
αf=σpl(˙ε2)−σpl(˙ε1)σpl(˙ε1) (5) 当变化前后的应变率
和
位于中、低应变率范围内时, 注意到此时式(4)中系数n保持不变, 式(5)变为:
αf=a(˙ε2)−a(˙ε1)a(˙ε1) (6) 式中:
和
为应变率分别为
和
时, 对应的式(4)中比例系数a。
同样, 由式(4)可以得到应变率从
增加到
时, 泡沫铝基体材料屈服强度的相对变化量:
αs=σdy(˙ε2)−σdy(˙ε1)σdy(˙ε1)=(˙ε2)1/P−(˙ε1)1/PC1/P+(ˉε1)1/P (7) 基于上式计算可以得到, 当应变率从10 s-1增加到100和1 000s-1时, 泡沫铝基体材料屈服强度分别增加12.87%和35.74%。将表 1中的数值拟合结果代入式(6)则可以得到, 当应变率从10 s-1增加到100和1 000 s-1时, 泡沫铝材料塑性坍塌强度分别增加12.32%和34.96%。可见, 在中、低应变率下, 泡沫铝与基体材料对应变率的敏感程度接近, 此时, 基体材料的率相关性是泡沫铝应变率效应的主要影响因素。
分别基于式(5)和(7)可以得到, 应变率从2 500 s-1增加到5 000 s-1和10 000 s-1时, 相对密度分别为35%、50%和65%泡沫模型及其基体材料塑性坍塌强度的相对变化量, 见表 2。可见, 泡沫模型塑性坍塌强度的相对变化量比基体材料大; 随着相对密度升高, 泡沫模型塑性坍塌强度的相对变化量逐渐减小。(1)在高应变率下, 微结构惯性的影响逐渐增大, 使泡沫材料对应变率的敏感性进一步提高, 导致泡沫铝比铝合金材料具有更大的应变率敏感性。(2)在高应变率下, 随着相对密度的降低, 泡沫模型的应变率敏感程度提升。这表明相对密度较低泡沫铝的微结构惯性效应更加显著。
表 2 泡沫模型及其基体材料强度的变化Table 2. Variation of flow stress for foams and its cell material˙ε2/s-1 αf/% αs/% ρ=0.35 ρ=0.50 ρ=0.65 5 000 12.04 11.24 10.66 8.34 10 000 28.59 26.77 25.44 18.25 3. 结论
基于三维随机球形泡孔模型, 采用有限元方法, 对35%~65%相对密度范围内泡沫铝在10~10 000 s-1应变率下的动态压缩变形进行了数值模拟。计算结果表明:
(1) 在中、低应变率(10~1 000 s-1)下, 基体材料的率相关性是泡沫铝应变率效应的主要影响因素; 在高应变率(2 500~10 000 s-1)下, 泡沫铝应变率效应受基体材料率相关性能和微结构惯性共同作用。
(2) 在中、低应变率下, 泡孔结构的变形主要表现为准静态变形模式; 在高应变率下, 动态压缩变形机制逐渐起重要作用。
(3) 在中、低应变率下, 泡沫铝与基体材料对应变率的敏感程度接近; 在高应变率下, 泡沫铝比铝合金材料具有更大的应变率敏感性。
(4) 在高应变率下, 相对密度较低泡沫铝的微结构惯性效应更显著。
-
表 1 试验构件参数
Table 1. Parameters of test members
编号 原文编号 m/kg L/m v0/(m·s−1) Do/mm h/mm fs/MPa fc/MPa w*f/mm 来源 1 CC1 465 1.74 9.10 180 3.65 247 65.74 57.00 文献[7] 2 CC2 920 1.74 6.40 180 3.65 247 65.74 60.00 3 CC3 465 1.74 9.67 180 3.65 247 65.74 72.00 4 DBF14 229.8 1.2 3.96 120 1.70 232 33.73 19.44 文献[5] 5 DBF16 229.8 1.2 4.14 120 1.70 232 47.5 25.66 6 DZF22 229.8 1.2 7.67 120 3.50 298 47.5 39.42 7 DZF23 229.8 1.2 9.90 120 3.50 298 47.5 63.78 8 DZF24 229.8 1.2 10.19 120 3.50 298 47.5 65.40 9 DZF25 229.8 1.2 8.93 120 3.50 298 47.5 72.42 10 DZF28 229.8 1.2 11.54 120 3.50 298 47.5 79.42 11 DZF30 229.8 1.2 11.63 120 3.50 298 47.5 82.30 12 DHF35 229.8 1.2 10.84 120 4.50 290 47.5 33.06 13 DHF36 229.8 1.2 14.48 120 4.50 290 47.5 73.24 14 DHF37 229.8 1.2 14.00 120 4.50 290 47.5 56.20 15 DHF39 229.8 1.2 11.71 120 4.50 290 47.5 38.30 16 DHF40 229.8 1.2 12.52 120 4.50 290 47.5 48.10 表 2 有限元算例参数
Table 2. Parameters of the FE model
组别 m0/kg v0/(m·s−1) Do/mm h/mm fc/MPa ρc/(kg·m−3) B/mm Ⅰ 920 16~6.75 300~120 3.6~9.0 60 2 440 27 Ⅱ 920 8.1~11.2 180 1.8~10.8 60 2 440 27 Ⅲ 920 9.4 180 5.4 60 2 440 9~450 Ⅳ 920 9.4 180 5.4 40~80 2 300~2 600 27 Ⅴ 150~1000 2.8~14.3 180 5.4 60 2 440 27 -
[1] HAH L H, LI W, BJORHOVDE R. Developments and advanced applications of concrete-filled steel tubular (CFST) structures: members [J]. Journal of Constructional Steel Research, 2014, 100: 211–228. DOI: 10.1016/j.jcsr.2014.04.016. [2] 余同希, 朱凌, 许骏. 结构冲击动力学进展 (2010−2020) [J]. 爆炸与冲击, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113.YU T X, ZHU L, XU J. Progress in structural impact dynamics during 2010−2020 [J]. Explosion and Shock Waves, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113. [3] FAN W, XU X, ZHANG Z, et al. Performance and sensitivity analysis of UHPFRC-strengthened bridge columns subjected to vehicle collisions [J]. Engineering Structures, 2018, 173: 251–268. DOI: 10.1016/j.engstruct.2018.06.113. [4] XU M, GAO S, GUO L H, et al. Study on collapse mechanism of steel frame with CFST-columns under column-removal scenario [J]. Journal of Constructional Steel Research, 2018, 141: 275–286. DOI: 10.1016/j.jcsr.2017.11.020. [5] 王蕊. 钢管混凝土结构构件在侧向撞击下的动力响应及其损伤破坏的研究 [D]. 太原: 太原理工大学, 2008: 19–49.WANG R. Study on the dynamic response and damage failure of concrete filled steel tube under lateral impact [D]. Taiyuan, Shanxi, China: Taiyuan University of Technology, 2008: 19–49. [6] WANG R, HAN L H, HOU C C. Behavior of concrete filled steel tubular (CFST) members under lateral impact: experiment and FEA model [J]. Journal of Constructional Steel Research, 2013, 80: 188–201. DOI: 10.1016/j.jcsr.2012.09.003. [7] HAN L H, HOU C C, ZHAO X L, et al. Behaviour of high-strength concrete filled steel tubes under transverse impact loading [J]. Journal of Constructional Steel Research, 2014, 92: 25–39. DOI: 10.1016/j.jcsr.2013.09.003. [8] 王潇宇, DEMARTINO C, 徐金俊, 等. 侧向冲击作用下钢管混凝土柱动力响应试验研究及计算方法 [J]. 土木工程学报, 2017, 50(12): 28–36. DOI: 10.15951/j.tmgcxb.2017.12.004.WANG X Y, DEMARTINO C, XU J J, et al. Dynamic response of concrete filled steel tube column under lateral impact load: experimental study and calculation method [J]. China Civil Engineering Journal, 2017, 50(12): 28–36. DOI: 10.15951/j.tmgcxb.2017.12.004. [9] ZHU A Z, XU W, GAO K, et al. Lateral impact response of rectangular hollow and partially concrete-filled steel tubular columns [J]. Thin-Walled Structures, 2018, 130: 114–131. DOI: 10.1016/j.tws.2018.05.009. [10] 王文达, 陈振福, 纪孙航. 长期持荷工况下钢管混凝土构件的抗撞击性能 [J]. 爆炸与冲击, 2021, 41(8): 083106. DOI: 10.11883/bzycj-2020-0204.WANG W D, CHEN Z F, JI S H, Impact resistance of concrete-filled steel tubular members under long-term loading [J]. Explosion and Shock Waves, 2021, 41(8): 083106. DOI: 10.11883/bzycj-2020-0204. [11] 王路明, 刘艳辉, 赵世春, 等. 侧向低速冲击作用下钢管混凝土构件开裂评估模型及影响因素研究 [J]. 土木工程学报, 2022, 55(3): 7–17, 35. DOI: 10.15951/j.tmgcxb.2022.03.001.WANG L M, LIU Y H, ZHAO S C, et al. Study on evaluation and influencing factors for cracking of concrete-filled steel tubular members subjected to lateral low-velocity impact [J]. China Civil Engineering Journal, 2022, 55(3): 7–17, 35. DOI: 10.15951/j.tmgcxb.2022.03.001. [12] 瞿海雁, 李国强, 孙建运, 等. 侧向冲击作用下钢管混凝土构件的简化分析模型 [J]. 同济大学学报 (自然科学版), 2011, 39(1): 35–41. DOI: 0253-374X(2011)01-0035-07.QU H Y, LI G Q, SUN J Y, et al. Simplified analysis model of circular concrete filled steel tube specimen under lateral impact [J]. Journal of Tongji University (Natural science), 2011, 39(1): 35–41. DOI: 0253-374X(2011)01-0035-07. [13] BAMBACH M R, JAMA H, ZHAO X L, et al. Hollow and concrete filled steel hollow sections under transverse impact loads [J]. Engineering Structures, 2008, 30(10): 2859–2870. DOI: 10.1016/j.engstruct.2008.04.003. [14] WANG Y, QIAN X, LIEW J Y R, et al. Impact of cement composite filled steel tubes: an experimental, numerical and theoretical treatise [J]. Thin-Walled Structures, 2015, 87: 76–88. DOI: 10.1016/j.tws.2014.11.007. [15] 余同希, 华云龙. 结构塑性动力学引论 [M]. 合肥: 中国科技大学出版社, 1994: 88−89. [16] MARTIN J B, SYMONDS P S. Mode approximations for impulsively loaded rigid plastic structures [J]. Journal of the Engineering Mechanics Division, 1965, 92(5): 61. DOI: 10.1061/JMCEA3.0001036. [17] YU T X, STRONGE W J. Large deflections of a rigid-plastic beam-on-foundation from impact [J]. International Journal of Impact Engineering, 1990, 9(1): 115–126. DOI: 10.1016/0734-743x(90)90025-q. [18] 于博丽, 冯根柱, 李世强, 等. 横向爆炸载荷下薄壁圆管的动态响应 [J]. 爆炸与冲击, 2019, 39(10): 103101. DOI: 10.11883/bzycj-2018-0295.YU B L, FENG G Z, LI S Q, et al. Dynamic response of thin-wall circular tubes under transverse blast loading [J]. Explosion and Shock Waves, 2019, 39(10): 103101. DOI: 10.11883/bzycj-2018-0295. [19] WALTERS R M, JONES N. An approximate theoretical study of the dynamic plastic behavior of shells [J]. International Journal of Non-Linear Mechanics, 1972, 7(3): 255–273. DOI: 10.1016/0020-7462(72)90049-2. [20] 张煜航, 陈青青, 张杰, 等. 混凝土三维细观模型的建模方法与力学特性分析 [J]. 爆炸与冲击, 2019, 39(5): 054205. DOI: 10.11883/bzycj-2018-0408.ZHANG Y H, CHEN Q Q, ZHANG J, et al. 3D mesoscale modeling method and dynamic mechanical [J]. Explosion and Shock Waves, 2019, 39(5): 054205. DOI: 10.11883/bzycj-2018-0408. -