边坡稳定分析边坡大变形的流形元法.docx

边坡稳定分析边坡大变形的流形元法 §2.7.1数值流形方法的位移 流形方法的局部位移函数是建立在单个物理覆盖上的。在整个材料域内，流形单元上所有物理覆盖的位移函数的加权求和形成该单元的总体位移函数。对于有限流形方法，物理覆盖上的位移函数与材料边界无关。覆盖上的位移函数可以取为常数、线性或高阶多项等。在平面问题中，物理覆盖上的位移函数可取为 流形单元e是由q个覆盖的交集构成，则该单元的总体位移函数为 式中：m——是物理覆盖i上所定义的位移函数的项数； uij、vij——覆盖i上第j个沿x向和y向的位移变量； Sij(x，y)——覆盖i上的基本函数； We(i)——覆盖i上位移函数对单元e贡献权值。 将式(2.102)化简写为 其中， 式中，权函数We(i)表示组成该流形单元e间的各物理覆盖的位移函数对流形单元位移函数的贡献比率。 对于基于有限数值流形方法，因数学网络是有限网络，对于均质材料由各物理覆盖交集构成的流形单元正是有限单元，此时各物理覆盖的位移函数就相当于流形单元e的“结点”位移。流形单元e由几个物理覆盖组成j就有几个“结点”。故将各覆盖对流形单元e的贡献权函数可直接取为该单元的形函数。对于含有节理、裂隙等不连续面的岩体，流形单元的权函数不同于有限元的形函数。一般情况，在每一单元中，权函数应满足 并且 式中：Vi——构成流形单元e的物理覆盖。 §2.7.2大变形流形方法的基本方程 (1)应力应变分析 流形方程的应力应变分析与有限元法类似，即本构关系为 式中：{σ}——应力矩阵； [D]——弹性矩阵； {ε}——应变矩阵。 考虑大变形以初始构形为参考构形，格林应变矩阵可表示为 式中第一项为应变分量的线性部分，即 其中， 第二项为应变分量的非线性部分，即 对于平面问题 而 由式(2.110)、式(2.111)、式(2.114)及式(2.116)有 （2.119） 式中， (2)大变形分析的平衡方程 数值流形方法的平衡方程的建立可借助于最小势能原理。设分析域内有N个流形单元，n个物理覆盖，每个覆盖有2m个待求量，则整个系统的总势能为 （2.120） 对于覆盖i （1.121） 于是，可导出岩石大变形分析的流形方法整体平衡方程 （1.122） 式中：Di——覆盖i内待求位移变量{di1，di2，…，di2m}T； Fi——覆盖i内分布到2m个位移变量上的等效荷载{Fi1，Fi2，…，Fi2m}T； Kij——切线刚度系数子矩阵，其为2m×2m阶矩阵，且具有；。 将式(2.122)简写为 （2.123） 式中：{F}——总体荷载矩阵； {U}——总体位移变量列阵，即 （2.124） {KT}——总体切线刚度矩阵，对于流形单元e，切线刚度矩阵可直接表示为 （2.125） 式中： （2.126） （2.127） （2.128） （2.129） 以上公式中的积分均是在物理覆盖的交集(即流形单元e)上进行的。