CN115358099A - 一种全长黏结锚索的数值模拟方法及系统 - Google Patents
一种全长黏结锚索的数值模拟方法及系统 Download PDFInfo
- Publication number
- CN115358099A CN115358099A CN202211286678.1A CN202211286678A CN115358099A CN 115358099 A CN115358099 A CN 115358099A CN 202211286678 A CN202211286678 A CN 202211286678A CN 115358099 A CN115358099 A CN 115358099A
- Authority
- CN
- China
- Prior art keywords
- anchoring
- anchor
- anchor cable
- shear
- numerical simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004873 anchoring Methods 0.000 claims abstract description 156
- 238000010008 shearing Methods 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 12
- 238000010276 construction Methods 0.000 claims description 49
- 239000000463 material Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 4
- 239000003795 chemical substances by application Substances 0.000 description 17
- 238000002474 experimental method Methods 0.000 description 13
- 239000011435 rock Substances 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Piles And Underground Anchors (AREA)
Abstract
本发明公开了一种全长黏结锚索的数值模拟方法及系统,该方法包括:创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;定义锚索与锚固剂间锚固界面的剪切滑移关系;创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。本发明能够解决原始的锚结构单元和桩结构单元无法模拟锚索与锚固剂间锚固界面失效后的非线性剪切软化和摩擦脱黏问题以及表格函数使用较为繁琐的问题,以准确模拟全长黏结锚索的锚固性能,并提高数值模拟的计算效率。
Description
技术领域
本发明涉及岩石力学技术领域,具体来说,涉及一种全长黏结锚索的数值模拟方法及系统。
背景技术
全长黏结锚索是煤矿巷道支护和硐室围岩加固中常用的支护工具,在维护巷道和硐室围岩稳定及安全使用方面发挥着至关重要的作用。因此,研究全长黏结锚索的锚固性能是岩石力学领域的一项热点内容。为了研究全长黏结锚索的锚固性能,研究人员往往采用室内或现场的锚索拉拔实验,检验全长黏结锚索的锚固力及峰值荷载对应的拉拔位移。但锚索拉拔实验存在很多缺陷。首先,为了进行锚索拉拔实验,研究人员需要购置相应的锚索拉拔设备及仪器,因此需要耗费大量的金钱。其次,为了开展锚索拉拔实验,研究人员需要携带锚索拉拔设备到达实验地点并在实验地点调试和组装锚索拉拔设备,而这又需要消耗大量的人力资源和时间资源。
随着计算机技术的快速进步,利用计算机开展数值模拟研究成为岩石力学领域行之有效的方法。为了对全长黏结锚索进行数值模拟,美国Itasca公司在三维有限差分数值模拟软件FLAC3D内开发了锚结构单元和桩结构单元,利用其可以对全长黏结锚索的轴向承载行为进行模拟。但这两种结构单元均存在明显缺点。主要缺点在于这两种结构单元将锚索与锚固剂间的锚固界面处理为理想弹塑性模型,忽略了锚索与锚固剂间锚固界面破坏后的力学行为,因此并不能准确模拟全长黏结锚索的轴向承载行为。
这种缺点具体表现在以下两方面。首先,物理实验结果表明,锚索与锚固剂间的锚固界面在达到其剪切强度后会表现出非线性剪切软化行为,而原始的锚结构单元和桩结构单元均不能模拟该行为。其次,当锚索与锚固剂间发生非线性剪切软化行为后,随着拉拔位移继续增加,锚固界面的剪切应力会逐渐趋于稳定并保持不变,表现出摩擦脱黏行为。但原始的锚结构单元和桩结构单元均无法模拟该行为。先前的研究人员曾经使用桩结构单元内的表格函数,实现对原始的桩结构单元进行改进。但该方法要求用户从锚固界面的剪切滑移曲线中手动提取数据点,并将所提取的数据点输入到表格函数中。该方法需要耗费大量的时间,使用较为繁琐,给实际使用造成不便。
因此,提出一种能够考虑锚索与锚固剂间锚固界面剪切软化和摩擦脱黏的且能够方便使用的数值模拟方法,对准确模拟全长黏结锚索的锚固性能具有重要意义。
发明内容
本发明提出一种全长黏结锚索的数值模拟方法及系统,能够解决原始的锚结构单元和桩结构单元无法模拟锚索与锚固剂间锚固界面失效后的非线性剪切软化和摩擦脱黏问题以及表格函数使用较为繁琐的问题,以准确模拟全长黏结锚索的锚固性能,并提高数值模拟的计算效率。
本发明的技术方案是这样实现的:
根据本发明的一方面,提供了一种全长黏结锚索的数值模拟方法。
该全长黏结锚索的数值模拟方法,包括:
创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
定义锚索与锚固剂间锚固界面的剪切滑移关系;
创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
其中,创建三维数值模拟模型包括:根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置。
可选的,所述模型配置的计算模式为大变形计算模式。
其中,在所述三维数值模拟模型中,创建锚索包括:根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。
可选的,每个锚固构建的长度需小于等于200mm。
可选的,所述基本力学性质参数包括:锚固界面剪切刚度、锚固界面法向刚度、锚固界面沿剪切方向的摩擦角、锚固开关、锚索杨氏模量、锚索泊松比、锚索横截面积、锚索周长、锚索抗拉极限、锚索极惯性矩、锚索沿Y轴的二次惯性矩和/或锚索沿Z轴的二次惯性矩。
其中,定义锚索与锚固剂间锚固界面的剪切滑移关系包括:定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。
其中,创建锚固界面剪切行为FISH函数包括:查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;为锚固界面剪切强度;为锚索直径;所述第二公式为:;其中,c为非线性剪切软化系数;所述第三公式为;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
根据本发明的另一方面,提供了一种全长黏结锚索的数值模拟系统。
该全长黏结锚索的数值模拟系统,包括:
模型创建模块,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
其中,所述模型创建模块在创建三维数值模拟模型时,可根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;并定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置。
可选的,所述模型配置的计算模式为大变形计算模式。
其中,所述模型创建模块在所述三维数值模拟模型中,创建锚索时,可根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。
可选的,每个锚固构建的长度需小于等于200mm。
可选的,所述基本力学性质参数包括:锚固界面剪切刚度、锚固界面法向刚度、锚固界面沿剪切方向的摩擦角、锚固开关、锚索杨氏模量、锚索泊松比、锚索横截面积、锚索周长、锚索抗拉极限、锚索极惯性矩、锚索沿Y轴的二次惯性矩和/或锚索沿Z轴的二次惯性矩。
其中,所述关系定义模块在定义锚索与锚固剂间锚固界面的剪切滑移关系是包括:定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。
其中,所述函数处理模块在创建锚固界面剪切行为FISH函数时,可查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;为锚固界面剪切强度;为锚索直径;所述第二公式为:;其中,c为非线性剪切软化系数;所述第三公式为;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
有益效果:
本发明能够模拟出锚索与锚固剂间锚固界面的弹性、非线性剪切软化和摩擦脱黏行为,与实际情况更为相符;改善了原始的桩结构单元无法模拟全长黏结锚索破坏后力学行为的缺点,从而能够更为准确的模拟全长黏结锚索的轴向承载行为;而FISH函数能够根据锚固界面的剪切滑移量自动修正锚固界面单位长度的剪切强度输入参数,从而避免了用户使用表格函数时需手动采集锚固界面剪切滑移关系数据点的工序,简化了计算工序,提高了计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本发明实施例的一种全长黏结锚索的数值模拟方法的流程示意图;
图2是根据本发明实施例的一种全长黏结锚索的数值模拟系统的结构框图;
图3是根据本发明实施例的锚索与锚固剂间锚固界面的剪切滑移关系示意图;
图4是根据本发明实施例的物理实验结果、原始模拟结果与本发明模拟结果之间的对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。
根据本发明的实施例,提供了一种全长黏结锚索的数值模拟方法及系统。
如图1所示,根据本发明实施例的一种全长黏结锚索的数值模拟方法,包括:
步骤S101,创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
步骤S103,定义锚索与锚固剂间锚固界面的剪切滑移关系;
步骤S105,创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
步骤S107,在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
在具体应用时,上述拉拔速度需小于等于1×10-6m/s。而创建三维数值模拟模型是,则可根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置。所述模型配置的计算模式为大变形计算模式。
在所述三维数值模拟模型中,创建锚时,则根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。每个锚固构建的长度需小于等于200mm。所述基本力学性质参数包括:锚固界面剪切刚度、锚固界面法向刚度、锚固界面沿剪切方向的摩擦角、锚固开关、锚索杨氏模量、锚索泊松比、锚索横截面积、锚索周长、锚索抗拉极限、锚索极惯性矩、锚索沿Y轴的二次惯性矩和/或锚索沿Z轴的二次惯性矩。
实际应用时,上述锚固界面剪切刚度可根据公式(1)进行计算;锚索极惯性矩则可根据公式(2)进行计算;锚索沿Y轴的二次惯性矩则可根据公式(3)进行计算;锚索沿Z轴的二次惯性矩则可根据公式(3)进行计算;而锚固界面法向刚度和锚固界面沿剪切方向的摩擦角统一设置为零。
其中,S k为锚固界面剪切刚度;τ max为锚固界面剪切强度;d 1为锚固界面剪切强度对应的剪切滑移量;D b为锚索直径;I为锚索极惯性矩;I y为锚索沿Y轴的二次惯性矩;I z为锚索沿Z轴的二次惯性矩。
在具体应用时,在定义锚索与锚固剂间锚固界面的剪切滑移关系时,可定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。而创建锚固界面剪切行为FISH函数时,可按照如下逻辑关系进行计算:查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;为锚固界面剪切强度;为锚索直径;所述第二公式为:;其中,c为非线性剪切软化系数;所述第三公式为;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
对应的,如图2所示,根据本发明实施例的一种全长黏结锚索的数值模拟系统,包括:
模型创建模块201,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块203,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块205,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块207,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
在具体应用时,所述模型创建模块201在创建三维数值模拟模型时,可根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;并定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置,所述模型配置的计算模式为大变形计算模式;所述模型创建模块201在所述三维数值模拟模型中,创建锚索时,可根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。
而所述关系定义模块203在定义锚索与锚固剂间锚固界面的剪切滑移关系是包括:定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。所述函数处理模块205在创建锚固界面剪切行为FISH函数时,可查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;为锚固界面剪切强度;为锚索直径;所述第二公式为:;其中,c为非线性剪切软化系数;所述第三公式为;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
为了方便了解本发明的上述技术方案,以下通过具体对比例,对本发明的上述技术方案进行进一步的描述。
为了验证本发明的有效性,采用本发明提出的数值模拟方法对Aoki et al(2003)开展的全长黏结锚索拉拔实验进行模拟,并将实验结果与数值模拟结果进行对比。
Aoki et al (2003)对直径为15.2mm的锚索进行了拉拔实验,锚索的杨氏模量为195GPa,整个锚索采用全长黏结的方式进行锚固且锚固长度为4m。为了对该拉拔实验进行模拟,采用如下的步骤:
S1:根据实际工况在FLAC3D内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高分别为4m、1m和1m;同时确定沿着三维数值模拟模型长、宽和高方向单元体的数量分别为10、50和10;定义三维数值模拟模型中单元体的本构模型为各向同性弹性模型,材料参数分别为杨氏模量20GPa,泊松比为0.25;
S2:定义该三维数值模拟模型的边界条件为靠近拉拔端的单元体面为滚筒支撑;设置三维数值模拟模型的初始条件为密度为2700kg/m3;定义该三维数值模拟模型的模型配置为大变形计算模式;
S3:在三维数值模拟模型内,利用桩结构单元创建锚索;在该模拟中,锚索的长度为4m,与实际工况中锚索的长度相同;同时,定义桩结构单元内锚固构建的数量为20,即每个锚固构建的长度为200mm,以符合本发明对每个锚固构建长度的要求;定义全长黏结锚索的基本力学性质参数,具体包括:锚固界面剪切刚度为1.53GPa、锚固界面法向刚度为0、锚固界面沿剪切方向的摩擦角为0、锚固开关为“true”、锚索杨氏模量为195GPa、锚索泊松比为0.25、锚索横截面积为181.46mm2、锚索周长47.75mm、锚索抗拉极限为260kN、锚索极惯性矩为5240.52mm4、锚索沿Y轴的二次惯性矩为2620.26mm4和锚索沿Z轴的二次惯性矩为2620.26mm4。
S4:定义锚索与锚固剂间锚固界面的剪切滑移关系,包括锚固界面的剪切强度为1.42MPa、锚固界面剪切强度对应的剪切滑移量为0.5mm、非线性剪切软化系数为200m-1、锚固界面残余剪切强度对应的剪切滑移量为4mm,该锚固界面剪切滑移关系如图3所示;
S5:创建本发明开发的锚固界面剪切行为FISH函数“interface_shear”,并将定义的锚固界面剪切滑移关系输入到该FISH函数“interface_shear”内;
S6:在全长黏结锚索的拉拔端施加恒定的拉拔速度且拉拔速度为1×10-6m/s,以实现对全长黏结锚索拉拔的模拟;
S7:记录全长黏结锚索拉拔端的拉拔力和拉拔位移,并以拉拔力为纵坐标,拉拔位移为横坐标,绘制全长黏结锚索拉拔力与拉拔位移间的关系曲线,如图4所示。通过该图可以看出,在物理实验中,全长黏结锚索的峰值荷载约196kN,利用本发明提出的数值模拟方法获得的锚索峰值荷载约194kN,与物理实验结果十分接近。此外,本发明提出的数值模拟方法获得的锚索承载性能曲线斜率及整体趋势与物理实验结果也非常接近,充分证明了本发明提出的数值模拟方法的有效性和准确性。与之相反,桩结构单元原始的模拟方法获得的结果与物理实验结果明显不符,充分说明桩结构单元原始的模拟方法并不能准确反应全长黏结锚索的锚固性能,而本发明提出的模拟方法可以有效解决桩结构单元的不足。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种全长黏结锚索的数值模拟方法,其特征在于,包括:
创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
定义锚索与锚固剂间锚固界面的剪切滑移关系;
创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
2.根据权利要求1所述的一种全长黏结锚索的数值模拟方法,其特征在于,创建三维数值模拟模型包括:
根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;
定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置。
3.根据权利要求2所述的一种全长黏结锚索的数值模拟方法,其特征在于,所述模型配置的计算模式为大变形计算模式。
4.根据权利要求3所述的一种全长黏结锚索的数值模拟方法,其特征在于,在所述三维数值模拟模型中,创建锚索包括:
根据实际工况,确定桩结构单元内锚固构建的数量;
根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;
其中,锚索长度与实际工况中的锚索长度相同。
5.根据权利要求4所述的一种全长黏结锚索的数值模拟方法,其特征在于,每个锚固构建的长度需小于等于200mm。
6.根据权利要求5所述的一种全长黏结锚索的数值模拟方法,其特征在于,所述基本力学性质参数包括:
锚固界面剪切刚度、锚固界面法向刚度、锚固界面沿剪切方向的摩擦角、锚固开关、锚索杨氏模量、锚索泊松比、锚索横截面积、锚索周长、锚索抗拉极限、锚索极惯性矩、锚索沿Y轴的二次惯性矩和/或锚索沿Z轴的二次惯性矩。
7.根据权利要求6所述的一种全长黏结锚索的数值模拟方法,其特征在于,定义锚索与锚固剂间锚固界面的剪切滑移关系包括:
定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。
8.根据权利要求7所述的一种全长黏结锚索的数值模拟方法,其特征在于,创建锚固界面剪切行为FISH函数包括:
查找桩结构单元的头指针并将其赋值给指针变量pnt;
利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;
判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;
在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;
在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
10.一种全长黏结锚索的数值模拟系统,其特征在于,包括:
模型创建模块,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211286678.1A CN115358099B (zh) | 2022-10-20 | 2022-10-20 | 一种全长黏结锚索的数值模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211286678.1A CN115358099B (zh) | 2022-10-20 | 2022-10-20 | 一种全长黏结锚索的数值模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115358099A true CN115358099A (zh) | 2022-11-18 |
CN115358099B CN115358099B (zh) | 2023-01-13 |
Family
ID=84007595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211286678.1A Active CN115358099B (zh) | 2022-10-20 | 2022-10-20 | 一种全长黏结锚索的数值模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115358099B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115795924A (zh) * | 2023-02-07 | 2023-03-14 | 中国矿业大学(北京) | 一种针对压花锚索的数值模拟方法 |
CN116502466A (zh) * | 2023-06-19 | 2023-07-28 | 中国矿业大学(北京) | 一种锚杆结构单元中锚固构件数量的确定方法 |
CN117113512A (zh) * | 2023-09-20 | 2023-11-24 | 中国矿业大学(北京) | 一种考虑锚固界面破坏后力学行为的数值计算方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130195368A1 (en) * | 2012-01-31 | 2013-08-01 | Synopsys, Inc. | Scalable pattern matching between a pattern clip and a pattern library |
CN111931383A (zh) * | 2020-08-24 | 2020-11-13 | 长江水利委员会长江科学院 | 一种锚杆/锚索剪切屈服破坏断裂模拟计算方法 |
CN113153315A (zh) * | 2021-04-26 | 2021-07-23 | 中国矿业大学 | 一种深部大断面切眼成巷方式及支护时机的确定方法 |
CN113153435A (zh) * | 2021-05-17 | 2021-07-23 | 中国矿业大学 | 一种双巷布置系统复采扰动下煤柱加固参数的确定方法 |
CN114547895A (zh) * | 2022-02-24 | 2022-05-27 | 中国矿业大学(北京) | 一种全长注浆锚索力学传递性能的确定方法 |
CN114810163A (zh) * | 2021-12-22 | 2022-07-29 | 晋能控股煤业集团有限公司 | 深部高应力薄层坚硬顶板底抽巷支护工艺 |
CN115081302A (zh) * | 2022-07-15 | 2022-09-20 | 中国矿业大学 | 支护构件与硐室围岩接触及相互作用的模拟方法和系统 |
-
2022
- 2022-10-20 CN CN202211286678.1A patent/CN115358099B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130195368A1 (en) * | 2012-01-31 | 2013-08-01 | Synopsys, Inc. | Scalable pattern matching between a pattern clip and a pattern library |
CN111931383A (zh) * | 2020-08-24 | 2020-11-13 | 长江水利委员会长江科学院 | 一种锚杆/锚索剪切屈服破坏断裂模拟计算方法 |
CN113153315A (zh) * | 2021-04-26 | 2021-07-23 | 中国矿业大学 | 一种深部大断面切眼成巷方式及支护时机的确定方法 |
CN113153435A (zh) * | 2021-05-17 | 2021-07-23 | 中国矿业大学 | 一种双巷布置系统复采扰动下煤柱加固参数的确定方法 |
CN114810163A (zh) * | 2021-12-22 | 2022-07-29 | 晋能控股煤业集团有限公司 | 深部高应力薄层坚硬顶板底抽巷支护工艺 |
CN114547895A (zh) * | 2022-02-24 | 2022-05-27 | 中国矿业大学(北京) | 一种全长注浆锚索力学传递性能的确定方法 |
CN115081302A (zh) * | 2022-07-15 | 2022-09-20 | 中国矿业大学 | 支护构件与硐室围岩接触及相互作用的模拟方法和系统 |
Non-Patent Citations (3)
Title |
---|
JIANHANG CHEN 等: "Studying the performance of fully encapsulated rock bolts with modified structural elements", 《INT J COAL SCI TECHNOL》 * |
言志信 等: "锚固参数对地震作用下岩质边坡锚固界面剪切作用影响的数值分析", 《煤炭学报》 * |
郭金刚 等: "基于残余剪切强度的全长黏结锚杆拉拔模拟", 《岩土力学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115795924A (zh) * | 2023-02-07 | 2023-03-14 | 中国矿业大学(北京) | 一种针对压花锚索的数值模拟方法 |
CN115795924B (zh) * | 2023-02-07 | 2023-04-18 | 中国矿业大学(北京) | 一种针对压花锚索的数值模拟方法 |
CN116502466A (zh) * | 2023-06-19 | 2023-07-28 | 中国矿业大学(北京) | 一种锚杆结构单元中锚固构件数量的确定方法 |
CN116502466B (zh) * | 2023-06-19 | 2023-09-05 | 中国矿业大学(北京) | 一种锚杆结构单元中锚固构件数量的确定方法 |
CN117113512A (zh) * | 2023-09-20 | 2023-11-24 | 中国矿业大学(北京) | 一种考虑锚固界面破坏后力学行为的数值计算方法 |
CN117113512B (zh) * | 2023-09-20 | 2024-02-09 | 中国矿业大学(北京) | 一种考虑锚固界面破坏后力学行为的数值计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115358099B (zh) | 2023-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115358099B (zh) | 一种全长黏结锚索的数值模拟方法及系统 | |
Nandwana et al. | Detection of the location and size of a crack in stepped cantilever beams based on measurements of natural frequencies | |
Otani | Inelastic analysis of R/C frame structures | |
CN111666699A (zh) | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 | |
CN113111559A (zh) | 一种钢混组合结构粘结滑移的有限元模拟方法 | |
CN107688691A (zh) | 一种木结构用齿板连接节点受力性能的数值模拟方法 | |
Chaudhari et al. | Effect of piled raft design on high-rise building considering soil structure interaction | |
Carpinteri et al. | The effect of contact on the decohesion of laminated beams with multiple microcracks | |
Ile et al. | Shaking table tests of lightly rc walls: Numerical simulations | |
Hu et al. | Multi-scale modelling of flexible end plate connections under fire conditions | |
Hu et al. | Modelling of flexible end plate connections in fire using cohesive elements | |
Kuijvenhoven | A design method for timber grid shells | |
Ronagh et al. | A rational model for the distortional buckling of tapered members | |
Costa et al. | Numerical modelling approaches for existing masonry and RC structures | |
Mougaard et al. | A partly and fully cracked triangular XFEM element for modeling cohesive fracture | |
TORKI et al. | Effect of opening dimensions on the relative flexural operation of coupled shear walls | |
CN113361099B (zh) | 一种裂隙岩体模拟方法及系统 | |
Tzamtzis et al. | A 3d model for non-linear ‘microscopic’fe analysis of masonry structures | |
Barrar | Structural optimization using the principle of virtual work and an analytical study on metal buildings | |
WO2023058148A1 (ja) | 地盤のn値推定方法及びプログラム | |
Roca et al. | Micro-mechanical modeling of masonry-Parametric study | |
El-Hafez et al. | Direct design of reinforced concrete skew slabs | |
Ackroyd | Nonlinear inelastic analysis of flexibly-connected steel frames | |
Oumer et al. | Assessment of distributed non-linear fiber models | |
Jansen | Influence of creep on concrete beam deflection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |