Loading [MathJax]/jax/output/SVG/jax.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

基于弱有限元法的渗流分析

何朝葵, 速宝玉, 盛金昌, 罗玉龙

何朝葵, 速宝玉, 盛金昌, 罗玉龙. 基于弱有限元法的渗流分析[J]. 岩土工程学报, 2023, 45(7): 1526-1532. DOI: 10.11779/CJGE20220635
引用本文: 何朝葵, 速宝玉, 盛金昌, 罗玉龙. 基于弱有限元法的渗流分析[J]. 岩土工程学报, 2023, 45(7): 1526-1532. DOI: 10.11779/CJGE20220635
HE Zhaokui, SU Baoyu, SHENG Jinchang, LUO Yulong. Seepage analysis based on weak finite element method[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(7): 1526-1532. DOI: 10.11779/CJGE20220635
Citation: HE Zhaokui, SU Baoyu, SHENG Jinchang, LUO Yulong. Seepage analysis based on weak finite element method[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(7): 1526-1532. DOI: 10.11779/CJGE20220635

基于弱有限元法的渗流分析  English Version

详细信息
    作者简介:

    何朝葵(1974—),男,硕士,讲师,主要从事有限元方法方面的教学和科研工作。E-mail: hzk@hhu.edu.cn

  • 中图分类号: TU139.14

Seepage analysis based on weak finite element method

  • 摘要: 在坝、闸和堤防等水工建筑物的设计和运行管理中,渗流计算占有非常重要的地位。目前渗流分析的方法以一般有限元方法为主。在渗流分析时往往希望通过加密网格来追求更高精度的数值结果,但是由于一般有限元法属于协调有限元方法,它不适合处理混合剖分网格,加密网格容易导致计算格式不稳定、计算结果不收敛。而弱有限元方法是一种非协调有限元方法,用弱函数的弱梯度算子替代一般有限元方法变分式中的经典梯度算子,并在变分式中增加稳定子项,保证了计算格式的绝对稳定性,并且它适用于混合网格。对于同一剖分网格,弱有限元法的总体代数方程组的自由度远远比一般有限元方法的要大很多,为了降低弱有限元法的自由度,探讨根据弱有限元方法适用于混合剖分网格的特点,采用有针对性的混合剖分网格,运用弱有限元法求解渗流自由面和进行有防渗帷幕的闸基渗流分析,数值结果表明弱有限元方法能灵活处理混合剖分网格,并且计算结果具有高精度。
    Abstract: In the design and operation management of hydraulic structures such as dams, gates and embankments, the seepage calculation plays a very important role. The current methods for seepage analysis are mainly based on the general finite element method. It is often hoped, by refining the mesh, to pursue the calculated results with higher-precision in the seepage analysis. However, the general finite element method, which belongs to the coordinated finite element method, cannot handle hybrid grids, and the refined mesh may lead to unstable calculation formats and the divergence of the calculated results. The weak finite element method is a non-coordinated one, which replaces the classical gradient operator in the variational formula for the general finite element method with the weak gradient operator of the weak function, and adds a stable sub-term in the variational formula, to obtain an absolute stability of the calculation format, and it can handle hybrid grids. For the same mesh grid, the degree of freedom of the overall algebraic equations for the weak finite element method is much larger than that of the general finite element method. In order to reduce the degree of freedom of the weak finite element method, the targeted mixed meshing grid is used to establish the weak finite element method which is used to solve the free surface of seepage and to analyze the seepage field of the gate foundation with anti-seepage curtain. The numerical simulation shows that the weak finite element method can be used to process flexibly hybrid grids, and it is of high accuracy.
  • 近年来,盾构法在城市地下工程修建中得到了广泛应用,而随着城市地下空间的开发及交通需求的增长,盾构隧道的建设趋于大断面化和大埋深化,所穿越地层的地质条件愈加复杂,施工中常面临地层软弱、高外水压等一系列问题。传统的单层管片衬砌结构难以满足盾构隧道安全性和耐久性的要求,为此国内武汉地铁8号线越江隧道、狮子洋水下隧道等盾构隧道工程开始采用双层衬砌作为支护结构。相比于单层管片衬砌,双层衬砌结构在补强、防水和防侵蚀等方面具有更多优势,在盾构隧道建设中的应用越来越广泛[1]

    目前,国内外已有许多学者采用理论分析、数值模拟、模型试验以及现场监测的方法对盾构隧道双层衬砌结构进行了研究。早期,村上博智等[2]开展了双层衬砌模型试验并提出了复合、叠合两种衬砌结构形式及特点。张厚美等[3]依托穿黄隧洞工程提出了双层衬砌叠加受力模型,得到了双层衬砌联合承载的内力分布曲线。王士民等[4]建立了基于塑性损伤的双层衬砌三维实体非连续接触模型,真实模拟了盾构隧道管片衬砌与二次衬砌接合面的接触情况。刘洋等[5]、高鑫等[6]通过数值模拟的方式,对双层衬砌结构横向力学特性、非协调变形、应力集中等力学行为进行了分析。

    封坤[7]、于清洋[8]完善了盾构隧道模型试验加载装置,实现了不均匀水土荷载的模拟施加,为双层衬砌相似模型试验打下了基础。何川等[9]依托狮子洋水下盾构隧道工程,开展双层衬砌纵向力学性能模型试验,研究了在软硬交替地层下双层衬砌隧道纵向内力与变形的变化规律。王士民等[10]开展了二次衬砌合理施作时机模型试验,对管片衬砌与二次衬砌的主从关系进行研究,认为二次衬砌在盾构隧道双层衬砌中仅作为辅助承载结构,应在管片衬砌承受部分设计荷载之后施作。郭文琦等[11]通过设计模型试验,探究了二次衬砌厚度对双层衬砌纵向力学性能的影响,综合考虑得到了大断面水下盾构隧道二次衬砌合理厚度的合理取值范围。

    由于上述双层衬砌结构的研究重点多集中在内外衬砌结构体系分析方面,实际工程中隧道洞身所穿越的地层条件复杂多变,而针对不同地层条件下双层衬砌结构力学特性变化的研究国内还鲜有报道。郭瑞等[12]认为地层侧压力系数的变化对隧道整体稳定有较大影响,且提出了一种基于极限位移与隧道直径变化率的衬砌结构失稳判据。但是目前,部分考虑侧向土压力对盾构隧道影响的研究[13-15]主要集中在管片衬砌结构,对于二次衬砌施作之后到破坏阶段的双层衬砌整体力学特性无详细讨论。以侧压力系数为指标,为进一步探明侧压力系数对盾构隧道双层衬砌结构力学特性的影响,确保结构满足安全运行要求。本文依托广深客运专线狮子洋盾构隧道工程,通过相似模型试验,针对不同侧压力系数下盾构隧道双层衬砌结构的力学特性进行了相关研究,揭示了侧压力系数对双层衬砌结构力学特性及破坏形态的影响规律,同时也对管片衬砌单独承载与双层衬砌联合承载性能差异进行了分析讨论。

    试验依托狮子洋水下铁路大直径盾构隧道工程,盾构段隧道长10.8 km,在处于软弱地层及软硬不均地层的进出洞口段采用双层衬砌结构。如图 1所示,其中管片衬砌采用“5+2+1”分块形式,外直径10.8 m,内直径9.8 m,管片厚度50 cm,幅宽2 m,混凝土强度等级为C50;二次衬砌环外直径为9.8 m,内直径9.2 m,厚度为30 cm,采用C30素混凝土。

    图  1  狮子洋盾构隧道双层衬砌分块示意
    Figure  1.  Schematic diagram of double-layer lining of Shiziyang tunnel

    相似试验按照“地层-结构模式”设计模拟。如图 2所示,试验采用“盾构隧道-地层复合体模拟试验系统”作为加载装置,中央装载隧道模型并装填模型土[16]。在平行于隧道横断面方向,该装置通过控制地层内部应力场调整地层侧压力系数,其中Ⅰ方向布置4台千斤顶施加竖向土压力,Ⅱ方向布置4台千斤顶施加侧向土压力;垂直隧道横断面方向通过控制千斤顶施加纵向约束力保证双层衬砌在试验加载过程中保持平面应变状态,其量值大于Ⅰ,Ⅱ方向荷载值。

    图  2  模型试验加载装置示意
    Figure  2.  Schematic diagram of loading devices for model tests

    参考模型试验相似理论及方法,根据相似理论,在本研究中采用几何相似比CL=20作为整个相似模型试验的基础相似比,表 1为其余各物理参数的原型与模型值之比[17]

    表  1  相似关系表
    Table  1.  Ratios of similarities
    物理量 符号 相似比 单位
    强度 R CR=20 Pa
    弹性模量 E CE=20 Pa
    应力 σ Cσ=20 Pa
    黏聚力 c Cc=20 Pa
    内摩擦角 φ Cφ=1 (°)
    下载: 导出CSV 
    | 显示表格

    根据相似关系,按照一定的比例将河砂、粉煤灰、石英砂、机油、重晶石粉及松香均匀混合模拟实际工程中狮子洋隧道洞身穿越段地层的土体[18],模型土与原型土的物理力学参数如表 2所示。

    表  2  土体材料物理力学参数表
    Table  2.  Physical and mechanical parameters of soil materials
    名称 γ/(kN·m-3) E/MPa c/kPa φ/(°)
    原型值 18.7~20.3 15.0~25.0 0 20.0~32.0
    模型值 20.0 1.0 0 28.0
    下载: 导出CSV 
    | 显示表格

    管片衬砌与二次衬砌均由可塑性较强的雕塑石膏按照不同的水灰比调制而成,双层衬砌模型尺寸按照相似比关系确定,如图 3所示。同时为保证管片衬砌的等效抗弯刚度,选择在管片衬砌中内嵌钢丝网模拟衬砌的环向主筋。

    图  3  模型结构及加载示意
    Figure  3.  Schematic diagram of model structure and loading

    图 4所示,根据盾构隧道管片衬砌不同位置的弯矩典型分布规律和环向接头的拉压特性,采用内外分区割槽的方式模拟管片衬砌环向接缝,弱化管片衬砌在环向接头处的等效刚度,从而更好地模拟实际工程中接头位置的强度削弱[19],其中正弯矩区深1.5 cm,负弯矩区深1.6 cm。此外,为约束管片衬砌环间错动,精确表征双层衬砌结构横向力学特性,通过细钢棒模拟管片衬砌纵向接头螺栓。

    图  4  盾构隧道管片衬砌接头处理示意
    Figure  4.  Treatment of segmental sheet lining joint of shield tunnel

    在模型制作过程中,叠合结构由于在管片衬砌内表面接头位置预制凹槽,使得二次衬砌与管片衬砌在该位置互相咬合,进而改变结构的弯曲刚度,复合结构则无此影响。此外,试验中为了减少切向剪力的干扰,使结构仅传递管片衬砌与二次衬砌间的法向压力,在模型管片衬砌施作完成之后,选择复合的方式施作二次衬砌[20]。如图 5所示,具体为管片衬砌结构加载至设计值时,在管片衬砌与二次衬砌之间放置一层塑料薄膜,用以模拟工程实际中的防水层,然后通过现浇的方式施作预制好的二次衬砌,实现二者的复合。

    图  5  接合面为复合类型的双层衬砌结构示意
    Figure  5.  Schematic diagram of double-lining structure with composite type at joint surface

    为了获得不同侧向土压下双层衬砌结构的各项试验数据,以中间整环管片及对应的二次衬砌作为主要量测对象,在二次衬砌浇筑前,将电阻应变片每隔15°粘贴在拼装好的管片环及预制石膏二次衬砌的内外侧,合计4环共96个测点,利用测得的应变数据计算管片衬砌与二次衬砌的弯矩和轴力[21]。除此之外,在内侧相同方位布设6个电阻式位移计,记录双层衬砌横向位移;在双层衬砌内侧拱顶、拱底及左右拱腰处安装声发射探头记录声发射信息。3种不同测量元件的排布方式如图 6所示。

    图  6  测试元件排布方式示意
    Figure  6.  Schematic diagram of arrangement of test components

    试验采用分级加载方式,在垂直于管片横断面方向先控制竖向千斤顶将模型纵向约束力设为18 MPa,然后同时施加竖向(Ⅰ方向)与横向(Ⅱ方向)土压力,最后通过环箍钢丝绳施加管片衬砌外侧的水压力。竖向土压加至工程实际埋深的设计荷载值(第8荷载步)后,在稳载状态下施作二次衬砌并继续加载直至双层衬砌结构整体损伤破坏(第24荷载步),加载工况如表 3所示。其中,每一个荷载步Ⅰ方向土压力按照千斤顶油压刻度0.4 MPa逐级加载,Ⅱ方向土压力值通过相应的侧压力系数计算得出。

    表  3  试验加载参数表
    Table  3.  Loading parameters for tests
    加载步 千斤顶油压/MPa 模型拱顶地层压力/kPa 模拟等效覆土厚度/m
    0 0.0 0.00 0
    1 0.6 1.44 3
    2 1.0 5.52 11
    3 1.4 8.67 17
    4 1.8 11.86 25
    5 2.2 14.46 30
    6 2.6 18.70 37
    7 3.0 21.79 43
    8 3.4 24.95 50
    9 3.8 29.29 60
    10 4.2 33.96 70
    11 4.6 41.88 85
    12 5.0 47.31 95
    13 5.4 52.83 105
    14 5.8 57.34 115
    15 6.2 61.33 125
    16 6.6 65.37 135
    17 7.0 70.04 145
    18 7.4 75.96 155
    19 7.8 81.39 165
    20 8.2 86.91 175
    21 8.6 91.42 185
    22 9.0 95.41 195
    23 9.4 99.45 205
    24 9.8 103.58 215
    下载: 导出CSV 
    | 显示表格

    大型跨江、海盾构隧道通常会穿越多种地层,地质条件复杂多变,地层性质呈现很大的不均匀性,广深客运专线狮子洋盾构隧道洞身主要穿越地层有淤泥质土层、砂层、岩石强风化带、黏性土层,根据相关地勘资料,本次试验主要探究侧压力系数分别为0.5,0.6,0.7,0.8的情况下,双层衬砌结构的力学特性变化规律以及对高侧压力地层的适应性﹐试验工况见表 4

    表  4  试验工况表
    Table  4.  Test conditions
    试验组号 侧压力系数 二衬厚度/cm 水头高度/m
    1 0.5 30 30
    2 0.6
    3 0.7
    4 0.8
    下载: 导出CSV 
    | 显示表格

    (1)管片衬砌受力特性分析

    图 7为不同侧压力系数条件下管片衬砌轴力弯矩变化曲线,其中红色虚线为管片衬砌弹塑性分界点(结构内力与声发射事件率突变),粉色虚线为管片衬砌失稳破坏临界点(结构轴力突变)。试验加载过程中,二次衬砌在管片衬砌承受100%设计荷载值即第8荷载步时施作[22]

    图  7  不同侧压力系数条件下管片衬砌内力随加载变化曲线
    Figure  7.  Variation curves of internal force of segmental lining with loading under different lateral pressure coefficients

    根据图 7可知,随着侧压力系数的增大,管片衬砌的弹塑性分界点和失稳破坏临界点逐渐后移。当侧压力系数为0.5时,管片衬砌内力呈线性变化,在前6级荷载下处于弹性变形阶段,二次衬砌施作前管片衬砌已进入塑性变形阶段,当荷载等级达到14级时,拱底的轴力急剧增加,表现出失稳破坏;而当侧压力系数增至0.8时,管片衬砌的弹性变形阶段延伸到12级荷载,并且拱顶的弯矩在第20级荷载时才表现出突增,发生失稳破坏,其中管片衬砌的失稳荷载在侧压力系数由0.7增至0.8时变化最大,为4个荷载步。

    在试验加载过程中,盾构隧道管片衬砌所承受的弯矩在拱底和拱顶表现为内侧受拉,而在左右拱腰则为外侧受拉。轴力以全环受压的形式存在。随着侧压力系数的增大,围岩土体对管片衬砌变形的约束作用增强,导致管片衬砌的弹塑性临界荷载及失稳破坏的临界荷载增大。管片衬砌的失稳位置主要集中在拱底处,拱腰在低侧压力系数条件下更易发生失稳破坏。

    (2)二次衬砌受力特性分析

    为纵向探究不同侧压力系数下二次衬砌结构力学特性及横向对比二次衬砌与管片衬砌力学特性差异。提取加载试验的应变数据,计算相应内力,绘制出二次衬砌弯矩和轴力随加载步变化曲线,如图 8所示,虚线为二次衬砌失稳破坏临界点。

    图  8  不同侧压力系数条件下二次衬砌内力随加载变化图
    Figure  8.  Variation of internal force of secondary lining with loading under different lateral pressure coefficients

    图 8中可以看出,二次衬砌在施作完成后即开始承受一定的弯矩和轴力,在其施作过程中石膏的水化作用导致二次衬砌内部出现初始应力。当侧压力系数为0.5时,根据弯矩和轴力变化,加载至第14级荷载,二次衬砌的右拱腰和拱顶位置轴力突变,发生失稳破坏;当侧压力系数增至0.8时,加载至22级荷载,二次衬砌左、右拱腰和拱底位置发生失稳破坏,侧压力系数从0.6增至0.7时,二次衬砌失稳荷载变化最大,为7个荷载步。

    二次衬砌与管片衬砌所受的弯矩和轴力表现相似,但所承担的荷载水平更低。随着侧压力系数的增大,二次衬砌失稳破坏临界荷载值也增加,较高的侧压力系数使得二次衬砌失稳破坏位置主要集中在拱顶和拱底。相对于管片衬砌而言,二次衬砌的失稳破坏发生时间较晚。因此,二次衬砌能够有效协助管片衬砌承担外部荷载。增大侧压力系数有助于提高二次衬砌的承载能力,确保盾构隧道双层衬砌结构的安全性。

    图 9为不同侧压力系数对应的双层衬砌位移变化曲线(位移数据以径向向外为正)。从图 9中可以看出,在外部荷载的作用下,隧道拱顶发生沉降变形,拱底隆起,左右拱腰变形径向向外,衬砌整体椭变呈“横鸭蛋”状,最大位移出现在双层衬砌拱底处,最大位移随荷载步增加而不断增大。二次衬砌施作后,双层衬砌变形速率减缓,衬砌整体刚度加强。随侧压力系数的增大,双层衬砌变形由拱腰外凸向拱顶和拱底外凸趋近,最后变形趋于均匀,结构椭圆度下降。为进一步分析侧压力系数变化带来的影响,提取了衬砌二次失稳时对应的位移数据进行分析,如表 5所示。

    图  9  不同侧压力系数条件下衬砌变形随加载变化图
    Figure  9.  Variation of lining deformation with loading under different lateral pressure coefficients
    表  5  双层衬砌失稳变形统计
    Table  5.  Instability deformation statistics of double-layer lining
    侧压力系数 失稳荷载等级 失稳位置 失稳变形/mm
    0.5 14 拱底 15.17
    0.6 14 拱底 10.68
    0.7 21 拱底 10.59
    0.8 22 拱底 9.45
    下载: 导出CSV 
    | 显示表格

    随着侧压力系数的变化,双层衬砌失稳破坏时的变形结果有所差异。由表 5中数据可得,侧压力系数从0.5增大至0.8,双层衬砌失稳荷载等级增大,双层衬砌承载能力增强。双层衬砌失稳临界变形值由15.17 mm降至9.45 mm减小了37.71%,其中最大降幅为29.60%,在侧压力系数为0.5~0.6时取得,可以认为地层侧压力的增大限制了双层衬砌的横向变形。

    为探明双层衬砌中二次衬砌与管片衬砌的相互作用关系,进一步分析地层侧压力系数对双层衬砌结构受力及变形特性的影响,需对双层衬砌内力与变形结果进行联合分析。基于此,将试验中二次衬砌失稳破坏时双层衬砌的关键数据汇总于表 6

    表  6  双层衬砌结构试验结果汇总
    Table  6.  Summary of test results for double-layer structures
    二次衬砌失效时衬砌结构试验数据 侧压力系数
    0.5 0.6 0.7 0.8
    竖向外荷载/kPa 57.34 57.34 91.42 95.41
    最大径向位移/mm 15.17 10.68 10.59 9.45
    最大弯矩值/(kN·m) 954.26 885.52 845.73 750.54
    最大轴力值/kN 4322.43 5699.15 8030.46 9421.30
    最大偏心距/mm 220.77 155.54 105.32 79.66
    二次衬砌轴力分配比例/% 24.5 26.4 19.9 32.7
    下载: 导出CSV 
    | 显示表格

    表 6中各项数据可以看出,失稳破坏阶段双层衬砌整体的内力与变形量值均处于较高水平,二次衬砌轴力占比低于40%,说明在二次衬砌的辅助下,管片衬砌的承载能力得到了强化,管片衬砌为主要承载结构。

    地层侧压力系数增加使得双层衬砌能够承受更大的外部荷载。侧压力系数从0.5增长至0.8,二次衬砌失效时的竖向外荷载增加了66.39%,管片衬砌所承受的最大轴力增加117.96%,最大弯矩减小21.35%。此时双层衬砌的最大径向位移反而从15.17 mm减小至9.45 mm,最大偏心距从220.77 mm减小至79.66 mm,表明结构的失稳变形减小、受力状态有利,结构可靠性增加。其中,地层侧压力系数由0.6变为0.7时变化最为明显,双层衬砌结构承载能力提高了59.43%,侧压力系数为0.6时,双层衬砌结构所能承受的最大竖向外荷载虽然没有改变,但最大径向位移与最大偏心距均大幅减小,很好地改善了双层衬砌结构的受力状态与变形。

    试验过程中,存在管片衬砌结构内力突变先于二次衬砌的现象,究其原因,是二者的刚度存在差异,两者共同承载的过程中,刚度较低的管片衬砌通过变形将部分外部荷载传递给二次衬砌(管片衬砌与二次衬砌的相互作用),完成了其自身的内力调整;二次衬砌则由于刚度较大,通过变形进行内力调整的能力相对较弱,在其荷载持续增加的情况下,内力随之增加,达到其承载极限时会出现损伤及局部破坏,进而释放变形,导致其内力突变。两者发生内力的突变的本质不同,管片衬砌结构的内力突变通过与二次衬砌接触传递荷载实现,二次衬砌结构的内力突变则通过结构损伤破坏实现。

    试验全过程利用声发射仪器监测不同地层侧压力系数下随加载步变化的声发射事件率及累计声发射事件数,从能量角度揭示双层衬砌的损伤过程以及地层侧压力系数改变带来的影响。如图 10,为侧压力系数0.5,0.6,0.7,0.8时声发射事件率及累计声发射事件数随加载步变化曲线。

    图  10  不同侧压力系数下加载全过程声发射信息
    Figure  10.  Acoustic emission information of whole loading process with different lateral pressure coefficients

    图 10可知,侧压力系数为0.5时,管片衬砌在第6~7级荷载下从弹性阶段过渡到了弹塑性阶段,施作二次衬砌后声发射信号趋于平稳,直到施加第14级荷载后,双层衬砌整体失稳破坏,此时累计声发射事件数达836次。随着侧压力系数的增大,管片衬砌单独承载情况下进入弹塑性阶段需要的荷载值不断增大,同时双层衬砌失稳破坏的临界荷载也逐渐增大。当侧压力系数增至0.8时,整体失稳破坏累计声发射事件数达到1236次,增加了47.8%,认为破坏前双层衬砌整体吸收了更多能量而产生了更多裂纹,二次衬砌的施作抑制了管片衬砌进入塑性阶段的进程,双层衬砌结构承载能力得到提升。

    为研究不同地层侧压力系数下双层衬砌结构破坏形态的变化规律,在加载试验结束之后对辅助承载的二次衬砌进行清理,观察作为主要承载结构的管片衬砌的宏观破坏裂纹,结合加载全过程声发射信息最终绘制出管片衬砌破坏形态素描图。如图 11所示,素描图以拱底为起点沿顺时针绘制,数字表示破坏区域编号,红色实线代表管片衬砌外侧破坏裂纹,绿色虚线代表内侧破坏裂纹。

    图  11  管片衬砌破坏形态素描
    Figure  11.  Sketch of damage pattern of segmental lining

    图 11可知,在侧压力系数为0.5的条件下,管片衬砌破坏裂纹早期主要出现在拱顶位置,后向拱腰位置延伸,裂纹具体表现为拉伸破坏和剪切破坏以及局部掉块现象。随着地层侧压力系数的增长,二次衬砌破坏裂纹数量少量增长,裂纹整体发育程度下降,局部掉块现象减少,剪切与拉伸破坏产生的裂纹由左右拱腰向拱顶转移。在侧压力系数为0.8条件下,管片衬砌局部掉块现象消失,仅在左右拱腰处有少数剪切破坏及拉伸破坏裂纹出现。可见地层侧压力系数的增加,延缓了管片衬砌局部裂纹的产生,提高了双层衬砌的承载能力以及结构的整体安全性。

    依托狮子洋水下铁路盾构隧道工程,通过模型试验对不同地层侧压力系数下的双层衬砌力学特性进行分析,得到3点结论。

    (1)在竖向土荷载恒定的情况下,地层侧压力系数增加对双层衬砌结构的极限承载能力有显著的提升效果。高地层侧压力下,双层衬砌的极限轴力值增大,弹塑性临界点与内力突变点滞后,二次衬砌补强了管片衬砌的承载性能,双层衬砌整体更加趋于静水压力场状态。

    (2)内外衬砌联合承载的条件下,二次衬砌与地层的侧向压力共同抑制了管片衬砌的侧向变形,使得管片衬砌并具备更高的抗变形能力。地层侧压力系数增大使得双层衬砌最大失稳位移减小,结构横椭变现象减弱,提高了衬砌整体结构的可靠性。

    (3)声发射率和声发射事件数的突变时机与的双层衬砌结构的破坏时机相对应。在高侧向土压力条件下,二次衬砌能够更好地协助管片衬砌承担外部荷载,延缓结构声发射率和声发射事件数的突变,减少管片衬砌大裂纹与局部掉块的产生。

  • 图  1   三角形单元图

    Figure  1.   Triangular element

    图  2   剖分网格图

    Figure  2.   Diagram of partial grids

    图  3   单元边界编号

    Figure  3.   Number of element boundary

    图  4   坝体剖分网格Ⅰ

    Figure  4.   Diagram of dam body meshing Ⅰ

    图  5   自由面

    Figure  5.   Free surface

    图  6   坝体剖分网格Ⅱ

    Figure  6.   Diagram of dam body meshing Ⅱ

    图  7   闸基模型

    Figure  7.   Model for brake base

    图  8   闸基剖分网格图

    Figure  8.   Meshes of gate base

    图  9   防渗墙深5 m时水头等值线分布图

    Figure  9.   Distribution of water head contour line at depth of anti-seepage wall of 5 m

    图  10   防渗墙深10 m时水头等值线分布图

    Figure  10.   Distribution of water head contour line at depth of anti-seepage wall of 10 m

    图  11   防渗墙深15 m时水头等值线分布图

    Figure  11.   Distribution of water head contour line at depth of anti-seepage wall of 15 m

    表  1   剖分单元信息

    Table  1   Information of elements

    单元编号 单元类型 节点编号
    1 3 3 2 4
    2 3 3 1 2
    3 4 3 4 5 6
    4 5 3 6 7 8 1
    下载: 导出CSV

    表  2   节点坐标信息

    Table  2   Information of node coordinates

    节点 1 2 3 4 5 6 7 8
    x 1 2 1 2 2 1 0 0
    y 0 0 1 1 2 2 2 0
    下载: 导出CSV

    表  3   计算结果对比Ⅰ

    Table  3   Comparison of calculated results Ⅰ

    横坐标 解析解 FEM 误差1 WG 误差2
    0 10.00 10.00 0.00 10.00 0.00
    1 9.59 9.73 0.14 9.62 0.03
    2 9.17 9.39 0.22 9.26 0.09
    3 8.72 8.99 0.27 8.79 0.07
    4 8.25 8.53 0.28 8.30 0.05
    5 7.75 8.03 0.28 7.78 0.03
    6 7.21 7.46 0.25 7.23 0.02
    7 6.63 6.83 0.20 6.63 0.00
    8 6.00 6.11 0.11 6.04 0.04
    9 5.29 5.21 -0.08 5.34 0.05
    10 4.47 4.50 0.03 4.50 0.03
    下载: 导出CSV

    表  4   计算结果对比Ⅱ

    Table  4   Comparison of calculated results Ⅱ

    横坐标 解析解 WG 误差
    0 10.00 10.00 0.00
    1 9.59 9.59 0.00
    2 9.17 9.18 0.01
    3 8.72 8.77 0.05
    4 8.25 8.23 -0.02
    5 7.75 7.74 -0.01
    6 7.21 7.16 -0.05
    7 6.63 6.61 -0.02
    8 6.00 6.00 0.00
    9 5.29 5.34 0.05
    10 4.47 4.50 0.03
    下载: 导出CSV

    表  5   水力梯度值

    Table  5   Values of hydraulic gradient

    防渗墙
    深/m
    防渗墙底端处水力梯度 闸基出口处
    水力梯度
    5 0.323 0.332
    10 0.333 0.262
    15 0.358 0.190
    下载: 导出CSV
  • [1]

    BASSI F, REBAY S. A high-order accurate discontinuous finite element method for the numerical solution of the compressible navier-stokes equations[J]. Journal of Computational Physics, 1997, 131(2): 267-279. doi: 10.1006/jcph.1996.5572

    [2]

    COCKBURN B, SHU C W. The local discontinuous Galerkin method for time-dependent convection-diffusion systems[J]. SIAM Journal on Numerical Analysis, 1998, 35(6): 2440-2463. doi: 10.1137/S0036142997316712

    [3]

    COCKBURN B, GOPALAKRISHNAN J, LAZAROV R. Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems[J]. SIAM Journal on Numerical Analysis, 2009, 47(2): 1319-1365. doi: 10.1137/070706616

    [4] 何朝葵, 速宝玉, 盛金昌. 稳定渗流分析的局部间断伽辽金有限元法[J]. 河海大学学报(自然科学版), 2012, 40(2): 206-210. https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201202018.htm

    HE Zhaokui, SU Baoyu, SHENG Jinchang. Local discontinuous Galerkin finite element method for steady seepage analysis[J]. Journal of Hohai University (Natural Sciences), 2012, 40(2): 206-210. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201202018.htm

    [5]

    WANG J P, YE X. A weak Galerkin finite element method for second-order elliptic problems[J]. Journal of Computational and Applied Mathematics, 2013, 241: 103-115. doi: 10.1016/j.cam.2012.10.003

    [6] 王军平, 叶秀, 张然. 弱有限元方法简论[J]. 计算数学, 2016, 38(3): 289-308. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSX202001001.htm

    WANG Junping, YE Xiu, ZHANG Ran. Basics of weak Galerkin finite element methods[J]. Mathematica Numerica Sinica, 2016, 38(3): 289-308. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSX202001001.htm

    [7] 朱弘泽, 林莉, 周晨光, 等. 弱Galerkin有限元法的稳定性[J]. 吉林大学学报(理学版), 2018, 56(6): 1427-1430. https://www.cnki.com.cn/Article/CJFDTOTAL-JLDX201806025.htm

    ZHU Hongze, LIN Li, ZHOU Chenguang, et al. Stability of weak Galerkin finite element method[J]. Journal of Jilin University (Science Edition), 2018, 56(6): 1427-1430. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JLDX201806025.htm

    [8]

    WANG J P, YE X. A weak Galerkin finite element method for the stokes equations[J]. Advances in Computational Mathematics, 2016, 42(1): 155-174. doi: 10.1007/s10444-015-9415-2

    [9]

    MU L, WANG J P, YE X, et al. A numerical study on the weak Galerkin method for the Helmholtz equation[J]. Communications in Computational Physics, 2014, 15(5): 1461-1479. doi: 10.4208/cicp.251112.211013a

    [10]

    MU L, WANG J P, YE X, et al. A weak Galerkin finite element method for the maxwell equations[J]. Journal of Scientific Computing, 2015, 65(1): 363-386. doi: 10.1007/s10915-014-9964-4

    [11] 张然. 弱有限元方法在线弹性问题中的应用[J]. 计算数学, 2020, 42(1): 1-17. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSX202001001.htm

    ZHANG Ran. Weak Galerkin finite element method for linear elasticity problems[J]. Mathematica Numerica Sinica, 2020, 42(1): 1-17. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSX202001001.htm

    [12] 张然, 翟起龙. 偏微分方程特征值问题的弱有限元方法[J]. 中国科学: 数学, 2019, 49(12): 1979-1994. https://www.cnki.com.cn/Article/CJFDTOTAL-JAXK201912015.htm

    ZHANG Ran, ZHAI Qilong. Weak Galerkin finite element methods for PDE eigenvalue problems[J]. Scientia Sinica (Mathematica), 2019, 49(12): 1979-1994. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JAXK201912015.htm

    [13] 吴梦喜, 张学勤. 有自由面渗流分析的虚单元法[J]. 水利学报, 1994, 25(8): 67-71. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB408.009.htm

    WU Mengxi, ZHANG Xueqin. Imaginary element method for numerical analysis of seepage with free surface[J]. Journal of Hydraulic Engineering, 1994, 25(8): 67-71. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB408.009.htm

    [14] 速宝玉, 沈振中, 赵坚. 用变分不等式理论求解渗流问题的截止负压法[J]. 水利学报, 1996, 27(3): 22-29, 35. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB603.003.htm

    SU Baoyu, SHEN Zhenzhong, ZHAO Jian. The cut-off negative pressure method for soling filtration problems based on the theory of variational inequalities[J]. Journal of Hydraulic Engineering, 1996, 27(3): 22-29, 35. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB603.003.htm

    [15] 王贤能, 黄润秋. 有自由面渗流分析的高斯点法[J]. 水文地质工程地质, 1997, 24(6): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG706.000.htm

    WANG Xianneng, HUANG Runqiu. Gavss point method for seepage analysis with free surface[J]. Hydrogeology and Engineering Geology, 1997, 24(6): 1-4. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG706.000.htm

    [16] 党发宁, 王晓章, 郑忠安, 等. 有自由面渗流分析的变单元渗透系数法[J]. 西北水力发电, 2004, 20(1): 1-3. https://www.cnki.com.cn/Article/CJFDTOTAL-SXFD200401001.htm

    DANG Faning, WANG Xiaozhang, ZHENG Zhongan, et al. Variable element seepage coefficient method for seepage numerical analysis with free surface[J]. Journal of NOTHUCEST RYDOOELECTNICPOWER, 2004, 20(1): 1-3. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SXFD200401001.htm

    [17] 孙伟建, 侯兴民, 李远东, 等. 基于虚单元法求解渗流自由面的曲线拟合法[J]. 水力发电, 2016, 42(11): 50-53, 67. https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201611013.htm

    SUN Weijian, HOU Xingmin, LI Yuandong, et al. A curve fitting method based on virtual element method for solving seepage free surface[J]. Water Power, 2016, 42(11): 50-53, 67. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201611013.htm

  • 其他相关附件

图(11)  /  表(5)
计量
  • 文章访问数:  0
  • HTML全文浏览量:  0
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-16
  • 网络出版日期:  2023-02-23
  • 刊出日期:  2023-06-30

目录

/

返回文章
返回