利用智能节段技术的物联网监测方法对盾构隧道纵向接头进行监控:通过可测量的节段旋转角度推断出隐藏的结构状态
《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》:Smart segment-enabled IoT monitoring of shield tunnel longitudinal joints: Inferring hidden structural states from measurable segment rotations
【字体:
大
中
小
】
时间:2026年08月29日
来源:TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY 8.7
编辑推荐:
**摘要:**
在盾构隧道中,接头是脆弱的组成部分。就纵向接头而言,诸如开口角度和混凝土损伤等关键状态在运行过程中仍然是隐藏的。直接测量这些状态需要在接头上放置传感器,但这些传感器在恶劣条件下可能会出现故障,甚至可能损坏防水系统。本研究考虑了从位于受保护位置的段内传感器推断隐
**摘要:**
在盾构隧道中,接头是脆弱的组成部分。就纵向接头而言,诸如开口角度和混凝土损伤等关键状态在运行过程中仍然是隐藏的。直接测量这些状态需要在接头上放置传感器,但这些传感器在恶劣条件下可能会出现故障,甚至可能损坏防水系统。本研究考虑了从位于受保护位置的段内传感器推断隐藏接头状态的问题,而不是在接头上放置传感器。这种方法基于物理基础的传递函数,将可测量的段旋转与隐藏的接头状态联系起来。整个框架包括四个步骤:首先,全尺寸弯曲试验为表征纵向接头的非线性力矩-旋转行为提供了基础,揭示了从弹性响应到失效的不同阶段;这些结果为数值校准提供了物理真实的依据。其次,针对试验结果校准段级有限元模型,以建立经过验证的数值模型,包括材料参数、接触设置和边界条件。第三,环级参数化模拟继承这些校准参数,并改变表面荷载,以建立从嵌入式传感器位置的旋转到接头开口角度和混凝土损伤的传递函数。第四,可靠性分析框架结合了荷载不确定性,以量化不同安全状态的概率。上述受保护的传感器是通过嵌入无线倾斜传感器的智能段实现的。这些物联网(IoT)支持的智能段与无线网关和云平台集成,作为核心监测元件,能够从浇筑到运行期间持续采集数据。一年期的现场监测数据被用来评估该框架的有效性。推断出的接头开口角度约为-0.13°,相当于73.3 kPa的等效表面荷载;在假设的荷载分布下,健康、警告和临界状态的概率分别为88%、8%和4%。该框架将隧道监测从被动的变形观察转变为主动的结构状态推断,为预防性维护提供了实用的智能解决方案。
**1. 引言**
段接头被广泛认为是盾构隧道中最脆弱的组成部分。特别是,纵向接头的状态直接影响结构安全(Ma等人,2025年;Jiang等人,2024年)。由于接头开口以及螺栓腐蚀和结构逐渐破坏,隧道可能会发生渗漏,这可能影响隧道的使用性和安全性(Liu等人,2024年;Gong等人,2019年)。不利的荷载条件,如衬砌背后的灌浆空洞引起的条件,也可能促进段结构的开裂和损伤(Bao等人,2025a)。大多数现有的接头状态监测技术都需要在接头上放置传感器。在高湿度等恶劣的接头条件下,传感器可能会随时间损坏,使得长期监测变得困难。此外,在接头上安装传感器还可能损坏防水系统。因此,可靠地监测接头状态对于隧道预防性维护至关重要。
盾构隧道的结构健康监测最近受到了重视,已经开发了几种针对接头的技术(Chen等人,2025年)。传统的监测技术,如收敛度和沉降监测,记录了整体结构反应,但无法单独量化各个接头的贡献或提供有关接头内部应力状态的信息(Zhang等人,2017年)。近年来,专门为接头开发了先进的监测技术。对于变形监测,基于布里渊散射的分布式光纤传感技术能够实现沿接头的连续应变测量,具有高空间分辨率,并能检测到开合变形和位移变形(Zhang等人,2022年;Wang等人,2018年)。类似的分布式光纤传感器系统也已被开发用于浸没式隧道接头,能够检测到亚毫米级的水平接头开口和垂直不均匀沉降(Zhang和Broere,2023年)。为了监测老化隧道中的环向变形,已经开发了优化的数据处理方法,将离散的分布式光纤传感测量转换为连续的隧道收敛度记录(Xiao等人,2024年)。多传感器数据融合技术结合了旋转和应变测量,提高了变形监测的准确性,而无需参考框架(Ding等人,2024年)。对于应力和损伤监测,已经开发了压电传感器,用于在施工期间测量接头处的接触应力(Jiang等人,2024年),并通过环形接缝压力分布识别接头开口(Xu等人,2024年)。还开发了一种使用基于深度学习的图像分析方法从隧道检查图像中自动检测和量化接头渗漏的方法(Ma等人,2025年)。在系统层面,Li等人(2026年)开发了集成多源数据和早期预警系统的智能监测平台,以实现全面的隧道安全管理。此外,Afrazi等人(2026年)表明,数字孪生和人工智能技术对实时隧道监测和主动维护有益。同时,物联网(IoT)和物理信息学习已被引入隧道工程:Zhang等人(2021年)构建了一个基于人工智能物联网(AIoT)的系统,实时传输盾构运行参数和施工引起的沉降,以支持闭环控制;物理信息模型如用于段衬砌抬升的Transfer Learning DeepONet(Elbaz等人,2026年)和用于正向及反向隧道问题的三保真TF-PINN(Elbaz等人,2025年)将机械原理嵌入数据驱动的预测中。然而,这些IoT和物理信息方法针对的是整体施工响应,如沉降和抬升,而不是运行过程中单个纵向接头的隐藏状态。
为了评估这些关键接头状态,目前的监测系统仍然主要依赖于通过接触传感器、干涉相机等直接测量接头级别的量。然而,如前所述,这些直接方法要么需要在接头位置进行侵入式安装,要么只能在损伤发生后才能检测到。另一种方法是借助基于物理的传递函数,从段级变形中获取接头状态,这需要理解段变形与接头行为之间的机械关系。物理模型中的弯曲试验表明,纵向接头的力矩-旋转关系表现出高度的非线性,通常可以分为几个阶段,包括弹性行为、软化和失效(Feng等人,2026年;Wang等人,2020年;Gong等人,2017年;Li等人,2015年)。通过使用受混凝土损伤塑性(CDP)模型的数值模拟,可以捕捉接头的渐进性损伤演变(Li等人,2024年;Liu等人,2022年)。这些机械研究主要集中在设计方面,但如何利用这些知识进行实时监测尚未研究。换句话说,目前还没有基于物理的框架可以从可测量的段变形中推断隐藏的接头状态,而不需要在接头上放置传感器。这一差距推动了基于物理的推断方法的发展,使得可以从段级测量中评估接头状态。
为了减少施工干扰并促进全生命周期监测,提出了智能段的概念,即配备传感器的预制混凝土元件,这些传感器可以在其整个使用寿命期间跟踪其性能(Tshimbombo等人,2023年;Tshimbombo,2023年)。这些段避免了在施工期间安装监测设备的需要,实现了从浇筑到运行期间的连续监测。在各种传感技术中,无线倾斜传感器具有特殊优势,如最小的操作干扰、可靠性、易于安装以及无需清晰视线即可提供实时反馈(Zhang等人,2024年;Bennett等人,2010年;Ko等人,2010年;Lewis,2004年)。已经开发了基于无线倾斜传感器的实时水平收敛度监测系统,并优化了传感器嵌入位置,以从段旋转计算收敛度(Chang等人,2022年;Wang等人,2021年)。然而,智能段的现场应用仍然有限,它们作为间接传感工具来推断隐藏接头状态的研究尚未进行。
为了支持地下基础设施的预防性维护,本研究提出了一个基于物联网的框架,可以从可测量的段旋转中推断隐藏的接头状态。该框架结合了全尺寸弯曲试验来研究非线性接头行为,校准的有限元(FE)模型来建立物理传递函数,以及基于物联网的智能段进行实时现场监测。还开发了可靠性分析框架,以计算在荷载不确定性存在下的安全概率。与现有的需要在接头上放置传感器或只能在损伤发生后才能检测到损伤的接头监测方法不同,当前框架通过实验校准的物理传递函数,从受保护的段内位置测量的旋转中推断出运行中盾构隧道的隐藏纵向接头状态,并进一步将推断的状态转换为荷载不确定性下的概率安全状态。这项工作的主要贡献包括:(1)一个基于物理的推断框架,它将可测量的段旋转与隐藏的接头状态联系起来,而无需在接头位置放置传感器;(2)一年的现场验证表明其适用于隧道预防性维护;以及(3)一个考虑荷载不确定性的可靠性分析框架,用于概率安全评估。
**2. 框架概述**
本研究提出了一个基于物联网的智能监测框架,将可测量的段旋转与隐藏的纵向接头状态联系起来。如图1所示,该框架包括四个相互连接的组件,具有渐进的信息流:(a)全尺寸测试,(b)段级FE模型校准,(c)环级参数化模拟,以及(d)结合可靠性分析的基于物联网的智能段。这些组件之间的信息传递建立了一条从实验表征到实时现场推断的完整链:
- **全尺寸测试**通过表征非线性接头行为建立了物理基础,如图1(a)所示。对全尺寸纵向接头进行负弯曲测试表明,它们会经历六个阶段,从弹性到最终失效。测试研究展示了开口角度与混凝土损伤之间的关系的发展。力矩-旋转(M-φ)关系是主要输出。实验确定的M-φ关系作为真实依据,并传递(图1中的链接I)到段级FE模型进行数值校准。
- **段级FE模型校准**过程包括使用全尺寸测试中的M-φ关系来建立和验证接头行为的数值模型,如图1(b)所示。一个整合了CDP本构行为和固体元件螺栓建模的FE模型被校准,以复制实验观察到的M-φ关系。这种校准验证了模型捕捉接头刚度退化和失效机制的能力。然后,校准后的模型参数(包括本构关系、接触设置和边界条件)被传递(图1中的链接II)到环级模型进行参数化模拟。
- **环级FE模型参数化模拟**利用段级校准的验证建模方法,并基于荷载-结构方法将其扩展到三环配置,如图1(c)所示。使用截断的对数正态分布对表面荷载进行建模,以考虑荷载不确定性,并在操作荷载范围内模拟多种情景。每种情景生成包含传感器位置段旋转作为可测量输入和接头状态(具体为开口角度和加权平均混凝土损伤)作为隐藏输出的数据对。通过回归分析,建立了将可测量的段旋转与隐藏接头状态联系起来的定量传递函数。这些基于物理的传递函数随后被应用于现场监测系统(图1中的链接III)进行实时推断。
- **物联网监测和可靠性评估**将参数化模拟获得的传递函数直接应用于现场测量,如图1(d)所示。带有嵌入式无线倾斜传感器的智能段在隧道运行期间持续测量旋转。物理传递函数使得可以直接从这些测量中推断隐藏的接头状态:将嵌入式传感器位置的旋转转换为段边缘旋转,计算接头开口角度,并评估相应的损伤状态。可靠性分析框架进一步结合荷载不确定性,以评估不同安全状态(健康/警告/临界)的概率,提供概率置信水平以支持基于风险的维护决策。最终输出包括推断的接头开口角度和被监测隧道段的安全状态概率。
**3. 全尺寸测试和段级FE模型校准**
所考虑的段尺寸如下:外径为8.8米,内径为8.0米,段宽度为1.8米。每个环由七个块组成,包括一个关键块(F)、两个相邻块(L)和四个标准块(B)。使用M30斜向螺栓进行连接,其中14个周向螺栓将七个块连接在一起,19个纵向螺栓连接周向接头。本节报告的满尺寸测试以及与之校准的节段级有限元模型(图2、图3、图4、图5)最初由Feng等人(2026年)提出,他们利用这些模型来量化水平轴向力和顶升力对纵向接头弯曲性能的耦合效应。在本研究中,这些模型具有不同的用途,为第4节建立的传递函数提供了实验和数值基础。
下载:下载高分辨率图像(409KB)
下载:下载全尺寸图像
图2. 满尺寸测试概述(Feng等人,2026年)
下载:下载高分辨率图像(432KB)
下载:下载全尺寸图像
图3. 满尺寸测试的加载方案和仪器布置
下载:下载高分辨率图像(299KB)
下载:下载全尺寸图像
图4. 纵向接头开角与负弯矩之间的关系(Feng等人,2026年)
下载:下载高分辨率图像(294KB)
下载:下载全尺寸图像
图5. 测试与有限元模型之间弯矩-旋转关系及失效模式的比较(Feng等人,2026年)
3.1. 测试配置和加载协议
测试需要选择三个载荷分量:单个环段内块与块之间的水平载荷N、环段之间的纵向相互作用力LT,以及通过接头两侧的垂直载荷P施加在纵向接头上的负弯矩。在本研究中,弯矩的符号约定基于接头变形:负弯矩导致外缘闭合而内缘打开,正弯矩则产生相反的效果(外缘打开而内缘闭合)。满尺寸测试装置以及详细的载荷和边界条件分别显示在图2(a)和(b)中。虽然纵向顶升力和垂直载荷不会直接引起末端块E1和E2的水平位移,但仍可能发生微小的水平移动。因此,施加水平载荷N以抑制这种移动。在本测试中,单个环段内块与块之间的水平载荷为80 kN。
环段之间的纵向相互作用力是根据Liao等人(2017年)报告的应力松弛现象确定的,他们发现环段之间的接触应力松弛经历了四个发展阶段:(1) 盾构隧道施工期间的动态调整;(2) 隧道内部施工期间的临时平衡;(3) 隧道内部施工后的逐渐减小;(4) 最终稳定。在这四个阶段之后,环段之间的接触应力从大约2 MPa逐渐减小到1 MPa,平均应力松弛程度约为50%。基于这种应力松弛行为,通过将1 MPa的稳定状态降低50%来设定失效测试的接触压力为0.5 MPa。与完全失效状态(接触为零)不同,这个值代表了一个保守的临界情景,在这种情况下,环段仍然通过摩擦、榫槽互锁和螺栓预紧等仍然活跃的耦合机制保持接触。在0.987 m2的受力面积上,这对应于493.5 kN的纵向顶升力??T。纵向顶升力分三个阶段施加(200 kN、200 kN和93.5 kN,即总计493.5 kN),施加速率为100 kN/min。对于垂直加载,每侧的力量P以15 kN/min的速率递增,直至270 kN,然后以10 kN/min的速率递增,直至峰值载荷约300 kN。总共,纵向顶升力施加了3个水平,垂直力施加了11个水平,如图3(a)所示。关于仪器,图3(b)详细说明了传感器布置,其中面板(b-I)和(b-II)分别描述了环段内缘和外缘的布置。该系统使用四个点来监测接头开角:两个在外缘测量开角,两个在内缘记录闭合情况。位移数据以1 Hz的频率记录。
弯矩计算遵循Gong等人(2017年)提出的方法,但在此研究中未对纵向接头施加水平轴向力N。考虑到中间环段的加载条件,弯矩分布如图3(a)所示,纵向接头的弯矩使用公式(1)计算:
(1)
?{ {?{ {???A=???(?????1)??B=????/2??=??A+??B
其中P表示作用在接头一侧的垂直载荷。e是从纵向螺栓中心到纵向接头的距离,等于693 mm。标准块M1和M2的水平尺寸为1300 mm。从两侧的垂直载荷到纵向接头的距离??1为400 mm。标准块M1和M2的弹性模量G为21.75 kN。需要注意的是,公式(1)仅考虑了中间环的外部载荷,并保守地忽略了由摩擦、榫槽互锁和螺栓预紧等环间耦合机制提供的弯曲约束,这些机制的作用在第3.3节中进一步研究。
3.2. 纵向接头的弯矩-旋转关系
基于满尺寸测试,分析在施加负弯矩下的纵向接头的弯矩-旋转关系。从弯矩-旋转分析中排除了负荷保持阶段的数据。然后使用公式(2)(Li等人,2015年)计算纵向接头的旋转:
(2)
??=??2???1?
其中??2表示图3(b)中位移传感器??2,1和??2,2测得的平均开角,??1表示图3(b)中位移传感器??1,1和??1,2测得的平均闭合角,所有单位均为毫米;h是环段厚度,等于400 mm。纵向接头开角与负弯矩之间的关系如图4所示,并根据六个不同阶段进行解释:
(1) 初始弹性阶段:弯矩-旋转关系大致呈线性,弯曲刚度很高。混凝土和钢螺栓均处于线性弹性状态,接头界面完好无损,没有可见裂纹。
(2) 非线性过渡阶段:从线性行为向非线性行为转变。随着接触压力分布的变化,切线刚度开始减小,局部材料非线性开始显现。
(3) 第一个塑性平台阶段:曲线出现明显的平台,弯曲刚度降至几乎为零(类似于完全塑性响应)。这归因于混凝土压缩区的局部压碎和螺栓的屈服,导致旋转增加但弯矩增益可以忽略不计。
(4) 刚度恢复阶段:曲线的斜率再次增加,表明弯曲刚度恢复。这种现象是由于螺栓的应变硬化和混凝土压缩区的压实,使接头在初始屈服平台后能够承受更高的载荷。
(5) 第二个塑性平台阶段:刚度再次显著减小,接近零,曲线进入第二个平台。螺栓发生广泛的塑性流动,混凝土达到其极限压缩应变,导致无限制的塑性旋转。
(6) 最终失效阶段:由于螺栓的最终应变硬化,弯矩略微增加至峰值。该阶段在图4中的红色圆圈和内嵌照片所示的失效点突然结束,这是由于螺栓达到其极限抗拉强度和混凝土的广泛压碎造成的。
3.3. 基于满尺寸测试结果验证节段级有限元模型
基于满尺寸测试结果,开发了一个节段级有限元模型来模拟在493.5 kN顶升力作用下的纵向接头弯曲行为,如图2(b)所示。表1中的材料参数采用来自C55混凝土的现有文献(GB 50010–2010,2010年,Li等人,2024年,Tian等人,2025年)。通过将模拟的M-φ关系与实验结果进行比较,确认了这些参数的有效性,如图5(a)所示,模拟曲线再现了测试的整体演变趋势和弯矩能力。螺栓锚固端连接到节段上。节段之间以及节段与螺栓之间采用表面接触,摩擦系数为0.3,法向为硬接触。
表1. 材料参数(Feng等人,2026年)
| 材料 | 属性 | 值 | |
|--------------|-----------------|----------------------|--------------------|
| C55混凝土 | 标准模量 | E_c (MPa) | 2.4 |
| | 泊松比 | ?_c (f) | 0.2 |
| | 密度 | ?_c (t/m3) | 3.3 |
| | 抗压强度 | f_ck (MPa) | 35.5 |
| | 抗拉强度 | f纺织品 (MPa) | 27.4 |
| | 扩张角 | ??/° | 0.0 |
| | 黏度系数 | ? | 0.6 |
| | 屈服表面参数 | E_c (MPa) | 0.6 |
| | 偏心率 | f_b?0?/??_c?0 | 3.8 |
| | 螺栓 | E_b (MPa) | 21.0 |
| | 泊松比 | ?_b (f) | 0.6 |
| | 抗拉强度 | f тек (MPa) | 35.5 |
| | 极限强度 | f_u (MPa) | 64.0 |
| | | | |
FE模型施加了与满尺寸测试一致的边界条件和载荷。如图5(a)所示,节段级FE模型捕获了测试的整体响应。然而,观察到与实验结果的某些偏差。首先,在弹性阶段(第1阶段),FE模型的旋转量略小于测试在相同时刻的水平。相反,在弹性极限之后(第2至4阶段),FE模型产生的旋转量大于测试,表明数值响应略软。这种转变主要是由于物理测试不可避免地涉及支撑和加载系统中的寄生摩擦约束,这些约束抵抗了大变形,而理想化的数值边界条件并未捕捉到这些因素。此外,CDP模型的连续介质损伤公式导致弹性极限后的刚度减小更为渐进。其次,FE分析的第(5)阶段没有显示出测试中观察到的急剧刚度下降;相反,弯矩-旋转关系大致保持线性,弯曲刚度几乎恒定为322.7 kN·m/rad,与测试中接近零的刚度(塑性平台)相反,因为CDP模型模糊了局部混凝土压碎现象,而螺栓遵循理想的弹性-塑性规律。应当注意的是,这种比较的目的是为了验证传递到环级别模型的本构关系、接触设置和螺栓建模,而不是逐阶段再现测试曲线。第4节中的环级别模拟得到的接头开角幅度不超过0.323°,对应于曲线早期偏差最小的部分,而第(5)阶段的平台出现在旋转量约为1.8°以上。因此,第4.3节推导出的传递函数和随后的接头状态推断在研究的操作范围内是可靠的。
图5(b)比较了M1块和M2块中靠近E2块的榫槽部分的混凝土损伤模式。由于测试考虑了环与环之间的相互作用,纵向接头经历弯曲,而周向接头在垂直载荷下抵抗剪切;周向接头中的榫槽连接部分的混凝土主要受到压缩。CDP模型中的压缩损伤指数用于量化局部微裂纹的传播和刚度下降。这个标量损伤变量随着非弹性应变的累积而变化,范围从0(完整材料)到1(完全失去压缩能力)。数值结果与满尺寸测试观察结果一致,表明周向接头在纵向接头弯曲期间抵抗剪切,榫槽混凝土中观察到显著的压缩-剪切破坏。这表明周向接头的剪切性能直接影响纵向接头的弯曲能力。
验证后的节段级FE模型建立了纵向接头行为的本构关系、接触设置和边界条件。这些建模参数随后被转移到第4节的环级别FE模型中,用于在不同表面附加荷载情况下进行参数模拟。
4. 隐藏状态推断的物理传递函数
本节建立了将可测量的节段旋转与隐藏接头状态联系起来的物理传递函数。开发过程分为两个阶段:首先,介绍环级别FE模型和在不同表面附加荷载情况下的参数模拟;其次,推导节段旋转与接头状态之间的定量关系,包括从传感器位置到节段边缘的传递函数、接头开角计算和混凝土损伤演变。
4.1. 监测截面的地质条件
所研究的盾构隧道的埋深范围为11.6 m至38.8 m。在监测截面,隧道顶部覆盖有淤泥、泥泞土壤和从表面向下沉积的固结泥质粘土层,盾构隧道主要穿过中粗砂层。在监测横截面处的地质剖面如图6(a)所示,从现场调查报告中获得的各土层的物理和力学参数总结在表2中。下载:下载高分辨率图像(521KB)下载:下载全尺寸图像图6. 基于荷载-结构方法的环形FE模型概览。表2. 监测横截面各土层的物理和力学参数。编号 土层 淹没单位重量 ???? (kN/m3) 厚度 ???(m) 内聚力 ????(kPa) 内摩擦角 ???? (°) 周围岩石分类 1 泥沙 6.3 18.4 6 5.0 5.0 6 2 泥质土壤 6.7 3.1 8 7.0 6.0 3 固结泥质粘土 7.9 7.7 5 10.0 8.0 4 中粗砂 9 1.3 3 0.0 3 2.0 5 隧道主体 4.2. 基于荷载不确定性的环形FE模型参数化仿真本小节开发了用于参数化仿真的环形FE模型,涵盖了模型配置(包括几何形状、材料属性、边界条件及荷载方案),并使用截断的对数正态分布来表示隧道运行期间表面附加荷载的变化性。4.2.1. 环形FE模型概览环形FE模型采用了第3节中节段级校准验证的建模方法,包括接触设置和材料参数。监测横截面处的地质剖面如图6(a)所示,施加在结构上的荷载如图6(b)所示。在图6(b)中,活荷载(即未受扰动的原位应力)包括:??v?1 是顶部处的垂直土压力和表面附加荷载;??h?1 和 ??h?2 分别是顶部和底部处的水平土压力和表面附加荷载;??0 是结构自重;??w?1 和 ??w?2 分别是顶部和底部处的静水压力。地面反力包括:??r?1 是由垂直土水压力引起的岩层反力;??r?2 是岩层的水平阻力;??g 是由节段自重引起的岩层反力(Wei等人,2025年,Wei等人,2024年,Liu等人,2022年)。如图6(c)所示,??r?1、??r?2 和 ??g 通过弹簧连接衬砌与地面,其中正常弹簧仅承受压缩荷载,而切向弹簧提供摩擦阻力。由于隧道主体位于中粗砂中,并根据静态圆锥贯入试验,地面反力系数K取为20 MPa/m。正常弹簧的刚度由网格节点数确定。切向弹簧的刚度设置为正常弹簧的三分之一(Arnau和Molins,2012年)。环形FE模型包含三个节段环,以考虑环与环之间的相互作用效应,如图6(d)所示。节段细节,包括齿槽连接,如图6(e)所示。对于FE分析,屏蔽隧道节段的纵向位移(沿环形FE模型的Z轴)受到限制。数值模拟包括两个分析步骤:步骤1涉及对螺栓施加100 kN的预紧力。步骤2涉及施加土水压力并激活地面弹簧以完成分析。4.2.2. 荷载不确定性为了捕捉基于可靠性的预警系统中表面附加荷载的变化性,使用截断的对数正态分布来表征荷载不确定性(Ang和Tang,2007年,Baecher和Christian,2003年)。选择对数正态分布是因为其非负支持和右偏特性,反映了从健康状态到预警状态再到临界荷载范围的概率递减。第5.2.2节检验了结果对分布参数的敏感性。表面附加荷载S的概率密度函数表示为:(3)???(??)=1/??·??ln???√2πexp?[?(ln??????ln???)2/2??2ln???]·1/???(??max),0 ????max,其中:S是表面附加荷载,单位kPa;??ln???和??ln???分别是ln???的均值和标准差;??max是表面附加荷载S的截断上限。???(??max)是在??max处评估的对数正态分布的累积分布函数(CDF),它对截断范围(0, ??max)内的密度进行归一化。给定目标均值????和变异系数COV??,这些参数分别由公式(4)和(5)确定:(4)??ln???=√ln?(1+COV2/??)(5)??ln???=ln?(????)???2ln???2在本研究中,表面附加荷载的均值取为???? = 75 kPa。该值对应于大约4米高的临时填料堆,单位重量约为19 kN/m3,这是软土中屏蔽隧道的代表性偶然附加荷载情况,并且与报告的附加荷载事件和全尺寸堆放试验中的荷载水平一致(Wei等人,2024年,Huang和Zhang,2016年)。变异系数设为COV?? = 0.5,因为材料堆放、交通和施工活动引起的临时附加荷载变化很大。第5.2.2节量化了这些参数对可靠性结果的影响,因此结论不依赖于单一参数的选择。因此,得到的分布参数为??ln??? = 4.206和??ln??? = 0.472。上限设定为??max = 150 kPa,基于以下考虑:(1)现场规定通常在隧道运行期间对表面附加荷载进行限制以确保结构安全;(2)所提出的框架侧重于运营预警而非极限承载能力评估,因此超出实际管理范围的极端超载情况不太相关。图7展示了截断的对数正态分布的概率密度函数。第4.3.3节的工作定义了预警阈值(107.46 kPa)和临界阈值(129.79 kPa),从而确定了三个荷载范围,即健康(0–107.46 kPa)、预警(107.46–129.79 kPa)和临界(129.79–150 kPa),分别对应健康、预警和临界结构状态。大约88%的概率质量位于健康范围内,8%位于预警范围内,4%位于临界范围内。进行了参数化仿真,以评估表面附加荷载从0到150 kPa(以5 kPa为增量)的结构响应,共得到31种荷载条件,覆盖健康、预警和临界荷载范围。重要的是要认识到,在实际中,这些荷载状态描述的发生概率将显著较低(即远小于1),对于预警和临界状态而言。由于仿真的均匀采样(5 kPa增量),转移函数对每种荷载条件都有效。下载:下载高分辨率图像(184KB)下载:下载全尺寸图像图7. 用于荷载不确定性表征的表面附加荷载的截断对数正态分布。4.3. 节段旋转与关节状态之间的转移函数本节建立了节段旋转与关节状态之间的转移函数。从环形FE模型中提取节段中跨和边缘的旋转角度并进行分析。然后进行回归分析,以确定节段中跨的旋转角度与其对应的节段边缘旋转之间的关系,以及关节开口角度与这些边缘旋转和表面附加荷载之间的关系,还有混凝土损伤的演变。这些转移函数构成了第5节中介绍的基于物理的IoT监测系统和可靠性分析框架的基础,使得可以从现场测量数据实时推断隐藏的关节状态。为了从FE模型中提取旋转数据,计算每个网格元素的旋转量作为变形前后的坡度差(Wang等人,2021年)。元素表示环形FE模型中的网格划分结果。每个元素的起点和终点按顺时针方向在拱顶处编号,从拱顶开始沿圆周方向编号。逆时针旋转定义为正,而顺时针旋转定义为负(Wang等人,2021年)。图8(a)展示了从零水平到150 kPa的31个表面附加荷载水平下节段的周向旋转。为了展示图8(a)中节段旋转的物理意义,图8(b)展示了0 kPa时的旋转模式,该模式下有四个零旋转点位于圆周角度0°、79.95°、170.64°和281.78°。这将环分为四个旋转区域,第一个区域为正,第二个区域为负,第三个区域为正,第四个区域为负。应该指出的是,零旋转点关于拱顶-底部轴不对称。这是由于节段的不对称组装造成的。尽管这些零旋转点的确切位置会随着附加荷载水平的变化而略有变化,但总体模式保持不变。这种变形模式反映了盾构隧道在垂直荷载作用下的典型水平椭圆形化模式,即隧道横截面的垂直直径减小而水平直径增加。随着表面附加荷载的增加,旋转幅度逐渐增大,而空间分布基本保持不变,只有零旋转点位置发生轻微移动。这种自相似行为表明,单一的功能关系可以描述整个运行范围内的荷载-旋转行为。150 kPa时的旋转幅度约为50 kPa时的2.5–3倍,表明表面附加荷载对结构变形有显著影响。对于第5.1节中描述的监测系统,选择B1和L1块作为最佳传感器位置,因为它们表现出中等至较大的旋转幅度,并且对荷载变化具有较好的敏感性。L1和B1块的位置范围由图8(a)中的两个垂直平面表示。下载:下载高分辨率图像(564KB)下载:下载全尺寸图像图8. 表面附加荷载下的节段旋转特性。为了全面描述节段旋转与影响因素之间的关系,为B1和L1块构建了响应面,以表示旋转随表面附加荷载和节段位置角的变化情况,见图9。投影在顶平面上的等高线有助于在2D平面上可视化旋转模式。对于图9(a)中显示的B1块,它位于弹簧线上方,位置角大约在295°到351°之间,所有旋转均为负,表示顺时针旋转。旋转幅度随表面附加荷载的增加而单调增加,在150 kPa时达到约?0.68°的峰值。响应面显示出平滑的梯度,没有不连续性,表明适合基于回归的转移函数开发。响应面相对于表面附加荷载的梯度代表了旋转对荷载的敏感性。这为从节段旋转推断关节状态提供了物理基础。下载:下载高分辨率图像(489KB)下载:下载全尺寸图像图9. 不同表面附加荷载和位置角下的节段旋转。对于图9(b)中显示的L1块,它跨越弹簧线区域,位置角大约在237°到295°之间,旋转行为更为复杂。L1的上部表现出与B1相似的负旋转,而下部则表现出正旋转。这种转变反映了靠近弹簧线的隧道变形模式的拐点。最大正旋转在150 kPa时达到约0.40°。B1和L1的对比行为对监测策略有重要意义,因为相邻节段边缘之间的旋转差异直接对应于纵向关节开口角。4.3.1. 节段边缘与传感器位置之间的旋转关系为了实际应用,无线倾斜传感器嵌入在节段内的特定位置,而不是关节所在的边缘。放置在节段中跨的传感器受到混凝土覆盖的保护,从而更好地免受施工损坏。然而,这种设计需要建立转移函数,将此类传感器位置的测量旋转与计算关节开口角所需的边缘旋转联系起来。图10展示了通过本研究的31个参数化仿真案例获得的B1和L1块的转移函数。对于图10(a)中显示的B1块,根据以下公式存在边旋转??Be和传感器位置旋转??Bs之间的强线性关系:(6)??Be?(??Bs)=0.014+0.654×??Bs,??2=0.999下载:下载高分辨率图像(282KB)下载:下载全尺寸图像图10. 从传感器位置到节段边缘的转移函数。斜率系数0.654表明,边缘旋转的幅度大约是传感器旋转的65.4%,反映了关节边缘相对于中跨嵌入传感器的旋转减少。对于图10(b)中显示的L1块,边旋转??Le和传感器旋转??Ls之间的关系具有非线性特征,如下所示:(7)??Le?(??Ls)=0.172×exp?(???Ls0.121)?0.199+0.493×??Ls,??2=0.992L1的非线性特征源于其跨越弹簧线的位置,此处变形模式在正旋转和负旋转之间转变。以0 kPa的表面附加荷载为例,B1和L1块的边缘旋转分别为?0.072°和?0.055°。由于这两个值都是负的,因此表示顺时针旋转。B1相对于L1的较大旋转角度,相差-0.017°,表明B1的边缘相对于L1顺时针旋转,在纵向接头的凸缘处张开,在凹缘处闭合。4.3.2. 纵向接头开口角:与节段边缘旋转和表面荷载的关系纵向接头开口角??J定义为相邻节段边缘之间的相对旋转,如图11(a)所示。具体来说,??J是通过计算块B1的边缘旋转与块L1的边缘旋转之间的差值来得出的:(8)??J=??Be???Le下载:下载高分辨率图片(394KB)下载:下载全尺寸图片图11. 地下盾构隧道在表面荷载作用下的接头开口行为。??J的负值表示接头在凸缘处张开,在凹缘处闭合,对应于接头位置的负弯矩。这种符号约定与第3.2节中的全尺寸试验观察结果一致。图11(b)显示了接头开口角随表面荷载的变化,其关系遵循指数函数:(9)??J?(??)=0.212?0.223×exp?(??170.281),??2=0.999在表面荷载为零时,FE模型给出的基线接头开口角约为-0.017°(通过方程(8)从边缘旋转得出),代表了仅由自重和土压力引起的变形;拟合的指数方程(9)在整个范围内都很好地捕捉到了这一趋势,R2 = 0.999。随着表面荷载的增加,接头开口角逐渐增大,在150 kPa时达到-0.323°。指数形式反映了纵向接头在不断增加的弯矩下的非线性刚度降低,这与第3.2节中描述的六阶段弯矩-旋转行为一致。图11(b)中的插图展示了在150 kPa荷载下接头处的相应拉伸损伤分布。4.3.3. 混凝土损伤随表面荷载的演变虽然接头开口角提供了关于接头变形的见解,但它并不能直接量化材料的退化程度。为了全面评估结构安全性,还必须评估混凝土的损伤状态。在模拟中,定义了一个加权平均混凝土损伤??Wa:(10)??Wa=∑????=1????·????∑????=1????其中????是拉伸损伤变量,????是元素i的体积,n是元素的总数。这种体积加权公式提供了一个标量指标,总结了每个节段的整体损伤状态。图12展示了块B1和L1的加权平均损伤随表面荷载的演变。损伤演变遵循S形曲线,可以用一般形式表示为:(11)??Wa=??1???21+exp?[(?????0)/??]+??2其中??Wa是加权平均混凝土损伤,??0是拐点,k是增长率参数,??2和??1分别是上渐近值和下渐近值。基于环形FE模型结果的回归分析,块L1和B1的具体拟合关系分别由方程(12)、(13)给出:(12)??WaL=?0.016?0.0611+exp?[(???54.879)/39.928]+0.061,??2=0.999(13)??WaB=?0.014?0.0641+exp?[(???70.183)/45.261]+0.064,??2=0.999其中??WaL和??WaB分别表示块L1和B1的加权平均混凝土损伤。下载:下载高分辨率图片(240KB)下载:下载全尺寸图片图12. 块L1和B1的加权平均混凝土损伤随表面荷载的演变。为了建立安全状态分类的客观标准,通过导数分析确定了S形损伤演变曲线的特征点。对方程(11)中的一般S形损伤函数关于表面荷载S的一阶和二阶导数进行计算,得到:(14)??′Wa=d???Wad???=(??2???1)·exp?[(?????0)/??]??·{1+exp?[(?????0)/??]}2(15)??″Wa=d2???Wad???2=(??2???1)·exp?[(?????0)/??]·{1?exp?[(?????0)/??]}??2·{1+exp?[(?????0)/??]}3二阶导数??″Wa代表了损伤增长的加速。将??″Wa = 0设定为拐点S = S0,此时损伤增长率达到最大值。为了确定从快速增长到饱和的过渡点,通过设定???Wa = 0来确定最小二阶导数的位置,得到:(16)??Decel=??0+ln?(2+√3)·??其中ln(2+√3)≈1.317。这个最大减速点标志着损伤增长从减速到渐近饱和的过渡。在这一点上,损伤已经达到了其最终值的大约79%,表明结构储备能力已经显著减少。如图12所示,块L1和B1的二阶导数曲线??″WaL和??″WaB与损伤演变曲线一起绘制。由于块L1的位置更靠近弹簧线,应力集中更为严重,因此在整个加载范围内其损伤略高于B1。因此,它较早达到了损伤饱和阈值。L1的最大减速点出现在S = 107.46 kPa,而B1的最大减速点出现在S = 129.79 kPa。因此,L1的最大减速点(107.46 kPa)被采用作为警告阈值,而B1的最大减速点(129.79 kPa)被采用作为临界阈值。这些数学推导出的特征点提供了物理意义上的界限:在这些值以下,损伤继续发展,结构仍具有承载能力;在这些值以上,损伤接近饱和状态,结构接近其承载极限。图13所示的混凝土损伤的空间分布证实了这些阈值对应的物理状态。如图13(a)至(d)所示,在低表面荷载水平下,拉伸损伤主要局限于螺栓孔和接头边缘附近。随着表面荷载接近警告阈值(107.46 kPa),沿接头界面形成了拉伸损伤带(例如,在图13(e)中的100 kPa),表明裂纹扩展,降低了接头的使用性能并增加了渗漏风险。当加载接近并超过临界阈值(129.79 kPa)时,整个接头区域形成了广泛的拉伸损伤区(例如,在图13(f)和(g)中的125–150 kPa),表明结构严重退化。此外,图13(h)显示,在极端加载(150 kPa)下,接头凹缘处的最大压缩损伤达到了大约0.85,表明局部混凝土发生破碎。在整个考虑的加载范围内,计算出的螺栓应力始终低于屈服强度,因此这两个阈值分别提前预警了螺栓屈服和最终失效。这些阈值在第5.2节基于可靠性的安全状态评估框架中被采用。下载:下载高分辨率图片(895KB)下载:下载全尺寸图片图13. 混凝土损伤随表面荷载的演变:(a)-(g) 150 kPa时的拉伸损伤DAMAGET;(h) 压缩损伤DAMAGEC。本节建立的传递函数提供了从可测量的节段旋转到隐藏接头状态的完整推理链。方程(6)、(7)中的传递函数将传感器测量值转换为节段边缘旋转。然后使用方程(8)计算接头开口角,其在不同表面荷载水平下的变化由方程(9)描述。混凝土损伤的演变由方程(12)、(13)量化。这些关系构成了第5节中提出的基于物理的可靠性分析框架的基础,使得能够从实时监测数据中概率推断隐藏的接头状态。应当注意的是,这些传递函数中的系数是针对本研究考虑的隧道配置和地质条件特定的。对于具有不同节段尺寸、土壤属性或加载条件的其他隧道项目,需要重新计算这些系数。然而,这里提出的方法论提供了一个通用框架,可以应用于不同情况以获得相应的系数。5. 物联网监测和基于可靠性的安全评估本节介绍了所提出框架的现场实施,并展示了如何将第4节中建立的物理传递函数应用于实时监测数据。实施过程如下:无线倾斜传感器嵌入在智能节段中,连续测量块B1和L1内预定义位置的旋转;然后将这些测量值转换为节段边缘旋转;随后计算接头开口角;最后,使用S形函数评估相应的损伤状态。可靠性分析框架进一步整合了加载不确定性,以根据第4.3.3节中得出的阈值评估不同安全状态的概率。5.1. 智能节段设计和现场监测由于地下环境的复杂性和整个生命周期内持续数据采集的需求,监测盾构隧道具有挑战性。为了应对上述挑战,开发了一种基于智能节段的监测系统,该系统包括带有传感器的智能节段、无线网关和云平台,如图14所示。智能节段是基本的监测单元。无线网关负责实现数据的双向传输。通过网关,智能节段上传物理数据,包括温度、旋转角和电池电压。网关还发送下行命令以调整监测周期、睡眠/唤醒状态和其他传感器参数。WIFI设置在无线网关中,安装到隧道后自动连接。云平台可视化监测数据,执行多项任务,包括远程调整系统参数。下载:下载高分辨率图片(735KB)下载:下载全尺寸图片图14. 带有智能节段、无线网关和云平台的监测系统。基于第4.3节建立的传递函数,选择块B1和L1进行传感器嵌入,因为它们对加载变化敏感。如图14(a)所示,在隧道的选定环中组装了带有无线倾斜传感器的智能节段。图14(b)展示了单个由七个块组成的环中两个传感器的嵌入位置,角度分别为271.8°(块L1)和309.6°(块B1),数据采集从节段浇筑开始并持续到整个操作阶段。图14(c)展示了具体的浇筑细节,显示了传感器嵌入情况。为了确保信号传输顺畅,无线网关(图14(d))被安装在嵌入式传感器附近,如图14(e)所示。然后收集的数据被上传到云平台(图14(f)),通过图14(g)中的门户访问,在那里可以选择监测项目并远程设置测量单元和采集参数。图15展示了智能节段的实施过程,涵盖了从节段浇筑(图15(a))、固化(图15(b))和堆放(图15(c)),到盾构隧道施工(图15(d)),直至盾构隧道运行(图15(e))的整个生命周期。本文中的智能节段侧重于数据采集和传输能力,代表了隧道监测技术的基础进步。未来的智能节段将扩展这些能力,包括边缘计算,实现现场数据分析、模式识别以及通过机器学习算法的自主决策。下载:下载高分辨率图片(600KB)下载:下载全尺寸图片图15. 智能节段的生命周期,从制造到运行条件。传感器元件是增强型全方位倾斜传感器节点6306,它能够实现低功耗、高精度(0.002°)的X/Y/Z轴相对水平面的角度测量。支持远程调整监测频率,间隔从一分钟到几小时不等,该设备能够克服与钢筋混凝土结构相关的信号传输挑战。为了确保长期运行,节点在预定测量之间处于睡眠模式,采用低功耗设计。表3总结了无线倾斜传感器和无线网关的详细性能参数。表3. 无线倾斜传感器和无线网关的性能参数(Chang等人,2022年)。增强型全方位倾斜传感器节点6306无线网关参数类型参数类型电源方法内置电池电源方法内置工业级1型电池尺寸200 × 120 × 75 mm尺寸180 × 140 × 60 mm重量3.5 kg重量≤2.0 kg监测范围-90°至+90°本地存储8 GB存储空间,100个传感器测量点的监测周期为一分钟,最小数据存储容量为一年监测精度0.002°无线模块兼容2G/2.5G/3G/4G(Micro SIM卡类型)监测分辨率0.0001°传感器网络协议WiSenMeshWAN?协议5.1.1. 盾构隧道健康监测数据盾构机的尺寸表明,盾构壳可以容纳几乎两个环的节段。组装完成后,智能节段在随后的隧道施工过程中逐渐从盾构尾部推出。到第二环施工结束时,智能节段完全暴露在土壤和水中;此时,它们不再受到盾构尾部刷子的圆周限制。因此,第二环施工完成后记录的数据被用作分析节段旋转的初始参考。块B1的初始倾斜角度分别为X轴-0.998°、Y轴49.598°、Z轴-40.384°,而块L1的初始倾斜角度分别为1.272°、88.645°和-0.467°。图16展示了2024年12月至2025年12月一年内智能节段的监测数据。对于块B1(图16(a)),Y轴和Z轴的倾斜角度呈现类似的趋势,前三个月有逐渐变化,随后趋于稳定。X轴的倾斜角度在整个监测期间相对稳定。对于块L1(图16(b)),三个轴的倾斜角度都保持稳定,仅有轻微波动。下载:下载高分辨率图片(472KB)下载:下载全尺寸图片图16. 智能节段的一年监测数据。5.1.2.从监测数据推断关节开口角度
第4节使用有限元(FE)建模分析了纵向关节开口角度与段体旋转之间的关系,而本节重点关注现场监测数据。为了与第4节中的参数区分开来,从监测数据中得到的参数用上标“m”表示(表示“测量值”):??mBs和??mLs表示B1块和L1块中嵌入传感器位置的旋转;??mBe和??mLe表示B1块和L1块的边缘旋转;??mJ表示B1块和L1块之间推断出的关节开口角度。
基于第4.3节建立的传递函数,可以从??mLs和??mBs推断出??mJ:首先使用方程(6)、(7)将嵌入传感器位置的旋转转换为边缘旋转,然后使用方程(8)计算??mJ。图17展示了整个监测期间??mLs、??mLe、??mBs、??mBe和??mJ的时间历史。逆时针旋转的??mLs和顺时针旋转的??mBs与第4.3节中的数值结果一致,验证了监测策略的可靠性。在整个监测期间,??mLe相对稳定在大约?0.08°。相比之下,??mBe在前三个月内表现出显著变化,从接近零的值变化到大约?0.21°,然后在剩余时间内趋于稳定。??mJ也遵循类似的趋势,在初始调整期后稳定在大约?0.13°。根据方程(9),这相当于大约73.3 kPa的等效表面超载,表明监测的隧道段在预期的运行范围内(0–150 kPa)处于健康状态。前三个月测量到的旋转变化反映了时间依赖的过程,而不是瞬时加载,主要是隧道上方软沉积物的施工后固结以及Liao等人(2017年)报告的环间应力松弛。经过这种调整后,推断出的关节开口角度稳定在大约?0.13°,因此图16和图17构成了全年测量倾斜角度和推断关节状态的时间依赖记录。
下载:下载高分辨率图像(209KB)
下载:下载全尺寸图像
图17. 测量到的段体旋转、边缘旋转和推断出的关节开口角度的时间历史。
5.2. 在负载不确定性下的可靠性分析
虽然第5.1.2节中的确定性推断提供了当前安全状态的实时评估,但它仅基于可用测量值考虑了瞬时条件。为了进行全面的安全评估,必须考虑整个运行期间的负载变化,因为表面超载会因交通负载、临时施工活动和附近的城市发展而变化。为此,第4.2.2节中建立的负载不确定性表征将表面超载S建模为一个截断的对数正态分布,其平均值为???? = 75 kPa,变异系数COV?? = 0.5,上限为Smax = 150 kPa。该分布代表了隧道在其运行期间可能经历的负载条件范围,参数是根据现场条件和工程经验估计的。在此基础上,本节开发了一个可靠性分析框架,首先介绍了用于计算整体安全状态概率的加权积分方法,然后进行了参数敏感性分析(第5.2.2节),以量化这些负载分布参数对可靠性评估结果的影响。
5.2.1. 加权积分可靠性方法
可靠性分析通过对负载分布进行积分来计算每种安全状态的概率。根据第4.3.3节中得出的损坏阈值,采用警告阈值??Warning = 107.46 kPa和临界阈值??Critical = 129.79 kPa来定义三种安全状态:健康状态(概率??Healthy),对应于表面超载低于??Warning;警告状态(概率??Warning),对应于表面超载在??Warning和??Critical之间;以及临界状态(概率??Critical),对应于表面超载超过??Critical。对于给定的负载概率密度函数f(S),每种安全状态的概率由方程(17)、(18)、(19)给出,其中截断的对数正态密度f(S)在方程(3)中定义。
(17) ??Healthy = ∫??Warning????(??)d??
(18) ??Warning = ∫??Critical??Warning???(??)d??
(19) ??Critical = ∫??max??Critical???(??)d??
预期的关节开口角度由以下公式给出:
(20) ???(??J) = ∫??max???J?(??)·???(??)d??
其中??J?(??)是作为表面超载函数的关节开口角度,如方程(9)中所定义。
5.2.2. 参数敏感性分析
可靠性评估结果取决于假设的负载分布参数。为了量化这种依赖性并为不同的运行场景提供指导,通过改变平均值????和变异系数COV??来进行参数敏感性分析。平均表面超载????的影响反映了基线负载条件的差异:较低的值表示隧道位于表面负载控制严格的区域,而较高的值对应于负载更重的城市环境。COV??的影响捕捉了负载条件的变化:较小的值表示稳定、控制良好的负载,而较大的值表示高度可变的条件。这个参数框架使工程师能够:(1) 评估特定地点的负载条件下的可靠性;(2) 评估负载控制措施的有效性;(3) 识别需要加强监测的关键参数组合;以及(4) 支持隧道运营和维护的风险知情决策。
图18展示了????和COV??对??Healthy、??Warning和??Critical概率的综合影响。如图18(a)所示,随着????的增加,??Healthy在所有COV??值下都降低。在低平均表面超载(???? = 15 kPa)时,??Healthy接近1.0,无论变异系数如何,表明处于健康负载范围内的隧道保持了良好的安全裕度。随着????增加到75 kPa,??Healthy在COV?? = 0.5时降低到0.88,在???? = 135 kPa时进一步降低到0.60。与????相比,COV??对??Healthy的影响是次要的,尽管较高的变异性倾向于在中等表面超载水平下降低健康状态的概率。图18(b)展示了警告状态??Warning的概率。在低COV??值时,??Warning与????的关系是非单调的,先增加然后在高表面超载水平时减少。例如,在COV?? = 0.1时,??Warning在???? = 115 kPa时达到约0.63的峰值,然后下降。这种行为反映了警告状态的过渡性质:在非常高的表面超载水平下,概率质量直接转移到临界状态,绕过了警告状态。在更高的COV??值下,这种非单调趋势变得不那么明显。在基准条件(???? = 75 kPa,COV?? = 0.5)下,??Warning约为0.08。图18(c)显示了临界状态??Critical的概率。这种概率随着????的增加而增加,特别是在负载分布集中的低COV??值时。在COV?? = 0.1时,??Critical从???? = 75 kPa时的接近零显著增加到???? = 135 kPa时的0.58。在基准条件(???? = 75 kPa,COV?? = 0.5)下,??Critical约为0.04。较高的COV??值倾向于通过将概率质量分散到所有状态来缓和??Critical的增长。
图18进一步展示了????和COV??对??Healthy、??Warning和??Critical概率的综合影响。
图19进一步展示了????和COV??对状态概率和预期关节开口角度???(??J)的影响。如图19(a)所示,在固定COV?? = 0.5的情况下,??Healthy从???? = 15 kPa时的接近1.0逐渐减少到???? = 75 kPa时的0.88,再到???? = 135 kPa时的0.60。同时,??Warning和??Critical都随着????的增加而增加,在???? = 135 kPa时分别达到0.23和0.17。在整个负载范围内,??Healthy始终占主导地位,尽管在高表面超载水平下显著减少。预期的关节开口角度???(??J)从???? = 15 kPa时的?0.03°增加到???? = 75 kPa时的?0.13°,再到???? = 135 kPa时的?0.19°,反映了在更高负载下的关节逐渐恶化。图19(b)展示了在固定???? = 75 kPa时对COV??的敏感性。概率??Healthy与COV??的关系是非单调的,从COV?? = 0.1时的接近1.0减少到中等变异性(COV??= 0.5)时的大约0.88,然后在更高的COV??值时略有恢复。这种行为反映了对数正态分布的特性:在低COV??时,负载紧密集中在平均值附近,这处于健康范围内;随着COV??的增加,概率质量分布在所有状态之间。预期关节开口角度???(??J)的幅度首先增加,然后随着COV??的增加而减少,在COV?? = 0.25时达到峰值?0.135°,然后在COV?? = 0.9时恢复到?0.107°。
参数敏感性分析为隧道运营和维护提供了实际指导。对于可能在高表面负载(????>100 kPa)下的城市环境中的隧道,建议采取严格的负载控制措施和增加监测频率。对于负载变化较小的隧道(COV??<0.3),从健康状态到临界状态的转变可能突然发生,需要更加保守的安全阈值。本研究中采用的基准条件(????= 75 kPa,COV?? = 0.5)得出??Healthy = 0.88,??Warning = 0.08,??Critical = 0.04,表明隧道在典型运营条件下保持了足够的安全裕度。
6. 讨论
6.1. 所提出框架的优势
本框架的主要贡献是将可测量的量与在隧道运营期间不可观察的安全关键参数联系起来。通过收敛或沉降测量捕捉和监测结构响应的标准技术无法识别可能开始损坏的个别纵向关节的状态。该框架通过段体的旋转生成关节级别的状态,从而将监测方法从全局变形的被动观察转变为局部结构状态的主动推断。
严格的基于物理的校准链确保了这种推断的可靠性。首先,全尺寸弯曲测试表征了实际的力矩-旋转行为,揭示了六个不同的阶段,提供了物理上的真实依据。FE模型根据实验结果进行了验证,随后建立了传递函数(见方程(6)、(7)、(8)、(9)、(10)、(11)、(12)、(13))。使用基于物理的验证方法为设计工程师提供了可解释的关系,而基于物联网的智能段体和可靠性分析框架将这些关系转化为对运行中隧道的连续概率安全评估。
6.2. 局限性和未来工作
所提出的框架将物理传递函数视为确定性的,并且只考虑了负载不确定性。实际上,推断还受到传感器漂移、温度对读数的影响以及简化土壤-结构相互作用的模型不确定性的影响,数值模型没有考虑时间依赖效应,如软土的固结和混凝土的蠕变。因此,从现场数据反计算的等效表面超载73.3 kPa可能会高估实际负载。此外,段级别FE模型没有完全捕捉到全尺寸测试中观察到的急剧的刚度下降,因此当前的传递函数仅适用于运行负载范围,而不是接近失效的条件,需要进一步完善对平台状态的本构建模,以便将框架扩展到最终状态。通过考虑所有不确定性来源以及时间依赖的土壤行为的综合框架可以获得更可靠的估计。
负载-结构模型采用了水平分层的土壤剖面,因此施加的土压力和地面反作用力遵循相当均匀的模式,而在现实中遇到的不均匀土壤条件、渗流下空间变化的水力特性(Chen等人,2026年)和不对称界面(Razgordanisharahi等人,2026年,Moritz,2011年)会在衬砌上产生特征性的正负摩擦力模式,这会改变内部力和传递函数系数。从实施的角度来看,单个仪器化的横截面展示了该框架,但不足以对隧道进行完整评估。传递函数是基于监测横截面的地质条件开发的,需要为具有不同土壤类型、地质条件或负载模式的隧道截面重新生成。此外,激光扫描-based点云重建的补充监测技术可以实现段体几何形状的自动化数字重建,并直接测量关节位移和旋转(Lin等人,2025年),这可以独立验证推断出的关节状态。将这种光学技术与分布式光纤传感技术结合使用,后者可以沿隧道连续测量关节变形(Zhang和Broere,2023年),将是跨验证的多技术监测关节状态的有希望的方向。
在未来的研究中,我们计划将监测系统扩展到隧道沿线覆盖不同地质单元的多个横截面,这可以通过动态更新3D地质模型来支持,以减轻地质不确定性(Jiang等人,2026年)。在那里,我们将测量和模型不确定性纳入可靠性框架,并开发机器学习方法,这些方法明确考虑这些不确定性以更稳健地预测关节状态。一个有希望的途径是通过贝叶斯更新进行连续模型细化。时间依赖效应,如固结和蠕变,也将被纳入环级别模拟中,以便传递函数同时依赖于载荷和时间。此外,由于环级别模型已经为每个超载等级提供了监测接头处的法向力和弯矩,因此可以进一步将这些内部力的代码兼容交互包(Razgordanisharahi等人,2026年)与推断的状态进行对比,从而提供更接近常规施工规范的联合指标。在结构方面,该框架也可以应用于具有增强抗裂性的衬砌材料,例如那些采用混合纤维钢筋系统加固的衬砌材料(Bao等人,2025b年),在这种情况下需要重新推导传递函数,以反映修改后的接头刚度和损伤演化。
7. 结论
本研究提出了一种基于物联网的智能监测框架,该框架能够从盾构隧道中可测量的节段旋转数据中推断出隐藏的纵向接头状态。全尺寸试验、有限元建模、基于物联网的监测和可靠性分析相结合,填补了通常可获得的测量数据与安全关键参数之间的差距,为隧道的预防性维护提供了实际可行的解决方案。本研究的主要贡献如下:
(1)本研究采用的全尺寸弯曲试验(Feng等人,2026年)表明,当纵向接头受到负弯曲时(即外拱展开而内拱闭合),会表现出六个行为阶段。从初始弹性响应到最终破坏的实验性弯矩-旋转关系为模型校准提供了基础。从校准后的有限元模型中得出的传递函数基于物理原理,将传感器测量值与接头开口角度和混凝土损伤联系起来。这保证了推断状态的一致性和可解释性。这些函数能够表征随着表面荷载变化而演变的损伤过程,从而客观评估预警和临界阈值,减少了对主观经验判断的依赖。
(2)在一座正在运营的地铁隧道中进行的一年现场监测验证了所提出框架的实际适用性。带有嵌入式无线倾斜传感器的智能节段能够在节段浇筑期间持续采集数据,无需在接头位置安装侵入式传感器。推断出的接头开口角度在初始变化后稳定在约-0.13°,相当于大约73.3千帕的表面荷载,表明隧道处于预期的正常运行范围内。需要注意的是,这个等效值可能会高估实际的表面荷载,因为当前模型没有考虑时间依赖的土壤固结效应。
(3)开发了一个可靠性分析框架,用于量化在荷载不确定性下的这些安全状态的概率,这些不确定性由截断的对数正态分布表示。在平均表面荷载为75千帕、荷载变异系数为0.5的典型运行条件下,健康状态、预警状态和临界状态的概率分别约为88%、8%和4%。参数敏感性分析显示,平均表面荷载水平对安全状态的概率具有主导作用,而变异系数对健康状态的概率的影响是非单调的。
**作者贡献声明:**
冯浩蓝:撰写——原始草稿、可视化、验证、软件开发、资金筹集、正式分析、数据管理。
张东明:撰写——审稿与编辑、监督、资源协调、项目管理、资金筹集、概念构思。
斯特凡诺·马里亚尼(Stefano Mariani):撰写——审稿与编辑、监督。
加布里埃莱·德拉韦基亚(Gabriele Della Vecchia):撰写——审稿与编辑、监督。
黄宏伟(Hongwei Huang):监督、资源协调、项目管理、概念构思。
陈刚(Gang Chen):资源协调、概念构思。
生物通微信公众号
生物通新浪微博
今日动态 |
人才市场 |
新技术专栏 |
中国科学人 |
云展台 |
BioHot |
云讲堂直播 |
会展中心 |
特价专栏 |
技术快讯 |
免费试用
版权所有 生物通
Copyright© eBiotrade.com, All Rights Reserved
联系信箱:
粤ICP备09063491号