求助各位老師，大神，我做的是一個(gè)ls-dyna模擬隧道爆破，采用流固耦合計(jì)算。圍巖采用拉格朗日算法，炸藥 、洞內(nèi)外空氣以及ALE重復(fù)網(wǎng)格采用ale算法，三者采用多物質(zhì)耦合（ALE_MULTI-MATERIAL_GROUP）；圍巖和流體采用流固耦合（CONSTRAINED_LAGRANGE_IN_SOLID）；采用方形炮孔等效計(jì)算。2號(hào)巖石乳化炸藥量為360kg，

有以下問題：

1 沖擊波超壓數(shù)值（較實(shí)測）計(jì)算偏大，振動(dòng)速度偏小一個(gè)數(shù)量級(jí)；計(jì)算距掌子面50m處洞內(nèi)超壓峰值為40 kpa,邊墻振動(dòng)速度為 1.1cm/s;距掌子面100m處洞內(nèi)超壓峰值為 25 kpa,邊墻振動(dòng)速度為 0.75cm/s.

2 空氣沖擊波超壓值沒有隨著爆心距增加遞減 ?

圖1 模型網(wǎng)格圖

            圖2   t=0.05s時(shí)沖擊波超壓云圖

            圖3   t=0.1s時(shí)沖擊波超壓云圖

            圖4   t=0.25s時(shí)沖擊波超壓云圖

            圖5   t=0.5s時(shí)沖擊波超壓云圖

           圖6 距掌子面50m、100m、150m、200m處超壓值時(shí)程曲線

           圖7 距掌子面50m、100m、150m、200m處邊墻振速值時(shí)程曲線

k文件

請(qǐng)各位大神抽空多指點(diǎn)指點(diǎn)，不勝感激！實(shí)在折騰不出來 搞了一個(gè)月了。