1.一种考虑弯矩贡献因子的二维时效破裂模型的构建方法，其特征在于：包括如下步骤：
步骤1：设定岩体细观颗粒平行粘结接触的几何参数量包括平行粘结面积和平行粘结惯性矩，Ra、Rb分别为二维平行粘结接触a端的颗粒半径、b端的颗粒半径，为岩体细观颗粒二维平行粘结直径乘数，在二维情况下，平行粘结单位厚度为1时的岩体细观颗粒二维平行粘结面积A和平行粘结惯性矩I分别通过公式(2)、公式(3)来确定：
其中：为岩体细观颗粒二维平行粘结半径，为岩体细观颗粒二维平行粘结直径乘数，A为岩体细观颗粒二维平行粘结面积，I为岩体细观颗粒二维平行粘结惯性矩；
步骤201：利用岩体细观颗粒二维平行粘结时效衰减劣化的初始时间步长增量Δt，通过指数型函数计算岩体细观颗粒二维平行粘结时效衰减劣化的直径，公式(4)来确定：
其中： 为考虑弯矩贡献因子的二维平行粘结法向应力， 为判断岩体细观颗粒二维平行粘结开始时效劣化衰减时的应力阀值， 为岩体细观颗粒二维平行粘结的拉伸强度，为考虑弯矩贡献因子的二维平行粘结应力比，β1、β2分别为岩体细观颗粒平行粘结时效劣化的第一控制参数、第二控制参数， 为岩体细观颗粒二维平行粘结随时间劣化衰减的直径， 为岩体细观颗粒二维平行粘结未衰减时的平行粘结直径，e为自然常数，Δt为岩体细观颗粒二维平行粘结时效衰减劣化的初始时间步长增量；
步骤202：根据步骤201中的公式(4)，设定岩体细观颗粒二维平行粒粘结直径的指数型时效衰减因子β，见公式(5)：
其中： A'、I'、分别表示为岩体细观颗粒二维平行粘结随时间劣化衰减的平行粘结直径、平行粘结半径、平行粘结面积、平行粘结惯性矩、平行粘结直径乘数，A、I、为岩体细观颗粒二维平行粘结未衰减时的平行粘结直径、平行粘结
半径、平行粘结面积、平行粘结惯性矩、平行粘结直径乘数；
步骤203：根据步骤1中的公式(1)～公式(3)以及步骤202中的公式(5)，设定岩体细观颗粒二维平行粘结几何参数时效劣化衰减模式；岩体细观颗粒二维平行粘结直径随着时间增加而不断劣化衰减，二维平行粘结单位厚度为1时的平行粘结面积和平行粘结惯性矩也随着时间增加而不断劣化衰减，分别见公式(6)、公式(7)：
其中： A'、I'分别表示为岩体细观颗粒二维平行粘结随时间劣化衰减的平行粘结半径、平行粘结面积、平行粘结惯性矩，A、I为岩体细观颗粒平行粘结未衰减时的平行粘结面积、平行粘结惯性矩；
步骤204：依次计算第j个至第k个的岩体细观颗粒包含时间效应的二维平行粘结法向弯矩增量；在二维情况下，由平行粘结两端颗粒的速度、角速度和给定的循环计算时间步长增量Δtc，通过公式(8)、公式(9)、公式(10)确定第i个岩体细观颗粒二维平行粘结相对转角 二维平行粘结法向增量位移 以及二维平行粘结切向增量位移 再结合
步骤203中的公式(7)和步骤202中的公式(5)，确定第i个岩体细观颗粒包含时间效应的二维平行粘结弯矩增量，具体见公式(11)：
其中，ff、j、k是自然数，且2≤j≤ff≤k，j为每次循环计算中，包含时间效应的岩体细观颗粒二维平行粘结衰减后未破裂的初始标号值，ff为中间某一个标号值，k为每次循环计算中，包含时间效应的岩体细观颗粒二维平行粘结衰减后未破裂的最末标号值，i为第一个至最后一个平行粘结颗粒标号值， 分别为第i个岩体细观颗粒二维平
行粘结接触的a端和b端的绝对运动速度和角速度，nn和ns为岩体细观颗粒二维平行粘结接触面的法向单位向量和切向单位向量， 和 分别为岩体细观颗粒二维平行粘结法向位移增量和切向位移增量， 为岩体细观颗粒平行粘结法向刚度， 为岩体细观颗粒二维平行粘结弯矩增；
步骤205：根据步骤203中的公式(6)和公式(7)、步骤204中的公式(8)、公式(9)和公式(11)以及步骤202中的公式(5)，依次更新计算第j个至第k个岩体细观颗粒平行粘结未破裂并包含时间效应的二维平行粘结法向力、切向力和切向弯矩；通过公式(12)、公式(13)、公式(14)计算第i个岩体细观颗粒接触的二维平行粘结法向力、切向力和切向弯矩；在二维情况下，通过公式(15)来确定岩体细观颗粒平行粘结法向弯矩：
法向力：
切向力：
切向弯矩：
法向弯矩：
其中： 分别为第i个岩体细观颗粒包含时间效应的
平行粘结法向力、平行粘结切向力、包含时间效应的平行粘结法向弯矩、平行粘结切向弯矩、平行粘结法向位移增量和平行粘结切向位移增量， 为岩体细观颗粒二维平行粘结切向刚度，+＝为加法自反运算符，-＝为减法自反运算符；
步骤206：设置弯矩贡献因子 考虑弯矩对岩体细观颗粒平行粘结法向正应
力的贡献程度，根据二维平行粘结正应力计算公式 和二维平行粘结剪应
力计算公式 同时将这两个公式中A、I以及 用A'、I'及 替换，然后将步骤203中的公式(6)和公式(7)以及步骤202中的公式(5)代入，获得包含指数型时间效应和弯矩贡献因子的二维平行粘结正应力 和包含指数型时间效应的二维平行粘结剪应力 计算公式，分别见公式(16)和公式(17)；
步骤207：将步骤206中包含时间效应的 代入公式(18)，可确定考虑弯矩贡献因子且带拉伸截止限的摩尔-库伦时效破裂准则，并且依次计算第j个至第k个的二维平行粘结应力，用于判断岩体细观颗粒平行粘结是否破裂以及破裂模式；在该准则的岩体细观颗粒二维平行粘结应力中包含了指数型时间效应和考虑弯矩贡献因子；
其中：fs、fn分别为岩体细观颗粒二维平行粘结的时效剪切破裂准则、时效拉伸破裂准则， 为第i个接触的含指数型时间效应的二维平行粘结剪应力， 为第i个接触的含指数型时间效应且考虑弯矩贡献因子的二维平行粘结正应力， 分别为岩体细观颗粒二维平行粘结的拉伸强度、抗剪强度，为岩体细观颗粒二维平行粘结的粘聚力，为岩体细观颗粒二维平行粘结的内摩擦角；fs大于等于0表示平行粘结剪切破裂，小于0表示平行粘结未发生剪切破裂；fn大于等于0表示平行粘结拉伸破裂，小于0表示平行粘结未发生拉伸破裂；
步骤208：如果步骤207中的公式(18)中的fs或fn大于等于0，表明岩体细观颗粒的平行粘结发生了破裂，此后岩体细观颗粒的运动模式采用考虑阻尼效应的二维线性接触模型来表达；如果步骤207中的公式(18)中的fs和fn都小于0，表明平行粘结未破裂，继续循环步骤
201至207，计算、更新、判断岩体细观颗粒接触的平行粘结状态，直至岩体不产生新的平行粘结破裂或者平行粘结破裂加速发展而形成宏观破坏，循环终止。
2.根据权利要求1所述的考虑弯矩贡献因子的二维时效破裂模型的构建方法，其特征在于：所述岩体细观颗粒二维平行粘结时效衰减劣化的初始时间步长增量Δt的确定，是采用考虑弯矩贡献因子的平行粘结时效劣化衰减的二维指数型模式，由每个时间步内的二维平行粘结首次衰减破裂所损耗的时间来确定，也即通过第一个平行粘结按指数型模式进行衰减破裂所历时的时间除以直至第一个平行粘结破裂所需要的计算循环次数来估算初始时间 步 长 增量 Δ t ，见 公 式 其 中 ，
为第i个接触的岩体细观颗粒二维平行粘结直径乘数，nc为第一个
岩体细观颗粒二维平行粘结破裂所需的循环计算的次数，βσ、βτ分别为岩体细观颗粒二维平行粘结拉伸强度对应的时效劣化因子、二维平行粘结剪切强度对应的时效劣化因子，i依次为第一个至最后一个岩体细观颗粒平行粘结数，∞为无穷大。
3.根据权利要求2所述的考虑弯矩贡献因子的二维时效破裂模型的构建方法，其特征在于：所述岩体细观颗粒二维平行粘结拉伸强度对应的时效劣化因子βσ和岩体细观颗粒二维平行粘结剪切强度对应的时效劣化因子βτ的确定包括如下步骤，其中，这些步骤中包含的公式下标1代表第一个按指数型模式进行时效衰减劣化的二维平行粘结破裂标号；
步骤211：在二维情况下，由岩体细观颗粒平行粘结两端颗粒的速度、角速度和给定的循环计算时间步长增量Δtc，通过公式 确定平行粘结接触的相对转
角 通过公式 确定平行粘结法向增量位移 通过公
式 确 定 平行 粘 结 切 向 增 量 位 移 通 过 公式
确定平行粘结接触的弯矩增量；
步骤212 ：根据步骤211 中的公式 通过公式
确 定 平 行 粘 结 法 向 力 ；根 据 步 骤 2 1 1 中 的 公 式
通过公式 确定平行粘结切向力；根据步骤
2 1 1 中 的 公 式 和 公 式 通 过 公 式
确定平行粘结切向弯矩；通过公式 确定平行粘结法向弯矩，
其中，+＝为加法自反运算符，-＝为减法自反运算符；
步骤213：在二维情况下，通过公式 确定平行粘结正应力，通
过公式 确定平行粘结剪应力，将这两个公式中A、I以及 用A'、I'及 替换，然后将步骤203中的公式(6)和公式(7)以及步骤202中的公式(5)代入，获得包含指数型时间效应和弯矩贡献因子的二维平行粘结正应力计算公式 和包含
指数型时间效应的二维平行粘结剪应力计算公式
步骤214：将 代入公式 并令β＝βσ，将
代入公式 并令β＝βτ，据此，分别得到所述岩体细观颗粒二维平行
粘结拉伸强度对应的时效劣化因子 以及岩体细
观颗粒二维平行粘结剪切强度对应的时效劣化因子
4.根据权利要求1所述的考虑弯矩贡献因子的二维时效破裂模型的构建方法，其特征在于：所述岩体细观颗粒平行粘结发生破裂后，岩体细观颗粒的运动模式采用考虑阻尼效应的二维线性接触模型来表达，用于描述岩体细观颗粒平行粘结时效破裂后细观颗粒的应力、变形及运行规律，考虑阻尼效应的二维线性接触模型的构建包括如下步骤：
步骤301：通过Monte Carlo搜索算法，遍历寻找岩体细观颗粒每个二维线性接触端a、二维线性接触端b，颗粒与颗粒或者颗粒与墙的中心坐标，在二维情况下，通过公式(19)计算两者中心距离：
其中：d为二维线性接触两端颗粒与颗粒或者颗粒与墙之间的中心距离， 为二维线性接触端a的坐标， 为二维线性接触端b的坐标；
步骤302：二维平面状态下岩体细观颗粒间每个接触点的单位向量通过公式(20)计算，如果是颗粒与颗粒之间的接触，利用步骤301中得到二维线性接触两端的中心点坐标及距离计算，如果是颗粒与墙接触，直接利用墙体的法向量等效替换来计算，确定每个接触段的单位向量：
其中：ni为接触的单位矢量， 为接触端b的方向向量， 为接触端a的方向向量，nwall为约束墙的方向向量；
步骤303：岩体细观颗粒平行粘结破裂后，二维线性接触段的接触重叠量U，通过步骤
301计算颗粒间距d以及二维线性接触两端的颗粒半径Ra、Rb，再利用公式(21)来确定；通过设定颗粒二维线性接触参考距离gr，并结合公式(22)，确定颗粒二维线性接触的距离gs：
gs＝|U|-gr  (22)
步骤304：确定岩体细观颗粒接触点法向、切向等效刚度，利用接触两端颗粒实体或是墙体的刚度ka，kb等效代替为接触点的刚度，按公式(23)计算：
其中：Kn、Ks为等效的法向刚度和切向刚度， 为颗粒与颗粒或者颗粒与墙的接触a端的法向和切向刚度， 为颗粒与颗粒或者颗粒与墙的接触b端的法向和切向刚度；
步骤305：确定岩体中接触两端颗粒间的相对运动速度，利用公式(24)、公式(25)来计算；其中epqz为Ricci指标置换符号，按照公式(26)来计算：
其中：Vp与Vq等价，Vp与Vq为岩体中接触两端颗粒间的相对运动速度，p、q为指标等价符号，p＝1，q＝1表示颗粒与颗粒接触，p＝2，q＝2时表示颗粒与墙接触， 为颗粒与颗粒或是颗粒与墙的接触b端单元的速度， 为颗粒与颗粒或是颗粒与墙的接触a端单元的速度， 为颗粒与颗粒或是颗粒与墙的接触a端单元的角速度，为颗粒与颗粒或是颗粒与墙的接触b端单元的角速度， 为颗粒与颗粒或是颗粒与墙的接触a端的位移， 为颗粒与颗粒或是颗粒与墙的接触b端的位移， 为位移指标变换的中间过渡符号， 表示指标符号为p时颗粒-颗粒或者颗粒-墙的接触a端的速度，表示指标符号为q时颗粒-颗粒或者颗粒-墙的接触a端的速度， 表示指标符号为p时颗粒-颗粒或者颗粒-墙的接触b端的速度， 表示指标符号为q时颗粒-颗粒或者颗粒-墙的接触b端的速度，只有a端和b端两个接触端；
步骤306：对于线性接触模型的初始时间步长增量Δt的取值，通过公式(29)估计最小的时间步长Δt，确保所构建模型的计算时间步长小于该值，即可保证系统积分计算趋于稳定；通过公式(30)、公式(31)、公式(32)确定每个线性接触的总位移增量、法向位移增量和切向位移增量：
R＝min(Ra,Rb)  (27)
ΔUp1＝Vp1Δt  (30)
a b
其中：R为岩体细观颗粒的等效半径，R、R 分别为二维平行粘结接触a端的颗粒半径、b端的颗粒半径，m为岩体细观颗粒质量，J1为岩体细观颗粒的转动惯量；k平为岩体细观颗粒系统平移刚度，k转为岩体细观颗粒系统旋转刚度；ΔUp1为岩体细观颗粒二维线性接触的总位移增量，Δδn、 物理意义相同，均表示岩体细观颗粒二维线性接触的法向位移增量，Δδs、 物理意义相同，均表示岩体细观颗粒二维线性接触的切向位移增量，Vp1与Vq1为岩体细观颗粒接触两端的相对运动速度，n为单位法向量，p1,q1为张量指标变换符号；
步骤307：由公式(22)判定岩体细观颗粒表面接触允许存在的最大距离，计算法向和切向位移更新因子，另外，岩体细观颗粒二维线性接触法向位移增量的更新是采用前一步的法向位移增量与更新因子α的乘积获得，岩体细观颗粒二维线性接触切向位移增量的更新是采用前一步的切向位移增量与更新因子α的乘积获得；
其中：(gs)0为模型计算初始时刻的表面接触距离，gs为岩体细观颗粒接触的距离，α为位移更新因子；
步骤308：岩体细观颗粒二维线性接触的法向线性力采取相对矢量累加(Ml＝1)和绝对矢量累加(Ml＝0)模式，通过公式(34)、(35)计算获得；岩体细观颗粒二维线性接触的切向线性力采用岩体细观颗粒接触滑动来表示，通过公式(36)计算获得：
其中： 分别为岩体细观颗粒间受力变形的二维线性接触法向线性力、切向线性力，kn、ks分别为岩体细观颗粒间受力变形的二维线性接触法向、切向线性刚度，Δδn、Δδs分别为岩体细观颗粒二维线性接触的法向位移增量、切向位移增量， 分别为岩体细观颗粒二维线性接触的初始法向力增量值、切向力增量值， 为岩体细观颗粒未滑动时的静摩擦力， 为岩体细观颗粒滑动摩擦力，其值可通过摩擦系数μ与 乘积得到；
步骤309：岩体细观颗粒线性接触的法向阻尼采用全法向模式Md＝{0,2}和无拉伸模式Md＝{1,3}两种，通过公式(37)计算，其中mc为等效颗粒质量，按公式(38)计算，岩体细观颗粒线性接触的切向阻尼采用全剪切模式Md＝{0,1}和滑-剪模式Md＝{2,3}，按照公式(39)来计算，
其中： 分别为岩体细观颗粒线性接触的法向线性阻尼力、切向线性阻尼力，βn为岩体细观颗粒线性接触的法向阻尼系数，βs为岩体细观颗粒线性接触的切向阻尼系数，kn为岩体细观颗粒线性接触的法向线性刚度，ks为岩体细观颗粒线性接触的切向线性刚度，分别为岩体细观颗粒线性接触的法向速率和切向速率，F*为岩体细观颗粒线性接触的全法向阻尼力，表达式为 mc为岩体细观等效颗粒质量，m(1)为岩体细观颗粒接触端1的细观颗粒质量，m(2)为岩体细观颗粒接触端2的细观颗粒质量，Fd为岩体细观颗粒线性接触的总阻尼力。