一种模拟巷道围岩蠕变失稳的离散元方法及系统

技术特征：
1.一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，包括：获取待模拟巷道围岩的赋存状况、待模拟巷道设定范围内岩层的物理力学参数和待模拟巷道围岩的变形量随时间变化的曲线；所述赋存状况包括：巷道尺寸、埋深、巷道周围采掘情况和设定范围内的岩体参数；所述岩体参数包括：岩体倾角、岩体厚度、岩体岩性和岩体的层理结构；所述物理力学参数包括：单轴压力试验应力
–
应变曲线、巴西劈裂实验应力
–
应变曲线和岩石蠕变曲线；利用udec离散元模拟软件建立压缩劈裂模型并根据所述单轴压力试验应力
–
应变曲线和所述巴西劈裂实验应力
–
应变曲线对所述压缩劈裂模型的接触参数和块体参数进行反演得到合理接触面参数值和合理块体参数值；所述压缩劈裂模型包括单轴压缩数值模型和巴西劈裂数值模型；根据所述合理块体参数值和所述合理接触面参数值利用udec离散元模拟软件建立块体单元为cvisc蠕变模型的第一单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第一单轴蠕变数值模型的块体蠕变参数进行反演得到合理块体蠕变参数值；根据所述合理接触面参数值和所述合理块体蠕变参数值利用udec离散元模拟软件构建块体单元为改进cvisc蠕变模型的第二单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第二单轴蠕变数值模型的块体蠕变参数退化方程进行反演得到合理块体蠕变参数退化方程；所述第二单轴蠕变数值模型的块体单元为改进cvisc蠕变模型；所述改进cvisc蠕变模型为退化方程为块体蠕变参数退化方程的cvisc蠕变模型；根据所述待模拟巷道围岩的变形量随时间变化的曲线、所述赋存状况和所述合理块体蠕变参数退化方程构建巷道离散元数值仿真模型并采用离散法对所述巷道离散元数值仿真模型进行网格划分得到经离散网格化处理的巷道离散元数值仿真模型；根据所述待模拟巷道围岩的赋存状况、所述待模拟巷道设定范围内岩层的物理力学参数和所述待模拟巷道围岩的变形量随时间变化的曲线对所述经离散网格化处理的巷道离散元数值仿真模型进行模拟开挖得到各网格的累计张拉裂隙长度和累计剪切裂隙长度；所述累计张拉裂隙长度为所述网格内所有的张拉裂隙长度和；所述累计剪切裂隙长度为所述网格内所有的剪切裂隙长度和；根据各所述网格的累计张拉裂隙长度、累计剪切裂隙长度和各所述网格的临界蠕变损伤值确定所述待模拟巷道围岩蠕变失稳的范围；根据各所述网格的岩石蠕变损伤和临界蠕变损伤值确定所述待模拟巷道围岩蠕变失稳的范围。2.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述利用udec离散元模拟软件建立压缩劈裂模型并根据所述单轴压力试验应力
–
应变曲线和所述巴西劈裂实验应力
–
应变曲线对所述压缩劈裂模型的接触参数和块体参数进行反演得到合理接触面参数值和合理块体参数值，具体包括：利用udec离散元模拟软件构建单轴压缩数值模型和巴西劈裂数值模型；使用迭代试错法根据所述单轴压力试验应力
–
应变曲线对所述单轴压缩数值模型的块体参数和接触面参数进行反演以使所述单轴压缩数值模型与所述单轴压力试验应力
–
应变曲线相匹配，根据所述巴西劈裂实验应力
–
应变曲线对所述巴西劈裂数值模型的块体参数和接触面参数进行反演以使所述巴西劈裂数值模型与所述巴西劈裂实验应力
–
应变曲线相
匹配，得到合理接触面参数值和合理块体参数值。3.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述根据所述合理块体参数值和所述合理接触面参数值利用udec离散元模拟软件建立块体单元为cvisc蠕变模型的第一单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第一单轴蠕变数值模型的块体蠕变参数进行反演得到合理块体蠕变参数值，具体包括：利用udec离散元模拟软件建立块体单元为cvisc蠕变模型、块体参数为所述合理块体参数值，接触面参数为所述合理接触面参数值的第一单轴蠕变数值模型；使用迭代试错法根据所述岩石蠕变曲线对所述第一单轴蠕变数值模型的块体蠕变参数进行反演以使所述第一单轴蠕变数值模型与所述岩石蠕变曲线相匹配得到合理块体蠕变参数值。4.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述根据所述合理接触面参数值和所述合理块体蠕变参数值利用udec离散元模拟软件构建块体单元为改进cvisc蠕变模型的第二单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第二单轴蠕变数值模型的块体蠕变参数退化方程进行反演得到合理块体蠕变参数退化方程，具体包括：利用udec离散元模拟软件建立块体单元为改进cvisc蠕变模型、接触面参数为所述合理接触面参数值，块体蠕变参数为所述合理块体蠕变参数值的第二单轴蠕变数值模型；使用迭代试错法根据所述岩石蠕变曲线对所述第二单轴蠕变数值模型的改进cvisc蠕变模型中的块体蠕变参数退化方程中的蠕变常数进行反演以使所述第二单轴蠕变数值模型与所述岩石蠕变曲线相匹配得到合理块体蠕变参数退化方程。5.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述根据所述待模拟巷道围岩的变形量随时间变化的曲线、所述赋存状况和所述合理块体蠕变参数退化方程构建巷道离散元数值仿真模型，具体包括：利用udec离散元模拟软件根据所述赋存状况构建接触面参数为所述合理接触面参数值、块体参数为所述合理块体参数值、块体蠕变参数为所述合理块体蠕变参数值、退化方程为所述合理块体蠕变参数退化方程的初始仿真模型；对所述初始仿真模型进行模拟开挖并采用迭代试错法对所述初始仿真模型的接触面参数、块体参数和块体蠕变参数进行反演以使所述初始仿真模型与所述待模拟巷道围岩的变形量随时间变化的曲线相匹配得到巷道离散元数值仿真模型。6.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述临界蠕变损伤值的确定过程具体包括：根据所述巷道离散元数值仿真模型的接触面参数值、块体参数值、块体蠕变参数值和块体蠕变参数退化方程构建第三单轴蠕变数值模型；采用分级加载方式，根据所述第三单轴蠕变数值模型进行单轴蠕变数值模拟，并以单轴蠕变数值模拟过程中试件发生加速蠕变时的临界蠕变损伤值作为所述临界蠕变损伤值。7.根据权利要求1所述的一种模拟巷道围岩蠕变失稳的离散元方法，其特征在于，所述根据各所述网格的累计张拉裂隙长度、累计剪切裂隙长度和各所述网格的临界蠕变损伤值确定所述待模拟巷道围岩蠕变失稳的范围，具体包括：根据各所述网格的累计张拉裂隙长度和累计剪切裂隙长度确定各所述网格的岩石蠕
变损伤；根据各所述网格的岩石蠕变损伤和各所述网格的临界蠕变损伤值确定所述待模拟巷道围岩蠕变失稳的范围。8.一种模拟巷道围岩蠕变失稳的离散元系统，其特征在于，包括：获取模块，用于获取待模拟巷道围岩的赋存状况、待模拟巷道设定范围内岩层的物理力学参数和待模拟巷道围岩的变形量随时间变化的曲线；所述赋存状况包括：巷道尺寸、埋深、巷道周围采掘情况和设定范围内的岩体参数；所述岩体参数包括：岩体倾角、岩体厚度、岩体岩性和岩体的层理结构；所述物理力学参数包括：单轴压力试验应力
–
应变曲线、巴西劈裂实验应力
–
应变曲线和岩石蠕变曲线；合理接触面参数值和合理块体参数值确定模块，用于利用udec离散元模拟软件建立压缩劈裂模型并根据所述单轴压力试验应力
–
应变曲线和所述巴西劈裂实验应力
–
应变曲线对所述压缩劈裂模型的接触参数和块体参数进行反演得到合理接触面参数值和合理块体参数值；所述压缩劈裂模型包括单轴压缩数值模型和巴西劈裂数值模型；合理块体蠕变参数值确定模块，用于根据所述合理块体参数值和所述合理接触面参数值利用udec离散元模拟软件建立块体单元为cvisc蠕变模型的第一单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第一单轴蠕变数值模型的块体蠕变参数进行反演得到合理块体蠕变参数值；合理块体蠕变参数退化方程确定模块，用于根据所述合理接触面参数值和所述合理块体蠕变参数值利用udec离散元模拟软件构建块体单元为改进cvisc蠕变模型的第二单轴蠕变数值模型并根据所述岩石蠕变曲线对所述第二单轴蠕变数值模型的块体蠕变参数退化方程进行反演得到合理块体蠕变参数退化方程；所述第二单轴蠕变数值模型的块体单元为改进cvisc蠕变模型；所述改进cvisc蠕变模型为退化方程为块体蠕变参数退化方程的cvisc蠕变模型；巷道离散元数值仿真模型确定模块，用于根据所述待模拟巷道围岩的变形量随时间变化的曲线、所述赋存状况和所述合理块体蠕变参数退化方程构建巷道离散元数值仿真模型并采用离散法对所述巷道离散元数值仿真模型进行网格划分得到经离散网格化处理的巷道离散元数值仿真模型；裂隙长度确定模块，用于根据所述待模拟巷道围岩的赋存状况、所述待模拟巷道设定范围内岩层的物理力学参数和所述待模拟巷道围岩的变形量随时间变化的曲线对所述经离散网格化处理的巷道离散元数值仿真模型进行模拟开挖得到各网格的累计张拉裂隙长度和累计剪切裂隙长度；所述累计张拉裂隙长度为所述网格内所有张拉裂隙长度的和；所述累计剪切裂隙长度为所述网格内所有剪切裂隙长度的和；失稳范围确定模块，用于根据各所述网格的累计张拉裂隙长度、累计剪切裂隙长度和各所述网格的临界蠕变损伤值确定所述待模拟巷道围岩蠕变失稳的范围。9.根据权利要求8所述的一种模拟巷道围岩蠕变失稳的离散元系统，其特征在于，所述合理接触面参数值和合理块体参数值确定模块，具体包括：单轴压缩数值模型和巴西劈裂数值模型构建单元，用于利用udec离散元模拟软件构建单轴压缩数值模型和巴西劈裂数值模型；合理接触面参数值和合理块体参数值确定单元，用于使用迭代试错法根据所述单轴压
力试验应力
–
应变曲线对所述单轴压缩数值模型的块体参数和接触面参数进行反演以使所述单轴压缩数值模型与所述单轴压力试验应力
–
应变曲线相匹配，根据所述巴西劈裂实验应力
–
应变曲线对所述巴西劈裂数值模型的块体参数和接触面参数进行反演以使所述巴西劈裂数值模型与所述巴西劈裂实验应力
–
应变曲线相匹配，得到合理接触面参数值和合理块体参数值。10.根据权利要求8所述的一种模拟巷道围岩蠕变失稳的离散元系统，其特征在于，所述合理块体蠕变参数值确定模块，具体包括：第一单轴蠕变数值模型构建单元，用于利用udec离散元模拟软件建立块体单元为cvisc蠕变模型、块体参数为所述合理块体参数值，接触面参数为所述合理接触面参数值的第一单轴蠕变数值模型；合理块体蠕变参数值确定单元，用于使用迭代试错法根据所述岩石蠕变曲线对所述第一单轴蠕变数值模型的块体蠕变参数进行反演以使所述第一单轴蠕变数值模型与所述岩石蠕变曲线相匹配得到合理块体蠕变参数值。

技术总结
本发明涉及一种模拟巷道围岩蠕变失稳的离散元方法及系统。所述方法包括：根据应力


技术研发人员：李国栋 刘洪林 杭银建 马述起 刘金虎 王宏志 张卫东 李剣峰
受保护的技术使用者：中国矿业大学 徐州矿务集团有限公司 西南交通大学
技术研发日：2022.02.25
技术公布日：2022/5/31