基坑开挖下倾斜长短组合桩的受力变形特性
2022-05-14王际凯石钰锋刘伟煌詹刚毅蒋亚龙徐长节
王际凯 石钰锋 刘伟煌 詹刚毅 蒋亚龙 徐长节
摘 要:倾斜长短组合桩兼具变形控制能力良好和造价低的优点。目前,对倾斜长短组合桩的研究较少,尚缺乏对其工作机理的认识。利用模型试验结合数值模拟的方法,研究開挖过程中倾斜长短桩在不同倾角和桩长组合下的变形受力特性,继而探讨其工作机理。结果表明:对于倾斜长短组合桩而言,其支护效果与倾角成正比,与排桩的桩长成反比,且倾角的影响最大,桩长次之;在总桩长一定且短桩满足一定嵌入深度时,倾角存在一个临界值,当倾角小于临界值时,最佳桩长分配受倾角大小的影响,而当倾角大于临界值时,斜桩越长,支护效果越好。
关键词:基坑开挖;倾斜长短组合桩;土体与支护结构位移;桩身弯矩
中图分类号:TU473.1 文献标志码:A 文章编号:2096-6717(2022)04-0035-10
Mechanical and deformation characteristics of inclined long short composite pile under foundation pit excavation
WANG Jikai1, SHI Yufeng1, LIU Weihuang1, ZHAN Gangyi2, JIANG Yalong1, XU Changjie1,3
(1. Jiangxi Key Laboratory of Infrastructure Safety Control in Geotechnical Engineering, East China Jiaotong University, Nanchang 330013, P. R. China; 2. China Railway Shanghai Design Institute Group Co., Ltd., Shanghai 200070, P. R. China;3. Research Center of Coastal and Urban Geotechnical Engineering Zhejiang Unirersity, Hangzhou 310058, P. R. China)
Abstract: The inclined long short composite pile has the advantages of good deformation control ability and low cost. Currently, there are limited studies on this aspect and understanding of its working mechanism is still lacking. In view of this, through the model test and numerical method, the deformation and mechanical characteristics of long and short inclined piles under different inclination and pile length combinations during excavation are studied, and then the working mechanism is discussed. The results show that: For inclined long short composite piles, the supporting effect is directly proportional to the inclination angle and inversely proportional to the pile length of row piles, and the influence of inclination angle is the greatest, followed by the pile length; When the total pile length is fixed and the short pile meets a certain embedded depth, the inclination angle has a critical value. When the inclination angle is less than the critical value, the optimal pile length distribution is affected by the inclination angle. When the inclination angle is greater than the critical value, the longer the inclined pile, the better the supporting effect.
Keywords:foundation pit excavation; inclined long short composite pile; displacement of soil and supporting structure; pile bending moment
基坑工程的排桩围护结构一般采用直桩,而将直桩倾斜一定的角度,则为倾斜桩。在相同条件下,倾斜桩能够减小排桩的变形与内力,进而无需水平支撑[1-2]。近年来,陆续有学者对倾斜桩在基坑工程中的应用展开了研究。
在现场试验方面,孔德森等[3]通过模型试验的方法对倾斜悬臂桩的桩顶水平位移、桩后土体沉降和桩身弯矩进行了研究,从而探讨其受力特性。刘涛[4]将倾斜桩拓展使用至基坑支护工程,通过开展现场试验和有限元数值模拟分析,着重研究了倾斜支护桩在基坑中的作用机理。Han等[5]根据基础板的垂直应力和倾斜角度,探讨了倾斜基础板对埋入式桩端承载力和桩底土性的影响。Goit等[6]在1g条件下进行了模型物理缩尺试验,获得了埋置在干砂中的倾斜单桩的竖向响应。孔德森等[7]根据基坑支护倾斜悬臂桩模型试验,对比分析了倾斜角度、布桩方式、桩顶约束方式、桩位置4种因素对基坑支护倾斜悬臂桩受力变形特性的影响。
在倾斜桩竖向承载特性研究方面,郑刚等[8]提出了基坑斜直交替支护桩的3个工作机理效应,即刚架效应、斜撑效应和重力效应。孔德森等[9]为了研究动力荷载作用下倾斜桩的动力p-y曲线规律,采用数值模拟的方法构建了液化场地倾斜桩土数值分析模型。Misir等[10]采用室内模型试验和数值分析相结合的方法,研究了水平荷载作用下均质/层状砂土中垂直/倾斜放置情况下短桩的抗拔性能。刁钰等[11]将主动控制技术与斜直交替倾斜桩支护技术相结合,提出了一种新型的可以有效控制基坑变形的主动式无支撑支护技术。王耀鹏[12]通过模型试验研究了竖向荷载作用下无桩周土侧向位移和有桩周土侧向位移的3种倾斜率下桩的竖向承载变形性状。
在傾斜桩水平向承载特性研究方面,Wang等[13]建立了Timoshenko梁Pasternak地基(T-P)模型,以估算水平荷载作用下倾斜桩的水平动力阻抗和相互作用系数。周越洲等[14]结合珠海市某工程项目,利用有限元程序模拟了土方开挖过程中桩侧不平衡土压力对桩受力和变形的影响。Zhou等[15]借助数值模拟技术,研究了斜交PCP桩基础的变形特性,提出了斜向桩的加固方案。
上述对于倾斜桩的研究多以等长桩为研究对象,但由弹性地基梁法可知,基坑开挖面以下,围护结构受到被动土体的反向作用力,桩体下部承受较小弯矩,因此,可以通过减小部分桩长,形成长短桩,来实现降低造价的目的。基于此,笔者对倾斜长短组合桩的工作性状、支护能力进行了研究探讨。
1 倾斜长短桩试验研究
1.1 试验目的与内容
为探究基坑开挖下倾斜长短组合桩的支护效果,通过室内模型试验研究基坑开挖对冠梁侧向位移、桩身内力、土体沉降等方面的影响。已有研究表明,全斜桩支护能力弱于斜直交替桩型[11],故试验设计为一长一短组合(组合1~3)及对照组试验,具体见表1。
倾斜桩的倾斜方式以及采集装置的布设见图1。模型基坑分4次开挖,根据对称原则,取一半基坑进行模拟,故开挖范围为1.74 m×0.74 m,模型桩以空心PVC管模拟,其外径32 mm、壁厚2 mm、桩距为53.7 mm,冠梁规格为1.73 m×0.06 m×0.06 m(长×宽×高),分为冠梁顶板及冠梁底板,其大小尺寸相同。
1.2 试验概况
模型箱主体框架为防锈钢材,四面布置厚度10 mm的透明钢化玻璃以便于观察试验现象,并承受水平土压力以保证试验安全进行,如图2所示。
1.3 冠梁与模型桩
1.3.1 冠梁
根据试验需要,冠梁底板的预留桩孔设置0°、10°、20°的倾斜角度。底板预留桩孔直径为32 mm、孔位中心间距为53.7 mm。冠梁顶板预留桩孔起到限制桩顶穿过冠梁以及便于导线穿过的作用,其直径为15 mm。冠梁如图3所示。
1.3.2 模型桩
因为基坑开挖模拟试验中支挡结
构所承受的土压力较小,模型桩材料的弹性模量不宜太大,故参考文献[3],模型桩采用PVC(聚氯乙烯)管作为主体材料。桩体通过冠梁底板的预留桩孔来达到不同倾斜角度。此外,为了采集试验数据,在选定的模型桩上的特定位置(管材内壁)布设应变片。
1.4 试验用砂
采用试验室内储量颇丰的河砂作为地基填筑材料。为了去除砂土中的杂物,保证土样的均匀性,需要对河砂进行筛分,本试验采用了5 mm筛对河沙进行筛分,其主要物理力学参数见表2,颗粒级配曲线见图4。
1.5 试验过程
试验过程见图5,具体如下:
1)模型桩采用预先埋入法成桩,即先将桩固定在模型槽中一定深度后,在直桩迎土面布置止砂帷幕,再用沙雨法填筑预定高度。
2)对基坑内部的土体进行分次开挖,并记录相应的试验数据。
3)试验结束后静置24 h,观察试验现象。
2 试验结果分析
2.1 冠梁侧向位移
图6给出了倾角为0°、10°、20°情况下的一长一短组合桩与等长长桩的冠梁侧向位移。
根据冠梁侧向位移变化图可知:
1)随着基坑开挖深度的增加,冠梁侧向位移增大。当开挖至坑底时,除倾角为20°外,其余桩型的冠梁侧向位移的增长速率突然加快,被动区土体产生了塑性变形,土体的约束变小,对变形的约束较小。
2)倾角为0°的长短组合桩与等长长桩的冠梁侧向位移曲线极为接近,且二者的差值随着开挖深度的增加而增大。当开挖至坑底时,等长长桩的冠梁侧向为位移3.74 mm,而一长一短组合桩倾角为0°、10°、20°时,冠梁侧向位移分别为4.16、2.46、0.95 mm,相比等长长桩,其冠梁侧向位移增长率为11.2%、-34.2%、-74.6%。说明在整个支护体系中,斜桩起到斜撑的作用,进而影响整体支护作用,其中,倾斜的角度为传递此力的关键因素。
综上可知,在长短组合桩中,减小桩长会增大冠梁侧向位移,但变化的程度较小;增大长桩倾角会减小冠梁侧向位移,且变化的程度较大。
2.2 地表沉降
图7给出了当基坑开挖深度H=40 cm时,不同倾角下的一长一短组合桩与等长长桩的地表沉降结果。
根据地表沉降结果,可知:
1)开挖至坑底时,等长长桩的最大地表沉降为8.21 mm,而倾角为0°时的地表沉降最大值为9.12 mm,增长率为11.0%,与冠梁侧向位移增长率接近。由此可见,当长短组合桩倾角为0°时,减小其桩身长度都会增大冠梁侧向位移和最大地表沉降,且对二者的影响程度相当。
2)对比等长长桩可知,当倾角为10°、20°时,地表沉降最大值分别为6.22、2.67 mm,其地表沉降最大值增长率为-24.4%、-67.5%。由此可知,长桩倾角增大时,地表沉降减小幅度小于冠梁侧向位移,但相差较小,效果同样显著。这是因为基坑支护导致坑外土体沉降,基坑开挖导致坑底土体回弹,二者使支护桩两侧产生差值,导致桩顶水平位移增大,也因此使桩后土体产生沉降。
2.3 桩身弯矩
图8给出了开挖至坑底时0°、10°、20°倾角下的一长一短组合桩与等长长桩的桩身弯矩结果。图8中弯矩分析只考虑大小,不考虑正负。
根据桩身弯矩变化规律可知:
1)由刚架效应可知,冠梁将长短支护桩连接成一个整体,形成一种“三角形刚架支护体系”,其中,排桩中的斜桩承受较大的土压力,进而分担了直桩受到的作用力,使支护体系的抗倾覆性增强。
2)桩身上下弯矩符号相反,变化处为反弯点。在反弯点上部的桩身受压,下部的受拉。其中,斜桩的负弯矩值大于直桩的负弯矩值,这是因为在开挖过程中前排斜桩受到的被动土压力更大,类似于普通双排桩结构的前后排桩受力机理。
当倾角为0°时,对于开挖面以上的桩身弯矩最大值而言,等长长桩最小,最大弯矩为1.14 N·m,长短组合桩略大,其长桩与短桩的最大弯矩分别为1.42、1.39 N·m,长桩、短桩的弯矩增长率分别为24.6%、21.9%。这主要是因为桩长减小降低了嵌固长度。
4)当倾角为10°、20°时,斜桩在基坑开挖面以上的最大弯矩分别为1.29、1.22 N·m,基坑开挖面以下分别为1.20、0.94 N·m。与倾角为0°时的长桩相比,在开挖面以上斜桩桩身最大弯矩分别减小9.2%、14.1%;开挖面以下则分别减小15.5%、33.8%。故由此可知,长短组合桩中长桩的弯矩随倾角增大而减小,且在开挖面以下弯矩减小的趋势显著。
5)当倾角为10°、20°时,其短桩在基坑开挖面以上的最大弯矩分别为1.33、1.27 N·m,基坑开挖面以下分别为0.82、0.66 N·m。与倾角为0°时的短桩相比,在开挖面以上短桩的桩身最大弯矩分别减小4.3%、8.6%;开挖面以下分别减小21.9%、37.1%。故由此可知,长短组合桩中短桩的弯矩随倾角增大而减小,且在开挖面以下弯矩减小的趋势显著。
6)当倾角为10°、20°时,在基坑开挖面以上,斜桩的最大弯矩值稍小于直桩,而在基坑开挖以下,斜桩弯矩明显大于直桩弯矩。可以看出,开挖面以上,斜桩分担的弯矩略小于直桩;而开挖面以下,斜桩分担的弯矩更多。
3 数值模拟研究
借助有限元分析方法,按照相似比还原试验原型并拓展研究,在长桩、短桩的总桩长一定的条件下合理分配长桩、短桩的桩长。长桩、短桩总桩长为31.2 m,其中,CD1-1和CD1-2的桩长比值与试验对应,具体设计桩长分配方式如表3所示。
3.1 模型建立
由于基坑具有对称性,故取一半的基坑进行模拟,同时模型试验结果表明,在长桩倾角达到20°时,基坑的极限开挖深度有所增加。因此,模型在研究桩长分配影响规律的同时,基坑将分7次开挖,以此来探究极限开挖深度。为减小边界效应,模型横向、竖向分别取71.8、50 m;为对桩变形受力特性进行精细化模拟,桩和冠梁采用梁单元,土体采用实体单元,网格采用混合网格,且对桩与基坑周围局部加密,模型共计23 175个单元,如图9所示。
有限元软件Midas-GTS可以通过设置施工工况对基坑施工的全过程进行模拟,数值分析的工况如表4所示。
3.2 材料参数的确定
以南昌地区典型的上软下硬地層为基础,选用经典模型本构修正MC本构,土层参数如表5所示,由于在地勘报告中,土体弹性模量E并未给出,故根据经验公式E0=(2~5)Es进行调整。在数值分析中不考虑地下水影响。支护桩、冠梁均采用C30混凝土。
3.3 有限元模拟结果分析
3.3.1 桩身深层水平位移分析
将不同桩长分配下的桩身深层水平位移与开挖深度的变化关系绘制成图10。
由图10可知:
1)在同一开挖深度下,倾角越大,桩身深层水平位移最大值越小,说明倾角增大可以有效控制桩身深层水平位移。此结论符合模型试验研究中倾角对冠梁侧向位移和地表沉降的影响趋势,故数值结果与模型试验可相互验证。
2)倾角为10°时,桩身深层水平位移曲线变化规律为从桩顶到桩底,位移逐渐减小,这是因为斜桩桩底的支点力在水平方向的分量并不大,对限制桩顶水平位移的作用相对较小。
3)倾角为20°时,桩身深层水平位移曲线变化趋势呈现“中间大,两头小”的规律,这是因为倾角较大时,冠梁、直桩、斜桩组成的空间超静定结构更加稳定,斜桩可以看成斜撑,将桩顶牢牢固定住。
4)在模型试验和数值分析中,开挖至坑底时,倾角为20°时的冠梁侧向位移与地表沉降值均未产生大幅度增长,说明被动区土体并未发生塑性变形,故得出倾角增大可以提高基坑极限开挖深度的结论。
5)倾角不同,桩身深层水平位移曲线趋势也并不相同,CD1-2组的桩顶侧向位移在开挖深度为9.6 m时为负值,即向坑外移动,且斜桩下部侧向位移皆大于直桩。这是因为倾斜越大,斜桩端部支点力水平分量越大,且斜桩下部越靠近基坑中部位置,受到坑内土体卸荷回弹的作用更大。
为了便于对比不同倾角、不同桩长分配方式的倾斜长短组合桩桩身深层水平位移曲线变化规律,以开挖深度10.8 m为例,将不同桩长分配下斜桩与直桩的桩身深层水平位移随倾角的变化曲线分别绘制于图11。
由图11可知:
1)在开挖深度相同的条件下,当倾角为10°时,不同桩长分配下斜桩、直桩的桩身位移从大到小依次为SCD3-1>SCD1-1>SCD2-1;当倾角为20°时,不同桩长分配下斜桩、直桩的桩身位移从大到小依次为SCD3-2>SCD2-2>SCD1-2。
2)当长短组合桩总桩长一定时,倾斜长短组合桩在控制支护结构位移方面是要明显优于等长倾斜桩的。
3)当长短组合桩总桩长一定时,桩长分配对支护结构位移的影响与倾角有关。
3.3.2 土体与支护结构位移分析
各组在不同开挖深度下的支护结构最大侧向位移变化情况如表6所示。
由表6可知:
1)同等条件下,倾角在0°~20°范围内,2种倾斜长短组合桩支护效果各有强弱,当倾角为10°时,在斜桩桩长增加的过程中,支护效果先增强后减弱,而当倾角为20°时,支护效果随斜桩桩长增加而增强。由以上可知,在总桩长一定且短桩满足一定嵌入深度时,倾角存在一个临界值,当倾角小于临界值时,最佳桩长分配受倾角大小的影响,而当倾角大于临界值时,斜桩越长,支护效果越好。
2)对于倾斜长短组合桩而言,极限开挖深度随着倾角增大而明显增加,而等长倾斜桩的极限开挖深度受桩身嵌入深度的影响更大,倾角的增加更可能导致斜桩的嵌入深度不足。
3.3.3 桩身弯矩分析
图12为不同开挖深度时的桩身弯矩。弯矩分析只考虑大小,不考虑正负。表7~表10分别为以上6组在不同开挖深度下的基坑开挖面上下斜桩与直桩的最大弯矩值。
根据桩身弯矩变化规律并由表7~表10可知:
1)在开挖面以上,随着开挖深度的增加,斜桩、直桩的最大弯矩值增加;在开挖面以下,随着开挖深度的增加,斜桩、直桩最大弯矩值先增加后减小,直至面向基坑开挖的一侧桩身由受压转为全部受拉。
2)在倾斜长短组合桩中,随着倾角的增大,斜桩、直桩的最大弯矩值呈减小的趋势;且开挖面以下变化显著。该结论与模型试验所得弯矩结论可相互验证。
3)在等长倾斜桩中,随着倾角的增大,斜桩、直桩的最大弯矩值呈减小的趋势;且开挖面以下变化显著。
图13为不同桩长分配时的桩身弯矩。
由图13可知:
1)在总桩长一定的条件下,在开挖面以上,斜桩、直桩的桩身最大弯矩值基本相近;在开挖面以下,当倾角为10°时,斜桩的桩身最大弯矩值从小到大依次为MCD3-1 2)当倾角为20°时,斜桩的桩身最大弯矩值从小到大依次為MCD3-2 4 结论 基于室内模型试验探究基坑开挖下倾斜长短组合桩的支护效果,并通过数值手段研究总桩长一定时,长桩与短桩的桩长比值变化对支护效果的影响。得到如下主要结论: 1)长短组合桩的支护效果与桩的倾角和桩长分配有关。 2)倾角增大,基坑抗倾覆性增强,极限开挖深度提高;在基坑开挖面以上,斜桩分担弯矩总是小于直桩,而在开挖面以下,则反之。 3)在长短组合桩中,随着倾角的增大,斜桩、直桩的最大弯矩值和最大水平位移呈减小的趋势,且主要体现在基坑开挖面以下。 4)由模型试验可知,减小部分桩长会造成长短组合桩支护效果的减弱,增大倾角会增强长短组合桩支护效果,且后者的影响大于前者。 5)在总桩长一定且短桩满足一定嵌入深度时,倾角存在一个临界值,当倾角小于临界值时,最佳桩长分配受倾角大小的影响,而当倾角大于临界值时,斜桩越长,支护效果越好。 参考文献: [1] 徐江, 龚维明, 张琦, 等. 大口径钢管斜桩竖向承载特性数值模拟与现场试验研究[J]. 岩土力学, 2017, 38(8): 2434-2440, 2447. XU J, GONG W M, ZHANG Q, et al. Numerical simulation and field test study on vertical bearing behavior of large diameter steel of inclined piles [J]. Rock and Soil Mechanics, 2017, 38(8): 2434-2440, 2447. (in Chinese) [2] CHEN Y D, DENG A, WANG A T, et al. Performance of screw-shaft pile in sand: Model test and DEM simulation [J]. Computers and Geotechnics, 2018, 104: 118-130. [3] 孔德森, 张杰, 王士权, 等. 基坑支护倾斜悬臂桩受力变形特性试验研究[J]. 地下空间与工程学报, 2020, 16(1): 160-168. KONG D S, ZHANG J, WANG S Q, et al. Experimental study on stress and deformation characteristics of cantilever inclined retaining pile for foundation pit support [J]. Chinese Journal of Underground Space and Engineering, 2020, 16(1): 160-168. (in Chinese) [4] 刘涛. 倾斜桩在基坑支护中的作用机理和简化计算方法[D]. 天津: 天津大学, 2018. LIU T. The function mechanism and simplified calculation method of inclined pile in foundation pit support [D]. Tianjin: Tianjin University, 2018. (in Chinese) [5] HAN K, SEO M J, HONG W T, et al. End-bearing capacity of embedded piles with inclined-base plate: Laboratory model tests [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2020, 146(8): 04020063. [6] GOIT C S, SAITOH M, IGARASHI T, et al. Inclined single piles under vertical loadings in cohesionless soil [J]. Acta Geotechnica, 2021, 16(4): 1231-1245. [7] 孔德森, 张杰, 王士权, 等. 倾斜悬臂支护桩受力变形特性模型试验[J]. 工业建筑, 2019, 49(3): 117-121,70. KONG D S, ZHANG J, WANG S Q, et al. Model test on bearing and deformation characteristics of inclined cantilever pile in deep excavation support [J]. Industrial Construction, 2019, 49(3): 117-121,70. (in Chinese) [8] 郑刚, 何晓佩, 周海祚, 等. 基坑斜直交替支护桩工作机理分析[J]. 岩土工程学报, 2019, 41(Sup1): 97-100. ZHENG G, HE X P, ZHOU H Z, et al. Working mechanism of inclined-vertical retaining piles in excavations [J]. Chinese Journal of Geotechnical Engineering, 2019, 41(Sup1): 97-100. (in Chinese) [9] 孔德森, 李文胜, 常龙龙. 液化场地倾斜桩动力p-y曲线研究[J]. 地震工程与工程振动, 2019, 39(4): 41-56. KONG D S, LI W S, CHANG L L. Study on the dynamic p-y curves of inclined pile in liquefied ground [J]. Earthquake Engineering and Engineering Dynamics, 2019, 39(4): 41-56. (in Chinese) [10] MISIR G, LAMAN M. Numerical analysis of laterally loaded battered rigid piles [J]. Teknik Dergi, 2019,30(5):9507-9531 [11] 刁钰, 苏奕铭, 郑刚. 主动式斜直交替倾斜桩支护基坑数值研究[J]. 岩土工程学报, 2019, 41(Sup1): 161-164. DIAO Y, SU Y M, ZHENG G. Numerical study on oblique-straight alternate piles for active support of excavations [J]. Chinese Journal of Geotechnical Engineering, 2019, 41(Sup1): 161-164. (in Chinese) [12] 王耀鹏. 有桩周土侧向位移的倾斜桩竖向承载性能研究[D]. 济南: 山东建筑大学, 2017. WANG Y P. Study on lateral displacement of the pile with soil pile vertical bearing capacity of the deflection of accident and reinforcement research [D]. Jinan: Shandong Jianzhu University, 2017. (in Chinese) [13] WANG J, ZHOU D, JI T J, et al. Horizontal dynamic stiffness and interaction factors of inclined piles [J]. International Journal of Geomechanics, 2017, 17(9): 04017075. [14] 周越洲, 李重陽, 李国清, 等.软土地基倾斜桩竖向承载力研究[J/OL]. 建筑结构学报.https://doi.org/10.14006/j.jzjgxb.2020.0349. ZHOU Y Z, LI C Y,LI G Q, et al. Research on vertical bearing capacity of inclined pile in soft soil foundation [J/OL]. Journal of Building Structures, https://doi.org/10.14006/j.jzjgxb.2020.0349. (in Chinese) [15] ZHOU D Q, FENG C X. Engineering characteristics and reinforcement program of inclined pre-stressed concrete pipe piles [J]. KSCE Journal of Civil Engineering, 2019, 23(9): 3907-3923. (编辑 胡玲) 收稿日期:2021-03-02 基金项目:江西省自然科学基金(20192BAB206043);国家地区科学基金(51768020) 作者简介:王际凯(1996- ),男,主要从事基坑工程研究,E-mail:512797846@qq.com。 石钰锋(通信作者),男,博士(后),E-mail:270347713@qq.com。 Received:2021-03-02 Foundation items:Jiangxi Natural Science Foundation (No. 20192BAB206043); National and Regional Science Foundation Project (No. 51768020) Author brief:WANG Jikai (1996- ), main reseatch interest: foundation pit engineering, E-mail: 512797846@qq.com. SHI Yufeng (corresponding author),postdoctor, E-mail: 270347713@qq.com.