1.一种计算爆破作用下隧道关键块体稳定性的方法，其特征在于，包括：S1：通过弹性波动力学及Snell定理对P波入射关键块体传播路径进行推导计算，获得P波入射关键块体造成块体节理的正应力值及剪切应力值；S2：构建基于弹性地基梁理论的喷射混凝土支护力学模型，通过所述喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力最大值q max ；S3：通过所述P波入射关键块体造成块体节理的正应力值及剪切应力值和所述喷射混凝土与围岩间的接触应力最大值q max ，计算获得关键块体的安全系数K；步骤S1具体为：入射J1节理的入射波P 1 表示为：P 1 (t-L 1 /c p )其中，t为应力波在关键块体内传播的时间，L 1 为应力波从出发点传播至J1节理所走的距离，C p 为纵波波速；入射波P 1 的应力波从J1节理传播至J2节理的距离表示为：L 2 ＝x i (tanα 1 +tanα 2 )其中，X i 为入射应力波的空间坐标，α 1 为J 1 节理与水平方向的夹角，α 2 为J2节理与水平方向的夹角；则入射J2节理的入射波P 1 可表示为：P 1 (t-L 1 /c p )(T P (t-L 2 /c p ))通过计算获得入射波P 1 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值，具体的计算公式为：其中，i为虚数单位，k n 为应力波波数，λ n 为拉梅常数，μ n 为岩体剪切模量， 为应力波振动方向矢量， 为应力波传播方向矢量，A n 为应力波振动幅值，η n 为k n (x·p (n) -c n t)，x为应力波空间坐标，c n 为应力波横波或纵波波速；重复上述步骤，获取入射波P 2 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值，和入射波P 3 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值；所述P波入射关键块体造成块体节理的正应力值及剪切应力值，为P 1 、P 2 和P 3 对应的正应力值及剪切应力值的叠加值；步骤S2中，所述通过所述喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力最大值q max ，具体为：喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力q表示为：q＝ky其中，k为围岩弹性抗力系数；y为弹性地基梁在重力作用下的位移；y可表示为：其中，弯矩M 0 与剪力Q 0 为初参数，α为弹性地基梁的弹性特征值，x为喷层梁与围岩地基接触面上任一点到荷载的距离；弹性地基梁的弹性特征值α的计算公式为：其中，E为喷射混凝土的弹性模量；h为喷射混凝土厚度；I为地基梁刚度；喷射混凝土与围岩间的接触应力q的最大值发生在x＝o处，此时M 0 ＝0，Q 0 ＝P，则所述喷射混凝土与围岩间的接触应力最大值q max 的计算公式为：其中，u m 为块体底部周长，W为关键块体重量；步骤S3具体为：计算喷射混凝土支护的关键块体受爆破振动作用的法向力F′ n 和剪切力F′ τ ，计算公式为：其中，σ 0 为围岩初始应力，σ J 为应力波造成的正应力，τ J 为应力波造成的剪应力，ρ为岩体的密度，s为楔形块体的面积，g为重力加速度，l 2 为J2节理的长度；根据极限平衡理论，所述关键块体的安全系数K的计算公式为：其中， 为节理内摩擦角。
2.一种计算爆破作用下隧道关键块体稳定性的系统，其特征在于，包括以下模块：应力计算模块，用于通过弹性波动力学及Snell定理对P波入射关键块体传播路径进行推导计算，获得P波入射关键块体造成块体节理的正应力值及剪切应力值；接触应力最大值计算模块，用于构建基于弹性地基梁理论的喷射混凝土支护力学模型，通过所述喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力最大值q max ；安全系数计算模块，用于通过所述P波入射关键块体造成块体节理的正应力值及剪切应力值和所述喷射混凝土与围岩间的接触应力最大值q max ，计算获得关键块体的安全系数K；应力计算模块的工作流程具体为：入射J1节理的入射波P 1 表示为：P 1 (t-L 1 /c p )其中，t为应力波在关键块体内传播的时间，L 1 为应力波从出发点传播至J1节理所走的距离，C p 为纵波波速；入射波P 1 的应力波从J1节理传播至J2节理的距离表示为：L 2 ＝x i (tanα 1 +tanα 2 )其中，X i 为入射应力波的空间坐标，α 1 为J 1 节理与水平方向的夹角，α 2 为J2节理与水平方向的夹角；则入射J2节理的入射波P 1 可表示为：P 1 (t-L 1 /c p )(T P (t-L 2 /c p ))通过计算获得入射波P 1 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值，具体的计算公式为：其中，i为虚数单位，k n 为应力波波数，λ n 为拉梅常数，μ n 为岩体剪切模量， 为应力波振动方向矢量， 为应力波传播方向矢量，A n 为应力波振动幅值，η n 为k n (x·p (n) -c n t)，x为应力波空间坐标，c n 为应力波横波或纵波波速；重复上述步骤，获取入射波P 2 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值，和入射波P 3 入射关键块体造成块体J1与J2节理的正应力值及剪切应力值；所述P波入射关键块体造成块体节理的正应力值及剪切应力值，为P 1 、P 2 和P 3 对应的正应力值及剪切应力值的叠加值；接触应力最大值计算模块中，所述通过所述喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力最大值q max ，具体为：喷射混凝土支护力学模型计算获得喷射混凝土与围岩间的接触应力q表示为：q＝ky其中，k为围岩弹性抗力系数；y为弹性地基梁在重力作用下的位移；y可表示为：其中，弯矩M 0 与剪力Q 0 为初参数，α为弹性地基梁的弹性特征值，x为喷层梁与围岩地基接触面上任一点到荷载的距离；弹性地基梁的弹性特征值α的计算公式为：其中，E为喷射混凝土的弹性模量；h为喷射混凝土厚度；I为地基梁刚度；喷射混凝土与围岩间的接触应力q的最大值发生在x＝o处，此时M 0 ＝0，Q 0 ＝P，则所述喷射混凝土与围岩间的接触应力最大值q max 的计算公式为：其中，u m 为块体底部周长，W为关键块体重量；安全系数计算模块的工作流程具体为：计算喷射混凝土支护的关键块体受爆破振动作用的法向力F′ n 和剪切力F′ τ ，计算公式为：其中，σ 0 为围岩初始应力，σ J 为应力波造成的正应力，τ J 为应力波造成的剪应力，ρ为岩体的密度，s为楔形块体的面积，g为重力加速度，l 2 为J2节理的长度；根据极限平衡理论，所述关键块体的安全系数K的计算公式为：其中， 为节理内摩擦角。
