岩体工程中的反分析方法概述(PPT).pptVIP

令 则 b 体积力势能 设体积力 则体积力势能为： c 表面力势能 设面力 则面力势能为： d 集中力势能 设集中力 则集中力势能为： e 单元总势能 将所有单元势能叠加得系统总势能 根据最小势能原理，真实解应使系统势能最小 即 由此得系统总平衡方程 引入位移边界条件，得最终求解方程： 5.3 线弹性位移反分析 5.3.1 反分析基本公式 以隧道开挖平面应变问题为例。设开挖边界上的初始地应力为： 将该初始地应力转化开挖边界上的等效节点力： 或写成： 在整个求解域上： 根据有限元求解基本方程： 令： 则： 假定量测点与有限元网格节点重合，则可把节点位移分成已知和待求量部分： 相应的平衡方程写为： 将未知位移消去： 其中： 上式可简写为： 其中： 因为测量都是两点相对位移，绝对位移与相对位移之间转化关系： [ T ] 为转换矩阵： 则： 其中 上式中待求量 为3个，若量测值 刚好为3个，则可从上式中求出唯一的 若量测值 多于3个，则通过最小二乘法求得 ，构造以下目标函数 令 得 的最小二乘估计为： 5.3.2 实例 根据围岩内部位移求得： 根据围岩收敛变形求得： 隧道埋深400m，可求得： 则： 5.4 粘弹性位移反分析 5.4.1 反分析基本公式 实际测试的位移是一组随时间变化的值，是由于释放荷载随时间的改变及围岩的蠕变造成，采用粘弹性反分析可以很好反映时间因素。 任意时刻 t 的应变可以看做由瞬时弹性应变和蠕变应变两部分构成： 其中弹性应变： 粘弹性应变： 则： 得： 其中： 为弹性矩阵 任意时刻 t 的有限元平衡方程： 将 代入上述平衡方程： 其中： 称之为综合模量。 常用流变模型的粘弹性模量和综合模量 粘弹性模型 粘弹性模量 综合模量 Maxwell Kelvin Kelvin-Voigt Poynting-Thomson Burgers 将位移 U( t ) 分为已知和未知两部分 仿照线弹性位移反分析公式推导，得粘弹性反分析基本方程 其中： 5.4.2 参数回归与优化 实际量测过程中，设置测点时间为 t0，到 t 时刻量测的位移是 t0~t 之间发生的位移，为此把综合模量也用t0~t 之间的变化来代替 设置测 点时间t0 时间t 设从 t0~t 之间共进行了n次量 ，每次量测有m个测点，量测结果为： 则利用每次量测的m个测点已知位移，可反求该时刻的综合模量 ，n 次量测可以反求出不同时刻的综合模量： 根据综合模量与模型参数之间的关系，比如P-T模型： 式中： 将上式变换为以下形式： 两边取对数： 令： 得： 用矩阵表示： 其中： 最小二乘法解为： 根据求得的a、b、yi 再求E0、E1、E2、 设 时刻的量测值为： 5.5 工程应用 5.5.1 量测数据预处理 量测数据带有误差，甚至有奇异数据，为此需对其进行一定的数学处理，改善其规律性。最常用的就是多项式拟合。 令： 建立如下目标函数： 求解方程组 得 ai 从而得到： 利用回归公式求得任意时刻的u(ti)进行反分析。 5.5.2 反分析结果检验 方法一：将反分析得到的参数作为输入，在相同模型下进行计算，通过计算结果与量测结果的比较来衡量反分析结果的正确性。 方法二：将反分析得到的参数作为输入，对尚未发生的位移进行预测，根据预测的精度来衡量反分析结果的正确性。 5.5.3 测点布置 （a） 3条收敛测线 （b） 6条收敛测线 （c）收敛测线加内部位移 5.5.4 围岩变形预测 利用小规模实验工程量测结果进行反分析，再用反分析结果对实际工