高速撞击充气压力容前壁损伤数值模拟.pdfVIP

76 中国空间科学技术 2010年8月 ChineseScienceand 第4期 Space Technology 高速撞击充气压力容器前壁损伤数值模拟 庞宝君盖芳芳 管公顺 (哈尔滨工业大学空间碎片高速撞击研究中心，哈尔滨150080) 摘要 针对空间碎片超高速撞击充气压力容器前壁损伤问题，应用非线性动力学分析 软件AUTODYN采用拉格朗日方法对球形弹丸撞击球形压力容器前壁穿孔进行了数值模 拟研究。在建模过程中通过对容器壁内侧施加压力边界条件来模拟由于内充气体的作用在 容器壁内产生的应力场，并通过与试验结果的比较验证了数值模拟方法的有效性。在此基 础上针对容器的内充气体压力、球形弹丸直径及撞击速度对充气压力容器前壁穿孔的影响 进行了研究。结果表明：在一定的气体压力下，气体压力对压力容器前壁穿孔直径与穿孔 形态的影响可以忽略不计；而撞击速度及弹丸直径对穿孔直径及穿孔形态有着较大的影响， 当撞击速度大于3km／s时，撞击穿孔边缘开始有裂纹产生，并且穿孔直径与裂纹直径随着 弹丸直径及撞击速度的增加而增大。利用压力容器前壁穿孔的数值模拟结果进行计算可以 得出当容器受到撞击速度大于3km／s的弹丸撞击后比撞击速度不大于3km／s时更易发生破 坏。 关键词 空间碎片 超高速撞击 穿孔压力容器数值模拟 1 引言 航天器上的某些压力容器是航天器用来储存液体或含能高压气体的部件，是受空间碎片／微流 星体威胁最大的关键部件之一。空间碎片／微流星体对压力容器超高速撞击造成的损伤从简单的容 器穿孔、撕裂到整个容器发生裂纹失稳破裂，是一个复杂的过程，其危害是很大的[1-z]。压力容器 受空间碎片／微流星体超高速撞击后，在压力容器前壁产生了穿孔及一些裂纹；在气体压力的作用 下，当裂纹尖端的应力强度因子超过压力容器材料的断裂韧性时，则裂纹发生失稳扩展，导致容器 发生破坏[3-4]。也就是说，在一定的撞击条件下，容器的内充气体压力存在一个门槛值(临界压 力)，当气体压力超过了这个门槛值，则容器发生裂纹失稳破坏。这种容器在气体压力的作用下发 生的断裂破坏，为了分析的方便认为是一种准静态过程。由于在不同的撞击条件下，受撞击后的压 力容器前壁产生的穿孔及裂纹的尺寸不同[5]，导致容器发生破坏的临界压力不同，因此对在不同撞 击条件下的压力容器前壁的穿孔及裂纹尺寸进行研究是有必要的。本文应用非线性动力学分析软件 AUTODYN对球形弹丸撞击充气球形压力容器的前壁穿孔进行了数值模拟研究，考察了气体压 力、弹丸直径及撞击速度对穿孔的影响；并对在不同的撞击参数下的球形压力容器发生准静态断裂 的临界应力进行了计算。 收稿日期：2009—11—04。收修改稿日期：2010—02—12 2010年8月 中国空间科学技术 2压力容器前壁穿孔的数值模拟 2．1计算模型的建立 针对高速撞击压力容器前壁穿孔及裂纹问题，由于Lagrange方法能够精确地描述结构边界的运 过程中通过对压力容器壁内侧施加压力边界条件的方法来实现内充气体对容器壁的作用。 图1初始几何模型 在压力容器壁内侧施加的压力产生的应力波在容器壁内传播，应力波在自由边界反射后在容器 壁内产生了没有衰减的应力波震荡。为了得到静应力场，可以通过AUTODYN的静态阻尼模块来 完成容器壁内的应力波衰减¨]。在笛卡尔坐标系下，z，Y方向的速度衰减公式为： 主井1／2一(1—2翘d)主，r垅／+(1～7浓d)≯At—I 多井V2一(1—2兀Rd)多”V2+(1一枳d)岁”At”l (】) R。：器 J 一 式中衰减因子风是由最小模态周期丁和计算时间步长础确定，稽为时间步数。 况进行了数值模拟，并将容器壁内应力的理论结果与数值模拟结果、前壁穿孔直径的