特厚板弯曲成形理论与三辊弯卷成形工艺研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
厚度>100mm的特厚板卷焊筒节是制造大型中高压厚壁容器和反应器的关键部件。近年来,在石油、化工、核电、造船等行业市场刚需与“十二五”重大装备国产化方针的合力推动下,重型三辊卷板装备研制和大型厚壁筒节成形工艺控制成为卷板行业研究重点。由于特厚板弯曲机理复杂、弯卷工艺影响因素多,因此,在重型卷板机设计阶段和厚壁筒节弯卷成形阶段存在诸多问题,如特厚板弯矩计算误差大、人工操控依赖度大、弯卷道次多和效率低等。本文在山西省科技重大专项——重型三辊卷板成套设备研发(NO.20111101032)的资助下,以提高重型卷板装备设计水平与厚壁筒节弯卷成形控制水平为目标,采用理论建模、有限元模拟和实验研究相结合的方法,对特厚板三辊弯卷成形工艺进行了系统研究;建立了特厚板弯矩模型、三辊预压弯和滚弯模型,制定了特厚板弯卷成形工艺规划方案,确定了工艺参数算法并提出了厚壁筒节成形质量控制策略。
     论文主要研究内容包括:
     1.阐述了大型厚壁中高压容器和重型卷板机的市场背景和研制现状;对厚壁卷焊筒节和重型三辊机的关键工艺参数、多辊弯卷成形理论模型和板材弯卷成形有限元模拟关键技术及国内外研究现状进行了分析和评述。
     2.以Q235特厚板为研究对象,在厚向不同位置截取圆棒拉伸试样并进行静拉伸试验,得到了特厚板厚向强度和塑性分布规律;建立了厚向强度非均匀分布的特厚板纯弯曲工程模型及回弹模型;研究了特厚板弯曲机理、弯矩影响因素及其影响比重随曲率的变化规律。
     3.在考虑厚向强度分布、材料硬化特性、回弹、板-辊接触点偏移等因素的基础上,建立了特厚板三辊滚弯精确模型,制定了上辊压下量与滚弯半径的工艺曲线及近似公式;分析了厚向强度分布和下辊中心距等对滚弯成形力学参数和几何参数的影响规律。
     4.建立了特厚板三辊预弯成形精确模型,获得了上辊压下量与预弯角工艺曲线及理论计算方法;分析了厚向强度分布等对预弯工艺参数的影响,并优化了预弯工艺;采用中厚板不对称压弯试验验证了所建预弯模型的正确性。
     5.将特厚板单道次弯卷成形流程规划为多个工步和动作步,建立了工作辊速度算法;针对重型卷板机的上辊驱动系统的严重时滞特性,建立了以分段恒速为特点的低阶速度响应算法;采用有限元方法对单道次弯卷成形和回弹过程进行了数值模拟;研究了悬置板材的动态惯性、卷板速度、滚弯道次和上辊变形对厚壁筒节质量的影响。
     6.建立了多道次滚弯成形控制算法;采用四种厚壁筒节的多道次滚弯成形模拟结果验证了新算法的精度和适用性;分析了厚壁筒节多道次滚弯成形质量的影响因素,提出了相应的控制策略。
Cylindrical segment with a thickness of more than100mm is one of the keycomponents of large, medium-high pressure, thick-wall vessels and reactors. Inrecent years, the research and production (R&P) of heavy three-roller bendingequipment and the forming process control of cylindrical segment have becomethe key research points in the field of plate bending, under the promoting ofmarket requirement of oil, chemical, nuclear power, shipping industries, etc. andthe “Twelfth Five Plan” policy on the demostic of heavy machine equipments.Because the mechanism of super-thick plate bending is complexy, and there aremany factors influencing the bend-rolling process, many problems exist in thedesign of heavy bender and the bend-rolling forming of thick-wall cylindricalsegment, such as obvious calculation error of thick plate bending moment,excessive dependence on manual operation, more roll-bending passes and lowproduction efficiency, etc.. Under the financial support by the major scientificproject of shanxi province (Research and Development of heavy three-rollerbending equipment NO.20111101032), in this paper, aims to prompt the designlevel of heavy bending equipment, the bend-rolling forming process ofsuper-thick plate was systematically studied with combination methods oftheoretical modeling, finite element simulation and experimental study. Thebending moment model, prebending and roll-bending models of super-thickplate were build up, planning pattern and parameter algorithm of bend-rollingforming process were proposed, and effective control strategies of improving theforming quality of cylindrical segment were proposed.
     The main contents are as follows:
     1. The marketing background and R&P status of large-and-thick vessel andheavy three-roller bender was discussed. The key parameters of thick-wallcylindrical segment and heavy three-roller bender, the theoretical model ofmulti-roller bend-rolling forming process, and FEM simulation key technologiesand research status at home and abroad of plate bend-rolling forming wereanalyzed and reviewed.
     2. Q235thick plates were selected as research objects, tensile testing ofbars cut from different positions through the plate thickness were carried out,strength and ductility distribution through the thickness of the super-thick platewas obtained. The engineering bending model and springback model ofsuper-thick plate considering nonuniform strength distribution through thicknesswere built up. The bending mechanism of super-thick plate, factors influencingthe bending moment and its variation rules with relative curvature were studied.
     3. The precise three-roller bending model of super-thick plate was set upconsidering strength distribution through thickness, material hardening property,springback, and plate-roller contact shift, etc.. The relationship betweentop-roller displacement and roll-bending radius and the approximation formulawas obtained, and the effect of strength distribution and bottom-rollers distance,etc. on the mechanical and geometrical parameters was analyzed.
     4. The precise model of super-thick plate prebending was proposed. Therelationship between top-roller displacement and prebending angle and thetheoretical model was obtained. The influence of strength distribution throughthickness on the prebending process parameters was analyzed, and theprebending process was optimized. The presented prebending model wasverified by the asymmetrical bending test of moderate-thick plate.
     5. The single-pass bend-rolling forming process of super-thick plate wasplaned into multiple working steps and movement steps, and the velocityalgorithm of the work roller was presented. According to the obvious time-lagproperty of the top roller driving system, the speed control algorithm oflow-order response was presented with a characteristic of section constant speed.The single-pass bend-rolling forming and springback process was simulatedwith finite element method. The effect of plate dynamic inertia, roll-bendingvelocity and top roller deformation on the bending quality of thick-wallcylindrical segment was studied.
     6. The control algorithm of multi-pass roll-bending forming process ofsuper-thick plate was set up. The accuracy and availability of the new algorithmwas verified with the simulation result of multi-pass roll-bending forming of four kinds of thick-wall cylindrical segments. Factors influencing the multi-passroll-bending quality of thick-wall cylindrical segment were analyzed, andcorresponding control strategies were proposed.
引文
[1]李文健,杨松,雷万庆.电站锅炉压力容器行业现状与未来[C].中国机械工程学会焊接学会四十周年、中国焊接协会十五周年纪念文集.2002,11:55-70.
    [2]陈建俊.我国压力容器制造技术发展方向的探讨[C].第六届全国压力容器学术会议压力容器先进技术精选集.2005,9:48-63.
    [3]刑伟荣.卷板机的现状与发展[J].锻压装备与制造技术,2010,2:10-16.
    [4]哈尔滨锅炉厂技术科.国外几家原子能压力容器制造厂的制造能力及其它有关情况[J].化工与通用机械,1973,7:27-40.
    [5]张振戎,张文辉,卢庆春.加氢反应器的发展历史[J].一重技术,2004,1:1-3.
    [6]陈晓玲,李多民,段滋华.加氢反应器的发展现状[J].化工装备技术,2009,30(1):29-32.
    [7]陈国理.新型高压容器的研究与实验[J].石油化工设备技术,1991,12(6):8-12.
    [8]朱国辉,陈志平,郑传祥等.新型绕带式压力容器技术的突出发展优势[J].化工机械,2000,27(3):162-167.
    [9]刘栋,宋鹏云,胡明辅等.整体多层包扎式高压容器应用于加氢反应器的前景[J].机械,2010,27(1):1-4.
    [10]杜敏.会员之窗[J].大型铸锻件,2013,6:49-50.
    [11]丁振君,陈贵卿.国内首台最大壁厚板焊厚壁加氢反应器应用分析[J].中国设备工程,2010,2:48-50.
    [12]陈崇刚.我国加氢反应器制造面临的课题[J].炼油技术与工程,2011,41(7):16-20.
    [13]刘振堂.我国锻压机械行业现状概况[J].锻压装备与制造技术,2011,4:9-15.
    [14]刑伟荣.水平下调式结构在大型三辊卷板机上的应用[J].锻压装备与制造技术,2006,5:20-23.
    [15]李云良,张汉谦,彭碧草等.核电压力容器用钢的发展及研究现状[J].压力容器,2010,27(5):36-42.
    [16]王元清,周晖,胡宗文等.钢结构厚板力学性能的低温试验研究[J].土木建筑与环境工程,2011,33(5):7-12.
    [17] Ludwik P. Technologische Studie uberBlechbiegung[J]. Technische Blatter,1903:133-159.
    [18]余同希,章亮炽.塑性弯曲理论及应用[M].北京:科学出版社,1992.
    [19] Hill R. The Methematical Theory of Plasticity [M]. New York: Oxford University PressOxford.1950.
    [20] Wollter K. Freies Biegen von Blechen [M]. Dusseldorf: VDI-Forschungsh,1952.
    [21] Crafoord R. Plastic sheet bending [D]. Sweden, Gothenburg: Chalmers University ofTechnology,1970.
    [22] Proksa F. Zur Theorie des Plastischen Blechbiegens [D]. Germany: TechnischeHochschule Hannover,1958.
    [23]莫施宁.弯卷机和矫正机[M].北京:机械工业出版社,1958.
    [24]陈宜周.幂强化时无限宽板的纯弯曲[J].力学学报,1962,5(2):107-116.
    [25] Spencer R M, Shapcott R. Bend forming of strain hardending plates [C], AustralianConference on Manufactring Engineering.1977,8:249-254.
    [26] Chuantao Wang, Gary Kinzel, Taylan Altan.Mathematical modeling of plane-strainbending of sheet and plate [J]. Journal of Materials Processing Technology,1993,(39):279-304.
    [27]章亮炽,余同希.宽板弹塑性纯弯曲的一个精华理论[J].北京大学学报(自然科学版),1988,24(1):65-72.
    [28] Tan Z, Persson B, Magnusson C. Plastic bending of anisotropic sheet metal [J].International Journal of Mechanical Sciences,1995,3:405-421.
    [29] Chakrabarty J, Leew B, Chank C. An exact solution for the elastic/plastic bending ofanisophopic sheet metal under conditions of plane strain [J]. International Journal ofMechanical Sciences,2001,43:1871-1880.
    [30]官英平,赵军.板材弹塑性弯曲应力应变中性层位置关系探讨[J].塑性工程学报,2002,9(2):39-41.
    [31]茅云生,侯磊,王呈方.板材弯曲伸长的理论研究[J].武汉理工大学学报,2002,26(5):592-595.
    [32]彭炎荣,江莹,段继承等.理想塑性及线性硬化板材弯曲成形的变薄规律[J].塑性工程学报,2003,10(3):22-25.
    [33]彭炎荣,段继承,李兆飞等.宽板弯曲过程中板厚的变化规律[J].模具技术,2003,12(1):10-12.
    [34]贺广零,罗文波,卢晋福等.宽板弯曲成形过程中板厚变化规律[J].塑性工程学报,2006,13(6):48-51.
    [35]贺广零,李倩妹,洪芳等.宽板弯曲成形过程中板厚变化精化模型[J].塑性工程学报,2007,14(6):205-208.
    [36] Gardiner F J. The Springback of Metals [J]. Transactions of ASME,1957,79:1-9.
    [37] Crandall S H. Discussion of Gardiner’s paper [J]. Transactions of ASME,1957,79:7-8.
    [38] Queener C A, Deangelis R J. Elastic springback and residual stresses in sheet metalformed by bending [J]. Transactions of ASME,1968,61:757-768.
    [39] Johnson W, Yu T X. On the range of applicability of results for the springback of anelastic/perfect plastic rectangular plate after subjecting it to biaxial pure bending-II [J],International Journal of Mechanical Sciences,1981(23):631-637.
    [40] Yu T X, Johnson W. Influence of axial force on the elastic-plastic bending and springbackof a beam [J]. Mechanical Working Technology,1982,6:5-21.
    [41] Leu D K. A simplified approach for evaluating bendabilitu and springback in plasticbendingof anisotropic sheet metals [J]. Journal of Materials Processing Technology,1997,66:9-17.
    [42] Gau J T,.Kinzel G L. A new model for springbacl predition in whick the bauschingereffect is considered [J]. International Journal of Mechanical Sciences,2001,43:1813-1832.
    [43] ZHANG Dong-juan, CUI Zhen-shan, RUAN Xue-yu, et al. An analytical model forpredicting springback and side wall curl of sheet after U-bending [J]. Computational MaterialsSciences,2007,38:707-715.
    [44] Morestin F, Boivin M, SILVA C.Elasto plastic formulation using a kinematic hardeningmodel [J]. Journal of Materials Processing Technology,1996,56:619-630.
    [45] Tekiner Z. An experimental study on the examination of springback of sheet metals withseveral thickness and properties in bending dies [J]. Journal of Materials ProcessingTechnology,2004,145:109-117.
    [46] ZHANG L C, LIN Z. An analytical slution to springback f sheet metals stamped by arigid punch and elastic die [J]. Journal of Materials Processing Technology,1997,63:49-54.
    [47] ASNAFI N. Springback and fracture in V-die air bending of thick stainless steel sheets[J]. Materials and Design,2000,21:217-236.
    [48] Yi H K, Kim D W, Tyne C J V, Moom Y H. Analytical predition of springback based onresidual differential strain during sheet metal bending [J]. Journal of Materials ProcessingTechnology,2008,222:117-129.
    [49] Kim H, Nargundkar N, Altan T. Predition of bend allowance and springback in airbending [J]. AMSE,2007,129:129-351.
    [50] Panthi S K, Ramakrishnan N. Semi analytical modeling of springback in arc bending andeffect of forming load [J]. Transaction of Nonferrous Merals Society of China,2011,21:2276-2284.
    [51] Fenoglietto F. Springback analysis in orthotropic sheet metal forming with anelastoplastic formulation wsing a kinematic hardening model [C]. Proc NUMIFORM5, Shen&Dawson(eds), Balkema, Rotterdam,1995:699-704.
    [52] Zang S L, Liang J, Guo C. A constitutive model for spring-back prediction in which thechange of Young’s modulus with plastic deformation is considered [J]. International Journalof Machine Tools&Manufacture,2007,47:1791-1797.
    [53]胡卫龙,王仲仁.各种卷板成形工艺的辊筒受力分析[J].锻压机械,1992,4:43-45.
    [54]胡卫龙.水平下调式卷板机最佳弯边位置的确定[J].机械工程师,1988,9:23-25.
    [55] Basset B M, Johnson W. The bending of plate using a three roller pyramid type platebending machine [J]. Journal of Strain Analysis.1966,1(5):396-398.
    [56] Hansen, Jannerup. Modeling of elastic bending of beams using a roller bending machine[J]. Trans of ASME,1979,104:304-310.
    [57]曾田,小西.塑性加工,1962,3(8):74.
    [58]曾田,小西.塑性加工,1963,4(31):515.
    [59]山川俊夫.塑性加工,1972,18(93):106.
    [60]徐辅仁.船体板成形后曲率与三辊弯板机中心辊进给量的关系[J].中国造船,1991,1:82-85.
    [61]舒东海.三辊弯板机自动控制模型的建立与实验研究[D].武汉:武汉水运工程学院,1993.
    [62]茅云生.板材弯曲时上辊下降值与成形曲率间关系的理论分析[J].武汉造船,1995,(1):22-27.
    [63] Lysov M I. Teoriya i raschet protsessov izgotovleniya detalei metodami gibki (Theoryand calculation of part production processes by bending)[M]. Moscow: Mashionostroenie,1966.
    [64] Lysov M I, Sosov N V. Formoobrazovanie detalei gibkoi(shaping of parts by bending)
    [M]. Moscow: Mashionostroenie,2001.
    [65] Bodunov N M. Calculation of setup variables for the process of bending and rollingthin-walled components using the finite difference method [J]. Russian Aeromatics,2011,54(1):64-67.
    [66] Hua M. The mechanics of continuous roller bending of plates [D], University of Aston,1986.
    [67] Hua M, Lin Y H. Effect of Strain Hardening on the Continuous Four-roll Plate EdgeBending Process [J]. Journal of Materials Processing Technology,1999,90:12-18.
    [68] Hua M, Lin Y H. Influence of Strain Hardening on the Continuous plate roll-BendingProcess [J]. International Journal of Non-Linear Mechanics,2000,35:883-896.
    [69] Yu T X, Johnson W. The large elastic-plastic deflection of a structure [J]. InternationalJournal of Non-linear Mechanics,1982,17:195-209.
    [70] Hua M, Lin Y H. Large deflection analysis f elastoplastic plate in steady continousfour-roll bending process [J]. Inernational Journal of Mechanical Science,1999,(41):1461-1483.
    [71]胡卫龙.板材弯卷变形的理论分析与实验研究[D].哈尔滨:哈尔滨工业大学,1992.
    [72] Shin Lee, Yim. Mechanics-Based Determination of the Center Roller Displacement inThree-Roll Bending for Smoothly Curved Rectangular Plates [J]. KSME International Journal,2001,15(12):1655-1663.
    [73]舒东海,王呈方.三辊弯板机自动控制数学模型的探讨[J].造船技术,1994,6:19-21.
    [74]茅云生,王呈方,杨明国.重力对板材滚弯成形半径的影响[J].武汉理工大学学报,2002,26(3):315-317.
    [75]龚发云,王训响.非对称式三辊弯板机卷板数学模型的研究[C].十二省区市机械工程学会学术年会.2005:315-317.
    [76] Gandhi A H, Raval H K. Analytical and Empirical Modeling of Top Roller Position forThree-roller Cylindrical Bending of Plates and its Experimental Verification [J]. Journal ofMaterials Processing Technology,2008,197:268-278.
    [77] Yang M, Shime S. Simulation of Pyramitical type Three-roller Bending Process [J].International Journal of Mechanical Sciences,1988,(30-12):877-886.
    [78]杨建国,方洪渊,万鑫.大尺寸板材滚弯过程应力场特点分析[J].材料科学与工艺,2008,16(5):684-689.
    [79] Zeng J, Liu Z, Champliaud H. FEM dynamic simulation and analysis of the roll-beningprocess for forming a conical tube [J]. Journal of Materials Processing Technology,2008,198(1-3):330-343.
    [80] Tailor V K, Gandhi A H, Moliya R D, et al.Finite element analysis of deformed geometryin three-roller plate bending process [J]. Proceedings of the ASME internationalmanufacturing science and engineering conference, MSEC2008,2009,1:609-616.
    [81] Zhao G. Y, Liu Y L, et al. Three-dimenstional finite-elements modeling and simulation ofrotary-draw bending process for thin-walled rectangular tube [J], Materials Science andEngineering A,2009,(499):257-261.
    [82] Feng Zhengkun, Henri Champliaud, Thien-My Dao. Numerical Study ofNon-kinematical Conical Bending with Cylindrical Rolls [J]. Simulation Modelling Practiceand Theory,2009,(17):1710-1722.
    [83] Feng Zhengkun, Henri Champliaud. Modeling and Simulation of AsymmetricalThree-roll Bending Process [J]. Simulation Modelling Practice and Theory,2011,(19):1913-1917.
    [84] Feng Zhengkun, Henri Champliaud. Three-stage Process for Improving Roll BendingQuality [J]. Simulation Modelling Practice and Theory,2011,(19):887-898.
    [85] Han X, Hua L.3D FE modeling of cold rotary forging of aring workpiece [J]. Jounal ofMaterials Processing Technology,2009,(209):5353-5362.
    [86] Ahmed Ktari, zied Antar, Nader Haddar. Modeling and Computation of the Three-rollerBending Process of Steel Sheets [J]. Journal of Mechanical Science and Technology,2012,26(1):123-128.
    [87] Hoggan E J, Scott R I, Barnett M R, et al. Mechanical properties of tension leveled andskin passed steels[J]. Journal of Materials Processing Technology,2002,125-126:155-163.
    [88] Hemmerich E, Rolfe B, Hodgson P D, et al. The effect of pre-strain on the materialbehavior and the Bauschinger effect in the bending of hot rolled and aged steel [J]. MaterialScience and Engineering A,2011,528:3302-3309.
    [89] Weiss M, Rolfe B, Peter D H, et al. Effect of residual stress on the bending of aluminium[J]. Journal of Materials Processing Technology,2012,212:877-883.
    [90]杨宏权,李宜男.国产SA-299特厚板性能试验研究[J].点焊机,2010,40(2):105-108.
    [91] Katashi FUJII, Shinsuke ISHIKAWA, Yasunori NAKAMO, et al. Materials propertiesand residual stresses of rolled thick plate and their effect on static strength [J]. DobokuGakkai Ronbunshuu A,2010,66(2):253-263.
    [92]刘瑞堂,刘文博,刘锦云.工程材料力学性能[M].哈尔滨:哈尔滨工业大学出版社,2001.
    [93] Perduijn A B, Hoogenboom S M. The pure bending of sheet [J]. Journal of MaterialsProcessing Technology,1995,51:274-295.
    [94] Kim Y M, Kim S K, Lim Y J, et al. Effect of microstructure on the yield ratio and lowtemperature toughness of linepipe steels [J]. ISIJ International,2002,42:1571-1577.
    [95] Hua M, Baines K, Cole I M.Continuous four-roll plate bending: A production process forthe manufacture of single seamed tubes of large and medium diameters [J]. InernationalJournal of Machine Tools&Manufacture,1999,(39):905-935.
    [96] Hua M, Sansome D H, Baines K. Mathmatical modeling of the internal bending momentat the top roll contact in muti-pass four-roll thin-plate bending [J]. Journal of MaterialsProcessing Technology.1995,(52):425-459.
    [97] Gandhi A H, Raval H K. Analytical Modeling of Top Roller Position for multiple pass(3-roller) Cylindrical forming of Plates [C]. Proceedings of international MechanicalEngineering Congress and Exposition. Chicago, IL, USA, November5-10, Paper No.IMECE2006-14279.
    [98]陈毓敏.板材与型材弯曲回弹控制原理与方法[M].北京:国防工业出版社,1990:23-27.
    [99] Wang Chuantao, Gary Kinzel, Taylan Altan. Mathematical Modeling of Plane-strainBending of Sheet and Plate [J]. Journal of Materials Processing Technology,1993,39:279-304.
    [100] Li Xuechun, Yang Yuying, Wang Yongzhi, et al. Effect of the Material-hardening Modeon the Springback Simulation Accuracy of V-free Bending [J]. Journal of MaterialsProcessing Technology,2002,123:209-211.
    [101] Narasimhan N, LovellM. Predicting springback in sheetmetal forming—an explicit toimplicit sequential solution procedure [J]. Finite Elements in Analysis and Design,1999,33:29-42.
    [102] Finn M J, Galbraith P C, Wu L, et al. Use of a coupled ex-plicit-implicit solver forcalculating spring-back in automotive body panels [J]. Journal of Materials ProcessingTechnology,1995,50:395-409.
    [103] Bambach M, Hirt G. Performance Assessment of Element Formulations andConstitutive Laws for the Simulation of Incremental Sheet Forming (ISF)[C]. Proceedings ofVIII International Conference on Computational Plasticity COMPLAS,2005. Barcelona:CIMNE (International Center for Numerical Methods in Engineering),2005.
    [104]王勖成.有限单元法[M].北京:国科科大出版社,2003.
    [105]李尚健,董湘怀.金属板材成形过程模拟[M].北京:机械工业出版社,1999.
    [106] Barlat F, Maeda Y, Chung K. Yield function development for aluminum alloy sheets [J].Journal of mechanical physical solids,1997,45(11-12):1727-1763.
    [107] Jackiewicz J. Optimization of the perforamance of prrous low-carbon steels means of anevolution strategy[J]. Materials and Manufacturing Process,2007,22:623-633.
    [108] Kjell Mattiasson, Per Thilderkvis. Simulation of Springback in Sheet Metal Forming [J].Simulation of Materials Processing: Theory, Methods and Application,1995,115-124.
    [109]曹金凤,石亦平.ABAQUS有限元分析常见问题解答[M].北京:机械工业出版社,2009.
    [110] Lee S W,Yang D Y. An Assessment of Numerical Parameters Influencing Springback inExplicit Finite Element Analysis of Sheet Metal Forming Process [J]. Journal of MaterialsProcessing Technology,1998,80-81:60-67.
    [111]付泽民.高强度金属板多道次渐进折弯成形及回弹研究[D].武汉:华中科技大学,2010.

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

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

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