Mesoscopic simulation of coal cylinders confined by CFRP with different sizes
-
摘要:
为探究碳纤维增强聚合物(CFRP)布约束煤圆柱的尺寸效应规律,利用三维PFC-FLAC构建碳纤维增强聚合物(CFRP)布约束煤圆柱细观模型。通过已有试验数据验证了该模型的有效性,并在此基础上开展0、1、2、3层CFRP布对不同尺寸(50、100、150 mm)煤圆柱尺寸效应的研究;得到了不同层数(L)下抗压强度演化规律的函数关系,以及不同层数下弹性模量演化规律的函数关系;并建立了CFRP布约束煤圆柱的修正Richart强度分析模型。结果表明:在CFRP布层数一定时,不同尺寸CFRP布约束煤圆柱的抗压强度,具有明显的尺寸效应,CFRP布约束煤圆柱的相对峰值强度增强比,随尺寸的增大而增大;直径50 mm煤圆柱随CFRP层数增加弹性模量逐渐增大,直径100、150 mm煤圆柱在0≤L≤3层弹性模量先增加后减小,具有明显的层效应规律;裂隙累积量随CFRP层数的增加而增大。
Abstract:In order to study the size effect law of carbon fiber reinforced polymer (CFRP) constrained coal cylinders, a numerical model of coal circular confined by carbon fiber reinforced polymer (CFRP) sheet is established by a three-dimensional PFC-FLAC coupling method. The validity of the model was verified with experimental data. On this basis, the study of the size effect of 0, 1, 2 and 3 layers of CFRP sheet on different sizes (50 mm, 100 mm, 150 mm) of coal cylinders were studied; the functional relationship of the evolution law of compressive strength under different layer number (L) is obtained; the functional relationship of elastic modulus evolution under different number of layers is obtained; a modified Richart strength analysis model of coal cylinder confined by CFRP sheet is established. The results show that when the number of CFRP sheets is constant, the compressive strength of coal cylinders is constrained by different sizes of CFRP sheets, and the size effect is obvious; the relative peak strength enhancement ratio of coal cylinder constrained by CFRP sheet increases with the increase of the size. The elastic modulus of a coal cylinder with a diameter of 50 mm increases with the number of CFRP layers. The elastic modulus of a coal cylinder with a diameter of 100 mm and 150 mm increases first and then decreases at 0≤L≤3 layers, showing an obvious layer effect law. The accumulation of cracks increases with the increase of the number of CFRP layers.
-
我国金属矿、煤矿等主要依靠井工开采,当发生冒顶、透水等井下灾害事故造成人员被困时,由于巷道坍塌严重,次生灾害频发,导致井下救援受限[1-3]。地面救援是一种首先施工大直径应急救援井,然后以提升舱为载体,使用提升设备将被困人员从井下提升至地面的救援方式,具有速度快、可保证救援人员安全等特点,在国内外已有多个成功应用实例[4-5]。国外矿难救援的成功,特别是2010年智利矿难救援的壮举,给予我国矿山救援装备制造企业很大的启发,国内外先后生产了不同规格型号的救援提升舱[6-8],但普遍存在功能单一、对被救人员和设备的保护性不足、缺乏指导救援系统及方法等问题。
中煤科工西安研究院(集团)有限公司在“十三五”国家重点研发计划课题“地面应急救援井安全提升装备”的资助下,研制了多功能柔性救援提升舱,从救援指挥系统、柔性转动、导向、缓冲、钢丝绳与舱体分动、断绳保护、救援提升方法等多方面提出了解决方案,以加强对人员和设备的保护,提高地面大直径钻孔提升救援的成功率、安全性、规范性。
1. 柔性救援能提升舱总体方案
经过充分调研,对煤矿巷道结构进行实地考察,并根据人机工程学原理,确定了救生舱的设计参数。矿用应急救援提升舱外径550 mm,内部净高1 930 mm,每次可载1名人员,要求舱体对钻孔的适应性、通过性强,同时实现空气监测、二次逃生、实时通信等功能,二次逃生高度800 m。提升舱主要包括上导向、上舱段、柔性转动装置、人员舱、下舱段、下导向、分动装置、断绳保护装置、上侧缓冲、下侧缓冲、底部缓冲、救援指挥系统等方面。柔性救援提升舱总体方案如图1所示。
提升舱整体结构为“子弹形”,从上到下的连接顺序为上导向、上舱段、柔性转动装置、人员舱、下舱段、下导向。分动装置安装于上导向与上舱段的结合面,上侧缓冲和断绳保护装置安装于上舱段内部,下侧缓冲安装于下舱段内部,底部缓冲安装于下舱段与下导向的结合面。上导向、下导向均为锥形结构,有利于提升舱下放和提升时进入救援孔[9]。救援钢丝绳同时具备提升舱体和救援指挥系统通信2方面的功能。
柔性救援提升舱实现可视化救援指挥、救生舱柔性转动、通信钢丝绳防扭转等功能,以加强提升舱的通过性和对人员和设备的保护,提高地面大直径钻孔提升救援的成功率、安全性、规范性。
2. 分动装置
分动装置实现了通信钢丝绳与提升舱的分动,避免救援过程中提升舱带动通信钢丝绳旋转,干扰内部信号传输,损坏钢丝绳编织结构。不同于普通电信号滑环的分动结构[10],分动装置实现了通信钢丝绳与提升舱转动分动的同时,承受钢丝绳轴向拉力,进而实现了通信钢丝绳提升和通信的双重功能。分动装置结构图如图2所示。
悬挂轴为阶梯式通轴结构,轴内部为中空,上部加工有紧固螺纹,悬挂轴下端法兰安装于上导向和上舱段的结合面,悬挂轴上端小直径段与上导向的固定板配合安装(悬挂轴安装时应从上导向底部向上安装);垫板1安装在悬挂台肩上,止推轴承1的2个端面分别与垫板1和接头下凹槽端面接触,止推轴承2的2个端面分别与接头上凹槽端面和垫板2接触;锁紧螺母通过悬挂轴上的紧固螺纹将止推轴承、止推轴承、垫板1、垫板2、接头压紧至悬挂轴台肩上,实现接头与悬挂轴的分动(一方不随另一方转动)。滑环安装座固定在悬挂轴上端面,防爆滑环安装在滑环安装座中,滑环下部导线通过悬挂轴中空孔引至上舱段中的信号转换器,防爆滑环的上部导线通过拉板固定在接头上,滑环上部与下部通过防爆电刷分动;防爆滑环上部导线与通信钢丝绳内部电缆连接,将提升舱内的电信号通过通信钢丝绳传输至地面终端;销轴安装在接头上。通信钢丝绳通过锁具与销轴连接,实现提升舱的提升和下放。防爆滑环出线连接救援指挥系统在救生舱上安装的部分,实现通信系统的连通。
3. 柔性转动装置
地面大直径钻孔过程中,可能会出现钻孔偏斜、弯曲等情况,提升舱面临无法通过的风险[11]。通过在上舱段和人员舱中间设计柔性转动装置,使提升舱遭遇钻孔偏斜、弯曲时转动一定角度,增加其钻孔适应性和通过性。舱体小角度弯曲示意图如图3所示,柔性转动装置结构图如图4所示。
柔性转动装置包括柔性旋转座、旋转体、限转块、油杯。柔性旋转内部加工为球形面,旋转体外部加工为球形面,两者球形面的直径相同,柔性旋转座与旋转体通过内外球面配合,两者可绕球心在任意方向转动。柔性旋转座与上舱段连接,旋转体穿过柔性旋转座后通过可拆卸方式与人员舱连接;柔性旋转座圆周方向对称加工有2个方形凹台,限转块安装于凹台之上,限制柔性旋转与旋转体之间绕球心旋转的角度,避免提升舱过渡弯曲及旋转过程中剪断下舱段底部监控器的通信电缆,同时可根据钻孔轴线弯曲程度,安装不同大小的限转块,达到改变提升舱柔性角度α。
柔性转动装置通过油杯注射润滑液体,通过旋转体上的润滑孔道进入柔性旋转座与旋转体的配合球面,实现润滑。
柔性转动装置不同于机器人及仿生结构的柔性结构,主要承受轴向拉力的作用,且内部大直径中空,实现上舱段和下舱段的大空间连通。
提升救援过程中提升舱主要受力来源于旋转体及以下舱段的设备和人员质量,旋转体承受向下拉力,旋转座为壳体结构,可能出现外扩变形,导致旋转体滑落,出现人员坠落的二次事故,故重点对柔性旋转座与旋转体之间的接触球面进行受力分析。柔性转动装置应变云图如图5所示。
通过Workbench有限元分析可以看出,在静态加载10 000 N(人员质量及设备质量的6倍)情况下,最大变形量为0.024 mm,出现在内球铰体与下舱体连接块连接部位,远小于旋转件外直径与设备舱壳体内直径之间的距离,柔性旋转座与旋转体之间的接触球面强度满足迁都要求。
4. 缓冲装置
缓冲装置主要分为上侧缓冲、下侧缓冲、底部缓冲3个部分。上部缓冲安装于上舱段内部,下侧缓冲安装于下舱段内部,底部缓冲安装于下舱段与下导向的结合面。
上、下侧缓冲装置结构图如图6所示。上侧缓冲和下侧缓冲结构相同,均由6组缓冲体组成,在上舱段圆周方向均布,缓冲体由底座、弹簧、滚轮组成,自然状态下缓冲体的弹簧为伸长状态,6个滚轮外圆周构成的圆形直径大于救援钻孔直径;在提升舱下放、提升过程中,滚轮压缩弹簧,使6个滚轮外圆周构成的圆形直径等于救援钻孔直径,弹簧反作用力将滚轮压紧在救援钻孔孔壁上,通过弹簧压缩量随救援钻孔直径的变化而变化,由此达到下放、提升过程中对于救生舱和舱内人员的缓冲保护效果。
底部缓冲装置结构图如图7所示。底部缓冲由4组弹簧、套筒、导向销、底板组成,导向销穿过下导向底板的导向孔,一端与底板固定,底板和导向销整体可沿导向孔轴线上下自由滑动,上套筒、弹簧安装在导向销上,弹簧自然状态的长度大于套筒长度,提升舱到达人员被困巷道或提升至地面接触地面的过程中,底板最先接触地面压缩弹簧,至套筒支撑底板,实现提升舱着陆缓冲。其中,套筒长度应大于下舱段底部监控器的外露长度,避免着陆过程中底板损坏下舱段底部监控器。
5. 救援指挥系统
救援指挥系统如图8所示。救援指挥系统的上舱段顶部监控器、人员舱内部监控器、下舱段底部监控器分别对提升舱上部救援钻孔情况、提升舱内部人员情况、提升舱下部救援钻孔情况进行监控,气体浓度监测仪对舱内气体浓度进行监控,人员舱内部语音对讲器实现地面与舱内人员实时双向联系,气体浓度信号、监控视频和语音信号均通过信号转换器经通信钢丝绳传输至地面终端。地面终端通过通信监测软件界面对所有通信信号进行集中显示和指挥控制[12-15]。
6. 提升舱人机工程分析
提升救援过程中,人员舱内的人员在一定程度上会受到灾害带来的影响而受伤,从而出现行动不便和无力支撑等情况。因此,获救人员进入人员舱内,需要提供一定的身体支撑,对人员舱内部结构进行人机工程分析也是必不可少。
6.1 入舱情况分析
对人员入舱情况进行分析,根据救生舱结构尺寸,入舱门槛高600 mm,对人员入舱有阻碍,尤其对人体尺寸百分位数为P50以下的人群入舱行为产生明显的障碍。
为了解决被救人员入舱阻碍的问题,需要在人员舱舱门内部设置入舱扶手,人员入舱扶手高度应该设置在离地1 400~1 800 mm处,这样可对人体尺寸百分位数为P5~P95的人群给予借力的帮助,也方便被救人员顺利进入舱内。人员舱入舱扶手高度图如图9所示,
6.2 人员姿态情况分析
在被救人员入舱后,对其在人员舱内的姿态进行分析,加装一些舱内辅助人机工程组件,保障被救人员在人员舱内能够处于一个比较舒适的姿态。依、靠、扶、握可使保持的姿态省力、放松,也有助于抵抗救生舱舱体在救援提升过程中由于舱体晃动产生的心理不安因素,帮助安抚受救人员的紧张情绪。
一般人员舱内被救人员的舱内姿态分为靠坐握和站趴扶2种形式,舱内人员靠坐握姿态分析图如图10所示,舱内人员站趴扶姿态分析图如图11所示。
由图10可以看出,靠坐握姿态的舒适度高(绿色),可保持姿态较长时间而不会疲劳,可预测为被救援人员较长时间保持的姿态。由图11可以看出,站趴扶姿态舒适度略低(黄色),长时间坚持会产生肌肉疲劳,可预测为被救援人员暂时变换姿态短暂使用。
根据对人员舱内被救人员的姿态分析,选取靠坐握姿态为主要姿态依据,设计人机工程组件。需要在被救人员的臀部部位设计坐姿的支点,提供尺寸合适的坐姿点和扶手,应以人体尺寸百分位数为P50的人体模型双膝微屈姿态为基础原型,支点的高度740 mm,宽5 mm,形态圆润,防止碰撞对人体造成伤害,并且扶手高度应设置为890~1 160 mm。人员坐姿支点设计图如如图12所示。
当人员舱内被救人员坐姿调整为站立姿态时,提供尺寸合适的立姿扶手,应以人体尺寸百分位数为P50的人体模型身体前倾为基础原型,手扶舱壁的高度确定为1 515 mm,外形圆润,也要防止碰撞对人体造成伤害。人员立姿扶手支点设计图如图13所示。
6.3 安全防护情况分析
舱内需设置防撞装置,从而避免在救援提升过程中,被救人员在舱内发生碰撞造成伤害。防撞装置主要针对头部进行设置,沿舱壁设置橡胶条防撞,尺寸下线以人体尺寸百分位数为P5的人体尺寸模型坐姿状态下,受到晃动后头部撞击舱壁的位置为确定高度,上限以人体尺寸百分位数为P99的人体模型立姿头顶位置为确定高度。
舱体内壁需要设置安全挂钩,可将被救人员身上的安全带固定在安全挂钩上,也能起到在舱体碰撞过程中固定被救人员避免冲撞的危险。安全挂钩的安装高度应该以人体尺寸百分位数为P99的人体模型尺寸加上鞋与安全帽的修正尺寸为依据,确定为1 920 mm。安全挂钩安装示意图如图14所示。
7. 测试及试验
柔性救援提升舱与救援提升车共同构成了提升救援成套装备,在室内和施工现场分别对舱体弯曲角度和综合性能进行了测试和试验。
1)弯曲角度测试。将提升舱顶端吊起,底端着地,保持上下舱体弯曲至最大值,利用数显倾角仪测量上下舱体倾角度数得到舱体相对的弯曲角度。上下舱体2次测量的相对弯曲角度分别为2.2°和2.6°,均大于2°的指标,满足设计指标要求。
2)现场综合试验。通过宁夏梅花井矿施工的大直径救援井,对提升救援成套装备进行了综合检验。通过空载、重载、载人登试验验证了提升舱的综合性能。具体为:①提升舱内径满足人员正常或辅助站立,上下舱体间可摆动,配备有井径测量、气体测量和音视频通讯装置,能在580 mm以上的井筒内顺利通过,视频音频清晰度高,能清晰看到救援舱内和井壁内的运行变化,内壁圆整,套管与裸井段接缝处无明显凸起,提升舱可顺利通过井下通道;②在线气体测量功能正常,能实时显示气体浓度,从现场软件记录数据来看,测试孔段CO、H2S 2种气体体积分数为0,O2体积分数在20%~21.2%之间、CH4体积分数在0.2%~0.3%之间、CO2体积分数在0.01%~0.02%之间变化,大气压差在107~110 Pa 之间变化,温度在30~34 ℃之间变化;③提升舱在运行速度突变及载重情况下,结构稳固,音视频数据传输正常,试验验证了提升舱重载条件下的可靠性。
8. 结 语
研制了多功能矿用救援柔性提升舱。多功能救援柔性提升舱从救援指挥系统、柔性转动、导向、缓冲、钢丝绳与舱体分动等多方面提出了解决手段。实现了可视化救援指挥、通信钢丝绳防扭转、救生舱柔性转动等功能,提高地面大直径钻孔提升救援的成功率、安全性、规范性。
1)分动装置实现了通信钢丝绳与提升舱的分动,避免救援过程中提升舱带动通信钢丝绳旋转,干扰内部信号传输,损坏钢丝绳编织结构,为实现通信钢丝绳提升和通信的双重功能提供了技术支撑。
2)柔性转动装置实现了上下舱体的弯曲,增强了提升舱体对曲率较大、有局部凸起等复杂井壁的适应性,以加强提升舱的通过性。
3)通过人机工程学相关理论和技术,开发了人员舱内辅助入舱、靠坐握姿态等人机工程组件,对被救人员进行全方位保护。
4)基于井径测量、气体测量和音视频通讯装置的通信系统则从监测井壁、气体及环境参数、被救人员等多方面实现可视化救援,保证救援过程的安全。
-
表 1 煤圆柱细观参数
Table 1 Mesoscopic parameters of coal cylinder
参数 符号 数值 颗粒半径/mm r 1~1.9 摩擦因数 μ 0.1 颗粒间刚度比 kn/ks 0.08 平行黏结强度/MPa pb_coh 1×108 平行黏结抗拉强度/MPa pb_ten 1×108 表 2 CFRP布模型参数
Table 2 Model parameters of CFRP sheets
表 3 细观模拟方案
Table 3 Mesoscopic simulation schemes
试件名称 截面尺寸/
(mm×mm)CFRP层数 加载速率/
(mm·min−1)B-C-M 50×100 0 0.1 B-3-CC-M 50×100 3 C-C-M 100×200 0 C-1-CC-M 100×200 1 C-2-CC-M 100×200 2 C-3-CC-M 100×200 3 D-C-M 150×300 0 D-1-C-M 150×300 1 D-2-C-M 150×300 2 D-3-C-M 150×300 3 表 4 截面尺寸试验结果统计
Table 4 Statistics of section size test results
试件名称 截面尺寸/(mm×mm) CFRP层数 峰值强度/MPa B-C-M 50×100 0 24.771 C-C-M 100×200 0 24.467 D-C-M 150×300 0 24.186 B-1-CC-M 50×100 1 130.538 C-2-CC-M 100×200 2 175.032 D-3-CC-M 150×300 3 223.209 注:当试件的尺寸分别为50 mm×100 mm、100 mm×200 mm、150 mm×300 mm时,相应的CFRP布约束层数设为1、2、3层,能够确保未约束煤圆柱试件受到的侧向约束效果一致,即约束比相同。 表 5 不同尺寸CFRP布约束煤圆柱抗压强度尺寸效应拟合关系式具体参数
Table 5 Specific parameters of the fitting relation for the size effect of coal cylinder compressive strength constrained by different sizes of CFRP sheet
层数 aP bP cP 0 20.753 4.347 635.542 1 136.409 −4.927 285.423 2 170.843 48.987 −40.665 3 222.917 621.692 −19.571 表 6 不同尺寸CFRP布约束煤圆柱弹性模量尺寸效应具体参数
Table 6 Specific parameters of size effect of elastic modulus of coal cylinder constrained by CFRP sheet with different sizes
层数 aE bE cE 0 24.259 27.027 575.086 1 41.158 46.896 423.625 2 53.351 59.114 576.079 3 151.071 157.108 1834.053 表 7 CFRP布约束煤圆柱强度修正模型参数
Table 7 Parameters of strength modified models of coal cylinder restrained by CFRP sheet
参数 修正 Richart 修正 HoekBrown α β γ ψ η θ 数值 2.993 2.564 0.542 −1 1.991×106 0.883 -
[1] 李庆文,高森林,黄筱. 基于加载速率效应的煤样损伤演化研究[J]. 煤矿安全,2023,54(11):105−115. LI Qingwen, GAO Senlin, HUANG Xiao. Study on damage evolution of coal samples based on loading rate effect[J]. Safety in Coal Mines, 2023, 54(11): 105−115.
[2] 纪孙航,王文达,鲜威. CFRP加固火灾作用后圆钢管混凝土构件的侧向撞击性能研究[J]. 工程力学,2021,38(8):178−191. JI Sunhang, WANG Wenda, XIAN Wei. Lateral impact behavior of CFRP-reinforced circular concrete-filled steel tubular members[J]. Engineering Mechanics, 2021, 38(8): 178−191.
[3] DAS A J, MANDAL P K, GHOSH C N, et al. Extraction of locked-up coal by strengthening of rib pillars with FRP-A comparative study through numerical modelling[J]. International Journal of Mining Science and Technology, 2017, 27(2): 261−267. doi: 10.1016/j.ijmst.2017.01.024
[4] 李庆文,禹萌萌,高森林,等. 加载速率对碳纤维布被动约束煤能量演化影响研究[J/OL]. 煤炭学报,2023:1-14. [2023-06-20]. DOI: 10.13225/j.cnki.jccs.2023.0238 LI Qingwen, YU Mengmeng, GAO Senlin, et al. The effect of loading rate on energy evolution of coal confined passively by CFRP sheets[J/OL]. Journal of China Coal Society, 2023: 1-14. [2023-06-20]. DOI: 10.13225/j.cnki.jccs.2023.0238.
[5] XIA Z, YAO Q, LI X, et al. Acoustic emission characteristics and energy mechanism of CFRP-jacketed coal specimens under uniaxial compression[J]. Construction and Building Materials, 2022, 342: 127936 doi: 10.1016/j.conbuildmat.2022.127936
[6] HUANG Y H, YANG S Q, ZHAO J. Three-dimensional numerical simulation on triaxial failure mechanical behavior of rock-like specimen containing two unparallel fissures[J]. Rock Mechanics and Rock Engineering, 2016, 49(12): 4711−4729. doi: 10.1007/s00603-016-1081-2
[7] YAO W, CAI Y Y, YUJ, et al. Experimental and numerical study on mechanical and cracking behaviors of flawed granite under triaxial compression[J]. Measurement, 2019, 145: 573−582. doi: 10.1016/j.measurement.2019.03.035
[8] YANG S Q, YIN P, ZHANG Y, et al. Failure behavior and crack evolution mechanism of a non-persistent jointed rock mass containing a circular hole[J]. International Journal of Rock Mechanics and Mining Sciences, 2019, 114: 101−121. doi: 10.1016/j.ijrmms.2018.12.017
[9] 李庆文,禹萌萌,刘艺伟,等. GFRP布被动约束标准煤矸石混凝土圆柱轴压性能细观模拟[J]. 硅酸盐通报,2023(7):2458−2471. LI Qingwen, YU Mengmeng, LIU Yiwei, et al. Mesoscopic simulation on axial compression performance of GFRP passively confined standard coal gangue concrete circular-columns[J]. Bulletin of the Chinese Ceramic Society, 2023(7): 2458−2471.
[10] SHI Can, ZHAO Chunfa, ZHANG Xu, et al. Coupled discrete-continuum approach for railway ballast track and subgrade macro-meso analysis[J]. International Journal of Pavement Engineering, 2020, 22(13): 1−16.
[11] 何涛,王丽. 基于颗粒离散元模型的煤岩组合体损伤演化规律研究[J]. 煤矿安全,2018,49(7):205−208. HE Tao, WANG Li. Research on damage evolution laws of rock-coal-rock combination based on particle discrete element model[J]. Safety in Coal Mines, 2018, 49(7): 205−208.
[12] 黄世斌,邹辰皓,王家全,等. 基于离散-连续耦合的加筋包裹碎石桩承载机理分析[J]. 广西大学学报(自然科学版),2022,47(1):13−21. HUANG Shibin, ZOU Chenhao, WANG Jiaquan, et al. Load-bearing mechanism analysis of geosynthetic encased stone columns based on discrete-continuous coupling[J]. Journal of Guangxi University (Natural Science Edition, 2022, 47(1): 13−21.
[13] 李俊孟. 矸石固废充填材料承载压缩三维组构时空演化及其透明化表征[D]. 徐州:中国矿业大学,2020. [14] TAN Xin, HU Zhengbo, CHEN Changfu, et al. 3D DEM-FDM coupled analysis of the behavior of an isolated geogrid-encased stone column under axial loading[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(6). https://doi.org/10.1061/(ASCE)GT.1943-5606.0002516.
[15] 李庆文,曾杏钢,张向东,等. 碳纤维布层数对煤圆柱力学特性影响的细观研究[J]. 煤炭科学技术,2023(8):73-85. LI Qingwen, ZENG Xinggang, ZHANG Xiangdong, et al. Mesoscopic study on the effect of CFRP layers on the mechanical properties of coal circular-columns[J]. Coal Science and Technology, 2023(8): 73-85.
[16] 石崇,褚卫江,张一平,等. 岩土工程颗粒流数值模拟技术应用案例[M]. 北京:中国建筑工业出版社,2023. [17] 石崇,张强,王盛年. 颗粒流(PFC5.0)数值模拟技术及应用[M]. 北京:中国建筑工业出版社,2018. [18] 李庆文,胡露露,曹行,等. CFRP布均匀约束煤圆柱轴压性能[J]. 复合材料学报,2022(11):5611-5624. LI Qingwen, HU Lulu, CAO Hang, et al. Axial compressive behavior of CFRP uniformly wrapped coal in circular columns[J]. Acta Materiae Compositae Sinica, 2022(11): 5611-5624.
[19] HE P F, KULA TILAKE P H S W, YANG X X, et al. Detailed comparison of nine intact rock failure criteria using poly axial intact coal strength data obtained through PFC3D simulations[J]. Acta Geotechnica, 2018, 13(2): 419−445.
[20] ZHANG L, REN T, LI X C, etal. Acousticemission, damage and cracking evolution of intact coal under compressive loads: experimental and discrete element modelling[J]. Engineering Fracture Mechanics, 2021, 252: 107690. doi: 10.1016/j.engfracmech.2021.107690
[21] 王刚,王锐,武猛猛,等. 渗透压-应力耦合作用下煤体常规三轴试验的颗粒流模拟[J]. 岩土力学,2016,37(S1):537−546. WANG Gang, WANG Rui, WU Mengmeng, et al. Simulation of conventional triaxial test on coal under hydro-mechanical coupling by particle flow code[J]. Rock and Soil Mechanics, 2016, 37(S1): 537−546.
[22] 谭鑫,曹明,冯龙健,等. 土工织物包裹碎石桩力学特性的数值模拟研究[J]. 中国公路学报,2020,33(9):136−145. TAN Xin, CAO Ming, FENG Longjian, et al. Numerical study on mechanical behaviors of geotextile wrapped stone column[J]. Chinese Journal of Highways, 2020, 33(9): 136−145.
[23] 郭润兰,范雅琼,王广书,等. 基于PFC3D的机床床身用树脂矿物复合材料损伤性能细观研究[J]. 复合材料学报,2022,39(2):834−844. GUO Runlan, FAN Yaqiong, WANG Guangshu, et al. Meso-scale study on damage performance of resin mineral composite material for machine tool bed based on PFC3D[J]. Journal of Composites, 2022, 39(2): 834−844.
[24] EID R, PAULTRE P. Plasticity-based model for circular concrete columns confined with fibre-composite sheets[J]. Engineering Structures, 2007, 29(12): 3301−3311. doi: 10.1016/j.engstruct.2007.09.005
[25] EID R, DANCYGIER A N, PAULTRE P. Elastoplastic confinement model for circular concrete columns[J]. Journal of Structural Engineering, 2007, 133(12): 1821−1831. doi: 10.1061/(ASCE)0733-9445(2007)133:12(1821)
[26] 梁正召,张永彬,唐世斌,等. 岩体尺寸效应及其特征参数计算[J]. 岩石力学与工程学报,2013,32(6):1157−1166. LIANG Zhengchao, ZHANG Yongbin, TANG Shibin, et al. Size effect of rock messes and associated representative element properties[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(6): 1157−1166.
[27] 陈卫忠,杨建平,邹喜德,等. 裂隙岩体宏观力学参数研究[J]. 岩石力学与工程学报,2008,27(8):1569−1575. CHEN Weizhong, YANG Jianping, ZOU Xide, et al. Research on macromechanical parameters of fractured rock masses[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(8): 1569−1575.
[28] SHI Xudong, FENG Guorui, BAI Jinwen, et al. Numerical and experimental study of the mechanical behaviour for FRP-wrapped cement mortar-coal composite disc[J]. Computers and Geotechnics, 2023, 158: 105373. doi: 10.1016/j.compgeo.2023.105373
[29] LIM J C, OZBAKKALOGLUT. Lateral strain-to-axial strain relationship of confined concrete[J]. Journal of Structural Engineering, 2015, 141(5): 04014141. doi: 10.1061/(ASCE)ST.1943-541X.0001094
[30] FERROTTO M F, FISCHERO, CAVALERIL. Analysis-oriented stress–strain model of CRFP-confined circular concrete columns with applied preload[J]. Materials and Structures, 2018, 51(2): 116.
[31] WU Yufei, ZHOU Yingwu. Unified strength model based on Hoek-Brown failure criterion for circular and square concrete columns confined by FRP[J]. Journal of Composites for Construction, 2010, 14(2): 175−184. doi: 10.1061/(ASCE)CC.1943-5614.0000062
[32] ZHANG Yang, LU Zhifang, CAO Yugui. Unified strength model based on the Hoek-Brown failure criterion for fibre-reinforced polymer-confined pre-damaged concrete columns with circular and square cross sections[J]. Journal of Central South University, 2020, 27(12): 3807−3820. doi: 10.1007/s11771-020-4563-z