心墙土石坝渗流自由面的数值模拟

？心墙与坝壳材料渗透系数相差比较大,使用有限元方法选代计算心墙土石坝渗流自由面时,往往会在两种材料交界面处失真,难以正确模拟出心墙下游端的出渗点.为解决此问题,采用具有严格理论基础的signorini型变分不等式,确定渗流渗出点和自由面的位置,可消除出渗点的奇异性,采用罚heaviside函数克服网格依赖性和数值不稳定性.同时,根据毛细水上升高度来确定非饱和区负孔隙水压力,消除饱和-非饱和渗流分析中坝壳节点中的虚值,使得计算结果能准确地反映出坝体中的实际渗流场.