不同于地上混凝土结构，隧洞衬砌结构会受到碳酸性水环境的影响，其劣化机制与大气环境中碳化作用存在明显的差异。为揭示碳酸性水环境作用下隧洞衬砌结构的劣化机制，开展了大气环境和碳酸性水环境下水泥砂浆室内加速侵蚀试验研究，研究结果表明：碳酸性水环境中水泥砂浆的碳化深度明显小于大气环境中的碳化深度，其碳化系数仅为大气环境中的0.18倍；大气环境下碳化作用导致水泥砂浆中孔体积减小，其单轴抗压强度和弹性模量均有明显提高，而碳酸性水环境加剧Ca元素流失，其碳化作用导致碳化层中Ca元素含量较低，从而破坏了水泥砂浆孔隙结构，导致水泥砂浆孔体积增加，最终碳酸性水环境导致水泥砂浆单轴抗压强度和弹性模量均降低，其中弹性模量降低幅度较大。研究结果揭示了碳酸性水环境会加速隧洞衬砌结构中Ca元素流失，造成孔隙增加而强度降低，最终导致隧洞衬砌结构劣化，从而影响隧洞衬砌结构长期安全运行。

对非饱和土降雨入渗过程及其致灾机制的深入认识有赖于室内外试验及多物理量联合监测。基于离心机超重力环境下土体响应监测技术，开展了非饱和黄土地基降雨入渗离心模型试验，测试研究了超重力环境对新研制的微型TDR探针及张力计测量的影响规律，并利用TDR、张力计及弯曲元对降雨入渗过程中土体响应进行多物理量联合监测。研究结果表明：在不同离心加速度下微型TDR探针与给定含水率土体实测原始波形重合，说明超重力环境对TDR测试没有影响，含水率的测试误差在2%以内。在离心机加速过程中，张力计所测吸力下降约2.0～2.9 kPa，当离心加速度稳定在40g时，所测吸力在10 min内上升并接近常重力下土体初始吸力。在降雨入渗过程中，埋设在同一深度的TDR探针、张力计和弯曲元对湿润锋响应的时间点基本一致，降雨入渗导致土体含水率增加，基质吸力降低，剪切波速降低。这些多物理量监测数据有助于建立非饱和土含水率-吸力-剪切模量之间的关系。

隐晶质玄武岩是白鹤滩水电站超大规模地下厂房洞室群中广泛揭露出的主要围岩类别之一，且普遍赋存隐微裂隙，开挖过程中展现出特殊的力学响应特征。综合采用CT扫描、高速摄像、声发射、扫描电镜等手段，分析了隐晶质玄武岩的变形破坏特征、裂纹扩展及声发射演化规律等，并探讨了其细观破坏机制。结果表明：完整试样应力-应变曲线光滑，加载至峰值强度时瞬间爆裂，声发射信号异常集中，碎裂破坏；含隐微裂隙试样应力-应变曲线呈锯齿状，出现多次表面剥落，整个过程声发射信号非常活跃，以劈裂破坏为主；含宏观裂隙试样应力-应变曲线呈双峰或多峰状，峰前原生裂隙局部滑移，第1次峰值时宏观破裂面形成，声发射信号集中于峰值应力跌落瞬间，以剪切破坏为主。隐晶质玄武岩细观破坏机制主要为矿物颗粒沿晶断裂和穿晶断裂。研究成果为准确认识并科学掌握白鹤滩地下厂房洞室群围岩的力学响应及其破裂演化特征奠定了坚实基础，也可为硬脆岩体高应力破坏的认识和灾害控制提供参考和借鉴。

通过开展饱水花岗岩单轴压缩声发射试验，采用快速傅里叶变换提取声发射信号主频及次主频，进一步提出了综合表征声发射信号频谱特征的参数：主频比F，即主频与次主频的比值，研究岩石破裂过程声发射频谱特征。研究结果表明：饱水花岗岩破裂过程声发射主频、次主频频带数量呈现“3→6→3”先增大后减小的规律，与主频分布特征相比，次主频低频区间占比减小，中、高频区间占比增大，次主频总体平均值比主频高约5 kHz。声发射主频比F分布区间在弹性阶段达到最大：0

研究边坡在反复地震作用下的累积损伤对于库区岩质边坡的长期稳定性分析具有重要意义。结合振动台试验与理论分析方法，对水平层状岩质边坡在反复地震作用下的累积损伤与破坏模式进行了深入研究。结果表明：采用混凝土疲劳裂纹扩展的双K准则可以很好地解释边坡的累积损伤曲线的变化趋势，完整的累积损伤演化曲线表现出初始损伤、细观裂纹扩展、宏观裂缝扩展3个阶段；分别采用三次多项式和幂指数建立了微小地震和强震作用下边坡岩体累积损伤演化模型并取得了良好的拟合效果，微小地震作用下边坡岩体累积损伤曲线表现出初期略微下降、中期线性增长、后期平缓增长的正“S”型的3段式特征，强震作用下边坡岩体累积损伤相比微震作用下有较为快速的增长，且在边坡破坏前增速加快。频发微震作用下，水平层状岩质边坡在坡体上部形成阶梯型破坏面并向坡面出露，在坡体下部形成平直的剪切破坏面，坡体最终的破坏面呈阶梯状。

抽采负压是影响钻孔孔周煤体瓦斯渗流规律的重要因素之一，为了深入研究负压变化对煤体瓦斯渗流特性的影响，通过设计瓦斯渗流特性相似模拟试验平台对孔内瓦斯流动规律进行测定，并运用数值软件对孔周煤体内瓦斯渗流状态进行分析，结合上述结论对钻孔孔周煤体瓦斯渗流特性进行研究。结果表明：抽采负压提供孔周煤体瓦斯向钻孔内渗流的动力，负压值设定为25～35 kPa时，可保证雷诺数集中分布在10～20之间，其低速紊流状态有利于瓦斯高效抽采；随着抽采时间的增加，负压对瓦斯的引流作用逐渐减弱，瓦斯流量随时间呈现负指数衰减规律，煤体渗流速率饱和值处于7.41×10?6～1.30× 10?5 m/s之间；煤体内部黏滞阻力对瓦斯渗流存在抑制作用，随着抽采负压的增大，孔口负压上升趋势比孔周煤体内负压上升趋势更快，说明抽采负压变化对钻孔孔口负压的影响更为显著。

基于微生物诱导碳酸钙沉淀作用（MICP）的土体改性技术近年来在岩土工程领域引起了人们的广泛关注。该技术在改善岩土材料的强度、刚度、抗液化、抗侵蚀及抗渗透性等性能的同时，还能维持土体良好的透气性和透水性，改善植物的生长环境。由于微生物矿化作用涉及一系列生物化学和离子化学反应，固化过程中的反应步骤较多，因此，MICP固化效果受许多因素的制约与影响。基于大量文献资料，系统总结了细菌种类、菌液浓度、温度、pH值、胶结液配比及土的性质等关键因素对微生物改善岩土材料性能的影响，讨论了这些影响因素的优化方式和未来的研究方向，主要得到了以下几点结论：菌种类型、菌液浓度、温度、pH、胶结液性质会从微观上影响碳酸钙的晶体类型、形貌和尺寸，进而在宏观层面影响岩土体的胶结效果；菌液浓度尽可能高、温度在20～40℃间、pH值在7.0～9.5左右、胶结液浓度在1 mol/L以内的因素条件对微生物加固岩土体具有较好的效果。上述范围内的低温、较高的pH值、低浓度胶结液有助于提高土体的抗渗性，而高温、较低的pH值以及中高浓度胶结液有助于提高土体的强度；MICP加固土体的有效粒径范围为10～1 000 ?m，相对密度越大、级配越好则加固效果越好。分步灌浆法、多浓度相灌注法及电渗灌浆法有助于提高土体固化均匀性，0.042 (mol/L)/h以下的注浆速度有利于提高胶结液利用率，砂土试样的灌浆压力一般在10～30 kPa之间，粉黏土试样的灌浆压力不宜超过110 kPa，过高的灌浆压力会破坏土体结构，降低固化效果。

为合理反映颗粒破碎对堆石料力学特性的影响，基于试验结果分析，得出了堆石料在压缩和剪切作用下的颗粒破碎特性规律。通过引入压缩破碎和剪切破碎的相关参数，借鉴已有本构模型的合理定义，吸收临界状态理论和边界面理论的优点，发展了考虑颗粒破碎和状态相关的堆石料静动力统一弹塑性本构模型，并阐述了模型参数的确定方法。该模型不仅能够反映堆石料在静力荷载作用下的低压剪胀、高压剪缩、应变软化和硬化等特性，还能够反映在循环荷载作用下应力-应变的滞回特性和残余变形的累积效应。最后为验证模型的合理性，分别对堆石料的静力三轴和循环三轴试验进行了数值模拟预测，结果表明：模型预测与试验数据吻合良好，所提出的本构模型能够合理地描述颗粒破碎对堆石料静动力变形特性的影响。

颗粒分析试验被广泛应用于土木、水利水电和采矿等诸多工程中，以确定材料的颗粒级配。但对于工程中常见的超大粒径填料合理取样质量问题，尚缺少可依据的文献资料。将试样选取视为一种分层抽样过程，结合抽样分析理论，建立了能够确定样本容量的理论公式；考虑到颗粒材料平均粒径与其级配分布的关系，提出了一种为确定材料整体级配所需试样质量的通用方法。对不同级配分布情况开展了针对性试验，验证了所提出方法的适用性。并将该方法应用于确定重庆江北机场和贵州龙洞堡机场填料的颗粒级配，取得了良好效果。所得结果可为相关工程实践及规范的发展和完善提供参考。

玄武岩纤维复合材料（BFRP）锚杆与传统钢锚杆相比具有比强度高、耐腐蚀性强、与围岩协调变形性好等优点，是一种新型高性能纤维锚杆，在边坡加固领域的应用才刚刚起步。通过BFRP锚杆加固黄土边坡的现场拉拔试验，较系统地研究了BFRP锚固体系在不同锚杆直径、锚固长度下的工作性能，并通过现场开挖式剖析，分析了BFRP锚固体系的破坏模式。试验结果表明，破坏模式受控于锚固系统诸界面的相对强度，?12 mm和?16 mm锚杆体系为锚杆与灌浆体界面（第1界面）剪切破坏，?25 mm锚杆体系为灌浆体与土层界面（第2界面）滑移破坏；一定锚固条件下，增大锚杆直径可显著提高锚固体系的极限抗拔力；随着锚固长度的增加，极限抗拔力并非始终线性增大，而是增幅逐渐减弱，存在临界锚固长度；第1界面和第2界面平均黏结强度均随锚固长度的增大而减小，并给出了诸界面平均黏结强度的建议值，可供实际工程设计使用；杆体轴力沿锚固深度逐渐衰减，分布形态与受拉荷载大小、锚杆直径和锚固长度等有关；锚杆界面摩阻力分布服从随锚固深度先增大后减小的单峰形态，峰值多出现在锚固前端0.5 m范围内，同样受锚固长度和直径影响。建议今后进一步改善BFRP材料的抗剪性能以及BFRP锚杆表面形态设计和制作工艺。

针对深大基坑支护工程中存在的问题提出了一种新型的支护形式-剪力键支护体系，由竖直桩和斜向支撑的剪力键构成。通过模型试验，对剪力键组和单排悬臂桩组分别量测了不同开挖深度时的桩顶水平位移及桩身内力，同时进行了有限元的数值模拟。结果表明剪力键组的桩顶最大水平位移约为悬臂组的1/10，桩身最大弯矩值约为悬臂组的一半左右。剪力键支护体系具有更大的刚度，能够有效减小桩顶位移及桩身最大弯矩，且耗材相对较少，支护深度较大，占据空间较少，不影响主体地下室施工，可节约工期和成本，为深大基坑的支护提供了新的型式。

土体地震反应分析对结构的地震反应和抗震安全性评价具有重要的意义，以大型振动台模型试验为手段，研究一致激励和非一致激励下自由场土体的非线性地震反应规律及其影响因素，能够为结构的地震破坏机制分析提供支撑。通过对场地土体的动力特性、加速度响应、剪应力-剪应变曲线以及土体沉降的分析，研究了不同地震动、不同地震强度作用下自由场土体的地震动反应特性及其变化规律。研究结果表明：土体的非线性发展程度不仅与地震动记录有关，还与地震动输入方式、加载等级有关；在大震或纵向非一致激励作用下，土体运动的不一致性导致了更强的土体结构性变化，频率降低，阻尼比增大，土体刚度弱化，土体的非线性发展相对较快；土体剪应力-剪应变曲线和沉降曲线的变化规律一定程度上也反映了土体的塑性发展情况。所得结论和自由场地震反应宏观现象一致，彼此佐证了结论的合理性。

为探讨摩阻力增强效应影响下的水平受荷桩承载机制，建立了水平荷载作用下考虑增强效应影响的桩侧摩阻力?-s曲线模型，并在此基础上推导了桩身单位长度侧阻抗力矩数值解。随后开展了桩径、临界位移、极限侧摩阻力等因素影响下Ms,ini /Msu,ini-? /? ref,ini无量纲化曲线规律分析。结合摩阻力增强效应理论进而分别建立了线弹性-塑性、双曲线、API规范砂土与黏土? -s曲线模型作用下的桩身侧阻抗力矩简化理论解，并基于传递矩阵法原理，解得了考虑桩侧摩阻力作用时基桩水平承载力半解析解。通过多组验证工况和工程实例分析验证了所建立的侧阻抗力矩简化解以及考虑侧阻抗力矩影响的基桩水平承载力方法的正确性，同时也证明了摩阻力增强效应对基桩水平承载力的影响不能忽略，并得出如下结论：无试验数据参考时，砂土、黏土与桩接触面的刚度指数可分别取0.725和0.600；桩身侧阻抗力矩随着桩径、极限摩阻力的增加以及临界位移的减小而增加，其中桩径对侧阻抗力矩影响更为明显。

在寒区隧道工程中，冻胀的水分迁移作用加剧了隧道围岩的冻融破坏。耦合水分迁移作用的岩石本构模型对防治冻融渗漏、崩塌等隧道围岩灾害有着重要意义。基于内状态变量理论，将冻胀过程中的水分迁移量作为本征变量引入Helmholtz自由能，在热力学框架下建立了一个岩石力-热-水耦合本构模型。模型描述了温度和水分迁移对冻融后岩石损伤阈值、等向强化饱和值、等向强化速率等力学参数的影响。模拟了岩石冻融后力学性质的劣化。区别于全量经验公式，本模型以增量形式给出，为复杂应力历史条件下的数值模拟提供了便利。通过将模型模拟曲线与冻融后岩石常规三轴压缩试验曲线进行对比，初步验证模型的可靠性，为实际寒区工程的冻胀破坏预测提供参考。

取3类具有明显片理面的武当群片岩，采用偏光显微镜、扫描电镜研究其微观组构、断面形貌特征。构建了面积百分比、尺寸与周长-面积分形维数、定向分布系数等指标，利用IPP专业图像处理软件，对矿物和孔隙（微裂隙）两系统的含量、形态、分布特征等进行定量分析；并基于Matlab的数字图像处理技术，对断裂面粗糙度展开定量评价。偏光显微镜下显示，片岩具有矿物定向排列与粒、片状矿物近互层状分布的典型特征。矿物系统的微观定量分析结果表明：武当群片岩白云母的含量为11%～30%；小岛法计算的白云母矿物周长-面积分形维数为1.31～1.39；定向分布系数为0.61～0.85；微裂隙依附于白云母矿物边缘孕育、扩展，导致武当群片岩的孔隙（微裂隙）微观指标与对应的矿物微观指标有较强的相关性。沿片理面的片岩断面主要呈现沿晶破坏特征，断口较平坦，基于断面形貌图绘制的灰度信号波曲线的盒维数为1.48～1.60；斜交片理面破裂的片岩断面以穿晶破坏为主，有穿晶-沿晶破坏交替而生的迹象，断口粗糙起伏，对应的灰度信号波曲线的盒维数为1.65～1.69。武当群片岩的断面粗糙度各向异性约为1.03～1.14。

TBM掘进经常遇到复合岩层，在该地层中施工对滚刀破岩非常不利。为了提高掘进效率，降低工程造价，开展复合岩层滚刀破岩效率研究很有必要。鉴于此，采用滚刀岩机作用综合试验台对砂岩、花岗岩复合（复合比例为4? 6）而成的岩层进行5组刀间距下3种贯入度的3把滚刀同步旋转切割试验；试验过程中监测破岩总法向力、总扭矩，并分别收集两种岩石的岩渣进行筛分、称重；依据试验结果分析不同刀间距、贯入度下的法向力、扭矩及比能的关系。研究表明：不同刀间距时，平均法向力、平均扭矩随贯入度的增加而增大但两者的增加趋势不同，即平均法向力随着贯入度的增加几乎呈直线增大，而平均扭矩随着贯入度的增加其递增趋势减小；不同贯入度下存在不同的最优刀间距；砂岩、花岗岩切割轨迹长度比例为4：6时刀间距与贯入度的比值在14左右，破岩效率最高。

大型水力压裂是干热岩地热能开发中人工储留层建造的最有效手段，其核心力学问题为高温、高压下岩石的水力破岩机制。通过单轴应力下带钻孔花岗岩注入高温蒸汽破坏试验，研究固-热耦合作用下花岗岩的水力破岩机制。结果表明：高温对花岗岩破裂有很大的促进作用，热效应导致强度弱化，降低破裂压力。高速率注入430℃和350℃蒸汽破坏试验中，破裂压力比常温水压裂至少降低58%；低速率注入400 ℃和450 ℃蒸汽破坏试验中，花岗岩破裂压力比常温水压裂降低75％。注蒸汽破坏过程可分为热破裂损伤和宏观裂缝扩展两个阶段。高温蒸汽产生的热应力在钻孔周围随机发生热破裂，随着注入蒸汽时间的增加，热破裂范围由钻孔附近逐渐向远处扩展，热破裂分布密度增大，为宏观裂缝的产生提供便利条件。初始宏观裂缝首先出现在钻孔两侧，沿着最终形成的宏观裂缝轨迹扩展，直到试样破坏。与常温水压裂相比，低速率注蒸汽破坏是一个缓慢的延性拉破坏过程，裂缝相对钻孔不对称扩展，宽度小于水力压裂裂缝宽度。

岩体破坏实际上是在能量驱动下的一种状态失稳现象，其变形破坏过程中的能量演化规律能够反映岩体变形破坏的本质特征。为探求煤样破坏过程中能量的演化规律，对其进行了常规三轴及不同卸围压速率的三轴卸围压试验。试验结果表明：不同试验路径下，煤样在峰值应力前以能量的存储及耗散为主，在峰值应力后以能量的释放及耗散为主；煤样在较高的卸围压速率下，卸荷段经历时间较短，内部裂纹扩展及环向变形不充分，耗散能和克服围压做功耗能都较小，大量的能量在破裂时释放出来，破裂更加剧烈；同时，卸围压速率越快，煤样在破裂前后损伤量增加速率越快，损伤曲线急剧增高，因此，高应力煤岩体在开挖卸荷时，在岩体开挖施工速度较快时，极易引起大量积聚能量的突然释放，造成严重的灾害。

含水合物土是由水合物赋存在土体中形成的特殊土，其力学特性不仅与饱和度有关，和赋存模式也有很大的关系。不同赋存模式的水合物胶结结构不同，其力学特点也不尽相同。考虑含水合物土在荷载作用下胶结结构会发生破坏，采用扰动状态概念，建立含水合物土的应力-应变模型来描述胶结结构受载破坏的全过程。同时，采用水合物的有效饱和度表示的伤变量来描述含水合物土的内部结构的在受载条件下的发展规律。根据含水合物土的力学行为，采用理想弹塑性模型和摩尔-库仑模型来描述两个参考状态，形式相对简单，参数易于确定。最后通过三轴剪切试验对模型进行了检验，结果表明该模型能较为合理的模拟含水合物土在扰动过程中的变形特点。

针对混合粒径钙质砂中不同粒径颗粒绝对破碎量无法获得和现有破碎率难以考虑破碎重叠掩盖破碎量这两个问题开展研究。设计了粗砂、中砂、细砂颗粒集中分布的3种级配钙质砂试样，进行侧限压缩试验。对不同粒径区间钙质砂分别染成不同颜色，拍照获取各粒径区间钙质砂破碎信息；采用Image J软件进行彩色图像颗粒分割、二值化处理、统计各颜色颗粒面积，换算得各颜色颗粒破碎后含量；并提出考虑破碎重叠掩盖的试样累积破碎率指标Ba。结果表明，随压力增大及颗粒分布集中，试样的重叠掩盖破碎量增大。混合粒径钙质砂中的中间粒径（0.25～1.00 mm）颗粒易于破碎，各粒径颗粒破坏模式以颗粒边角破碎为主；累积破碎率Ba值较相对破碎率Br较大，与垂向压力对数值间满足线性关系，为颗粒破碎研究提供了新的思路。

为了研究岩石在加载-卸载过程中的应力-应变关系，以砂岩为例，对其进行常规三轴加卸载试验。分析了峰后卸载阶段岩石的非线性特性，对岩石的损伤变量进行定义，给出了峰后卸载过程中用于描述应力-应变关系的弹性模量模型。通过分析加载-卸载过程中的轴向应变与径向应变的关系，得到了卸载过程中泊松比模型。引入D-P塑性模型，针对砂岩的塑性硬化特性，对硬化函数进行修正，建立了与等效塑性应变相关联的损伤模型。将计算模型矩阵化后进行数值计算。在此过程中得到如下结论：多孔隙岩石在加载过程中表现出明显的非线性特征，随着体应力的增大，岩石的弹性模量逐渐增大。岩石峰后卸载过程中，当轴向应力大于围压时，应力-应变可以利用峰前弹性阶段的弹性模量模型乘以连续性因子进行描述。随着等效塑性应变的增大，泊松比先增大后减小，最终趋于稳定。峰后卸载过程中，等效塑性应变不发生变化，此时泊松比保持不变。利用提出的本构模型进行了数值计算，数值计算结果与试验结果进行对比，结果表明，提出的模型能够反映出岩石在峰后卸载过程中的应力-应变规律。

常规室内模型试验中，土体侧向位移作用下桩身变形特性能够获得，而桩周土体内部的位移很难量测。因此，结合透明土材料、粒子图像测速（PIV）技术和光学测量系统，设计了侧向位移作用下被动桩桩周土体变形测量试验系统。通过非插入式测量，得到桩不同埋置深度下桩周土体的位移场。试验结果表明：相同侧向位移模式下，桩埋置深度越大，桩对桩前土体的遮拦效应越明显，中间与两侧土体的相对位移越大。基于模型试验的结果建立了数值模型，参数敏感分析表明土体侧向位移形状对桩前土体位移遮拦效应的影响程度依次为矩形、抛物线形和三角形，且遮拦效果随着桩径的增加而更加明显。

确保海床场地的动力稳定性是近海工程安全运行的前提。目前针对复杂应力路径对海域环境饱和粉土动力学行为特性的影响研究尚属少见。利用GDS空心圆柱扭剪仪，开展了轴向-扭转耦合循环加载的重塑饱和粉土不排水试验，模拟波浪的波幅大小和海床深度变化，以广义剪应变?g = 5%为液化标准，研究了均等固结条件下循环应力路径（循环应力比CSR和循环加载幅值比?）对饱和粉土不排水动力特性的影响。结果表明：CSR与? 值的大小对饱和粉土的超静孔压与变形发展特性影响显著，圆形应力路径（? = 1）循环加载时饱和粉土最易液化；当CSR ≤0.050时，饱和粉土不会发生液化；当CSR≥ 0.065且? = 1时，饱和粉土会发生液化；CSR > 0.150时，粉土易发生液化。?g与孔压比ru的相关性受CSR与? 值大小的影响较小，且?g可表示以ru为变量的正切函数。以等效循环应力比ESR作为表征复杂应力路径下动应力大小的指标，饱和粉土达到?g = 5%所需的液化振次NL与施加的ESR值具有事实上的唯一性关系，ESR随NL的增大而减小。

由于煤矿再生岩体围岩结构和成分进行了重组，其力学特性将发生本质变化。根据再生岩体压实成岩过程中破碎岩体多次受压的实际情况，采用广西州景煤矿再生顶板破碎围岩进行有侧限条件的分次压实试验。结果表明：分次压缩过程可分为压实破碎、压实应力记忆和压密固结3个阶段。在第1次压缩时，压实特性曲线呈线性变化，后续压实特性曲线逐渐向指数形式转变，而“压实应力记忆”阶段表现为多项式变化。在侧限压缩条件下，提出以最小势能原理作为判定不规则岩样挤压破碎判据。破碎岩样压缩特性受压实次数和装岩高度的影响，破碎岩样分次压实次数越多，累积压缩应变越大，后续压缩时的应力增量越大。装岩高度越大，压缩至相同压应力时，累积总应变越大。压实应力达到一定程度，不同粗粒径破碎岩样的粒径分布规律趋于一致，据此建立了破碎岩样粒径分布与工程中破碎岩体压实程度之间的相关关系。

针对层状地基，采用改进的比例边界有限元算法（SBFEM）与有限元算法相结合，求解层状地基内刚性埋置管道动力响应问题。改进的比例边界有限元算法引入一个新的坐标变换关系，以一条相似轴代替传统SBFEM的相似中心，并利用加权余量方法得到层状地基动力刚度矩阵方程，最后与有限元法在远场和近场交界面进行耦合，由此求得层状地基埋置管道的动力响应。数值算例验证了算法的准确性，并利用该算法分析了层状地基的非均质特性、泊松比以及埋置深度对管道动力响应的影响，以期为工程实际提供必要的数值依据。结果表明：地基的非均质特性和泊松比对埋置管道的动力响应均有显著的影响；随着埋置深度的加大，管道在相同荷载作用下的动力响应均减小。

确定统一的表征单元体（REV）尺寸是研究土体细观孔隙结构时首先需要解决的问题。使用6.5 μm分辨率同步辐射显微CT扫描南京粉细砂试样，从土样三维重构模型的5个代表性部位提取5组立方体孔隙REV，对每个REV采用最大球算法分析，建立孔隙网络模型，从中提取孔隙率、单位体积孔隙数、孔隙平均体积、孔隙最小体积、孔隙最大半径、孔隙最小半径、孔隙平均半径、孔隙截面积平均形状因子等8个孔隙结构参数，建立其与REV尺寸间的相关性。利用假设检验T检验和F检验，最终确定样品孔隙结构参数的统一REV边长尺寸为400像素，即2.60 mm。该方法可用于砂土、粉土等颗粒土体细观孔隙结构分析。

针对弹性波三维散射和动应力集中问题，提出一种基于球面波势函数的基本解方法（SWP-MFS）。方法基于单层位势理论，采用膨胀波和矢量剪切波球面波势函数构造散射波场，根据边界条件建立边界积分方程并配点求解。精度检验表明，该方法具有良好的数值精度及稳定性。以无限空间中三维夹杂体及空洞对平面P、SV波的散射为例，进行方法展示，并揭示了三维夹杂体周围弹性波散射的一些重要规律。结果表明：三维夹杂体随其内部介质刚度降低，位移谱曲线震荡越加强烈。三维球形空洞在P波和SV波水平入射下应力集中规律不同，前者在顶部和底部及附近更明显，后者在纵截面两45°角线附近更明显。与以往的集中力源函数相比，新的波场构造基本解更为简洁易用，为三维弹性波动分析提供了一种新型无网格数值方法。

提出了一种考虑地震波散射效应的沉积河谷空间相关多点地震动模拟方法。基于弹性波传播理论，采用合理功率谱、相干函数和传递函数模型，借助原型谱表示法得到沉积河谷空间相关多点地震动，并验证了方法的可行性。其中，传递函数应用动力边界元法计算，可在很宽的频带内精细模拟沉积河谷中波型转换及聚焦效应。针对常见V型河谷开展了计算模拟，结果表明，河谷地形及松软沉积层对地震波传播具有显著影响，河谷地表地震动空间差异性十分显著，一维模型难以考虑地震波的散射及聚焦效应，从而会明显低估沉积内部部分区域的地震动峰值加速度。沉积河谷中大跨结构抗震设计宜采用考虑地震波散射效应的空间相关多点地震动输入。

岩土工程中的不确定性是工程风险的重要来源。在我国7个主要地区的20个工程场地上开展了国内外首次大范围现场剪切波速不确定性专项试验，共由47家单位用11种工程常用仪器采用单孔法完成，获取了600组土层剪切波速Vs实测数据。依据实测结果拟合了场地时间平均剪切波速Vs, z的变异系数COV与计算深度Z的对应公式，揭示了土层剪切波速不确定性对场地刚性判断的影响，进而分析了我国场地分类的可能潜在误判区域。结论如下：在目前我国现场剪切波速测试水平下，Vs, z的变异系数与深度负相关——在近地表处达到极大值15%，并随着深度增加明显减小；地表至10 m深度区间内，Vs, z的变异系数迅速减小；10～20 m深度区间内，Vs, z的变异系数减小趋势变缓；超过20 m深度后，Vs, z的变异系数小于5%且基本不变；我国III、IV类场地分类采用的等效剪切波速Vse和欧美场地分类指标Vs, 30的变异系数基本相等；依据我国建筑抗震设计规范GB50011－2010，对I类场地以及覆盖层厚度小于5 m或等效剪切波速接近500 m/s的II类场地分类时需考虑剪切波速不确定性的潜在影响，其他类型的II类场地、III类和IV类场地在分类时基本可忽略其影响。

分析国际上现有液化场地剪切波速473组数据进行参数特征后，以此为基础，采用分区方法，利用成熟的Logistic模型，提出了以地表峰值加速度、剪切波速、地下水位、可液化层埋深等4个参数表达的土体液化概率计算公式和不同概率下液化临界值计算公式，研究了不同概率水平下公式的适用性，并与现有主要方法进行比较。研究表明：地震动强度为液化判别首要影响参数，液化层与非液化层剪切波速区分度不显著，采用以往一个公式构造液化判别公式的方式难以要达到基本要求；现有Andrus公式和石兆吉公式会将很多明显为非液化的场地误判为液化场地，违背了现有认识，达到了不可接受的程度，需要改进。文中公式取50%概率时液化点和非液化点回判成功率符合对等原则，不同烈度下成功率均在70%左右；总体看，公式表现良好，可为中国工程建设提供一个合理、可操作性强的剪切波速液化概率判别方法。

长春地铁1号线卫星广场站下穿轻轨3号线工程中，下穿段负一层顶板与轻轨隧道底板密贴接触，形成“叠合板结构”。施工中，先注浆加固下穿段土体，然后采用六导洞PBA法配合千斤顶顶升支撑修建下穿段，具体方法为：于1#+4#导洞冠梁上布置液压千斤顶，对上层板形成两端临时支撑；再于3#+5#、2#+4#导洞内各施作一组钢管混凝土立柱及纵梁，对下层板形成永久支撑，纵梁兼做负一层顶板。然后逐步移除千斤顶、拆除导洞初支、浇筑负一层顶板，实现上下层板密贴；最后修建剩余下穿段结构。为评估施工对轻轨隧道的影响，建立数值模型验算轻轨隧道底板弯矩，结果为安全；数值结果表明，底板沉降经历了先下沉、后抬升、再下沉的复杂过程，最终沉降曲线为单凹槽状。实测数据表明：沉降曲线可划分为振荡型、掉落型两种；曲线类型与测点位置有很大关系，震荡型曲线测点位于下穿段两端，掉落型曲线测点位于下穿段中部，分别受千斤顶、5#及6#导洞开挖控制；底板中线最终沉降为单凹槽状，最大沉降量与数值结果相符；变形缝差异沉降在千斤顶作用下一直处于可控状态。

不完整井在实际工程中应用广泛，且存在某些承压不完整井，其滤管顶、底均不在含水层层面，不能直接按常规的方法进行计算。可根据流线图中的水平流线对渗流场分层后进行简化计算，但分层位置的确定非常困难。因此，提出了简便可行的渗流场分层方法，包括假设上下部渗流场渗径相等的等渗径改进方法、假设上下部渗流场中井滤管长度与其含水层厚度比值相等的等比例改进方法和以等比例方法为基础的等效阻力法，并以此为基础对计算简便、精度较高的沙金煊公式进行完善，得到了适用于滤管顶、底均不在含水层层面的承压不完整井近似计算方法。经过算例验证和与有限元比较，该方法计算精度较高，实用性较强。

结合GIS（geographic information systems，地理信息系统）的空间数据处理能力，将潘家铮法从二维扩展到三维，提出一个基于栅格柱体单元的库岸滑坡滑速计算方法。建立了基于栅格柱体单元的三维滑坡滑速计算分析模型，给出了模型参数在GIS中的空间计算表达式。其次，通过对计算模型进行受力分析，结合牛顿运动定律，选取了沿滑体滑动方向和垂直方向来建立动力平衡方程，从而求解三维滑坡滑速。在GIS中开发出一个三维滑坡滑速计算扩展模块，并通过实例与潘家铮法所得结果进行对比。结果表明，该方法比潘家铮法计算的最大速度结果大15.2%，比启动加速度结果大32.8%，比达到最大速度的滑动时间短1 s。

桩端下伏溶洞顶板厚度是影响嵌岩桩竖向承载力的一个重要影响因素，在工程实践中，我国相关规范要求桩端平面以下的顶板厚度不小于3倍桩径。根据规范确定溶洞顶板厚度具有一定的经验性和笼统性，没有考虑溶洞顶板岩体质量，亦忽略了溶洞自身尺寸的影响。基于桩端岩体的冲切破坏模式，并结合广义Hoek-Brown准则和极限定理上限法建立了溶洞顶板临界厚径比（ ，即顶板厚度h/基桩直径d）的计算方法，给出了不同岩体基本质量级别的临界厚径比建议值。结果表明：极限端阻力、岩石坚硬程度和岩体质量均对临界厚径比产生影响，极限端阻力越小、岩体质量越好、岩石越坚硬，则溶洞顶板的临界厚径比越小；以《工程岩体分级标准》（GB/T50218－2014）中岩体基本质量分级为参考依据，I～IV级岩体的临界厚径比建议值分别是I级岩体 、II级岩体 、III级岩体 、IV级岩体 。

极限平衡方法计算锚固力作用下边坡稳定性时，通常将锚固力作为集中力处理，得到的滑面正应力和条间力分布曲线极不合理。除此之外，若按常规方法处理锚固力的作用，采用临界滑动场法分析锚固边坡稳定性时，搜索的滑动面会产生突变。为克服以上固有缺点，将半无限体受法向力作用的弹性力学解答近似作为锚固力在边坡体内产生的附加应力等效模型，并在此基础上建立了锚固力作用下的边坡临界滑动场计算方法，且通过算例分析表明了这种集中力近似等效处理方式的合理性，进一步丰富了锚固边坡稳定性分析方法。算例分析与工程应用表明，锚固力的作用会改变边坡临界滑动面位置，该方法能够搜索到锚固力作用下的合理临界滑动面，且滑面正应力和条间力分布曲线也较合理，相应的安全系数也较可靠。同时，该法能够全面评价边坡的整体和局部稳定性，可为实际边坡工程提供欠稳定区域的空间分布情况，使其得到有效的治理。

地下水是滑坡变形与失稳的重要影响因子，目前关于天然边坡水分运移过程的研究多借助传感器监测或数值模拟，但尚不能充分满足无扰、快捷的测量要求。针对该类问题，以天然边坡为研究对象，利用高密度电法结合时域反射技术，通过测量边坡土体含水率与电阻率，构建两者的对应关系，据此反演并分析边坡介质地下水的分布特征及迁移规律，验证高密度电法技术在探究天然边坡水分运移规律的有效性。现场试验结果显示：研究边坡土体含水率与电阻率呈明显对数关系；该边坡有效入渗或蒸发深度为2 m，大于该深度范围内，边坡土体含水率随深度增加逐渐减小；中雨强度下埋深为1.3～2.0 m范围内土体含水率激增的主要原因在于后缘裂缝的优势流补给；边坡高含水率（低电阻区）与低含水率（高电阻区）过渡区与后期滑坡滑动面位置高度吻合。

基于青藏铁路沿线P32和P33监测断面连续10年的含融化夹层路基的地温和变形场地实测数据，分析了该两处监测断面左路肩下多年冻土人为上限、季节冻结最大深度、融化夹层厚度及多年冻土人为上限附近地温的年变化过程；同时分析了P32和P33监测断面左右路肩的总沉降年变化过程、P32监测断面左路肩地温场对变形的影响及P33监测断面左右路肩地温场差异对左右路肩差异沉降的影响。结果表明：P32和P33监测断面左路肩下多年冻土人为上限逐年下降、季节冻结最大深度基本不变、融化夹层厚度逐年增厚及多年冻土人为上限附近地温逐年升高；观测期内，P32和P33监测断面左右路肩变形均以沉降为主，且P32监测断面左右路肩的总沉降变形量均小于P33监测断面；其中P32监测断面左路肩暖季沉降变形明显，冷季发生轻微的冻胀变形，且发生沉降和冻胀的时间略滞后于路基下部温度场的变化；P33监测断面左右路肩地温场的差异导致左右路肩存在差异沉降，且其差异沉降值随时间逐年变大。

为全面探究土岩组合岩体中抗拔桩的极限承载力，结合工程岩土参数及试验数据，运用Flac3D数值分析软件对其进行数值模拟分析，即可得到土岩组合岩体中抗拔桩的极限承载力。基于被动状态下的Kotter极限平衡方程式求解土层提供的抗拔力，根据岩石强度，基于Hoek-Brown破坏准则求解抗拔桩嵌岩端岩体的抗拉强度，从而可计算得到嵌岩端岩体的抗拔力；由静力平衡原理，叠加土层及嵌岩端岩体提供的抗拔力及破坏锥体重量，即可得到土岩组合岩体中嵌岩抗拔桩的极限承载力理论解析式。在嵌岩深度较小的情况下，该解析式的理论计算值与数值模拟分析值相接近，但随着嵌岩深度的增加，理论计算值会偏离数值计算值。故结合数值模拟试验值，对提出的极限承载力理论解析式作进一步的修正，得到修正后的极限承载力解析式能反映嵌岩端岩石风化程度、嵌岩深度、土层厚度、桩长对极限承载力的影响。运用修正后的解析式对该地质条件下不同抗拔桩的极限承载力计算表明：数值模拟结果与理论计算结果相吻合，说明所建立的抗拔桩极限承载力解析式的方法是可行的。同时，运用该方法可确定类似工程中嵌岩抗拔桩的极限承载力。

不同的桩-土剪切模式对桩-土剪切位移量的大小和桩侧摩阻力的发挥有重要影响。对抗压桩、顶拔桩和负摩擦桩等3种模式的6根钻孔灌注桩的剪切特性进行了现场加载试验与长期监测。通过对不同桩型的桩身应变、桩顶位移和桩周土沉降等试验资料的系统分析，建立了“侧阻-位移”曲线确定各模式桩-土剪切位移量和不同深度的极限侧阻力。根据不同桩-土剪切模式极限侧阻力的发挥特征，分析了桩-土相对运动方向、剪切边界约束条件、剪切位移量、加载速率和剪切历时等5个因素对极限侧阻力发挥的不同影响。讨论了利用现场顶拔桩试验结果对桩负摩阻力增大段的极限负摩阻力进行预测的可行性，供同行参考借鉴。

川藏铁路拉林段控制性工程巴玉隧道，施工过程中岩爆频发且呈现邻近区域岩爆位置不一现象。为此，通过微震监测系统原位观测了岩爆位置偏转过程中的岩体破裂演化规律，进一步探讨了岩爆位置发生偏转的原因。研究结果表明：岩爆位置偏转段两侧微震信息呈不均匀分布特征，微震活动（事件数、能量和视体积等参数）随岩爆偏转向岩爆发生位置集核偏转，两者具有良好的空间相关性。岩爆位置偏转段岩体损伤区及视应力的演化特征说明岩爆位置偏转与局部应力集中区变化相关，应力集中程度大的区域更易发生岩爆。相关研究成果对岩爆位置偏转预警及防治具有一定借鉴和参考意义。

基于地理信息系统（GIS）和数值模拟软件，提出了一套适用性强，运行流畅的三维边坡建模和模拟计算技术方案。此方案能充分利用复杂边坡的勘察历史资料，真正实现了从勘察数据处理－建立三维边坡模型－稳定性分析－结果反馈的边坡模拟全过程。为建立满足精度要求的模型，深入探讨了地表和土层的建模理论和方法，借助GRASS完成边坡的三维建模，并对已建模型进行详尽的误差分析，为量化建模质量提出了相应方法。用Python语言实现GIS和数值模拟软件的耦合，通过多个自主编写的程序和数据接口，完成不同平台间的数据传递和格式转换，避免了传输过程中数据的丢失或失真，还可深度挖掘已建模型的数据。借助杭口岭隧道工程实例，对此方案进行验证。分析结果表明，利用此方案建立的三维边坡模型，建模速度快，精度高，能准确地反映真实边坡的情况，并搜索到边坡潜在的最危险区域。

提出BP神经网络的分布式区域滑坡预测方法，算法设计在大数据分布式处理平台Spark下实现，通过构造包含均方误差和L2正则化的代价函数，提高运算实时性和算法泛化能力。统计影响滑坡评价因子的量化指标和定义监测剖面危险级别评价值，并进行评价因子特征选取，用于三峡库区忠县区域9个滑坡11年月监测海量数据挖掘，对研究区所有滑坡监测剖面每月进行危险级别评价，实现以月为周期的区域滑坡危险程度空间预测。试验表明，采用所述方法得到的拟合精度、准确度、效率均比梯度提升决策树、随机森林算法好，预测的滑坡危险级别准确，该方法可作为区域滑坡空间预测的一种新思路。

水力压裂技术广泛使用于页岩气开采工程中。为了分析压裂过程中多裂缝扩展形成复杂裂缝网的机制，尝试将态型近场动力学理论引入页岩水平井水力压裂过程的力学建模与数值仿真，在物质点间相互作用力模型中加入等效水压力项以实现在新生裂缝面上跟踪施加水压力，建立了水力压裂过程的近场动力学分析模型。通过模拟页岩储层的水力压裂过程，可得到复杂的裂缝扩展路径、裂缝网络的形成过程以及裂缝扩展受射孔间距及页岩天然裂缝和层理的影响。研究结果表明：射孔间距过小会造成起裂干扰，使中心射孔的裂缝扩展受到抑制；在压裂压力一定的情况下适当增大射孔间距，可以显著增强页岩压裂形成裂缝网的能力；压裂过程中水平层理面可能张开形成水平裂缝，且天然裂缝会诱导形成更复杂的垂直裂缝。模型和方法可为页岩水力压裂过程和机制研究及工程实践提供参考。

基于FLAC3D软件，以含软弱层的锚固岩体边坡为研究对象，通过修正的cable单元建模和改进剪应力提取方法，分别模拟分析了地震作用下含软弱层岩体边坡两锚固界面剪切作用及其演化。研究发现：砂浆-岩体界面上的剪应力远比锚杆-砂浆界面上的小，并均以锚杆上的中性点为界方向相反，分布很不均匀，且均在中性点附近发生突变；随边坡地震响应增强危岩中的锚固界面率先脱黏，基岩中的锚固界面紧随脱黏，且脱黏分别向锚头和锚根发展，直至拉拔段或锚固段全部脱黏锚固破坏。获得了地震波作用下边坡两锚固界面上剪应力及其分布和脱黏破坏过程，揭示了锚固界面上的剪切相互作用和锚固破坏机制，并得到了试验印证，为边坡锚固设计和施工提供了参考。