爆轰计算JWL状态方程参数的不确定度

针对爆轰计算JWL状态方程参数,开展不确定度量化研究。首先从JWL状态方程性质出发,结合数学、物理分析确定JWL状态方程中主要参数R1 的概率分布范围,然后利用非嵌入多项式混沌展开方法对PBX9404炸药驱动飞层外界面速度进行了不确定度量化分析。给出的分析方法与分析流程对爆轰计算及其他工程数值模拟不确定度量化有一定参考作用。     

一种与爆轰参数封闭的JWL方程参数确定方法

利用数值模拟和理论计算分析了炸药JWL状态方程参数与爆轰参数封闭的重要性;获得了利用圆筒实验测试结果计算炸药爆轰产物绝热等熵指数和爆压的方法;建立了与爆轰参数封闭的JWL状态方程参数的确定方法,并依据公布的圆筒实验数据,应用于两种典型炸药JWL状态方程参数的确定,获得的参数与炸药爆轰参数严格封闭,数值模拟结果与实验结果一致性好,表明炸药JWL参数确定方法合理可靠。     

含铝炸药爆轰产物JWL状态方程研究

对ＲＯＴ ９０１含铝炸药进行了爆轰产物ＪＷＬ状态方程研究。对ＲＯＴ ９０１炸药进行了２５ｍｍ和５０ｍｍ两种装药直径的圆筒试验,测试了圆筒膨胀过程的Ｒ ｔ关系。由冲击Ｈｕｇｏｎｉｏｔ关系和参考的ＰＯＰ关系拟合了ＲＯＴ ９０１炸药的反应速率。对两种装药直径的圆筒试验进行了综合分析和二维数值模拟计算,综合评估了ＲＯＴ ９０１的作功能力和圆筒试验的特性。通过对两种装药直径圆筒试验圆筒膨胀过程的解析分析计算和二维数值模拟,确定得到了ＲＯＴ ９０１含铝炸药爆轰产物ＪＷＬ状态方程参数。     

基于自适应和投影Wiener混沌的圆筒实验不确定度量化

由于炸药爆轰现象的复杂性和人们对它的认知缺陷,其表征爆轰流体力学过程的物理数学模型具有较强的不确定性,要降低基于爆轰建模与模拟的数值结果做出决策的风险,量化和评估不确定输入对爆轰系统输出结果的影响尤为重要。本文中针对具有高维随机变量的爆轰问题的不确定度量化,使用自适应基函数的Wiener混沌方法、耦合旋转变换和投影方法,减少截断空间的长度。针对输入变量相关性,使用Rosenblatt变换使其相互独立。针对不符合标准正态分布的变量使用等概率原则,将它化为标准正态分布。最后,使用自主研发的具有完全知识产权的爆轰数值模拟软件LAD2D,研究了具有高维不确定参数的圆筒实验的不确定度量化,给出期望、标准差、置信区间等统计信息,所得问题与实验数据比对,从而确认了模型的有效性。     

RDX基PBX的做功能力及JWL状态方程参数确定

为研究RDX基PBX炸药的做功能力并确定其爆轰产物的JWL状态方程参数,对RDX基PBX炸药和TNT炸药进行?50 mm标准圆筒实验,获得了圆筒膨胀位移和速度的时程曲线,对比得出RDX基PBX炸药的做功能力明显高于TNT炸药;基于能量守恒对实验数据进行非线性拟合,得到2种炸药爆轰产物的JWL状态方程参数。TNT炸药的拟合参数和通过AUTODYN软件计算得到的结果符合较好;将采用上述方法得到的RDX基PBX炸药爆轰产物JWL状态方程参数用于数值模拟,计算结果与实验值吻合较好,符合数值模拟标定JWL状态方程参数的要求。     

基于贝叶斯方法的模态参数识别的不确定度分析

振动模态分析和模态参数识别是动态测试的一个重要研究方向。模态参数在模型的修正、响应的预测、系统的健康检测及控制等方面有着重要的作用。但动态测试的不确定度分析,尤其是模态参数的不确定度的研究还十分缺乏。本文主要基于贝叶斯方法,通过傅立叶变换(FFT)建立时域数据和频域数据之间的对应关系。根据共振频带内的多个模态的响应数据得到相对应的模态参数,优化得到模态参数的最佳估计值,评定模态参数识别的不确定度。在固支梁的模态实验中,加速度传感器采集环境激励中的振动数据,运用贝叶斯法进行处理得到模态参数的最佳估计值。在此基础上,通过模态参数的最佳估计值,以及仪器的检定报告数据,结合合成不确度分析方法,系统分析了模态参数识别的不确定度。     

非嵌入式多项式混沌法在爆轰产物JWL参数评估中的应用

介绍了非嵌入多项式混沌法的数学模型,给出了非嵌入式多项式混沌法进行不确定度量化的主要步骤。采用此方法研究了平面、散心爆轰问题数值模拟中, JWL模型参数R₁、R₂服从均匀分布的随机变量时所引起的爆轰过程计算结果的不确度性,着重分析了爆轰传播过程中压力与密度的统计特性。研究结果表明,非嵌入式多项式混沌法可以为模型输入参数不确定性的传播对输出结果响应量的影响建立一种有效不确定度评估方法,为使用JWL模型时选取参数提供参考。     

模态参数不确定性分析的贝叶斯方法研究

在结构损伤诊断和参数识别中,实测结构模态参数不可避免地存在误差。本文将模态参数视为随机变量,采用贝叶斯方法对模态参数的不确定性进行分析。分析中选用高斯联合概率密度函数作为先验密度函数,通过多次独立的模态参数测试,得到传递函数的条件概率密度函数和模态参数的后验估计表达式,再利用拉普拉斯渐近方法求解边缘概率密度函数,得到模态参数的最大后验估计。在钢筋混凝土框架结构的模态试验中,利用本文方法给出了结构模态参数的估计值,结果表明,本文方法具有良好的收敛性。     

基于高斯过程分类的结构贝叶斯可靠性分析

贝叶斯可靠性方法是处理不完备信息条件下结构可靠性问题的有效途径之一。在实际应用中,由于可靠性分析的计算量较大,常须采用各种近似替代模型以提高计算效率。传统的替代模型方法是对结构的功能函数予以近似建模。这种方法不易定量考虑模型误差对可靠性分析的影响,且难以应用于诸如功能函数不连续和失效域不连通等情况。为此,本文提出一种基于高斯过程分类的替代模型,直接辨识结构的极限状态曲面,并将其应用于结构贝叶斯可靠性分析之中。分析了替代模型不确定性对可靠性预测结果的影响,给出了失效概率分布参数的方差算式,进而提出了改善模型精度的补充采样准则。通过算例验证了方法的适用性和有被性．     

基于多项式混沌方法对C-J爆轰参数不确定度的分析

Chapman-Jougeut理论是预测波后爆轰物理量状态的有力工具,但以往的研究未考虑模型中的不确定因素及其影响。事实上,不确定度会影响数值模拟的预测能力和可靠性。首先,通过剖析爆轰机理,深入挖掘爆轰建模与模拟中的不确定因素。假设PBX-9502的初始密度和爆速服从对数正态分布,结合真实的试验数据,通过参数估计和Anderson-Darling假设检验法标定初始密度和爆速的概率密度函数。Beta分布用以定量刻画没有物理意义的、唯象参数的不确定度,形状参数和支集源于工程经验。Rosenblatt变换将相关的、非Gauss随机变量转化成相互独立的标准正态分布。然后,使用非嵌入多项式混沌研究高维爆轰不确定度传播。具体而言,针对一元多项式混沌,正交多项式通过Gauss-Hilbert空间中的Gram-Schmidt方法导出,六点Gauss求积方法用以计算多项式混沌的系数。使用权重和Gauss求积点的全张量积计算多元多项式混沌。最后,通过多元多项式混沌得到感兴趣量的概率密度函数以及对应的期望、标准差和置信区间等Gauss统计量。研究结果表明：波后压力波动较大,置信区间较宽,与孙承纬的“爆轰压力...     

RHT-902和Octol炸药爆轰产物JWL状态方程研究

本文对RHT-902和Octol两种炸药做了标准圆筒试验,运用能量守恒原理分析处理试验结果,用解析关系式确定了这两种炸药爆轰产物的JWL状态方程参数。这些参数与国外发表的相应参数基本一致。     

基于圆筒实验的RDX/Al炸药反应进程

对RDX炸药和2种铝粉质量分数分别为15%、30%的RDX基含铝炸药进行?50mm圆筒实验,研究铝粉含量对炸药做功能力的影响,根据格尼公式分析铝粉与爆轰产物的反应进程。结果表明:在圆筒实验记录的时间范围内,铝粉质量分数为15%的含铝炸药做功能力最强,RDX炸药次之,铝粉质量分数为30%炸药做功能力最弱;34μs时刻,铝粉质量分数为15%的炸药,铝粉的反应度为0.49,而铝粉质量分数为30%炸药铝粉的反应度仅为0.21,含铝炸药中铝粉的反应时间在50~200μs之间。     

爆轰产物状态方程的水下爆炸反演理论研究

由于水下爆炸过程中爆轰产物的信息以水中压缩波形式向外传递,本文旨在探讨如何利用水下爆炸试验数据确定爆轰产物的状态方程。相较于常规圆筒试验,水下爆炸试验具有装置简单成本低、装药尺寸限制少、测定压力范围更广的特点,更适用于大药量或非理想炸药的现场测试。本文从水中冲击波轨迹和波后压力时程曲线出发,发展了由冲击波及其波后流场还原水气界面的逆特征线算法,以及根据水气界面确定爆轰产物状态方程的遗传算法。与水下爆炸正演结果对比,发现逆特征线法可以较准确地还原水气界面的位置和压力参数,有效压力下限可达2 MPa,远低于圆筒试验的测试下限0.1GPa。根据水气界面的反演结果,遗传算法也能稳定地优化JWL方程参数,8种常用炸药的等熵衰减压力误差在爆压至0.01 GPa的区间内都小于3%。结果表明,利用本文的逆特征线算法和遗传算法,理论上可以反演出压力范围较宽、准确性较高的爆轰产物状态方程。     

混凝土高压状态方程实验与数值模拟研究

为研究高压下混凝土冲击压缩特性以及确定HJC本构模型状态方程参数,采用$\varnothing $58 mm火炮加载技术和多普勒探针系统（Doppler probe system,DPS）测速技术,对抗压强度为26.5 MPa（C25）和42.1 MPa（C40）的混凝土进行反向平板撞击实验研究与数值模拟。DPS探针记录TU1无氧铜靶自由面粒子速度历史,根据一维应变弹塑性波理论,计算撞击压力,拟合得到2～11 GPa高压条件下混凝土波速与粒子速度（u_s-u_p）、压力与体积应变（p-μ）关系。实验结果表明：高压条件下,混凝土波速-粒子速度呈线性关系;两种初始密度、孔隙率相近,强度不同的混凝土波速-粒子速度、压力-体积应变关系存在明显差异,相同压力下,混凝土试件强度越高,体积应变越小。基于实验结果,确定了两种强度混凝土HJC本构模型状态方程参数,利用LS-DYNA动力有限元分析软件对平板撞击实验进行了数值模拟,靶板自由面粒子速度历史与实验曲线吻合较好,仿真结果表明混凝土中冲击波的追赶卸载现象仅存在于低速撞击条件下。

     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。