摘要
我国是一个多山、多雨雪天气的国家,滑坡灾害分布范围很广且频率很高。此外我国处在基础建设飞速发展的阶段,交通设施、水利水电设施的建设在全国范围内全面展开,极大的改变了地球浅表层的自然状态,打破了原有的平衡状态,使得滑坡灾害越发严重。因此,快速准确的边坡稳定性评价对指导边坡工程的设计施工、保证边坡工程的安全,具有极大的意义。
目前,二维方法仍是边坡稳定分析的最主要手段。但是自然界发生的绝大多数滑坡的三维效应十分明显,三维边坡稳定分析可以更加真实地评价边坡的稳定性。以垂直条分为基础的三维极限平衡分析方法已经有了充分的发展,并且开始用于边坡工程的稳定性分析。但是由于三维极限平衡分析方法本身的复杂性,其在工程应用中仍难以大规模推广。本文针对该问题,进行了三维极限平衡分析软件开发技术的研究,主要内容有以下几个方面:
(1)选取水科院三维极限平衡分析方法作为应用技术研究的基础,详细推导了该方法的求解过程,并提出了迭代初值的确定方法。
(2)提出了Delaunay三角网数字地层模型加下覆土层规则构建坡体模型的方法,并设计了一种构建Delaunay三角网的新型生长法——壳外插入法。
(3)设计了一种新型的三维滑裂面——双二次Bezier滑裂面。基于垂直条分法的特点,提出了基于格网数据的三维滑体构造方法。
(4)建立了三维临界滑裂面求解的目标函数,根据其特点确定了以遗传算法为基础的搜索方法。收集整理了改进遗传算法性能的研究成果。对遗传算法的实现技术进行了研究,设计了内存二进制编码技术和基于多线程的并行遗传技术,同时以面向对象方法设计了遗传算法实现框架。
(5)研究了三维边坡稳定分析应用技术的集成方法,设计了三维边坡稳定分析软件的开发构架,并以该框架设计了Ramtil-SLOPE应用程序。
In our country which is mountainous rainy and snowy landslide is widely and high-frequency happened. Moreover our country is at the stage of high-speed developing. Constructions of transportation establishments and hydropower stations are progressing all over the nation. These constructions caused more and more landslide disaster because of greatly destroying the natural state of the land surface and the balance of slopes. So rapidly and exactly evaluating the stability of slopes is meaningful for supervising the construction and design and ensuring the safety of slope projects.
At present, 2D method is still the main means for analysis of slope stability. But 3D effect is very evidence for mass landslides in the nature, so 3D analysis of slope stability can be more trustful. Nowadays, 3D limit equilibrium analysis methods which based on vertical slice method are highly developed and start used in analysis of slope stability. But 3D limit equilibrium analysis methods are hard to spread widely in slope project for the sake of complexity. For solving this problem, software developing techniques researches are presented in this paper. The main research contents can summarized as follows:
(1) IWHR method had been chosen as the base method for research. The solution procedure of the method had been derived. The method for sampling initial value of Newton-Raphson method which is for solving the limit equilibrium formulas was also presented.
(2) The modeling method using Delaunay Triangulation and surface coverage rule for describing 3D slope was presented. More over a new growth method for creating Delaunay Triangulation was presented.
(3) A new type of slip surface called two dimensional quadric Bezier slip surface was designed. For the characteristic of vertical slice method a method based on grid data for creating 3D slide mass model was designed.
(4) Objective function was setup for searching the critical slip surface and Genetic Algorithm (GA) was chosen as the searching method based on its characteristic. The main achievement of GA research work was summarized. The research works for application techniques of GA were done. Memory binary code and parallel GA based on Multi-Thread techniques were researched. And the architecture of GA application was designed using the Object-Oriented Design method.
(5) The methods for integrating the application techniques of 3D slope stability analysis were researched. The architecture was designed for developing 3D slope stability software and the Ramtil-SLOPE application was developed based on the architecture.
引文
[1]Duncan,J.M.State of the art:Limit equilibrium and finite element analysis of slopes[J].Journal of Geotechnical Engineering.1996,122(7):577-596.
[2]Stark,T.D.,Eid H.T.Performance of three-dimensional slope stability analysis method in Practice[J].ASCE Journal of Geotechnical Engineering.1998,(124):1049-1060.
[3]Lam.,L.,D.G.Fredlund.A general limit equilibrium model for three-dimensinal slope stability analysis[J].Canadian Geotechnical Journal.1993,30:905-919.
[4]Chen,R.H.,Chameau J.L."Three-Dimensional Slope Stability Analysis," presented at Numerical Methods in Geomechanics Edmonton,1982.
[5]Leshchinsky,D.,Baker R.,Silver M.L.Three Dimensional Analysis of Slope Stability[J].International Journal for Numerical and Analytical Methods in Geomechanics.1985,3(9):199-223.
[6]陈祖煜.土质边坡稳定分析--原理、方法、程序[M].北京:中国水利水电出版社,2003.
[7]Hungr,O.,Salgado E M.,Byme P.M.Evaluation of a three-dimensional method of slope stability analysis[J].Canadian Geotechniacal Journal.1989,26:679-686.
[8]ZhangXing.Three-dimensional stability analysis of concave slopes in plan view[J].Jouornal of Geotechnical Engineering,ASCE.1988(114):658-671.
[9]Hutchison,J.N.,Sarma S.K.Discussion on three-dimensional limit equilibrium analysis of slopes[J].Geotechnique.1985,35:215.
[10]Huang,C.C.,Tsai C.C.New method for 3D and asymmetrical slope stability analysis[J].Journal Geotechical and Engineering.2000,10(126):917-927.
[11]冯树仁,丰定祥,葛修润,谷先荣.边坡稳定性的三维极限平衡分析方法及应用[J].岩土工程学报.1999,21(6):657-661.
[12]谢谟文,周国云,江崎哲郎,三谷浩泰.基于GIS空间数据库的三维边坡稳定性分析[J].岩石力学与工程学报.2002,10(21):1494-1499.
[13]姜清辉,王笑海,丰定祥,冯树仁.三维边坡稳定性极限平衡分析系统软件SLOPE3D的设计及应用[J].岩石力学与工程学报.2003,7(22):1121-1125.
[14]孙平.基于非相关流动法则的三维边坡稳定极限分析[D].中国水利水电科学研究院,北京,2005.
[15]顾晓强.边坡稳定分析方法及其应用研究[D].上海交通大学,上海,2007.
[16]曹亚星.饱和-非饱和土坡稳定性的三维极限平衡分析[D].天津大学,天津,2006.
[17]张常亮.边坡最危险滑动面的三维搜索方法研究[D].长安大学,西安,2005.
[18]谢谟文,蔡美峰.信息边坡工程学的理论与实践[M].北京:科学出版社,2005.
[19]Chen,Z.,Wang J.,Wang Y.,Haberfield C.A three-dimensional slope stability analysis method using the upper bound theorem,Part II:Numerical approaches,applications and extensions.[J].International Journal of Rock Mechanics and Mining Sciences.2001,(38):379-397.
[20]张浩,王常明,马文东.边坡最危险滑动面全局搜索的模拟退火算法及改进[J].吉林大学学报(地球科学版).2007,1(37):129-133.
[21]吴余生,陈胜宏.并行组合模拟退火算法在边坡稳定分析中的应用[J].岩土力学.2006,9(27):1554-1558.
[22]张发明,郭丙跃,庞正江,周祖权,周民权.基于模拟退火算法的边坡稳定分析方法[J].岩土力学.2006,4(27):514-518.
[23]李亮,迟世春,林皋.禁忌模拟退火复合形法及其在边坡稳定性分析中的应用[J].岩石力学与工程学报.2005,18(24):3342-3349.
[24]吕文杰,李晓军,朱合华.基于遗传算法的边坡稳定分析通用算法[J].岩土工程学报.2005,5(27):595-599.
[25]李亮,迟世春,郑榕明,林皋.一种新型遗传算法及其在土坡任意滑动面确定中的应用[J].水利学报.2007,2(38):157-162.
[26]何则干,陈胜宏.遗传模拟退火算法在边坡稳定分析中的应用[J].岩土力学.2004,2(25):316-319.
[27]包太,刘宝琛.遗传算法在边坡稳定系数计算中的应用[J].重庆建筑大学学报.2007,1(29):49-51.
[28]陈昌富.混沌扰动启发式蚁群算法及其在边坡非圆弧临界滑动面搜索中的应用[J].岩石力学与工程学报.2004,20(23):3450-3453.
[29]陈昌富,龚晓南,王贻荪.自适应蚁群算法及其在边坡工程中的应用[J].浙江大学学报(工学版).2003,5(37):566-569.
[30]李亮,迟世春,林皋.基于蚁群算法的复合形法及其在边坡稳定分析中的应用[J].岩土工程学报.2004,5(26):691-696.
[31]陈云敏,魏新江,李育超.边坡非圆弧临界滑面的粒子群优化算法[J].岩石力学与工程学报.2006,7(25):1443一1449.
[32]Voronoi,G Nouvelles Applications des Parameters Continus,a la Theorie des Formes Quadratiques,Deuxieme Memorie:Recherches sur les Parrallelloedres Primitifs[J].Jounal fur die Reine und Angewandte Mathematik.1908,(134):198-287.
[33]Green,P.J.,Sibson.R.Computing Dirichlet Tesselations in the Plane[J],The Computer Journal.1978,2(21):168-173.
[34]Brassel,K.E.,Reif D.Procedure to Generate Thiessen Polygons[J].Geophysical Analysis.1979,11:289-303.
[35]McCaullagh,M.T.,Ross C.G.Delaunay Triangulation of a Random Data Set For lirarithrn ic Mapping[J].The Cartographic Journal.1980 17:93-99.
[36]Maus,A.Delaunay Triangulat ion and the Convex Hull of n Points in Expected Linear Time[J].BIT.1984,24:151-163.
[37]Lawson,C.L.Software for C' Surface Interpolation Mathematical Software Ⅲ[M].New York:Academic Press,1977.
[38]Watson,D.F.Computing the n-dimension Delaunay Tesselation with Application to Voronoi Polytopes [J].Computer Journal.1981,24:167-172.
[39]Lewis,B.A.,Robinson J.S.Triangulation of Planar Regions with Applications[J].The Computer Journal.1978,4(21):332-332.
[40]Lee,D.T.,Schachter B.J.Two A Igorithms for Constructing a Delaunay Triangulation[J].Int.J.of Computer and Information Sciences.1980,3(9).
[41]Fortune,Steven.Sweepline Algorithm For Voronoi Diagrams[J].Algorithrnica.1987,2(2):153-174.
[42]Bowyer,A.Computing Dirichlet Tesselations[J].Computer Journal.1981 24:162-166.
[43]蒋红斐,涂鹏,李国忠.基于生长算法构建Delaunay三角网的研究[J].公路交通科技.2004,12(21):38-41.
[44]凌海滨,吴兵.改进的自连接Delaunay三角网生长算法[J].计算机应用.1999,12(19): 10-12.
[45]徐青,马旭东,常歌,杨力.TIN建立及其向GRID转换优化算法[J].解放军测绘学院学报.1999,4(16):263-266.
[46]宋占峰,蒲浩,詹振炎.快速构建Delaunay三角网算法研究[J].铁道学报.2001,5(23):85-91.
[47]詹振炎.道路规划与设计原理[M].长沙:中南大学出版社,2004.
[48]钱家欢,殷宗泽.土工原理与计算[M].北京:中国水利水电出版社,1996.
[49]谢谟文,江崎哲郎,周国云.基于边坡单元的三维滑坡灾害评价的GIS方法[J].岩石力学与工程学报.2004,6(22):969-976.
[50]David,E R.Procedureal Elements for Computer Graphics(2E)[M].New York:McGraw-Hill Companies,1998.
[51]Holland,J.H.Adaptation in Nature and Artificial System[M].MIT Press,1992.
[52]De Jong,K.A.An Analysis of the Behavior of a Class of Genetic Adapive System[D].University of Michigan,Michigan,1975.
[53]Goldberg,D.E.Genetic Algorithms in Search,Optimization and Machine Learning[M].Reading:Addison-Wesley,1989.
[54]周明,孙树栋.遗传算法原理及应用[M].北京:国防工业出版社,2001.
[55]Srinivas,M.Genetic Algorithms:a Survey[J].Computer Journal.1994,7:17-26.
[56]Cavicchio,D.J.Adaptive Search Using Simulated Evolution[D].University of Michigan,1970.
[57]金菊良,杨晓华,丁晶.标准遗传算法的改进方案--加速遗传算法[J].系统工程理论与实践.2001,4:8-13.
[58]Hockney,R.W.,Jesshope C.R.Parallel Computers Architecture,Programming and Algorithms[M].Bristol:Adam Hilger Ltd.,1981.
[59]兰博,雅格布森,博什.UML参考手册[M].北京:机械工业出版社,2005.
[60]奥多彻蒂.面向对象分析与设计[M].北京:清华大学出版社,2006.
[61]刘艺.Delphi模式编程[M].北京:机械工业出版社,2004.
[62]Floudas,C.A.,Pardalos P.M.,"A Collection Test Problems for Constrained Global Optimaization Algorithms," in LNCS,vol.455,Springer-Verlag ed,1987.
[63]Hock,W.,Schittkowski K.,"Test Examples for Nonlinear Programming Codes," in Lecture Notes in Economics and Mathematical Systems,vol.187,Springer-Verlag ed,1981.
[64]黄云锴,黄艳华.GLScene在地形可视化应用中的研究[J].昆明理工大学学报(理工版).2007,2(32):7-11.
[65]欧训勇.基于GLScene三维图形引擎的教学资源开发与应用研究[D].华中师范大学,武汉,2006.