岩石力学反分析课件

岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论11绪论岩石力学中的反分析最早由Kavangh（1973）、Gioda和Maier（1980）等人提出，Sakurai（1983）完成了岩体弹性模量和初始地应力的线弹性有限元位移反分析，此后又发展了弹塑性、粘弹性、粘塑性等非线性位移反分析，并引入了误差分析、优化技术等一些手段，以求获得非线性反分析中的最佳值。1.1发展历史1绪论岩石力学中的反分析最早由21.2正分析与反分析1正分析2反分析①②1.2正分析与反分析1正分析2反分析①②3岩石力学中的反分析主要有以下几种类型：已知岩体的本构模型、初始地应力和位移量测值，求岩体物理力学参数;已知岩体的本构模型、物理力学参数和位移量测值，求初始地应力;已知岩体的本构模型、物理力学参数、初始地应力和位移量测值，求开挖空间最佳几何形状；已知初始地应力和位移量测值，求岩体的本构模型及模型参数，即系统辨识。3反分析问题的特点多解性、无解性、不稳定性。岩石力学中的反分析主要有以下几种类型：已知岩体的本构模型、初41.3反分析中的几个要素1模型模型是“原型”的一种“类似”，任何模型都不能反映出原型的一切特征。模型的表达形式可以是概念的、物理的或数学的，用数学描述形式建立的模型为数学模型。2参数和状态参数是系统的内部状态变量，反映了系统的本质，是不可测量的；状态是系统的外部表现，是可以测量的。1.3反分析中的几个要素1模型5在岩石力学数学模型中，因变量，如位移、应力、应变均为外部状态变量，弹性模量、泊松比、内粘结力等均为参数。3准则函数由于模型的近似性和量测误差的存在，在已知量和待求量之间对等的情况下，求出的结果往往不能很好地反映系统的本质。可行的方法就是增加已知量的数量，求待求量的最优值，为此需要引入一个准则函数。准则函数有两类：以量测值为基础的第一类准则函数；以量测误差及其统计特性为基础的第二类准则函数。常用准则函数。在岩石力学数学模型中，因变量，如位移、应力、应变6①常规最小二乘法②高斯——马尔可夫估计③最大似然估计④贝叶斯估计①常规最小二乘法②高斯——马尔可夫估计③最大似然估计④71.4反分析求解方法1逆法将模型输出表达成待求量的显函数，与量测值构成准则函数直接求解。2正法当模型输出不能表达成待求量的显函数时，先给出待求量的初值，计算出模型的输出，与量测值一起代入准则函数求出准则函数值，按一定的路径取待求量的值，可计算出一系列准则函数值，使得准则函数值达到最小的待求量值即为最优值。该方法是由一系列正算过程构成，故名正法。其适用范围较逆法更广。1.4反分析求解方法1逆法将模型输出表达成待8正法中要用到最优化方法，最常用的有模式搜索法、变量轮换法、单纯形法、鲍威尔法。正法中要用到最优化方法，最常用的有模式搜索法、变量9岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论102.1有限单元法的基本思路将连续求解域离散为有限个、按一定方式相互连接在一起的单元组合体，在每个单元内用一假设的位移函数来表示待求的未知位移场函数，而假设的位移函数用单元节点上的未知位移来表示，以此可导出单元内以未知节点位移所表示的应力、应变，最后通过最小势能原2有限元法正分析简要2.1有限单元法的基本思路将连续求解域离散11理导出单元每个节点上以未知节点位移表示的平衡方程，整个求解域所有节点平衡方程将构成一方程组，通过求解该方程组可求得各节点上的位移，从而求得单元内的位移、应力、应变。总之，有限单元法就是将连续无限自由度问题转变为求离散的有限自由度问题，将偏微分方程组的求解转化为代数方程组的求解。理导出单元每个节点上以未知节点位移表示的平衡方程，整个求解域12对求解连续体的边值问题，有限单元法是一种近似方法，近似的程度随单元划分加密而提高，但会带来计算工作量的增加。2.2有限单元法求解的一般过程2.2.1计算模型的建立确定计算域，以工程建设影响不到的地方作为计算模型的边界。划分单元，靠近开挖区域单元化分密些，远离开挖区域单元化分疏些。单元边长尽可能相近。以平面四节点的等参单元为例。对求解连续体的边值问题，有限单元法是一种近似方13岩石力学反分析课件142.2.2计算方法第一步计算初始地应力场；第二步将释放荷载转化成等效节点力施加在开挖边界上，计算相应的应力、位移；第三步将前两步应力叠加即为开挖引起的应力，第二步的位移即为开挖引起的位移。2.2.3平面有限元法基本公式①位移模式2.2.2计算方法第一步计算初始地应力场；2.15其中：②几何方程其中：其中：②几何方程其中：16（i=1、2、3、4）根据等参单元的坐标变换式：（i=1、2、3、4）根据等参单元的坐标变换式：17得：其中：得：其中：18③本构方程其中：③本构方程其中：19④单元势能分析a单元应变能④单元势能分析a单元应变能20令则b体积力势能设体积力则体积力势能为：令则b体积力势能设体积力则体积力势能为：21c表面力势能设面力则面力势能为：c表面力势能设面力则面力势能为：22d集中力势能设集中力则集中力势能为：d集中力势能设集中力则集中力势能为：23e单元总势能将所有单元势能叠加得系统总势能根据最小势能原理，真实解应使系统势能最小e单元总势能将所有单元势能叠加得系统总势能根据最小势能原24即由此得系统总平衡方程引入位移边界条件，得最终求解方程：即由此得系统总平衡方程引入位移边界条件，得最终求25岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论263线弹性位移反分析3.1反分析基本公式以隧道开挖平面应变问题为例。设开挖边界上的初始地应力为：将该初始地应力转化开挖边界上的等效节点力：或写成：3线弹性位移反分析3.1反分析基本公式以隧道27在整个求解域上：根据有限元求解基本方程：令：则：假定量测点与有限元网格节点重合，则可把节点位移分成已知和待求量部分：在整个求解域上：根据有限元求解基本方程：28相应的平衡方程写为：将未知位移消去：其中：相应的平衡方程写为：将未知位移消去：29上式可简写为：其中：因为测量都是两点相对位移，绝对位移与相对位移之间转化关系：上式可简写为：其中：因为测量都30岩石力学反分析课件31岩石力学反分析课件32[T]为转换矩阵：则：其中上式中待求量为3个，若量测值刚好为3个，则可从上式中求出唯一的。若大于3个，则的最小二乘估计：[T]为转换矩阵：则：333.2实例3.2实例34根据围岩内部位移求得：根据围岩收敛变形求得：隧道埋深400m，可求得：则：根据围岩内部位移求得：根据围岩收35岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论361绪论岩石力学中的反分析最早由Kavangh（1973）、Gioda和Maier（1980）等人提出，Sakurai（1983）完成了岩体弹性模量和初始地应力的线弹性有限元位移反分析，此后又发展了弹塑性、粘弹性、粘塑性等非线性位移反分析，并引入了误差分析、优化技术等一些手段，以求获得非线性反分析中的最佳值。1.1发展历史1绪论岩石力学中的反分析最早由371.2正分析与反分析1正分析2反分析①②1.2正分析与反分析1正分析2反分析①②38岩石力学中的反分析主要有以下几种类型：已知岩体的本构模型、初始地应力和位移量测值，求岩体物理力学参数;已知岩体的本构模型、物理力学参数和位移量测值，求初始地应力;已知岩体的本构模型、物理力学参数、初始地应力和位移量测值，求开挖空间最佳几何形状；已知初始地应力和位移量测值，求岩体的本构模型及模型参数，即系统辨识。3反分析问题的特点多解性、无解性、不稳定性。岩石力学中的反分析主要有以下几种类型：已知岩体的本构模型、初391.3反分析中的几个要素1模型模型是“原型”的一种“类似”，任何模型都不能反映出原型的一切特征。模型的表达形式可以是概念的、物理的或数学的，用数学描述形式建立的模型为数学模型。2参数和状态参数是系统的内部状态变量，反映了系统的本质，是不可测量的；状态是系统的外部表现，是可以测量的。1.3反分析中的几个要素1模型40在岩石力学数学模型中，因变量，如位移、应力、应变均为外部状态变量，弹性模量、泊松比、内粘结力等均为参数。3准则函数由于模型的近似性和量测误差的存在，在已知量和待求量之间对等的情况下，求出的结果往往不能很好地反映系统的本质。可行的方法就是增加已知量的数量，求待求量的最优值，为此需要引入一个准则函数。准则函数有两类：以量测值为基础的第一类准则函数；以量测误差及其统计特性为基础的第二类准则函数。常用准则函数。在岩石力学数学模型中，因变量，如位移、应力、应变41①常规最小二乘法②高斯——马尔可夫估计③最大似然估计④贝叶斯估计①常规最小二乘法②高斯——马尔可夫估计③最大似然估计④421.4反分析求解方法1逆法将模型输出表达成待求量的显函数，与量测值构成准则函数直接求解。2正法当模型输出不能表达成待求量的显函数时，先给出待求量的初值，计算出模型的输出，与量测值一起代入准则函数求出准则函数值，按一定的路径取待求量的值，可计算出一系列准则函数值，使得准则函数值达到最小的待求量值即为最优值。该方法是由一系列正算过程构成，故名正法。其适用范围较逆法更广。1.4反分析求解方法1逆法将模型输出表达成待43正法中要用到最优化方法，最常用的有模式搜索法、变量轮换法、单纯形法、鲍威尔法。正法中要用到最优化方法，最常用的有模式搜索法、变量44岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论452.1有限单元法的基本思路将连续求解域离散为有限个、按一定方式相互连接在一起的单元组合体，在每个单元内用一假设的位移函数来表示待求的未知位移场函数，而假设的位移函数用单元节点上的未知位移来表示，以此可导出单元内以未知节点位移所表示的应力、应变，最后通过最小势能原2有限元法正分析简要2.1有限单元法的基本思路将连续求解域离散46理导出单元每个节点上以未知节点位移表示的平衡方程，整个求解域所有节点平衡方程将构成一方程组，通过求解该方程组可求得各节点上的位移，从而求得单元内的位移、应力、应变。总之，有限单元法就是将连续无限自由度问题转变为求离散的有限自由度问题，将偏微分方程组的求解转化为代数方程组的求解。理导出单元每个节点上以未知节点位移表示的平衡方程，整个求解域47对求解连续体的边值问题，有限单元法是一种近似方法，近似的程度随单元划分加密而提高，但会带来计算工作量的增加。2.2有限单元法求解的一般过程2.2.1计算模型的建立确定计算域，以工程建设影响不到的地方作为计算模型的边界。划分单元，靠近开挖区域单元化分密些，远离开挖区域单元化分疏些。单元边长尽可能相近。以平面四节点的等参单元为例。对求解连续体的边值问题，有限单元法是一种近似方48岩石力学反分析课件492.2.2计算方法第一步计算初始地应力场；第二步将释放荷载转化成等效节点力施加在开挖边界上，计算相应的应力、位移；第三步将前两步应力叠加即为开挖引起的应力，第二步的位移即为开挖引起的位移。2.2.3平面有限元法基本公式①位移模式2.2.2计算方法第一步计算初始地应力场；2.50其中：②几何方程其中：其中：②几何方程其中：51（i=1、2、3、4）根据等参单元的坐标变换式：（i=1、2、3、4）根据等参单元的坐标变换式：52得：其中：得：其中：53③本构方程其中：③本构方程其中：54④单元势能分析a单元应变能④单元势能分析a单元应变能55令则b体积力势能设体积力则体积力势能为：令则b体积力势能设体积力则体积力势能为：56c表面力势能设面力则面力势能为：c表面力势能设面力则面力势能为：57d集中力势能设集中力则集中力势能为：d集中力势能设集中力则集中力势能为：58e单元总势能将所有单元势能叠加得系统总势能根据最小势能原理，真实解应使系统势能最小e单元总势能将所有单元势能叠加得系统总势能根据最小势能原59即由此得系统总平衡方程引入位移边界条件，得最终求解方程：即由此得系统总平衡方程引入位移边界条件，得最终求60岩石力学反分析1绪论2有限元法正分析简要3线弹性位移反分析岩石力学反分析1绪论613线弹性