土石坝渗流分析中Chebyshev谱元法.docVIP

土石坝渗流分析中Chebyshev谱元法 　　摘要：基于谱元法理论，对土石坝一维二维渗流问题进行了理论分析。首先建立了与土石坝渗流微分方程初边值问题等价的积分形式，选取Chebyshev多项式作为函数插值的基函数，得出渗流问题的离散矩阵方程。然后采用与有限元相似的方法集成整体系数矩阵，根据不同的边界条件，最终计算得到结点水头列阵{h}，通过谱插值点来求出整个土石坝各点的水头值。 　　关键词：谱元法；渗流；切比雪夫多项式；有限元 　　中图分类号：TV139.1；TV641 文献标识码：A 文章编号： 　　1672-1683（2011）01-0169-04 　　A Chebyshev Spectral Elements Method for Earth Dam Seepage Analysis 　　CHANG Ming-yuan，SONG Zi-heng，LI Nan-sheng，NIU Yong-chang 　　(Department of Hydraulic Engineering,Tongji University,Shanghai 200092,China) 　　 　　Abstract: Based on the theory of spectral elements method,this paper theoretically analyzed the issue of earth dam seepage under both one-dimensional and two-dimensional conditions. First,an integral form about the seepage problem was set up;then,the discrete matrix equations were formed by using Chebyshev polynomials as the basis function and the global coefficient matrix was formed in the way used in the finite-element method.Finally,according to different types of boundary conditions,water head linear array {h} was calculated and the whole dams water head was educed from this linear array. 　　Key words: Spectral elements method;seepage;Chebyshev Polynomial;finite element method 　　 　　 渗流问题是影响土石坝工程安全性的最重要的因素之一，因而渗流问题的研究具有重要的理论意义和工程应用价值。谱元方法是用谱近似技术结合有限元的思想，在单元内部仍然是使用加权参量法。谱元法和有限元两种方法最根本的区别是展开函数的选取，有限元是基于固定的低阶插值函数，是靠增加单元数目来提高精度；而谱方法是基于总体插值，是靠增加插值的阶数来提高精度，谱方法还有一个优点是收敛速度快（指数型）。这两个概念的成功结合，使谱元法即可达到较高的求解精度，又可适应复杂的几何边界［1］。因而本文基于谱元法对土石坝一维、二维渗流问题进行了理论分析。?? 　　1 谱方法求解一维渗流计算问题?? 　　1.1 一维稳定渗流?? 　　在坝体宽度较大下游坝坡防渗条件较好，而且没有坝顶水流补给情况下，应用Dupuit假设，可以将二维渗流问题简化为一维稳定渗流［2］。一维渗流方程有：?? 　　假设隔水底板水平，对上式直接积分得：?? 　　式（2）中C1，C2为积分常数，应用边界条件可确定：?? 　　故得一维稳定渗流解：?? 　　1.2 一维非稳定渗流问题?? 　　如果坝顶有水流补给w(x,t)，此时坝体渗流必须考虑时间因素，问题变成为非稳定问题［2］，即：?? 　　对于一维非稳定渗流方程（4），以下采用高精度谱元法对其进行近似数值计算。对未知水头函数h(x,t)用谱多项式对其进行数值逼近，也就是逼近展开的基函数采用完备的正交多项式――谱多项式。根据一维非稳定渗流方程的类型，此处应用Chebyshev多项式对水头函数h(x,t)进行谱逼近，即假定未知水头函数h(x,t)渐进解取如下展开形式：?? 　　hN(x,t)=??Nm=0am(t)Tm(x)（5） 　　式中：N-渐进展开项数，其整数值大小由计算精度和效率综合决定，Tm(x)=cos(marccosx)为m级Chebyshev多项式［3