花岗岩残积土研究进展论文

油新华1 何刚1 李晓2

（1.北京城建集团 北京 100025 2.中国科学院地质与地球物理研究所 北京 100029）

摘要 土石混合体边坡是由滑坡堆积、残坡堆积、崩坡堆积、冲洪堆积等形成的一种非均质边坡。由于其物质组成的不均质性，使其与其他边坡具有显著不同的变形性质和破坏机理。针对其物质结构特点，对土石混合体边坡应该采用细观处理技术，将之分别看成均质连续体、等效的均质连续体、非均质不连续体和结构面控制的非连续体4种情况来处理。

关键词 土石混合体边坡 细观处理技术

土石混合体边坡是由土夹碎石或碎块石以及碎石或碎块石夹土等组成的边坡［1，2］。它一般发育在第四系松散堆积层中，主要由滑坡堆积、残坡堆积、崩坡堆积、冲洪堆积等形成。土石混合体边坡在我国有着广泛的分布，如香港、广东、福建等地大面积分布的花岗岩残积土边坡，西南地区、长江三峡库区及黄河中上游广泛分布的古滑坡、崩积层边坡，西藏等地分布的冰碛土边坡等［3］。由于该类边坡通常规模较大、影响因素众多、失稳突发性强、滑移条件复杂，常给国民经济建设、人民正常生活与生命安全带来严重危害和巨大的财产损失。因此，对这类边坡进行系统分析、研究和总结具有特殊的理论和实用价值。

对于土体边坡以及岩体边坡来说，人们已经进行了长期、系统、卓有成效的研究工作，并取得了丰富的实践经验。而对于土石混合体边坡来说，由于物质组成和结构的不同而与其他边坡有着显著不同的变形特点，目前人们对于它的研究还只是处于一种定性分析的阶段，或是通过它与某种因素的相关性分析来探讨滑坡的机理［4，5］，或是通过模型实验来进行稳定性分析［6，7］，而缺少像研究土坡和岩体边坡那样的理论和技术方法。作者认为首先应该认清问题的本质，即土石混合体边坡与众不同的特点，然后以一种新的角度来分析和研究这些问题。

1 细观处理技术的概念

对岩土体的破坏机理研究依据所研究对象的尺寸大小，可分为宏观、细观和微观三类。关于这三类的划分，目前还没有一个统一的标准。按照谢强等人的观点［8］，将野外普遍发育、影响工程岩土体力学特征的断层、节理归为宏观级；将发育在岩土体结构中，用肉眼或显微镜观测，直接影响岩土体力学特征的裂纹归为细观级；将发育在矿物晶体内部，一般对岩土体的宏观力学性质没有直接影响的那些位错归为微观级。从研究实施的角度看，宏观研究多结合野外测绘及现场实验进行。由于观测手段及耗资等原因，所做的研究总是有限的。细观研究多以室内岩样力学实验与显微镜观测的方式进行。以现有的经济技术手段可行性而论，细观水平上岩土体破裂过程的研究是当前主要的研究方法。

传统的岩石力学研究采用的是宏观力学的实验与分析方法，已取得了许多瞩目的成果。随着众多重大岩土工程要求的不断提高和研究工作的逐步深入，人们逐渐认识到，由于岩土体固有的非均质性，仅从唯像学角度以宏观尺寸探讨其力学机理显然是不够的，1998年孙钧提出，应当进一步从细观甚至微观尺度，“三观”相结合，更加全面、深入地来探讨岩石力学问题。

所谓细观尺度，介于微观和宏观之间，它相对于分子、原子的微观尺度足够大，而相对于所涉及对象的宏观尺度又足够小；对于不同的材料及其不同的工程特征尺寸，它是一种相对的概念。对岩石来说，细观尺度可能落到厘米甚至米的范围内。人们认识自然现象的规律，一般总是要经历宏观—细观—微观—宏观的不断反复过程。高一级尺度的研究总是基于低一级尺度的研究之上。这就意味着宏观现象的发生能够在微观研究中得到解释，而微观现象的发生又能够在细观研究中得到解释。

2 细观处理的必要性

一些非均质土，譬如砾质沉积物和崩积物，可以含有包括漂石在内的范围比较广的不同大小的颗粒，从而给岩土工程师提出了一些特殊的问题。对于这些物质的采样和试验的困难已经有过很好的报道［9］，并且得到了一定程度的解决［10］。但是，由于评价这种边坡的稳定性而产生的方法论问题却没有引起足够的重视。常规土的边坡分析方法实际上是建立在级配相对良好的基础之上的，其性状假定为与单元尺寸无关。但是，对于土石混合体边坡来说，由于其物质组成的不均质性，使得有必要从单元尺寸即细观角度上来处理。

土石混合体的力学性质主要受其粒度组成的影响，即受其组成颗粒的大小、形状、均一情况、颗粒的磨圆度特别是孔隙度的影响。如有粘土、粉土或有机质参与时，性质将更加复杂。

图1 土石混合体边坡破坏的三种方式

当土石混合体边坡中含有少量大的颗粒时，一般可能发生穿过颗粒、偏离颗粒和包含颗粒的破碎区变宽3种破坏情况［11］（图1），其破坏机理符合最小抗剪强度原理。哪种现象发生取决于颗粒和周围基质的相对强度、颗粒和破碎区的相对大小、大颗粒的含量以及应力水平。例如穿过颗粒可能在颗粒较弱的情况下发生。如果仅有少数几个大颗粒，则倾向于偏离颗粒破坏，此时破裂区域的几何形状将由于颗粒的存在而在一定程度上略有改变。Hencher ＆ Martin［12］认为当大颗粒含量超过30%时，如香港的崩积物，通过填充物的破裂面一般都要受到大颗粒的干扰。如果存在许多大的颗粒，将会形成锯齿状破裂面（图2）。如果颗粒与边坡尺寸相比太大，或者填充物倾向于沿薄的破裂面破坏，则破碎区变宽的可能性要小于偏离颗粒式的破坏。偏离颗粒破坏在可能发生剪胀的低应力区比限制剪胀的高应力区更可能发生，而在高应力区则可能导致破碎区变宽。

图2 土石混合体边坡的锯齿状破裂面

在上述的讨论中，似乎大颗粒要么增加边坡的稳定性，要么没有大的影响。然而，大颗粒可能降低边坡的稳定性。当沿大颗粒与充填物的交界面或大颗粒之间的剪切强度小于充填物本身的强度时，这种情况就会发生，而且占有相当大的比例。如果大颗粒占据不太理想的位置，则这种作用更加明显。另外还应注意到大颗粒的存在将十分显著地增加坡体的排水性。

从以上的分析可以看出，为了更准确地研究边坡的稳定性，应从细观角度来进行处理。对于重大的工程，应尽可能地进行详细勘察，以确定土石混合体边坡的不均质地质剖面，将其中影响边坡稳定的大的砾石块体找出来，并在稳定性分析中将这些块体考虑进去。

3 细观处理的方法

对于所要研究的具体工程来说，首先要在相对充分的勘察资料基础上建立尽可能反映工程实际情况的岩土工程模型。其中心思想就是首先应该把它看成一种什么样的材料或介质，因此它的内容应包括与边坡有关的基本地质条件，如边坡的尺寸、所含介质的颗粒大小及其均匀性、岩体的风化程度、节理等不连续面的发育程度、岩层的相互组合等。对于不同的工程条件，应该有不同的处理办法。下面是根据土石混合体材料的力学性质研究提出的针对不同边坡工程的处理办法，即根据工程地质条件将边坡材料看成是均质或非均质、连续还是非连续体。

3.1 均质连续体

图3 可看成均质连续体的边坡

最常用的方法是将土石混合体作为一种实质上的均质连续体，它把最主要物质的参数作为区域代表性的参数。这种方法适合于土石混合体中含有少量的砾石块体而且块体的大小和边坡的尺寸相比特别小不足以影响整体力学性质的情况（图3）。这时，将岩土体看成是均质材料是有效的，其工程参数和纯土体相比没有明显的变化，而且不可能发生块体与土体之间结构控制的行为。在这种情况下，没有必要进行野外原位试验，而只需实验室试验即可获得所需的参数。

3.2 等效的均质连续体

当土石混合体中虽然含有较多的砾石，但是砾石分布比较均匀且大小与边坡尺寸相比也比较小时，可以把它看成一种等效的均质连续体（图4）。由前面的分析可知，由于砾石的存在而使得其力学性质有所增强，但是因为试验方法的困难，许多工程师仍将最软弱成分的性质作为区域的代表性质，这种方法一般是保守的。

图4 可等效成均质连续体的边坡

尽管对之进行了整体简化，即将某种物质的性质代表大区域的性质，但是材料对于扰动的敏感性意味着实验室试验得出的数据是不可信的。因此，许多设计者宁愿依赖于原位试验和经验公式。另外值得注意的是，在土石混合体中，砾石与土体之间的接触面在受力变形时变成了不连续面。在许多情况下，不连续面的出现不能被忽略，即使不存在明显的机理，这时不连续面可能占有主要的地位，它对强度的影响是非常大的。

3.3 非均质、非连续体

当土石混合体边坡中含有和边坡相比尺寸较大的砾石块体，即使含量不高但位置较为重要时，则必须要把它看成是非均质、非连续体（图5）。这时，大的砾石块体的存在必然影响着边坡的破坏形态和机理，在进行边坡稳定性分析时，一定要考虑大块体的影响。

图5 非均质、非连续体的土石混合体边坡

在处理非均质岩土体时，由于软弱基质中含有强度较大的核岩（corestone），而使得确定具有实际意义的参数变得十分困难。这时，不可能以岩土体的代表尺寸来测试试样，因此从核岩到较弱的物质的岩土体强度和模量的分布，要么被忽略（假定在没有结构控制时一般是保守的），要么从理论和模拟上来处理［13］。在这方面，许多学者都对之进行了不同程度的研究。Anon以及West等人［11～14］在对较高砾石含量的堆积体边坡分析时，试图计算由不规则破裂面剪胀引起的额外的剪切强度，其几何形状由交错分布的碎屑来控制。Irfan和Tang［15］给出了由于堆积体中强度较高包裹体引起的附加强度的试验性指导原则。只要它们具有相同的特点，相同的方法或许可以适用于不均质岩土体中。Blight［16］以及随后出版的Jaros的工作［17］报道了在预测包含核岩的地基的沉陷时获取实际参数的困难。Jaros将一建筑的低沉陷归结于建筑物下残积土中石英核岩的影响，然而Blight却认为所观察到的沉陷同样可以用稍微不同的参数归因于母岩残积土。DeMello［18］在预测不同压缩性物质区域中建筑物的沉陷时使用相同的方法考虑了同样的问题，他将之归因于不同物质的比例和不同压缩性。

3.4 结构面控制的不连续体

当边坡中含有很多的块石，比如弱风化岩体、碎裂岩体时，其中土体充填物的含量较少，所起的作用不大，起主导作用的是块体与块体之间的节理等不连续面，这时要把它看成由结构面控制的不连续体（图6）。很清楚，在岩土体行为完全被不连续面控制的时候，可以使用正常的岩体力学方法，但是必须注意到风化对沿潜在滑动面剪切强度的影响。

图6 节理面控制行为的土石混合体边坡

由以上的分析可以看出，对于上述四种方法的区分主要是以土石混合体中砾石或块体的含量以及尺寸大小为标准的。但是，其含量多少以及大小如何还没有一个具体的界限，只能根据现场的具体情况和工程经验来确定。作者给出如下建议值：当含石量少于10%时，可以作为第一种情况考虑；当含石量在10%～25%时，属于第二种情况；当含石量为25%～70%时，属第三种情况；当含石量大于70%特别是达到90%时，则可以作为节理岩体来考虑。对于第一和第四种情况，可以使用土力学和岩体力学的方法来处理。而对于第二种情况，可以根据文献中的研究成果来确定土石混合体的强度参数，以进行相应的计算和分析。对于第三种情况，具体方法有两种：一是经过详细的勘察，将大的砾石块体的位置探明出来，利用土石混合体的实测结构模型进行数值模拟，当然这是一种理想方法，要想达到是非常困难的；二是根据局部的统计资料利用蒙特-卡洛方法模拟砾石的分布情况，在随机结构模型的基础上进行数值模拟，这是一种比较实用、切实可行的方法。

致谢 在本文的写作过程中，得到了许多现场和科研机构工作人员的大力协助，在此向他们，特别是原中国地质环境监测院殷跃平副总工程师、中国科学院地质与地球物理研究所张年学研究员、曲永新研究员表示衷心的感谢！

参考文献

［1］郭俊仃，夏季华.超压密土坡稳定分析的非线性有限元法.见：计算机方法在岩石力学及工程中的应用国际学术讨论会论文集.武汉：武汉测绘科技大学出版社，1994

［2］湖南省水利水电勘测设计院.边坡工程地质.北京：水利电力出版社，1983

［3］孙广忠.中国典型滑坡.北京：科学出版社，1988

［4］王发读.浅层堆积物滑坡特征及其与降雨的关系初探.水文地质工程地质，1995，（1）

［5］谢守益，徐卫亚.降雨诱发滑坡机制研究.武汉水利电力大学学报，1999，32（1）

［6］刘光华.河岸松散堆积层滑坡机制及防治对策研究——以重庆市为例.见：三峡库区地质环境暨第二届中日地层环境力学国际学术讨论会文集.北京：煤炭工业出版社，1996

［7］新杰罗·科玛达.大坝水库运行诱发的边坡不稳定性及其产生的机理（内部资料）

［8］谢强等.岩石细观力学试验与分析.成都：西南交通大学出版社，1997

［9］油新华.土石混合体的随机结构模型及其应用研究.北方交通大学博士学位论文，2001

［10］ Watkins D J.Problems in measuring the properties of glacial tills.M.Sc.Thesis.University of London，Imperial College，1969

［11］ West L J，Hencher S R.Assessing the stability of slopes in heterogeneous soils.Landslides Balkema：Rotterdam，1991

［12］ Hencher S R，Martin R P.The description and classification of weathered rocks in Hong Kong for engineering purposes.Proceedings of the 7thSouth East Asia.Geotechnical Conference.Hong Kong.1982，1

［13］ Lindquist E S，Goodman R E.Strength and deformation properties of a physical model mélange.Proc.1st NARM Symp.Austin Balkema Rotterdam，1994

［14］ Anon.Effect of the coarse fractions on the shear strength of colluvium.Geotechical Engineering Office Report 23.Hong Kong，1993

［15］ Irfan R Y.Potential construction problems for bridge foundations in areas underlain by complex marble formation in Hong Kong.Engineering Geology of Construction Geological Society.London Engineering Geology Special Publication，1995

［16］ Bligth R D.Construction in tropical soils.Proceedings of 2ndInternational Conference on Geomechanics in Tropical soils Singapore.Balkema，1988

［17］ Jaros M B.The settlement of two multi-storey buildings on residual Venterdorp lava.South African Institution of Engineers Symposium on Soil-Structure Interaction.Durban，S.Africa，1978

［18］ Demello V F B.Thoughts on soil engineering applicable to residual soils.Proceedings of 3rdSoutheast Asian Conference on Soil Engineering.1972

作为主编或副主编出版学术会议论文集4部；主审全国教材1部；在国内外重要学术刊物或学术会议上发表论文110余篇；其中7篇论文获中国土木工程学会、中国岩石力学与工程学会等有关学会的优秀论文奖，被EI、ISTP等检索系统收录达40余次。 1. 周晶，栾茂田等主编. 中国水利水电工程未来与发展（第九届全国水利水电工程青年学术交流会论文集，大连，2002年9月25-28日）. 大连理工大学出版社，2002年9月，pp.600（120万字).2. 栾茂田主编 （2001）. 岩土力学数值分析与解析方法（第七届全国岩土力学数值分析与解析方法讨论会文集，大连，2001年9月26-28日）. 大连理工大学出版社，2001年9月，pp.564（106万字）.3. 栾茂田主编 （1998）. 土动力学理论与实践（第五届全国土动力学学术会议论文集，大连，1998年6月25-27日）. 大连理工大学出版社，1998年6月，649pp.（100万字）.4. M.T. Luan（栾茂田）and Keizo Ugai (Editor) （1998）. Performance and Evaluation of Soil Slopes under Earthquakes and Rainstorms (Proceedings of the Special Sino-Japanese Forum on Performance and Evaluation of Soil Slopes under Earthquakes and Rainstorms,Dalian,June 28-29）. Dalian University of Technology Press,June 1998，191pp. （30万字）.5. 栾茂田主审 （2001）. 土力学及地基基础 （陈兰云主编）. 机械工业出版社，2001（28.4万字）. 2.1 国际杂志上发表的学术论文2篇（均为第一作者）1. Maotian Luan,Toyoaki Nogami （1997）. Variational Analysis of Earth Pressure on A Rigid Earth-Retaining Wall. Journal of Engineering Mechanics,ASCE,123⑸：524-530. (SCI).2. Maotian Luan and Chongpan Jin （1988）. Dynamic Stability and Deformation of Offshore Gravity Platform Foundation. China Ocean Engineering,2⑴：31-42.2.2 国际与国外会议上发表的学术论文21篇（其中第一作者12篇，非第一作者9篇）1. Maotian Luan （栾茂田） and Dong Wang （2001）. Dynamic Response of Seabed under Wave-Induced Loading. Proceedings of the Fourth International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics,March 26-31,2001,San Diego.2. Xinzheng Wu,Maotian Luan （栾茂田） and Junxia Xin （2001）. Effects of Dynamic Properties of Rockfill Materials on Seismic Response of Concrete-Faced Rockfill Dams. Proceedings of the Fourth International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics,March 26-31,2001,San Diego.3. Fuming Wang,Yigong Ji and Maotian Luan （栾茂田） （2001）. A Simplified Method of Vibration Analysis of Layered Foundation and Applications in Pavement Parameter Identification. Proceedings of the Fourth International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics,March 26-31,2001,San Diego.4. Maotian Luan （栾茂田），Rong Tian and Qing Yang （2001）. A New Numerical Method-Finite-Cover Element-Free Method. Proceedings of the 10th International Conference on Computer Methods and Advances in Geomechanics,January 7-12,2001,Tuscon,Arizona.5. Maotian Luan （栾茂田），Yong Li （2001）. Discontinuous-Nonlinear Deformation Analysis Method and Applications in Geotechnics. Proceedings of the 10th International Conference on Computer Methods and Advances in Geomechanics,January 7-12,2001,Tuscon,Arizona.6. Qing Yang,Jiankui Jiao and Maotian Luan （栾茂田） （2001）. The Experimental Study on Swelling Properties of Swelling Rock. Proceedings of the 10th International Conference on Computer Methods and Advances in Geomechanics,January 7-12,2001,Tuscon,Arizona.7. Qing Yang,Maotian Luan and Jiankui Jiao （2001）. Reliability Analysis of Progressive Failure of Slope and Its Application. Proceedings of the 15th International Conference on Soil Mechanics and Geotechnical Engineering,August 27-31,2001,Istanbul,Turkey.8. Maotian Luan,Tingkai Nian,K.T. Law,C.F. Lee,Keizo Ugai and Qing Yang （1999）. Lower-Bound Solution of Earth Pressure of Cohesive Backfill with Inclined Surface. Proceedings of International Symposium on Slope Stability Engineering: Geotechnical and Geoenvironmental Aspects (IS-Shikoku’99），Matsuyama,Shikoku,Japan,Nov. 8-11,281-286. (ISTP)9. Maotian Luan,Yongzhe Cao and Keizo Ugai （1999）. Shakedown Analysis of Soil Foundations under Varied Loads. Proceedings of International Symposium on Slope Stability Engineering: Geotechnical and Geoenvironmental Aspects (IS-Shikoku’99），Matsuyama,Shikoku,Japan,Nov. 8-11,663-668. (ISTP)10. Qing Yang,Jiankui Jiao,Maotian Luan,Dazhai Shi （1999）. Reliability Analysis and Risk Evaluation of the Slope of Open Pit Mine. Proceedings of International Symposium on Slope Stability Engineering: Geotechnical and Geoenvironmental Aspects (IS-Shikoku’99），Matsuyama,Shikoku,Japan,Nov. 8-11,1095-1100. （ISTP）11. Maotian Luan,Lili Yang,Shangmu Lu,Jiying Yin and Wai-Fah Chen （1999）. Numerical Elasto-Plastic Analyses of Stress and Deformation of Embankment on Soft Soils. Proceedings of International Conference on Enhancement and Promotion of Computational Methods in Engineering and Science (EPMESC Ⅶ），August 1999,Macau.12. Maotian Luan,Chongpan Jin and Gao Lin （1998）. Evaluation Technique for Earthquake-Induced Movements and Dynamic Stability of Embankments. Proceedings of the Special Sino-Japanese Forum on Performance and Evaluation of Soil Slopes under Earthquakes and Rainstorms,Dalian,June 28-29）. Dalian University of Technology Press,June 1998,116-131.13. K.T. Law,C.F. Lee,Maotian Luan and Y. Zhai （1998）. Engineering Behaviour of Loose Fill of Grantic Soil. Proceedings of the Special Sino-Japanese Forum on “Performance and Evaluation of Soil Slopes under Earthquakes and Rainstorms”，Dalian,June 28-29）. Dalian University of Technology Press,June 1998,18-34.14. Maotian Luan and Zhen Zhang （1996）. FEM-Based Formulations of Discontinuous Deformation Analysis Method with Applications in Geomechanics. The Advances in Computational Mechanics (Edited by Wanxie Zhong,Gengdong Cheng and Xikui Li),International Academic Publishers,429-438.15. Max Y. Ma,Austin D.-E.Pan,Maotian Luan （1996）. Seismic Analysis of Stone Arch Bridges Using Discontinuous Deformation Analysis. Proceedings of 11th World Conference on Earthquake Engineering,Elsevier Scientific Ltd.,Paper No.1551,8pp.16. Maotian Luan,Gao Lin and W.F. Chen （1995）. Lumped Parameter Model and Nonlinear DSSI Analysis. Proceedings of the 3rd International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics,St. Louis,April 2-7,1995,355-360.17. Max Y. Ma,Austin D. Pan,Maotian Luan and Joseph M. Gebara （1995）. Stone Arch Bridge Analysis by the DDA Method. Proc. of the 1st International Arch Bridge Conference—Bolton,U.K.,Sept.3-6,1-10pp. (ISTP)18. Austin D. Pan,Joseph M. Gebara and Maotian Luan （1994）. A Finite Element Interpretation of the Method of Discontinuous Deformation Analysis. Proceedings of the Third World Congress on Computational Mechanics.19. Maotian Luan,Gao Lin and Hirokazu Takemiya （1992）. 2-D Computational Model and Numerical Method for Nonlinear Seismic Response of Soil Site. Proc. of the EPMESC'92,Dalian,606-613.20. Maotian Luan,Chongpan Jin,Gao Lin and Hirokazu Takemiya （1990）. Parametric Studies on Dynamic Stability and Permanent Deformation of Cohesive Embankments. Proceedings of the Eighth Japanese Symposium on Earthquake Engineering,Tokyo,975-980.21. Maotian Luan,Guocheng Han （1990）. Dynamic Stability and Sliding Displacements of Rockfill Dams. Proceedings of the Eighth Japanese Symposium on Earthquake Engineering,Tokyo,969-974.2.3 国内杂志上发表的学术论文45篇 （其中第一作者15篇，非第一作者30篇；核心期刊40篇1. 栾茂田，田荣，杨庆 （2002）. 平面广义四节点等参元GQ4及其性能探讨. 力学学报，34⑷：578-585.2. 王栋，栾茂田 （2002）. 波浪作用下粘弹性海床动力响应的数值分析. 工程力学，19⑷：130-134.3. 王栋，栾茂田 （2002）. 广义Biot理论隐－隐式解法的推广及其应用. 世界地震工程，⑵：23-26.4. 黎勇，冯夏廷，栾茂田，王泳嘉 （2002）. 非连续变形计算力学模型中的广义有限单元. 东北大学学报，2002,23⑽：1004-1007.5. 栾茂田，王栋 （2001）. 波浪作用下海床动力反应的数值分析. 海洋工程，19⑷：40-45.6. 栾茂田，黎勇 （2001）. 非连续变形计算力学模型弹塑性分析方法及其应用. 岩土工程学报，23⑸：530-535.7. 姬亦工，王复明，栾茂田 （2002）. 层状粘弹性地基动力响应的样条半解析法. 岩土工程学报，24⑴：34-37.8. 栾茂田，田荣，杨庆 （2002）. 有限覆盖无单元法在岩土类弱拉型材料摩擦接触问题中的应用. 岩土工程学报，24⑵：137-141.9. 栾茂田，辛军霞，赵尚传 （2002）. 面板堆石坝边坡稳定性可靠性分析及其应用. 大连理工大学学报，2002,42⑷：467-471.10. 张大林，栾茂田，杨庆，田荣（2002）. 基于流形方法的动态应力强度因子数值算法. 大连理工大学学报，2002,42⑸：590-593.11. 年廷凯，栾茂田 （2002）. 均布荷载作用下挡土墙后粘性填土的土压力计算. 岩土力学，2002,23⑴：17-22.12. 郭莹，董秀竹，栾茂田 （2002）. 《建筑地基基础设计规范》承载力的确定方法研究. 岩土力学，23⑷：474-477.13. 栾茂田，姬亦工，王复明 （2001）. 非均质粘弹性地基的动力柔度系数. 大连理工大学学报，41⑸： 621-625.14. 王栋，栾茂田，郭莹 （2001）. 波浪作用下海床动力反应有限元数值模拟与液化分析. 大连理工大学学报，41⑵：216-222.15. 王飞，栾茂田，朱菊芬，屠晟 （2001）. 基于弹塑性有限元分析土工安定性计算方法及其应用. 大连理工大学学报，41⑶：349-354. (EI)16. 栾茂田，黎勇，林皋 （2001）. 非连续变形计算力学模型及其在有缝重力坝静力分析中的应用. 水利学报，⑷：40-46.17. 黎勇，栾茂田，林皋 （2001）. 水位循环变化作用下有纵缝重力坝的非线性特性. 水力发电学报，⑵：26-35.18. 栾茂田，吴兴征，阴吉英，辛军霞 （2001）. 堆石料动力特性参数对面板堆石坝三维非线性地震响应的影响. 水力发电学报，⑴：7-18.19. 吴兴征，栾茂田，李相崧 （2001）. 复杂加载路径下堆石料动力本构模型及数值模拟. 世界地震工程，17⑴：9-14.20. 栾茂田，吴兴征，李相崧 （2001）. 堆石料的亚塑性边界面模型及其验证. 岩石力学与工程学报，20⑵：164-170.21. 田荣，栾茂田，杨庆，Keizo Ugai （2001）. 高阶流行方法及其应用. 工程力学，18⑵：21-26.22. 黎勇，栾茂田 （2001）. 有纵缝重力坝的非连续与弹塑性耦合数值模拟. 岩石力学与工程学报，20（增刊2）：1555-1560.23. 杨庆，田荣，栾茂田，张大林 （2001）. 有限覆盖－无单元法及其在孔边应力集中问题中的应用. 岩石力学与工程学报，20（增刊2）：1561-1564.24. 吴兴征，栾茂田，辛军霞 （2001）. 修正邓肯模型及在面板堆石坝应力与变形分析中的应用. 岩石力学与工程学报，20（增刊1）：1098-1102.25. 王飞，栾茂田，朱菊芬 （2001）. 弹塑性有限元法在土工安定性分析中的应用. 岩石力学与工程学报，20（增刊1）：1111-1115.26. 王栋，栾茂田，郭莹 （2001）. 波浪作用下弹性海床的动力反应及其影响因素分析. 岩石力学与工程学报，20（增刊1）：1132-1136.27. 阴吉英，栾茂田，吴兴征 （2001）. 面板堆石坝地震反应特性的比较分析. 岩石力学与工程学报，20（增刊1）：1158-1162.28. 杨庆，焦建奎，栾茂田 （2001）. 膨胀岩土侧限膨胀试验新方法与膨胀本构关系. 岩土工程学报，23⑴：49-52.29. 姬亦工，栾茂田，王复明 （2000). 粘弹性非均质地基动力柔度的样条半解析解法. 振动工程学报，13⑷：535-543.30. 姬亦工，杜朝伟，王复明，栾茂田 （2000). 空间地基有限元分析及反问题. 郑州工业大学学报，21⑵：25-27.31. 栾茂田，罗锦添，李焯芬，翟阳 （2000). 不排水条件下全风化花岗岩残积土的工程特性与本构关系. 大连理工大学学报，大连理工大学学报创刊50周年特刊，40 （特刊）:S82-89,S94.32. 吴兴征，栾茂田，阴吉英 （2000). 面板堆石坝应力与变形的弹塑性有限元分析与计算. 大连理工大学学报，40⑸：602-608.33. 栾茂田，黎勇，杨庆 （2000). 非连续变形计算力学模型其在岩体边坡稳定性分析中的应用. 岩石力学与工程学报，19⑶：289-294. (EI)34. 黎勇，栾茂田 （2000). 非连续变形计算力学模型基本原理及其数学规划解. 大连理工大学学报，40⑶：351-357. (EI)35. 田荣，栾茂田，杨庆 （2000). 高阶形式广义节点有限元法及其应用. 大连理工大学学报，40⑷： 492-495.36. 栾茂田，田荣，杨庆 （2000). 广义节点有限元法. 计算力学学报，17⑵：192-200.37. 吴兴征，栾茂田 （2000). Nelton土工网格CE131拉伸特性的试验研究. 大坝观测与土工测试，⑴.38. 韩国城，栾茂田，刘文廷，肖白云，汪荣大 （2000). 瀑布沟心墙堆石坝地震反应三维分析. 大连理工大学学报，40⑵：218-223.39. 杨庆，焦建奎，栾茂田，师大寨 （2000). 边坡可靠性与经济风险性分析及其应用. 工程地质学报，⑵.40. 栾茂田，Keizo Ugai （1999）. 关于岩土工程研究中若干基本力学问题的思考. 大连理工大学学报，大连理工大学建校50周年特刊，39⑵：309-317.41. 栾茂田，李锦，杨庆 （1999）. 地震和暴雨滑坡中日研讨会学术总结. 国际学术动态，⑶：69-74.42. 杨庆，栾茂田，崇金著，薛永生，庞海臻，杨军 （1999）. 混凝土底板与碎石垫层室内水平滑动试验研究”. 工程勘察，⑹：5-8.43. 杨庆，栾茂田 （1999）. 地下水易污性评价方法－DRASTIC指标体系. 水文地质工程地质，⑵：4-9.44. 杨庆，栾茂田，崇金著，王国利，周慧成，周集体，王栋 （1999）. DRASTIC指标体系法在大连市地下水易污性评价中应用的研究. 大连理工大学学报，39⑸：684-688.45. 栾茂田，周晶，康海贵，宋玉普 （1998）. 浅议博士研究生课程建设和素质培养. 高等教育研究，1998年，⑴：39-44.2.4 国内会议上发表的主要学术论文45篇（其中大会主题报告4篇，第一作者17篇，非第一作者28篇）1. 刘占阁，栾茂田，许成顺，王栋，郭莹. 成层非均质土的一维固结分析.2. 许成顺，栾茂田，刘占阁，郭莹，王栋. 不同应力状态下土的强度参数之间关系探讨.3.孔德森，栾茂田，杨庆. 桩－土－结构动力相互作用分析研究.4. 年廷凯，栾茂田，郑德凤，崇金著. 卫生填埋场压实粘土衬垫渗透性的试验研究.5. 栾茂田,田荣，杨庆，张大林（2002）. 有限覆盖无单元方法在岩质材料宏观裂纹扩展分析中的应用.6. 张大林，栾茂田，杨庆（2002）. 数值流形方法的网格生成技术.7. 娄峰，杨庆，栾茂田，张大林（2002）. 时间序列AR模型定阶准则的改进方法.8. 栾茂田 （2002）. 岩土工程学科发展动态. 山东省岩土工程学科发展战略研讨会 （泰安，2002年7月8-9日）. （大会主题报告）9. 栾茂田，王栋，郭莹，王飞，刘占阁，许成顺，冯晓莹，赵少飞，李波 （2002）. 海床与海洋地基动力分析理论与设计方法研究进展评述. 第六届土动力学学术会议论文集 （南京，2002年5月18-21日），北京：中国建筑工业出版社. （大会主题报告）10. 栾茂田，吴兴征 （2002）. 多向地震动作用下面板堆石坝动力响应非线性数值分析.11. 栾茂田，姬亦工，王复明，冯晓莹 （2002）. 非均质地基上刚性条形基础动力阻抗的样条半解析法.12. 王飞，栾茂田，刘占阁，许成顺 （2002）. 波浪荷载作用下海洋平台地基承载力评价方法及其应用.13. 栾茂田，郭莹，李木国，王静，王栋，崇金著 （2002）. 土工静力－动力液压三轴－扭转多功能剪切仪研制中若干问题探讨.14. 郭莹，董秀竹，栾茂田，王栋，李木国，王静 （2002）. 波浪荷载作用下砂土动力变形特性的试验研究.15. 栾茂田,吴兴征 （2001）. 堆石料静动力耦合本构模型及参数确定与实验验证.16. 吴兴征，栾茂田 （2001）. 面板堆石坝三维耦合非线性数值计算方法与实例分析.17. 翟阳，李焯芬，罗锦添，栾茂田 （2001）. 全风化花岗岩松散填土稳定状态线的归一化分析.18. 吴兴征，栾茂田 （2001）. 关于堆石料加卸载体积模量的探讨.19. 屠晟，栾茂田，王飞 （2001）. 非均质地基的安定性分析与计算.20. 栾茂田，王飞，屠晟，朱菊芬 （2001）. 非均质地基的极限荷载与安定荷载比较分析.21. 年廷凯，栾茂田，杨庆 （2001）. 倾斜表面挡土墙后粘性填土的地震土压力计算.22. 姜韶华，栾茂田，吴兴征，杨庆 （2001）. 河谷岸坡对面板堆石坝三维非线性地震响应的影响.23. 栾茂田，张晨明，王栋，郭莹 （2001）. 波浪作用下缓坡海床弹性动力响应数值分析.24. 贺杰，杨庆，栾茂田，张慧珍 （2001）. 膨胀土抗剪强度的试验研究.25. 张慧珍，杨庆，栾茂田 （2001）. 膨胀力与含水量关系的试验研究.26. 辛军霞，栾茂田，吴兴征 （2001）. 单桩竖向承载力的可靠性分析.27. 栾茂田，张大林，杨庆，田荣 （2001）. 有限覆盖无单元法及其在模拟裂纹扩展中的应用.28. 娄峰，杨庆，栾茂田 （2001）. 应用概率分析改进分维数测量精度的方法.29. 武亚军，栾茂田，卢文阁，杨亚洲 （2001）. 深基坑支护设计方法的探讨.30. 黎勇，栾茂田 （2000). 非连续变形计算力学模型在坝基岩体弹塑性分析中的应用.31. 田荣，栾茂田，杨庆 （2000). 一种新的数值方法－有限覆盖无单元法.32. 栾茂田，曹永哲 （1999）. 基于温度参数法的土工安定分析方法及其应用.33. 栾茂田，年廷凯，杨庆 （1999）. 倾斜表面情况下挡土墙后粘性填土的侧向土压力.34. 田荣，栾茂田，杨庆 （1999）. 连续与非连续变形分析的流形方法及其在土力学中的应用.35. 姬亦工，王复明，栾茂田 （1999）. 非均质粘弹性半空间地基动力柔度系数的样条半解析法.36. 杨庆，栾茂田，焦建奎，师大寨 （1999）. 边坡渐进破坏可靠性分析及其应用.37. 栾茂田，林皋，杨庆 （1998）. 岩土地震工程与土动力学的若干最新进展评述.38. 李锦，栾茂田，郭莹 （1998）. 一种砂土动力弹塑性本构模型的探讨.39. 栾茂田，邵宇，林皋 （1998）. 土体地震反应非线性分析方法比较研究.40. 王栋，栾茂田，文立 （1998）. 波浪作用下海床动力学研究进展.41. 翟阳，K.T.Law，栾茂田，C.F.Lee （1998）. 一种全风化花岗岩残积土稳定状态的不排水三轴压缩试验研究. 第五届全国土动力学学术会议论文集《土动力学理论与实践》，大连理工大学出版社，大连，119-123.42. 郭莹，栾茂田 （1998）. 能量分析方法在土动力学中的应用.43. 栾茂田，罗锦添，李焯芬，翟阳 （1998）. 不排水条件下松散砂质土的本构模型. 中国科协第三届青年学术年会，中国科学技术出版社，北京，431-433.44. 栾茂田，黎勇 （1998）. 基于有限块体模型的下限极限分析方法及其应用. 第六届全国岩土力学数值方法与解析讨论会论文集，广东科技出版社，广州，212-217.45. 栾茂田 （1997）. 关于岩土力学与岩土工程研究中若干问题的思考. 中国科协青年科学家论坛《二十一世纪土木工程学科的发展趋势》，科学出版社，北京，247-260. （大会主题报告）