神经网络反演研究法在基坑支护变形数值模拟与预测中应用

1、神经网络反演研究法在基坑支护变形数值模拟与预测中应用【摘要】：本文将神经网络与数值模拟相结合，用 神经网络反演分析被支护体“综合力学性能参数”，建立 flac 3d模型，通过神经网络表达从基坑变形特征值到土体 物性参数的非线性映射关系，创建神经网络结构，并将其应 用于数值模拟中来实现对基坑支护变形的预测。【关键词】：神经网络、反演分析、flac 3d、数值模拟abstract： in this article, authors combine neural network with numerical simulation, back analyse the "combinetion

2、 property parameter” of supported soil by neural network, create flac 3d model, express the nonlinear relationship from deformation character of foundation pit to property parameter of soil, establish neural network structure, and use it in deformation numerical simulation of foundation pit in order

3、 to realize forecast for deformation .keywords： neural network, back analysis, flac 3d, numerical simulation中图分类号：tn711文献标识码：a文章编号：基坑的变形对于开挖所引起的环境效应评价以及信息 化施工都有着十分重要的意义。然而，在目前的各种数值计 算方法中，很难对基坑的变形做出较为准确的模拟和预测。 其原因就在于对岩土体的初始参数难于确定，这也是目前数 值计算方法应用于岩土体模拟中遇到的最大困难之一【1, 2o初始参数不易确定的原因是多方面的，主要有：岩土体 在空间分布上的复杂性

4、与多变性，室内土工试验不能准确测 定某些参数或不能全面反映岩土体力学性能，初始地应力的 影响，地下水的影响，土体本构模型的不准确性以及数值模 拟中因模型的简化与实际情况的差异等。鉴于此，本文在对 基坑变形进行数值模拟时，尝试引入神经网络反演分析法， 以消除上述因素的影响，同时也希望借此方法打开初始参数 不易确定这一阻碍数值模拟在岩土工程中发展的"瓶颈” 问题。1、神经网络反演分析应用于数值模拟中的思路神经网络几乎可模拟任何复杂的非线性系统，特别适用 于参数变量和目标函数之间无数学表达式的复杂工程问题， 一个简单的神经网络模型可以反映一个非常复杂的映射关 系【3, 4】。因此，用神经网

5、络模型来表达岩土工程中岩土 体物性参数与岩土体位移值之间的映射关系是有理论可行 性的。神经网络反演分析的思路：首先建立较为符合工程实际 的数值模拟分析模型，通过不断调整土体物性参数，由数值 模型计算分别得到与之对应的基坑变形特征值。然后将上述 依依对应的土体物性参数和变形特征值作为训练神经网络 的样本，以基坑变形特征值作为输入向量，以土体物性参数 作为目标输出向量，对建立的神经网络模型进行学习、记忆 和训练，利用神经网络的黑箱处理方法，建立起输入向量和 目标输出向量之间的内在关系。最后，利用施工阶段实际监 测到的变形数据，输入训练好的网络，得到能够反映基坑变 形特征的土体“综合力学性能参数”，

6、并以此作为正演数值 模拟分析程序的输入参数，来实现对下一步开挖基坑的变形 以及基坑最终变形的预测。该方法的优点在于：通过反演分析得到的土体“综合力 学性能参数”包含了由于岩土体分布不均匀性、以及初始地 应力、地下水影响以及数值模拟模型简化所引起的误差修正 等综合信息，以此进行正演模拟分析，可以得到较为贴近工 程实际的分析结果。2、工程实例验证本文以某广场深基坑支护工程为例，运用上述方法进行 数值模拟分析，以此来验证该方法的应用可行性。2.1、工程背景某基坑支付工程，基坑东西向长308. 2m,南北向宽 46. 5m82. 5m,总开挖面积近20000平方米。基坑开挖深度 东西北三侧为17.62

7、1.0m,基坑南侧开挖深度为13.81& 7mo场地地层分布：人工填土层(qml)；第四系坡洪积 层(qdl+pl)；第四系残积层(qel)；燕山期基岩(丫53)； 基坑南侧主要为粒质粉质粘土。基坑南侧采用预应力锚索复合土钉支护，其他三面采用 桩锚式复合土钉支护。场地地层分布主要为粒质粉质粘土， 少量粉土和燕山期花岗岩。2.2、建立flac 3d模型【5】由于flac在岩土工程数值模拟中的众多优势，因此本 文以美国itasca公司开发的flac 3d,即三维快速拉格郎日 差分分析程序对某广场基坑南侧的预应力锚索复合土钉支 护进行数值模拟分析。在建立flac 3d模型时有以下几个问 题值

8、得注意：(1) 合理选择结构单元。尽管flac 3d程序中已经建立 的梁(beam)、索(cable)、壳(shell)以及桩(pile)结构单元, 但使用时要根据实际情况注意区分，尤其注意各结构的连接 以及接触面的相互作用。本文建立模型时，对于土钉和预应 力锚索采用索结构单元，对面层喷射混凝土采用壳结构单 元，腰梁采用梁结构单元。(2) 土体本构模型的选取。目前土体的本构模型已有200 多种，其中能较好地模拟实际土体的受力变形性能的本构模 型大多比较复杂，参数太多，很难在实际工程分析中加以应 用；较简单的模型虽然考虑因素不是很全面，但能够说明最 主要的问题，因此，要根据数值模拟的主要目的，选

9、择既简 单又能反映主要问题的模型。本文侧重于基坑变形的分析， 因此选择了考虑土体的剪胀性和静水压力影响的drucker-prager 模型。(3) 边界条件的确定。duncan和goodman认为，边坡的 计算范围可取边坡坡脚以下h深、坡脚左右水平2h3h的 范围，其中h为边坡高度。另外，在模型的两个侧面，认为 土体不发生y方向的位移，其应力边界条件按静止土压力计 算。结合上述注意事项的考虑，本文建立flac 3d数值模拟 分析模型如图1所示。图1、flac 3d分析模型fig 1 analysis model of flac 3d2.3、生产训练样本由建立好的flac 3d模型，通过输入不断

10、改变的土体物 性参数，得到与之对应的基坑边坡位移和沉降值，从而产生 多组由土体物性参数与位移和沉降值依依对应的多组训练 样本。土体物性参数主要包括土体弹性体积模量k、弹性剪切模量g、粘聚力c和内摩擦角，其他参数基本都可由c、值 导出。目前室内土工试验对土体的粘聚力c和内摩擦角测量 的比较准确，因此本文采用c、的室内土工试验值，不将其 作为神经网络的输出向量。而土体的弹性体积模量k、弹性 剪切模量g通常由土体的弹性模量求得，对此室内土工试验 不易把握，同时k、g值又是影响基坑边坡土体变形的两个 主要因素，因此，本文把k、g值作为神经网络的目标输出 向量，把通过flac3d计算得到的基坑变形特征值

11、，即位移 (x-disp)和沉降(z-disp)作为神经网络的输入向量，这样就 建立了 x-disp、z-disp与k、g对应的样本集合。据此，建立了 23组样本，如下表所示。表1、神经网络训练样本tabl training samples of neural network2.4、创建神经网络结构如果说flac 3d程序计算过程表明了从土体物性参数到 基坑变形特征值的基于有限差分理论的映射关系，那么，通 过神经网络所表达的则是从基坑变形特征值到土体物性参 数的高度非线性的映射关系。鉴于此映射关系的高度复杂 性，本文采用两层前馈bp网络结构，第一层选用8个神经 元和双曲正切s形激励函数，第二层

12、选用2个神经元和正线 性激励函数，所建立的bp网络模型结构如图2所示。图2、bp网络结构模型fig 2 strueture model of bp neural network2.5、训练网络、输出结果对bp网络的训练可以由mat lab程序实现，主要可分为 两个阶段【3】：第一阶段是输入已知学习样本，通过设置的 网络结构和前一次迭代的权值和阈值，从网络的第一层向后 计算各神经元的输出；第二阶段是对权值和阈值进行修正， 从最后一层向前计算各权值和阈值对总误差的影响（剃度）， 据此对各权值和阈值进行修改。以上两个过程反复交替，直 到达到收敛为止。当误差性能达到预先设定的标准后，即认 为对网络的训

13、练已经完成。此时网络的输出结果已经十分接 近期望输出目标值，达到可以接受的范围，这一点可由误差 性能曲线体现（图3）。图3、误差性能曲线fig 3 error performance curve从图3看到，总误差已经比较接近目标值0.01,说明输 入与输出之间的映射关系已经建立起来。下面就以基坑开挖 完第八层后s22监测点的数据（x-disp=36. l,y-disp=20. 6）, 输入训练好的网络，得出土体的“综合力学性能参数”k = 6. 6847e6pa, g=3. 9487e6pao2.6、模拟结果分析以神经网络反演出的k、g值输入flac 3d正演分析程 序来实现对第九步开挖和最后

14、一步开挖（第十步）基坑支护 变形的预测，并与监测点s22和s23的实测值比较，其结果 如下。表2、模拟值与实测值对比表tab2 comparison of between simulation resuits and actual measurement resuits表2中可以看出，模拟值与实测值差别不是很大，吻合 程度较高，基本说明应用神经网络反演分析土体“综合力学 性能参数”，再用于数值模拟对后续施工过程中基坑支护变 形的预测是切实可行的。神经网络反演分析法溶入基坑支护变形数值模拟中，不 仅仅能够对某些点的水平位移和沉降做出较为准确地预测， 而且还能够较为准确地反映被支护边坡的总体变形趋

15、势以 及产生最大变形的位置，如图4、图5所示。图4、x方向位移等值线(水平位移)fig 4 contour of x displacement(horizontai displacement)图5、z方向位移等值线(沉降)fig 5 contour of z displacement(settlement)从上面的两幅图中可以看到：预应力锚索复合土钉支护 结构的水平位移最大值发生在基坑边坡中下部；沉降最大值 并不是在靠近基坑开挖边线处，而是在距开挖边线一定距离 的位置；同时基坑底靠近坡脚的位置发生较大的隆起变形。上述特点都与实际观测到的总体变形趋势是一直的。为了进一步了解对被支护边坡总体变形趋

16、势预测的准 确程度，通过flac3d对不同深度土体的水平位移进行预测, 并与实际测斜仪的测量结果比较，发现吻合程度相当好（见 图6）。这说明预测结果有较高的可信程度，可以应用于工程 实际以指导现场施工和信息化施工。fig 6 comparison between simulation resuits and actual measurement results of soil along depth为了对比溶入神经网络反演分析法的数值模拟与传统 数值模拟二者的差异，本文还将由室内土工试验得出的k= 4.917e6pa、g=2. 472e6pa（k、g值均通过土体弹性模量由公 式求得，具体过程见

17、文献6）,输入flac3d分析程序，并保 证在相同的模型下进行计算，结果如下。图7、s22、s23监测点水平位移模拟曲线fig 7 horizontal displacement simulation curve of s22 and s23图8、s22、s23监测点沉降模拟曲线fig 8 settlement simulation curve of s22 and s23图7、图8表示了基坑开挖完成后两个监测点水平位移 和沉降随时间步长的变化，从中可以看出，达到稳定后的水 平位移和沉降值明显偏大于表2中最后一步开挖完成后基坑 变形的实际监测结果，其详细对比结果见表3。表3、传统模拟值与实测值

18、对比表tab3 comparison of between traditional simulationresults and actual measurement results由表3知，传统数值模拟方法所分析的结果与实测值出入较大，误差均在30%左右，用此结果对基坑变形的预测意 义不大。这也从另一方面说明了神经网络反演分析法在数值 模拟中所发挥的积极作用，它克服了因土体初始参数不易确 定所带来的影响，在实际工程应用中有着重要的现实意义。3、结论通过实例验证，得到如下结论：（1）将神经网络反演分析法应用于基坑支护变形的预测 中，具有理论可行性和应用可行性。（2）神经网络反演分析法可以解决数值模拟中，岩土体 初始参数不易确定的难题，但需要建立较为合理的神经网络 结构和足够的训练样本对网络进行训练。（3）溶入了神经网络反演分析法的数值模拟，能够对后 续施工过程中基坑的变形、基坑最终变形，以及基坑变形的 总体趋势进行较好的模拟和预测，有利于信息化施工。4、深化研究方向受时间关系本文未对神经