岩质高边坡稳定性分析是岩土工程中经典的研究课题，其假定不同，分析方法就不同。从现有稳定性分析方法中存在的问题入手，指出其存在的优劣性，在刚体极限平衡和有限元等传统边坡稳定性分析方法的基础上，提出了用有限单元法直接求解边坡稳定性。该方法的思路是用潜在滑动面上在开挖、支护完成后尚未发生滑动的所有滑面单元模拟边坡在各工况下的真实应力场和变形场，根据这个真实应力场直接计算边坡安全系数。结合现场即时监测手段和计算机仿真技术可逐步模拟边坡在分步开挖期间的稳定性。算例分析表明，该方法完善了刚体极限平衡法的不足，克服了强度折减法强度参数c、?值折减比例的难题，弥补了超载法力学方向不明确的缺陷。计算所得的边坡安全系数和稳定状态与工程实际相当接近，表明提出的方法是可靠有效的新方法，为复杂岩质边坡稳定性分析提供了新的思路和途径。

针对比利时HADES地下实验室PRACLAY现场加热试验，应用温度-渗流-应力耦合弹塑性模型，模拟现场加热过程中泥岩核废料处置库的水力学响应特征。采用单因素分析法，就泥岩热、水、力学参数对核废料处置库围岩孔压、温度、有效应力的影响进行了三维有限元分析。并基于参数敏感性分析结果，就温度、渗流、应力三场两两耦合作用对处置库围岩水力学响应的影响程度进行了系统分析。研究结果表明：泥岩热、水、力学参数中，渗透系数、弹性模量以及导热系数对加温所导致的超孔压的值影响较大；凝聚力、内摩擦角以及热膨胀系数对孔压的影响较小，但会显著影响围岩的有效应力；导热系数对围岩温度场的分布有决定性影响，温度传递的差异会显著影响围岩的孔压和有效应力；不同的热、水、力学参数对孔压、温度以及有效应力的影响机制是不同的，温度、渗流、应力三场两两耦合作用对围岩水力学响应的影响程度也存在显著的差异性。温度场对应力场、温度场对渗流场的耦合效应十分显著，加热后，围岩超孔压的产生以及热膨胀导致的有效应力变化会显著影响处置库的稳定。该研究结果在一定程度上可以为核废料处置库泥岩的热、水、力学参数的确定及耦合机制分析提供科学依据。

为了客观重复制作与原岩结构面的岩性、表面形态和力学性质完全一致的试样，集成三维扫描和三维雕刻技术，提出了原岩结构面试样重构新方法。采用该方法制作了具有相同表面形态的大理岩结构面试样，开展了不同法向力下的直剪试验，并采用声发射技术对结构面剪切破坏过程进行了监测。结果表明：重构结构面试样与原始结构面试样具有高度相似性；结构面剪应力-剪切位移曲线可分为两类：剪胀滑移型和剪断跌落型；剪切破坏后结构面可明显地划分为两个区域：剪胀脱开区和剪断磨损区；提出了结构面剪切强度新公式，理论公式与试验结果更接近；结构面剪切过程中剪应力、能量率、计数、撞击率随时间的变化规律具有较好的一致性，声发射定位位置与剪切过程中结构面损伤破坏区域基本吻合。

比利时放射性核废料地质备选场址Boom clay是一种典型的横观各向同性材料，在Drucker-Prager帽盖模型的基础上，构建了适用于Boom clay的横观各向同性特点的热-水-力耦合弹塑性本构模型，该模型可反映温度对其强度、弹性模量、渗透性等的影响，并在ABAQUS中进行了二次开发。为验证所建立模型的合理性，结合比利时HADES地下实验室ATLAS III现场加热试验结果，应用所提出的模型对加热过程中围岩的温度和孔隙水压力的变化规律进行了数值仿真分析。结果表明：所建立的模型能够正确地描述现场加热试验过程中围岩温度场和孔压场所呈现出的各向异性特征，主要表现为热源水平面内测点的孔压在加热功率升高时先略为下降后才升高，在加热功率下降时先略为升高后才下降，而竖直面内测点的孔压在加热功率升高时立即升高，在加热功率下降时立即下降。研究成果表明，考虑各向异性的THM耦合分析能更好地反映加热过程中泥岩温度场和孔压场的实际分布情况，研究结果对类似工程或现场试验的设计、安全运行提供重要的决策依据。

为研究吸力及压实度对非饱和压实填土压缩变形特性的影响，并建立脱湿（吸力增长）状态下的填方土体工后沉降变形修正计算应用模型，开展了控制吸力和压实度的一维非饱和土侧限压缩试验。结果表明：压缩曲线平缓段随压实度和吸力的提高而增长，表明土样的结构屈服强度同步得到提高；在比容v变化差值序列上，饱和土压缩土样最大、常规压缩土样次之、脱湿土样最小，并且吸力越大，比容v变化差值越小。经历脱湿（吸力增长）后的压实土压缩性降低，定义和建立了吸力压缩系数 及其经验模型，用以表征和度量吸力和压实度对压缩特性的影响规律及程度，发现 随竖向应力增加呈现指数型衰减。对模型参数与压实土相关参量间的关联性进行了分析探讨，同一压实度土样不同吸力所对应的参数 均值随压实水平提高而线性减小，参数λ总体上随吸力的增加而增大，但随压实度的提高吸力对于土体抗压性增强的贡献水平降低且参数 试验数值点靠拢线性趋近线，土样的压缩性随压实度和吸力大小变化而动态调整。基于分层总和法的基本原理初步构造建立了高填方非饱和填土压缩变形修正计算模型（MS），其应用途径应建立在进一步对非饱和压实填土受荷状态下脱湿土-水特征关系的研究基础之上。

为揭示深埋碳质千枚岩的遇水软化损伤特性和微层理结构的各向异性，开展了不同层理角度和不同含水状态下的千枚岩力学试验。并基于能量损伤演化机制分析了层理角度和含水状态对千枚岩储能机制和释能机制的影响。研究结果表明：深埋碳质千枚岩的各向异性特性和遇水软化效应是千枚岩内部的沉积层理、层状薄片矿物和高黏土矿物成分导致的固有特性；其破坏机制及其强度各向异性主要受控于层理面间的弱胶结作用，这种弱面效应会随着围压的增大而降低；千枚岩遇水后其力学特性和脆性被弱化，宏观破坏角增大，表现出张拉破坏减弱，剪切破坏增强；基于弹性能和耗散能比例曲线的S型变化规律与能量演化规律，将岩石的渐进破裂演化过程划分为4个阶段，并采用能量损伤演化理论探讨了I型和II型两种典型岩石破坏方式的能量机制及其损伤演化机制；层理方向和含水状态对千枚岩的储能能力具有较大影响，其中层理方向对千枚岩的释能机制和损伤破裂演化机制影响不大，而千枚岩的储能与释能机制则对含水状态的变化更为敏感。

为了进一步理解岩爆的机制并解释岩爆发生的过程，通过对抛掷型岩爆能量来源、去处及大小的分析，提出了抛掷型岩爆的震源两体模型。抛掷型岩爆本质上是释能体-岩爆体相互作用下岩石裂纹动态扩展和能量快速释放的过程。从岩爆体张拉断裂机制出发，采用脆性内聚裂纹阵列扩展动力学模型分析了裂纹间距与应变率对岩爆体能量释放速率的影响。据此得出释能体对岩爆体的作用可归结为应变率效应，释能体能量释放速率越大，对岩爆体的动态加载速率就越大，使得岩爆体能量释放速率加快。岩爆体破裂能量释放速率越大，作用于释能体的动态卸载速率就越大，又使得释能体的能量释放速率增大，两体相互作用形成的动态加、卸载效应构成了岩爆能量释放的正反馈机制。基于上述理论推导得出岩爆震源体的几何尺度是岩爆体的几倍到几十倍，理论结果与现场实际监测分析基本吻合。

珊瑚礁灰岩层中桩基承载特性的研究是岩土工程的热点问题。基于某跨海大桥工程开展的珊瑚礁灰岩层中桩端后压浆桩现场静载荷试验，对比分析了压浆前后的实测结果，研究了桩端后压浆的预压作用对桩阻力的影响，并在桩端后压浆试验结果分析的基础上，进一步研究了后压浆桩增强效应作用机制。结果表明，桩端后压浆技术可应用于珊瑚礁灰岩地层中，能有效地提高桩基承载力和减小沉降量；钻孔取芯试验明确了压力浆液在桩端以下一定范围内的分布情况，并证实了压力浆液对珊瑚礁灰岩的孔隙有充填作用；桩端压浆后，在桩端未能消散的压力产生负摩阻力，使负摩阻力在竖向荷载作用下提前发挥，并增大了压浆桩的侧阻力。此外，桩端后压浆的预压作用提高了桩端阻力，并减小了桩端位移，改善了桩基的承载性状。

滑面的摩擦特性是高速岩质滑坡动力学研究的关键问题之一。以典型脆性灰岩为研究对象，借助滚动磨损试验设备，开展不同压力和转速条件下干燥和湿润灰岩的滚动摩擦试验，研究岩石动态摩擦特性（摩擦系数和磨损率）。试验结果表明，干燥和湿润岩石的滚动摩擦系数分别与摩擦速率呈负相关和正相关，与法向压力关系不明显。对于干燥灰岩，滚动试件（滑体）和静止试件（滑床）的磨损率与摩擦速率分别呈正相关和负相关。对于湿润灰岩，磨损率与摩擦速率呈正相关。基于Hertz弹性接触理论，推导出摩擦斑处的最大压应力和拉应力，进而提出相应条件下脆性灰岩的摩擦碎裂机制。结论可为灰岩高速滑坡防灾减灾提供设计参数。

颗粒形态是影响砂土力学特性的重要因素，特别是影响砂土在低应力状态下的抗剪强度、剪胀效应和临界状态行为，以及高应力状态下的颗粒破碎行为。因此，准确地重构砂粒的三维形态，并量化计算其形态表征参数是研究砂粒形态效应的前提工作。借助于高精度的CT扫描技术和图像处理技术，获得近海石英砂和风化花岗岩残积砂这两类砂土颗粒的三维形态信息。采用球谐函数序列实现两种砂颗粒三维形态的准确重构，并通过球谐函数分析计算砂土颗粒的体积来验证该方法的有效性。基于球谐重构的三维砂粒表面，提出了实用性的方法来计算砂粒的表面积、表面曲率和三维尺寸等，进而计算砂粒的三维球度、圆度和伸长率等形态表征参数。结果表明，当球谐函数阶达到15时，其重构的砂粒基本形状和表面纹理均与真实砂粒非常接近；近海石英砂在水流搬运和磨蚀的作用下颗粒形态较为规则和圆滑，球度和圆度较大，而风化花岗岩残积砂则在物理风化和剥蚀作用下颗粒形态较为复杂和粗糙，球度和圆度较小；而这两种地质作用对砂土颗粒的伸长率则没有明显的影响。

钙质砂是富存于热带海洋环境（包括中国南海海域）中的一种特殊岩土介质，具有不同于陆源砂的水理性质。采用传统常水头渗流试验，首先探究了不同不均匀系数和曲率系数条件下级配对钙质砂渗透性的影响。针对钙质砂的颗粒形状特性，采用扫描电镜（SEM）与图像处理技术，引入球度 和圆度X的比值 从三维空间角度上对颗粒形状进行了定量描述。并在粒径区间、相对密实度相同的条件下，通过与福建标准砂、玻璃珠进行对比试验来考察钙质砂渗透特性，从而进一步探究颗粒形状对钙质砂渗透性的影响。试验结果表明，钙质砂渗透性随着颗粒粒径不均匀系数和曲率系数的增大而增大，符合级配对一般砂土渗透性的影响规律；钙质砂颗粒形状具有较强的结构性和不均匀性，同等密实度下钙质砂的渗透性小于陆源的石英砂。本研究获得的钙质砂渗透规律对今后南海岛礁填筑和海上平台基础工程设计具有一定的指导意义。

锚固浆液的类型及其与杆体、土体的兼容协调性问题是土遗址锚固领域的研究焦点。由糯米浆为胶凝主材和由黏土、粉煤灰为填充料组成的改性糯米灰浆与木锚杆、玻璃纤维增强塑料锚杆、钢筋锚杆组成了3种锚固系统。通过对3种锚固系统的原位锚固、拉拔试验和界面应力-应变监测，分析比较了3种锚固系统的破坏模式、极限荷载、荷载-位移特征，剖析了杆体-浆体界面应力分布传递规律和界面测点应变对荷载时步的响应程度。在初步分析其锚固机制的基础上，最终得出了此类基于改性糯米灰浆的3种锚固系统的锚固性能的优劣完全取决于杆体-浆体的受力机制、变形和强度特征所决定的力学兼容性的结论。该研究成果为以糯米浆为主材的改性浆液在土遗址锚固领域中的应用提供了依据和参考。

粗粒料的力学特性不仅取决于应力状态，还与粗粒料本身的松密程度密切相关，即与粗粒料的材料状态相关。通过不同初始干密度的粗粒料常规大型三轴各向等压固结排水剪切试验研究了干密度对粗粒料力学特性的影响。试验结果表明：粗粒料应力-应变曲线的形态取决于密度和围压的共同作用，破坏状态之前密度和围压共同决定了粗粒料的强度和变形，而当应变足够大处于渐进状态或临界状态时，粗粒料的应力和体积应变受初始干密度的影响逐渐减小直至消失。对于软化型曲线，相变状态的应力小于渐进状态的应力，而渐进状态的应力小于破坏状态的应力。无论围压大小如何，随着初始干密度的增大，应力-应变曲线的硬化性逐渐减弱，而软化性则逐渐增强。

为观测和分析安装中拖曳锚在土中的运动特性，自制了由玻璃水槽、拖曳系统和拉力传感器组成的板锚拖曳试验装置。将硅凝胶与水混合制备成了透明土，依据有效重度强度比确定模型锚材料，按照1:50缩尺比例，在自制板锚拖曳试验装置上开展固定锚腚拖曳锚拖曳试验，研究了偏心系数和锚腚角对固定锚腚拖曳锚运动特性的影响规律。结果表明：该试验装置可以直接观测拖曳板锚嵌入土中运动的全过程；锚眼法向偏心距和切向偏心距不是独立参数，偏心系数能综合描述锚眼偏心对板锚在土中运动的影响；锚腚角是影响板锚在土中运动的重要参数，随着锚腚角增加，板锚最终嵌土深度先增加后减小；随着锚眼偏心系数增加，板锚最终嵌土深度先增加后减小；固定锚腚拖曳锚在土中的运动轨迹近似服从负指数函数曲线。

扩散是低渗透系数情况下溶质迁移的主要方式，目前尚没有测定高塌落度防渗墙填料扩散系数的试验方法。根据沥出试验原理，将处于流动状态的填料用半透膜包裹进行透析试验，通过有限圆柱中溶质迁移解析解对溶质迁移过程进行拟合求其扩散系数，探讨了试验和数据处理的相关问题，研究了膨润土含量和盐浴浓度对扩散系数的影响规律。结果表明，随拟合天数增加，填料中NaCl有效扩散系数减小，而误差增大；高塌落度填料中溶质迁移明显比压实土样中快，按拟合误差小于0.5%作为试验数据取舍依据，这一误差要求实际上是对溶质迁移以扩散为主导这一条件的量化，保证了采用纯扩散解析解进行数据拟合的合理性；扩散系数随填料中膨润土含量增加而减小，随盐浴浓度增加而减小，但在试验范围内数值变化不大，透析试验是快速测定防渗墙填料中溶质有效扩散系数的可行方法。

分别开展砂土和粉质黏土两种典型土质条件下基坑悬臂式开挖离心模型试验，详细叙述试验过程中所要解决的关键问题，并提出合理的解决方案。通过对比分析两组试验结果，得到以下结论：非饱和土地基制备中参数控制困难，分层夯实法有待进一步改进，而砂雨法制备的砂土地基参数可控性更好；两组试验的结果有差异，砂土地基试验所呈现的土压力、地基变形、支护弯矩的变化规律更好，因此，岩土离心试验可适当考虑以砂土代替非饱和土；对于采用悬臂式支护结构的基坑，开挖引起的地表沉降曲线在砂土中呈指数型，而在粉质黏土中呈直线型；开挖引起的粉质黏土地基土体位移范围较砂土地基更大；开挖引起的砂土中挡墙弯矩较粉质黏土更大，砂土和粉质黏土中最大弯矩位置都随开挖逐渐下移；在砂土试验中开挖引起主动区土压力各处均减小，而在粉质黏土试验中开挖引起土压力在挡墙底有增大趋势。该基坑工程离心模型试验过程及数据处理方法可为进一步试验提供参考。

为研究列车振动荷载作用下盾构隧道结构及周围土体的动力响应特性，采用模型试验方法，通过布置在盾构隧道底部的激振器施加扫频激振荷载和列车振动荷载，采用频率响应函数FRF与最大加速度分析了盾构隧道衬砌结构与周围土体不同位置处的动力响应及其衰减规律。研究结果表明：FRF是隧道衬砌结构和周围土体自身的动力响应特性的体现，与激振力的大小、扫频方向及扫频时间无关；在隧道管片衬砌结构的底部和顶部均体现出高频响应大于低频响应的特性，隧道顶部加速度响应沿隧道纵向衰减幅度明显小于隧道底部；隧道周围土体的动力响应沿深度有明显变化，但均表现出沿隧道轴向衰减的规律。隧道结构上部第1层测点土体的动力响应在全频域内随频率的增加逐渐增大。但在第2层和地表的第3层测点，土体的动力响应在30～90 Hz区段线性增长，在90～300 Hz区段出现波动变化，并无明显增大趋势；与扫频激振荷载引起的动力响应规律一致，由列车运行振动荷载引起的隧道管片衬砌结构和周围土体的振动沿隧道轴向逐渐衰减，隧道底部的加速度响应大于顶部，随着列车车速的增大，在隧道内部引起的加速度响应将显著增大。同时，在列车振动荷载作用下发现地表存在加速度放大效应，地表第3层测点的加速度响应均大于隧道结构上部第1层测点。

由于黏性土表面带有丰富的负电荷，孔隙水溶液化学状态的变化对黏性土的物理力学特性存在明显影响。随着化学-力学耦合的相关岩土工程问题日益突出，进行有效的化学-力学耦合行为的数值分析评价显得尤为重要。因此，建立一个简单有效的考虑化学-力学耦合的本构模型是非常关键的。基于传统的修正剑桥模型，提出了一个简单的化学-力学耦合模型。该模型采用渗透吸力π描述孔隙水的化学状态，建立了前期屈服应力，临界状态线斜率M和弹性刚度与渗透吸力π之间的关系式，从而实现了模型对盐溶液饱和黏性土的变形和强度特性的有效模拟。通过与试验数据的对比和分析，说明该模型能有效地模拟孔隙盐溶液饱和黏性土的等向压缩行为、 状态下压缩行为以及 状态下化学-力学循环加载行为。此外，通过对黏性土三轴压缩试验的模拟，说明该模型能反映黏性土三轴应力状态下的基本力学特征。

多年冻土地区的排水沟渠经常遭受水平冻胀力的破坏，严重影响构筑物的服役性能。根据青藏高原多年冻土区路基坡脚U型槽排水沟渠的现场试验，并通过数值计算与理论分析，分析了结构埋深、粗颗粒换填范围以及结构形式对其温度场和水平冻胀力的影响。结果表明：改变工况，对土体温度场分布的形状影响较小，仅对其大小有一定影响。不同工况下水平冻胀力均沿深度非均匀分布，最大水平冻胀力主要出现在结构中部，而其上部和下部较小，结构侧壁在1/2至1/3处易发生冻胀破坏。梯形结构所受的水平冻胀力较U型结构增大13%～15%左右，但其分布形式基本相同。因此，U型结构在降低水平冻胀力方面优于梯形结构。随着换填范围增大，排水沟渠的水平冻胀力最大值逐渐减小；对于埋深1.7 m的排水沟渠来说，其侧边的0～2.8 m是影响水平冻胀力的主要换填宽度范围，而当换填宽度超过2.8 m后，水平冻胀力几乎不再降低。

活动层水热状态直接影响多年冻土和寒区工程的稳定性。已有研究大多基于附面层理论研究冻土温度场变化，较少研究液态水和水汽运移过程及其对冻土温度场的影响。结合多年冻土活动层包气带水-热耦合迁移的物理过程和内在机制，以温度和含水率为基本变量，建立了考虑液态水和水汽相变、水分对流传热和水汽运移的土壤-地表-大气能量平衡及土壤内部水热变化的耦合模型，分析了真实野外气象条件下活动层液态水和水汽运移规律。计算结果表明：白天温度梯度水分向土壤内部运移，夜间温度梯度水分向土壤表层运移；暖季温度梯度水分以向土壤内部运移为主，冷季温度梯度水分以向地表运移为主；就全年而言，活动层各个深度处水汽运移作用大于15%，水势梯度水汽运移极小可以忽略不计，特别是浅层土壤在无降雨状态下，水分运移以温度梯度水汽运移为主；水势梯度液态水通量受降雨影响明显，在降雨事件期间和之后，液态水和水汽下渗占主导地位，降雨降低表层土壤温度、减小土壤热传导通量，有利于土壤热稳定性。

为研究不同初始破损程度、不同干湿循环次数及两者耦合工况对膨胀土孔洞及裂隙发育、演化规律的影响，以及相应的土体变形、力学行为，以合肥膨胀土为研究对象，对制备的不同初始破损程度（圆柱孔直径分别为0、2.5、5 mm）的3组20个试样进行干湿循环试验及三轴剪切试验（控制吸力为50 kPa，有效围压分别为100、200、400 kPa）。并解释了孔洞-裂隙的产生及发育新机制。结果表明：干湿循环次数为0时，一定程度的初始破损会提高土体的强度；与常规认识不同，干湿循环次数为1～2次时，初始破损的存在会诱导土体裂隙开裂的方向和贯通模式，形成的结构块较无破损土体形成的更完整，因而具有相对于无初始破损土体更高的强度；干湿循环次数为3次时，裂隙完全扩展，土体结构性被严重破坏，有无初始破损的土体强度趋于一致；孔洞-裂隙的演化机制是：脱湿-增湿过程中产生的水力梯度使土体具有不同的软硬性质和胀缩变形，产生的拉压应力在微裂纹、尖端形成应力集中而开裂，更易与具有的初始孔洞破损联合、交汇、贯通，萌生的裂隙增大了水汽交换界面面积，进一步平衡拉压应力而最终完全碎裂。研究可为普遍存在初始破损的膨胀土工程研究提供新的思路。

对哈尔滨原状粉质黏土进行共振柱试验，通过扭转及弯曲试验得出相应的弹性模量E、剪切模量G、泊松比ν，旨在分析低幅应变下粉质黏土泊松比随应变发展的非线性特征。试验结果表明：随着应变的增大，粉质黏土剪切模量与弹性模量均出现衰减，但二者衰减并不同步，剪切模量先于弹性模量发生衰减；泊松比随着应变的增大而增大，并随着围压的升高而降低；在半对数坐标系下，泊松比与归一化剪切模量 近似呈线性关系，基于此提出低幅应变下泊松比预测模型。理论推导证明：低幅应变下ν- 曲线呈线性具有其必然性；随着应变的发展和 值的降低，ν- 曲线由斜率一定的线性段过渡到斜率逐渐降低的非线性段，并最终趋近于一定值。基于理论推导建立了泊松比预测模型参数与Hardin-Drnevich模型参数的关系，由此可依据低幅应变下ν- 曲线形态对较高应变下泊松比进行预测，实例验证表明该方法预测结果较好，具有一定应用价值。

依据室内试验所得的石膏岩的基本力学性质随着大气中相对湿度变化的规律，并结合石膏矿采空区的破坏特征，将护顶层简化为铅直方向受均布荷载q作用的三铰拱分析模型，建立了相对湿度作用下的石膏矿采空区护顶层的尖点突变模型，研究了相对湿度作用下护顶层突发破坏的失稳机制。试验结果表明：石膏岩的力学性质在采空区中环境湿度的作用下，其强度和变形参数会大大降低。分析和工程计算结果表明：随着采空区中环境湿度的增加，支撑体系越不稳定，越容易跨越分叉集而发生突发失稳破坏；相空间的控制参数a随着采空区中环境湿度的增加，其代数值是逐渐增大的（绝对值减小，由负值向0变化），而控制参数b逐渐增大，系统突发失稳控制参数（a，b）的变化路径为向东北方向变化，该路径是系统较为容易发生突发失稳的路径。

基于Young-Laplace方程利用土-水特征曲线估算孔隙分布，借鉴非饱和渗流统计模型，建立了土-水特征曲线的滞后模型。在某一基质吸力下对应的毛细半径上限的毛细水对更大孔隙半径的毛细水具有阻塞作用，阻塞概率与孔隙分布函数直接相关。孔隙分布函数自身属性体现了土体孔隙空间分布的非均匀性。该模型显示：在高基质吸力与低基质吸力阶段，脱湿曲线与吸湿曲线趋近相等；中等基质吸力阶段，脱湿曲线饱和度高于吸湿曲线饱和度，两者的差值存在一个明显的峰值。并采用实例验证了该模型的可靠性，发现该模型对中、细粒土预测效果较好，而砂类土因孔隙分布不符合假设，导致存在较大误差。对于砂类土等特征尺寸较大的土质，引入阻塞概率修正系数，发现最佳修正系数与土-水特征曲线半对数坐标下的土-水特征曲线最大斜率呈反比例关系。

自三峡水库蓄水后，库区地质灾害频发，研究库水位升降条件下的库岸边坡稳定性与破坏模式日益重要。以三峡库区云阳县凉水井滑坡为原型，基于长江科学院CKY?200土工离心机，结合滑坡区工程地质概况，模拟制作了这一典型库岸滑坡的简化模型。采用高速相机、激光位移传感器、孔隙水压力传感器、土压力传感器等监测手段，在模型试验中实现了对边坡前缘水位高程及升降的控制，获取了试验过程中坡体水平位移和沉降、孔隙水压力、土压力值的变化规律，得到以下结论：试验从开始至加水阶段结束，模型位移以竖直方向为主，而在排水阶段则水平方向位移远大于竖直方向；试验过程中，在库水位上升与骤降作用下，首先在前缘产生裂隙，并逐渐向中部、后缘扩展，随后在中部产生较大变形，出现局部滑移，模型总体破坏形式呈牵引式破坏模式；试验结束时，模型自后缘起体积含水率逐渐增大；在水位上升期间，模型变形相对较小，主要滑移出现在水位下降期间，因此，该滑坡为动水压力型，产生变形破坏主要是由于动水压力效应所致。试验揭示了库岸牵引式滑坡在库水位升降条件下的破坏模式和变形失稳机制，为库区滑坡的防治提供了重要依据。

为构建能够反映砂土高低压下不同力学特性的统一模型，分析了砂土在较大的压力范围内的力学试验结果并获取其强度、等向压缩以及临界状态特性参数。通过引入应力路径相关因子来修正塑性应变增量中与应力路径相关的部分，从而使得模型硬化参量能够反映密实砂土在常压下的剪胀特性。同时，基于砂土的临界状态特性提出能够与砂土内部状态相对应的潜在状态面概念，由屈服面与潜在状态面间的动态关系确定加载过程中的动态密实参数及潜在强度，进而使得硬化参量也能够反映密实砂土在常压下的软化特性及高压下的剪缩、硬化特性。分析模型屈服面及潜在状态面间的演化规律并对不同压力等级下的砂土受荷力学行为进行模拟预测，证实了该模型能够反映密实砂土常压下剪胀软化及高压下剪缩硬化的特性。

为了模拟深部圆形洞室三维受力情况下岩爆的发生过程，以具有极强岩爆倾向性的花岗岩作为试验材料，利用自行研制的TRW?3000岩石真三轴电液伺服诱变试验机，对含直径50 mm贯穿圆形孔洞的100 mm×100 mm×100 mm立方体花岗岩试样进行了三向不等压加载试验。试验首先模拟1 000 m深度的初始应力环境，以垂直于洞室轴向的水平方向作为最小主应力方向，通过增加竖直方向应力模拟应力调整过程诱发岩爆情况，并借助三维加载岩样内部结构破坏实时视频监控系统监控洞壁破坏情况。试验结果表明，岩爆过程可以划分为平静阶段、颗粒弹射阶段、岩片剥落和爆裂阶段，洞壁两侧中间部位应力集中系数最高，最先产生破坏，然后沿径向向深部发展，最终形成两个对称型的V型槽。圆形洞室洞壁两侧破坏的3个典型受力节点，符合深部工程中统计得到的无支护圆形巷道的远场应力状态与破坏模式之间的关系，达到了模拟的效果。利用现有的4种岩爆判据进行了分析，判别结果与试验过程中拍摄到的洞壁破坏情况基本一致。洞壁的破坏伴随着声发射计数增加和能量升高，岩爆越剧烈，声发射越活跃。

煤是自然界中一种常见的沉积岩，它具有显著的各向异性特征。然而目前针对煤体渗透特性研究，多数学者为了简化问题，多假设煤体为各向同性材料，提出了相应的各向同性渗透模型。这类模型并不能完全反映含瓦斯煤气-固耦合真实工程和室内试验的实际情况。假设煤体为横观各向同性，推导出以不同方向模量损失率为关键参数的煤体各向异性渗透模型，在此基础上，推导出含瓦斯煤的气-固耦合控制方程，并植入Comsol计算平台，系统研究煤体各向异性对气体扩散和渗透的影响。理论和数值研究结果表明：不同方向的模量损失率 反映出煤体结构各向异性变化程度，若 不同，其煤体各方向渗透特性也不相同；煤体渗透率的改变主要受解吸附效应和有效应力作用双重影响， 反映了这两种效应对于渗透率的影响程度；单轴应变或位移控制边界条件下，水平方向的模量损失率 对于垂直方向的煤体渗透率改变量 的影响程度大于对水平方向的煤体渗透率改变量 的影响程度，垂直方向模量损失率 对 的影响则弱于对于 的影响。

针对三轴压缩试验中含夹层盐岩特殊的破坏形式和机制开展研究。基于夹层单元的受力平衡和界面剪应力的传递机制，建立了夹层径向位移的微分方程，通过求解该微分方程得到了界面剪应力、夹层应力和应变的理论公式，并采用数值模拟验证了应力公式的准确性，进而结合多种破坏准则分析了含夹层盐岩的两种破坏形式及其机制。结果表明：界面剪应力和夹层应力分布不均匀，其最大值一般分别在边缘和中心位置出现；对于试样中夹层纵向劈裂，计算显示夹层存在横向的拉应力或拉应变；在一些试样界面边缘发现了环状白色裂纹，这是由于界面边缘剪应力最大，该处容易发生界面滑移破坏。夹层和界面的破坏会导致含夹层盐岩的密闭性降低，研究成果可在这一方面提供参考。

基于热力学原理及二维液桥计算模型，以3个不等径土颗粒为研究对象，推导出基质吸力与体积含水率之间的关系，得到微观模型下的土-水特征曲线。将3个不等粒径颗粒概化为限制粒径、中值粒径和有效粒径，研究了粒径大小、颗粒级配和土颗粒接触角的大小等因素对土-水特征曲线影响的规律。计算结果表明：体积含水率相同，颗粒粒径越大，相应的基质吸力越小；在同一基质吸力下，不均匀系数 越大即土颗粒越不均匀，其体积含水率就越低；接触角对土-水特征曲线的影响是明显的，即接触角越小土颗粒亲水性越好，同一基质吸力下的体积含水率也就越大，且曲线的进气值会随着接触角的减小而减小。

以石膏、浮石、硅藻土等材料为相似材料对多孔含水岩层进行相似材料模拟。试验中采用四因素（硅藻土水泥石膏比、标准砂浮石比、骨胶比、重晶石骨料比）、五水平的正交试验方法，分别对多孔含水岩层的相似材料试件进行密度、单轴抗压强度、孔隙率、弹性模量以及软化系数的测定，得出相应的物理力学性质指标，采用直观分析法分析各因素对物理力学性质的影响规律。分析结果表明：由于各配比度不同，相似材料强度分布范围较广，在与多孔岩层的岩石参数对比中发现，在一定的相似条件下，基本符合多孔岩层岩石的相似要求；硅藻土含量与密度、抗压强度、弹性模量、孔隙率呈负相关，与软化系数呈正相关；浮石含量与密度、抗压强度、弹性模量、软化系数呈负相关，与孔隙率呈正相关。所加入的两种材料可较好地模拟样品的多孔特性。在允许波动范围之内，所建立的模拟多孔含水岩层的经验方程可用于工程实践。

白鹤滩水电站左岸坝基边坡地质构造错综复杂，坡内多条陡倾断层和顺坡向错动带是控制边坡稳定性的主要因素。结合现场地质资料、监测数据及勘察认识，利用通用离散元软件UDEC建立左岸边坡开挖变形分析模型。首先在无支护开挖工况下研究主要结构面F17、LS331、LS3319上下盘岩体关键点应力及位移变化，并指出潜在失稳破坏区；同时采用预应力锚索单元在不同施工时序下对潜在不稳定块体进行加固计算，对比研究了不同方案下左岸边坡的变形破坏机制和整体稳定性。计算结果和现场监测表明：左岸边坡变形破坏与结构面密切相关，在开挖卸荷作用下主要表现为沿错动带的压剪破坏和陡倾断层的张拉破坏，采用预应力锚索加固处理能够有效提高边坡稳定性，并且在开挖后及时跟进支护可以有效抑制边坡剪切变形。研究结果对白鹤滩左岸坝基边坡后续工程加固和施工工序具有一定的参考价值。

为考虑位移因素对于滑坡渐进性发育的影响，基于简布条分法基本原理，提出一个新的能够考虑滑带土抗剪能力与剪切位移关系的滑坡稳定性简化评价方法。该方法从以下两个方面对简布条分法进行了改进：（1）将滑动面剪应力-剪切位移本构模型与沿滑动面的剪切位移模型相结合，以考虑滑坡随着变形发展渐进破坏的过程；（2）提出局部安全系数的概念并对整体安全系数进行重新定义，物理意义更加明确。推导了该方法的计算公式，给出了其迭代计算思路，并编制Matlab计算程序。通过计算实例表明，应变软化边坡的安全系数不仅取决于土的抗剪强度指标，还与其剪应力-剪切位移关系密切相关，基于峰值强度的极限平衡条分法应用于实际工程是偏于危险的。同时该方法还可以有效地模拟边坡从开始破坏到最终滑动面贯通的发育过程，实现了坡体位移与安全系数之间一一对应的定量关系，进而实现边坡的稳定性的预测预报。

地震动作用下小净距隧道双洞之间存在动力放大效应，依托某强震区小净距盾构隧道工程，基于傅里叶贝塞尔级数展开法，给出了平面P波垂直入射下地下半空间弹性介质中圆形双洞复合式衬砌洞室动应力解析解，深入研究了洞室间距、隧道全环注浆加固参数对衬砌动应力集中的影响。研究表明：洞室间距对衬砌的动应力有显著影响，随着洞室之间距离减小，左右洞室衬砌动应力明显放大；2.5倍洞径可以作为双洞洞室动力临界净距，当洞室间距小于2.5倍洞径时，应进行抗减震设计；注浆全环加固隧道周边土体，可以明显降低衬砌动应力，随着注浆加固区土体刚度增加，衬砌动应力幅值减小；结合工程实际，建议注浆加固区厚度取0.5倍的洞室半径值，加固区岩体和围岩剪切波速度之比取1.2。

工作面长度的增加是导致采场矿压强度增加的原因之一。通过相似模拟试验证明了覆岩能破断成梯形台形态，并指出了梯形台形成机制及梯形台的参数计算原理。利用梯形台分析了工作面长度对覆岩破断规律的影响，指出了梯形台与关键层理论的关系。基于薄板理论推导了岩层破断步距的计算公式，分析了工作面长度对梯形台空间结构演化及采场矿压显现特征的影响。结果表明：由于覆岩按梯形台破断，工作面长度较短时，加载层厚度较小，来压强度小。对比研究了赵庄2号井1305大采高工作面（工作面长度为85 m）和1302大采高工作面（工作面长度为180 m）矿压特征。结果表明：1305工作面基于梯形台和薄板理论进行选型计算，支架计算工作阻力为4 738 kN，现场选用额定阻力5 500 kN能够满足顶板控制的要求；1302工作面基于梁理论进行选型计算，支架计算工作阻力为7 623 kN，现场需选用额定阻力7 800 kN才能满足顶板控制的要求。

针对土坡稳定体系可靠度分析问题，提出了一种分析方法框架，包括采用缩减方差抽样技术生成随机变量样本值、采用全局优化算法搜索边坡最小安全系数、采用Monte-Carlo法计算边坡体系可靠度3个主要部分。在此框架下，建立了一种较为简便实用的高土石坝坝坡稳定体系可靠度分析方法。该方法采用拉丁超立方抽样技术生成随机变量的样本值，再用商业软件STAB搜索相应的坝坡最小安全系数，最后用可靠指标法或Monte-Carlo法计算坝坡体系可靠度。工程算例表明，筑坝材料强度参数的随机不确定性对坝坡临界滑弧的位置影响较大，坝坡稳定体系可靠度小于单一滑动面的坝坡稳定最小可靠度，提出的方法可用于实际工程坝坡稳定体系可靠度分析。

采前煤柱稳定性研究是工作面冲击危险性评估和开采方案设计的关键。以山东某矿深井巨厚砾岩条件工作面开采遗留煤柱为背景，采用案例调研、理论分析、数值模拟和工程实践等方法，对巨厚岩层-煤柱协同变形机制及其煤柱稳定性进行了研究，建立了巨厚岩层-煤柱协同变形的简化力学模型，探讨了引起煤柱变形的主要应力来源和变形形式，推导了在协同变形条件下煤柱的应力-应变关系。以此为基础，综合煤柱煤体应力、围岩稳定性和变形特征等条件，提出了煤柱整体失稳的力学判据。研究结果表明：巨厚岩层-煤柱失稳诱发冲击与煤柱的位置、尺寸和上覆岩层运动或变形关系密切，上覆岩层运动或变形是诱发煤柱失稳的动力因素；巨厚岩层-煤柱的变形主要包括受集中力F压迫的协同挠曲压缩变形和受集中力G作用的重力沉降变形，二者保持内在协调性；巨厚岩层下煤柱整体失稳的工程判据为煤柱煤体平均支承应力p超过其平均极限支承强度Rc（p≥Rc）；评估得到遗留的50 m煤柱具有强冲击危险性，并通过优化开采设计，取得了良好的效果。该研究成果对相似条件煤柱留设及其稳定性分析具有参考意义。

为进一步改进和完善空心包体应变计地应力计算方法，基于地应力现场实测过程中空心包体应变计的安装方式，采用线性参数的最小二乘拟合方法对地应力分量计算公式进行了推导，得出了6个地应力分量的改进算法及其标准误差的计算公式。采用实现完全温度补偿并考虑岩体非线性的地应力测量技术，对弓长岭井下矿-160、-220、-280 m 3个水平3个测点进行了地应力现场实测，得到了3组孔壁应变数据。分别使用常规算法和改进算法对3个测点的地应力分量进行计算并分析其标准误差。结果表明：采用改进算法计算得到的地应力分量其标准误差普遍小于常规算法计算地应力分量产生的标准误差，说明改进算法比常规算法具有更高的可靠度。在此基础上，对弓长岭井下矿地应力赋存规律进行了研究。结果显示：矿区地应力场以水平构造应力为主导，最大水平主应力的走向总体上为南东东-北西西方向，最大水平主应力、最小水平主应力和垂直主应力均随深度呈增长关系。

海上风力发电作为一种绿色能源越来越受到重视，单桩基础局部冲刷的研究受到广泛关注。建立了不同冲刷深度的海上风机系统模型，对各模型进行静力、动力分析。结果与前人冲刷模型试验吻合，随着桩基冲刷深度的增大，支撑系统的最大水平位移、最大水平应力和最大弯矩值逐渐增加。冲刷深度为两倍桩径时，最大水平位移、应力及弯矩分别比无冲刷条件下增加了3.61%、12.7%、10.3%。风机系统的自振频率随冲刷深度的增加而降低，高阶频率更易受影响。桩基冲刷对风机系统的动力特性有显著影响，塔顶水平位移时程曲线随冲刷深度不断上移。海上风机设计时必须考虑冲刷对支撑系统承载性能的影响。

岩体中普遍存在着断层﹑节理和裂隙等结构面，这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此，研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法，模拟裂隙时独立于网格，因此，在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究，建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型，并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。

裂隙渗流会引起裂隙周围岩体中的温度场变化，在低温岩体中其影响更为明显；此外，裂隙水与周围低温岩石介质发生热交换会引起裂隙中的水冰相变过程发生，而裂隙水冻结将阻碍裂隙渗流，引起裂隙渗流场的变化。因此，低温下的裂隙岩体水-热相互作用是一个强耦合过程。考虑裂隙中的水冰相变过程和渗流作用，建立了低温冻结条件下裂隙岩体水-热耦合模型；以冻结法施工为例，考察了低温冻结过程中裂隙水渗流对裂隙冻结交圈的影响。研究结果表明：由于裂隙渗流的存在，距裂隙较远处岩石先冻结，裂隙冻结所需时间远大于周围岩石；裂隙宽度和裂隙水压力差都会影响冻结交圈时间，裂隙越宽、水压力差越大，裂隙冻结需要时间越长；随着冻结时间的推进，裂隙水渗流速度逐渐降低，当裂隙冻结后裂隙渗流停止。最后通过构建随机裂隙网络模型，利用所建立的水-热耦合模型考察了裂隙网络渗流对冻结交圈的影响，说明了在冻结法施工中考虑裂隙的重要性。

建立尾砂的本构模型是开展尾矿坝数值模拟和安全评价的重要基础，而目前尾砂的本构模型研究多集中于非线性弹性模型，如Duncan-Chang模型，关于弹塑性模型的研究较少。结合尾砂的应力-应变特性，提出了一种适用于描述尾砂力学特性的改进广义塑性模型。基于用户自定义材料子程序UMAT，将提出的模型在ABAQUS中二次开发实现，应力积分采用Runge-Kutta显式积分。通过三轴试验模拟验证，偏差应力曲线表明，有限元计算结果可以反映围压对应力-应变关系的影响，抗剪强度随应变逐渐硬化，达到峰值强度，随后发生应变软化，抗剪强度有所下降。体应变曲线结果表明，广义塑性模型可以很好地描述体应变曲线的体缩-体胀发展过程，与试验曲线吻合较好。同时有限元数值模拟结果和理论值误差很小，且和试验结果拟合较好。该研究成果可进一步用于尾矿坝的应力变形计算和安全评价。

浮力计算是地下结构抗浮设计的关键，而弱透水黏性土中浮力大小的确定一直存在争议。对此，设计一套免受摩擦力影响的模型试验装置，对埋置于软黏土中的地下结构模型进行浮力测试。试验采用卸重上浮的思路，力学模型简单明确，浮力即等于浮起瞬间的剩余配重，并通过水压和位移的实时监测来判定结构上浮的临界状态。试验结果表明，软黏土中的水压是一个动态变化量，受渗流、温度、卸重回弹以及土-水相互作用的影响，且往往小于理论水压；实测浮力是理论浮力的0.93倍，存在少量折减；水压和位移时程曲线显示上浮过程是持续动态变化的，最终浮起时具有突发性。

应用光纤传感技术测量模型相似材料内部应变，必须解决光纤与相似材料耦合的问题。传统的研究中普遍存在应变传递理论分析结果无法通过试验检验的问题。设计实现了一种基于光纤-胶体-相似材料3层应变传递结构的光纤Bragg光栅（FBG）相似材料一维传感器。通过对该传感器应变传递特性和标定方法的研究探索解决光纤与相似材料耦合性的问题。该传感器采用无涂覆层的光纤Bragg光栅作为传感元件，将光纤埋入相似材料应变试件内部直接测量轴向中心应变。通过应变传递理论分析、数值模拟以及全新设计的标定试验研究3种方法研究了该相似材料光栅传感器的应变传递特性。平均应变传递率结果显示3种方法所得到的数值非常接近，表明该传感器的光纤、胶体、相似材料具有良好的耦合性以及应变测量的准确性。该方法为解决模型相似材料内部三维应变测量提供了可行方案。

一系列深埋长大隧道工程在我国西部涌现，深埋隧道围岩的岩爆成因机制及监测预警成为亟待解决的课题。微震监测作为岩体微破裂空间监测技术在隧道工程岩爆预警上发挥了重要作用。由于微震数据存在的波动性（时间序列）、离散性（空间分布）和误差，利用震源参数随施工循环的演化特征评估岩爆发育过程并非足够稳定、有效。基于3种微震监测易于获取且紧密联系的震源参数（地震能量、地震矩和视应力），扩展建立了用于岩爆评估与预警的EMS方法。该方法定义了地震能量、地震矩、视应力3种评价指标，并获得了3种评价指标间的关系；将3种评价指标用于岩爆预测，并进行了岩爆等级评估；同时还给出了微震路径与震源参数空间的概念。最后，将EMS方法应用于巴陕高速米仓山隧道的岩爆灾害实例，对岩爆里程段微震事件簇（微震事件聚集区域）内岩体累进破裂和岩爆发育过程进行了评估，实现了简易有效的微震数据日常处理与岩爆灾害评估预警。