Bearing characteristics of ultra-long piles based on distributed fiber optic sensing technology
-
摘要: 为探究不同地质条件下超长桩承载变形特性及数值模拟中摩擦系数对桩-土相互作用的影响,以深圳湾超级总部某项目系列试桩为背景,采用FBG密集分布式光纤传感技术进行基桩内力测试,并采用ABAQUS有限元软件进行数值模拟,对比分析桩、土界面采用不同摩擦系数取值方式时土层平均侧阻力模拟结果。结果表明:①灌注桩承载变形特性与桩长和持力层性质密切相关,超长桩轴力沿深度方向基本呈线性递减,端阻比明显减小;持力层为强风化花岗岩层时,试桩轴力沿深度方向呈“大肚形”递减,端阻比低于17.1%。②数值模拟中摩擦系数取值方式会对土体侧阻力模拟结果产生较大影响,根据土层性质差异,对桩土界面采用多个摩擦系数能够提高模拟结果准确性。
-
关键词:
- FBG密集分布式光纤传感技术 /
- 超长灌注桩 /
- 承载变形特性 /
- 数值模拟 /
- 侧阻力
Abstract: Based on a series of pile foundation tests from a project in the Super Headquarters Base in Shenzhen Bay, the load deformation characteristics of super-long piles under different geological conditions and the impact of friction coefficient on pile-soil interaction in numerical simulations are investigated. The FBG (fiber bragg grating) densely distributed fiber optic sensing technology is employed for internal force tests of the piles and the ABAQUS for numerical simulations. The results of the simulated average lateral resistance of the soil layers are compared and analyzed with different friction coefficient values applied to the pile-soil interface. The results indicate that: (1) The load deformation characteristics of bored piles are closely related to the pile length and the properties of the bearing stratum. The axial force of super-long piles decreases linearly with depth, and the end resistance ratio significantly reduces. In cases where the bearing layer consists of moderately weathered granite, the axial force of the test pile decreases in a "bulging" shape along the depth direction, with the end resistance ratio lower than 17.1%. (2) The method for determining the friction coefficient in numerical simulations significantly impacts the simulated lateral resistance of the soil. Multiple friction coefficients at the pile-soil interface can improve the accuracy of the simulated results, depending on the nature of the soil layer. -
0. 引言
钙质砂是一种分布于热带海洋的特殊岩土介质,由于其成因和组构上的特殊性导致其工程力学性质与普通的陆源砂有明显区别。目前国内外针对钙质砂地基应力变化及破坏模式等研究很少,其工程特征及地基性能尚无明确统一的说法,这给基础设施的设计、施工建设带来很大的困难,因此需要研究钙质砂类岩土的应力特性和承载机理。
钙质砂的应力特性,国内外已作了一些研究。文哲等[1]采用中国南沙某岛礁吹填珊瑚钙质砂,在级配和组分分析的基础上,研究不同密实度、含水率条件下轴向荷载钙质砂的剪切力学特性。王刚等[2]进行钙质砂不破碎三轴剪切过程的数值试验,以确定临界状态与固定级配的对应关系。相盈盈[3]采用三轴固结不排水剪切试验,对南海海区岛礁饱和钙质砂进行剪切特性分析。张季如等[4]为了研究不同应力路径对钙质砂的颗粒破碎和力学性质的影响,对不同固结压力的钙质砂进行了5种应力路径下的排水三轴压缩试验。Lade等[5]针对珊瑚砂的时间变量开展了三轴压缩试验,提出了应力下降原因。Yamamuro等[6]开展了钙质砂三轴排水试验,考虑了不同的围压状态和相对密实度,抗剪强度与轴向应变的变化规律。其他学者[7-9]对钙质砂的工程承载特性进行了部分研究。
由以上研究可以发现,针对钙质珊瑚砂的力学工程特性性质已经已有一定基础,与其他砂类试验对比分析,结果显示钙质砂是一种天然特殊砂,钙质砂在高压条件下发生颗粒破碎,其破碎性质也会引起钙质砂的力学性质发生改变,其性质与普通砂存在很大差异,工程中不能使用普通砂的成果分析。目前有关三轴试验剪切速率的影响因素相关研究较少,本文在不同围压下采用不同剪切速率对高低围压下钙质砂进行三轴固结排水试验,分析围压、剪切速率对钙质砂应力影响分析。
1. 三轴排水剪切试验
1.1 试验材料
试验采用的试样取自中国海南省三沙市西沙群岛,试样中存在大块的砂砾和尖锐的贝壳,见图 1,为避免戳破橡胶膜影响试验,对原试样进行烘干、筛分、按一定比例进行配比后使用见图 2,试验筛除大于2 mm颗粒的钙质砂,将钙质砂分为9个粒径范围,不同粒径含量见表 1。颗粒级配见图 3,试验用钙质砂基本物理参数见表 2。
表 1 试样粒径范围Table 1. Grain-size ranges of sample粒径/mm 1.45~2 1~1.45 0.71~1 0.5~0.71 0.355~0.5 百分比/% 7 25 25 12 7 粒径/mm 0.25~0.355 0.18~0.25 0.125~0.18 0.075~0.125 百分比/% 5 4.5 4.5 10 表 2 试样基本物理参数Table 2. Basic physical parameters of sampleGS emax emin D60/mm D30/mm D10/mm 2.80 0.58 1.07 0.9 0.49 0.125 1.2 试验仪器
试验采用英国GDS公司生产的自动三轴试验系统GDSTAS完成(图 4)。应变控制速率从0.00001~9.99999 mm/min。压力室最大承受压力为3.5 MPa,2 MPa孔隙水压力传感器精度为全量程的0.15%,内置水下荷重传感器精度为全量程的0.1%。3 MPa/ 200 mm3标准围压/体积控制器(围压控制),具有自动压力和体积溢出保护,测量精度为量测值的0.25%,压力精度为全量程的0.1%,体积变化测量和显示分辨率为1 mm3。3 MPa/200 mm3标准反压/体积控制器(反压控制),具有自动压力和体积溢出保护,测量精度为量测值的0.25%,压力精度为全量程的0.1%,体积变化测量和显示分辨率为1 mm3。
1.3 试验方案
试验方案如表 3所示。试验为三轴固结排水试验,试样相对密实度为70%,有效围压为300,600 kPa。剪切应变速率分别为0.15%,0.035%,0.0045%/min,研究钙质砂在不同围压、不同剪切速率下的应力变化特性。试样直径为50 mm,高度100 mm。
表 3 试验方案Table 3. Test schemes相对密实度Dr 有效围压/kPa 剪切应变速率/(%·min-1) 0.70 300,600 0.0045,0.035,0.15 试验采用干法制样,将配比好的试样按一定质量分层装入模具,因钙质砂的特性导致内部空隙较多,试验中利用CO2、水头饱和、反压饱和相结合的方法对试样进行饱和,提高试样的饱和度。首先从试样底部缓慢通入CO2,从试样顶部排出,通气时间为2 h,此过程置换试样内部的空气。然后水头饱和,从试样底部缓慢进水,从试样顶部排水,控制水流速度,时间为2 h,此过程置换试样中的CO2。最后进行反压分级饱和,并监测试样饱和度,直到饱和度大于95%。接下来进行试样固结,实时监测孔隙水压力变化,直到孔隙水压力消散,有效围压达到预定值。按照设定的剪切速率进行三轴排水试验,轴向位移达到20%试验停止。
2. 结果分析
2.1 围压与剪切速率对偏应力影响
不同围压下各剪切速率的偏应力应变曲线,见图 5。偏应力曲线趋势相同,剪切初期偏应力上升较快,当应变超过5%时偏应力增长变缓,低围压下应变达到10%左右时偏应力出现峰值,高围压下应变达到12%左右时偏应力出现峰值,随后偏应力开始下降,高围压偏应力峰值出现的较晚。
不同围压下剪切速率对峰值的影响规律相同,偏应力峰值随剪切速率的增大而增大,剪切速率越快应力峰值越大,应力下降也越快应力下降比越大,可见剪切速率对应力的峰值大小及残余剪切强度有一定影响。低围压下不同剪切速率最大偏应力峰值差比为11%,且在低围压约束力下剪切速率减小时试样颗粒发生重组后再剪切,因此速率0.035%,0.0045%/min的偏应力峰值软化后又会出现硬化趋势;高围压下不同剪切速率最大偏应力峰值差比为4%,可见低围压下剪切速率对偏应力峰值影响更大。这与福建标准砂的性质不同,周杰[10]对不同应力下福建标准砂进行5种剪切速率试验,表明在低应力下砂的抗剪强度与剪切速率基本无关,高围压下,砂的抗剪强度随剪切速度增大而减小。
2.2 围压与剪切速率对偏应力比影响
不同围压下各剪切速率的偏应力比曲线,见图 6。偏应力比曲线趋势相同先增大后减小趋于稳定。围压越小偏应力比上升越快、偏应力比越大,可见低围压下的剪切效果更明显,偏应力和峰值的变化情况相对于平均有效应力的变化更大,因为围压越低对颗粒束缚较小,钙质砂的越容易发生错动导致偏应力上升。应力应变分析中低围压下剪切速率越小时偏应力出现峰值然后减小有又硬化的趋势,因此在低围压下剪切速率越小偏应力比峰值后减小有又上升趋势。
相同围压下剪切速率越大,偏应力增长越快、偏应力比越大,剪切时砂粒之间重组排列,高剪切速度使得砂粒之间迅速咬合偏应力快速上升导致试样破坏,钙质砂的结构特征使得对剪切速率更加敏感。但是相对于0.15%/min的剪切速度,0.035%/min和0.0045%/min的速度偏应力比的变化几乎相同,可见当剪切速率低于一定速率时,偏应力比的变化并不受剪切速度的影响。
2.3 围压与剪切速率对应力路径影响
固结排水试验的偏应力-平均有效应力曲线,见图 7。不同围压条件下,偏应力-平均有效应力关系曲线都是斜率为3的直线,是因为剪切时排水条件使得孔压保持稳定值不变,且围压为一定值,根据平均有效应力为p′=[2(σ3−μ)+(σ1−μ)]/3,σ3为围压,σ1为轴向应力,μ为孔隙水压力,经计算相邻点有效应力变化为p′=(σ11−σ12)/3,即为轴向应力差除以3,相邻点偏应力为q′=σ11−σ12,即为轴向应力差,所以曲线斜率为q′/p′=3。
不同剪切速率下的破坏线q′=Mp′趋势相同,由应力应变分析可以看出偏应力峰值随剪切速率的增大而增大,破坏线是关系偏应力峰值与圆点的连线,可见剪切速率越大破坏线角度越大,三条线的M值近似为1.75,由公式φ=arcsin(3M/(6+M))=42°,钙质砂的内摩擦角近似为42°。这由钙质砂的结构特性决定,钙质砂形状不规则且尖锐,颗粒之间咬合错动,剪切破坏角较大,内摩擦角较大。
2.4 围压与剪切速率对孔隙比影响
不同围压及剪切速率条件下,孔隙比-平均有效应力关系曲线,见图 8。偏应力峰值以前平均有效应力逐渐增大,偏应力峰值以后平均有效应力逐渐变小,且偏应力峰值以前不同围压下孔隙比相对于平均有效应力的下降速率相近呈线性变化,可见孔隙比与平均有效应力的变化趋势在剪切初始阶段不受剪切速率的变化影响。偏应力峰值以后孔隙比相对于平均有效应力变化呈曲线增大状态。
剪切过程中试样体积先减小后增大,孔隙比表现为先减小后增大,拐点发生在偏应力峰值的前后,剪切过程表现为密砂的特征。不同围压孔隙比下降速率相近,低围压剪缩阶段剪切速率对孔隙比的影响不明显,剪缩阶段颗粒之间位置调整砂粒发生相互错动重新排列压实,试样体积减小,孔隙比减小。剪胀阶段剪切速率对孔隙比的影响较大,剪切速率越大,孔隙比增大越快,随着剪切阶段颗粒间的作用力逐渐小于颗粒的抗剪强度试样破坏,尖利不规则的砂粒重新组合排列试样体积增大,孔隙比增大。
3. 结论
本文开展了钙质砂三轴排水剪切试验速率对应力特性影响研究,分析不同围压、不同剪切速率状态下钙质砂剪切过程中应力的变化规律,得出以下4点结论。
(1)偏应力应变曲线趋势相同,偏应力曲线均出现峰值,高围压偏应力峰值出现比低围压偏应力峰值出现的晚。剪切速率越快,偏应力峰值越大,偏应力下降越快,剪切速率对偏应力峰值大小及残余剪切强度有一定影响,低围压下剪切速率对偏应力峰值影响更大。
(2)偏应力比应变曲线趋势相同先增大后减小趋于稳定。围压越小偏应力比上升越快,偏应力比越大,相同围压下剪切速率越大,偏应力比增长越快偏应力比越大。当剪切速率低于一定速度时,偏应力比的变化并不受剪切速度的影响。
(3)不同围压下偏应力平均有效应力的关系是斜率为3的直线,不同速率下的破坏线趋势相同,偏应力峰值随剪切速率的增大而增大,破坏线是关系偏应力峰值与圆点的连线,因此剪切速率越大破坏线角度越大,破坏线的斜率近似为1.75,钙质砂内摩擦角近似为42°。
(4)偏应力峰值以前平均有效应力逐渐增大,偏应力峰值以后平均有效应力逐渐变小,且偏应力峰值以前不同围压下孔隙比相对于平均有效应力的下降速率相近。低围压剪缩阶段剪切速率对孔隙比下降速率的影响不明显,剪胀阶段剪切速率越大,孔隙比的增长越大。
-
表 1 灌注桩参数
Table 1 Parameters of cast-in-place piles
桩号 桩径/mm 桩长/m 混凝土强度等级 Z5 1000 67.8 C50 Z8 1200 51.6 C50 Z10 800 50.9 C50 表 2 灌注桩加载情况
Table 2 Loading conditions of cast-in-place piles
桩号 首级荷载/kN 分级荷载/kN 最终荷载/kN 加载级数 Z5 4000 2000 14000 6 Z8 6200 3100 32000 9 Z10 2000 1000 13000 12 表 3 桩周土实测极限平均侧阻力
Table 3 Measured average ultimate side resistances of soils
单位: kPa 土层 Z5 Z8 Z10 规范建议值[8] ①1杂填土 17.95 23.08 30.05 20~28 ②3淤泥 34.31 19.90 21.91 12~18 ④1砾质黏性土 41.52 132.98 37.38 38~96 ⑤1全风化花岗岩 55.56 150.61 105.41 120~140 ⑤2强风化花岗岩 60.80 150.62 191.24 160~240 ⑤4中风化花岗岩 73.42 218.11 — — 表 4 土层参数
Table 4 Soil parameters
土层 厚度/m 重度γ/(kN·m-3) 压缩模量Es/MPa 黏聚力c/kPa 内摩擦角φ/(°) 泊松比μ ①1 9.1 18.0 3.0 12.0 15 0.38 ②3 3.5 17.0 2.5 10.0 4 0.35 ④1 19.2 19.0 8.0 22.5 20 0.33 ⑤1 4.7 20.5 15.0 33.0 29 0.35 ⑤2 — 25.0 24.0 37.0 31 0.35 表 5 土层极限平均侧阻力
Table 5 Ultimate average lateral resistances of soil layers
土层 摩擦系数μ 模拟值/kPa 实测值/kPa 规范建议值[8]/kPa 比值 ①1 0.22 23.86 30.05 20~28 0.80 ②3 0.50 21.87 21.91 12~18 1.00 ④1 0.20 41.17 37.38 38~96 1.10 ⑤1 0.29 150.75 105.41 120~140 1.43 ⑤2 0.27 193.71 191.24 160~240 1.01 -
[1] 朴春德, 施斌, 魏广庆, 等. 分布式光纤传感技术在钻孔灌注桩检测中的应用[J]. 岩土工程学报, 2008, 30(7): 976-981. doi: 10.3321/j.issn:1000-4548.2008.07.005 PIAO Chunde, SHI Bin, WEI Guangqing, et al. Application of distributed fiber optic sensing techniques in bored pile detection[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(7): 976-981. (in Chinese) doi: 10.3321/j.issn:1000-4548.2008.07.005
[2] 杨建平, 朴春德, 常鸿飞, 等. 水平荷载下灌注桩变形分布式检测及承载机制研究[J]. 岩石力学与工程学报, 2014, 33(增刊1): 2983-2988. YANG Jianping, PIAO Chunde, CHANG Hongfei, et al. Distributed detection of deformation and bearing mechanism on bored pile under horizontal load[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S1): 2983-2988. (in Chinese)
[3] 贺嘉, 陈国兴. 基于ABAQUS软件的大直径桩承载力-变形分析[J]. 地下空间与工程学报, 2007, 3(2): 306-310. doi: 10.3969/j.issn.1673-0836.2007.02.025 HE Jia, CHEN Guoxing, Analysis of bearing capacity- deformation of the large diameter pile by the ABAQUS software[J]. Chinese Journal of Underground Space and Engineering, 2007, 3(2): 306-310. (in Chinese) doi: 10.3969/j.issn.1673-0836.2007.02.025
[4] 高鲁超, 叶生洋. 基于ABAQUS软件的单桩静载试验数值模拟[J]. 河北工程大学学报(自然科学版), 2015, 32(3): 51-54. GAO Luchao, YE Shengyang. Numerical simulation of the single pile static load test based on ABAQUS software[J]. Journal of Hebei University of Engineering (Natural Science Edition), 2015, 32(3): 51-54. (in Chinese)
[5] 秦仕伟, 高磊, 韦兵兵, 等. 基于光纤技术和有限元计算的灌注桩承载特性[J]. 科学技术与工程, 2023, 23(13): 5685-5692. QIN Shiwei, GAO Lei, WEI Bingbing, et al. Bearing characteristics of cast-in-place piles based on optical fiber technology and finite element calculation[J]. Science Technology and Engineering, 2023, 23(13): 5685-5692. (in Chinese)
[6] HUANG X D, WANG Y, SUN Y Y, et al. Research on horizontal displacement monitoring of deep soil based on a distributed optical fibre sensor[J]. Journal of Modern Optics, 2018, 65(2): 158-165. doi: 10.1080/09500340.2017.1382594
[7] 深圳市建筑基桩检测规程: SJG 09—2020[S]. 2020. Testing Specification of Building Foundation Piles in Shenzhen City: SJG 09—2020[S]. 2020. (in Chinese)
[8] 建筑基桩检测技术规范: JGJ 106—2014[S]. 2014. Technical Code for Testing of Building Foundation Piles: JGJ 106—2014[S]. 2014. (in Chinese)
[9] 曹文昭, 杨志银, 王雪阳. 大直径超长桩静载试验中试桩-锚桩相互影响分析[J]. 工业建筑, 2022, 52(9): 206-213. CAO Wenzhao, YANG Zhiyin, WANG Xueyang. Analysis of interaction between test piles and reaction piles of super-long large diameter piles in static loading tests[J]. Industrial Construction, 2022, 52(9): 206-213. (in Chinese)