本发明属于岩土加固数值计算领域，具体涉及一种考虑锚固界面破坏后力学行为的数值计算方法。
背景技术：
：：1、锚杆(索)支护是岩土加固领域的重要支护技术。与传统砌碹支护、型钢支护等方式相比，锚杆(索)支护可深入岩土体内部，并与岩土体形成一个整体，共同协调变形。此外，锚杆(索)支护还可在外端头施加预紧力，并根据实际工况调节预紧力大小，从而对岩土体形成主动约束作用，提高支护效果。2、工程实践结果表明，在锚杆(索)支护过程中，锚杆(索)与锚固剂间锚固界面剪切破坏是一种常见破坏形式。为了分析锚固界面剪切破坏机理，研究人员对锚杆(索)支护系统开展了大量的室内物理实验及现场工程实践。这些室内物理实验及现场工程实践为揭示锚杆(索)与围岩间力学传递原理提供了理论依据，也阐明了锚杆(索)支护系统中锚固界面破坏机理。3、但室内物理实验与现场工程实践的缺点是往往很难研究各种不同工况条件下的锚杆(索)支护系统锚固性能。以室内物理实验为例，开展室内物理实验需要准备相应的实验设备、耗材、材料及实验平台等。若开展各种不同工况条件下的室内物理实验，需要投入大量的人力、物力、财力，可操作性及经济性较差。若依托现场工程实践开展各种不同工况条件下的锚杆(索)支护系统实验难度更大。这是由于现场工程环境是由地质条件决定，并不能严格符合研究人员预期的实验环境。此外，在现场开展不同工况条件下的锚杆(索)支护系统实验，可能诱发安全事故并威胁实验人员的安全。4、为了弥补室内物理实验和现场工程实践的不足，数值计算广泛应用于锚杆(索)支护系统锚固性能研究。与室内物理实验、现场工程实践相比，数值计算依赖计算机及数值计算软件。用户可在计算机中设定各种复杂的工况，并利用数值计算软件分析各种不同复杂工况条件下锚杆(索)支护系统锚固性能及破坏机理。5、在各种数值计算软件中，由itasca公司开发的flac3d软件得到广泛应用。在该软件中，itasca公司开发了锚单元、桩单元。这两种单元均可模拟岩土加固中的锚杆(索)支护。此外，itasca公司为这两种单元设定了本构方程及材料参数。因此，用户只需输入材料参数，便可分析锚杆(索)支护锚固性能，可操作性强且使用方便。6、但锚单元、桩单元存在明显缺陷。其中一项是忽略了锚固界面破坏后的力学行为。锚单元、桩单元采用理想弹塑性模型模拟锚固界面剪切破坏行为。当锚杆(索)支护系统受载后，锚固界面剪切应力随剪切位移增加而增加。当锚固界面剪切应力达到黏结强度后，锚固界面剪切应力保持不变。但实际工程中，锚固界面剪切应力达到黏结强度后会逐渐下降，表现出软化特性。因此，原始锚单元、原始桩单元无法准确模拟出锚固界面破坏后的力学行为。7、为了克服该问题，前人曾提出利用锚单元、桩单元中的内聚力参数调整锚固界面黏结强度，模拟出锚固界面峰前强化和峰后软化力学行为。这种方法需设定锚固界面摩擦角为零。此时，用户可利用fish函数或table函数将预定义黏结滑移关系嵌入锚单元、桩单元，模拟出锚固界面破坏后的力学行为。但这种方法忽略了岩土体中围压对锚杆(索)锚固系统的影响。由于锚固界面摩擦角为零，岩土体中围压不影响锚杆(索)锚固性能。8、前人室内物理实验结果表明，摩擦力直接影响锚固界面黏结强度。因此，设定锚固界面摩擦角为零，不能准确反映锚杆(索)支护系统锚固作用。9、此外，由于忽略了锚固界面摩擦作用，在不同围压条件下，锚固界面剪切应力达到黏结强度时对应的位移相同。但这与实际情况不符，因为，随锚杆(索)支护系统围压增大，锚固界面剪切应力达到黏结强度时对应的位移会明显增加。因此，这种方法无法准确模拟出锚固界面在有围压条件下的力学性能。10、针对上述问题，本发明提出一种考虑锚固界面破坏后力学行为的数值计算方法。该方法克服了原始锚单元无法模拟锚固界面破坏后力学行为的缺陷，也克服了仅调整内聚力参数无法考虑锚固界面摩擦作用的缺陷。本发明提出的方法对揭示锚杆(索)支护锚固机理及破坏后力学行为具有重要意义。技术实现思路1、本发明的目的在于提供一种考虑锚固界面破坏后力学行为的数值计算方法。该方法改善了现有数值计算方法无法准确模拟锚固界面破坏后力学行为的缺陷。2、本发明采用如下技术方案，提供一种考虑锚固界面破坏后力学行为的数值计算方法，包括：在flac3d中创建网格并划分单元体形成数值模型；定义数值模型的本构模型、材料参数；定义数值模型的位移边界条件、围压边界条件；在数值模型内部安装锚单元并对锚单元划分锚固构件；设置锚单元材料参数；定义fish函数results并利用其记录拉拔力、拉拔位移；定义fish函数postfailure并利用其实现锚固界面内聚力软化、锚固界面摩擦角软化；在锚单元拉拔端施加恒定的拉拔速度；设置整个数值模型的largestrain为true；定义时间步数并完成时间步计算。3、作为上述技术方案的进一步描述：4、数值模型的本构模型设置为弹性模型；数值模型的材料参数包括杨氏模量、泊松比。5、作为上述技术方案的进一步描述：6、数值模型的位移边界条件为：与锚单元安装方向垂直且穿过锚单元拉拔端的面为滚筒支撑。7、作为上述技术方案的进一步描述：8、数值模型的围压边界条件为：与锚单元安装方向平行的面为压应力边界条件。9、作为上述技术方案的进一步描述：10、对锚单元划分锚固构件时，锚固构件数量大于1。11、作为上述技术方案的进一步描述：12、锚单元材料参数包括：锚固界面内聚力初始值ci、锚固界面摩擦角初始值fi、锚固界面剪切刚度、滑动参数、锚固构件杨氏模量、锚固构件周长、锚固构件横截面积。13、作为上述技术方案的进一步描述：14、所述fish函数results逻辑结构如下：查找锚单元拉拔端位置处节点指针；将其赋值给变量pn；查找变量pn所在地址沿第1自由度方向的不平衡力；将其取绝对值后赋值给变量force，变量force即为拉拔力；查找变量pn所在地址沿第1自由度方向的位移；将其取绝对值后赋值给变量disp，变量disp即为拉拔位移。15、作为上述技术方案的进一步描述：16、所述fish函数postfailure逻辑结构如下：定义变量a、变量b、变量d，其中变量a设置为零，变量b大于零且小于1，变量d大于零；查找结构单元头指针并赋值给变量pb；设置循环，判断变量pb是否为空；若不为空，检查变量pb所在地址第1号节点滑移状态是否大于零；若滑移状态大于零，检查变量a是否等于零；若变量a等于零，取出变量pb所在地址第1号节点的位移，将其取绝对值后赋值给变量sp，变量sp即为锚固界面破坏时对应的位移；设置变量a为1；变量pb指向下一结构单元并继续循环；若变量a不等于零，取出变量pb所在地址第1号节点的位移，将其取绝对值后赋值给变量s；将变量s代入第1公式并将计算结果赋值给pb所在地址的锚固界面内聚力函数；将变量s代入第2公式并将计算结果赋值给pb所在地址的锚固界面摩擦角函数；变量pb指向下一结构单元并继续循环；若滑移状态不大于零，变量pb指向下一结构单元并继续循环；当变量pb为空时，结束循环。17、作为上述技术方案的进一步描述：18、所述第1公式为公式中：c为锚固界面内聚力函数；ci为锚固界面内聚力初始值；s为pb所在地址第1号节点位移；sp为锚固界面破坏时对应的位移。19、作为上述技术方案的进一步描述：20、所述第2公式为公式中：f为锚固界面摩擦角函数；fi为锚固界面摩擦角初始值。21、作为上述技术方案的进一步描述：22、拉拔速度小于等于1×10-6m/s。23、本发明有益效果主要包括以下六个方面：24、1.本发明提供了一种考虑锚固界面破坏后力学行为的数值计算方法。该方法能模拟锚杆(索)支护系统锚固性能及锚固界面破坏后的力学行为。为用户研究锚杆(索)支护系统锚固性能及锚固界面破坏机理提供了新思路、新方法。25、2.设计了自定义fish函数results。利用其能获得锚杆(索)拉拔过程中的拉拔力、拉拔位移。为用户分析锚杆(索)拉拔过程中的锚固性能提供了基础数据。26、3.设计了自定义fish函数postfailure。利用其能获得锚固界面剪切破坏时对应的位移，并实现锚固界面内聚力软化、锚固界面摩擦角软化；能有效模拟锚固界面破坏后的力学行为。克服了原始锚单元无法模拟锚固界面破坏后力学行为的缺陷。27、4.在本发明中，锚固界面摩擦角不为零。在应用本发明时，用户首先设定锚固界面摩擦角初始值。随后，fish函数postfailure会根据锚固界面是否破坏，自动修正和软化锚固界面摩擦角，且软化后的摩擦角仍大于零。本发明在有效模拟锚固界面破坏后力学行为的同时，又能考虑岩土体中围压对锚杆(索)锚固性能的影响。克服了前人将锚固界面摩擦角设置为零且无法考虑岩土体围压对锚杆(索)锚固性能影响的缺陷。28、5.本发明设计了第1公式、第2公式。基于第1公式、第2公式，用户可利用上述2个公式中的变量b、d调整锚固界面内聚力软化程度、锚固界面摩擦角软化程度。最终利用不同软化程度的锚固界面内聚力、锚固界面摩擦角模拟不同工况条件下锚杆(索)支护系统锚固性能。29、6.所述第1公式、所述第2公式中变量b数值大小可以不同。所述第1公式、所述第2公式中变量d数值大小也可以不同。用户应用本发明时，可根据实际情况对锚固界面内聚力、锚固界面摩擦角实现不同程度的软化。能分别研究锚固界面内聚力、锚固界面摩擦角在不同软化程度条件下对锚杆(索)锚固性能的影响。当前第1页12当前第1页12