In the wind-water erosion crisscross region on the Loess Plateau,the existence of stone fragments influences soil water hydraulic properties and restricts plant water efficiency.Therefore,results of this investigation for seepage flow of soil-stone mixtures could be useful to vegetation restoration in this region.Although there are a lot of researches on the seepage flow of underground water,few of them are for soil stone mixture because of its largely increasing scale of FEM solution.By employing rectangular element of sub-domain technique of finite element model and assuming that the increase of stone content in the mixture has no effects on porosity of soil and that stone is a material of no seepage,this work has analyzed the seepage flow of soil stone mixture.The computational results of this model show that the main factor which affects saturated water conductivity of soil stone mixture is the stone contents,while the stone size almost has no influence on it.However the distribution layout of stone in soil of this model is the second affecting factor for the average seepage flow velocity if the stone layout increases the length of seepage flow path.And the numerical achievements are in good agreement with some author's theoretical results and experimental results.The sub-domain method is of high computational performance and accuracy.And it can be used for analyzing large scale seepage flow in soil stone mixture.In addition,the method could be generalized for 3D analysis of soil stone mixture and for larger scale simulation of seepage by multi-sub-domain technique.