基于能量突变的滑动面搜索方法.docVIP

基于能量突变的滑动面搜索方法 张伟 潘建平 江西理工大学建筑与测绘工程学院 X 关注成功！ 加关注后您将方便地在 我的关注中得到本文献的被引频次变化的通知！ 新浪微博 腾讯微博 人人网 开心网 豆瓣网 网易微博 摘????要： 安全系数和滑动面位置能够准确分析边坡的稳定性, 虽然有限元强度折减法可以得到边坡的安全系数, 但是对滑动面位置的确定无法达成共识。我们基于能量理论, 分析了边坡能量组成, 边坡破坏时的能量变化, 从单元体的弹性应变能入手, 发现处在滑动面处的单元体弹性应变能发生突变, 我们以坐标定位的形式划分一个个单元体, 比较相邻单元体弹性应变能是否发生突变, 通过一系列离散点来刻画滑动面。将本文方法和Spencer法进行比较发现误差很小, 证明该方法是可行的。实际操作过程中还发现, 单元体体积大小、离散点密度都影响滑动面的搜索。 关键词： 边坡; 稳定性; 弹性应变能; 突变; 滑动面; 边坡的稳定性一直是岩土力学的研究方向, 如何获得准确的安全系数和滑动面位置学者们进行了大量研究。近年来随着有限元分析的长足发展, 有限元强度折减法为我们获得安全系数提供了有效的方法。但是如何确定滑动面位置目前还不成熟。Cheng[1]等指出最大剪应变和最大剪应变云图都可以描绘滑动面, 且两者给出的滑动面位置极为接近。Lane[2]等发现边坡破坏后的网格会出现一个畸变带, 认为这条畸变带就是滑动面位置。岩土体破坏是由于某一面上的剪应力达到了岩土体的剪切强度, 此时剪切面上必然发生剪切变形, 据此李剑、陈善雄[3]等提出了根据最大剪应变增量的边坡潜在滑动面搜索方法。根据边坡破坏的特征可以知道, 边坡破坏时滑动面上节点位移和塑形应变将产生突变, 滑动面就在水平位移和塑形应变突变的地方。依据滑动面位置附近位移等值线最为密集, 袁维、胡叶江[4]等提出了基于位移场的滑动面搜索方法。近年来对岩土体能量研究的深入, 涂义亮[5]等认为边坡破坏时伴随着能量变化, 边坡释放能增加, 弹性变形能先增大后减小有个突变的过程。赵阳生[6]等认为岩土体弹性变形与集聚能量是一个稳态过程, 岩土体破坏与能量释放是一个失稳过程, 能量转换遵循岩土体破坏的最小能量理论。郭建强、刘新荣[7]等提出了基于弹性应变能的岩石强度准则, 引入弹性应变能释放分散系数, 建立了基于弹性应变能的强度准则。通过对比岩石破坏时的能量转换和释放, 同样在边坡问题中也存在这种现象。边坡滑动面处的土体能量释放显著, 弹性应形能变化显著。 1 边坡破坏过程中的能量变化 根据能量守恒原理: 式中:U为外部环境对边坡输入的总能量;U为边坡耗散能;U为岩石的弹性应变能。 其中为研究方便U按下式计算: 边坡中的应力应变并非均匀分布, 所以边坡单元体储存的弹性应变能都不一样。当边坡破坏时, 滑动面周围的土体应力分布开始调整, 应力状态迅速转变为二向, 最终转变为单向应力状态, 发生能量跃迁, 主要就是弹性应变能的改变。当破坏完成时能量释放主要发生在滑动面处, 滑动面周围土体发生塑形变形, 弹性应变能得到释放, 一部分转化为边坡的动能, 在滑动面处材料弹性应变能相较于周围不处滑动面上的材料会变小。滑动面必然出现在某个单元体中, 该单元体的弹性应变能释放, 就相当于原子的裂变释放能量。在我们依据有限元强度折减法中的材料强度参数, 利用FLAC3D内嵌的FISH语言编制单元体弹性应变能计算程序。 2 建立一种边坡分析模型 边坡单元体在同一纵向上储存的弹性应变能应当是连续性的, 在滑动面处单元体能量发生突变。为了衡量这种突变性质, 我可以用在纵向上相邻两个单元体的能量差, 差值最大的地方 (发生能量突变) 就是滑动面所在的位置, 我们取两个单元体坐标的中点为滑动面所经过的点。边坡单元体弹性应变能进行同一纵向比较, 记坐标单元体弹性应变能为Ui, 坐标单元体弹性应变能为Ui+1。可以得到相邻单元体能量差值ΔU。 图1 ??下载原图 如果某个纵向有i个单元体就会得到 (i-1) 个差值ΔU, 比较这些ΔU其中取得最大时, 相邻两个坐标的中点就是我们所求滑动面经过的点。等间距的在横向上布置一系列纵向单元体, 这样每个纵向方向上都有与滑动面相交的的离散点, 连结这些离散点就可以得到滑动面所处的位置。 3 滑动面曲线拟合 由于在有限元分析中边坡单元体中的所有应力、应变数据都是以单元体中心点的数据, 这样两个单元体中心点的中点所得的离散点组成的滑动面不一定是光滑面, 需要进行光滑处理。我们运用最小二乘法对离散点所组成的曲线进行拟合。最小二乘法就是所有离散点到拟合函数曲线的距离达到最小或者满足某个精度要求。我们将拟合函数写成向量的形式: 选取一组线性无关的函数组定义: 离散点到拟合曲线距离 依据最