边坡稳定分析边坡工程的反分析法.docx

边坡稳定分析边坡工程的反分析法 §2.9.1反分析法简述 由确定的力学模型，给定岩土介质物理参数以及受力状态而求边坡岩土体的位移、应变、应力等，称之为正问题，我们在前面各节中介绍的方法可以求解此类问题。通过工程现场的位移等量测信息作为“计算参数”，推求边坡岩土体的初始地应力以及介质参数，这类分析称之为反分析，或反问题。 反分析的基本思想最先由Kavangh(1973)、Gioda和Maier(1980)等人提出，Sakurai(1983)首次给出了均匀地应力与岩体弹性模量的有限元反分析数值解。近年来，国内许多学者已做了大量的研究工作。刘维宁(1993)对岩土工程反分析方法的信息问题进行了系统的研究，杨林德(1996)、吕爱钟和蒋斌松(1998)对岩土工程和岩石力学的反演理论和分析方法进行了系统的研究和阐述。李宁、崔政权等(1996)对边坡安全监测进行了仿真反演分析(2000)，对碎裂块体围岩安全监测与仿真进行了反演分析。谭昌明(1999)提出了动态反演予测分析法和迭代反演综合预测法。郑颖人院士和高玮(2001)提出了用快速遗传算法进行岩土工程反分析的方法。刘迎曦等(2001)对边坡的地层参数提出了优化反演的有限元方法。 反分析法作为解决岩土工程介质本构模型及物性参数选求问题的有效方法，其分析方法主要有优化法和逆解法。反演参数求解应满足适定性，即解的存在性、稳定性。 §2.9.2反分析基本原理 (1)初始应力场 边坡开挖后，边坡的位移与开挖形状、边坡岩土力学状况、地质结构、初始地应力场等诸多因素有关。对于一个给定的边坡结构，边坡岩土中某点的位移总可表示为 （2.145） 式中：{u}——空间点(x，y，z)处t时刻的位移； x、y、z——边坡中测得位移的空间坐标点； t——观测位移的对应时刻； E、μ、c、、——分别为边坡岩土体介质的变形参数和强度参数； ——边坡岩土体的容重； T——代表温度因素； k——代表渗流因素； {σ0}——初始应力场； {ST}——代表其他地质构造因素。 若初始地应力场是一非均匀分布场，则{σ0}可表示为 （2.146） 式中：[X]——空间坐标的函数矩阵； {P}——初始地应力分布参数组成的矢量； {S0}——随已知因素变化而变化的量(如自重变化引起的应力)。 对平面问题，假定初始地应力场为位置坐标(x，y)的线性函数，则有： 式中：P1，…，P9——待定参数。 式(2.147)中的9个参数Pl，P2，…，P9应满足以应力表示的变形协调方程和平衡微分方程，即 （2.148） （2.149） 在Y轴向上的右手螺旋坐标中，如只计入自重引起的体力，则有 按式(2.147)分布的初始应力场将满足变形协调方程式(2.148)，由此可得 （2.150） 则可将式(2.147)改写为 与式(2.146)对应，就有 （2.152） （2.153） （2.154） 可证明在满足式(2.148)、式(2.149)后，式(2.153)中的7个参数P1，P2，P3，P4，P5，P6，P7，是相互独立的。 在式(2.145)中，如除{σ0}为未知数外，其他因素及参数都已知，则可由观测所得的位移u，根据有限元列式反推求得表征{σ0}的初始地应力场分布参数。 (2)初始应力场反推基本方程 在有限元计算中，有限元的平衡方程为 引入条件后，有 （2.155） 这里，参数{R}是初始地应力场经开挖后释放所产生的荷载，即 （2.156） （2.157） （2.158） 设任一点i在αi方向(αi为与x轴正向夹角)的位移ui为已知，并设i点在单元e内，则ui可表示为 （2.159） 式中：u、v——i点在坐标轴x、y方向位移矢量； {Tαi}T——位移转换矩阵，{Tαi}T=[cosαi，sinαi]。 由有限元列式的几何关系，有 （2.160） 式中：——单元e在测点i处的形函数矩阵； {δe}——单元e的位移矢量。 将式(2.160)代入式(2.159)，有 （2.161） 令 （2.162）（2.163） 则 （2.164） 如果有n个测点，则可将各测点位移ui(i=1，2，…，n)代入式(2.164)，从而得到n个方程组。可用矩阵表示为 （2.165） 式中：m——{P}的维数。 而有 （2.166） 上式即为用有限元列式建立的初始地应力场分布参数的绝对位移反推关系式，或反推基本方程。 §2.9.3参数反演的优化法 (1)优化法定义 参数反演求解要满足适定性。所谓适定性，是指除物理基本定律所确定的关系外，尚需知道多少有关解的局部信息，才足以在全局范围