滑移爆轰驱动钢管的层裂 spallation of steel tube driven by sliding detonation.pdfVIP

第28卷第2期 爆炸与冲．击 V01．28．No．2 WAVES 2008年3月 EXPLOSl0NANDSHOCK Mar．．2008 文章编号：1001·1455(2008)02一0105一05 滑移爆轰驱动钢管的层裂。 熊俊，周海兵，刘文韬，张树道，孙锦山 (北京应用物理与计算数学研究所，北京100088) 摘要：将VG损伤模型推广到二维情况，考虑了最大主应力方向对损伤演化的影响，并使用显式断裂算 法对20钢管在GI-920炸药滑移爆轰驱动下层裂的同题进行二维数值模拟。分析了一维内爆和滑移爆轰两 种加载方式下作用于钢管外表面的压力及钢管内部受力状态的区别，考察了滑移爆轰加载方式下钢管外表 面的受力随炸药厚度变化的规律．进而研究了钢管内损伤的分布和演化．以及裂纹的产生和扩展现象。计算 得到的层裂片初始厚度随炸药厚度的变化规律与实验结果符合较好． 关键词：爆炸力学；层裂f数值模拟；滑移爆轰}损伤 中图分类号：0381；0346国标学科代码：130·35 文献标志码：A 1 引 言 ． 动态断裂在理论和工程上都是重要而复杂的研究课题，其机理和特征因材料性质和加载方式的不 同而差异很大。滑移爆轰驱动金属圆管的层裂是典型的动态断裂现象，实验已观测到20钢管在GI一920 (固态920)炸药滑移爆轰驱动下的层裂图像，并给出了层裂厚度与炸药厚度的定量关系n咕3；文献[3—4] 在一维内爆假定下对该爆轰驱动圆管的层裂问题进行了数值模拟。 从爆轰加载条件来看，一维内爆时存在内聚效应，而滑移爆轰时存在侧向稀疏波的影响，因此同样 厚度的炸药在一维内爆时作用于圆管外表面的峰值压力通常比滑移爆轰高一些。同时，一维内爆情况 下压力比滑移爆轰时卸载速度快。这些因素都会影响炸药驱动金属圆管的能力。 从金属圆管的受力状态来看，在一维内爆下金属柱壳内各点在rz平面上只有正应力口，、％；而在 滑移爆轰下，还存在剪切应力分量％。同时，在一维内爆加载条件下仅需关注裂纹的产生，而在滑移爆 轰加载下，初始裂纹产生后沿圆柱壳轴向的扩展也是动态断裂过程的重要特征． 因此，有必要对滑移爆轰驱动金属圆管的层裂问题进行二维数值模拟，避免一维假定带来的误差。 本文中将在对二维动力学问题，特别是断裂问题数值模拟算法研究的基础上，对20钢管在GI-920炸药 滑移爆轰驱动下层裂的问题进行数值模拟，并与典型实验[21、一维数值模拟结果‘叼进行比较。 2数值方法 主体计算采用总能量守恒的二维拉氏算法，具体的计算格式这里不作介绍，可参见文献[5]。下面 依次介绍在爆轰、损伤、断裂处理方面的算法。 2．1爆轰模型 采用体积起爆方式，未反应炸药与爆轰产物的混合单元的平均压力取为p=即，，其中炸药反应份 为采用JWL状态方程得到的爆轰产物压力。 -收稿日期：2006—10-27I修回日期：2007一01—24 基金项目：国家自然科学基金委员会与中国工程物理研究院联合基金(NSAF)项目；国家重点基础研 究专项基金项目(2005CB321703) 作者简介：熊俊(1975一)，男。博士，副研究员． 万方数据 爆 炸 与 冲 击 第28卷 2．2损伤模型 在滑移爆轰驱动的层裂问题中钢管处于二维应力状态，因此实际的损伤状态并非严格的各向同性， 但是由于钢管所受应力张量中偏量相对较小，损伤演化仍然主要决定于压力，因此本文中仍然采用 Johnson(1981年)提出的适用于延性金属孔洞各向同性生长的VG(voidgrowth)模型[e3来描述断裂 前钢管的损伤演化。令多孔介质的单元宏观体积为V，孔洞体积为U，基体材料的体积K=y—U。 令a=y／K，其时间演化方程为 fO △夕≥O 占2{一!曼尝(口一1)V3△户△夕o