数值流形方法研究及其在岩土工程中的应用
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
数值流形方法是一种具有一般化意义的数值方法,能对连续变形与非连续变形问题进行统一求解,对问题具有较高的求解精度。同时,数值流形方法以覆盖技术为基础,为工程问题提供了一种更自然和简便的解决方法。数值流形方法是目前岩土工程数值分析方法的一个研究热点。
    本文针对流形方法及其在岩土工程中的应用,并结合西部交通建设科技项目“红层软岩填石路堤不均匀沉降机理研究” 开展了对流形方法的研究工作,开展的工作及取得的研究成果主要有以下几个方面:。 
    (1)研究了材料裂隙及分层界线对流形元覆盖系统形成的影响,指出流形元数学网格必须与材料物理网格(裂隙及分层界线)相协调来构造流形元的覆盖系统,并建议了按此原理形成流形元覆盖系统的基本过程;
    (2)基于流形单元上位移函数的组成,提出了流形方法中固定约束处理的新方法,在流形单元上构造能反映边界位移条件的位移函数使边界的固定约束条件得到自动满足,推导了相应的流形单元刚度矩阵及荷载矩阵的数值计算格式;
    (3)基于有限覆盖技术,建立了全长粘结杆件的数值流形方法模型。该模型是一种新的杆件数值模型,它采用局部解析的位移模式,有较高的计算精度。该模型中,规则的网格可以适应杆件的复杂布局,没有数值模拟中杆件布置与单元划分相互制约的矛盾,且可以考虑杆件中施加预应力的情况。该模型可用于岩土锚杆和钢筋混凝土问题的模拟计算; 
    (4)系统的提出了数值流形方法对开挖卸荷问题的模拟方法。提出了流形元模拟开挖时覆盖系统的处理方法、不平衡力计算方法、使用高阶覆盖函数时对开挖的模拟方法及流形元对开挖模拟的一般过程,提出了开挖作用的直接分析法和间接分析法的概念。流形方法对岩土工程开挖卸荷问题的模拟比传统数值方法有更多的优点:①覆盖技术的使用,使简单的数学覆盖可以适应任意的开挖过程;②开挖作用的直接分析法使流形方法对开挖的模拟可以放弃开挖释放荷载的概念,使计算变得更自然、简便和可靠;③高阶的流形方法对开挖问题的模拟具更好的求解精度;
    (5)针对“红层软岩填石路堤不均匀沉降机理研究”,推导了流形元与DDA的耦合计算方程,并提出了填石路堤工后变形的简化计算方法,首次使用完全非连续散体系统模型对填石路堤问题进行数值模拟分析,在模拟方法的研究上取得了一些有意义的成果,在一定程度上揭示了填石路堤变形的细观机理及规律。数值模拟的结果显示,填石路堤的不连续变形与不均匀变形明显,这与将填石体视
The numerical manifold method is a generalization numerical method, it is capable of uniformly dealing with the problem of continuous deformation and discontinuous deformation, and it has higher calculation precision for problem. At the same time, The numerical manifold method is based on finite cover, it is a simple and convenient resolvent for engineering problem. The numerical manifold method is pop method of research for the numerical method in geotechnical engineering.
    Aim at the theory and application of numerical manifold method, and combine the west traffic construction scientific and technical project “study on non-uniformed settlement mechanism of rock-fill embankment of red-bed”, have researched the theory and application of numerical manifold method, works and research findings mainly has following aspects mainly:
    (1) Have research the effects of crack and boundary line in rock mass to the generation of cover system for numerical manifold method, put forward that mathematical cover must suit with crack and boundary line in rock mass, and put forward corresponding method that form the cover system for numerical manifold method;
    (2) Base on the composition of displacement function on manifold element, a new method of constraints treatment of fixed boundary is proposed, fixation constraint terms is strictly met by using corresponding cover function on the physical cover made up of the manifold element, deduce the corresponding numerical calculation format of manifold element’s stiffness matrix and load matrix with fixation constraint terms;
    (3) Created numerical manifold model of fully bonded bar based on the limited cover technology of manifold method, this model is a new numerical model of bolt which use partially analytic displacement pattern. When use high-order cover function, it have feature of partial analytic solution and better calculation accuracy. At the same time, bolt and mathematics cover relatively independent, general grid can meet complex arrange of bolt, have eliminated the contradiction of mutual restriction between bolt disposal and grid partition in numerical simulation of bolt, and this model can consider to the case that exert prestressing force in anchor rod. This model can be used in simulating geotechnics bolt and analyzing the interaction of reinforcing bar and concrete;
    (4) Have research the simulation of excavation with numerical manifold method.
    Have research the treatment method of cover system and calculation method of unbalance force while simulate excavation with manifold method, and have research the simulation method for excavation while use high-order cover function, finally have made general course of simulation of excavation with manifold method. Have made the direct analysis and indirect analysis of excavation action, and have discussed their relation in detail. The manifold method have more merits than FEM while apply in excavation: ①based on cover technology, simple grid can meet the arbitrary course of excavation; ②the direct analysis of excavation can give up the concept of released load, the simulation course is simpler and more natural and more reliable; ③high-order manifold have better precision for excavation; (5) Aim at “study on non-uniformed settlement mechanism of rock-fill embankment of red-bed”, research the coupling problem of numerical manifold method and DDA, put forward the simplified calculation method of the settlement after construction of full-rock embankment, used this completely discontinuous system study for rock fill problem for the first time, have disclosure microscopic the mechanism and law of displacement of fill-rock embankment in some degree. The simulation result display that the uneven and discontinuity of embankment displacement is obvious, it have greatly difference from condition considered rock fill as equivalent continuum. Have carried out numerical simulated analysis for the sensitivity of various factors that affects deformation of embankment, the reduction of friction coefficient between rock block and reduction of deformation parameter of rock have greater influence for deformation of embankment after construction, but the deformation form and law caused by them is different, the influence of reduction of rock block is very little. (6) Programmed the Object-Oriented program with C++ language in VC++ platform, it is capable of calculating the common problem of manifold method, fully bonded bar and excavation etc. The program have better preprocess and postprocess function. And the program have been proved correct by example. At the same time, programmed calculation program EDAnalysis according to “study on non-uniformed settlement mechanism of rock-fill embankment of red-bed”, EDAnalysis have better preprocess and postprocess and the function of dynamic simulation.
引文
[1] Zienkiewicz O C. The finite element method. London Mc Graw-Hill, 1977.
    [2] 朱伯芳著. 有限单元法原理与应用[M]. 水利电力出版社,1979
    [3] 张有天,周维垣. 岩石高边坡的变形与稳定[M]. 北京:中国水利水电出版社,1999
    [4] 朱合华,陈清军,杨林德. 边界元法及其在岩土工程中的应用[J]. 上海: 同济大学出版社,1997
    [5] 王泳嘉. 拉格朗日元法及其在岩土力学中的应用. 第二届全国青年岩土力学与工程学术会议论文集(岩土力学与工程),大连:大连大连理工大学出版社,1995,22-33
    [6] 王泳嘉. 拉格朗日元法-一种分析非线性大变形的数值方法. 第三届华东地区岩土力学学术讨论会暨二十一世纪的岩土力学专题讨论会论文集,武汉华中理工大学出版社,1995:15-27
    [7] 寇晓东,周维垣,杨若琼. FLAC-3D 进行三峡船闸高边坡稳定分析[J]. 岩石力学与工程学报,2001,20(1):6-10
    [8] Cundall P A. A computer model for simulating progressive large-scale movements in blocky rock systems. Proceeding of the International Symposium on Rock Fracture. Nancy France, 1971
    [9] 黄润秋. 边坡治理工程的数值模拟研究[J]. 地质灾害与环境保护,1996,7(1):69-76
    [10] Kawai T. New discrete structural models and generalization of the method of limit analysis. Finite Elements in Nonlinear Mechanics. NITTronheim, 1977, 885-909
    [11] 钱令希,张雄. 结构分析中的刚性有限元法. 计算结构力学及其应用,1991,8(1):1-14
    [12] 钱令希,张雄. 刚性有限元的参变量变分原理及有限元参数二次规划解. 计算结构力学及其应用,1992,9(2):117-123
    [13] 张雄. 刚性有限元的数学理论基础及其在岩土工程中的应用[D]. 大连:大连理工大学,1992
    [14] 石根华著,任放等译,块体系统不连续变形数值分析新方法[M]. 北京:科学出版社,1993
    [15] 石根华著,裴觉民译. 数值流形方法与非连续变形分析[M]. 北京:清华大学出版社,1997
    [16] 周少怀,杨家岭. DDA 数值方法及工程应用研究[J]. 岩土力学,2000,21(2):123-125
    [17] 邬爱清,任放,董学晟. DDA 数值模型及其在岩体工程上的初步应用研究[J]. 岩石力学与工程学报,1997,19(5):411-417
    [18] 吴洪词. 用改进的DDA 法模拟岩体的破裂[J]. 岩石力学与工程学报,1996,15(增):556-558
    [19] 姜清辉. 三维非连续变形分析方法的研究[ D]. 武汉:中国科学院武汉岩土力学研究所,2000
    [20] 郑榕明,张勇慧,王可均. 耦合算法原理与DDA 的耦合[J]. 岩土工程学报,2000,22(6):727-730
    [21] 王书法,朱维申. 考虑空间影响的两种非连续变形分析方法[J]. 岩石力学与工程学报,2000,19(3):369-372
    [22] 姜清辉,丰定祥. 三维非连续变形分析方法中角-面接触模型的研究[J]. 岩石力学与工程学报,2000,19(增):930-935
    [23] 姜清辉,丰定祥. 三维非连续变形分析方法中的锚杆模拟[J]. 岩土力学,2001,22(2):176-178
    [24] 殷坤龙,姜清辉,汪洋. 新滩滑坡运动全过程的非连续变形分析与仿真模拟[J]. 岩石力学与工程学报,2002,21(7):959-962
    [25] 栾茂田,黎勇,林皋. 非连续变形计算力学模型及其在有缝重力坝静力分析中的应用[J].水利学报,2001,4:40-46
    [26] 黎勇,冯夏庭,栾茂田,王泳嘉. 非连续变形计算力学模型中的广义有限单元[J].东北大学学报(自然科学报),2002,23(10):1004-1007
    [27] 黎勇,栾茂田. 非连续变形计算力学模型基本原理及其线性规划解[J]. 大连理工大学学报,2000,40(3):351-357
    [28] 章青,卓家寿. 加锚岩体的界面应力元模型[J]. 岩土工程学报,1 998,2 0 (5):50-53
    [29] 章青,卓家寿. 三峡船闸高边坡稳定分析的界面元法与评判标准[J]. 岩石力学与工程学报,2000,19(3):285-288
    [30] 邓聚龙.灰色系统基本方法[M].武汉:华中理工大学出版社,1987
    [31] 陈新民,罗国煜. 基于经验的边坡稳定性灰色系统分析与[J]. 岩土工程学报,1999,21(5):638-641
    [32] 卢才金,胡厚田,徐建平等. 改进的BP 网络在岩质边坡稳定性评判中的应用[J]. 岩石力学与工程学报,1999(3): 303-307
    [33] 肖专文,张奇志.遗传进货算法在边坡稳定分析中的应用[J].岩土工程学报,1998,20(1):44-46
    [34] 崔政权,李宁. 边坡工程:理论与实践最新发展[M]. 北京:中国水利水电出版社,1999
    [35] Shi G.H. Manifold method of material analysis .Transaction of the Ninth Army conference on applied mathematics and computing .Minneapolis .Minnesota.USA.1991
    [36] Shi G.H. Modeling Rock Joints and Blocks by Manifold Method. Proceedings of 32nd U.S. Symposium on Rock Mechanics. 1992
    [37] 周维垣,杨若琼,剡公瑞. 流形元法及其在工程中的应用[J]. 岩石力学与工程学报,1996,15(3):211-218
    [38] Lin Dezhang , Mo Haihong. Manifold method of materiall analysis. The University Of Okahoma. 1994
    [39] Hideomi Ohtsubo , Katsuyuki Suzuki , Kenjiro Terada And Katsuyoshi Nakanishi. Utilization of Finite Covers in the Manifold Method for Accuracy Control. Proceedings of the 2nd International Conference on Analysis of Discontinuous Deformation. 1997, Kyoto , Japan: 317-322
    [40] 彭华. 岩石边坡稳定分析的数值流形法研究[D].武汉大学博士学位论文,2000
    [41] 蔡永昌. 流形方法的理论与应用研究[D].重庆大学博士学位论文,2000
    [42] 蔡永昌,张湘伟. 流形方法的矩形覆盖系统及其全自动生成算法[J] . 重庆大学学报(自然科学版),2001,1:42-45
    [43] Shyu K.and M.R.Salami.Manifold With Four-node isoparametric finite element method. Work Forum on Manifold method of material. 1995,California. USA.Vol.1:61-68
    [44] Yaw-Jeng Chiou, Ren-Jow Tsay and Wailin Chuang. Crack Propagation Using Manifold Method. Proceedings of the 2nd International Conference on Analysis of Discontinuous Deformation. 1997, Kyoto , Japan: 298-308
    [45] Guangqi Chen, Yuzo Ohnishi and Takahiro Ito. Development of High order Manifold Method. Proceedings of the 2nd International Conference on Analysis of Discontinuous Deformation. 1997, Kyoto , Japan: 132-154
    [46] Lin J.S. and Lee D. H. Manifold method using polynomial basis function of any order.Proc of the 1st Int. Forum on DDA Simulation of Discontinuous Media.1996, Berkeley, California. USA.: 365-372
    [47] Lu, M. Complete N-Order Polynomial Cover Function for Numerical Manifold Method. SINTEF report,2001
    [48] Y. M. Cheng, Y.H.Zhang, W.S.Chen. Wilson non-conforming element in numerical manifold method .COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING. 2002, 18: 877-884
    [49] 曹文贵,速宝玉. 流形元覆盖系统自动形成方法之研究[J]. 岩土工程学报,2001,23(2):187-190
    [50] Ming Lu. High-order manifold method with simplex integration. Fifth International Conference on Analysis of Discontinuous Deformation,ICADD-5 Oct 6-10, 2002 Beer Sheva, Israel, 75-83
    [51] Guangqi Chen, Yuzo Ohnishi, Takahiro Ito. Development of high-order manifold method. International Journal for Numerical methods in Engineering. 1998, Vol.43,no.4, 685-712
    [52] 陈刚,刘佑荣.流形元覆盖系统的有向图遍历生成算法研究[J].岩石力学与工程学报,2003,(5):711-716
    [53] Lin C TAmadei B and Sture S. Using an augmented Lagrangian method and block fracturing in the DDA method. Computer Methods and Advances in Geomech.1994, 837-842
    [54] 王水林. 数值流形方法与裂纹扩展的模拟[D].中国科学院博士学位论文,1998
    [55] 蔡永昌,张湘伟. 使用高阶覆盖位移函数的数值流形方法及其应力精度的改善[J]. 机械工程学报,2000,36(9):20-25
    [56] 蔡永昌,廖林灿,张湘伟. 高精度四节点四边形流形单元[J]. 应用力学学报,2001,18(2):75-80
    [57] 蔡永昌,张湘伟,骆少明. 数值流形方法在连续体数值分析中的应用[J]. 力学与实践,1999,6:53-54
    [58] 邹建德,章光,闵弘. 流形方法数值结果初步探讨[J]. 土工基础,1998,12(2):45-48
    [59] 何俊,陈倩. 一种高精度八节点流形单元[J]. 汕头大学学报(自然科学版),2001,16(2):68-73
    [60] 居炎飞,章光,王水林. 数值流形方法在P 型自适应分析中的初探[J]. 岩土力学,2001,22(1):88-91
    [61] 章光,王水林,闵弘. 数值流形方法及其应用介绍[J]. 岩土工程界,2000,3(12):44-46
    [62] 周维垣,杨若琼,剡公瑞. 流形元法及其在工程中的应用[J]. 岩石力学与工程学报,1996,15(3):211-218
    [63] Betting BP,Han RP.An objected-orented frame work for internective numerical analysis in a graphical user interface environment. Int J Num Meth Engng, 1996, 39:2945-2971
    [64] 张几,许晶月,阮雪榆.面向对象的有限元程序设计[J]. 计算力学学报,1999,16(2):216-225
    [65] 张湘伟,蔡永昌,廖林灿. 数值流形方法物理覆盖系统的自动剖分[J]. 重庆大学学报(自然科学版),2000,23(1):28-31
    [66] 刘欣,朱德懋,陆明万,张雄. 基于流形覆盖思想的无网格方法的研究[J]. 计算力学学报,1997,16(5):405-410
    [67] 骆少明,蔡永昌,张湘伟. 数值流形方法中的网格重分技术及其应用[J]. 重庆大学学报(自然科学版),2001,24(4):34-41
    [68] 王芝银,王思敬,杨志法. 岩石大变形分析的流形方法[J] . 岩石力学与工程学报,1 997, 1 6 (5):3 99-40 4
    [69] 朱以文,曾又林,陈明祥. 岩石大变形分析的增量流形方法[J] . 岩石力学与工程学报,1999,18(1):1-5
    [70] Chen G, Miki S and Ohnishi Y. Automatic creation of mathematical meshes in manifold method of material analysis. Working Forum on the Manifold Method of Material Analysis. California. USA, 1995
    [71] Sasaki T, Morikawa S,Ishii D and Yuzo O. Elastic-plastic analysis of jointed rock models by Manifold Method. Proceedings of 2nd International Conference Analysis of Discontinuous Deformation. Kyoto, 1997, 309-316
    [72] 王书法,朱维申,李术才,陈胜宏. 岩体弹塑性分析的数值流形方法[J]. 岩石力学与工程学报,2002,21(6):900-904
    [73] 王水林,葛修润. 数值流形方法在模拟裂纹扩展中的应用[J]. 岩石力学与工程学报,1997,16(5):405-410
    [74] Zhang G-X, Sugiuta Y and Hasegawa H. Crack propagation and thermal fracture analysis by Manifold Method. Proceedings of 2nd International Conference Analysis of Discontinuous Deformation. Kyoto, 1997, 282-297
    [75] Tsay Ren-JowChiou Yaw-Jeng and Chuang Wai-lin. Crack growth prediction by manifold method. Journal of Engineering Mechanics, 1999, 125(8): 884-890
    [76] 卓家寿,赵宁. 不连续介质静、动力分析的刚体-弹簧元法[J]. 河海大学学报,1993,21(5): 40-43
    [77] Kawai T.A new discrete model for analysis of solid mechanics problem. Seisan Kendyn, 1977, 29(4) : 208-210
    [78] 曹文贵,速宝玉. 岩体锚固支护的数值流形方法模拟及其应用[J]. 岩土工程学报,2001,5:581-583
    [79] 王书法,朱维申,李术才,陈胜宏. 加锚岩体变形分析的数值流形方法[J]. 岩石力学与工程学报,2002,21(8):1120-1123
    [80] 唐树名. 碎裂结构岩体路堑边坡锚固机理分析及其应用分析[D].重庆大学博士学位论文,2003
    [81] 唐树名,吕常新,邓安福. 预应力锚索加固路堑边坡的数值流形分析[J].公路交通科技,2002,19(6):10-12
    [82] 王芝银,李云鹏. 数值流形方法中的几点改进[J]. 岩土工程学报,1998,20(6):33-36
    [83] 朱爱军,邓安福. 岩体工程数值流形方法固定边界约束处理研究[J].岩石力学与工程学报,已录用,拟2005 年10 月刊出
    [84] 朱爱军,邓安福,唐树名等. 流形元法固定边界约束处理研究[J].计算力学学报,已录用待刊出
    [85] Cook R D. 有限元的概念和应用[M].北京:科学出版社,1989
    [86] 张君禄,陈胜宏. 加锚节理岩体时步自适应弹粘性有限元算法[J]. 水利学报,1997(增):168-175
    [87] 朱维申,李术才,陈卫忠. 节理岩体破坏机理和锚固效应及工程应用[M]. 北京:科学出版社,2002
    [88] 张玉军,朱维申. 三峡工程船闸高边坡锚固方案的平面有限元计算[J]. 岩土工程学报,1997,19(1):70-74
    [89] 吕西林,金国芳,吴晓涵. 钢筋混凝土结构非线性有限元理论与应用[M]. 同济大学出版社,1997
    [90] 潘昌实,隧道力学数值方法[M].中国铁道出版社,1995
    [91] 王书法. 岩体不连续非线性变形分析的数值流形方法研究[D].中科院武汉岩土力学研究所博士学位论文,2000
    [92] 曹文贵,唐学军. 岩石块体数值流形分析网格形成方法之研究[J].土木工程学报,2003,36(2):81-85
    [93] 朱爱军,邓安福,颜昌武. 岩体材料物理网格对流形元覆盖系统形成的影响[J].岩土力学,2004,25(12):1933-1936
    [94] Freitag L A, Ollivier-Gooch C.A cost. Benefit analysis of simplicial mesh improvement techniques as measured by solution efficiency[J]. Int J Comput Geom App,2000,10(4):361-382.
    [95] Cavendish J C, Field D A, Frey W H. Anapproach to automatic three-dimensional finite element mesh generation. Int J Numer Meth Engng,1985,21:329-347.
    [96] 彭宣茂,钱向东.有限元网格自动剖分的分区直接法[J]. 河海大学学报,2000,28(6):94-96
    [97] 杜群贵,黎启柏,罗立峰. 有限元网格自动剖分改进的结点连接法[J]. 华南理工大学学报,1999,27(6):20-26
    [98] REBAY S. Efficient unstructured mesh generation by means of Delaunay triangulation and bowyer-waston algorithm. Joural of computational physics, 1993, 106: 125-138.
    [99] 张国新,王光纶,裴觉民.基于流形方法的结构体破坏分析[J]. 岩石力学与工程学报,2001,20(3):281-287
    [100] 徐芝纶.弹性力学简明教程[M]. 北京:高等教育出版社,1980
    [101] 郑宏,葛修润,谷先荣等.关于岩土工程有限元分析中的若干问题[J].岩土力学,1995,16(3):7-11
    [102] 朱维申,李术才,陈卫忠. 节理岩体破坏机理和锚固效应及工程应用[J]. 北京:科学出版社,2002
    [103] 何满潮,景海河,孙晓明. 软岩工程力学[M]. 北京:科学出版社,2002
    [104] 肖明. 地下洞室施工开挖三维动态过程数值模拟分析[J]. 岩土工程学报,2000,22(4):421-425
    [105] 陈健,盛谦,高锋. 构皮滩工程高边坡开挖数值模拟分析及稳定性评价[J]. 长江科学院院报,1996,13(增):31-34
    [106] 吕凤梧,徐伟,刘建航. 深基坑开挖支护的弹道性有限元数值模拟与分析[J]. 建筑技术,2002,33(2):88-91
    [107] 朱爱军,邓安福,吴军.岩土工程开挖卸荷模拟的数值流形方法[J]. 岩土力学,已录用,待刊出
    [108] 唐学军,王建华,寇新建等. 自适应有限元法与开挖过程模拟分析[J]. 矿冶工程,1999,19(4):1-5
    [109] 王勇,殷宗泽. 有限元计算深开挖中挖方等效荷载的分析[J]. 河海大学学报,1998,26(5):71-74
    [110] 高俊合,赵维炳,李兴文. 深开挖有限元分析中释放荷载模拟-三种常用方法比较及改进的Mana 法[J]. 河海大学学报,1999,27(1):47-52
    [111] 章青. 有限元分析中开挖释放荷载的正确计算[J]. 河海大学学报,1999,27(3):112-115
    [112] 王敏强,许原. 有限元分析中开挖释放荷载计算的讨论[J]. 武汉大学学报,2001,34(1):56-59
    [113] Shi Genhua. Discontinuous deformation analysis: a new numerical model for the statics and dynamoics of deformable block structures. Engineering Computations, 1992(9):157-168
    [114] Lanru Jing, Yue Ma, Zulie Fang. Modeling of fluid flow and solid deformation for fractured rocks with discontinuous deformation analysis(DDA) method. International Journal of Rock Mechanics & Mining Sciences, 2001, 38:343-355
    [115] 梁军,刘汉龙,高玉峰. 堆石蠕变机理分析与颗粒破碎特性研究[J]. 岩土力学,2003,24(3):479-483
    [116] 王勇. 堆石流变的机理及研究方法初探[J]. 岩石力学与工程学报,2000,19(4):526-530
    [117] 卡西奥B.维奥蒂(巴西).堆石坝中堆石变形特性的讨论[J]. 水利水电快报,2000,21(6):15-18
    [118] 刘雄.岩石流变学概论[M].北京:地质出版社,1994
    [119] 邓安福,朱爱军,曾祥勇. 红层软岩填石路堤不均匀沉降机理研究.重庆大学科研报告,2004

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700