《兵工学报》论文投稿模板.docx

某大口径榴弹炮弹带挤进过程数值模拟研究二号,中文:小标宋/英文、数字：TimeSNeWRomaIi,居中）张三“，李四二王一五2（小四，中文：楷体/英文、数字：TimesNewRoman,居中）（L南京理工大学机械工程学院江苏南京210094：2.南京理工大学能源与动力工程学院.江苏南京210094）（六号，中文：宋体/英文、数字:TimesNewRoman收稿日期；（六号,黑体，粗体）202X-XX-XX （六号，Times New Roman）*W（六号，黑体，粗体）国家自然科学基金项目（12345678、23456789）（六号，中文*宋体/英文、数字：Tiaes New Kcxnan）通慎作者邮箱：居中）摘要：（五号，黑体，粗体）为探索某大口径榴弹炮发射条件下弹带挤进过程的力学机理，建立了弹带挤进坡膛的有限元模型。通过数值模拟研究弹带的动态挤进过程,分析弹带变形及刻槽形成过程，计算得到弹带动态挤进阻力、挤进压力和弹丸运动规律,给出了最大挤进阻力值、挤进压力值及对应时刻弹丸速度值。研究结果表明，挤进过程中弹带材料经历塑性变形流动，发生剪切失效占主导的韧性断裂并形成刻槽，弹带挤进后具有明显的层状特征，其内部区域的塑性变形量很小，弹带绝热变形产生的热量不足以使弹带材料熔化。（五号，中文：仿宋/英文、数字：TimesNewRoman）关键词：（五号，黑体，粗体）弹带挤进；弹带变形；挤进阻力；挤进压力；有限元模拟中图分类号，（五号，黑体，粗体）TJ301（五号，TimesNewRoman）文献标志码：（五号，黑体，粗体）A（五号，TimeSNewRoman）NumericalResearchonRotatingBandEngravingProcessOfaLarge-CaliberH。WitZer（四号，TimeSNeWRoman,粗体，居中）ZHANGSan12,LISiK,WANGYiw/（五号，TimesNewRoman,居中）（1.SchoolofMechanicalEngineering,NanjingUniversityofScienceandTechnology,Nanjing210094,Jiangsu,China;2.SchoolofenergyandPowerEngineering,NanjingUniversityofScienceandTechnology,Nanjing210094,Jiangsu,China）（小五，TimesNewRoman，居中）Abstract:（五号，黑体，粗体）Afiniteelementmodeloftherotatingbandengravingintotheforcingconesectionofguntubeisestablishedtoexplorethemechanismofengravingprocessofrotatingbandofalarge邺caliberhowitzerunderlaunchconditions.Thedynamicengravingprocessoftherotatingbandisstudiedthroughnumericalsimulation,andthemaximumresistance,engravingpressureandprojectilevelocityatthecorrespondingtimeareobtained.Thedeformationandgrooveformationprocessesoftherotatingbandareanalyzed.Thedynamicengravingresistanceoftherotatingband,theengravingpressureandmovementofprojectilearealsocalculated.Thecalculatedresultsshowthattherotatingbandundergoesplasticdeformationandmaterialflow,andformsthegroovesonitduetoductilefracture,whereshearfailureisdominant.Therotatingbandhasalayeredfeatureafterengraving,andtheplasticstraininthebandissmall.Theheatgeneratedbyadiabaticdeformationoftherotatingbandisnotenoughtomeltit.（五号，TimesNewRoman）Keywords:（五号，黑体，粗体）engravingofrotatingband;deformationofrotatingband;engravingresistance;engravingpressure;finiteelementsimulation（五号，TimesNewRoman,除特有名词外全小写）0引言（小四号，黑体，粗体）火炮发射条件下的弹带挤进过程具有高瞬态、高温强冲击、高速摩擦、大变形等复杂特点，经典内弹道理论忽略挤进过程，认为挤进过程瞬时完成,直接将挤进压力作为弹丸开始运动的起动压力，而不考虑挤进时期火药气体压力和弹丸速度的变化过程。随着现代火炮向高初速、远射程、高射速和高精度方向发展，弹带挤进过程对火炮膛压、初速、射击精度、身管寿命以及安全性等性能指标产生的影响更显著，瞬时挤进假设理论难以适用，需要研究弹带的动态挤进过程。（五号，中文：宋体/英文、数字：TimesNewRoman）文献1-3等经典理论工作中，针对弹带挤进阻力的研究受限于计算手段，对于弹带大变形弹塑性力学及材料损伤失效的计算过于简化。文献4从弹带挤进的角度研究了火炮内弹道峰现象的机理及其影响，其挤进压力和挤进阻力的计算仍然采用经典理论公式，有一定局限性。文献采用有限元法和Abaqus软件进行了挤进过程大变形和弹丸旋转弹带磨损分析，将身管和弹丸作为刚体，假设身管内膛是光滑对称的，铜弹带采用弹塑性模型，考虑摩擦因素，计算得到了挤进过程弹带和身管上的接触力。由于采用身管内膛光滑假设和不考虑弹带损伤失效的材料模型，计算模型与物理对象有较大差异。文献6-7就弹带构造对火炮身管受力的影响进行了非线性有限元数值计算和实验分析，但是略去膛线,将身管简化成光滑内膛，并只取周向一部分研究，使其分析与真实挤进过程偏差较大。文献网对两种大口径榴弹炮弹丸铜质弹带挤进坡膛过程进行了有限元计算，得到了变形、应力及相互作用力结果，但采用的是基于轴对称单元的二维模型分析，同样存在对物理几何过度简化的问题。文献9针对弹丸在枪管挤进过程中带有过盈摩擦大变形接触问题，运用动态显式算法以及网格自适应技术，通过有限元方法数值模拟研究了铅芯弹头的挤进过程，分析了弹头上压痕的形成过程和材料的流动情况，但是其所用膛压曲线数据是瞬时挤进假设下的计算值，载荷的施加存在合理性问题。文献10-11引入弹带材料的初始损伤及累计损伤模型，对在两种不同结构坡膛条件下的弹带挤进过程进行显式非线性有限元计算，得到了挤进过程中弹丸及弹带的动力学响应，分析了坡膛结构变化对挤进冲击力及内弹道性能的影响。但是，这些研究工作对初始内弹道中关注的动态挤进阻力和挤进压力等问题涉及较少，对弹带材料大变形及断裂失效模式的分析仍需补充完善。为探索某大口径榴弹炮发射条件下弹带挤进过程的力学机理，本文运用有限元方法，建立其弹带挤进线膛身管坡膛的三维有限元网格模型，采用考虑弹塑性大变形及断裂失效的弹带材料模型，以实测弹底压力作为载荷数据，通过数值模拟研究挤进过程中弹带塑性变形流动，由应力三轴度和Lode角参数判断弹带的断裂失效模式，并分析刻槽形成机理，由弹带变形后的温度分布考查弹带材料是否熔融，最终计算获取弹带动态挤进阻力、挤进压力和弹丸运动规律。1弹带挤进过程及模型假设1.1 弹带挤进坡Bi的过程（号，中文：黑体/英文、数字：TimeSNeWRoman,粗体）该大口径榴弹炮射击时，首先将弹丸装填到炮膛的正确位置。为了更好地密闭膛内火药气体，提供较为稳定的弹丸起动压力和良好的内弹道性能，弹丸采用带凸台的双弹带结构，且主弹带留有一定弹带强制量。装填后，弹带与坡膛紧密接触，使药室处于密闭状态。火炮发射时，膛内气体压力逐渐上升，当达到起动压力时，弹丸开始运动，弹带产生塑性变形逐渐挤进膛线。弹带的变形阻力随弹带挤进坡膛的长度而增加，直至最大值。之后随着弹带变形量减小，阻力减小。当弹带及其延伸部分全部进入身管膛线部，挤进终了。1.2 模型假设为简化计算模型，作如下假设：1）不考虑卡膛过程，初始时弹带凸台与药室前的坡膛密切接触而定位，不计弹带的初始应力和变形。2）忽略身管、弹丸本体的变形，假设它们为刚体；忽略弹丸的动不平衡，将弹丸前端和尾部截去，以等效质量代替。3）不考虑身管的后坐运动；不计弹丸前端空气动力的影响，也不计重力。4）不计温度应力场;假设挤进过程中弹带材料变形为绝热过程；不考虑摩擦产生的热量。2弹带挤进过程的有限元模型2.1 有限元网格弹带、弹丸本体、身管坡膛段的网格采用八节点六面体单元。弹带是挤进成形的关键部位,应采用细化网格。弹丸本体和身管坡膛段是刚性体,其网格尺寸可以稍大0根据模型试算的收敛性情况确定网格的基本尺寸，建立网格模型如图1所示,共有1748736个单元,其中弹带部分有1156800个单元。（五号,中文:宋体/英文、数字：TimesNewRoman）（a）弹带与弹丸本体六号，中文：宋体/英文、数字：TimeSNewRoman）（a）Rotatingbandandprojectile（六号，TimesNewRoman）(b)身管坡膛段(b) Forcing cone section of the gun tube图1有限元网格模型（小五，中文：宋体/英文、数字:TimesNewRoman）Fig.lMeshoffiniteelementmodel（小五，TimesNewRoman）所建立的有限元模型运用LS-DYNA软件U21的拉格朗日算法显式求解，采用单点积分和基于刚性的沙漏控制，既有利于大变形计算稳定可靠“叫又节省计算机资源。2.2 材料模型弹带材料为H90黄铜，弹丸本体材料为弹钢，身管材料为炮钢。弹带在挤进过程中经历弹塑性大变形及损伤，最终发生局部化韧性断裂，涉及到弹带材料的应变硬化、应变率硬化和温度软化，故采用Johnson-Cook塑性及断裂失效模型。2.2.1 JOhnSon-COok塑性模型（五号，中文：宋体/英文、数字：TimeSNewRoman）Johnson-Cook塑性模型”可中，VOnMises屈服应力是塑性应变、应变率以及温度的函数，即=（A+）l+Cln（l-T"）,m-£=上(3)式中：金为等效塑性应变；J为等效塑性应变率；备为参考应变率；T为温度；7；为室温；7；为熔化温度；A、B、C、?和为常量。表1弹带材料模型参数(五号，仿宋/英文、数字：TinICSNewRoman)TableIMaterialparametersofrotatingband(小五，TimesNewRoman)参数数值参数数值日MPa11500()Di0.54V0.31D24.89(kgm3)885-3.03(Jkg,K,)376Di0.014Z0.9DS1.124/MPa206c(cmsl)0.394BJMPa505S1