1.一种围岩力学性质动态损伤本构模型的构建方法，其特征在于，所述围岩力学性质动态损伤本构模型的构建方法包括以下步骤：S1：获取围岩性质的各项参数；S2：根据所述围岩性质的各项参数，计算渗流作用下对围岩裂隙拓展的损伤程度和流-力耦合作用下围岩岩土损伤程度；S3：根据所述渗流作用下对围岩裂隙拓展的损伤程度和所述流-力耦合作用下围岩岩土损伤程度，得到流-力耦合场围岩损伤程度；S4：根据所述流-力耦合场围岩损伤程度对所述各项参数中的部分参数进行折减，得到折减后的部分参数；S5：根据所述折减后的部分参数和所述流-力耦合场围岩损伤程度，计算热-力耦合场损伤后的围岩损伤程度；S6：根据所述热-力耦合场损伤后的围岩损伤程度和所述述流-力耦合场围岩损伤程度，得到由于热对流引起的水-力-热多场耦合作用下围岩损伤程度；S7：根据所述由于热对流引起的水-力-热多场耦合作用下围岩损伤程度和流变时效模型，得到围岩的应变增量；S8：根据所述围岩的应变增量，构建应变强度随着应变增量变化的围岩动态损伤本构模型；所述步骤S3中，所述流-力耦合场围岩损伤程度D w 为：D w ＝D 1 +D 2其中，D 1 表示渗流作用下对围岩裂隙拓展的损伤程度且 D 2 表示流-力耦合作用下围岩岩土损伤程度且 表示无损状态下围岩初始孔隙率，v表示泊松比，E s 表示岩体的压缩模量，σ i 表示围岩应力且σ i (i＝1、2、3)，p w 表示渗流压力， 表示内摩擦角，F pk 表示岩体峰值强度，acr表示岩体各向异性系数且acr＝ln[Eε pk /(σ 1 -2vσ 3 )]，ε pk 表示峰值应变，E表示岩体的弹性模量；所述步骤S4中，所述折减后的部分参数包括折减后的内摩擦角、折减后的岩体强度和折减后的粘聚力；所述折减后的内摩擦角 为：所述折减后的岩体强度F w 为：所述折减后的粘聚力c w 为：c w ＝c(1-D w )其中，D w 表示流-力耦合场围岩损伤程度， 表示内摩擦角，F pk 表示岩体峰值强度，v表示泊松比，σ i 表示围岩应力且σ i (i＝1、2、3)，p w 表示渗流压力， 表示无损状态下围岩初始孔隙率，λ表示岩体导热系数，h表示表面传热系数，c表示粘聚力；所述步骤S5中，所述热-力耦合场损伤后的围岩损伤程度D T 为：其中，λ表示岩体导热系数，v表示泊松比，C p 表示岩体的比热容，F w 表示折减后的岩体强度且，σ i 表示围岩应力且σ i (i＝1、2、3)， 表示折减后的内摩擦角，c w 表示折减后的粘聚力，K v 表示岩体完整性系数，ctg是三角函数余切的符号且 所述步骤S6中，所述由于热对流引起的水-力-热多场耦合作用下围岩损伤程度D为：其中，D w 表示流-力耦合场围岩损伤程度，D T 表示热-力耦合场损伤后的围岩损伤程度，e αT 表示，h表示表面传热系数，p w 表示渗流压力， 表示无损状态下围岩初始孔隙率；所述步骤S7中，所述围岩的应变增量ε t 为：其中，ε 1 为沿最小主应力方向的塑性应变增量且 ε 2 为沿最大主应力方向的塑性应变增量且 I 1 表示应力张量的第一不变量；J 2 表示应力偏量的第二不变量，t表示从隧道开挖卸荷开始的时间，E表示岩体的弹性模量，D表示由于热对流引起的水-力-热多场耦合作用下围岩损伤程度，v表示泊松比；所述步骤S8中，所述围岩动态损伤本构模型包括动态损伤的岩体弹性模量、动态损伤的粘聚力和动态损伤的内摩擦角；所述动态损伤的岩体弹性模量E t 为：所述动态损伤的粘聚力c t 为：所述动态损伤的内摩擦角为：其中，σ i 表示围岩应力且σ i (i＝1、2、3)，v表示泊松比，D表示由于热对流引起的水-力-热多场耦合作用下围岩损伤程度，ε t 表示围岩的应变增量，E表示岩体的弹性模量。
2.根据权利要求1所述的围岩力学性质动态损伤本构模型的构建方法，其特征在于，所述步骤S1中，所述围岩性质的各项参数包括：岩体的压缩模量E s 、无损状态下围岩初始孔隙率 泊松比v、渗流压力p w 、围岩应力σ i (i＝1、2、3)、岩体峰值强度F pk 、峰值应变ε pk 、内摩擦角 岩体各向异性系数acr，岩体导热系数λ、表面传热系数h、岩体的比热容C p 和粘聚力c。
3.根据权利要求1所述的围岩力学性质动态损伤本构模型的构建方法，其特征在于，所述步骤S2中，所述渗流作用下对围岩裂隙拓展的损伤程度D 1 为：其中， 表示无损状态下围岩初始孔隙率，v表示泊松比，E s 表示岩体的压缩模量，σ i 表示围岩应力且σ i (i＝1、2、3)，p w 表示渗流压力。
4.根据权利要求1所述的围岩力学性质动态损伤本构模型的构建方法，其特征在于，所述步骤S2中，所述流-力耦合作用下围岩岩土损伤程度D 2 为：其中， 表示内摩擦角，σ i 表示围岩应力且σ i (i＝1、2、3)，F pk 表示岩体峰值强度，acr表示岩体各向异性系数且acr＝ln[Eε pk /(σ 1 -2vσ 3 )]，ε pk 表示峰值应变，v表示泊松比，E表示岩体的弹性模量，p w 表示渗流压力。
