一种用于深厚软土地铁基坑土体参数反演分析的方法
【专利摘要】本发明涉及岩隧道施工及土工程参数计算【技术领域】，具体涉及一种用于深厚软土地铁基坑土体参数反演分析的方法，是将有限元数值计算结果与实际监测值之间差值的最小二乘函数乘以一个特定的加权系数作为目标函数，应用改进的高斯-牛顿方法对目标函数进行最优化求解，通过反复迭代计算使得有限元数值计算结果逐渐逼近真实的监测值，从而确定待定参数的最优解，建立起了监测值与基坑土体参数基本特征值之间的关系。本发明结合实际观测数值不断去反演与修正土体参数进而不断地预测下一阶段土体变形，预测结果随着监测值的数量增加而更加精确，为深厚软土地铁基坑工程的施工提供非常好的监测与预警方法，在土体参数确定的过程中具有较高的应用价值。
【专利说明】一种用于深厚软土地铁基坑土体参数反演分析的方法
[【技术领域】]
[0001]本发明涉及岩隧道施工及土工程参数计算【技术领域】，具体涉及一种用于深厚软土地铁基坑土体参数反演分析的方法。
[【背景技术】]
[0002]基坑工程作为一项岩土工程，其特点是很强的系统性和不确定性，若基坑开挖不当，会引发基坑工程自身的安全和稳定问题，对周边建筑物、地下管网等造成影响，特别地，当地质环境为深厚软土时，潜在的破坏风险更高。因此能否准确预测基坑在施工过程中的变形，对基坑工程的安全施工具有重要指导意义。有限元法由于其强大的适用性，正日益成为深基坑工程中预测基坑变形的强有力工具。有限元分析是利用数学近似的方法对真实物理系统进行模拟，利用简单而又相互作用的元素，就可以用有限数量的未知量去逼近无限未知量的真实系统。
[0003]现有的有限元分析时的土体材料本构采用修改剑桥本构模型(Modified Cam-ClayModel)，该本构从理论上与试验上都较好的阐明了土体的弹塑性变形特性，是应用最为广泛的软土本构模型之一。而对于本构模型参数的确定，目前有三种方法:理论方法、试验方法与参数反演法。理论方法依据力学理论与土体的性质进行理论求解，但计算过程中作了一定的假设，且其假设往往与实际工程不相符合，故理论解析无法正确地表现地基土的力学参数。试验方法一般是由室内土工实验结果选取或是由现场原位测试确定，由于土样受到扰动、试验方法的不完善及原位测试地点受人为主观因素的影响，使得试验结果无法全面地反映地基土体的性质。参数反演分析方法是根据基坑工程施工现场实测的数据(如位移、应力、孔隙压力)运用反演理论来计算、分析地基土体的力学参数，使力学参数更加接近土体的实际力学性质，以作为工程设计的依据而进行正分析变形预测，这种方法逐渐成为当前研究的热点。但参数反演分析的相关研究表明，其优化目标函数是一个高度复杂的非线性多峰值函数，目前，优化的方法主要是BP神经网络反分析、单纯形法反分析等优化反分析的方法，但这些方法要么目标函数容易陷入局部极小值，要么算法当复杂，优化参数过多而鲁棒性不强。
[
【发明内容】
]
[0004]本发明为了解决上述问题，设计了一种利用监测的位移值来反演得到基坑土体的基本参数特征值，将有限元数值计算结果与实际监测值之间差值的最小二乘函数乘以一个特定的加权系数作为目标函数，应用改进的高斯-牛顿方法对目标函数进行最优化求解，通过反复迭代计算使得有限元数值计算结果逐渐逼近真实的监测值，从而确定待定参数的最优解，建立起了监测值与基坑土体参数基本特征值之间的关系的用于分析深厚软土地铁基坑土体参数的反演分析方法。
[0005]为了实现上述目的，提供一种用于深厚软土地铁基坑土体参数反演分析的方法，该方法包括以下步骤:
[0006](I).采用ABAQUS大型通用有限元分析软件建立合理的计算模型，根据深厚软土的地质特性，土体的本构模型采用修正剑桥模型Modified Cam-Clay并且输入模型计算所需的已知参数与待反演参数的初始值；
[0007](2).编制土体参数反演分析新方法的核心计算程序并嵌入到ABAQUS有限元软件中，该反演分析新方法的核心算法如下:
[0008]a.确定最小二乘加权目标函数:S(b) = [y-y' (b)]Tco[y-y，(b) ] = eTcoe，其中，b为待反演参数的一个向量,该向量包含所有待反演的参数，y为实际监测值的向量，y’ (b)为与监测值相对应的有限元计算值，ω为权重系数矩阵，e为残差向量；
[0009]b.计算最小二乘加权目标函数是一个非线性回归迭代计算过程，采用修正高斯-牛顿法进行计算，该过程的数学描述为:
=CiX:'⑴(_>’-;/(/7,.))与\+1 = Prdr+br，其中，dr 是一个通过更新初始的反演输入参数值来优化b向量的一个向量，下标r是预计迭代的次数，&是一个用来评估br敏感性的系数矩阵，&中每个值根据计算得到，C是一个对角缩放矩阵，该矩阵中每个元素&，I为单位矩阵，πν是一个用来改进非线性回归的一个调整参数，P ^是一个迭代衰减参数，特别地，对于残差较大且非线性程度很高的问题，当迭代计算超过3次之后最小二乘加权目标函数变化值任然小于0.01时将X,7
替换为及来重新迭代计算:
[0010]C.敏感性系数矩阵Xr采用有限差分法进行计算，根据差分计算方式的不同，有:向前差分格式
【权利要求】
1.一种用于深厚软土地铁基坑土体参数反演分析的方法，其特征在于该方法包括以下步骤: (1).采用ABAQUS大型通用有限元分析软件建立合理的计算模型，根据深厚软土的地质特性，土体的本构模型采用修正剑桥模型Modified Cam-Clay并且输入模型计算所需的已知参数与待反演参数的初始值； (2).编制土体参数反演分析新方法的核心计算程序并嵌入到ABAQUS有限元软件中，该反演分析新方法的核心算法如下: a.确定最小二乘加权目标函数:S(b)= [y-y' (b)]Tco[y-y’ (b) ] =eTcoe，其中，b为待反演参数的一个向量，该向量包含所有待反演的参数，y为实际监测值的向量，I，(b)为与监测值相对应的有限元计算值，ω为权重系数矩阵，e为残差向量； b.计算最小二乘加权目标函数是一个非线性回归迭代计算过程，采用修正高斯-牛顿法进行计算,该过程的数学描述为:
与=其中，4是一个通过更新初始的反演输入参数值来优化b向量的一个向量，下标r是预计迭代的次数，Xr是一个用来评估br敏感性的系数矩阵，Xr中每个值根据Xii=SyiAbi计算得到，C是一个对角缩放矩阵，该矩阵中每个元素
为单位矩阵，πν是一个用来改进非线性回归的一个调整参数，P ^是一个迭代衰减参数，特别地，对于残差较大且非线性程度很高的问题，当迭代计算超过3次之后最小二乘加权目标函数变化值任然小于0.0i时将替换为ωΧν + R未重新迭代计算: C.敏感性系数矩阵I采用有限差分法进行计算，根据差分计算方式的
不同，有:向前差分格式
以及中心差分格式
两种格式，其中，I' i是第i次的有限元计算结果；bj是待反演参数向量中第j个参数；Abj是bj的一个扰动量；八2用来表示该差分格式为中心差分格式； d.为了确保反演分析算法的合理性与正确性，还需对优化算法进行拟合统计分析，从三个方面入手，综合考虑模型拟合统计:1)采用误差方差
来评估总体的加权残值，其中，S(b)是最小二乘加权目标函数，ND是监测值的数量，NP是待反演参数的数量；ii)每次迭代完成后反演效果采用FI = [S(b) W,-S(b) ttft]/S(b) 来评价，其中，s(b) 为初始的最小二乘加权目标函数值，S(b) 为优化后最小二乘加权目标函数值；iii)采用相关系数
来评估整体参数反演分析算法； e.实际监测值权重系数ω系数矩阵是一个对角矩阵，该矩阵的每个元素ω = l/σ2，其中，σ i为bi的标准偏差值，bi为第i个待反演的参数值；(3).联合参数反演算法程序与ABAQUS有限元程序进行迭代计算，若最小二乘加权目标函数达到计算收敛精度要求，则停止计算，此时反演所得参数即为所需的土体参数值，并将反演所得参数重新代入有限元模型中进行计算以预测深厚软土地铁基坑开挖过程中的土体变形。
2.如权利要求1所述的一种用于深厚软土地铁基坑土体参数反演分析的方法，其特征在于迭代衰减系数P r的取值范围在O~1.0之间。
3.如权利要求1所述的一种用于深厚软土地铁基坑土体参数反演分析的方法，其特征在于标准偏差值采用Oi = 2Ei/l.96进行估算，其中EiS仪器设置自身误差。
4.如权利要求1所述的一种用于深厚软土地铁基坑土体参数反演分析的方法，其特征在于待反演的参数值在优化过程中被限定在一定的区间之内，该区间采用数学表达式描述为:
其中，Yo含义为y1-b i空间中，bi取值变化曲线与Ji轴的交点值。
【文档编号】G06F17/50GK104077451SQ201410314277
【公开日】2014年10月1日 申请日期:2014年7月3日 优先权日:2014年7月3日
【发明者】贾春雷, 刘习生, 穆保岗, 竺明星, 曾旭, 钱琨 申请人:中铁上海工程局集团有限公司, 中铁上海工程局集团华海工程有限公司