[发明专利]一种用于减小后坐力的炮口制退器超压计算方法

权利要求书

1.一种用于减小后坐力的炮口制退器超压计算方法，其特征在于，包括如下步骤：

步骤一、获取前一时间步波阵面距炮口的距离波阵面的速度分量和冲击波球内气体总质量m_s(t)、火药气体质量m_g(t)、压缩空气质量m_a(t)；

步骤二、计算当前时间步的波阵面位置

当前时间步冲击波球内火药气体质量

m_g(t+Δt)＝m_g(t)+GΔt

冲击波球内压缩空气质量

当前时间步冲击波球内气体总质量

其中，G为流向此冲击波球的火药气体流量，d_s为冲击波球直径，ρ_a为空气密度；

步骤三、计算冲击波球内气体的平均密度

冲击波球内气体的平均压力

根据冲击波球内气体的平均压力，按经验公式计算冲击波波阵面上的平均压力

p_spj＝C_sp(p_pj-p_a)+p_a

其中，C_sp为符合系数；

步骤四、计算冲击波波阵面传播的马赫数

其中，v＝v_r+C_svv_c cosα，C_sv为修正系数；

步骤五、计算超压

。

2.根据权利要求1所述的用于减小后坐力的炮口制退器超压计算方法，其特征在于，步骤四中，根据正激波前后压力关系，计算出冲击波波阵面传播的平均马赫数

冲击波波阵面传播的平均速度为

当前时间步冲击波球的动量为

M(t+Δt)＝M(t)+Gλa^*Δt

其中，a^*为炮口截面的音速；

由此可得冲击波球心的速度

合成波阵面速度为

其中，α为波阵面点与冲击波球心的连线与x轴的夹角；C_sv为修正系数。

3.根据权利要求2所述的用于减小后坐力的炮口制退器超压计算方法，其特征在于，修正系数C_sv＝1.8。

4.根据权利要求2所述的用于减小后坐力的炮口制退器超压计算方法，其特征在于，符合系数C_sp＝10。