一种加卸载环境中锚索支护锚固性能确定方法技术

本发明专利技术属于锚索支护数值计算领域，具体涉及一种加卸载环境中锚索支护锚固性能确定方法。包括在FLAC3D中创建三维模型并划分网格；定义模型的本构模型、材料参数、边界条件、初始条件；大变形计算模式为true；利用锚索结构单元安装锚索并设置锚索结构单元材料参数；定义FISH函数recording，记录锚索拉拔端的拉拔力、拉拔位移；定义FISH函数incorporating；第一次加载，进行时间步计算直至第一次加载模拟位移；卸载，进行时间步计算直至卸载模拟位移；第二次加载，并进行时间步计算直至第二次加载模拟位移。本发明专利技术提出的方法能够克服传统加卸载拉拔实验获取结果单一且往往无法获取锚固界面剪切应力分布形态等缺点，有效揭示加卸载环境中的锚索支护锚固机制。境中的锚索支护锚固机制。境中的锚索支护锚固机制。

【技术实现步骤摘要】
一种加卸载环境中锚索支护锚固性能确定方法


[0001]本专利技术属于锚索支护数值计算领域，具体涉及一种加卸载环境中锚索支护锚固性能确定方法。

技术介绍

[0002]锚索支护是煤矿巷道支护中的一种重要支护元件。与锚杆相比，锚索的长度通常在3m以上，而且在实际施工中，可以根据需求动态调整锚索长度，从而具有锚固长度大、锚固深度可动态调节等优点。随着煤炭深井高应力赋存环境日趋普遍，锚索支护在煤矿巷道围岩加固中有着越来越广泛的应用前景。
[0003]在实际工程应用中，锚索周围岩土体中的应力可能突然处于升高的状态，此时锚索周围岩土体运动将对锚索形成加载效应。但随着应力环境改变，锚索周围岩土体中的应力可能又会处于降低状态，此时锚索处于卸载状态。因此，在实际工程环境中，锚索支护极有可能处于加卸载环境。
[0004]为了检验锚索支护在加卸载环境中的锚固性能，对锚索支护进行加卸载拉拔实验是一种重要方法。但加卸载拉拔实验首先需要配备加卸载拉拔设备，包括千斤顶、压力盒、位移传感器、泵站等一系列实验器材，因此需要投入大量的物力和资金，实验成本比较高。此外，开展加卸载拉拔实验需要多名施工人员协同工作，因此对人力资源的要求较高。
[0005]其次，开展加卸载拉拔实验往往只能得到单一的拉拔力、拉拔位移关系，往往很难获得实验过程中锚固界面剪切应力结果。这是由于锚固界面埋置在锚索锚固体内，不便于提取锚索锚固体内的数据。因此，开展加卸载拉拔实验，不便于分析加卸载过程中锚索支护内部的力学传递规律。最后，开展加卸载拉拔实验很难获得锚索周围岩土体内的应力分布形态，因此不便于分析加卸载环境中锚索支护周围岩土体内的应力演化规律。
[0006]与传统加卸载拉拔实验相比，数值模拟已经成为研究锚索支护锚固性能的重要手段。由于数值模拟能够基于时间步迭代的方式分析锚索支护锚固性能，因此，用户可以分析每个时间步时锚索支护锚固体内的力学传递特征及应力分布形态。
[0007]现有数值模拟技术在分析锚索支护锚固性能时，往往只局限于单调加载的受力环境，并不能反映锚索支护在加卸载环境中的锚固性能。因此，提出一种加卸载环境中锚索支护锚固性能确定方法对揭示锚索支护在加卸载环境中的锚固机理具有重要意义。

技术实现思路

[0008]本专利技术的目的在于提出一种加卸载环境中锚索支护锚固性能确定方法。该方法克服了传统加卸载拉拔实验实施成本高、获取结果单一等缺点，能够有效模拟加卸载环境中锚索支护的锚固性能。
[0009]本专利技术采用如下技术方案，提供一种加卸载环境中锚索支护锚固性能确定方法，包括如下步骤：在FLAC3D中创建三维模型并划分网格；定义模型的本构模型、材料参数、边界条件、初始条件；设置模型的大变形计算模式为true；利用锚索结构单元安装锚索并设置
锚索结构单元材料参数；定义FISH函数recording，并利用FISH函数recording记录锚索拉拔端的拉拔力、拉拔位移；定义FISH函数incorporating；第一次加载，在锚索拉拔端设置恒定的拉拔速度，进行时间步计算直至第一次加载模拟位移；卸载，在锚索拉拔端设置恒定的推入速度，并进行时间步计算直至卸载模拟位移；第二次加载，在锚索拉拔端设置恒定的拉拔速度，并进行时间步计算直至第二次加载模拟位移。
[0010]作为上述技术方案的进一步描述：所述边界条件为：与锚索安装方向垂直且通过锚索拉拔端的边界面为滚筒支撑。
[0011]作为上述技术方案的进一步描述：所述锚索结构单元材料参数包括：grout
‑
friction为0；grout
‑
cohesion为0；grout
‑
stiffness为1
×
10
11
；slide为false；young与锚索杨氏模量一致；grout
‑
perimeter与锚索周长一致；cross
‑
sectional
‑
area与锚索横截面积一致。
[0012]作为上述技术方案的进一步描述：所述FISH函数recording逻辑结构如下：在每一次时间步计算时取出锚索拉拔端位置处锚索结构单元节点的地址，将其赋值给变量location；取出该节点沿拉拔方向的不平衡力，将该不平衡力取绝对值后赋值给变量force；取出该节点沿拉拔方向的位移，将该位移取绝对值后赋值给变量displacement；所述变量force为所述锚索拉拔端的拉拔力，所述变量displacement为所述锚索拉拔端的拉拔位移。
[0013]作为上述技术方案的进一步描述：所述FISH函数incorporating根据用户设定的参数自动修正锚固界面黏结滑移关系，其逻辑结构如下：在每一次时间步计算时取出锚索结构单元头指针并赋值给变量address；基于变量address，遍历锚索结构单元所有构件；遍历至任意一构件时，提取该构件第一号节点沿拉拔方向的位移并将其传递给变量s_m，将s_m代入第一公式并将计算结果传递给变量stress_Nm；将变量stress_Nm传递给该构件的struct.cable.grout.cohesion函数；所述第一公式为，公式中：F为单位长度锚固界面剪切力；为锚索杨氏模量；为锚索横截面积；k和t为调节系数；s为锚固界面剪切滑移量。
[0014]作为上述技术方案的进一步描述：所述拉拔速度与所述推入速度大小相等，小于等于1
×
10
‑6m/s，拉拔与推入方向相反。
[0015]作为上述技术方案的进一步描述：所述拉拔速度和所述推入速度通过锚索拉拔端位置处锚索结构单元节点的component
‑
id进行定位。
[0016]作为上述技术方案的进一步描述：所述第一次加载模拟位移由用户自行设定，其与加卸载拉拔实验第一次加载阶段内拉拔力最大值对应的拉拔位移间相对差值为第一次加载模拟位移相对差值；采用第二公
式计算，所述第二公式为，公式中：为第一次加载模拟位移相对差值，小于等于5%；为第一次加载模拟位移；为加卸载拉拔实验第一次加载阶段内拉拔力最大值对应的拉拔位移。
[0017]作为上述技术方案的进一步描述：所述卸载模拟位移由用户自行设定，其与加卸载拉拔实验卸载阶段内拉拔力最小值对应的拉拔位移间相对差值为卸载模拟位移相对差值；采用第三公式计算，所述第三公式为，公式中：为卸载模拟位移相对差值，小于等于5%；为卸载模拟位移；为加卸载拉拔实验卸载阶段内拉拔力最小值对应的拉拔位移。
[0018]作为上述技术方案的进一步描述：所述第二次加载模拟位移等于加卸载拉拔实验第二次加载阶段最后时刻的拉拔位移。
[0019]作为上述技术方案的进一步描述：第一次加载阶段内，调节系数k和t采用如下方式调参获得：将调节系数k和t代入第四公式直至同时满足以下两个条件：第一次加载阶段内，数值模拟中拉拔力、拉拔位移关系曲线与加卸载拉拔实验中拉拔力、拉拔位移关系曲线整体趋势一致；第一次加载本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种加卸载环境中锚索支护锚固性能确定方法，其特征在于，包括如下步骤：在FLAC3D中创建三维模型并划分网格；定义模型的本构模型、材料参数、边界条件、初始条件；设置模型的大变形计算模式为true；利用锚索结构单元安装锚索并设置锚索结构单元材料参数；定义FISH函数recording，并利用FISH函数recording记录锚索拉拔端的拉拔力、拉拔位移；所述FISH函数recording逻辑结构为：在每一次时间步计算时取出锚索拉拔端位置处锚索结构单元节点的地址，将其赋值给变量location；取出该节点沿拉拔方向的不平衡力，将该不平衡力取绝对值后赋值给变量force；取出该节点沿拉拔方向的位移，将该位移取绝对值后赋值给变量displacement；定义FISH函数incorporating；所述FISH函数incorporating根据用户设定的参数自动修正锚固界面黏结滑移关系，其逻辑结构为：在每一次时间步计算时取出锚索结构单元头指针并赋值给变量address；基于变量address，遍历锚索结构单元所有构件；遍历至任意一构件时，提取该构件第一号节点沿拉拔方向的位移并将其传递给变量s_m，将s_m代入第一公式并将计算结果传递给变量stress_Nm；将变量stress_Nm传递给该构件的struct.cable.grout.cohesion函数；所述第一公式为，公式中：F为单位长度锚固界面剪切力；为锚索杨氏模量；为锚索横截面积；k和t为调节系数；s为锚固界面剪切滑移量；第一次加载，在锚索拉拔端设置恒定的拉拔速度，进行时间步计算直至第一次加载模拟位移；卸载，在锚索拉拔端设置恒定的推入速度，并进行时间步计算直至卸载模拟位移；第二次加载，在锚索拉拔端设置恒定的拉拔速度，并进行时间步计算直至第二次加载模拟位移。2.根据权利要求1所述的锚固性能确定方法，其特征在于，所述边界条件为：与锚索安装方向垂直且通过锚索拉拔端的边界面为滚筒支撑。3.根据权利要求1所述的锚固性能确定方法，其特征在于，所述锚索结构单元材料参数包括：grout
‑
friction为0；grout
‑
cohesion为0；grout
‑
stiffness为1
×
10
11
；slide为false；young与锚索杨氏模量一致；grout
‑
perimeter与锚索周长一致；cross
‑
sectional
‑
area与锚索横截面积一致。4.根据权利要求1所述的锚固性能确定方法，其特征在于，所述拉拔速度...

【专利技术属性】
技术研发人员：陈见行，张文博，马俊明，赵仪强，杨胜利，张俊文，赵洪宝，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：