脉冲压力作用下扁球壳屈曲的数值模拟研究.pdfVIP

第九届全国冲击动力学学术会议论文集 脉冲压力作用下扁球壳屈曲的数值模拟 陈刚陈忠富牛伟 91940l#．62】900 中国l程物理Ⅻ究院总体『程Ⅳ}究所．叫川绵m 的脉冲E力载荷作用T的动态Ⅶ曲过程进行T有限i教值模拟．并进一步丹析了加载脉冲持续时间和脉冲峰值等㈨ i嚣；?；；：罗：晶，数值模拟 1引言 无论在民用山面，还是在军用方面．薄壁球壳作为一种基本结构元件，广泛应用于各种运载工 具和压力容器，如航天器和深水丁程的压力容器等。随着现代科技的发展，球壳在各种荷载和边界 条件下的稳定性研究不断深入。zhong和Ru斟”研究了球壳在平头圆柱体冲击下的破坏行为．假定 球壳在冲击作用下形成一止多面体的窝陷，然后由能量平衡给出了中心残余挠度和弹体初始动能之 问的关系式。宁建国等“’’研究了球形薄壳在子弹撞山作用下的大变形动态过程和动力贯穿行为，给 山一个描述贯穿过程的壳体破坏模埠!。马春生等”5对薄壁扁球壳结构H不同撞击速度撞击刚性板的 的变形模式和动态响应进行了试验研究，研究了薄罐扁球壳带附加质量撞击刚性面的动态吸能特性。 上述研究怂针对平头弹体或刚性平面对先体的撞击研究，本文针对脉冲压力载荷作用r的扁球 壳结构响应．在材料力学性能女i试的基础上，运用Ls．DYNA对～种扁球壳结构的动态屈曲过程进 行了有限元数值模拟，jf迸一步分析了加载脉冲持续时间和脉冲峰值等冈素对结构变形的影响规律。 2计算模型 21结构与有限元网格 扁球壳材料为纯铝1060，结构经过模具冲压而成，球壳直径250㈣丁贞部球冠半径200㈣ 壁厚3mma根据结构的特点．提取模型中面，采用四边形壳单兀进行阿格离散，如图l所示。在壳 结构边缘采用固支边界条件。 图1扁球壳罔格模型 22材料参数 材料的力学性能参数是数值模拟的基础，为提高数值模拟的可靠性，加T了拉伸试件，通过开 展力学性能实验米获得材料特性数据。 拉伸实验，图3为一发试件的拉仲应力应变曲线，试什在整体约百分2五塑性变形时发生拉伸必稳； 凹2同时给出了拉伸断裂后的残余试件．试什的总体仲长较小．但断口处的断面收缩超过70％。这 种状况表明该材料具有较好的韧性，但没有明显的应变强化效麻。在数值模拟中采用实验获得的材 第九届全国冲击动力学学术会议论文集 料应力成变数据。对图3所示的材料应力应变曲线的塑性段采川J_c模型的第项进行拟合，得到 参数为A=80MPa．8=57蛳Pa，n=015。在数值模拟中即运用这些参数进行计算，没有考虑材料的应 变率效应和温升软化效麻。 崮2 1060锅板状拉忡试件 目3 1060纯目的＆山应变曲线 23载荷条件 数值模拟中，选取脉冲形式为单脉冲压力．施加于扁球壳结构的外表面。在空问上，考虑壳体 结构的结构特点，假设载荷幅值随位置的变化为高斯分布P=晶e“p(一2r2，督)．其中屹取为 75㈣找荷最大值在半球壳结构顶部中心部位；脉冲在时间上为半正弦分布，如图4所示。 0 n 0 n 0 逸缝 8)空问分布 b)时问分布 圈4脉冲载荷形式 3计算结果 31结构屈曲塌陷过程 仝一 第九届全国冲击动力学学术会议论文集 c1 t=4．7ms mt=5ms e1 t=5．5ms nt=6ms 曲产7ms h、t=8ms 图5脉冲载荷作用F扁球壳的屈曲过程 壳体在上述载荷作用下，顶部将向下变形。当峰值载荷大于壳体的屈曲载荷时，壳体将发生屈 曲变形，结构