基于FLAC3D改进的初始地应力场回归方法

摘要：

为提高地应力反演的精度，提出了FLAC3D中利用偏最小二乘回归方法拟合地应力场及将地应力场精确加载至计算模型并实现平衡的方法，首先对各应力分量分别进行回归，以提高局部异常应力区域的拟合精度。然后依据FLAC3D的计算原理，分析了单元应力与节点不平衡力的关系，将回归所得的单元初始地应力加载至模型并计算节点不平衡力，在节点施加与不平衡力反向的节点荷载，实现力的平衡，完成地应力场的加载。当采用弹塑性模型时，假定平均应力不变，根据屈服条件修正偏应力大小，对不满足强度条件的单元地应力进行修正，以满足弹性假设条件，修正后的单元地应力满足平衡条件及强度条件，并能精确满足地应力场回归估计值，接着介绍了一种将模型单元地应力张量表示成各单元坐标的函数，通过求取插值函数获得任意坐标的地应力张量的三维超曲面样条插值方法。最后以工程实例证实，改进方法可以有效提高局部地应力异常区域的拟合精度，解决了边界应力奇异分布的问题。

关键词：地应力反演、偏最小二乘回归、数值计算、FLAC3D

分类号：TU459;P642.3(土力学、地基基础工程)

资助基金：国家自然科学基金项目41072234,51274157；国家自然科学青年基金项目51109123；道路桥梁与结构工程湖北省重点实验室开放基金项目DQJJ201304；山东省优秀中青年科学家科研奖励基金项目BS2012NJ006

在线出版日期：2014-06-19（万方平台首次上网日期，不代表论文的发表时间）

页数：共7页

页码：892-898