三峡永久船闸高边坡开挖变形的离散元法分析
|
摘要
由于岩体中的结构面主导着岩石边坡的破坏机理，因此基于连续性理论的分析方法在分析岩石边坡的变形特征具有相当的局限性。而基于非连续介质力学的离散单元法因为不仅能够很好地表征岩土体材料的非均匀性和各向异性，而且能很好地模拟碎块的形成及运动。本文重点是以颗粒元为基础，建立离散元框架内的网格实体模型，并自行开发数值计算软件，结合典型工程背景进行岩石高边坡变形的研究和探索。
本文主要进行了如下工作：
(1)以颗粒元的离散单元法为基础，提出离散元框架内的网格实体模型。网格实体模型是指用颗粒单元将研究对象离散化，用接触发现算法确定每个颗粒单元的邻居单元，并在所有相邻单元间施加“梁”，这样在整个介质中就形成了一个“梁”的网络。通过所建立的强度准则，判别“梁”是否存在，利用矩阵位移法求解“梁”的受力，用离散单元法描述颗粒运动。
(2)采用C语言开发了基于离散元框架的网格实体模型程序DEM_LSM2D。计算软件的前、后处理模块由C调用OpenGL函数库实现。软件主界面和功能扩展模块(图片浏览器和批处理运算控制器)由VB6.0开发完成。
(3)应用DEM_LSM2D模拟三峡永久船闸高边坡的变形。在卸荷开挖条件下计算船闸高边坡的水平向和竖向的变形量。监测数据与模拟位移结果的对比表明该方法对卸荷条件下岩体的变形规律模拟是非常有效的。
Where a rock slope comprises multiple joint sets, which control the mechanism of failure, continuum modeling is not so best suited for the displacement analysis of slopes. A discontinuum modeling approach may be considered more appropriate, because it not only can token the nonuniformity and heterogeneity of the geo-materials, but also be good at modeling the generation and movement of the fragments. Based on particle element model, a lattice solid model used in distinct element method is established and applied to rock slope stability analysis in three gorges project in this paper.
The main research focuses are as followed:
(1) Based on particle element model, a lattice solid model (LSM) used in distinct element method is proposed. In the Lattice solid model, the elements are distinct and detected by collision detection method. The neighbor elements are connected by the beams, which help to establish a net of beams. Through the strength criterion, the existence of the beam is detected. Matrix Displacement Method is used to solve the force applied to the beams, and the movements of the particles are described by the distinct element method.
(2)A C program called DEM_LSM2D is developed, in which the Lattice Solid Model is adopted. The pre-process and post-process of the DEM_LSM2D is programmed by C and OpenGL function library. The interface and the of function extension module is finished by VB6.
(3)The deformation of high slope in the Three Gorges project is simulated by the DEM_LSM2D. Taking excavation into consideration, the displacements on the horizontal and vertical are calculated. Comparing the field data and the modeling result, it shows that the Lattice solid model method is a proper way to model the rock deformation under unloading conditions.
本文主要进行了如下工作：
(1)以颗粒元的离散单元法为基础，提出离散元框架内的网格实体模型。网格实体模型是指用颗粒单元将研究对象离散化，用接触发现算法确定每个颗粒单元的邻居单元，并在所有相邻单元间施加“梁”，这样在整个介质中就形成了一个“梁”的网络。通过所建立的强度准则，判别“梁”是否存在，利用矩阵位移法求解“梁”的受力，用离散单元法描述颗粒运动。
(2)采用C语言开发了基于离散元框架的网格实体模型程序DEM_LSM2D。计算软件的前、后处理模块由C调用OpenGL函数库实现。软件主界面和功能扩展模块(图片浏览器和批处理运算控制器)由VB6.0开发完成。
(3)应用DEM_LSM2D模拟三峡永久船闸高边坡的变形。在卸荷开挖条件下计算船闸高边坡的水平向和竖向的变形量。监测数据与模拟位移结果的对比表明该方法对卸荷条件下岩体的变形规律模拟是非常有效的。
Where a rock slope comprises multiple joint sets, which control the mechanism of failure, continuum modeling is not so best suited for the displacement analysis of slopes. A discontinuum modeling approach may be considered more appropriate, because it not only can token the nonuniformity and heterogeneity of the geo-materials, but also be good at modeling the generation and movement of the fragments. Based on particle element model, a lattice solid model used in distinct element method is established and applied to rock slope stability analysis in three gorges project in this paper.
The main research focuses are as followed:
(1) Based on particle element model, a lattice solid model (LSM) used in distinct element method is proposed. In the Lattice solid model, the elements are distinct and detected by collision detection method. The neighbor elements are connected by the beams, which help to establish a net of beams. Through the strength criterion, the existence of the beam is detected. Matrix Displacement Method is used to solve the force applied to the beams, and the movements of the particles are described by the distinct element method.
(2)A C program called DEM_LSM2D is developed, in which the Lattice Solid Model is adopted. The pre-process and post-process of the DEM_LSM2D is programmed by C and OpenGL function library. The interface and the of function extension module is finished by VB6.
(3)The deformation of high slope in the Three Gorges project is simulated by the DEM_LSM2D. Taking excavation into consideration, the displacements on the horizontal and vertical are calculated. Comparing the field data and the modeling result, it shows that the Lattice solid model method is a proper way to model the rock deformation under unloading conditions.
引文
[1] 高大钊．岩土工程的回顾与前瞻[M]．北京：人民交通出版社，2001．
[2] 高长胜，徐光明，张凌等．边坡变形破坏离心机模型试验研究[J]．岩土工程学报．2005，27(4)：478～481．
[3] 李桂荣，余成学，陈胜宏．层状岩体边坡的弯曲变形破坏试验及有限元分析[J]．岩石力学与工程学报．305～311．
[4] 石安池，赵明华，薛果夫．三峡工程永久船闸高边坡岩体变形特征与机理分析[J]．岩石力学与工程学报．2001，20(5)：638～642．
[5] 石豫川，冯文凯，刘汉超等．某水电站高边坡变形破坏模式及机制分析[J]．西南交通大学学报．2004，39(5)：609～613．
[6] 邓建辉，D．H．Chan，C．D．Martin等．一例由蓄水诱发德库岸边坡变形[J]．中国地质灾害与防治学报．2002，13(1)：21～24．
[7] 杨杰，胡德秀，关文海．李家峡拱坝左岸高边坡岩体变位与安全性态分析[J]．岩石力学与工程学报．2005，24(19)：3551～3560．
[8] 陆晓敏，任青文，盛芳．裂隙岩质边坡的弹粘塑性变形及稳定性分析[J]．岩石力学与工程学报．21(4)：493～497．
[9] 程桦，孙钧．三峡船闸及高边坡非线性大变形数值分析[J]．岩土力学．1998，19(4)：1～7．
[10] 郑东健，顾冲时，吴中如．边坡变形的多因素时变预测模型[J]．岩石力学与工程学报．2005，24(17)：3180～3184．
[11] 蒋刚，林鲁生，刘祖德等．边坡变形德灰色预测模型[J]．岩土力学．2000，21(3)：244～251．
[12] 关秦川，张志勇，冯浩．大型干坞边坡变形及其神经网络预测模型[J]．西南交通大学学报．2004，39(2)：157～161．
[13] 黄银，刘俊，潘翔．加权灰色与加权统计模型在边坡变形预测中的应用研究[J]．长江科学院院报．2005，22(4)：76～79．
[14] 杨学堂，哈秋黔，高希章等．厚层软岩与硬岩互层岩体高边坡卸荷变形与支护研究[J]．岩石力学与工程学报．2004，23(16)：2681～2686．
[15] 唐建新，陶洪久，章光．利用遥感图像对滑坡进行调查和监测分析研究[J]．岩石力学与工程学报．2002，21(增2)：2515～2520．
[16] 王旭春，何满潮，蒋宇静等．在GIS中实现基于力学原理的滑坡稳定性计算[J]．岩石力学与工程学报．2003，22(6)：977～980．
[17] 欧阳祖熙，张宗润，丁凯等．基于3S技术和地面变形观测的三峡库区典型地段滑坡监测系统[J]．岩石力学与工程学报．2005，24(18)：3203～3210．
[18] 周翠英，刘祚秋，董立国等．边坡变形破坏过程的大变形有限元分析[J]．岩土力学．2003，24(4)：644～652．
[19] 徐平，甘军，丁秀丽等．三峡工程船闸高边坡岩体长期变形及稳定性有限元分析．长江科学院院报．1999，16(2)：31～34．
[20] 刘文平，赵燕明，郑颖人．岩质边坡开挖应力与变形的有限元模拟[J]．后勤工程学院学报．2004，2：45～48．
[21] 刘凯欣，高凌天．离散元法研究的评述[J]．力学进展．2003，33(4)：483～490．
[22] 俞良群，邢纪波，徐建等．岩质高陡边坡失稳的离散元法与非连续变形法的对比分析[J]．烟台大学学报(自然科学与工程版)．1999，12(1)：63～66．
[23] 强天驰．不连续变形分析方法(DDA)在边坡稳定分析中的应用[J]．岩土工程界．2004，7(6)：40～42．
[24] 王书法，李树忱，李术才等．节理岩质边坡变形的DDA模拟[J]．岩土力学．2002，23(3)：352～354．
[25] 吴建宏，大西有三，石根华等．三维非连续变形分析(3D DDA)理论及其在岩石边坡失稳数值仿真中的应用[J]岩石力学与工程学报．2003，22(6)：937-942．
[26] Chert, G & Ohnishi, Y. (1999). Slope stability using discontinuous deformation analysis method. In Amadei et al. (eds.), Proceedings of the 37th U.S. Rock Mechanics Symposium, Vail. A.A. Balkema, Rotterdam, 535-541.
[27] Amadei, B., Lin, C.T., Sture, S. & Jung, J. (1994). Modeling fracturing of rock masses with the DDA method. In Nelson & Laubach (eds.), Proceedings of the 1st North American Rock Mechanics Symposium, Austin. A. A. Balkema, Rotterdam, 583-590.
[28] 杨天鸿，唐春安，郑雨天等．FLAC程序在抚顺西露天矿边坡变形治理工程中的应用[J]．地质灾害与环境保护．1999，10(3)：6～11．
[29] 虢曙安，贺春宁，芮勇勤．开挖边坡变形滑动的FLAC数值模拟与分析[J]．湖南交通科技．2003，29(3)：4～5，15．
[30] 罗海宁，焦玉勇．对三维离散单元法中块体接触判断算法的改进[J]．岩土力学．1999，20(2)：37～40，45．
[31] 王泳嘉，刘连峰．三维离散单元法软件系统TRUDEC的研制[J]．岩石力学与工程学报．1996，15(3)：201～210．
[32] 胥建龙，唐志平．离散元与有限元结合的多尺度方法及其应用[J]．计算物理．2003，20(6)：477～482．
[33] 金峰，贾伟伟，王光纶．离散元—边界元动力耦合模型[J]．水利学报．2001，1：23～27．
[34] 李世海，汪远年．三维离散元计算参数选取方法研究[J]．岩石力学与工程学报．2004，23(21)：3642～3651．
[35] 李本平．烟囱定向爆破仿真模拟[J]．爆破．2003，20(3)：11～13．
[36] 方新秋，万德钧，王庆等．离散元法在分析综放采场矿压中的应用[J]．湘潭矿业学院学报．2003，18(4)：11～14．
[37] 徐泳，李艳洁，李红艳．离散元法在农业机械化中应用评述[J]．2004，5：26～30．
[38] 焦玉勇，葛修润，刘泉声等．三维离散单元法及其在滑坡分析中的应用[J]．岩土工程学报．2000，22(1)：101～104．
[39] 郑文刚，刘凯欣．离散元法工程计算软件的前后处理系统[J]．计算机工程与科学．2000，22(6)：14，15，23．
[40] Ching-Lung Liao、Ta-Peng Chang and Dong-Hwa Young. Elasto-plastic constitutive modeling of granular assembly, Int. J. Solids & Structures. 1995, 32(8-9): 1121-1133.
[41] Ching S. Chang, Sao J. Chao Y. Chang. Estimates of elastic moduli for granular material with anisotropic random packing structure. Int. J. Solids & Structures. 1995, 32(14): 1989-2008.
[42] Ching S. Chang and Jian Gao. Second-gradient constitutive theory for granular material with random packing structure, Int. J. Solids & Structures. 1995, 32(16):2279-2293.
[43] Ta-Peng Chang and Dong-Hwa Young. Stress-strain relationship for granular materials based on the hypothesis of best fit, Int. J. Solids & Structures. 1997, 34(31-32):4087-4100.
[44] F. Emeriault and B. Cambou. Micromechanical modeling of anisotropic non-linear elasticity of granular medium, Int. J. Solids & Structures. 1996, 33(18):2591-2607.
[45] J. bardet and I. Vardoulakis. The asymmetry of stress in granular media, Int. J. Solids & Structures. 2001, 38:353-367.
[46] 黄松元．《散体力学》，机械工业出版社，1992．
[47] 赵彭年．《松散介质力学》，地震出版社，1995．
[48] F. Kun and H.J. Herrmarm: Transition from Damage to Fragmentation in Collision of Solids,Physical Review E, 59 (3/1999), p. 2623-2632.
[49] Kun, F, Herrmann, H. J. A study of fragmentation process using a discrete element method[J]. Computer Methods in Applied Mechanics and Engineering, 1996, 138:3-18.
[50] 邢纪波．梁-颗粒模型导论．北京：地震出版社，1999．
[51] G.A. D'Addetta, F. Kun, E. Ramm and H.J. Herrmann. From Solids to Granulates-Discrete Element Simulations of Fracture and Fragmentation Processes in Geomaterials, in P.A. Vermeer etal.(Eds.), Continuous and Discontinuous Modeling of Cohesive-Frictional Materials, Lecture Notes inPhysics, Springer (2001), p. 231-258.
[52] Cundall, P. A. and Strack, O. D. L. A Discrete Numerical Model for Granular Assemblies[J]. Geotechnique, 1979 (29): 47-65.
[53] Schinner, A. Fast algorithms for the simulation of polygonal particles[J]. Granular Matter, 1999, 2: 35～43.
[54] 谢祚水，王自力，吴剑国．计算结构力学[M]．武汉：华中科技大学出版社，2004．
[55] 长江勘测规划设计研究院，长江科学院等．长江三峡工程永久船闸边坡变形与稳定分析及预报研究．武汉：长江勘测规划设计研究院，2000．
[2] 高长胜，徐光明，张凌等．边坡变形破坏离心机模型试验研究[J]．岩土工程学报．2005，27(4)：478～481．
[3] 李桂荣，余成学，陈胜宏．层状岩体边坡的弯曲变形破坏试验及有限元分析[J]．岩石力学与工程学报．305～311．
[4] 石安池，赵明华，薛果夫．三峡工程永久船闸高边坡岩体变形特征与机理分析[J]．岩石力学与工程学报．2001，20(5)：638～642．
[5] 石豫川，冯文凯，刘汉超等．某水电站高边坡变形破坏模式及机制分析[J]．西南交通大学学报．2004，39(5)：609～613．
[6] 邓建辉，D．H．Chan，C．D．Martin等．一例由蓄水诱发德库岸边坡变形[J]．中国地质灾害与防治学报．2002，13(1)：21～24．
[7] 杨杰，胡德秀，关文海．李家峡拱坝左岸高边坡岩体变位与安全性态分析[J]．岩石力学与工程学报．2005，24(19)：3551～3560．
[8] 陆晓敏，任青文，盛芳．裂隙岩质边坡的弹粘塑性变形及稳定性分析[J]．岩石力学与工程学报．21(4)：493～497．
[9] 程桦，孙钧．三峡船闸及高边坡非线性大变形数值分析[J]．岩土力学．1998，19(4)：1～7．
[10] 郑东健，顾冲时，吴中如．边坡变形的多因素时变预测模型[J]．岩石力学与工程学报．2005，24(17)：3180～3184．
[11] 蒋刚，林鲁生，刘祖德等．边坡变形德灰色预测模型[J]．岩土力学．2000，21(3)：244～251．
[12] 关秦川，张志勇，冯浩．大型干坞边坡变形及其神经网络预测模型[J]．西南交通大学学报．2004，39(2)：157～161．
[13] 黄银，刘俊，潘翔．加权灰色与加权统计模型在边坡变形预测中的应用研究[J]．长江科学院院报．2005，22(4)：76～79．
[14] 杨学堂，哈秋黔，高希章等．厚层软岩与硬岩互层岩体高边坡卸荷变形与支护研究[J]．岩石力学与工程学报．2004，23(16)：2681～2686．
[15] 唐建新，陶洪久，章光．利用遥感图像对滑坡进行调查和监测分析研究[J]．岩石力学与工程学报．2002，21(增2)：2515～2520．
[16] 王旭春，何满潮，蒋宇静等．在GIS中实现基于力学原理的滑坡稳定性计算[J]．岩石力学与工程学报．2003，22(6)：977～980．
[17] 欧阳祖熙，张宗润，丁凯等．基于3S技术和地面变形观测的三峡库区典型地段滑坡监测系统[J]．岩石力学与工程学报．2005，24(18)：3203～3210．
[18] 周翠英，刘祚秋，董立国等．边坡变形破坏过程的大变形有限元分析[J]．岩土力学．2003，24(4)：644～652．
[19] 徐平，甘军，丁秀丽等．三峡工程船闸高边坡岩体长期变形及稳定性有限元分析．长江科学院院报．1999，16(2)：31～34．
[20] 刘文平，赵燕明，郑颖人．岩质边坡开挖应力与变形的有限元模拟[J]．后勤工程学院学报．2004，2：45～48．
[21] 刘凯欣，高凌天．离散元法研究的评述[J]．力学进展．2003，33(4)：483～490．
[22] 俞良群，邢纪波，徐建等．岩质高陡边坡失稳的离散元法与非连续变形法的对比分析[J]．烟台大学学报(自然科学与工程版)．1999，12(1)：63～66．
[23] 强天驰．不连续变形分析方法(DDA)在边坡稳定分析中的应用[J]．岩土工程界．2004，7(6)：40～42．
[24] 王书法，李树忱，李术才等．节理岩质边坡变形的DDA模拟[J]．岩土力学．2002，23(3)：352～354．
[25] 吴建宏，大西有三，石根华等．三维非连续变形分析(3D DDA)理论及其在岩石边坡失稳数值仿真中的应用[J]岩石力学与工程学报．2003，22(6)：937-942．
[26] Chert, G & Ohnishi, Y. (1999). Slope stability using discontinuous deformation analysis method. In Amadei et al. (eds.), Proceedings of the 37th U.S. Rock Mechanics Symposium, Vail. A.A. Balkema, Rotterdam, 535-541.
[27] Amadei, B., Lin, C.T., Sture, S. & Jung, J. (1994). Modeling fracturing of rock masses with the DDA method. In Nelson & Laubach (eds.), Proceedings of the 1st North American Rock Mechanics Symposium, Austin. A. A. Balkema, Rotterdam, 583-590.
[28] 杨天鸿，唐春安，郑雨天等．FLAC程序在抚顺西露天矿边坡变形治理工程中的应用[J]．地质灾害与环境保护．1999，10(3)：6～11．
[29] 虢曙安，贺春宁，芮勇勤．开挖边坡变形滑动的FLAC数值模拟与分析[J]．湖南交通科技．2003，29(3)：4～5，15．
[30] 罗海宁，焦玉勇．对三维离散单元法中块体接触判断算法的改进[J]．岩土力学．1999，20(2)：37～40，45．
[31] 王泳嘉，刘连峰．三维离散单元法软件系统TRUDEC的研制[J]．岩石力学与工程学报．1996，15(3)：201～210．
[32] 胥建龙，唐志平．离散元与有限元结合的多尺度方法及其应用[J]．计算物理．2003，20(6)：477～482．
[33] 金峰，贾伟伟，王光纶．离散元—边界元动力耦合模型[J]．水利学报．2001，1：23～27．
[34] 李世海，汪远年．三维离散元计算参数选取方法研究[J]．岩石力学与工程学报．2004，23(21)：3642～3651．
[35] 李本平．烟囱定向爆破仿真模拟[J]．爆破．2003，20(3)：11～13．
[36] 方新秋，万德钧，王庆等．离散元法在分析综放采场矿压中的应用[J]．湘潭矿业学院学报．2003，18(4)：11～14．
[37] 徐泳，李艳洁，李红艳．离散元法在农业机械化中应用评述[J]．2004，5：26～30．
[38] 焦玉勇，葛修润，刘泉声等．三维离散单元法及其在滑坡分析中的应用[J]．岩土工程学报．2000，22(1)：101～104．
[39] 郑文刚，刘凯欣．离散元法工程计算软件的前后处理系统[J]．计算机工程与科学．2000，22(6)：14，15，23．
[40] Ching-Lung Liao、Ta-Peng Chang and Dong-Hwa Young. Elasto-plastic constitutive modeling of granular assembly, Int. J. Solids & Structures. 1995, 32(8-9): 1121-1133.
[41] Ching S. Chang, Sao J. Chao Y. Chang. Estimates of elastic moduli for granular material with anisotropic random packing structure. Int. J. Solids & Structures. 1995, 32(14): 1989-2008.
[42] Ching S. Chang and Jian Gao. Second-gradient constitutive theory for granular material with random packing structure, Int. J. Solids & Structures. 1995, 32(16):2279-2293.
[43] Ta-Peng Chang and Dong-Hwa Young. Stress-strain relationship for granular materials based on the hypothesis of best fit, Int. J. Solids & Structures. 1997, 34(31-32):4087-4100.
[44] F. Emeriault and B. Cambou. Micromechanical modeling of anisotropic non-linear elasticity of granular medium, Int. J. Solids & Structures. 1996, 33(18):2591-2607.
[45] J. bardet and I. Vardoulakis. The asymmetry of stress in granular media, Int. J. Solids & Structures. 2001, 38:353-367.
[46] 黄松元．《散体力学》，机械工业出版社，1992．
[47] 赵彭年．《松散介质力学》，地震出版社，1995．
[48] F. Kun and H.J. Herrmarm: Transition from Damage to Fragmentation in Collision of Solids,Physical Review E, 59 (3/1999), p. 2623-2632.
[49] Kun, F, Herrmann, H. J. A study of fragmentation process using a discrete element method[J]. Computer Methods in Applied Mechanics and Engineering, 1996, 138:3-18.
[50] 邢纪波．梁-颗粒模型导论．北京：地震出版社，1999．
[51] G.A. D'Addetta, F. Kun, E. Ramm and H.J. Herrmann. From Solids to Granulates-Discrete Element Simulations of Fracture and Fragmentation Processes in Geomaterials, in P.A. Vermeer etal.(Eds.), Continuous and Discontinuous Modeling of Cohesive-Frictional Materials, Lecture Notes inPhysics, Springer (2001), p. 231-258.
[52] Cundall, P. A. and Strack, O. D. L. A Discrete Numerical Model for Granular Assemblies[J]. Geotechnique, 1979 (29): 47-65.
[53] Schinner, A. Fast algorithms for the simulation of polygonal particles[J]. Granular Matter, 1999, 2: 35～43.
[54] 谢祚水，王自力，吴剑国．计算结构力学[M]．武汉：华中科技大学出版社，2004．
[55] 长江勘测规划设计研究院，长江科学院等．长江三峡工程永久船闸边坡变形与稳定分析及预报研究．武汉：长江勘测规划设计研究院，2000．