[發明專利]一種用于計算網狀坑道內爆炸沖擊波傳播的方法

有效

• 申請號: 202011636189.5
• 申請日: 2020-12-31
• 公開（公告）號: CN113158404B
• 公開（公告）日: 2021-11-30
• 發明（設計）人: 鄧燁;陶西貴;趙健;奉祁林;陳一村;賈藝凡;樊成飛;李炳臻
• 申請（專利權）人: 中國人民解放軍軍事科學院國防工程研究院
• 主分類號: G06F30/20
• 分類號: G06F30/20
• 代理公司: 北京翔石知識產權代理事務所(普通合伙) 11816
• 代理人: 李勇

1.一種用于計算網狀坑道內爆炸沖擊波傳播的方法，其特征在于，包括：

步驟一：構建網狀坑道的坑道網絡圖；

步驟二：根據爆炸點位置對所述坑道網絡圖更新，進行初始節點參數初始化，確定爆炸沖擊波入射超壓峰值初始計算節點；

步驟三：確定所述坑道網絡圖中的任意一個未計算爆炸沖擊波入射超壓峰值的坑道網絡內部節點作為下一個計算節點，并轉入步驟四；

步驟四：預設表征衰減因素的全局靜態指標，根據所述全局靜態指標，均按照網絡最短路方法計算由所述初始計算節點到達所選定內部節點的沖擊波最不利傳播路徑；

所述全局靜態指標包括指標一、指標二以及指標三，所述指標一為傳播累計距離，所述指標二為傳播累計節點度，所述指標三為傳播累計體積；

步驟五：根據所述步驟三中得出的三條最不利傳播路徑以及步驟二中更新的坑道網絡圖，計算三條最不利傳播路徑各自對應的綜合衰減系數，并確定最大綜合衰減系數和全局最不利傳播路徑，根據所述步驟二中的爆炸沖擊波入射超壓峰值初始計算節點和所述全局最不利傳播路徑，計算坑道網絡節點的爆炸沖擊波全局最不利入射超壓峰值；

步驟六：判定是否還有內部節點未計算爆炸沖擊波入射超壓峰值，若有節點未計算爆炸沖擊波入射超壓峰值，則返回步驟三，重新確定計算節點；若沒有內部節點未計算爆炸沖擊波入射超壓峰值，則進行步驟七；

步驟七：計算坑道段(v_i,v_j)中任意位置處離散點v_x的沖擊波超壓峰值P_x：

其中，為坑道段(v_i,v_j)在v_i點的全局最大入射沖擊波超壓峰值；為坑道段(v_i,v_j)在v_j點的全局最大入射沖擊波超壓峰值；|v_iv_j|表示坑道段(v_i,v_j)長度，|v_iv_x|表示坑道段(v_i,v_x)的長度；

所述步驟五中，計算所述綜合衰減系數時，根據爆炸點所處的坑道位置，分別采用相應的綜合衰減系數計算方法，其過程包括：

若爆炸點為坑道口部外點時，所述綜合衰減系數為，

其中：n_i表示第i種傳播衰減在某次沖擊波傳播路徑R中重復出現的次數，i＝1,2,3,4,5,n_i為整數；ξ₁為口部傳入衰減系數；ξ₂為直線傳播衰減系數；ξ₃為拐點傳播衰減系數；ξ₄為分叉傳播衰減系數；ξ₅為斷面變化衰減系數；

若爆炸點位于坑道內部時，則不存在口部傳入衰減情況，所述綜合衰減系

數為，

其中：n_i表示第i種傳播衰減在某次沖擊波傳播路徑中重復出現的次數，i＝2,3,4,5,n_i為整數；ξ₂為直線傳播衰減系數；ξ₃為拐點傳播衰減系數；ξ₄為分叉傳播衰減系數；ξ₅為斷面變化衰減系數