飞片厚度对冲击波压力峰值衰减特性的影响分析

飞片厚度对冲击波压力峰值衰减特性的影响分析赵传荣;孔德仁;王胜强;杨凡【摘要】Incurrentstudies,theinfluenceofflyer'sthicknessontheattenuationcharacteristicsofshockwave'speakpressureisrarelyconsidered.Aimingatthisproblem,asimulationmodelwasbuilttosimulatetheflyersofdifferentthicknesshittinganoxygenfreecoppertargetboardandwassolvedbythenumericalmethod.Thesimulatedvaluesofpulsewidthattheinterfacebetweenflyerandtargetplatewerecomparedwiththetheoreticalvaluesandthecredibilityofthenumericalsimulationmodelaswellasthesimulationresultswasverified.Processingthesimulationdatawithleastsquaresmethodandestablishingthequantitativemathematicalrelationshipbetweentheflyer'sthicknessandtheattenuationcoefficientoftheexponentialdecaymodelofshockwave'speakpressure,itisshownthattheinfluenceofflyer'sthicknessontheattenuationcharacteristicsofshockwave'speakpressureisobvious,andtherelationshipisapproximatelylinear.Thesmallertheflyer'sthicknessis,thelargertheattenuationcoefficientaswellasthefastertheattenuationrateofshockwave'spressurepeakwillbe.Theresultsprovideatheoreticalreferencetothedesignandanalysisofrelatedexperiments.%针对现有研究没有考虑飞片厚度对冲击波压力峰值衰减特性影响的问题，采用数值仿真的方法，通过建立不同厚度的飞片撞击无氧铜靶板的仿真模型，比较飞片与靶板撞击面脉冲宽度的仿真值和理论计算值，验证了数值仿真模型和仿真结果的可信度。采用最小二乘法对仿真数据进行处理，建立了飞片厚度与冲击波压力峰值指数衰减模型衰减系数之间的定量数学关系式，结果表明飞片厚度对冲击波压力峰值衰减特性的影响比较明显，飞片厚度与衰减系数近似成线性关系，飞片厚度越小，衰减系数越大，冲击波压力峰值的衰减速率越快，为相关实验设计与分析提供了理论参考。【期刊名称】《振动与冲击》【年(卷)，期】2016(000)003【总页数】4页(P135-138)【关键词】冲击波;压力峰值;数值仿真;飞片厚度;衰减特性【作者】赵传荣;孔德仁;王胜强;杨凡【作者单位】南京理工大学机械工程学院，南京210094;南京理工大学机械工程学院，南京210094;西安近代化学研究所，西安710065浦京理工大学机械工程学院，南京210094【正文语种】中文【中图分类】O385平面冲击波特性是高能密度材料的合成［1］、炸药冲击起爆［2］及材料本构特性的获取［3］等研究中的重要课题。高速飞片撞击靶板是产生平面冲击波的有效方法之一。冲击波压力峰值和脉冲宽度是表征平面冲击波强度的两个重要参数。Erkman等［4］研究了平面冲击波在2024-T351铝中的衰减特性，结果表明，右行追赶稀疏波在没有追上靶板中右行冲击波之前，平面冲击波在靶板中的衰减主要表现为冲击波脉冲宽度的衰减；当追赶稀疏波追上右行冲击波以后，冲击波压力峰值将迅速衰减。汤文辉等［5］应用锰铜压阻计测量了峰压为4.92GPa的冲击波压力峰值在LY12-M铝中的衰减规律，飞片的厚度为2mm，材料为黄铜。程和法等［6］采用轻气炮装置驱动厚度为4mm、直径为55mm的飞片分别以101m/s和379m/s的速度撞击靶板，研究了不同冲击速度下，冲击波压力峰值在泡沫铝材料中的衰减特性，获得了相应的指数衰减模型。王永刚等［7］基于轻气炮装置驱动厚度为10mm的LY12铝飞片以670m/s的速度撞击厚度为5mm的LY12铝靶板，研究了冲击波压力峰值在C30混凝土中的衰减特性。然而，这些衰减模型或特性均是在特定的冲击条件下获得的，没有考虑飞片厚度的不同对冲击波压力峰值衰减特性的影响，目前也没有这方面的相关报道。由于现有的测量超高压力用传感器，如锰铜传感器或PVDF传感器，在10GPa左右的量程下，均是采用高速飞片撞击靶板，将测得的飞片撞击靶板的速度经流体力学模型换算得到的压力作为标准压力，并与传感器的输出电压或电荷量进行比较，从而实现对传感器的动态标定，获取其动态特性参数，如压阻系数或压电系数。因此，本论文采用数值仿真方法探讨飞片厚度对冲击波压力峰值衰减特性的影响，将数值仿真［8-9］结果与理论计算结果进行比较，验证仿真模型和结果的可信度。根据飞片与靶板撞击面冲击波压力的流体弹塑性模型［5］，将飞片撞击靶板过程看作为一维平面平行正碰撞，假设飞片直径远远大于飞片的厚度，则靶板撞击面中心处冲击波的压力峰值p0和脉冲宽度T0的理论计算模型为：式中，L是飞片的厚度，p01为飞片材料在零压下的密度，US1为飞片中冲击波速度，UP1为飞片中的粒子速度，C1为飞片材料在高压下的弹性声速，p1为飞片材料在P0下的密度。当追赶稀疏波以C1的速度到达飞片与靶板撞击面时，将以C2的速度继续向靶板中传播，这个波所到之处将会使该处的冲击波压力迅速衰减。故在靶板中传播距离X处冲击波脉冲宽度TX即是追赶稀疏波到达该处的时间与右行冲击波达到该处的时间之差，即：式中：p02为靶板材料在零压下的密度，p2为靶板材料在P0下的密度，C2为靶板材料在高压下的弹性声速，US2为靶板中冲击波速度。根据式(1)和式(2)，可将tx的计算模型简化为：式(3)中。令tx=0,则冲击波压力峰值保持不变的传播距离x1=t0/K。可见，飞片厚度L不同，脉冲宽度t0不同，会直接影响冲击波压力保持不变的传播距离x1,这势必将影响冲击波压力峰值的衰减。2.1飞片撞击靶板仿真模型的建立基于ANSYS/LS-DYNA软件分别建立厚度为2cm、3cm及4cm、直径为120cm的圆柱形平头飞片，撞击直径为120cm、厚度为60cm的靶板，飞片与靶板材料均为无氧铜。为了研究在10GPa峰值压力下无氧铜材料中冲击波的衰减特性，飞片的速度取520m/s。为了较好的反映靶板在撞击条件下的应力变化，飞片与靶板均采用Johnson-Cook材料模型和Gruneisen状态方程，无氧铜材料的Johnson-Cook本构参数和Gruneisen状态方程参数分别如表1和表2所列。采用Lagrange算法，用Solidl64单元进行划分，采用面面侵蚀接触算法，并在飞片与靶板对称面上施加对称边界约束，靶板采用无反射边界约束。考虑到飞片与靶板的结构具有对称性，为了减少计算规模和计算时间，只取1/4模型进行计算，数值模型单位采用g-cm-us单位制。2.2仿真结果可信度分析为了验证数值仿真模型和仿真结果的可信度，通过适当调整全局弹性刚度、网格尺寸等仿真参数，采用ANSYS/LS-DYNA专用后处理软件分别提取2cm、3cm及4cm厚飞片撞击靶板后，靶板撞击面中心处单元压力-时间历程曲线如图1所示。平台峰值压力均为10GPa，脉冲宽度分别为8.00卜、12.01ps及16.10ps,脉冲宽度是按照95%压力峰值的标准获取的。脉冲宽度的仿真值与理论计算值［5］的比较结果如表3所示。由图1和表3可见，脉冲宽度的仿真值与理论计算结果较吻合，其最大相对误差优于2%，二者随飞片厚度的变化趋势也一致，验证了数值仿真模型和仿真结果的可信度。2.3仿真结果分析为了研究飞片厚度对冲击波压力峰值衰减特性的影响，提取不同飞片厚度L下距靶板撞击面x处的冲击波压力峰值如表4所示。假设Px为距靶板撞击面x处冲击波压力峰值，则冲击波压力峰值随传播距离的衰减可表示为：式中，P0为靶板撞击面冲击波压力峰值；a为衰减系数，它体现了冲击波压力峰值的衰减速率；x1为冲击波压力峰值开始发生衰减的冲击波传播距离。为了综合评价在上述不同的冲击条件下冲击波压力峰值在无氧铜靶板中的衰减特性，将表4中的数据作归一化处理，可获得指数拟合关系曲线如图2所示。仿真数据的拟合结果列于表5中。根据飞片厚度L和衰减系数a的一一对应关系，采用最小二乘法对表5中的数据进行处理，可以得到其散点图和拟合曲线图，如图3所示。拟合曲线的函数关系式为：式中，L的单位为cm,a的单位为1,相关系数为0.9914。由图2和图3可见，飞片厚度不同，在相同的撞击速度下，无氧铜中冲击波压力峰值衰减特性的差异较大，飞片厚度与衰减系数近似成线性关系。当飞片速度及飞片与靶板的材料一定时，飞片厚度越小，衰减系数越大，冲击波压力峰值在靶板中衰减速率越快。采用数值仿真的方法分析了飞片厚度对冲击波压力峰值衰减特性的影响，结果表明飞片厚度对冲击波压力峰值衰减特性的影响比较明显，飞片厚度与衰减系数近似成线性关系，飞片厚度越小，衰减系数越大，对应的冲击波压力峰值衰减速率越快，为相关实验设计与分析提供了理论参考。【相关文献】FibbianiFP,PulhamCR.High-pressurestudiesofpharmaceuticalcompoundsandenergeticmaterials[J].ChemSocRev,2006,35(10):932-942.KimM,YooCS.HighlyrepulsiveinteractioninnovelinclusionD2-N2compoundathighpressure:Ramanandx-rayevidence[J].J.Chem.Phys,2011,134(4):44-51.PavlenkoAV,BalabinSI,KozelkovOE,etal.Aone-stagelightgasgunforstudyingdynamicpropertiesofstructuralmaterialsinarangeupto40GPa[J].InstrumentsandExperimentalTechniques,2013,56⑷：482-484.ErkmanJO,ChristensenAB.AttenuationofShockWavesinAluminum[J].JournalofAppliedPhysics,1967,38(13):5395-5403.汤文辉，张若棋,陈雪芳.LY12-M铝中冲击波衰减机理的实验研究[J].高压物理学报,1988,2(3):218-226.TANGWen-hui,ZHANGRuo-qi,CHENXue-fang.ExperimentalstudiesontheattenuationofshockwavesinLY12-Maluminum[J].ChineseJournalofHighPressurePhysics,1988,2(3):218-226.程和法，黄笑梅,薛国宪，等.冲击波在泡沫铝中的传播和衰减特性[J].材料科学与工程学报,2004,22(1):78-81.CHENGHe-fa,HUANGXiao-mei,XUEGuo-xian,etal.Propagationandattenuationcharacteristicofshockwaveinaluminumfoam[J].JournalofMaterialsScienceandEngineering,2004,22(1):78-81.王永刚，王礼立.平板撞击下C30混凝土中冲击波的传播特性[J].爆炸与冲击,2010,30(2):119-124.WANGYong-gang,WANGLili.ShockwavepropagationcharacteristicsinC30concreteunderplateimpactloading[J].ExplosionandShockWaves,2010,30(2):119-124.李顺波，东兆星,齐燕军，等.爆炸冲击波在不同介质中传播衰减规律的数值模拟[J].振动与冲击,2009,28(7):115-117.LIShun-bo,DONGZhao-xing,QIYan-jun,etal.Numericalsimulationforspreaddecayofblastingshockwaveindifferentmedia[J].JournalofVibrationandShock,2009,28(7):115-117.唐廷,朱锡,侯海量，等.高速破片在防雷舱结构中引起的冲击荷载的理论研究[J].振动与冲击,2013,32(6):132-136.TANGTing,ZHUXi,HOUHai-liang,etal.The