岩土体渗流自由面问题的重心插值无网格方法     

李树忱  王兆清  袁超 《岩土力学》2013,34(7):1867-1873

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。  相似文献   

含圆柱形弹性夹杂地基沉降量的解析解     

王明斌  李术才  李树忱  王兆清 《岩土力学》2008,29(6):1575-1579

研究各向同性地基中含一个圆柱形弹性夹杂的地基沉陷问题。考虑全场的应力、变形特征,运用复变函数的罗朗级数展开技术,构造合适的应力函数,结合夹杂的边界条件和无穷远的性态,获得了全场复势的解析表达。数值结果显示：由于两相材料剪切模量之比的不同,在某个范围内沉陷量大于基底中点沉陷;在这个范围之外,则小于基底中点沉陷。当夹杂埋深为夹杂半径的5倍以上时,就可以运用所得结果,而不致引起太大误差;含圆形隧洞和刚性夹杂的问题可作为特殊情况得到。  相似文献   

圆形隧道围岩附加应力场的解析解答     

王明斌  李术才  李树忱  王兆清 《岩土力学》2006,27(Z1):207-210

研究各向同性岩体中含一个圆形隧道围岩的附近应力问题。考虑全场的应力变形特征,运用复变函数的罗朗级数展开技术,构造合适的应力函数,结合隧道的边界条件和无穷远的性态,获得了全场应力复势的解析表达式。数值分析发现,应力集中程度最剧烈时,不一定总是发生在孔的附近;应力集中程度最弱时,不一定总是距离孔边较远;当隧道埋深为其半径的5倍以上时,可以作为无限大基体含一个圆形孔洞来处理。  相似文献   

多边形有限单元形函数的几何构造方法     

王兆清  李术才  李树忱  王明斌 《岩土力学》2006,27(Z1):65-68

根据非均质材料的细观结构,划分多边形的计算网格。采用多边形单元进行有限元分析,实现了基于材料真实结构的数值模拟。给出了多边形有限单元形函数的几何构造方法,构造了一个辅助多边形,采用散度定理推导出多边形单元形函数的表达式。分析总结了多边形Wachspress插值、Laplace插值和平均值插值的构造方法和性质。  相似文献   

无压渗流问题分析的多节点有限元方法   总被引：2，自引：1，他引：1  

王兆清  李术才  李树忱 《岩土力学》2008,29(10):2647-2650

采用基于平均值插值的多节点有限元方法分析有自由面渗流问题。在自由面附近采用多节点单元逼近自由面,利用平均值插值建立多节点单元的形函数,在远离自由面的区域采用四边形/三角形单元剖分。给定一个初始自由面位置,通过对渗流控制方程的多节点有限元求解,根据自由面上节点水头值判断自由面节点的调整方向和大小,最终迭代求出自由面的位置。土坝渗流问题的数值分析表明了所提方法的有效性和足够的计算精度。  相似文献   

某海底隧道最小岩石覆盖厚度影响因素分析     

李树忱  李术才  王兆清  王燕 《岩土力学》2006,27(Z1):261-265

目前国内正在修建或拟修建大量海底隧道,其中最小岩石覆盖厚度直接制约着海底隧道的造价。根据已确定某海底隧道最小岩石覆盖厚度,结合海底隧道所处的工程地质和水文地质条件,分析了裂隙岩体的断裂损伤、海水渗流和施工过程及爆破、地震等对海底隧道最小岩石覆盖厚度的影响分析。结合工程类比法和数值计算方法得到某海底隧道最小岩石覆盖厚度,考虑各种因素对最小岩石覆盖厚度的影响,最终分析出较优的最小岩石覆盖厚度。  相似文献