CN117133388B 一种含能材料三项式点火增长模型参数优化计算方法 （北京理工大学）

号一种含能材料三项式点火增长模型参数优本发明公开了一种含能材料三项式点火增于混合粒子群算法和初始参数及其适应度函数2采用非线性有限元计算软件，建立冲击起爆含能材料二维轴对称计算采用三项式点火增长模型，描述被测含能材料，该模型由未反应含能材料JWL状态方;;R12和;G2步骤三，在三项式点火增长模型中，针对未反应含能材料JWL状态方程和反应速率方和试验结果的误差；在标定反应速率方程参数时，打开含能材料的所有反应项，设计变量定义为Figmax、FG1maxG2minG2，适应度函数定义为最高点起跳速度和起跳斜率计算结果与试验结果的误差；3其中，被测含能材料由通过楔形表面配合的两部分药柱构成，示波器记录信号；当爆轰波依次到达单U型电磁粒子速度计和组合型电磁粒子速度计所在型电磁粒子速度计和组合型电磁粒子速度计的水平部分工作段长度切割均匀磁场的磁感的情况下，通过示波器测量出电磁粒子速度计运动所产生的感应电动势ε,根据Faraε=B×L×v；通过改变聚四氟乙烯隔板的厚度，调节冲击波入射被测含能材料时的在三项式点火增长模型参数优化计算中，调节三项式点火增在待优化参数的取值范围内，自定义划分参数取值区间，利用评估待优化参量对适应度函数值影响：针对待优化参量每个取值区间4参数交叉变化，是指基于待优化参量对适应度函数值影响的参数变异变化，是指基于待优化参量对适应度函数值影响的将粒子群中每个粒子的适应度函数值与自身历代最优适应度函数值进5采用该计算方法获得的三项式点火增长模型参数计算的有限元计算结果与试验结果相吻[0006]一种含能材料三项式点火增长模型参数优化计算方法，该计算方法包括以下步6触发示波器记录信号；当爆轰波依次到达单U型电磁粒子速度计和组合型电磁粒子速度计单U型电磁粒子速度计和组合型电磁粒子速度计的水平部分工作段长度切割均匀磁场的磁变的情况下，通过示波器测量出电磁粒子速度计运动所产生的感应电动势ε,根据Faraday7设定反应速率方程中反应度λ的最大值与最小值，以控制不同反应项的开始与关闭；当λ>[0032]在标定未反应含能材料JWL状态方程参数时，令点火增长模型的参数Figmax＝0、该待优化参量的取值区间；对适应度函数值变化幅度小且增大的待优化参量的取值区间，8[0043]判断全局最优适应度函数值是否满足计算阈值要求小于15若满足，则输出三[0051]图3为本发明含能材料三项式点火增长模型参数优化计算方法的一种具体流程增长模型包含未反应含能材料JWL状态方程和反应速率方程，针对三项式点火增长模型的9未反应含能材料JWL状态方程在实施例一中进行详细描述，针对三项式点火增长模型的反[0060]步骤一S1，开展三项式点火增长模型中未反应含能材料JWL状态方程参数标定试粒子速度计11所在位置处，单U型电磁粒子速度计10和组合型电磁粒子速度计11随着波后爆轰产物一起运动，单U型电磁粒子速度计10和组合型电磁粒子速度计11的水平部分工作的水平部分工作段长度L不变的情况下，通过示波器12测量出电磁粒子速度计运动所产生的感应电动势ε,根据Faraday电磁感应定律：设定反应速率方程中反应度λ的最大值与最小值，以控制不同反应项的开始与关闭；当λ>[0079]在标定未反应含能材料JWL状态方程参数时，令点火增长模型的参数Figmax＝0、该待优化参量的取值区间；对适应度函数值变化幅度小且增大的待优化参量的取值区间，[0089]判断全局最优适应度函数值是否满足计算阈值要求小于15若满足，则输出三算阈值要求小于15输出未反应含能材料JWL状态方程最优参数。图5给出了未反应含能冲击起爆含能材料计算模型，定义未反应含能材料JWL状态方程待优化参数、参数取值范法和初始参数及其适应度函数值，进行未反应JWL状态方程参数交叉和变异变化，迭代计[0093]由于将基于遗传变异的混合粒子群算法应用到三项式点火增长参数的标定问题子速度计11所在位置后，单U型电磁粒子速度计10和组合型电磁粒子速度计11会随着波后爆轰产物一起运动，单U型电磁粒子速度计10和组合型电磁粒子速度计11的水平部分工作量，点火项是冲击波强度和压力持续时间的函数；G1和d控制了点火后热点早期的反应生该待优化参量的取值区间；对适应度函数值变化幅度小且增大的待优化参量的取值区间，代最优适应度函数值，则粒子的自身最优适应度函数值修改为当前粒子的适应度函数值；适应度函数值，直到全局最优适应度函数值满足计算阈值要求小于15输出含能材料反