利用自主研发的水分迁移和冻胀试验装备，研究了不同初始含水率、冷端温度、干密度对砂土水分迁移规律的影响，分析了三因素作用下砂土水平冻胀力和冻胀量的变化规律，确定了冻结锋面位置。研究表明：初始含水率和冷端温度对砂土的水分迁移和冻胀效果影响明显；初始含水率从0%增加至10%，试样含水率峰值增大了5.00倍，水平冻胀力和冻胀量不断增大，冻结锋面位置上移至2.5 cm高度处。冷端温度从?5 ℃降低至?15 ℃，试样含水率峰值增大了4.38倍，水平冻胀力和冻胀量不断增大，冻结锋面位置上移至2.6 cm高度处。干密度对试样水分迁移和冻胀特性影响相对不明显，整体呈现出在较小干密度下，试样含水率、水平冻胀力和冻胀量增幅稍大的趋势，冻结锋面集中在2.2～2.5 cm高度处。针对不同影响因素提出了水平冻胀力和冻胀量预测公式，为认识水汽补给下砂土水分迁移规律及合理预防冻胀提供参考。

煤岩复合承载结构所处的应力边界条件不同，冲击地压在巷道中的显现特征和前兆规律亦不相同。采用高频振动采集及孔内成像三轴动静载试验系统，开展了高静载和动静载耦合作用下煤岩组合体真三轴单面临空试验，分析了煤岩组合体界面处力学特征和强度条件，探究了不同应力边界下煤岩组合体的破坏形态、动力显现特征和声发射信号的演变规律。研究结果表明：（1）受煤岩变形相互制约的影响，交界面处砂岩强度被“弱化”。当界面处煤体裂隙尖端的应力大于“弱化”后砂岩强度时，裂隙将穿过煤岩界面发育至砂岩中，砂岩呈现出屈曲层裂、劈裂成板的破坏形态。（2）高静载作用下，煤岩组合体变形破坏特征和声发射信号具有明显的前兆规律，组合体发生承载失效前煤体局部颗粒弹射动能增大、弹射颗粒块度降低，声发射信号由“高频低能”向“高频高能”转变，组合体的破坏形态以剪切?张拉复合破坏为主。（3）受冲击动载影响，顶底板砂岩夹持作用减弱，煤体裂纹尖端应力得不到有效积聚，裂纹扩展到煤岩交界面时被阻隔，组合体以煤样的张拉破坏为主，声发射信号呈现出“高频高能”的特点，但大多集中在冲击破坏之后，导致组合体动力破坏难以预测。（4）与纯静载作用相比，虽然动静载耦合作用下静载水平较低，但煤岩组合体的抛射质量和抛射碎块的分形维数较大、平均破碎块度较小，说明动载对煤岩组合体的破坏起到了正向激励的作用，静载为煤岩体动力破坏提供了应力和能量条件。

通过干湿循环效应下的数字图像三轴剪切试验、CT扫描试验及扫描电镜SEM（scanning electron microscopy）试验，研究了玄武岩纤维加筋黄土干湿循环过程的三轴剪切力学行为及微细观结构演化机制。结果表明：随干湿循环次数增加，纤维含量较高试样的三轴剪切鼓胀破坏形态转变为剪切带破坏；干湿循环早期阶段，剪切破坏形态随纤维含量增加，由剪切带破坏转变为鼓胀破坏。干湿循环作用和纤维含量对应力?应变曲线的类型及特征无明显影响，均表现为应变硬化型。破坏偏应力随干湿循环次数增加而逐渐减小，但衰减速率逐渐减小；破坏偏应力随纤维含量增大先增加而后减小，呈抛物线变化特征，存在一个最佳纤维含量为0.6%。CT数均值ME值呈现与破坏偏应力相似的变化规律。干湿循环作用下筋-土界面产生一定的开裂和松弛现象，弱化了纤维的加筋效应，但与素黄土相比，纤维加筋黄土的微观结构表现出显著的整体稳定性。构建了干湿循环效应下玄武岩纤维加筋黄土的宏细观损伤变量，其表现出一致的变化规律。

以滇中引水香炉山隧洞为背景工程，通过小尺度物理模型试验，对断层错动模拟过程中的隧洞上覆围岩破裂形态、衬砌破坏形态和裂纹发展、应变分布特征等关键响应特征开展试验分析，深入地研究了走滑断层错动影响下跨活断层铰接隧洞的破坏形式及破坏机制。在铰接隧洞抗错断设计参数作用机制方面，详细研究了衬砌节段长度、衬砌厚度、隧洞直径、隧洞轴线与断层带夹角、隧洞断面形式，衬砌材料力学特性等因素对铰接设计隧洞抗断性能的影响。研究结果表明：（1）当跨活断层隧洞不采用铰接设计时，在断层错动作用下隧洞破坏模式为剪切和弯曲组合破坏，衬砌破坏程度严重，存在剥离脱落现象。隧洞截面呈现椭圆化变形，整体坍塌趋势明显；在当前试验中，无铰接隧洞破坏范围达4Wf (Wf为断层带的宽度)。 （2）当跨活断层隧洞采用铰接设计时，在断层错动作用下隧洞整体变形呈现S形，衬砌结构破坏形式为节段间的转动和错台，衬砌节段相对完整，破坏程度较轻；当前试验中，铰接隧洞破坏范围为2.14Wf，相较于无铰接隧洞减少48%，表明铰接设计可改变活断层错断作用下隧洞的变形破坏模式，并减小隧洞结构的损伤范围。（3）铰接设计条件下，隧洞衬砌结构应变峰值主要分布在断层带内，隧洞易发生屈服破坏；与无铰接隧洞相比，铰接隧洞左右边墙纵向最大拉应变减少了56%，纵向最大压应变减少了68%，进一步表明铰接设计可以有效提高隧洞的抗错断性能。（4）在铰接隧洞抗错断设计参数作用机制方面，研究表明增加衬砌厚度、增加衬砌混凝土强度等级、减小节段长度、减小隧洞直径均可提高铰接隧洞的抗错断性能，隧洞穿越断层带的最佳角度为正交，圆形断面相较于三心圆断面可提升铰接隧洞的抗断能力。综上所述，研究成果可以为跨活断裂隧洞工程抗错断措施提供必要的理论参考和技术支持。

岩石节理剪切变形对岩体工程的安全性和稳定性具有重要影响。为研究常法向应力下岩石节理剪切变形本构关系，采用RDS-200型岩石直剪仪对非规则砂岩节理进行了不同法向应力下的直剪试验。根据岩石节理剪切应力?位移全程曲线形状特征，将其依次划分为峰前压密阶段、线性阶段、屈服阶段和峰后软化阶段；根据剪切应力在峰后软化阶段降低幅度和速率大小，将岩石节理剪切应力?位移曲线划分为3种类型：峰后平台型、峰后缓降型和峰后跌落型。基于岩石节理剪切应力?位移曲线各阶段剪切变形特征，采用分段函数建立了岩石节理剪切变形本构模型。与其他模型相比，新提出的岩石节理剪切变形本构模型对试验数据拟合精度更高，更好地描述了岩石节理剪切应力?位移全程曲线。另外，在通过岩石节理直剪试验由经验公式确定模型参数之后，所提出本构模型可在不同法向应力下实现对不同粗糙度岩石节理剪切应力?位移曲线的预测。研究结果对岩石节理剪切变形的数值模拟和工程估计具有一定的实用价值。

金属装配式基础因施工方便与综合造价低而被越来越多地应用于风积沙地区输电线路基础工程。针对中密风积沙开展全模和半模抗拔模型试验，研究了不同深宽比下基础的上拔荷载?位移特性和抗拔破坏模式。结果表明，抗拔荷载?位移关系曲线随着基础深宽比? 变化而改变。当? =1.0时，曲线呈缓变型；当? =2.0～5.0时，曲线呈软化型。随着? 的增大，抗拔拔承载力随之增大。在加载过程中，地表以基础顶部为中心逐渐形成近似为圆形的隆起区域，隆起程度随着上拔位移增大而不断增加，最终形成整体破坏，此时滑动面贯穿地表；通过半模试验观测的上拔土体滑动面子午线可用直线方程近似描述，且抗拔角随着? 的增大而减小。根据上述试验结果，基于Veesaert滑动面摩擦强度理论和极限平衡原理，提出了金属装配式基础的抗拔承载力改进计算方法，计算结果相较于土重法和剪切法更接近试验值，且离散程度更小。

在低品位、多夹层盐矿地层中建设地下储库，利用水采溶腔（盐穴）底部不溶物沉渣空隙储气有望较大幅度扩大储气库库容，但需要开展沉渣空隙储气可行性和适用性评估。首先，采用不溶物沉渣制备大尺寸CT扫描试样，构建了沉渣空隙三维重构模型，定性分析了沉渣空隙尺寸和空隙连通性，发现沉渣空隙尺寸大且连通性好；其次，采用自主研制的盐矿水采沉渣储气试验装置，利用不溶物小颗粒和夹层垮塌大尺寸块体混合堆积样品进行注气排卤试验，测定注气排卤前后沉渣空隙率，分析了沉渣空隙率变化的原因；最后，建立沉渣空隙体积预测模型，提出了盐腔采动空间和沉渣可排卤水体积的计算公式。通过考虑盐层不溶物小颗粒连续堆积，又考虑夹层垮塌大尺寸块体无序堆积，计算了沉渣堆积体积及可排卤储气体积，计算结果与现场采卤数据吻合较好。该研究可为评价低品位、多夹层盐矿沉渣空隙储气能力提供参考。

为了验证高能激光技术能够辅助隧道掘进机(tunnel boring machine，简称TBM)盘形滚刀（以下简称滚刀）实现高效破岩，以滚刀破岩过程中基本的运动形式——侵岩为研究对象，参照17 inch（432 mm）工程用滚刀，按1:8的比例制备了缩尺比例滚刀压头（以下引述为滚刀压头），开展了以孔孔距（2、3、4、5 mm）、刀孔距（3、4、5、6、7 mm）为变量的激光辅助滚刀压头侵岩L20(4×5)正交试验以及无激光作用的对照组试验。试验研究结果表明：随着刀孔距的增加，滚刀压头的破岩量与岩石破碎块度均增加，而垂直力和比能耗均降低，最优刀孔距为6 mm；随着孔孔距的增加，破岩量、垂直力和岩石侵入难度系数均降低，而岩石破碎块度和比能耗均增加，最优孔孔距为3 mm。与传统滚刀回转滚压破岩方式相比，在激光辅助作用下，可促进岩石张拉裂纹的产生，提高了滚刀侵岩效率。所提出的激光辅助滚刀耦合破岩模式在未来TBM极硬岩层隧道掘进中具有一定的应用前景。

为阐明高地应力地下工程开挖卸荷条件下错动带变形破坏的塑性规律与硬化特征，通过一系列不同围压、不同加卸荷应力路径下的三轴试验，着重研究了不同加卸荷应力路径下的错动带塑性变形规律；基于试验结果，在应力空间进一步探究了等效塑性功、等效塑性应变、塑性体应变等内变量作为硬化参数的应力路径相关性，并提出应力路径无关的硬化参量修正公式。研究结果表明：（1）高地应力下错动带的变形力学特性存在着明显的应力路径效应，在相同卸荷应力路径下，高围压会抑制错动带试样环向塑性变形的发展；在相同初始围压下，卸轴压卸围压应力路径下错动带试样的塑性体积残余应变（6%）明显大于其余卸荷应力路径（2%～4%）；不同应力路径对错动带塑性体积变形的促进作用为卸轴压卸围压应力路径＞恒轴压卸围压应力路径＞增轴压卸围压应力路径；（2）错动带在不同卸荷路径变形破坏过程中，其塑性体应变、等效塑性应变、等效塑性功等内变量均存在一定的应力路径相关性，因而在错动带的弹塑性分析中，直接将以上任何一个状态参量作为其硬化参数并假定其与应力路径无关是不准确的。为此，提出等效塑性功修正方法，发现当修正等效塑性功公式中反映错动带材料特性的参数ns = ?0.4时，修正等效塑性功具有明显的应力路径无关性，更适宜作为高地应力卸荷条件下错动带的硬化参数来描述其卸荷塑性应变硬化特征。揭示的错动带塑性力学特性可进一步深化高地应力卸荷条件下错动带的变形破坏认识，为实际工程中错动带的破坏机制分析和支护控制提供理论依据。

基于既有的三组元局域共振结构，考虑到岩土材料流塑性及弹塑性特征，提出一种由管桩、软体层、桩芯和基体组成的四组元局域共振桩结构，并基于周期性理论，利用频散曲线分析了管桩密度、弹性模量、厚度对四组元结构带隙特征影响。研究发现：相比于无管桩的三组元结构，带管桩的四组元结构更加容易出现完全带隙，且管桩只对带隙上边界产生影响；相比于管桩密度，管桩弹性模量可对四组元结构带隙产生质变影响，影响程度更大；通过四组元结构带隙边界点模态分析，确立了四组元结构带隙计算公式。通过时域分析发现，相比于三组元结构，四组元结构在带隙内振动衰减效果可提升近30%，其内在原因是四组元结构中管桩提高了基体整体刚度，从而有效地预防了基体“击穿”带及泄露波的出现。相关研究可为新型局域共振结构设计及隔振效果改善提供新思路。

经典的Richards入渗控制方程属于偏微分方程，具有强烈的非线性，难以求得解析解。以入渗时间为最小作用量，基于Richards方程建立关于入渗路径的时间泛函，将考虑重力项的非饱和土垂直入渗问题转化为泛函极值问题，并构造等价的Euler-Lagrange方程进行求解。计算结果表明，扩散系数D（?）与概化湿润锋距离具有函数关系，当扩散系数D（?）形式已知时，可求得最优路径下湿润锋处含水率、较远处湿润锋最小含水率、土壤含水率最大熵分布3个问题，并基于最优路径检验了本研究条件下，Boltzmann变换和线性变换求解Richards方程的精度。求解过程未引进新变量化简Richards方程，不改变原方程结构，因此其解具有普遍性，可作为非饱和土力学计算的一个补充。

地震过程中地下结构变形主要受周边土体变形控制。基于该思想，现已提出一种在土?结构有限元模型侧边施加土体变形的地下结构抗震简化分析方法——边界位移法，但缺乏相关试验研究。为探索土体侧边施加推覆位移的试验方法可行性，以自行研制的岩土综合试验模型箱为试验平台，以1:10缩尺的大开车站区间隧道为研究对象，开展了大型土?地下结构系统推覆试验。通过应变、位移和应力的分析，揭示了试验过程中地下结构及周边土体的反应特性。结果表明：由于土体材料的强非线性特征，土体侧边施加的倒三角形变形在传递过程中会产生衰减，地下结构受到剪切变形和挤压变形的耦合作用；中柱与底板交接处是整体结构中的抗震薄弱位置；水平基床系数与土层位移水平及结构侧壁的破坏阶段相关；结构整体刚度大于等代土体，两者的侧向变形比值小于1。随着推覆水平增加比值逐渐增大，通过变形特征能够量化土?结构相互作用，有效填补了土?结构相互作用系数试验研究的空白。试验方法与结论，对地下结构抗震分析及推覆试验可行性研究具有重要指导意义。

为了描述岩石强度的各向异性特征，通过引入材料的微结构张量，考虑单轴抗压强度 随层理倾角? 而变化，提出了基于微结构张量的修正Hoek-Brown屈服准则。开展不同层理倾角（0°、15°、30°、45°、60°、75°和90°）的页岩试样不同围压（0、5、10、15和20 MPa）下的三轴压缩试验，结果显示，页岩峰值强度与层理倾角? 呈“U”形曲线关系，在围压较低时呈“肩”形曲线关系，弹性模量Et和泊松比ν对围压不敏感，其平均值随层理倾角? 的变化关系分别呈倒“U”形、“肩”形。页岩的微观和宏观破坏图像表明，其破坏模式受围压和层理倾角双重控制。试验值与模型预测值的对比表明，修正的Hoek-Brown屈服准则可较好地对页岩强度的各向异性进行模拟。同时，采用Koteshwar千枚岩、Koteshwar板岩、绿泥石片岩等层状结构岩石的单轴和三轴压缩试验数据，进一步验证了提出的修正Hoek-Brown屈服准则对多数具有层理结构的横观各向同性岩石峰值强度模拟具有普遍适用性，并对屈服准则在试验数据较多和较少情况下的峰值强度模拟效果进行了评价。

从水平简谐振动作用下二维土?悬臂式刚性墙计算模型出发，基于波动力学理论，同时考虑土层的竖向应力和竖向位移，对二维场地中悬臂式刚性墙的动力响应特性进行了解析研究。首先对土层振动方程进行变换，得到关于体积应变? 的方程，通过分离变量法进行求解，再回代振动方程得到关于位移运动方程的非齐次方程，结合墙与土层的相互作用条件及远场边界条件得到振动方程定解，进而得到地下刚性墙墙上土压力、墙底剪力及弯矩的更为严格的解析解。将所得解与忽略竖向应力解、忽略竖向位移解进行了对比。研究表明，所得解能多反映出一个共振频率，且当土体泊松比大于0.45时，忽略竖向位移解失去意义。通过参数分析，表明激励频率与土体阻尼因子对墙体动力响应影响较大，考虑的振动模态阶数对墙体动力响应影响较小。

如何合理地描述超固结地层静止侧压力系数 分布规律对岩土工程有重要的理论意义和应用价值。通过超固结地层形成过程的应力发展分析，开展了侧应力松弛试验，研究了超固结粉质黏土水平应力释放规律，建立了残余超固结水平应力 与超固结比OCR的联系，提出了基于水平应力释放的超固结土 计算方法，并用原位土体水平应力KSB试验和室内试验进行了验证。提出的超固结土 计算方法物理意义清晰，并在中高OCR条件下实现了与前人归纳经验公式的统一。结果表明：超固结水平应力 经历回弹与松弛2个释放阶段后稳定于残余值 ，粉质黏土 释放主要发生在回弹阶段；超固结水平应力残余值与初始值之比 随OCR增大而增大，但在OCR达到一定值OCRr后趋于稳定。研究成果可供超固结地层中的挡土结构物土压力计算和数值分析提供参考和借鉴。

江苏近海属南黄海海域，是我国海上风电场最集中的地区，目前约占全国总装机容量的70%～75%，该区覆盖层厚度大，桩长深度范围（一般40～60 m）多为粉砂、粉土、粉质黏土地层，尤其是泥面以下20 m深度为易液化粉砂?粉土地层。而南黄海地区亦是我国的地震多发区，因此，研究了该区海上风电场水平场地的地震液化特性。首先基于南黄海某海上风电场50个机位的钻探资料统计特征建立了地层概化模型；随后通过动三轴试验和共振柱试验标定了土体动力分析参数；然后反演了3条地震波（EL-Centro波、Northridge波、Kobe波），将地表附近地震动峰值加速度PGA分别调整为0.05g、0.10g、0.20g、0.40g并进行场地液化分析，重点分析了地震作用下地层的超孔压比、总沉降、分层沉降等特性。研究发现，该区地层为可液化地层。当PGA=0.05g时，各层监测点的超孔压比均小于1.0，地层总沉降为1 cm左右；当PGA为0.10g和0.20g时，仅表层（12 m内）地层完全液化，地层总沉降分别为10 cm和17 cm；PGA=0.40g输入时，泥面下20 m地层均完全液化，地层总沉降为30 cm左右。不同地震作用下，表层土体沉降占比最大，当PGA为0.10g及0.20g时，表层土体沉降约占总沉降的95%以上。因此，该区海上风机桩基础设计应考虑地震作用下场地的上述液化特性。

块体塌落或滑动是岩质边坡工程的重要破坏形式之一，因此块体稳定性分析是岩质边坡工程中一项极为重要的内容。以浙江仙居县神仙居岩质边坡为研究背景，对块体稳定性分析及可视化进行研究。基于线性回归法和非均匀有理B样条法提出拟合结构面和临空面的新方法。基于坐标投影法提出单滑面和双滑面型块体的稳定系数计算方法，将无人机与坐标投影法结合，并基于Matlab开发出适于岩质边坡工程中平面多面体块体和曲面块体稳定性分析的CPG程序。该程序可实现结构面、临空面及不稳定块体的空间表示及可视化。工程实践表明，该研究可高效、准确地处理落石、崩塌等工程地质灾害问题。程序计算结果与坐标投影块体理论计算结果基本一致，块体稳定性分析方法可靠，开发的CPG程序可行，可大幅提高块体稳定性分析效率，具有实际工程意义。

城市地铁隧道施工不可避免地对围岩产生扰动引起地表沉降，动态预测隧道开挖引起的地表沉降是确保地表建构筑物与隧道施工安全的重要基础。针对隧道施工过程中地表沉降难以准确动态预测的问题，在定义纵向开挖度系数? 的基础上，建立横向地表沉降动态预测模型。该模型能够描述同一监测位置沉降曲线随掌子面推进而不断变化的规律，实现施工现场沉降动态预测。结果表明：在特定约束条件下，该模型能够退化为Peck模型以及随机介质理论预测模型；通过现场施工验证了该动态预测模型的准确性与适用性；基于纵向开挖度系数? 将隧道沿纵向分为强烈、中度和轻度共3个影响段，较好地反映了开挖掌子面在不同位置对同一监测截面的影响程度；通过分析建筑物和隔离桩对地表沉降曲线的影响可知，建筑物与其附近地层呈现出协同变形、共同承载的特征，在隧道一侧添加地质钻隔离桩对该侧地表沉降值的减小程度可达71.9%。研究成果对滇中引水工程现场施工与类似工程具有一定指导和借鉴意义。

为研究双巷布置工作面留巷受重复采动影响围岩的变形破坏规律及稳定控制问题，综合运用数值模拟、理论分析和现场观测等方法，研究双巷间距对重复采动下留巷围岩塑性区空间分布特征及偏应力的影响规律，揭示二次采动对留巷塑性区的诱导扩展机制，提出围岩塑性区二次扩展抑制方法。结果表明：重复采动巷道围岩塑性区在巷道轴向上均可划分为6个区域，整体表现为双巷间距越大，则塑性区破坏深度越小的特征，且其形态由非对称分布向对称分布状态转化，双巷距离增加造成的应力旋转程度和偏应力峰值低是主要原因。一次采动滞后影响稳定区维护距离及周期最长，且为二次采动塑性区叠加扩展的基础，是重复采动巷道围岩稳定的关键控制区域。据此提出通过调整巷道空间位置，达到改善应力环境、调整巷道围岩破坏状态的目的，并建立了分区补强支护相结合的调控方法和控制技术体系。现场实测结果表明，该技术体系能够达到巷道安全稳定目的，效果良好。

目前国内工程界对高填方地基的排水稳定性关注较多，然而针对地下水位较高、排水性能较差的原状地基，由快速填筑引起的不排水稳定性问题更为突出。综述了现行规范中关于边坡不排水稳定性分析的相关条文及国内外常用方法。通过分析总应力摩擦角的应力路径相关性以及计算原状地基典型点位的总应力加载路径，解释了采用三轴固结不排水（consolidated-undrained，简称CU）总应力强度指标进行不排水分析，会高估高填方地基的不排水稳定性，该方法存在理论缺陷和工程隐患。以简单假想边坡模型和某高填方机场实际工程为算例，利用简化Bishop法，选用5种不排水分析指标或模型，分别计算其不排水稳定性安全系数。结果表明：CU总应力法会高估高填方地基的不排水稳定性；其余4种方法计算得到的不排水稳定性安全系数比较接近，相对适用。

孔压静力触探（piezocone penetration test，简称CPTU）常被应用于土层划分、土体分类，其各种土体分类方法有一定的适用性和局限性。引用两组实际工程案例，对5种常用的基于CPTU的土体分类图进行了分析比较。分析表明，已有分类图虽然能够区分黏性土和无黏性土，但对于粉质黏土、粉土和粉砂之间的判别均存在较大误差，且多数分类图采用过渡带分类，导致分类结果的准确性偏低。基于Robertson土类指数分类图，根据实际CPTU测试数据分类结果，对土类指数Ic进行了修正，将其同心圆模型修正为平行线模型，由此提出一种基于CPTU的修正土类指数土体分类方法。根据对实测数据的应用表明，提出的修正土类指数土体分类方法可以更好地满足土体分类，且能够满足粉质黏土、粉土和粉砂之间的判别。

为了研究岩石张拉破坏的过程和机制，首先通过模拟计算，验证了加权反演法在进行巴西劈裂试验矩张量求解方面的优越性；其次基于花岗岩的巴西劈裂试验，进行声发射监测、事件定位及矩张量反演，并分析标准反演法和加权反演法对应矩张量结果的差别；使用聚类分析K均值算法对矩张量结果进行分簇研究。研究结果表明：模拟计算中加权反演法可以降低巴西劈裂试验矩张量反演误差；室内花岗岩巴西劈裂试验中，加权反演法较大程度地优化了声发射事件的剪切成分以及压/拉应力轴的分布，使得反演结果更加合理。针对加权反演结果，所有事件可分成3簇，同簇内事件的破裂类型、应力状态、破裂机制、辐射模式等均相近；不同簇中的事件在时空上交叉分布，试验后期在圆盘东?西方向上具有较为明显张拉特性的事件骤增，于圆盘中心附近聚集并占据主导，最终导致试样的宏观破裂。研究结果解释了花岗岩巴西劈裂试验的破裂过程和破裂机制，为岩石力学行为研究提供进一步指导。

巨石作为一种大体积、大质量的滚石，其崩塌失稳及高速、高能远程运动往往导致沿途建筑物和交通线路的毁灭性灾难。以西藏自治区G318国道K4580典型滑坡为工程背景，基于三维非连续变形分析（three-dimensional discontinuous deformation analysis，简称3D-DDA）方法研究巨石崩塌失稳及运动全过程的特征与现象。分别建立该边坡未滑坡、浅层滑坡后和深层滑坡后3种坡形的巨石崩塌3D-DDA数值模型。采用滚石运动横向偏移经验模型，验证3D-DDA巨石运动模拟的准确性。在此基础上，分析巨石崩塌失稳机制及破坏后沿3种不同坡形边坡的运动轨迹和动能演进等运动特征。结果表明，3D-DDA能够有效模拟巨石崩塌失稳、运动发展、剧烈冲击碰撞直至最终静止等整个动力学过程。巨石崩塌表现为滑动→倾倒?滑动→倾倒→翻转?下落的失稳模式转换；巨石运动表现为碰撞、弹跳、飞跃、滚动、滑动等多种运动形式以及横向偏移、侧向偏转等三维空间运动特征，经过道路并与高架桥发生碰撞，引发巨石灾害。不同坡面几何特征下的巨石运动偏移量、弹跳高度、运动至坡底碰撞时间、最终稳定时间等均随着未滑坡、浅层滑坡后、深层滑坡后3种坡形变化而减小。通过3D-DDA巨石崩塌运动分析，预测巨石运动全过程、影响范围、冲击能量、停积位置等，可为巨石防灾减灾对策或措施制定提供依据。

为研究挡墙变位模式以及墙后填土宽度对无黏性有限土体主动土压力的影响，在不同墙后填土宽度条件下分别开展了挡墙平动（T）模式、绕墙底转动（RB）模式以及绕墙顶转动（RT）模式的离散元模拟。根据离散元模拟结果对主动土压力、墙后土体破坏模式以及应力状态进行了分析。研究结果表明：挡墙变位模式和墙后填土宽度的变化使得土体破坏模式和应力状态发生变化，引起主动土压力大小及分布的差异。T与RB模式滑动土楔中内摩擦角调动值会相对初始值增加，且T模式滑动土楔中会出现小主应力拱。RT模式较为特殊，在填土宽度较小时，其应力状态与T模式相似；在填土宽度较大时，滑动土楔上部会出现内摩擦角调动值相对初始值减小的区域，并出现大主应力拱。

降雨入渗滑坡稳定可靠度分析往往仅考虑了土体参数的空间变异性，而忽略土体初始含水率非均匀性空间分布的影响。基于此，推导了不同降雨工况条件下初始含水率任一分布的Green-Ampt入渗模型，利用变步长复化Simpson法求解初始含水率不同空间分布条件下湿润锋深度与降雨历时的函数关系。结合多维正态累计分布函数Mvncdf，计算整个降雨评估基准期滑坡稳定的时变可靠度，并利用某一滑坡实例进行对比验证。计算结果表明，初始含水率空间分布的不同对降雨条件下滑坡稳定时变可靠度的影响十分明显，且相较于初始含水率为均匀和梯形分布；相同降雨条件下，指数分布时滑坡稳定时变可靠度的下降幅度最小。该改进模型适用于土体初始含水率的任意分布，有利于Green-Ampt模型在滑坡稳定性评价的普适性推广应用。

落石冲击棚洞结构作用过程复杂，缺乏统一的落石冲击力表达式。首先，将落石简化为刚性球体，基于Hertz接触理论，推导得到落石冲击力半正弦算法的理论表达式，考虑落石冲击下棚洞的非弹性特征，根据落石与材料碰撞过程中落石加速度曲线特征，采用函数拟合法推导得到落石法向冲击下其冲击力的理论计算方法；然后，基于ANSYS/LS-DYNA软件建立落石冲击棚洞数值计算模型，研究不同冲击速度下落石冲击棚洞动力特征；最后，与现存常见的多种方法进行对比，得出以下结论：Hertz半正弦法得到的落石冲击力远大于函数拟合法和数值法，而函数拟合法和数值法得到的落石冲击力时程曲线相接近，表明函数拟合法更能反映落石与棚洞接触碰撞动力关系；对比其他计算方法可以得到，Hertz算法适用于分析无能量损失下的弹性碰撞问题，而Logistic算法适用于材料大塑性变形的情况，弹塑性接触理论结果和动力有限元结果存在差异，而采用函数拟合推导的计算方法得到的落石最大冲击力和落石冲击作用时间与动力有限元法更接近，更能反映落石冲击棚洞动力响应特征，推导的落石冲击力计算方法可为工程实践中棚洞防护设计提供理论参考。

为研究天然气水合物的合成及分解特性、储层的开采特性以及开采过程中的出砂防砂等问题，自主研制了一套多维天然气水合物开采模拟试验系统。该试验系统主要由气体注入系统、稳压供液系统、蒸汽或热水注入系统、一维/二维/三维模型、回压控制系统、出口计量系统、数据采集控制处理系统等组成。该试验系统的主要特色如下：（1）将高压玻璃做成的一维反应釜外壁和高清摄像头结合达到可视化的目的，监测中低压力下水合物的合成、分解情况以及开采过程中砂粒的运移情况；（2）采用自主研发的双缸恒速恒压泵控制回压阀，使得整个降压开采过程中反应釜中的压力能够实现高精度的均匀减小，达到分段缓慢降压的目的，避免水合物的二次生成，使产气速度更加均匀。以不同粒径的石英砂、去离子水和纯度为99.9%的甲烷气为原材料，开展了一维相平衡试验和二维降压开采试验，通过以上初步应用验证了该试验系统的可靠性。

为了提高我国南海钙质砂地基的抗液化性能，提出利用微生物诱导碳酸钙沉积（MICP）技术联合纤维加筋技术对钙质砂进行改性处理。通过开展动三轴试验，对比分析了改性前后钙质砂试样的动应变、动孔压、应力−应变滞回曲线以及动弹性模量的发展规律和演化特征，并结合扫描电镜（SEM）试验探究了MICP和纤维加筋技术对钙质砂的联合改性机制。研究结果表明：（1）MICP技术可以明显改善钙质砂试样的抗变形与抗液化性能，相比于未胶结处理试样，仅MICP处理试样的动应变和动孔压分别降低了95.74% 和 92.46%；（2）纤维的掺入进一步提升了MICP的改性效果，相比于仅MICP处理试样，MICP和纤维加筋联合处理试样的动应变和动孔压分别降低了 74.32%和 74.18%；（3）MICP 和纤维加筋技术通过减轻试样在循环荷载作用下的循环活动强度和能量耗散、提高试样的动弹性模量和减小动弹性模量的衰减速率，从而实现试样抗变形与抗液化性能的显著提高；（4）SEM 试验分析结果表明，MICP 与纤维对钙质砂动力特性的改善具有协同作用。纤维的掺入为细菌提供了更多的附着场所，促进了碳酸钙晶体的生成量，该部分碳酸钙不仅增加了颗粒间的胶结强度，同时也将纤维固定在砂颗粒上增强了纤维网的约束作用。

为深入探究两淮矿区典型砂质泥岩劈裂注浆起裂机制，研制了常规三轴劈裂注浆试验装置，开展了类砂质泥岩浆压致裂起裂压力模型试验，基于试验结果分析了岩石强度与应力状态对注浆起裂压力、裂缝扩展形态影响规律，揭示了砂质泥岩劈裂注浆起裂机制。研究表明：起裂压力与岩石抗压强度呈正相关，且岩石抗压强度越高，劈裂路径越复杂；起裂压力对围压的敏感程度远高于轴压，且应力差 Δσ =σ _V −σ _H 越大，裂缝形态越规整；孔压三轴条件下，封闭裸孔段浆压致裂法确定的岩石抗拉强度值约为单轴抗拉强度的 2.5倍。该研究结果可为今后类似岩层劈裂注浆参数设计与施工提供参考。

针对天津地区滨海软黏土，开展了一系列动渗耦合三轴试验，系统分析了长期循环作用下软黏土渗透系数的变化规律。结合扫描电镜（SEM）和压汞试验（MIP）对不同变形阶段软黏土的微观结构特征进行分析，阐述了渗透系数的演化规律及机制。研究表明：在动渗耦合作用下，软黏土渗透系数呈现三阶段变化特征；渗透系数随渗流力呈现先增大后减小的变化规律，随动应力和频率呈现先减小后增大的特征。微观试验表明，渗透系数的演化是微结构孔隙大小、形态、分布特征渐进调整的宏观体现：加载初期，颗粒间的架空超大孔隙（孔隙直径 D >2.5 μm）大幅压缩，孔隙形态变化不大，渗透系数呈线性衰减；加载中期，大颗粒破碎生成的小颗粒密布填充于大颗粒间，生成大量“密小而曲折”的小孔隙（0.05 μm<D<0.1 μm），渗透系数呈现缓慢下降；加载末期，结构致密稳定，渗透系数基本不变。对于非稳定变形，颗粒团（聚）体间重新生成“条带”状孔隙，渗透系数反弹增大。该研究成果可为流固耦合分析中渗透系数的确定提供理论依据。

砂岩受到高温影响后，宏细观性质会发生不同程度的变化，400 ℃ 和1 000 ℃ 是砂岩宏细观性质变化的两个重要节点。为进一步研究升温速率对高温砂岩物理性质、力学特性以及内部微观破裂的影响，以350 ℃和 950 ℃为目标温度，开展了不同升温速率下砂岩的单轴压缩试验，单轴压缩全程采用声发射监测，采用扫描电镜分析了砂岩破坏后的细观形貌。研究结果表明：经不同升温速率处理后，350 ℃砂岩的质量、体积、密度以及纵波波速变化较小；950 ℃砂岩的质量、密度及纵波波速显著降低，体积明显增大，且升温速率越高变化率越小。350 ℃砂岩的应力−应变曲线、应力−体积应变曲线、抗压强度及弹性模量受升温速率的影响较小；950 ℃ 下砂岩试样应力−应变曲线随升温速率的增加向上偏移，抗压强度和弹性模量则先增加后达到恒定。350 ℃ 及 950 ℃砂岩的声发射振铃累计数都有随升温速率增大而减小的趋势，950 ℃ 下升温速率为1 ℃/min 时砂岩的振铃累计数最大，突发性声发射活动区域最多。对不同升温速率砂岩进行扫描电镜分析发现，升温速率对350 ℃ 砂岩的微观形貌影响较小，950 ℃ 砂岩试件则会随着升温速率的下降出现微裂纹和微孔隙数量增多、体积扩大的情况。

岩石三轴强度包络线通常为非线性，通过线性回归法（LRM）获得的抗剪强度参数对围压的变化具有很高的敏感性。为使 LRM能够考虑围压对抗剪强度参数估计的影响，定义了岩石三轴强度围压效应系数，构造一个指数函数表示该系数与围压的关系，并将其引入到 LRM 的修正中，提出考虑围压效应的线性回归分析法（CCPE-LRM）。同时，提出了CCPE-LRM 的合理性检验方法，并定义了距离系数作为指标来表征抗剪强度参数估计值与参数实际值的差异。通过以公开发表文献中多种类型岩石的三轴强度数据为例进行验证分析，结果显示各种岩石的距离系数均较小，CCPE-LRM 获得的抗剪强度包络线紧贴摩尔圆，且近似相切状态。表明通过CCPE-LRM 获得的抗剪强度包络线在一定程度上可以代替理想抗剪强度包络线，该方法估计的抗剪强度参数与理论抗剪强度参数的一致性较高，证明CCPE-LRM 具有较好的适用性。

棚洞等刚性防护结构中具有良好缓冲性能的垫层材料可以有效削减冲击荷载，提高结构抗冲击能力。为研究重复冲击下再生混凝土垫层宏观缓冲性能演化规律，开展了落锤冲击试验和离散元数值模拟研究。试验结果表明，相比于石英砂垫层，首次冲击下采用再生混凝土骨料垫层防护的棚洞顶板中心处传递力减小了83%，且应力分布更加均匀。随着落锤重复冲击次数增加，再生混凝土骨料和石英砂垫层缓冲性能均存在劣化现象，第6次冲击时棚洞顶板中心传递力峰值分别为首次冲击的11.2倍和1.4倍。此外，数值模拟研究结果发现，当长条状颗粒数量占比由0%增大到100%时，颗粒旋转角度和平移距离分别减小了40%和20%，落锤冲击力峰值增加了37%。颗粒形状越不规则，颗粒间咬合程度越大，限制了颗粒旋转和平移，增大了滚石冲击力和挡墙传递荷载。该研究成果可为再生混凝土骨料作为新型绿色垫层提供理论依据和工程指导。

海洋能源开发不断向深海迈进，大型深水海上平台对锚泊结构的承载力提出了更高的要求。针对传统鱼雷锚承载比不足的问题，提出了一种新型深水注浆锚。锚体带动连接其尾部的注浆管贯入海床，然后在外部注浆设备驱动下向海床注浆，浆液挤压土体并在锚体周围形成注浆固结块，从而大幅增加锚体抗拔力。采用自主设计的海洋锚注浆和拉拔试验系统进行了模型试验，研究了不同注浆量对浆液扩散特性以及锚体承载特性的影响规律。结果表明：浆液可以较好地包裹锚体，浆块呈倒锥形与锚体紧密黏结为一体，共同承受上部荷载，浆块承载比可达 24.6，远高于传统鱼雷锚的 2.4～4.1，从而验证了深水注浆锚的可行性。随注浆量的增加，浆块端部截面积和整体高度增加，注浆锚整体承载力也不断增大。

土−水特征曲线在研究非饱和土的水力与力学特性中发挥着重要的作用。生物炭具有多孔结构、高比表面积和强吸附的特性。将生物炭改性土应用于垃圾填埋场上覆盖层，因受自然环境因素的影响会使其水力特性发生改变。为了研究全吸力范围内生物炭掺量对生物炭−黏土混合土保水特性的影响，利用蒸汽平衡法（吸力范围 3～368 MPa）、滤纸法（吸力范围 0 ～40 MPa）和压力板法（吸力范围 0～1.5 MPa）控制土样的吸力，测定吸力平衡后土样的含水率和饱和度，得到全吸力范围内生物炭−黏土混合土的土−水特征曲线。试验结果表明：（1）3种吸力测试方法很好地表达了生物炭−黏土混合土全吸力范围内的土−水特征曲线。（2）生物炭能够影响黏土的保水性，但在一定的吸力范围内，生物炭−黏土混合土的保水性还与孔隙结构和孔隙中水的形态相关。（3）通过压力板法测得，试样的进气值随着生物炭掺量的增加而减小。当吸力值小于进气值时，曲线出现水平段，土样始终处于饱和状态，生物炭掺量越大，试样的保水性越好。（4）由生物炭−黏土混合土微观孔隙结构以及生物炭在黏土中的分布形态来解释生物炭改性黏土的保水能力随生物炭掺量的变化关系。

在渗流作用下土中细粒随着水流在粗粒骨架间运移的现象被称为内侵蚀。以往对于内侵蚀的研究主要集中在恒定的水力梯度下，缺乏对于波动水力条件下内侵蚀发展特征和细观机制的认识。设计了透明土渗流试验装置，开展了恒定和波动水力条件下不同内部稳定性土体的渗流试验，以探究波动水力条件下内侵蚀的发展特征。宏观试验现象表明，对于内部稳定性较差型土，当水力梯度达到临界水力梯度后，波动水力条件下的渗透系数比恒定水力条件下呈现出更快的增长趋势，说明水力波动会加剧细颗粒的迁移和流失。为了进一步揭示波动水力条件加剧土体内侵蚀发展过程的细观机制，利用一种基于平面激光扫描的二维截面图像开发的土体细观组构三维重建方法，建立了粗粒骨架、内部孔隙通道和细颗粒的三维可视化数字模型。通过观察细粒在孔隙通道中的分布情况，发现恒定水力条件下，细粒在运移至狭窄的孔隙喉道处时会发生堵塞堆积现象；而波动水力条件下的水流扰动会破坏这种暂稳态的堵塞堆积结构，继而重启细颗粒的迁移过程。

黄土边坡中竖直裂隙的发育往往会对边坡稳定产生影响。相对于平面应变机制，建立三维破坏机制下边坡稳定性分析方法更能接近实际边坡失稳情况。基于塑性极限分析上限法，考虑预先存在竖直裂隙的三维黄土边坡不同破坏机制（坡面破坏、坡脚破坏和坡底破坏），建立能量平衡方程及其无量纲临界高度值γH/c表达式，采用随机搜索法得到了临界高度的上限解。分析了约束宽度、边坡坡度、内摩擦角以及裂隙深度对三维竖直裂隙黄土边坡临界高度值的影响。结果表明：对于坡脚破坏机制，临界高度值随着裂隙深度的增加而减小，减小至临界裂隙深度 (δ /H) _min 后，裂隙深度的增加不再影响临界高度值；临界裂隙深度随着坡度 β  的增大而增大，随着内摩擦角 φ  的增大而减小。当约束宽度 B/H<0.8时，大多数破坏机制为坡面破坏。当约束宽度 B/H=0.8、内摩擦角 φ =10° 及约束宽度 B/H=0.6、内摩擦角 φ =15° 时，边坡的破坏从坡面破坏机制逐渐过渡到坡脚破坏机制。存在竖直裂隙的黄土边坡比完整边坡具有更小的临界高度，约束宽度及内摩擦角会对三维黄土边坡破坏机制产生影响。

温度的变化会导致土体的物理力学性质改变，且在一些实际工程中，饱和黏土会处于非等温分布状态。为此，针对非等温分布条件下饱和黏土的一维固结问题，考虑了更具普遍性的半透水边界，通过某些假定推导了单级线性加荷形式下饱和黏土一维固结控制方程，并利用分离变量法求解得到了控制方程的解析解。通过将所提解析解分别与已有解析解和有限差分解展开对比分析，验证了所提解答的正确性。基于所提解析解，利用某一算例分析了温度梯度、半透水边界参数及加荷时间对固结性状的影响。结果表明：温度梯度 M 越大，土体的渗透性越大，土体的固结速率越快；半透水边界参数 R1和 R2越大，相同时间内土体的超孔隙水压力越小，土体的平均固结度越大；土体的平均固结度随加荷时间 tc 的增大而减小，这主要是由于加荷阶段所施加的外荷载小于最终荷载，但加荷时间tc的延长可一定程度减小土体中产生的最大超孔隙水压力。

常规三轴压缩试验中具有较强结构性的黏性土在围压较低时其应力−应变关系会呈现应变软化现象，一般还伴有塑性变形，通常土体内部结构损伤是应变软化产生的主要原因。考虑到采用经典塑性理论描述材料的应变软化不仅会违背 Drucker 的稳定性假设，而且也不能描述卸载塑性。因此，基于修正剑桥模型及 Li 和 Meissner 提出的塑性硬化准则，建立了一个描述饱和黏性土不排水应变软化的弹塑性双面模型。该模型以应力−应变曲线的峰值点分界，将应变硬化和应变软化分别作为独立的加载事件进行分析，同时引入新的结构性参数表征剪切过程中土体结构损伤导致的塑性刚度衰退。对不同固结状态饱和结构性黏土的三轴固结不排水压缩试验结果的模拟表明，所建模型能够较好地描述饱和黏性土的不排水应变软化特性。

在能源岩土工程发展需求背景下，土体工程特性的温度效应研究受到了国内外学者的广泛关注。针对现有关于循环变温热固结理论研究的不足，提出考虑循环变温条件下土体先期固结压力随温度循环次数变化的关系式，计算多次循环变温后软土的固结压缩量。在此基础上，推导考虑循环变温的太沙基固结修正公式，分析整个循环变温过程中孔压和固结沉降的变化规律。最后通过与单元体测试结果进行对比，验证该计算模型的可靠性。研究结果表明，在循环变温和瞬时荷载的耦合作用下，土体的沉降变形和孔压随温度循环次数的变化而变化，且影响程度与循环次数相关；新型热固结模型为涉及循环变温的岩土工程问题分析提供了一定的理论支撑。