论文部分内容阅读
随着国家科学技术的快速发展,炸药在国防、建筑和生活等领域的应用越来越广泛,爆炸冲击波作为炸药所引起的主要效应,其传播规律的研究显得极其重要。早在诺贝尔发明炸药的时候,人们就开始研究炸药所引起的一系列现象,爆炸冲击波作为动力荷载直接作用在结构物上,会引起结构物的剧烈振动和损伤,甚至超过结构物的极限承载能力发生而倒塌等非常严重的灾难。虽然很多学者对爆炸冲击波传播问题进行了深入的研究,但在爆炸初期爆炸产物和冲击波相互作用,使得冲击波传播变得非常复杂,这对爆炸冲击波的模拟研究带来了很大的困难。总的来说对爆炸冲击波的研究主要有理论分析、实验、数值仿真三种方法,人们一般利用实验及物理现象来研究其传播规律,但耗资巨大且受实验条件的限制,并且得到的实验数据非常有限。随着计算机的发展及有限元理论的普及,利用商业有限元软件来模拟爆炸冲击波的传播成为可能。现有的很多商业软件有LS DYNA﹑AUTODYN等,研究成本低且能够得到广泛的数据,但在计算的过程中这些软件引入了人工粘性技术来处理冲击波的强间断,由于人为抹平冲击波阵面,故使计算的压力峰值产生很大的误差,从而不能真实的反映爆炸冲击波的传播过程。FCT(简称通量修正输运法)方法是一种高阶的计算方法,它在计算连续性方程时,通过添加限制耗散通量的方法,很好的处理强间断区域的解析解问题且无数值振荡。本文在研究FCT方法原理的基础上,通过泰勒公式展开分析一维方程的耗散项,发现关键参数△t、△x、u之间一定的取值关系能够在原有的FCT方法基础上在空间上提高一阶计算精度。本文用方波的精确解与二阶中心差分方法、FCT方法、新方法进行对比分析,验证了新改进的方法提高了FCT方法的计算精度。然后简单介绍了程序中四个主要子程序的变量及功能,并对边界条件的模拟进行讨论。最后结合爆炸冲击波的控制方程,利用FCT方法对三个方程进行耦合求解实现爆炸冲击波的数值模拟。由于冲击波的强间断特性,一般的高阶算法会产生很大的振荡,所以在计算之前必须对程序的有效性进行验证。由于激波管具有强间断性和精确解,本文用三个不同的激波管算例验证了程序的适用性。根据爆炸冲击波基本原理,对程序的初始参数进行简单的讨论和计算。最后应用计算程序模拟爆炸冲击波的传播,分别从下面三个方面进行分析:(1)对空中自由场的爆炸现象进行模拟,将计算结果与LS DYNA软件计算结果进行对比分析,发现压力时程曲线吻合较好但压力峰值偏大;(2)对结构物在空中爆炸荷载作用下的动态响应进行模拟,将计算结果与实验进行对比,发现程序计算结果与实验测量结果吻合较好,证明程序的有效性。(3)对不同比距离球形炸药爆炸的压力峰值进行模拟,将计算结果与经验公式、LS DYNA软件进行对比分析,发现在近场区域程序计算结果比软件值大但更接近于经验公式预测值,在远场区域三者计算结果基本吻合,验证了程序的适用性。通过上述的研究表明本文的计算程序能够较好的模拟爆炸冲击波的传播,具有一定的应用研究价值。