聚能射流失稳断裂.pdfVIP

聚能射流失稳断裂 秦承森 段庆生 林忠 李勇 (北京应用曲理与计算数学研究所．北京l呻O船) 摘要从哈密顿原理出发，给出描述射流运动的近似方程．数值解显示了影响射流失稳断裂的 因素。解析研究得出射流断裂时间及相邻断裂颗粒速度差的近似公式-与实验结果基本一致． i关t词聚脆射流失稳，断裂时问，速度差∑ 引言 广 聚能射流具有强大的破甲侵彻能力．在军事和工业中有着重要的应用。_般射流的侵彻 深度随射流拉伸长度的增长而增加，但是，由于射流失稳断裂，中止了射流的拉伸过程，降 低了射流的侵彻能力。因此，研究聚能射流断裂失稳的机理．延长射流断裂时间，是改善射 流质量，提高射流侵彻能力的重要课题。 I聚能射流失稳的数值模拟 1．I二雏敷值模拟 我们使用FCM(自由网格法)二维拉氏程序，采用J妇n￡00k屈服应力本构关系和 断裂应变公式，数值计算了聚能射流失稳断裂的过程。 计算初态采用美国105成型装药射流形成后某一时刻的一段射流形状和速度分布实测 p 值．计算到100s之后．该射流断裂成19个颗粒，与实验结果恰好一致；断裂颗粒的平均 速度差，断裂时间与实验平均值基本一致，颗粒半径计算值略小于实验值。 压力、应力、应变、速度、密度、温度等物理量的二维分布分析表明．在射流靠近轴线 部分，一直到射流断裂前，均可当成一维单向拉忡处理。 1．2一雏数值模拟 使用二维程序标定“等效”本构关系和“等效”断裂应变。可以使一维计算的断裂时间 与二维计算一致。一维数值解表明，粘度、剪应力是致稳因素；表面张力、气动力、材料强 um量级的微射流，表面张力效应才有必要 度、速度梯度均可导致射流失稳；只有直径为10 考虑；而当射流速度大于10km／s时，气动力的失稳效应开始显现，大于20km／g时．气动力 的影响将有可能超过材料强度．这与射流材料特性及射流头部形状等因素有关：屈服应力增 加．断裂时间缩短；应变强化效应增加，也使射流断裂时间缩短，但它的影响比屈服应力小 一个量级；热软化效应增加对射流时间缩短的影响介于屈服应力和应变硬化效应之间。 2射流运动方程 假设射流的运动为轴对称的，射流物质是不可压缩的，射流同一横截面上的轴向应力口， 和轴向速度r仅与时间f及轴向坐标=有关．与径向r无关。 则由哈密顿变分原理，可以导出描述射流运动的一维近似方程 一倒。鼍争+去嘉I鲁瓮善卜杀(打：+；p。“咄r)+ 磊1i磊I丽p船岬血J+—石一丽If刊裔：u cu(1’ 未{丢杀I啬扛船+p血)+等署I}一筹：o 鲁=粤一l 、7(2) 棚 |t 其中，¨：为轴向位移，』射流横截面积，月射流半径，r表面张力系数，dk剪切应力，p射 流表面气动压力，P。：2一孚，p置z：F!当，∥黏度系数，P密度，^拉氏坐标。 在郑哲敏先生所做的近似条件下，上述小扰动方程就是郑哲敏先生给出的小扰动方程。 3射流断裂无量纲时间^近似公式 一维矧辨晰冗发现，“满足的来系式为 一；蔷(剥幅¨碗志嘞 ∽ 舯岫舳髓皴龋靛脯戴毛。惫捂，毛2譬-芦。靠·本构关 系孑=矾t两可根据物质模型确定。公式(3)在特殊条件下，可以自然演化为以往文献中的 某些经验公式形式。 4相鄣囊艟速度薹公式 粒无量纲速度差dldr7J鲁J的近似公式 矿*等等【¨吉砉]n如I 一山J 卿～ 其中，二≥：墨；[：娶堕，胄为断裂射流颗粒的半径。 』^点：n+jto 使用二维程序计算的理论值为J矿=n873；使用计算的置-。，‰及近似公式可以得到 2RR．两者+替格：沂． ^矿一n 5结束语 在给出