Abstract：Seepage analysis has found wide application in many areas. In practice, the conventional numerical tools such as the finite element method are mainly used in the seepage analysis, which often demand large computational resources and therefore impose restrictions on the problem scale. The weak form quadrature element method is a simple and efficient numerical tool which has been applied to structural analysis. Based on the weak form description of a problem, it has the characteristics of global approximation and enjoys rapid convergence. The method is used for simulation of two- and three-dimensional confined and unconfined seepage. For the unconfined seepage, the adaptive mesh method is employed, and the free surface is expressed by polynomial interpolation at integration points. The results are compared with those of other methods and good agreement is reached. It is shown that a relatively small number of degrees of freedom are needed to attain convergence by the quadrature element formulation as compared with those of the finite element method. The weak form quadrature element method is expected to be an effective numerical tool for seepage analysis.
 AHMED A, BAZARAA A. Three-dimensional analysis of seepage below and around hydraulic structures[J] . Journal of Hydrologic Engineering, 2009, 14(3): 243-247.  BATHE K J, KHOSHGOFTAAR M R. Finite element free surface seepage analysis without mesh iteration[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1979, 3(1): 13-22.  KAZEMZADEH-PARSI M J, DANESHMAND F. Unconfined seepage analysis in earth dams using smoothed fixed grid finite element method[J] . International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(6): 780-797.  KAZEMZADEH-PARSI M J, DANESHMAND F. Three dimensional smoothed fixed grid finite element method for the solution of unconfined seepage problems[J] . Finite Elements in Analysis and Design, 2013, 64(1): 24-35.  SAKURAI H, KAWAHARA M. Three-dimensional groundwater flow analysis system using the element-free galerkin method[J]. International Journal of Computational Fluid Dynamics, 2004, 18(4): 309-315.  JIE Y, JIE G, MAO Z, et al. Seepage analysis based on boundary-fitted coordinate transformation method[J]. Computers and Geotechnics, 2004, 31(4): 279-283.  RAFIEZADEH K, ATAIE-ASHTIANI B. Three dimensional flow in anisotropic zoned porous media using boundary element method[J]. Engineering Analysis with Boundary Elements, 2012, 36(5): 812-824.  ZHONG H, YU T. A weak form quadrature element method for plane elasticity problems[J]. Applied Mathematical Modelling, 2009, 33(10): 3801-3814.  ZHONG H, YU T. Flexural vibration analysis of an eccentric annular Mindlin plate[J]. Archive of Applied Mechanics, 2007, 77(4): 185-195.  MO Y, OU L, ZHONG H. Vibration analysis of timoshenko beams on a nonlinear elastic foundation[J]. Tsinghua Science & Technology, 2009, 14(3): 322-326.  HE R, ZHONG H. Large deflection elasto-plastic analysis of frames using the weak-form quadrature element method[J]. Finite Elements in Analysis and Design, 2012, 50(1): 125-133.  SHU C. Differential quadrature and its applications in engineering[M]. Berlin, Hoidelerg: Springer, 2000.  POLUBARINOVA-KOCHINA P K. Theory of groundwater movement[M]. Princeton: Princeton University Press, 1962: 506-509.  BRESCIANI E, DAVY P, DE DREUZY J R. A finite volume approach with local adaptation scheme for the simulation of free surface flow in porous media[J]. International Journal for Numerical and Analytical Methods in Geomechanics. 2012, 36(13): 1574-1591.  ELLIOTT C, OCKENDON J R. Weak and variational methods for moving boundary problems[M]. Boston: Pitman Pub, 1982.  WESTBROOK D R. Analysis of inequality and residual flow procedures and an iterative scheme for free surface seepage[J]. International Journal for Numerical Methods in Engineering. 1985, 21(10): 1791-1802.