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

第 46 卷           贾永胜,等: 水下钻孔爆破孔底复合垫层装药结构的损伤控制模型                                第 11 期


                            Image input
                                            ①
                      Conversion to bipartite graphs

                              k=1

                           Min(w,h)/k>1                                            10
                                                                                   9      y=−1.347 1x+9.856 2
                                                                                          R 2 =0.99
                                                                                   8
                        Divide the image into                      r=4
                        2 k ×2 k  point subregions  Original image                 7
                                                                                  lgN  6
                   Count the number of grids containing                            5
                         pixel value 0 as N(r)  ②                                  4
                                                                                   3
                    Add ln(1/2 k ) with N(r) to the X, Y                            0   1  2   3   4  5
                                                                                             lgr
                             k=2k                                  r=10
                                                   Image binarization                        ③
                          Linear fit to X,Y:
                          Y=kX+b, D=−k                  ①               ②
                                               ③
                         Outputting the fractal
                           dimension D

                                             图 10    岩样爆生裂纹分形维数的计算流程
                             Fig. 10    Calculation procedure of fractal dimension for blast-induced cracks in rock samples


                        10                                       10

                                            y=−1.347 1x+9.856 2                       y=−1.195 2x+9.346 6
                         9                  R 2 =0.99             9                   R 2 =0.99
                         8                                        8
                         7                                        7
                      lgN                                       lgN
                         6                     (S1)               6                     (S2)

                         5                                        5
                         4                                        4
                         3                                        3
                          0     1     2     3     4     5          0     1     2     3     4     5
                                          lgr                                      lgr
                                        (a) S1                                    (b) S2
                        10                                       10
                                            y=−1.174 2x+8.316 1                       y=−1.142 1x+8.406 3
                         9                  R 2 =0.98             9                   R 2 =0.99
                         8                                        8
                         7                                        7
                      lgN                                       lgN
                         6                     (S3)               6                     (S4)
                         5                                        5
                         4                                        4
                         3                                        3
                          0     1     2     3     4     5          0     1     2     3     4     5
                                          lgr                                      lgr
                                        (c) S3                                    (d) S4
                                          图 11    岩样爆破顶面的爆生裂纹盒维数拟合曲线
                                Fig. 11    Fitting curves of box dimension for blast-induced cracks on the top surface
                                                   of rock sample after blasting


                                                         113902-12
   136   137   138   139   140   141   142   143   144   145   146