有自然面渗流分析的高斯点法

有自由面的渗流问题是比较复杂的非线性问题，这类问题的有限元求解一般采用的定网格法。本文根据饱和区和非饱和区的渗流特点，提出了高斯点法，该算法原理简单，易于实现。     

有自由面渗流分析的三维数值流形方法

提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。     

复合单元法在求解有自由面渗流问题中的应用

论文指出了在有限元分析过程中地下水自由面上、下的区别,提出可采用复合单元法处理常应变交截单元,并以四面体单元为例,推导出了自由面与单元交截的截面方程,详细说明了复合单元法的原理、系数矩阵的求解方法、推导出体积分和曲面积分的求解公式。     

基于Ansys的渗流自由面计算方法

对渗流自由面计算的单元传导矩阵调整法进行了改进。在局部坐标系下,将单元离散为一系列面元（三维为体元）,根据面元或体元形心处的压力水头值来确定面元或体元的权系数,然后按加权平均确定单元的等效渗透系数,该方法既可以避免复合单元中自由面上下面积或体积的计算,面元或体元的离散密度也可以任意选择,简化了计算程序。通过预设的收敛参数来控制迭代过程,可避免相邻两次计算的渗透系数变化太大而难于收敛。以Ansys为开发平台,用APDL语言编写了相应的计算程序,扩展了Ansys的计算功能,为渗流自由面计算提供了新途径。算例表明：该方法的收敛性较好、计算精度较高,能满足工程需要。     

三维非恒定渗流自由面边界积分项的高精度数值积分

在三维非恒定渗流有限元计算中,不可避免地需计算自由面上的边界积分项。建议一种高精度数值积分的方法求取自由面边界积分项。其基本思路是,基于8节点空间等参单元,根据压力为零的边界条件确定自由面满足的曲面方程,将自由面边界积分转化为ξη平面上的二重积分;然后再计算自由面与三维等参单元的交点,将交点投影到局部坐标平面ξη上,并根据点与线的相对关系,确定投影点所围成的局部坐标积分区域;再将总积分区域划分成若干个三角形子区域,并利用变步长Simpson方法计算各三角形子区域上的二重积分,从而实现了高精度的自由面边界积分。该方法避免了单元中自由面为平面的假设,可提高计算精度,特别对于单元内自由面变化剧烈情况,更为显著。将该方法应用于砂槽模型和土坝的渗流分析中,计算结果与实际情况较为接近。     

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

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

渗流自由面的数值分析

在岩土边坡、土坝、地下洞室及地下水运动等渗流分析中,均存在有渗流自由面问题。渗流自由面的确定,对于正确认识地下水在边坡稳定以及变形中地应用,了解地下水在边坡中的赋存特性及其运动过程等具有重大意义。本文用FLAC数值软件对影响渗流自由面的形态的相关参数,包括渗透系数、孔隙率、水头压力的大小、渗流路径长短、岩土体的各向异性程度等分别建模进行计算分析,得到相关结论。     

杂交边界点法求解带自由面的渗流问题

将杂交边界点法与迭代法相结合,求解有自由面的渗流问题。杂交边界点法基于杂交位移变分原理和移动最小二乘近似,利用基本解插值域内的场函数,而边界上的变量则用移动最小二乘近似,是一种纯边界类型的无网格方法。利用该方法只需在边界上布点而不需要划分任何网格的特性,先假定自由面的初始位置,再进行迭代求解。数值算例表明,该方法精度较高、计算量较小,适合于求解各种具有自由面的渗流问题。     

无压渗流分析的复合单元法

针对无压渗流问题,在传导矩阵调整法的基础上,对穿越自由面的单元用复合单元来处理,每个复合单元具有两套结点水头,一套用于饱和区的水头插值;另一套则用于无水区。通过变分原理,推出了控制方程并在程序中实现。应用复合单元法可以较好地解决自由面的插值问题,消除了穿过自由面单元的水上部分对自由面附近结点水头的贡献,同时对逸出面边界条件的处理也更加方便和准确。算例表明了该方法的有效性。     

改进的有自由面渗流问题的Bathe算法

建议了一个新的有自由面渗流问题的变分不等式提法,该提法通过将潜在出渗面上的边界条件提为Signorini型条件,从而从理论上消除了出渗点的奇性,解决了出渗点的定位问题。同时在离散求解时,通过引进依赖于网格参数的连续型Heaviside函数,克服了Bathe算法中所固有的网格依赖性,提高了这类方法的数值稳定性。     

饱和模型非定常渗流的数值求解方法

提出了饱和模型非稳定流的有限元算法，并对自由面流量补给项提出一种简单的处理方法，从而避免了在迭代中求解自由面的位置，最后给出了计算实例。     

不透水地基上均质土坝渗流自由面的有限元计算

为了保证土坝安全,水利工程中需要对土坝进行渗流计算,首先需要确定自由面。根据达西定律和渗流连续方程,以及水头和流量边界条件,建立了不透水地基上的均质土坝渗流的计算模型。利用Geoslope软件,采用有限元离散的方法计算了土坝的渗流自由面,并分析了土坝在不同的渗透系数情况下的自由面位置的差异。最后,解析解计算的自由面结果验证了有限元计算自由面方法是正确的。     

裂隙网络无压渗流分析的初流量法

借鉴初流量法的思想,引入初流速来抵消在裂隙网络干区实际不存在的流速,将Darcy定理扩展到整个区域,定义了在整个区域上的非线性边值问题,并将潜在溢出面边界条件归纳为Signorini型边界条件,建律了等价的变分不等式提法。通过结合连续型的Heaviside函数,给出了基于变分不等式的初流量有限元迭代算法,研制开发了相应的计算程序,通过两个典型算例验证了本文算法在求解复杂裂隙网络渗流自由面的有效性和鲁棒性。