Today underground structures are increasingly being built and exploited, due to the expansion and development of cities. These structures play a key role in the urban traort system. Considering that the excavation and construction of these structures is time consuming and expensive, and important equipment and facilities may be maintained in those, consideration of their sustainability issues, especially under dynamic loadings, is very important. Deep underground structures in many cases located in rock masses and the discontinuities of these rocks have a significant effect on the response of the system during blast loading. On the other hand, the correct estimate of the blast load is one of the most complex and sensitive issues in engineering sciences. Regarding these issues, in this research, the response of the urban shallow tunnels located in rocky masses of the southern part of the Isfahan City was studied using FLAC3D software. For the Correct modeling of the rock masses, the Discrete Fracture Network (DFN) Advanced method was used and the required data were extracted by scanline method from outcrops located at a project site in south of Isfahan which were processed using DIPS software. The Coupled Euler-Lagrange method was also used to estimate the blast load in ANSYS AUTODYN software. Considering the available experimental data and sensitivity analysis on the geomechanical properties, it was observed that increasing the rock mass strength properties and reducing the modulus of elasticity of the intact rock lead to increasing the PPV and residual displacement at the roof of tunnel. According to this sensitivity analysis, three types of M, S and C media were considered with different geomechanical properties and for each of these media a design experiment with two factor and three level were considered. The regression analysis showed that according to PPV criterion and Sakurai's critical strain criterion the explosive charge weight that need for destruction of tunnel in S Medium is 22% and 33% larger than explosive charge need in M medium, respectively and this values for C medium is 17% and 24% lower than explosive charge need in M medium, respectively. The design experiment using taguchi method and ANOVA analyzes showed that among the operating parameters, the depth of penetration of explosive charge and among the geomechanical parameters, the modulus of elasticity of intact rock had the highest percentage contribution and impact on the response of system. Finally, it was observed that an increase in the thickness of the tunnel lining from 30 to 70 cm in the medium M with an explosive charge weight of 300 kg and a depth of 8 m, lead to a 26.2% and 13.8% reduction in residual displacement and PPV, Respectively.