ARTICLE · 1249110
重磅!用AI来操作PLAXIS3d 建模并计算以及出报告
本来有些事情是不想承认的，但是也不得不承认。接着之前那一篇写报告的，这一个搞计算。附带我跟ai对话过程。
顶管下穿既有铁路桥梁
PLAXIS 3D 数值模拟计算报告
—— 大港电厂至南港工业区蒸汽管网顶管下穿南港铁路桥梁工程
（依据安全评估报告，采用 PLAXIS 3D 有限元软件远程脚本自动化建模与计算）
大港电厂至南港工业区蒸汽管网工程需下穿既有南港铁路桥梁，穿越段采用顶管法施工。为了评估顶管施工及运营对既有铁路桥梁（墩台）的影响，采用岩土工程有限元软件 PLAXIS 3D（2024 版）建立三维数值模型，通过远程脚本接口（Python + plxscripting）完成从建模、网格划分、分阶段计算到结果提取的全流程自动化，据此分析墩台附加位移与顶管上方地表沉降，为施工与监测控制提供依据。
依据《安全评估报告》确定的关键设计参数如下：
•顶管：内径 1.60m，外径 1.98m，壁厚 0.19m，管中心高程 -3.14m；穿越地层主要为淤泥质黏土（⑥1 层），顶进方向沿 X 轴，建模长度 80m（x∈[-40, +40]m）。
•桥梁：双线 32m 简支 T 梁铁路桥，133#/134# 墩中心距 32m（y=±16m），承台底高程分别为 -7.53m / -7.11m，墩顶（梁底）高程 +13.5m。
•列车荷载：双线中-活载换算均布 60kPa，作用于桥面板。
二、三维数值模型
2.1 模型范围与地层
模型平面范围 100m×80m（x∈[-50,50]m，y∈[-40,40]m），底面高程 -38.25m（深度 40m）。按安评报告钻孔资料自上而下建立 5 层地层，土体本构均采用摩尔-库仑（Mohr-Coulomb）模型：
层号 | 土层名称 | 顶底高程 (m) | 重度γ (kN/m³) | 弹性模量 E (kPa) | c (kPa) | φ (°) |
① | 填土 | 1.75 ~ -0.71 | 17 / 18 | 8000 | 12 | 8 |
⑥1 | 淤泥质黏土 | -0.71 ~ -13.5 | 16.5 / 17.5 | 3500 | 15 | 5 |
⑥2 | 粉质黏土 | -13.5 ~ -15.5 | 17.5 / 18.5 | 7000 | 20 | 10 |
⑧1 | 粉质黏土 | -15.5 ~ -19.0 | 18 / 19 | 14000 | 25 | 14 |
⑧2 | 粉砂 | -19.0 ~ -38.25 | 19 / 20 | 22000 | 2 | 32 |
2.2 结构材料
构件 | 材料 | 力学参数 | 模拟方式 |
桥墩 / 承台 | C30 混凝土 | 线弹性、非多孔，γ=25kN/m³，E=3.0×10⁷kPa，ν=0.2 | 实体单元（Volume） |
顶管管壁 | C50 混凝土 | 弹性板，厚 0.19m，γ=25kN/m³，E=3.45×10⁷kPa | 板单元（Plate，圆柱面） |
桥面板 | C30 混凝土 | 弹性板，厚 1.5m，E=3.0×10⁷kPa | 板单元（Plate） |
2.3 几何建模要点
（1）顶管管壁采用板单元模拟：沿管轴线在半径 0.895m 处建立直线，绕管中心轴线旋转 360° 生成圆柱面并转换为 plate 板单元，管壁直接嵌入土层（不建立开挖实体），分为东段（0~+40m）与西段（-40~0m）两段，对应顶进方向自东向西。
（2）桥梁墩台采用实体单元：133#/134# 承台与墩身为挤压成型的实体（Volume），桥面板为板单元；承台底深入土体、墩身高出地面。
（3）列车荷载：在桥面板上施加 60kPa 均布面荷载（双线中-活载换算）。
图1Phase_4 运营阶段三维整体竖向位移云图（PLAXIS Output 直接导出，管-土-桥整体模型）
图2沿管轴纵向剖切（y=0 竖直面）竖向位移云图（PLAXIS Output 直接导出，显示管-土相对位置与沉降槽）
2.4 网格划分
采用 1.5m 目标单元尺寸生成网格，全部由六面体主导单元组成，网格质量检查（checkgeometry）通过，单元数满足三维精细化计算要求。模型网格与变形显示如下图所示（PLAXIS 模型直接导出）。
图3初始阶段模型网格与变形显示（PLAXIS Output 直接导出）
2.5 施工阶段（分阶段施工）
阶段 | 名称 | 阶段内容 |
0 | 初始阶段 | 仅地层，地应力平衡 |
1 | Phase_1 既有铁路现状 | 激活桥梁墩台（显式赋予 C30 材料）、桥面板与列车荷载；位移清零（作为基准） |
2 | Phase_2 顶管顶进·东段 | 东段管壁板（0~+40m）就位 |
3 | Phase_3 顶管顶进·西段 | 西段管壁板（-40~0m）就位 |
4 | Phase_4 运营阶段 | 顶管全部完成，沉降稳定 |
各施工阶段的模型状态（Phase_0~3 为竖向位移云图，Phase_4 为总位移云图，均为 PLAXIS Output 直接导出）如下：
图4初始阶段（地应力平衡，仅地层）
图5Phase_1 既有铁路现状（激活桥梁与列车荷载，位移清零）
图6Phase_2 顶管顶进·东段（0~+40m）
图7Phase_3 顶管顶进·西段（-40~0m）
图8Phase_4 运营阶段（顶管全部完成，沉降稳定）
三、建模与计算中的关键技术问题
本次建模与计算通过 PLAXIS 3D 远程脚本接口全流程自动化完成，过程中发现并解决了一系列技术难点：
1.必须先创建钻孔再添加土层：PLAXIS 3D 中必须先创建钻孔（borehole）才能添加土层，否则触发参数错误。
2.材料赋值采用对象属性方式：直接使用 setmaterial 命令会导致参数错误，改为对象属性赋值（obj.Material = mat）可正常完成土层与管壁材料的绑定。
3.管壁板单元生成：polycurve 命令不接受点列表，采用 line + revolve(360°) 生成圆柱面后转换为 plate 板单元。
4.网格生成后对象代理失效：结构模式创建的对象在进入网格/阶段模式后 GUID 失效，必须在 gotostages 后从集合重新获取阶段对象。
5.板单元按创建索引重取：g_i.Plates 集合顺序稳定，按创建时记录的索引取回对应管壁与桥面板。
6.承台被网格吸收为土体：完全埋入土体的承台在网格化后被 PLAXIS 归类为 BoreholeVolume 土壤体，名称与索引均不可靠；改用创建时原始代理的 get_equivalent() 找回全部 4 个桥墩体，并在激活阶段用 setmaterial 显式赋予 C30 混凝土材料，确保阶段模型中桥墩按混凝土计算（这是本次桥墩材料问题修复的关键）。
四、计算结果
以 Phase_1（既有铁路现状，位移清零基准）为基准，提取 Phase_4（运营阶段）的相对增量位移，评价顶管施工对既有桥梁与地表的影响。
4.1 墩台顶部位移
测点 | Uz(Phase_1) mm | Uz(Phase_4) mm | Uz 增量 mm | Utot(Phase_1) mm | Utot(Phase_4) mm | Utot 增量 mm |
133# 墩顶 | -87.26 | -87.64 | -0.38 | 87.28 | 87.71 | +0.43 |
134# 墩顶 | -90.93 | -91.30 | -0.37 | 90.96 | 91.48 | +0.52 |
顶管施工引起的 133#、134# 墩顶附加竖向位移仅约 0.38mm、0.37mm，附加总位移约 0.4~0.5mm，对既有桥梁的影响很小，远低于一般桥梁墩台沉降控制值。
4.2 顶管上方地表沉降
x (m) | 沉降增量 (mm) | x (m) | 沉降增量 (mm) | x (m) | 沉降增量 (mm) | x (m) | 沉降增量 (mm) |
-40 | -2.97 | -15 | -5.58 | 10 | -6.18 | 35 | -4.77 |
-35 | -4.78 | -10 | -5.29 | 15 | -6.10 | 40 | -3.16 |
-30 | -5.38 | -5 | -5.22 | 20 | -5.78 | ||
-25 | -5.68 | 0 | -6.24 | 25 | -5.65 | ||
-20 | -5.76 | 5 | -6.56 | 30 | -5.33 |
顶管施工引起管轴线上方（y=0）地表附加沉降最大约 6.6mm（位于 x≈0~5m 区段），小于一般地表沉降控制值（20~30mm）；沉降影响范围沿管轴约 ±30m。横断面（x=0）沉降槽底约 6.2mm，位于管正上方 y=0 处，向两侧衰减，影响范围约 ±30m。
图9顶管施工引起的地表沉降曲线（管轴线上方，y=0，计算结果数据）
图10横断面沉降槽曲线（x=0，计算结果数据）
图11地表沉降盆平面分布（顶管施工增量，mm，计算结果数据）
图12纵向剖面（y=0）沉降增量等值线（计算结果数据，虚线框为管道位置）
4.3 位移云图（PLAXIS Output 导出）
以下云图由 PLAXIS Output 直接渲染导出。需要说明：累计竖向位移（Uz）包含既有桥梁自重与列车荷载引起的沉降，最大约 93mm 位于桥墩区域；而顶管施工引起的附加沉降（增量）仅约 6~7mm，见 4.1~4.2 节。
图13Phase_4 变形网格 |u| 云图（PLAXIS Output 导出，累计，放大 50 倍）
图14横向剖切（x=0，垂直管轴）Uz 云图（PLAXIS Output 导出，沉降槽横断面）
图15水平剖切（地表附近）Uz 云图（PLAXIS Output 导出，沉降盆平面分布）
五、结论与建议
7.顶管施工引起的既有桥梁墩顶附加竖向位移约 0.4mm（133# 墩 0.38mm、134# 墩 0.37mm），附加总位移约 0.4~0.5mm，影响甚微，既有铁路桥梁结构安全有保障。
8.顶管施工引起顶管轴线上方地表附加沉降最大约 6.6mm，位于管正上方 x≈0~5m 区段；横断面沉降槽底约 6.2mm，影响范围约 ±30m，均处于常规控制标准之内。
9.板单元模拟管壁、分阶段施工（东段→西段）的计算方法可有效模拟顶管下穿过程；网格化后桥墩材料通过 get_equivalent() + setmaterial 显式赋值得到正确传递。
10.建议施工期对 133#、134# 墩及顶管轴线上方地表进行自动化监测，重点控制管正上方地表沉降，并结合注浆压力、顶进速度等参数动态调整施工。
六、附：本次建模计算的技术迭代历程
本次工作从零开始，通过 PLAXIS 3D 远程脚本逐步打通“建模—网格—分阶段计算—结果提取—成图—报告”的全流程，主要历程如下：
11.先跑通最小示例：实现“新建项目→创建土体轮廓→钻孔与土层→材料→赋值→网格→计算”全流程，解决“先建钻孔再建土层”与“setmaterial 命令参数无效改用对象属性赋值”两个问题。
12.按安评报告参数建立完整三维模型：不断迭代调整模型结构，先后删除圆柱形桩基、移除顶管开挖区实体、仅以板单元模拟管壁，修复网格生成失败与集合索引重排导致的对象引用错误。
13.修复桥墩材料问题：发现网格化后埋入土体的承台被 PLAXIS 吸收为土壤体（BoreholeVolume），且属性赋值无法可靠传递到阶段模型；改用创建时原始代理 get_equivalent() 找回全部桥墩体，并在阶段中对每个桥墩显式调用 setmaterial 赋 C30 混凝土材料。
14.完成全部分阶段计算（初始阶段→既有铁路现状→顶管东段→顶管西段→运营阶段）。
15.输出结果提取：解决 Output 远程脚本端口占用问题（10001 被无关服务占用，改用 10002），并完成 Output 2024 新 API 的探测（getresults + 结果集转数组 + 最近节点匹配），提取墩台位移与地表沉降数据。
16.成果可视化：导出 PLAXIS Output 渲染云图与剖切图，绘制沉降曲线、沉降盆等值线与管-土相对位置剖面图，形成本报告。
—— 完 ——
