
Citation: | ZHANG Jiandong, WU Haijun, LI Wei, LI Jinzhu, PI Aiguo. Experiment and Numerical Simulation on Oblique Penetrating Concrete Targets by a Special-Shaped Projectile with Ribbed Head[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 065103. doi: 10.11858/gywlxb.20210723 |
由于爆炸焊接过程具有瞬时性、复杂性,人们在通过实验及相关设备来观察、研究爆炸焊接过程时,效果十分有限。随着计算机技术的发展,数值模拟技术的突破为爆炸焊接研究提供了帮助。迄今为止,国内外诸多学者利用数值模拟软件对爆炸焊接过程及参数开展了大量研究,取得了较理想的研究成果。Nassiri等[1]分别利用光滑粒子流体动力学(Smoothed particle hydrodynamics,SPH)法、任意拉格朗日-欧拉(Arbitrary Lagrange-Euler,ALE)法模拟了爆炸焊接的高速撞击过程,并研究了结合界面的剪切应力和速度分布,结果显示只有SPH法再现了碰撞时的射流现象。Abe[2]使用二维有限差分法研究了爆炸焊接波的形成机理,计算结果与模拟结果保持一致。Yuan等[3]采用SPH法研究了6061Al/AZ31B的爆炸焊接,模拟得到波形界面和射流现象,但是并没有详细指出波的形成机理。Tabbataee等[4]使用有限单元法模拟了射流。Mousavi等[5-6]采用欧拉法模拟了射流和类似流体状的波状界面,并成功预测了波形尺寸和射流速度。王宇新等[7]应用无网格物质点法(Material point method,MPM),对爆炸荷载作用下基板和复板的动态变形过程进行了三维数值模拟,并分析了碰撞点压力和速度的变化情况。刘江等[8]用SPH法进行二维建模,重现了实验中观察到的射流和界面波现象,并揭示了碰撞压力、有效塑性应变、温度等参数随时间的变化关系。周春华等[9-10]研究了均匀布药下结合界面波的变化规律后,提出了新的不等药厚布药工艺。
从众多学者的研究成果来看,数值模拟软件能够帮助人们全面、系统地了解爆炸焊接过程,降低实验成本,优化实验参数。在平板爆炸焊接中,通过炸药爆炸产生的爆轰能量推动复板与基板发生斜碰撞,碰撞界面在微秒级的时间内产生数千兆帕高压,发生塑性变形,从而实现基板和复板的固相焊接。由此可以看出,炸药是整个爆炸焊接过程的能量来源,因此对炸药药量和布药方式进行控制具有重要的研究意义。本研究应用LS-DYNA软件,采用SPH法进行二维建模,研究炸药与爆炸焊接界面波的关系及波的形成机理。
SPH法是一种无网格法,适合处理爆炸、冲击等大变形问题,但是其搜索算法较为耗时,粒子越多,该问题越突出,计算效率越低[11]。为了提高计算结果的精确性,同时兼顾计算效率,炸药采用粗网格划分,基板和复板采用细网格划分。先采用LS-DYNA软件建立二维平面计算模型,再利用LS-prepost软件将网格全部生成SPH粒子(见图1)。复板和基板材料选取Cu和Q235钢,爆炸焊接选用硝铵炸药。模型尺寸:复板150 mm × 4 mm,基板150 mm × 15 mm,间隙6 mm,基板和复板网格尺寸为0.2 mm × 0.2 mm,炸药网格尺寸为0.4 mm × 0.4 mm。采用cm-g-μs单位制。
硝铵炸药选用Jones-Wilkins-Lee(JWL)[12]状态方程(EOS),计算公式为
p=AJ(1−ωR1V)e−R1V+BJ(1−ωR2V)e−R2V+ωE0V |
(1) |
式中:AJ、BJ、R1、R2和ω为材料系数;p为爆轰产物压力,GPa;E0为初始比内能,kJ·cm–3;V为爆轰气体产物的相对比容。炸药的具体JWL状态参数见表1,其中ρ为密度,D为爆速。
ρ/(kg·m−3) | D/(m·s−1) | AJ/GPa | BJ/GPa | R1 | R2 | ω |
800 | 2 800 | 132.75 | 0.423 | 5.3 | 1.2 | 0.21 |
基板、复板均采用 Mie-Grüneisen状态方程[13]和Johnson-Cook材料模型[14]。Grüneisen状态方程被广泛应用于解决材料的大变形问题,它描述了粒子速度up和冲击速度D之间的基本关系,计算公式如下
ρ=ρH+Γρ(e−eH) |
(2) |
其中
Γρ=Γ0ρ0=Constant |
(3) |
ρH=ρ0c0μ(1+μ)[1−(s−1)μ]2 |
(4) |
eH=12ρHρ0μ1+μ |
(5) |
μ=ρρ0−1 |
(6) |
式中:Γ0为Grüneisen系数,ρ为当前密度,ρ0为初始密度,c0为体积声速,s为压缩比。
Johnson-Cook材料模型经常用来反映金属的力学行为,计算公式如下
σ=(A+Bεnp)(1+Cln˙ε∗p)(1−T∗m) |
(7) |
其中无量纲温度T*表示为
T∗=(T−Tr)/(Tm−Tr) |
(8) |
式中:εp为有效塑性应变;
Material | ρ/(g·cm−3) | G/GPa | A/GPa | B/GPa | n | C | m | Tm/K | Tr/K |
Cu | 8.96 | 46 | 0.090 | 0.292 | 0.31 | 0.025 | 1.09 | 1 356 | 294 |
Q235 | 7.83 | 77 | 0.792 | 0.510 | 0.26 | 0.014 | 1.03 | 1 793 | 294 |
Material | c/(km·s−1) | S1 | Γ0 | a |
Cu | 3.940 | 1.489 | 2.02 | 0.47 |
Q235 | 4.569 | 1.490 | 2.17 | 0.46 |
爆炸焊接质量的好坏极其依赖选取的爆炸焊接参数,如药量比、间距、碰撞速度、碰撞角等。参考文献[15]中的爆炸复合窗口,选取复板厚度δ1 = 4 mm时药量比(单位面积炸药与复板质量之比R)分别为R1 = 1.0、R2 = 1.5的两组参数开展研究。经计算得到对应的炸药厚度分别为44.8 mm和67.2 mm。由图2可知,R1靠近可焊窗口下限,R2位于可焊窗口内部,两组参数均位于爆炸复合窗口内部。
为了研究炸药量与焊接界面波形的关系,首先进行均匀布药模拟。分别建立炸药厚度为44.8 mm和67.2 mm的计算模型,输出模拟图像,如图3所示。
为了更方便地了解爆炸焊接过程参数的变化规律,在结合面上等间距区域选取7个关键点,选取位置如图4所示。输出关键点的碰撞压力值,如表4所示。硝铵炸药的爆轰压力计算公式为
Key point | Pressure/GPa | Key point | Pressure/GPa | |||
R1 = 1.0 | R2 = 1.5 | R1 = 1.0 | R2 = 1.5 | |||
A1 | 0.581 | 1.602 | A5 | 4.999 | 7.086 | |
A2 | 1.672 | 2.002 | A6 | 5.754 | 8.576 | |
A3 | 4.507 | 5.289 | A7 | 0.526 | 1.031 | |
A4 | 4.654 | 6.191 |
p=ρ0D21+γ |
(9) |
式中:γ为硝铵炸药的多方指数,ρ0为炸药密度,D为炸药爆速。
由式(9)可以看出,模拟中使用的硝铵炸药的密度为常数,所以爆轰压力p只与炸药爆速D有关。由于炸药的爆速与单位面积的药量成正比[16],且R1 < R2,因此R2对应的爆轰压力要高于R1的爆轰压力。炸药的爆轰压力作用于复板,复板又传递到焊接界面,最终导致R2对应的基板、复板碰撞压力高于R1。根据图5的模拟结果可知,理论分析结果与模拟结果一致。
焊接界面波长的计算公式[17]为
λ=πvpd2c(cvc−1)2 |
(10) |
式中:vp为碰撞速度;d为复板厚度;c为材料声速;vc为碰撞点移动速度,等于炸药的爆轰速度。
从式(10)可以看出,给定金属板材的材料参数后,复板厚度d和材料的声速c为常数,碰撞速度vp和碰撞点移动速度vc为变量,界面波的波长λ与vp成正比,与vc成反比。
从图3可以看出,复板上部均出现了打伤现象,这是由于前期炸药和板材划分的网格尺寸不同所致。R1和R2两组模拟均得到了波状焊接界面,这被认为是焊接良好的标志[18],模拟结果与文献[19]的金相实验结果吻合较好。从图3(a)和图3(b)看出,焊接界面在起爆端区域没有出现波状界面,随着距起爆端距离的不断增大,界面先后出现了无波到微波、微波到大波的转变。图3(a)中射流现象最先出现在4.299 6 μs时刻,而图3(b)中的射流最先出现在4.198 2 μs,早于前者,这是由于药量比越大,爆轰能量越大,导致射流出现较早。焊接界面的放大图像显示,波状界面先于射流出现,可以认为波的形成不一定需要射流,印证了Godunov提出的应力波机理[20]。从图3(c)、图3(d)可以清晰地看出,碰撞点前方在基板和复板碰撞作用下,板材表面喷射出大量射流,可以清洗掉金属表面的杂质和氧化层,提高焊接质量[21]。喷射出的射流中,基板的成分明显高于复板,基板材料是射流成分的主要来源,这是由于Cu的密度高于Q235钢,密度较小的材料更容易产生射流[22]。图3(e)和图3(f)为模拟图像,表明基板与复板已经完全复合。从复合板中央部分的取样结果来看,图3(f)的波状尺寸大于图3(e),由于小波增加了焊接面积,有助于提高焊接质量,波形越大,越容易出现裂缝和孔洞等缺陷,因此可以认为药量比R1优于R2的焊接质量,药量比尽量取在可焊窗口的下限。
图5是均匀布药下各个关键点的碰撞压力折线。考虑到边界效应的影响,去掉首尾两个关键点后,折线显示沿着爆轰方向焊接界面的碰撞压力逐渐增大,由于碰撞压力过大会产生熔化层,产生的稀疏波可能将已结合的复合板拉开,降低焊接质量,因此均匀布药不是最终选择。
根据均匀布药方式下界面波的分布特点,起爆端区域的结合呈平直界面,这是由于炸药起爆时能量不稳定[23],碰撞压力小。在爆炸焊接中,碰撞压力过小会导致界面未复合[16],即起爆端附近极易出现未复合现象,降低金属板材的利用率,这种现象被称为边界效应,在实验和工程实践中应尽量避免。碰撞压力要超过某个临界值,复合板才能获得良好的焊接质量,这也证明了爆炸焊接是一种压力焊[24]。沿着爆轰方向,结合界面逐渐由无波发展成微波、大波界面,这是爆炸产物的堆集以及基板和复板碰撞引起的振动能在待复合区叠加并共同作用的结果[25-26]。结合式(9)、式(10)可知,波状界面的尺寸与碰撞点的移动速度有关,碰撞压力与炸药爆速成正比,在平行法爆炸焊接装置中爆速等于碰撞点的移动速度[27],由此得出,炸药的爆轰能量赋予了复板巨大的压力,碰撞压力决定界面波形。因此可以通过调整炸药厚度来控制爆轰能量,使界面波形维持在微波状,进而提高复合板的焊接质量。
设计如图6所示的呈梯形装药结构,通过调整炸药起爆端和末端厚度a、b的值,来确定一个合适的倾斜装药角度,沿着爆轰方向控制碰撞点的压力基本保持不变,避免因大波状界面而产生涡旋和裂缝,提高复合板的结合强度,使整个焊接界面都保持在微波状结合状态。布药方案见表5。
Scheme | a/mm | b/mm |
Ⅰ | 67.2 | 58.8 |
Ⅱ | 67.2 | 50.4 |
Ⅲ | 67.2 | 42.0 |
Ⅳ | 67.2 | 33.6 |
由于网格划分较大,导致波形尺寸不易测量。当网格划分较小时,可以较清晰地观察界面的形貌特征。由前述分析得知,碰撞压力决定界面波形,因此将界面上的碰撞压力作为研究对象,使界面上各点处碰撞压力基本保持一致。
不同梯形布药方案下界面关键点的碰撞压力如图7所示,具体压力值列于表6。由图7中折线的波动情况可知,当采用梯形布药方式后,基板、复板之间的碰撞压力不再呈均匀布药下逐渐增长的趋势,而是呈较小幅度的上下波动,因此梯形布药方式可以消除不利因素,使焊接界面的碰撞压力基本保持不变,以获得波长较小、波幅较小的微波界面。观察图7发现,起爆端和末端的碰撞压力最小,这是由稀疏波引发的边界效应所致。在不考虑首尾两端的压力时,方案Ⅲ的碰撞压力波动最小,因此最终焊接界面的波形状态好于其他方案。而方案Ⅱ中界面上的碰撞压力起伏最大,最终可能导致基板与复板未复合。因此确定合适的首尾装药高度是梯形布药方案的关键。
Key point | Pressure/GPa | |||
Scheme Ⅰ | Scheme Ⅱ | Scheme Ⅲ | Scheme Ⅳ | |
A1 | 1.738 | 1.158 | 1.694 | 0.351 |
A2 | 2.549 | 2.548 | 2.141 | 1.936 |
A3 | 3.839 | 3.830 | 3.621 | 2.776 |
A4 | 5.230 | 7.547 | 5.439 | 5.665 |
A5 | 7.890 | 8.330 | 3.413 | 4.537 |
A6 | 6.423 | 3.505 | 3.731 | 2.064 |
A7 | 0.780 | 0.431 | 0.397 | 0.435 |
综上所述,采用梯形布药方式可以消除焊接界面波形不均匀现象,而且减少了炸药用量,同时也可以减小爆炸焊接实验现场引发的噪声。通过分析各个关键点的碰撞压力可知,方案Ⅲ的复合效果最好,即炸药起爆端和末端的厚度分别为67.2 mm和42.0 mm时装药倾斜角度最佳。在实际的工程实践中,梯形布药方式的应用并不广泛,原因是其影响因素较多,参数确定比较复杂。尽管减小炸药厚度会使碰撞压力减小,但碰撞压力的变化与炸药厚度并非成严格的正比关系,所以才会得到图7的模拟结果。如何优化梯形布药参数使界面波波形基本保持不变,还需要大量深入的研究。
爆炸焊接过程具有瞬时性,因此很难用实验方法来研究界面波的形成与焊接参数之间的关系。虽然有关界面波形成机理的研究已经进行了数十年,但迄今为止仍然没有达成一个统一的认识[28]。期间众多学者提出了一些波形形成机制,主要有复板流侵彻机理[29]、Helm-holtz不稳定流机理[30-31]、涡脱落机理[32-33]和应力波机理[34]。随着数值模拟软件的发展,SPH法能够清晰地模拟出波形成过程,为解释界面波的形成机理提供了帮助。
图8展示了采用SPH法得到的爆炸焊接界面波形成过程。基板与复板在P点高速碰撞,此时P点可看作低黏塑性流体状态,巨大的碰撞压力使P点产生塑性变形。如图8(a)所示,复板与基板碰撞后形成初始射流和再入射流,再入射流在前方复板的阻碍下向斜上方喷出。沿着焊接方向向前发展,碰撞点P转移到波峰上,见图8(b),先前的再入射流被基板捕获,形成在波峰的前后方,此处极易形成漩涡、空洞等缺陷,从而降低复合板的焊接质量。碰撞点P转移到波峰上后,再入射流与波峰碰撞,沿斜下方喷出。图8(c)中碰撞点P越过波峰后形成新的碰撞点,在巨大的爆轰压力下基板再次发生塑性变形,形成新的波形,由此产生周期性的波状界面。由图8可见,SPH法模拟的波形成过程与复板流侵彻机理的一致性较好,证明了SPH法在解释爆炸焊接界面波形成过程中复板流侵彻机理的有效性。
(1)采用LS-DYNA软件和SPH法,模拟了整个爆炸焊接过程,得到典型的波状界面和射流现象,理论分析与模拟结果具有一致性。
(2)在两种药量比、均匀布药方案中,忽略边界效应后的结果显示,沿着爆轰方向基板与复板之间的碰撞压力均逐渐增大,而且药量比越大,碰撞压力越大,界面波形越大。
(3)在梯形布药方案中,通过改变炸药起爆端和末端的高度设计了4种方案,结果显示梯形布药可以有效消除爆炸焊接界面波的不均匀现象,使界面波形尺寸基本保持一致,而且节省了炸药用量。当起爆端和末端的高度分别为67.2 mm和42.0 mm时效果最佳。
(4)介绍了几种波形成机理,展示并研究了界面波形成过程,模拟的波形成过程与复板流侵彻机理的一致性较好,证明了复板流侵彻机理解释界面波形成过程的有效性。
[1] |
FORRESTAL M J, ALTMAN B S, CARGILE J D, et al. An empirical equation for penetration depth of ogive-nose projectiles into concrete targets [J]. International Journal of Impact Engineering, 1994, 15(4): 395–405. doi: 10.1016/0734-743X(94)80024-4
|
[2] |
FORRESTAL M J, FREW D J, HICKERSON J P, et al. Penetration of concrete targets with deceleration-time measurements [J]. International Journal of Impact Engineering, 2003, 28(5): 479–497. doi: 10.1016/S0734-743X(02)00108-2
|
[3] |
CHEN X W, FAN S C, LI Q M. Oblique and normal perforation of concrete targets by a rigid projectile [J]. International Journal of Impact Engineering, 2004, 30(6): 617–637. doi: 10.1016/j.ijimpeng.2003.08.003
|
[4] |
马兆芳, 段卓平, 欧卓成, 等. 弹体斜侵彻多层间隔混凝土靶实验和数值模拟 [J]. 北京理工大学学报, 2016, 36(10): 1001–1005. doi: 10.15918/j.tbit1001-0645.2016.10.003
MA Z F, DUAN Z P, OU Z C, et al. Experimental and simulative research on projectile oblique penetration into concrete targets with multi-layered space structure [J]. Transactions of Beijing Institute of Technology, 2016, 36(10): 1001–1005. doi: 10.15918/j.tbit1001-0645.2016.10.003
|
[5] |
马兆芳. 动能弹斜侵彻有限厚混凝土靶体的弹道规律研究 [D]. 北京: 北京理工大学, 2016: 12-26.
MA Z F. Investigation on trajectory regularity of kinetic energy projectile oblique penetration into concrete targets of finite thickness [D]. Beijing: Beijing Institute of Technology, 2016: 12-26.
|
[6] |
段卓平, 李淑睿, 马兆芳, 等. 刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型 [J]. 爆炸与冲击, 2019, 39(6): 063302. doi: 10.11883/bzycj-2018-0411
DUAN Z P, LI S R, MA Z F, et al. Analytical model for attitude deflection of rigid projectile during oblique perforation of concrete targets [J]. Explosion and Shock Waves, 2019, 39(6): 063302. doi: 10.11883/bzycj-2018-0411
|
[7] |
KONG X Z, FANG Q, HONG J, et al. Numerical study of the wake separation and reattachment effect on the trajectory of a hard projectile [J]. International Journal of Protective Structures, 2014, 5(1): 97–117. doi: 10.1260/2041-4196.5.1.97
|
[8] |
邓佳杰, 张先锋, 刘闯, 等. 头部非对称刻槽弹体侵彻混凝土目标性能研究 [J]. 兵工学报, 2018, 39(7): 1249–1258. doi: 10.3969/j.issn.1000-1093.2018.07.001
DENG J J, ZHANG X F, LIU C, et al. Research on penetration of asymmetrically grooved nose projectile into concrete target [J]. Acta Armamentarii, 2018, 39(7): 1249–1258. doi: 10.3969/j.issn.1000-1093.2018.07.001
|
[9] |
LIU J C, PI A G, HUANG F L. Penetration performance of double-ogive-nose projectiles [J]. International Journal of Impact Engineering, 2015, 84: 13–23. doi: 10.1016/J.IJIMPENG.2015.05.003
|
[10] |
AMON J, SCHWARTZ A, BRANDEIS Y. Missile warhead: US20150059610A1 [P]. 2015-03-05.
|
[11] |
李金柱, 吕中杰, 张宏松, 等. 弹丸非正侵彻贯穿混凝土靶实验研究 [J]. 弹箭与制导学报, 2013, 33(5): 86–90. doi: 10.3969/j.issn.1673-9728.2013.05.021
LI J Z, LÜ Z J, ZHANG H S, et al. Non-ideal perforation experiments of concrete targets with steel projectiles [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2013, 33(5): 86–90. doi: 10.3969/j.issn.1673-9728.2013.05.021
|
[12] |
JENA P K, JAGTAP N, KUMAR K S, et al. Some experimental studies on angle effect in penetration [J]. International Journal of Impact Engineering, 2010, 37(5): 489–501. doi: 10.1016/j.ijimpeng.2009.11.009
|
[13] |
FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive-nose steel projectiles [J]. International Journal of Impact Engineering, 1996, 18(5): 465–476. doi: 10.1016/0734-743X(95)00048-F
|
[14] |
FREW D J, HANCHAK S J, GREEN M L, et al. Penetration of concrete targets with ogive-nose steel rods [J]. International Journal of Impact Engineering, 1998, 21(6): 489–497. doi: 10.1016/S0734-743X(98)00008-6
|
[1] | CUI Niansheng, WEI Jianlin, YUAN Zengsen, XU Zhenyang, LIU Xin, WANG Xuesong. Simulation Analysis of Mesoscale Characteristics in the Dynamic Fracture Damage of Heterogeneous Rock[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 044204. doi: 10.11858/gywlxb.20230638 |
[2] | XIU Chengdong, WANG Changfeng, LI Bing, GUAN Renguo. Numerical Simulation Study on the Influence of Hard Phase Shape on the Fracture Behavior of Ti-Al3Ti Bionic Composites[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 044201. doi: 10.11858/gywlxb.20230629 |
[3] | ZHANG Haiguang, WANG Yu, AN Lianhao, WANG Ke, WU Xiaodong. Experimental Study and Numerical Simulation of Dynamic Fracture Behavior of Branch Staggered Laminated Biomimetic Composites under Impact Loading[J]. Chinese Journal of High Pressure Physics, 2022, 36(1): 014101. doi: 10.11858/gywlxb.20210776 |
[4] | MENG Xiangsheng, WU Xiaodong, ZHANG Haiguang. Numerical Simulation on Interlaminar Fracture Toughness of 3D Printed Mortar Laminated Composites[J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 044206. doi: 10.11858/gywlxb.20190827 |
[5] | YAO Pengfei, HAN Yang, YAO Fen, LI Zhiqiang. Simulation of the Impact Fracture Behavior of Double Laminated Glass Based on Intrinsic Cohesive Model[J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 064105. doi: 10.11858/gywlxb.20190718 |
[6] | YU Yin, LI Yuanyuan, HE Hongliang, WANG Wenqiang. Mesoscale Lattice Model for Dynamic Fracture of Brittle Materials[J]. Chinese Journal of High Pressure Physics, 2019, 33(3): 030106. doi: 10.11858/gywlxb.20190707 |
[7] | XIAO Jiaxin, BAI Jingsong, WANG Tao. Numerical Study of Shock Wave Impacting on the Double-Mode Interface in Nonuniform Flows[J]. Chinese Journal of High Pressure Physics, 2018, 32(1): 012301. doi: 10.11858/gywlxb.20170501 |
[8] | GUO Chun-Hai, ZHANG Wen-Wu, RU Hao-Lei. Numerical Simulation of High Pressure Micro Water Jet Modulation with the Constraint of Gas Flow[J]. Chinese Journal of High Pressure Physics, 2017, 31(5): 585-592. doi: 10.11858/gywlxb.2017.05.012 |
[9] | PAN Jian-Hua, CHEN Xue-Dong. Numerical Simulation of Dynamic Fracture Toughness Test Using Three-Point Bending Specimen in SHPB[J]. Chinese Journal of High Pressure Physics, 2013, 27(6): 856-862. doi: 10.11858/gywlxb.2013.06.010 |
[10] | LIU Wei, LI Huo-Kun, LIU Cheng-Mei, LIU Wei-Lin. Numerical Simulation of Microchannel of Dynamic High-Pressure Microfluidization Based on FLUENT[J]. Chinese Journal of High Pressure Physics, 2012, 26(1): 113-120. doi: 10.11858/gywlxb.2012.01.017 |
[11] | MA Zhong-Liang, LIU Lin-Lin, XIAO Zhong-Liang. Mathematical Model and Numerical Simulation of Two-Phase Flow Interior Ballistics of Variable Burning-Rate Propellant[J]. Chinese Journal of High Pressure Physics, 2012, 26(1): 102-106. doi: 10.11858/gywlxb.2012.01.015 |
[12] | LIU Fang, LIU Yong-Gang, XIE Hong-Sen. Numerical Simulation of Temperature Field in Sample Assembly of Cubic Press[J]. Chinese Journal of High Pressure Physics, 2012, 26(2): 135-140. doi: 10.11858/gywlxb.2012.02.003 |
[13] | WANG Gang, ZHANG De-Liang, LIU Kai-Xin. Study on Chemical Reaction Models in Gaseous Detonation Numerical Simulation[J]. Chinese Journal of High Pressure Physics, 2008, 22(4): 350-356 . doi: 10.11858/gywlxb.2008.04.003 |
[14] | WANG Yong-Gang, HE Hong-Liang, WANG Li-Li, JING Fu-Qian. Percolation Description for the Early Stage of Void Coalescence during Dynamic Tensile Fracture in Ductile Materials[J]. Chinese Journal of High Pressure Physics, 2006, 20(2): 127-132 . doi: 10.11858/gywlxb.2006.02.003 |
[15] | WU Yu-Yu, HE Yuan-Hang, LI Jin-Zhu. Application of the Coupling Method in Simulating the Hypervelocity Impact[J]. Chinese Journal of High Pressure Physics, 2005, 19(4): 385-389 . doi: 10.11858/gywlxb.2005.04.019 |
[16] | CAO Yu-Zhong, LU Ze-Sheng, GUAN Huai-An, ZHANG You-Ping. Numerical Simulations of Blast Flow-fields in Closed Blast-Resistant Containers[J]. Chinese Journal of High Pressure Physics, 2001, 15(2): 127-133 . doi: 10.11858/gywlxb.2001.02.009 |
[17] | ZHANG Wan-Jia, YANG Zhong-Zheng. Studies on the Fracture Behaviour for 93 Tungsten Alloy[J]. Chinese Journal of High Pressure Physics, 1995, 9(4): 279-288 . doi: 10.11858/gywlxb.1995.04.007 |
[18] | WANG Ze-Ping, HUANG Feng-Lei, DING Jing, HOU Min. Studies of Dynamic Damage Behavior in LY12 Aluminum Alloy under the Condition of High Loading Rate[J]. Chinese Journal of High Pressure Physics, 1993, 7(1): 23-32 . doi: 10.11858/gywlxb.1993.01.003 |
[19] | WANG Ze-Ping, HUANG Feng-Lei, YUN Shou-Rong. A Model for Numerical Calculations of Spallation in Brittle Solids[J]. Chinese Journal of High Pressure Physics, 1991, 5(2): 90-97 . doi: 10.11858/gywlxb.1991.02.002 |
[20] | DONG Yu-Bin, ZAHNG Wan-Jia, JING Fu-Qian, HAN Jun-Wan, CHEN Da-Nian, SU Lin-Xiang, FENG Jia-Bo. Numerical Analysis for Dynamic Damage Processes and LY-12 Aluminum Spallations[J]. Chinese Journal of High Pressure Physics, 1988, 2(4): 305-312 . doi: 10.11858/gywlxb.1988.04.003 |
ρ/(kg·m−3) | D/(m·s−1) | AJ/GPa | BJ/GPa | R1 | R2 | ω |
800 | 2 800 | 132.75 | 0.423 | 5.3 | 1.2 | 0.21 |
Material | ρ/(g·cm−3) | G/GPa | A/GPa | B/GPa | n | C | m | Tm/K | Tr/K |
Cu | 8.96 | 46 | 0.090 | 0.292 | 0.31 | 0.025 | 1.09 | 1 356 | 294 |
Q235 | 7.83 | 77 | 0.792 | 0.510 | 0.26 | 0.014 | 1.03 | 1 793 | 294 |
Material | c/(km·s−1) | S1 | Γ0 | a |
Cu | 3.940 | 1.489 | 2.02 | 0.47 |
Q235 | 4.569 | 1.490 | 2.17 | 0.46 |
Key point | Pressure/GPa | Key point | Pressure/GPa | |||
R1 = 1.0 | R2 = 1.5 | R1 = 1.0 | R2 = 1.5 | |||
A1 | 0.581 | 1.602 | A5 | 4.999 | 7.086 | |
A2 | 1.672 | 2.002 | A6 | 5.754 | 8.576 | |
A3 | 4.507 | 5.289 | A7 | 0.526 | 1.031 | |
A4 | 4.654 | 6.191 |
Scheme | a/mm | b/mm |
Ⅰ | 67.2 | 58.8 |
Ⅱ | 67.2 | 50.4 |
Ⅲ | 67.2 | 42.0 |
Ⅳ | 67.2 | 33.6 |
Key point | Pressure/GPa | |||
Scheme Ⅰ | Scheme Ⅱ | Scheme Ⅲ | Scheme Ⅳ | |
A1 | 1.738 | 1.158 | 1.694 | 0.351 |
A2 | 2.549 | 2.548 | 2.141 | 1.936 |
A3 | 3.839 | 3.830 | 3.621 | 2.776 |
A4 | 5.230 | 7.547 | 5.439 | 5.665 |
A5 | 7.890 | 8.330 | 3.413 | 4.537 |
A6 | 6.423 | 3.505 | 3.731 | 2.064 |
A7 | 0.780 | 0.431 | 0.397 | 0.435 |
ρ/(kg·m−3) | D/(m·s−1) | AJ/GPa | BJ/GPa | R1 | R2 | ω |
800 | 2 800 | 132.75 | 0.423 | 5.3 | 1.2 | 0.21 |
Material | ρ/(g·cm−3) | G/GPa | A/GPa | B/GPa | n | C | m | Tm/K | Tr/K |
Cu | 8.96 | 46 | 0.090 | 0.292 | 0.31 | 0.025 | 1.09 | 1 356 | 294 |
Q235 | 7.83 | 77 | 0.792 | 0.510 | 0.26 | 0.014 | 1.03 | 1 793 | 294 |
Material | c/(km·s−1) | S1 | Γ0 | a |
Cu | 3.940 | 1.489 | 2.02 | 0.47 |
Q235 | 4.569 | 1.490 | 2.17 | 0.46 |
Key point | Pressure/GPa | Key point | Pressure/GPa | |||
R1 = 1.0 | R2 = 1.5 | R1 = 1.0 | R2 = 1.5 | |||
A1 | 0.581 | 1.602 | A5 | 4.999 | 7.086 | |
A2 | 1.672 | 2.002 | A6 | 5.754 | 8.576 | |
A3 | 4.507 | 5.289 | A7 | 0.526 | 1.031 | |
A4 | 4.654 | 6.191 |
Scheme | a/mm | b/mm |
Ⅰ | 67.2 | 58.8 |
Ⅱ | 67.2 | 50.4 |
Ⅲ | 67.2 | 42.0 |
Ⅳ | 67.2 | 33.6 |
Key point | Pressure/GPa | |||
Scheme Ⅰ | Scheme Ⅱ | Scheme Ⅲ | Scheme Ⅳ | |
A1 | 1.738 | 1.158 | 1.694 | 0.351 |
A2 | 2.549 | 2.548 | 2.141 | 1.936 |
A3 | 3.839 | 3.830 | 3.621 | 2.776 |
A4 | 5.230 | 7.547 | 5.439 | 5.665 |
A5 | 7.890 | 8.330 | 3.413 | 4.537 |
A6 | 6.423 | 3.505 | 3.731 | 2.064 |
A7 | 0.780 | 0.431 | 0.397 | 0.435 |