CN115358099A - 一种全长黏结锚索的数值模拟方法及系统 - Google Patents

一种全长黏结锚索的数值模拟方法及系统 Download PDF

Info

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
Application number
CN202211286678.1A
Other languages
English (en)
Other versions
CN115358099B (zh
Inventor
陈见行
王世纪
胡松松
张汉
孙国新
张俊文
杨胜利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202211286678.1A priority Critical patent/CN115358099B/zh
Publication of CN115358099A publication Critical patent/CN115358099A/zh
Application granted granted Critical
Publication of CN115358099B publication Critical patent/CN115358099B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force 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对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:
Figure DEST_PATH_IMAGE001
;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;
Figure DEST_PATH_IMAGE002
为锚固界面剪切强度;
Figure DEST_PATH_IMAGE003
为锚索直径;所述第二公式为:
Figure DEST_PATH_IMAGE004
;其中,c为非线性剪切软化系数;所述第三公式为
Figure 100002_DEST_PATH_IMAGE005
;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
根据本发明的另一方面,提供了一种全长黏结锚索的数值模拟系统。
该全长黏结锚索的数值模拟系统,包括:
模型创建模块,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
其中,所述模型创建模块在创建三维数值模拟模型时,可根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;并定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置。
可选的,所述模型配置的计算模式为大变形计算模式。
其中,所述模型创建模块在所述三维数值模拟模型中,创建锚索时,可根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。
可选的,每个锚固构建的长度需小于等于200mm。
可选的,所述基本力学性质参数包括:锚固界面剪切刚度、锚固界面法向刚度、锚固界面沿剪切方向的摩擦角、锚固开关、锚索杨氏模量、锚索泊松比、锚索横截面积、锚索周长、锚索抗拉极限、锚索极惯性矩、锚索沿Y轴的二次惯性矩和/或锚索沿Z轴的二次惯性矩。
其中,所述关系定义模块在定义锚索与锚固剂间锚固界面的剪切滑移关系是包括:定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。
其中,所述函数处理模块在创建锚固界面剪切行为FISH函数时,可查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:
Figure DEST_PATH_IMAGE006
;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;
Figure 100002_DEST_PATH_IMAGE007
为锚固界面剪切强度;
Figure DEST_PATH_IMAGE008
为锚索直径;所述第二公式为:
Figure DEST_PATH_IMAGE009
;其中,c为非线性剪切软化系数;所述第三公式为
Figure DEST_PATH_IMAGE010
;其中,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)进行计算;而锚固界面法向刚度和锚固界面沿剪切方向的摩擦角统一设置为零。
具体的,公式(1)为
Figure DEST_PATH_IMAGE011
;公式(2)为
Figure DEST_PATH_IMAGE012
;公式(3)为
Figure DEST_PATH_IMAGE013
;公式(4)为
Figure DEST_PATH_IMAGE014
其中,S k为锚固界面剪切刚度;τ max为锚固界面剪切强度;d 1为锚固界面剪切强度对应的剪切滑移量;D b为锚索直径;I为锚索极惯性矩;I y为锚索沿Y轴的二次惯性矩;I z为锚索沿Z轴的二次惯性矩。
在具体应用时,在定义锚索与锚固剂间锚固界面的剪切滑移关系时,可定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。而创建锚固界面剪切行为FISH函数时,可按照如下逻辑关系进行计算:查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:
Figure DEST_PATH_IMAGE015
;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;
Figure DEST_PATH_IMAGE016
为锚固界面剪切强度;
Figure DEST_PATH_IMAGE017
为锚索直径;所述第二公式为:
Figure DEST_PATH_IMAGE018
;其中,c为非线性剪切软化系数;所述第三公式为
Figure DEST_PATH_IMAGE019
;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
对应的,如图2所示,根据本发明实施例的一种全长黏结锚索的数值模拟系统,包括:
模型创建模块201,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块203,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块205,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块207,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
在具体应用时,所述模型创建模块201在创建三维数值模拟模型时,可根据实际工况,在FLAC3D软件内创建三维数值模拟模型,确定三维数值模拟模型的长、宽和高以及沿着三维数值模拟模型长、宽和高方向单元体的数量;并定义单元体的本构模型和材料参数以及三维数值模拟模型的边界条件、初始条件和模型配置,所述模型配置的计算模式为大变形计算模式;所述模型创建模块201在所述三维数值模拟模型中,创建锚索时,可根据实际工况,确定桩结构单元内锚固构建的数量;根据桩结构单元内锚固构建的数量,在所述三维数值模拟模型内利用桩结构单元创建锚索;并根据实际工况中全长黏结锚索的基本情况,输入锚索的基本力学性质参数;其中,锚索长度与实际工况中的锚索长度相同。
而所述关系定义模块203在定义锚索与锚固剂间锚固界面的剪切滑移关系是包括:定义锚固界面的剪切强度、锚固界面剪切强度对应的剪切滑移量、非线性剪切软化系数、锚固界面残余剪切强度对应的剪切滑移量。所述函数处理模块205在创建锚固界面剪切行为FISH函数时,可查找桩结构单元的头指针并将其赋值给指针变量pnt;利用所述指针变量pnt对桩结构单元的全部锚固构建进行遍历;每遍历一个锚固构建时,取出该锚固构建前端节点的剪切滑移量,将其取绝对值后赋值给剪切滑移变量slip;判断剪切滑移变量slip是否小于等于锚固界面剪切强度对应的剪切滑移量;在判断结果为是的情况下,根据预定第一公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,继续判断剪切滑移变量slip是否小于等于锚固界面残余剪切强度对应的剪切滑移量;并在判断结果为是的情况下,根据预定第二公式计算锚固构建单位长度的剪切强度输入参数,并将计算结果赋值给指针变量pnt对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
其中,所述第一公式为:
Figure DEST_PATH_IMAGE020
;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;
Figure DEST_PATH_IMAGE021
为锚固界面剪切强度;
Figure DEST_PATH_IMAGE022
为锚索直径;所述第二公式为:
Figure DEST_PATH_IMAGE023
;其中,c为非线性剪切软化系数;所述第三公式为
Figure DEST_PATH_IMAGE024
;其中,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对应的锚固构建单位长度的剪切强度输入参数;在判断结果为否的情况下,根据预定第三公式计算锚固构建单位长度的剪切强度输入参数。
9.根据权利要求8所述的一种全长黏结锚索的数值模拟方法,其特征在于,所述第一公式为:
Figure 240959DEST_PATH_IMAGE001
;其中,F s为锚固构建单位长度的剪切强度输入参数;d 1为锚固界面剪切强度对应的剪切滑移量;d为锚固界面的剪切滑移量;
Figure 20696DEST_PATH_IMAGE003
为锚固界面剪切强度;
Figure DEST_PATH_IMAGE005
为锚索直径;所述第二公式为:
Figure 415905DEST_PATH_IMAGE006
;其中,c为非线性剪切软化系数;所述第三公式为
Figure DEST_PATH_IMAGE007
;其中,d 2为锚固界面残余剪切强度对应的剪切滑移量。
10.一种全长黏结锚索的数值模拟系统,其特征在于,包括:
模型创建模块,用于创建三维数值模拟模型,并在所述三维数值模拟模型中,创建锚索;
关系定义模块,用于定义锚索与锚固剂间锚固界面的剪切滑移关系;
函数处理模块,用于创建锚固界面剪切行为FISH函数,并将定义的锚固界面剪切滑移关系输入到该FISH函数中;
速度测试模块,用于在锚索拉拔端施加恒定的拉拔速度;记录锚索拉拔端的拉拔力和拉拔位移并绘制拉拔性能曲线。
CN202211286678.1A 2022-10-20 2022-10-20 一种全长黏结锚索的数值模拟方法及系统 Active CN115358099B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 中国矿业大学 支护构件与硐室围岩接触及相互作用的模拟方法和系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
JIANHANG CHEN 等: "Studying the performance of fully encapsulated rock bolts with modified structural elements", 《INT J COAL SCI TECHNOL》 *
言志信 等: "锚固参数对地震作用下岩质边坡锚固界面剪切作用影响的数值分析", 《煤炭学报》 *
郭金刚 等: "基于残余剪切强度的全长黏结锚杆拉拔模拟", 《岩土力学》 *

Cited By (6)

* Cited by examiner, † Cited by third party
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