高等岩石力学-岩石力学有限元法概念.ppt

高 等 岩 石 力 学;岩石力学有限元法; 第一节 概 论;有限元法基本方程 就数学概念来说，有限元法是通过变分原理(或加权余量法)和分区插值的离散化处理把基本支配方程转化为线性代数方程，把求解待解域内的连续场函数转化为求解有限个离散点（节点）处的场函数值。显然，这种离散化的处理是一种近似，因面中有当单元划分得充分小时，才能保证满意的求解精度。 由于单元充分小，在一个微小的单元内，未知场函数可以采用十分简单的代数多项式近似地表述。通常可取为如下的插值形式： (4-1);下标m表示单元的节点数目，下标e表示单元的序号。 有限元法即是以所有节点处的1值作为基本未知量。对于二维的工程物理问题，如热传导、渗流等，其基本支配方程为如下形式的准谐方程： (4-2) 边界条件为： ① ，边界 在上 ② 边界在 上 ;; 对于一个单元体的子域内，把方程式（4-1） 带入方程式 （4-3)并进行变分后即可得到单元的基本方程： (4-5) 式中 即单元特性矩阵： 即为节点未知量组成的矢量： 是由在 上的边界条件给定的a（x,y）有关的矩阵，故仅对a≠0的二类边界上的单元才予考虑。对整个求解域A由变分驻值条件式（4-4）可得 或写为 (4-6) ; 式中，[k]为由各单元特性矩阵 及 按节点号组集得到的总体特性矩阵：{u}为所有节点的待求值组成的矢量：{p}为与Q（x,y）有关的矢量及q(x,y)有关的矢量。 解线性方程式（4-6）即可求得场函数 在各单元节点处的值。 对于弹性力学问题，可通过最小势能原理或虚功原理导出有限元法的基本方程。有限元法求解弹性力学问题通常以位移作为基本未知量。位移在直角坐标系中沿坐标轴X, Y, Z 的分量分别表示为 ，故单元位移模式这时可写为 (4-7) 式中，l为3×3阶单位矩阵；Nt 称为单元的位移插值函数或形函数（对二维问题l为2×2阶，位移分量 ）；{ }为由单元各节点位移分量形成的矢量称为移矢量， 利用弹性力学的几何方程及物理方程可导出单元的应变及应力表达式： ; （4-8） （4-9） 应用最小势能原理或虚功原理可以推导出单元刚度矩阵的表达式 （4-10，a） （4-10，b） （4-10，c）; 式中，{q}为分布体力