论文题目：基于复合材料细观力学方法的复杂岩体力学问题研究：综述

关键词： 复杂岩体，复合材料细观力学方法，Eshelby等效夹杂理论，自洽模型，均匀化方法

引言

岩体的复杂性与研究意义： 岩体是含有孔洞、层理、节理、弱面等多种缺陷的复杂地质体，其微观结构中的微裂纹、多相包裹体及矿物晶体形态对岩石的物理力学性能有显著影响。因此，研究这些缺陷在力、温度、湿度的外部因素下的演化规律，对于裂隙岩体、节理岩体等特殊岩体结构具有重要意义。

研究方法的演进与交叉： 早期主要采用弹性力学，随后塑性力学、断裂损伤力学等方法被越来越多地应用。然而，鉴于岩石的复杂性，岩体力学需要发展自身理论，并与其他基础学科进行跨学科交叉协作。复合材料力学正是这样一门在多学科基础上发展起来的科学。

复合材料力学方法的应用缘起与发展： 将复合材料力学方法应用于岩石力学研究的想法始于1956年（Einstein），后有Hashin、Salamon、Eshelby等学者跟进。1979年，中国王景涛教授指出该方法能更真实地反映岩体力学性能。1990年，赵平皋教授认为这是一种有望推动岩石力学发展的新思路。

本文目的： 本文旨在总结复合材料细观力学方法在岩石力学领域的应用，通过分析和讨论相关问题，为相关研究者提供该领域研究现状的更全面认识，并推动该方法在地岩工程问题中的进一步应用。

复合材料力学基本理论及其在岩石力学中的应用

复合材料力学的研究方法分类： 主要分为宏观力学方法和细观力学方法。宏观方法将复合材料视为均质整体，忽略内部组分相互作用；细观力学则旨在建立材料宏观性能、组分性能与微观结构之间的定量关系。

核心细观力学理论： 自1957年Eshelby提出等效夹杂理论后，相继发展出Mori-Tanaka方法、自洽理论、广义自洽理论等较为成熟的细观力学理论。20世纪70年代发展起来的均匀化方法，则用于寻找等效均质材料的有效特性，主要包括代表性体积单元法、多尺度渐进均匀化方法和基于夹杂理论的均匀化方法。

Eshelby等效夹杂理论及其应用

理论核心与优势： Eshelby利用三维弹性问题的格林函数法，求解了椭球夹杂问题的解析解，并分析了夹杂周围的弹性应力场。

在岩石力学中的应用：

宏观性能研究： 学者们利用该理论研究岩石宏观等效弹性模量、流变特性及加固工程（Hu等，2013；Xiao和Huang，2013；Nguyen等，2016）。

细观结构演化： 结合细观损伤力学，研究岩石细观结构的变化（Bai，2016）。

损伤模型构建： Pensee和Zhu等将Eshelby理论与能量法结合，推导了脆性材料的损伤模型。例如，Zhu（2006）通过考虑基体弹性应变和裂纹引起的非弹性应变，得到了代表性单元自由能的一般形式（公式1）。Chen等（2011）利用M-T法和Taylor法建立了热力耦合条件下岩石的细观损伤力学模型。

弹塑性本构关系： Weildle等在宏观损伤力学框架下，利用Eshelby等效夹杂法建立了考虑受损岩石塑性变形的弹塑性本构关系。但Eshelby理论解不适用于夹杂比例大的情况，导致模型模拟结果与后期试验结果存在较大差异。

自洽模型及其应用

基本思想： 在均匀边界条件下，通过自洽模型获得夹杂相中的平均应变，从而求得有效弹性刚度张量。

改进与应用： 针对自洽法在多相（夹杂）复合材料等效模量计算中的缺陷，Kerner（1956）提出了广义自洽模型，其预测结果更为合理。但总体而言，自洽法在岩石力学中应用不多。例如，温建华和周彩英等（2011）基于广义自洽理论研究了岩石材料的损伤效应，建立了反映荷载下岩石微裂纹发展变化的损伤本构模型。Wang等（2015）发展了一个用于粗糙表面弹性接触的自洽模型，以考虑微凸体相互作用。

均匀化方法及其应用

均匀化理论广泛应用于具有细观非均匀连续介质的物理和工程领域，在岩石力学中，如节理岩体、锚固岩体、层状岩体等均有广泛应用。

节理岩体中的应用：

问题与思路： 对于单节理岩体，可直接建立节理模型（如Goodman单元）。但对于多组节理或非连续面问题，逐一建模十分困难。因此，学者们将节理岩体等效为广义复合材料（节理为基体，岩块为夹杂），选取代表性单元，利用宏观-细观应力应变协调条件和细观力学方法，获得节理岩体的宏观等效力学参数和本构模型。

研究成果： 张武等（1987）获得了节理材料和完整岩体的弹性常数及节理产状参数；杨海天等（1991）建立了节理岩体的复合蠕变模型。研究者们还针对台阶状、单向/双向/多向、倾斜、层状等多种节理形式的岩体，利用复合材料细观均匀化方法获得了其宏观等效本构模型和参数。

另一种思路： 假设节理岩体具有细观周期性特征，首先通过渐近均匀化方法建立等效弹性参数与细观结构的关系（公式2、3），并基于此对工程问题进行数值分析。

层状岩体中的应用：

定义与重要性： 指由至少两种不同工程行为与性质的岩性单元组成的互层岩体（如软硬岩交替），具有显著的各向异性和非均质性，确定其力学性质对工程建设中的设计选值和稳定性分析至关重要。

研究方法： 常采用层状复合岩石试件进行试验，其力学参数比单一岩石更复杂。在建立本构模型时，常利用均匀化方法，通过代表性单元（如图3所示，可取作平面问题或空间问题）将不同厚度的岩体等效为均质介质。关键在于根据岩层结构特征（如应变协调原理、多相混合法则）选择合适的混合率，从而建立层状岩体的等效弹性、弹塑性和粘性本构模型。例如，杨春和研究组建立了层状盐岩的宏观等效弹性、弹塑性及非线性蠕变本构模型；刘考学等建立了互层和节理岩体的粘弹性模型。

两步均匀化方法及其应用：

方法介绍： 近年来，两步均匀化法被用于解决非均质材料的多尺度问题。例如，Shen和Shao等将粘土岩视为嵌入孔隙粘土基质中的线弹性矿物颗粒复合体。如图4所示，首先在细观尺度上对粘土固相和球形微孔进行第一次均匀化，得到粘土基质的宏观均匀应力-应变关系；在此基础上，对粘土基质和矿物颗粒进行第二次均匀化，建立粘土岩的细观力学模型，评估其宏观塑性性能。

应用范围： 均匀化方法的适用性与岩体结构特征密切相关。对于由基质包裹颗粒构成的各向同性多孔岩石（如石灰石、粘土岩等），复合材料力学方法（如Mori-Tanaka方法、Hill的增量法、Shao方法等）在均匀化过程中应用效果良好。

用复合材料细观力学原理研究其他复合岩土工程问题

复合路基与复合地基： 这类结构通常由连续分布的土体（基体）和被土体包围的桩或土工合成材料（增强体）组成。由于其非均匀结构，研究其本质特性较为困难。叶观宝等（2002）提出了复合地基的细观力学观点，核心是建立复合地基宏观性能与其组分性能及细观结构之间的定量关系。许多学者基于细观力学均匀化方法预测了桩-土复合地基的等效弹性参数。任回平等（2009）将土工格栅和土体视为广义复合材料，建立了预测其弹性等效本构参数的数值计算模型。

用复合材料力学原理研究含杂质盐岩体的力学问题

问题背景： 含杂质盐岩体除了纯盐岩外，还含有石膏、硬石膏、钙芒硝等其他共生矿物组分（杂质）。这些杂质对盐岩体的整体力学性能有显著影响。

研究模型（两种思路）：

均匀介质模型： 将含杂质盐岩体视为均匀介质，通过连续应力-应变场推导连续介质本构模型。从宏观角度将其视为均匀材料，可采用代表性单元，利用各组分应力-应变关系建立其宏观等效应力-应变关系。这符合复合材料细观力学的研究目的，即建立宏观性能、组分性能和微观结构之间的定量关系。
非均匀介质模型： 将含杂质盐岩体视为具有明显结构特征的多相材料（盐岩基体、杂质增强相、界面相）。考虑到各相力学性能的差异，可采用复合材料的细观有限元方法，考虑不同的微观结构和成分（杂质含量、尺寸、形状、分布及力学性能），研究其内部细观应力场、位移场和整体力学性能。揭示杂质盐岩体微观特征对其性能影响的关键在于如何利用细观力学方法预测和分析其性能，揭示其损伤和破坏的本质。

关于用复合材料细观力学方法研究岩石力学关键问题的讨论

复杂岩体与复合材料在结构组成上的相似性

相似性前提： 运用复合材料力学方法研究岩石力学问题，关键在于岩体与复合材料在结构组成上必须存在相似性。

复合材料视角： 复合材料是由两种或多种组分通过优化微观结构而制成的新材料。

岩体视角： 岩体是一种天然地质材料，其本质是非均质的，由颗粒、孔隙、裂隙和胶结物组成。这种非均质性对其强度和变形有显著影响。这与复合材料“由不同组分构成”的思路相反，但本质一致。

结构对比： 从组成上看，岩石可视为包含基体和缺陷（或微裂纹）的两相材料。此外，常见的含胶结物、杂质颗粒、软弱夹层的岩体，如软硬互层岩体、层状盐岩、含杂质盐岩等，其复杂结构与复合材料的宏观复合（如层压）和微观复合（如颗粒增强）模型高度相似。图5展示了双相复合材料的微观结构模型，它们与柱状节理岩体、复合地基、互层岩体、含杂质盐岩等常见岩土工程对象的微观结构非常相似。

复杂岩体与复合材料在力学分析上的相似性

复合材料力学分析特征： 复合材料力学是一种双尺度力学理论，关注细观结构特征（增强相的体积分数、分布、形状、界面性质等）对宏观性能的影响，旨在建立细观力学与宏观力学之间的定量关系，通过微观组分的应力-应变关系和体积分数，推导出复合材料的宏观应力-应变关系（公式4）。

岩体力学分析特征： 岩体力学的观点认为，岩石类型的差异主要在于造岩矿物体积分数的不同。只要确定了主要造岩矿物的物理力学性质和体积分数，就能通过正确的混合律计算出任何类型复合多矿物岩石的物理力学性能。这一观点与复合材料定量分析中细观力学向宏观力学过渡的关系一致。

其他相似性示例： 有学者比较了复合材料与岩块，指出两者都存在清晰的界面，区别在于界面上的受力与位移连续性不同。锚杆支护系统可视为特殊的复合材料（锚杆为增强材料，岩体为基体）。含杂质盐岩可视为颗粒增强复合材料（杂质为增强相，盐岩为基体），杂质对盐岩层的“增强和抑蠕”作用在本质上与颗粒对基体的作用一致。

讨论与展望

未来研究应重点关注以下几个方面：

界面的力学性能： 岩体内部边界（有胶结与无胶结）的力学行为差异显著。界面过渡区（ITZ）是颗粒-基体类复合材料的关键问题。应加强岩石力学试验和微观试验，理解界面的力学性能和粘结作用，研究界面应力传递、损伤和滑移机制，探索界面性能对宏观和细观力学性能的影响。

宏观破坏准则： 复合材料的破坏模式与均质材料不同，通常发生在增强体、基体或界面中的最薄弱环节。其破坏形式与规律不仅与应力状态有关，还与组分性能和组合情况有关。类似地，复杂岩体的宏观变形和破坏是内部组分微观效应的宏观体现。应从细观结构层面研究其宏观破坏机制，这是建立复杂岩体宏观破坏准则的关键问题。

岩土工程的多相复合问题： 与双相复合材料相比，多相材料具有更广泛的物质组成和微观结构，力学行为更复杂。例如，岩土工程中的锚固岩土体可视为由锚杆（索）、水泥浆和岩土体组成的三相复合材料，这可以成为复合材料细观力学方法的进一步研究对象。

数字图像处理技术与有限元方法的结合： 现有的微观行为检测方法无法提供应力应变的微观信息，而纯有限元模型又难以反映真实的微观特征。未来可将数字图像处理技术用于岩石试样切片，通过扫描电镜数字图像获取真实的二维微观结构，然后导入有限元软件进行网格划分，从而更真实地反映岩体中颗粒的形状和分布，获得更精确的计算结果。

结论

利用复合材料细观力学方法研究复杂岩体的力学问题具有重要的现实意义。本文详细总结了Eshelby夹杂理论、自洽模型，特别是均匀化方法在节理岩体、锚固岩体、层状岩体和复杂盐岩体等组合研究中的应用。在此基础上，作者认为应用该方法研究岩体复合问题时，要求复杂岩体在结构上与复合材料相似，在力学分析上与复合材料一致。下一步可借助或借鉴复合材料力学方法，研究复合岩体界面力学特性、破坏准则、多相岩土复合问题以及数字图像处理技术与有限元方法的结合。

参考文献
(此处列出原文中所有参考文献，格式已整理，内容略)