第五章变形驱动力反演.pptVIP

第五章 变形驱动力反演 §5.1 引言 §5.2 弹性力学基本方程 §5.3 弹性力学平面问题的有限单元法 §5.4 基于弹性力学平面问题的变形驱动力反演 §5.5 基于弹性力学空间问题的变形驱动力反演 §5.1 引言 1970年开始的国际地球动力学10年计划将探讨构造运动的力源作为该计划的重点目标。这从一个侧面反映了对变形体所受到的外力来源和分布的探讨的重要性。 变形驱动力的反演在构造应力场研究中已经做了许多工作，1983年王仁提出了构造应力场反演的叠加法，许多学者纷纷以该方法或以该方法为基础进行驱动力的变化信息，力源变化信息无疑对变形解释和预报有着不可替代的作用。 §5.1 引言 叠加法思想 将区域边界划分为若干段，首先计算对各段施加单位力时在区域内产生的应力场或位移场，然后将这些与单位边界力对应的场乘以待定系数作线性叠加而得出区域内的应力场。再利用他们去拟合区域内应力或位移的观测结果，并根据最佳拟合准则求解出一串待定系数，而所求边界作用的大小和方向即由这些系数乘以单位力确定。 本章提出一种与上述方法不同的方法，他同样是基于有限元方法，根据大地测量测定的位移，可同时得到变形体的外界力和底部驱动力。 §5.2 弹性力学基本方程 定义 弹性力学基本方程描述了弹性物体物理力学性质、几何关系、应力应变关系、力的平衡关系和边界条件。它们是求解弹性力学问题的基础。 基本公式 设（x y z）T为弹性体中任意一点的坐标，其位移、应变、应力分别表示为： 位移： 应变： 应力： §5.2 弹性力学基本方程 任一点的位移与应变关系描述 弹性体内任一点处的位移与应变之间的关系由如下几何方程描述： 其矩阵形式为： 此处L为微分算子矩阵。 §5.2 弹性力学基本方程 任一点的应力-应变关系描述 弹性体内任一点处的应力与应变之间的关系由如下物理方程描述： 这里D称为刚性矩阵： 作用在弹性体内任意一点处的体积力f=（fx fy fz）T与该点处的应力满足平衡方程： 其中L是几何方程中定义的微分算子矩阵。 §5.2 弹性力学基本方程 边界条件 弹性问题定解的求取需要给出边界条件。边界条件可以以力的边界条件或位移边界条件的形式给出。力的边界条件是给出作用在边界上的外力，位移边界条件是给出边界上的位移。对于一个弹性体，边界S可以分为力的边界Sσ和Sμ位移边界。有S=Sσ+Sμ。 原理上，弹性力学问题可有上述基本方程和边界条件得到解答。但是由于实际问题所涉及的弹性体比较复杂，解析解一般难以得到，只能通过数值方法求解。有限元方法是解这种力学问题最常用的一种数值方法。