Page 11 - 《爆炸与冲击》2026年第5期
P. 11

第 46 卷       黄    阳,等: 基于物理信息及数据融合驱动的复杂街区爆炸荷载快速计算方法                             第 5 期

               道三维张量,因此,其输入为三通道三维张量。网络输出为单通道三维张量,作为目标时刻的非起爆街
               道压力场预测结果。在后处理部分,基于                  3D-UNet 输出的三维场景压力场,可获取指定时刻的压力分布
               及指定测点的压力时程。


                2D-UNet                                                                        Pre-processing
                                                                     Boundary
                                                                      pressure              Rigid wall
                    [16,16,16]  [16,16,16]  [16,32,32]  [16,32,32]  [16,32,32]  [16,32,32]  [48,32,32]  [48,32,32]  [8,32,32]  [8,32,32]  [1,32,32]  [1,32,32]  [1,21,21]  [1,21,21]

                          [32,16,16]  [32,16,16]  [32,16,16]  [32,16,16]  [96,16,16]  [96,16,16]  [32,16,16]  [32,16,16]  [16,16,16]  [16,16,16]  [16,16,16] [16,32,32] [16,16,16] [16,32,32]  L 1  z

                              [64,8,8]  [64,8,8]  [64,8,8]  [64,8,8]  [192,8,8]  [192,8,8]  [64,8,8]  [64,8,8]  [32,8,8]  [32,8,8]  [32,8,8]  [32,8,8]  Concat  d 2  y  h Ground
                                                                                                 x
                                              2D Conv+ReLu
                                              2D Max pool
                                              2D Up-conv
                                    [128,4,4]  [128,4,4]  [64,4,4]  [64,4,4]  [64,4,4]  [64,4,4]  2D Conv  Location-channel  d 1/2  d 1/2  L 2
                                                                         ECF-channel
                                                  0.16                     Deep  Pressure/  Numerical  Pressure/  Relative
                              3D-UNet                 The proposed method  learning  MPa  model  MPa  error  Error
                                                                                                      1.0
                                                                                           0.14
                                                                                 0.14
                                                      Numerical model
                                                  0.14                           0.12      0.12       0.8
                          …
                                                  0.12
                   [3,16,16,16]  [3,32,32,32]  [8,32,32,32]  [1,16,16,16]  Pressure/MPa  0.10  0.10  0.10  0.6
                                                                                                      0.4
                                                                                           0.08
                                                                                 0.08
                                                                                                      0.2
                                                  0.08
                                                                                           0.06
                                                           4
                                                       2
                                                    0
                                                           Time/ms 6  8  10 y    0.06                 0
                                                         Pressure histories  x       Pressure contour
                Calculation model                                     O                     Post-processing
                                              图 5    非起爆街道预测模型的推理过程
                                    Fig. 5    Inference procee of the non-detonation street prediction model
                1.3.1    物理机制引入
                   非起爆街道模型基于在           3D-UNet 网络数据输入端的能量密度因子(energy concentration factor, ECF)
               引入物理机制。能量密度因子是一种与炸药当量无关的物理量,用于量化能量在空间中的近似分布。
               由于镜像爆源法难以处理非点源爆炸工况,且                     Silvestrini 等  [29]  和钟巍等  [30]  的研究指出,在半约束空间中
               爆炸产生的冲击波超压与场景的能量密度因子表现出显著相关性,因此本研究采用能量密度因子作为
               非起爆街道模型的输入,以降低网络的数据映射复杂度。
                   在指定炸点位置及冲击波传播半径后,该半径处的能量密度因子                             f  F  的计算公式为:
                                                                             EC
                                                              V sph
                                                      f ECF =                                           (5)
                                                           V semi-confined
               式中:V  sp h  为自由场爆炸时冲击波围成的空间体积,V               semi-confine d  为冲击波在受限空间内围成的体积。针对
               各典型街道的      f EC F  计算如图  6  所示。其中,f  EC  F  计算时依据的炸点中心被设置在边界条件处的地面中心,
               淡绿色球体代表自由场爆炸下传播到当前位置时的爆炸波体积。对                                f  F  的计算过程可基于三维体积布
                                                                            EC
               尔运算实现。







                             r                             r                             r


                         (a) L-intersection            (b) T-intersection             (c) Crossroads

                                            图 6    典型街道能量密度因子     f EC 计算示意图
                                                                    F
                          Fig. 6    Schematic diagrams of energy concentration factor calculation for typical street configurations



                                                         051411-6
   6   7   8   9   10   11   12   13   14   15   16