• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

钝感炸药点火增长模型的欧拉数值模拟

郝鹏程 冯其京 洪滔 王言金

王元博, 王肖钧, 卞梁, 余育苗. CDM模型及在纤维增强层合材料侵彻数值模拟中的应用[J]. 爆炸与冲击, 2008, 28(2): 172-177. doi: 10.11883/1001-1455(2008)02-0172-06
引用本文: 郝鹏程, 冯其京, 洪滔, 王言金. 钝感炸药点火增长模型的欧拉数值模拟[J]. 爆炸与冲击, 2012, 32(3): 243-250. doi: 10.11883/1001-1455(2012)03-0243-08
WANG Yuan-bo, WANG Xiao-jun, BIAN Liang, YU Yu-miao. CDM model and its application to numerical simulation on fiber-reinforced laminate under penetration[J]. Explosion And Shock Waves, 2008, 28(2): 172-177. doi: 10.11883/1001-1455(2008)02-0172-06
Citation: HAO Peng-cheng, FENG Qi-jing, HONG Tao, WANG Yan-jin. Euleriansimulationoninsensitiveexplosives withtheignition-growthreactivemodel[J]. Explosion And Shock Waves, 2012, 32(3): 243-250. doi: 10.11883/1001-1455(2012)03-0243-08

钝感炸药点火增长模型的欧拉数值模拟

doi: 10.11883/1001-1455(2012)03-0243-08
基金项目: 

国家自然科学基金项目(10926188)

中国工程物理研究院科学技术发展基金项目(2009B0202024)

计算物理实验室基金项目(9140C6901011001)

爆炸科学与技术国家重点实验室开放基金项目(KFJJ11-5M)

详细信息
    通讯作者:

    郝鹏程

Euleriansimulationoninsensitiveexplosives withtheignition-growthreactivemodel

Funds: 

NationalNaturalScienceFoundationofChina(10926188)

More Information
    Corresponding author: HAO Peng-cheng
  • 摘要: 在自主研发的二维多介质欧拉弹塑性流体力学程序中,通过引入点火增长的反应率模型以及炸药 减敏模型,借助网格自适应技术,研究钝感炸药的冲击点火、直径效应以及死区形成等爆轰现象。数值模拟结 果表明,该程序能够正确模拟平面爆轰波的爆速、CJ状态、vonNeumann尖点等爆轰参数;并能够较好模拟炸 药的直径效应。另外,通过引入考虑减敏效应的反应率模型,能较好地模拟钝感炸药的死区形成过程。
  • 金属材料的宏观力学性能与微细观结构密切相关。随着晶粒尺寸减小,晶界数量增大。晶界会阻碍位错运动,导致金属材料的屈服强度增大[1-3],屈服强度与晶粒尺寸之间满足Hall-Petch关系[4-5]。晶粒尺寸不仅影响材料的屈服行为,而且对材料的其他力学性能也有明显影响。

    层裂是材料在卸载波作用下发生的一种典型的动态拉伸断裂形式,细观上归结于材料中微空穴的成核、增长与聚集,强烈依赖于材料的微观结构。针对晶粒结构对多晶金属层裂特性的影响已开展了较多的研究,部分研究结果显示,层裂强度与晶粒尺寸之间服从反Hall-Petch关系,即晶粒越大,层裂强度越大[6-7];还有的研究结果显示晶粒尺寸对金属层裂特性几乎没有影响[8]。针对高纯铝和商业1050纯铝,Trivedi等[9]的研究结果显示:在低冲击应力下,高纯铝和1050合金铝的层裂响应基本不依赖于晶粒尺寸;在强冲击应力下,高纯铝的层裂强度随晶粒尺寸的增大而增大,而1050合金铝则保持不变。Chen等[10]则发现,在弱冲击应力下,晶粒尺寸对1060铝的层裂强度有影响,但在强冲击应力下,晶粒尺寸的影响不明显。除了金属铝,学者们对其他延性金属也开展了研究。Escobedo等[11]和Chen等[12]发现,高纯铜的层裂强度不依赖于晶粒尺寸,但损伤分布、平均孔洞尺寸、断裂模式强烈依赖于晶粒尺寸。近年来,也有一些学者关注轧制过程引起的材料微观结构变化对材料层裂特性的影响[13-14],发现当冲击加载方向垂直于轧制方向时,材料的层裂强度较低,而当冲击加载方向沿板材轧制方向时,材料的层裂强度增大。总之,晶粒结构对多晶金属层裂特性的影响非常复杂,有待于进一步研究。

    本文中,以轧制高纯铝(纯度大于99.99%)板材为研究对象,通过不同热处理工艺改变晶粒尺寸,采用一级轻气炮加载技术,开展一维应变冲击加载下高纯铝板材的层裂实验研究,重点关注晶粒尺寸对高纯铝板材宏观层裂强度的影响;对软回收靶板内部的损伤演化和断口进行金相显微分析,讨论晶粒尺寸对高纯铝板材损伤演化和断裂微观机制的影响规律。

    高纯铝粉由贵州铝厂提供,纯度大于99.99%。委托北京有色金属研究院制备成10 mm厚度板材,通过3种不同热处理工艺(300、450和550 ℃下保温2 h退火处理)改变板材的晶粒尺寸。采用线切割技术,沿板材厚度方向切割出直径为50 mm、厚 6 mm的圆片状层裂靶板,靶板取样如图1所示,其中RD为轧制方向,TD为板材厚度方向,ND为RD的垂直方向。飞片材料同为高纯铝,直径为53 mm,厚3 mm。图2(a)(c)给出了3种晶粒尺寸高纯铝靶板纵截面(TD和ND轴组成的平面)的电子背散射衍射(electron back scattering diffraction,EBSD)表征照片,可以看出,随着退火温度升高,晶粒尺寸逐渐增大,晶粒尺寸分布较均匀,但在550 ℃条件下晶粒尺寸非常大,在EBSD有效表征区内难以见到完整的晶粒。为此,图2(d)给出了由8个图像拼接而成的550 ℃条件下的高纯铝晶粒分布,其平均晶粒尺寸统计值约为500 μm。300和450 ℃条件下晶粒尺寸统计值分别为60和100 μm。图3给出了300 ℃条件下高纯铝板微结构EBSD[10]表征的三维照片,显示高纯铝板材在纵截面和横截面上晶粒尺寸分布较均匀,但是在轧面上晶粒被明显拉长,晶界沿轧制方向取向。

    图  1  板材轧制方向指示和层裂靶板取样示意图
    Figure  1.  Schematic showing the three principal orthogonal directions ofrolled plate and sample preparation
    图  2  不同热处理工艺得到的高纯铝微结构EBSD表征照片
    Figure  2.  Microstructures of high-purity (HP) aluminum by EBSD at different heat-treatment temperatures
    图  3  300 ℃条件下高纯铝微结构EBSD表征的三维照片
    Figure  3.  3D microstructure of HP aluminum by EBSD at the heat-treatment temperature of 300 °C

    采用内径为57 mm的一级轻气炮开展一维应变条件下高纯铝板材的层裂实验,实验装置如图4所示。飞片安装在铝合金弹托上,在高压氮气驱动下沿着炮管加速,通过调整驱动气体的压力控制飞片速度,由安装于炮口位置的测速探针测量飞片击靶速度。层裂实验基本原理为:飞片高速撞击靶板,在飞片和靶板中分别产生相背而行的两束冲击压缩波进入飞片和靶板中,当两束冲击压缩波分别到达飞片和靶板的自由表面时反射成为两束稀疏波,两束稀疏波在靶体中间位置相遇,从而产生拉伸应力,使得材料发生损伤断裂。实验中,采用全光纤激光位移干涉测量仪(displacement interferometer system for any reflector,DISAR)[15]测量靶板后自由面的粒子速度时程曲线。实验采用的飞片直径大于靶板直径,飞片撞击靶板后被靶环挡住,避免了飞片和靶板之间发生二次碰撞。通过合理设计回收装置,实现靶板软回收。对软回收靶板内部的微损伤演化和断口形貌进行显微分析,讨论晶粒尺寸的影响。为了单一研究晶粒尺寸对高纯铝板材层裂特性的影响,一组实验中尽量保持飞片撞击速度不变,实验条件及相关结果列于表1,其中:D为晶粒直径,v为撞击速度,σp为冲击加载的峰值应力,Δu为回跳(pull-back )速度,σs为材料层裂强度,˙uf1˙uf2分别为卸载段的自由面速度变化率和回跳速度斜率,˙ε为拉伸应变率。

    图  4  层裂实验装置示意图
    Figure  4.  Schematic diagram of spall experiments
    表  1  高纯铝板材层裂实验条件和结果
    Table  1.  Conditions and results of spall experiment for high-purity aluminum plate
    No.D/μmv/(m·s−1)σp/GPaΔu/(m·s−1)σs/GPa˙uf1/(m·s−2)˙ε/s−1˙uf2/(m·s−2)
    1 601421.0556.80.411.05×1080.98×1042.70×107
    21001391.0358.80.421.23×1081.14×1044.16×107
    35001401.0478.90.571.30×1081.21×1047.35×107
    4 602501.8785.10.622.68×1082.49×1042.80×107
    51002461.8486.90.632.61×1082.43×1045.16×107
    65002521.8986.10.622.58×1082.40×1047.92×107
    下载: 导出CSV 
    | 显示表格

    采用低速飞片撞击靶板,在靶板中产生初始层裂状态,中间剖切的高纯铝靶板如图5(a)所示,靶板整体发生弯曲变形,但没有观察到断裂面和飞片发生二次碰撞痕迹。随着飞片击靶速度提高,靶板发生完全层裂,横截面上观察到连贯的断裂面,如图5(b)所示。图6(a)(b)分别给出了不同飞片撞击速度下3种晶粒尺寸靶板的自由面粒子速度时程曲线。自由面速度时程曲线包含了丰富的波传播和波系相互作用信息,可以间接反映材料内部发生的微损伤演化动力学过程。对比图6(a)(b) 发现:无论是初始层裂还是完全层裂,自由面速度时程曲线都出现了明显的回跳层裂信号。传统观点一直把回跳信号作为判断层裂是否发生的重要依据。事实上,对回跳速度信息的解读一直存在争议[16]。以晶粒尺寸为500 μm靶板的实测自由面速度时程曲线为例,定量讨论自由面速度波形的典型特征,如图6(a)所示。

    图  5  软回收靶板照片
    Figure  5.  Images of soft recovery HP aluminum samples
    图  6  不同撞击速度下3种晶粒尺寸的高纯铝靶板实测自由面速度曲线
    Figure  6.  Free-surface velocity profiles of HP aluminum with different grain sizes at different impact velocities

    利用自由面速度曲线上的粒子速度峰值,可以计算冲击加载的峰值应力[17]

    σp=12ρ(c0+12λu)u
    (1)

    式中:ρ为材料的密度,c0λ分别为材料的Hugoniot参数,u为自由面粒子速度的峰值。高纯铝材料参数为ρ=2.7 g/cm3c0 = 5 386 m/s,λ=1.34。

    利用卸载段的自由面速度变化率˙uf1来估算层裂过程中的拉伸应变率[17]

    ˙ε=˙uf12cb
    (2)

    式中:cb为体波声速。

    利用回跳速度Δu可以近似计算出材料的层裂强度[17]

    σs=12ρcbΔu
    (3)

    冲击应力峰值、拉伸应变率和层裂强度的计算结果见表1图7给出了层裂强度随晶粒尺寸的变化曲线,图7中包含了Trivedi等[9]的实验结果。本实验结果显示:高纯铝板材在低冲击应力下的层裂强度随晶粒尺寸的增大而增大,体现了反Hall-Petch关系;在高冲击应力下,晶粒尺寸对层裂强度的影响较弱。Trivedi等[9]则给出了相反的结果:在高冲击应力下,晶粒尺寸和层裂强度之间服从反Hall-Petch关系;而在低冲击应力下,层裂强度与晶粒尺寸之间没有明显的依赖关系。无论是本文中的实验结果还是Trivedi等[9]的实验结果都显示,随着冲击应力增大,层裂强度增大。事实上,强冲击加载下,材料经历严重塑性变形,从而引起材料微观结构变化。材料初始微观结构及冲击应力都会对材料抵抗拉伸层裂破坏的能力产生显著影响,并且这种影响相互耦合。后续将通过微损伤演化分析及EBSD微观表征来进一步讨论层裂强度与材料微观结构之间的复杂关系。

    图  7  不同冲击应力条件下高纯铝层裂强度随晶粒尺寸的变化关系
    Figure  7.  Relationship between spallation strength and grain size of HP aluminum under different impact stress conditions

    除回跳速度Δu外,回跳速度斜率˙uf2也是反映材料内部损伤演化过程的重要物理量。已有研究指出,˙uf2与内部损伤演化速率之间是强关联的[18]。基于特征线方法,Kanel等[18]建立了回跳速度斜率˙uf2与材料损伤演化速率˙vv的关系, 具体形式为:

    ˙uf2=A(˙vvB)
    (4)

    式中:AB为材料参数。

    式(4)表明回跳速度斜率与材料损伤演化速率成正比。˙uf2随晶粒尺寸增大而增大,表明大晶粒尺寸的高纯铝靶板中损伤演化速率大。微孔洞长大过程与位错运动引起塑性变形区的发展密切相关[19]。屈服应力或流动应力对微孔洞增长率起着重要作用。图8给出了不同晶粒尺寸高纯铝板材在应变率为10−2 s−1时的应力应变曲线,由图8可以看出:随着晶粒尺寸增大,材料的屈服应力显著减小,满足Hall-Petch关系。大晶粒尺寸的高纯铝板材因具有较低的屈服应力,微孔洞长大较容易,导致损伤演化速率较高,在自由面速度曲线上则表现为回跳速度斜率随着晶粒尺寸的增大而增大。

    图  8  不同晶粒尺寸高纯铝板材的准静态应力应变曲线
    Figure  8.  Stress strain curves of HP aluminum with different grain sizes

    为了研究初始层裂靶板内部的微损伤演化特征,首先采用线切割方式将靶板(如图5所示)从中间分开,然后对损伤区域进行局部切割,再进行表面研磨和抛光,最后利用超景深光学显微镜观察靶板横截面上的损伤分布特征。图9(a)(c)给出了3种晶粒尺寸靶板内部微损伤分布的光学照片,靶板内部没有形成连续的断裂面,还处于微损伤状态,但损伤分布特征明显依赖于材料晶粒尺寸。微损伤带的宽度随晶粒尺寸的增大而增大,部分微孔洞之间发生相互聚集现象,聚集优先方向为垂直于冲击方向。这与祁美兰等[20]、Qi等[21]给出的高纯铝棒材层裂实验结果(见图9(d))明显不同,图9(d)清晰地显示微孔洞沿冲击方向聚集。对于高纯多晶金属,微孔洞成核、长大于晶界(见图9(e)),微孔洞之间易于沿着晶界发生聚集。在冲击加载方向与晶粒取向一致时,微孔洞沿着晶界聚集,不容易形成断裂带,因而材料具有较强的抗层裂破坏能力。此外,从图9中还可以看出:随着晶粒尺寸增大,微孔洞的数量减少而尺寸增大;对于晶粒尺寸较小的靶板,材料内部晶界较多,导致微孔洞成核点多,在相同冲击能量作用下,微孔洞的数量多但尺寸小,其中大孔洞是由于微孔洞之间聚集形成的;而对于晶粒尺寸大的靶板,由于内部晶界较少,微孔洞成核点少,因而不能排除微孔洞会成核于晶粒内部,如图9(f)所示,在相同冲击能量作用下,横截面上显示的微孔洞数量少但尺寸大,并且分布不均匀。3种晶粒尺寸靶板内部的微孔洞尺寸统计结果如图10所示。

    图  9  软回收不同晶粒尺寸靶板中微孔洞分布光学照片
    Figure  9.  Micro-void distribution in the HP aluminum plate samples with different grain sizes
    图  10  不同晶粒尺寸靶板中微孔洞尺寸统计结果
    Figure  10.  Statistical distribution of micro-void size in the HP aluminum plate samples with different grain sizes

    为了进一步揭示微孔洞长大和聚集的微观机理,图11给出了不同晶粒尺寸高纯铝靶板横截面损伤分布的EBSD表征照片。承受强冲击的靶板由于发生严重的塑性变形(特别是在微孔洞周围),同时靶板内残余应力增大,导致EBSD的标定率下降,给EBSD表征带来困难。对于标定率差的靶板,只能通过减小表征面积来提高标定率,因而图11中各图的图形标尺不同。从图11中可以看出:原始晶粒尺寸为60和100 μm的靶板,晶粒被压扁,且晶粒细化现象(尤其是微孔洞周围)非常严重;原始晶粒尺寸为500 μm的靶板,晶粒变化不明显,晶界和晶粒内部都观察到了微孔洞。众所周知,在层裂实验中,总是靶板先受到冲击压缩,然后自由表面反射的稀疏波相互作用形成冲击拉伸。目前,关于冲击压缩对材料微观结构影响的研究较少,Sencer等[22]研究了冲击压缩加载波形对金属铜和316L不锈钢微观结构和性能的影响,发现冲击压缩后金属铜的屈服应力明显提高,而316L不锈钢的应变硬化效应得到加强。迄今为止,关于晶粒尺寸对材料损伤特性的影响已开展了许多研究工作,但基本都忽略了冲击压缩过程对材料原始微观结构的影响。为了进一步验证冲击压缩对高纯铝板材微观结构的影响,开展了一维应变加载条件下的纯冲击压缩实验,飞片和靶板厚度比为1∶1,飞片撞击速度为200 m/s。

    图  11  不同晶粒尺寸高纯铝损伤后靶板微观结构的EBSD表征照片
    Figure  11.  EBSD characterization of HP aluminum samples with different grain sizes after damage

    图12(a)(b)分别给出了原始晶粒尺寸为 60和100 μm的靶板经历纯冲击压缩变形之后内部微观结构的EBSD表征照片。图13给出了平均晶粒尺寸为60 μm靶板的原始晶粒尺寸和压缩之后晶粒尺寸分布的对比,结果显示:冲击压缩后的靶板晶粒细化现象非常严重,再次证实冲击压缩过程会导致材料微观结构出现显著变化,从而对材料层裂特性产生严重影响。图6图7给出的宏观自由面速度剖面上的不同特征和层裂强度与晶粒尺寸之间的关系,本质上都与材料内部微观结构演化密切相关。裴晓阳等[16]指出,层裂实验中自由面速度回跳点本质上对应于微孔洞的成核。在低速冲击下,晶粒尺寸越大,回跳点的位置越低,这可以归因于成核率的影响。晶粒尺寸大的金属,内部成核点少,成核率低。张凤国等[23]采用数值模拟方法讨论了晶粒尺寸对高纯铜层裂损伤演化过程的影响,也发现随着材料晶粒尺寸增大,自由面速度回跳点降低,而回跳后速度曲线的斜率增大,损伤材料内部的孔洞数减少,平均孔洞尺寸增大。张凤国等[23]的计算结果与本文中的实验结果一致。在高速冲击下,前期的冲击压缩过程已经引起晶粒细化,从而导致初始晶粒尺寸对层裂损伤演化的影响被弱化。

    图  12  纯冲击压缩后不同初始晶粒尺寸的高纯铝靶板微观结构EBSD表征照片
    Figure  12.  EBSD characterization of HP aluminum sampleswith different grain sizes under shock compression
    图  13  高纯铝靶板原始的和冲击压缩后的晶粒尺寸分布
    Figure  13.  Grain size distribution of original and shocked high-purity aluminum plate.

    高速冲击加载下,不同晶粒尺寸的高纯铝靶板均发生完全层裂破坏,图14(a)(c)分别为不同晶粒尺寸靶板的层裂断口SEM照片。由图14可知:(1)3种晶粒尺寸靶板的断口表面都存在大量的韧窝,晶粒尺寸对韧窝的几何特征有显著影响;(2)晶粒尺寸小(60和100 μm)的靶板,断口上呈现高密度小韧窝,表现出明显的韧性沿晶断裂特征;(3)晶粒尺寸为500 μm的靶板,断口上除分布有一些小韧窝外,还有一些大韧窝,韧窝非常浅且表面较光滑,这些大韧窝并不是孔洞扩展聚集形成的,而是晶界断裂的痕迹,表现出明显的脆性沿晶断裂特征;(4)在高纯铝层裂断口还发现了一些金属小圆球,其局部放大如图14(d)所示,通过能谱分析可知,这些金属小圆球的成分为铝,由此推测这些小金属球可能是由于微孔洞长大和聚集过程中,局域塑性变形引起温升而导致材料局域熔化形成的。Liao等[24]采用分子动力学计算方法研究了金属铝层裂过程中微观结构和热耗散对微孔洞成核、长大的影响,发现微孔洞长大与热熔化之间相互促进。微孔洞长大过程中的热耗散导致微孔洞附近的温度升高,并促进微孔洞周围的材料熔化。此外,材料熔化导致材料的强度下降,有利于微孔洞长大。计算结果显示,微孔洞周围的局域温度可以达到1 200 K,远高于高纯铝的熔化温度(933 K)。

    图  14  不同晶粒尺寸靶板的断口SEM表征照片
    Figure  14.  SEM micrographs of fracture surface of HP aluminum with different grain sizes

    通过平板撞击层裂实验研究了晶粒尺寸对高纯铝板材层裂行为的影响,获得了不同飞片击靶速度条件下的自由面速度时程曲线,并对软回收靶板内部损伤演化特性和断口形貌进行了显微分析与讨论,获得以下3点认识。

    (1)基于实测自由面速度曲线,发现晶粒尺寸对高纯铝板材层裂行为的影响依赖于冲击加载应力幅值,在1.04 GPa冲击应力条件下,层裂强度与晶粒尺寸之间表现出反Hall-Petch关系;而在1.87 GPa 冲击应力下,晶粒尺寸对层裂强度几乎没有影响。

    (2)基于软回收样品金相显微分析,发现晶粒尺寸对微损伤分布特性和断口形貌有明显影响,随着晶粒尺寸增大,微孔洞分布范围扩大、数量减少、尺寸增大,断口形貌呈现由韧性向准脆性转变的沿晶断裂特征。

    (3)首次在断口上观察到随机分布的小金属圆球,归因于微孔洞长大、聚集过程中严重塑性变形引起的热效应,热效应不仅会促进微孔洞增长,还可能导致材料发生局部熔化,值得进一步关注。

  • 加载中
计量
  • 文章访问数:  2836
  • HTML全文浏览量:  111
  • PDF下载量:  239
  • 被引次数: 0
出版历程
  • 刊出日期:  2012-05-25

目录

/

返回文章
返回