大型高尾矿库隐患治理是“十三五”期间遏制重特大事故的重要任务和重点工程之一。高应力下尾矿材料的细观结构和宏观力学特性演化、大型尾矿坝沉积规律及复杂条件下高尾矿坝的劣化理论是其中的关键问题。为揭示尾矿细观结构特性，对尾矿材料进行了光学显微镜观察、电镜扫描和CT扫描处理，建立了尾矿材料的三维重构模型；为揭示高压下尾矿材料的力学特性，对尾矿进行了高应力三轴试验、高应力渗透固结联合试验，研究了高应力下尾矿材料的强度、渗透和固结特性；为揭示尾矿的沉积特性，进行了尾矿材料的大型模型试验，得到了矿浆流动、淤积形态、尾矿料物理力学性质的空间分布以及夹层结构在空间内的分布规律；为揭示高尾矿坝的失稳机制，通过理论分析与数值模拟，对高压下尾矿的强度准则进行了讨论，高尾矿坝的稳定性分析建议使用强度折减法；为揭示尾矿的排渗设施淤堵机制与服役性能，进行了尾矿?土工布渗透试验以及微观观测，发现细粒尾矿更容易产生物理淤堵，且化学淤堵程度与溶液中离子的种类和浓度有关；为了推进尾矿坝风险管控技术，搭建了尾矿库重大风险评估监测系统平台，以便实现尾矿库基础数据的提取、实时风险预警等智能化功能。通过以上坝体劣化理论研究、溃坝灾害评估与预警，提升了我国大型高尾矿库筑坝与灾害防控水平。

法向应力对土与结构接触面循环单剪力学特性的影响研究具有重要意义。运用80 t大型三维接触面试验机，进行了不同法向应力下粗粒土与结构接触面单剪试验，分析了法向应力对接触面循环单剪力学特性的影响规律。试验结果表明，法向应力对接触面变形位移、滑动位移、切向应力、可逆性剪切体变和不可逆性剪切体变以及抗剪强度等接触面性能参量的数值有显著影响，但对参量间关系形式影响较小。法向应力越大，接触面变形位移在初始几个循环内越大，随循环剪切的进行，其减小速率越快，最终稳定值越小；切向刚度越大，切向应力稳定值也越大，初始几个循环内切向刚度系数越小；不可逆性剪切体变越大，可逆性剪切体变峰值和稳定值越小，相变位移和相变应力比越小；抗剪强度越大，其异向程度则先增大后减小。不同法向应力下，接触面循环抗剪强度与法向应力符合摩尔?库仑准则；切向应力比?切向位移关系、不可逆性剪切体变和剪切功密度关系均具有良好一致性，基本不受法向应力影响，可用双曲线模型描述，这将大大简化接触面本构描述。

基于可视化三轴压缩伺服控制试验系统和三维数字图像相关技术（简称3D-DIC技术），开展了不同流变方向系数α下的广义应力松弛试验，对广义应力松弛过程中砂岩表面的应变演化规律进行了探讨。研究结果表明：在流变过程中若砂岩有发生破坏的趋势，其表面的轴向、径向应变场的应变集中区域逐步显现，轴向应变集中区域由层状离散分布逐渐向即将出现裂纹的位置集中，而径向应变场中在裂纹即将形成区域附近应变增长速率最快。不同表面位置轴向、径向应变差异的演化速率与其整体演化速率的变化趋势呈现正相关，在即将形成裂纹附近的轴向、径向应变变化量发展经历减速、等速、加速3个阶段，而远离裂纹处的应变则可能先增大后减小，不同位置应变值的等时曲线演化则预示裂纹即将形成的区域。随着流变方向系数?由0.3、0、?3.0、∞、3.0依次变迁，流变过程中轴向、径向应变场整体演化的均值与方差均上升，说明其演化速率加快，应变集中现象更显著，且砂岩破坏前后的应变场演化差异增大。

交叉节理普遍存在于实际岩体中，但其对岩体力学破裂特性的影响至今未有全面深入的研究。制作了含3D打印交叉节理试样，采用声发射、CT扫描技术和数字图像相关技术（DIC）分析手段对单轴条件下含交叉节理试样的力学及破裂特性进行了研究。研究结果表明：（1）主节理倾角较小时裂隙的分形维数DB与单轴强度、弹性模量和泊松比之间的变化关系比主节理倾角较大时的一致性更好，前者破裂模式的规律性好于后者。（2）主节理倾角较大时压密和弹性阶段节理和块体之间在张拉作用下会产生聚集性微破裂声发射事件，主节理倾角较小时不会产生这种现象。DIC应变云图中裂隙扩展路径上的微应变先于宏观破裂，可通过计算压缩条件下微应变最大区域预测破裂模式。（3）主节理对起裂角和裂隙扩展模式起主要作用，且声发射和分形特性具有统一性。次节理平行最大主应力时，总损伤主要由主节理扩展造成，累积声发射事件数量AAEE和DB均最大。AAEE和DB越大，破裂模式越复杂。（4）三维破裂模式以节理端部扩展产生的2～3个张拉翼裂隙面为主，内部三维宏观裂隙面上存在由不完全张拉破裂引起的非连续结构。

沿海沿江城市地下开挖产生的工程渣土含泥量大、含水率高且较松散，主要运往渣土场进行堆填处置。由于产量巨大而处置场地有限，许多渣土场在运营过程中存在堆填速度快、缺乏排水设施、超高超库容堆填等问题，容易引发堆填体失稳事故。目前对于非饱和工程渣土堆体在快速堆填过程中的失稳机制认识尚不清晰，尤其是对这一过程中的高饱和度工程渣土强度增长规律缺乏足够的认识。以深圳红坳渣土场填料——花岗岩风化料（CDG）填土为研究对象，对不同初始饱和度土样进行三轴不排水不排气等向压缩和剪切试验，结果表明：非饱和CDG填土不排水抗剪强度随围压增大呈非线性增长，增长速率与试样初始饱和度密切相关；当土样压缩后的饱和度超过0.7，不排水强度随围压的增长速率显著降低。基于有效固结应力法的原理，结合Hilf孔压公式和修正剑桥模型，提出了一种工程渣土不排水抗剪强度估算方法，并通过与试验结果对比验证了该方法对初始饱和度高于0.6的CDG填土的适用性。利用该方法确定的不排水强度 与正应力 的关系可应用于高饱和度工程渣土快速堆填中的稳定性分析。

交叉裂隙作为裂隙网络的基本结构是控制场地尺度非饱和渗流特性的关键。随着流量的减小，交叉裂隙非饱和渗流将从连续、稳定的线状流转变为非连续、非稳定的滴状流，但目前对后者的研究还很少。在滴状流条件下，开展了交叉裂隙非饱和流动特性的理论和试验研究，发现了交叉裂隙滴状流存在毛细力驱动以及重力与毛细力联合驱动两种分流驱动模式。基于瞬态静力平衡方法，提出了一套液滴分流理论及计算模型，实现了动态分流过程的定量化描述与预测。结合模型预测与可视化试验，详细探究了两种模式下液滴的分流过程，阐明了动态分流行为受重力、毛细力、黏滞力共同控制的细观机制，并揭示了液滴长度、裂隙倾角、通道中液体的累积等因素对界面流速、分流体积比例的影响。研究成果为低流量、低饱和度条件下非饱和裂隙岩体渗流的预测与控制提供了理论及试验支撑。

堆石料的应力?应变曲线具有应变硬化和软化的特性，并且体积变形会伴随着剪缩和剪胀现象的发生；此外，应力水平、应力路径和密度等因素均会对堆石料的力学特性产生显著的影响。基于临界状态土力学的理论框架，建立了适用于堆石料的分数阶下加载面模型，该模型为典型的双屈服模型。在平均主应力p?剪应力q平面，剪切过程中堆石料当前应力状态点和参考应力点分别位于下加载面和参考屈服面上，两者之间的相对位置与孔隙比之差? 相关。相比于状态参量?，? 能够额外考虑应力路径对堆石料剪胀特性的影响。新模型的另一个显著特点在于能够考虑堆石料塑性流动方向与加载方向之间的差异性。在不引入塑性势函数的情况下，新模型对屈服函数进行Caputo微分并得到分数阶塑性流动法则，进而统一描述相关联和非相关联的塑性流动法则。新模型具有形式简单且参数较少等优点，所包含的7个材料参数均具有明确的物理意义，并且可以采用常规的室内试验结果进行标定。通过将模型计算结果与Tacheng堆石料在不同初始孔隙比和围压条件下的三轴压缩排水试验结果进行对比分析，结果表明新模型能够准确地描述堆石料的应力?应变曲线和体积变形特点；同时，新模型也能够合理地描述初始孔隙比对堆石料在e-lnp平面（e为孔隙比）临界状态线的影响。

将氯氧镁水泥（MOC）创新性地引入淤泥固化，通过无侧限抗压强度、固化体含水率、扫描电子显微镜（SEM）和X射线能谱测定（EDS）等试验，研究了初始含水率、MgO/MgCl2摩尔比、养护龄期和MgO活性等复杂因素对淤泥固化效能影响及其驱动机制。结果表明：初始含水率越高，试样强度越低，中等含水率时MOC固化淤泥中出现针杆状的3相或5相晶体；随MgO/MgCl2摩尔比升高，试样抗压强度随之增加，水化产物从无定形凝胶逐渐转化为3相、5相和Mg(OH)2晶体；随养护龄期延长，固化淤泥强度总体呈上升趋势，28 d前强度增长相对较快，28 d后强度趋于稳定；养护后期，MgO/MgCl2摩尔比较高的试样表面易出现泛霜现象；MgO活性提高使得MOC固化淤泥试样含有更多有效活性组分，无侧限抗压强度更高，但MgO活性高低对水化产物演变规律并无显著影响。研究成果可为绿色、低碳MOC基胶凝材料研发及其在淤泥固化等土体加固领域中应用提供理论支撑。

冻融作用可以改变土的结构，引起其力学特征的变化，从而影响到冻土工程的稳定性。由于冻融周期的设定不尽相同，因此大量的试验结果不能有效地进行对比分析。另外，对于冻融作用下土体力学特征的预测也成为研究的难点。在冻土遗传蠕变理论的基础上，利用球型模板压入仪作为测试方法，提出了一种冻融次数?物理时间比拟的方法。把冻融周期（次数）转化成物理时间（min），将不同冻融周期土的长期强度曲线簇映射到同一坐标中，利用该归一化曲线进行长期抗剪强度的预测，并选取了两种土样进行了测试，分别得到了相关的长期抗剪强度的预测方程。该方法对冻融作用下土体力学行为的成果对比研究以及长期强度的预测有重要的理论意义，对寒区工程稳定性预测分析也具有工程实践价值。

基于无黏性可压缩理想流体介质波动理论和Biot流体饱和多孔介质波动理论，考虑水下饱和土的流固耦合，借助Hankel函数积分变换法(HFITM)给出入射平面P1波在海底洞室周围散射问题的解析解。相比传统研究中的“大圆弧假定”，Hankel函数积分变换法可以较好地处理半空间表面边界条件。利用该解析解，计算分析了洞室表面透水条件、入射角度、入射频率、海水水深和饱和土的孔隙率等因素对水?土交界面处水平位移、竖向位移和洞室表面动水压力、环周总应力的影响。结果表明：洞室表面透水条件对水?土交界面处水平位移和竖向位移影响较小；随着斜入射角的增加，水?土交界面处竖向位移减小；随着入射频率的增加，水?土交界面处水平位移随之增加；海水水深为2.5倍SV波的波长时，水?土交界面处水平位移及洞室表面动水压力最大值最大；随着孔隙率的增加，水?土交界面处水平位移、竖向位移和洞室表面环周总应力减小，而洞室表面动水压力随之增加。

红黏土成团现象突出，导致石灰处治红黏土时难以拌和均匀；同时，红黏土呈弱酸性，石灰与黏土会发生互损作用，最终影响处治效果。首先，从水分子与黏土矿物的电荷分布特征角度，解释了红黏土呈弱酸性的原因；然后，通过酸、碱溶液浸泡石灰处治红黏土，模拟石灰?红黏土的互损行为；最后，提出应用偏高岭土抑制“石灰?红黏土”互损的方法。为此，选择两种团粒尺寸的红黏土，掺入不同比例的偏高岭土和石灰，并评价其无侧限抗压强度。结果表明，与仅用石灰处治的红黏土相比，掺入偏高岭土过多或过少都不利于其强度提高，掺量为5%效果最佳，验证了偏高岭土抑制“石灰?红黏土”互损行为的可行性。偏高岭土与石灰快速反应，不仅生成了以离子键结合的硅、铝酸胶结物，还生成了以共价键结合的网状胶结物，而且后者具有明显的抗酸侵蚀能力，从而解释了偏高岭土能够减少石灰?红黏土互损的内在原因。

目前通过对软土地基预加固处理来提高桩基水平承载力已被工程界认可，但如何在工程前期设计过程中估算软土地基预处理后桩基水平承载力提高值仍是技术难点。基于此，参考Bowles[1]的地基土水平抗力计算式，同时考虑成层软土地基预排水固结处理影响，通过数学推导，推求出根据原状软土室内土工试验抗剪强度指标及预加固处理时间，估算软土地基预处理后桩基水平承载力提高值的实用计算方法。考虑桩侧土弹塑性屈服影响，推导出成层软土中水平受荷桩弹塑性解析解及塑性区深度的计算式，给出了桩顶水平位移、桩身最大弯矩的无量纲计算式及相关计算源代码。依托于浙江省某水闸桩基工程案例，根据提出的计算方法对桩基水平承载力、桩顶水平位移及桩身最大弯矩等性状进行预估计算，并与地基预处理前、后现场试桩检测值进行验证对比，认为桩基水平承载力、桩顶水平位移及桩身最大弯矩等预估计算成果与工程现场试桩的检测值较接近，对类似工程设计具有较好的参考价值。

宽级配土作为堆积层滑坡的主要物源，其渗透性研究是开展降雨型堆积层滑坡机制研究的前提和基础。宽级配土的渗透是一个包括水分运移和土体细颗粒迁移的复杂过程，但在研究其渗透性时通常只考虑了水分的运移而忽略了细颗粒的迁移。为此，采用自制大型渗透仪对3组不同D15 /d85值（D15为粗粒组中小于该粒径的颗粒质量分数为15%的粒径；d85为细粒组中小于该粒径的颗粒质量分数为85%的粒径）的土样进行了饱和渗流试验，研究了宽级配土的水分运移特征和细颗粒迁移规律。研究结果表明：D15 /d85值对宽级配土的渗透系数和细颗粒迁移有重要影响，D15 /d85值越小，则土体渗透系数越小，细颗粒不易发生迁移；D15 /d85值越大，渗透系数越大，试验过程中渗透系数变化越剧烈，迁出细颗粒的量也更大。根据渗透系数的变化也可判定土体内部细颗粒的运动情况，据此提出了3种宽级配土颗粒迁移模式。该研究成果加深了对宽级配土渗透特性的认识，为完善降雨型堆积层滑坡机制研究提供了新思路。

非饱和土是地球表层土体一种常见的存在形式。与经典的饱和土Biot固结理论相比，非饱和土固结理论还亟待发展。基于Fredlund非饱和土双应力变量固结理论，放弃传统理论固结过程中骨架总应力不变的假设，推导出变荷载作用下非饱和土全耦合轴对称固结理论的控制方程组。通过Laplace-Hankel变换处理变量t和r，所得的控制方程被处理为常微分方程组；扩展的精细积分法将被进一步运用求解方程组以得到层状非饱和土地基在变换域内的固结解答，而最终的解答将通过Laplace-Hankel逆变换技术求得。将非饱和土地基退化为饱和土地基，与现有文献结果进行对比，验证所提方法结果的可靠性；最后，提供3个数值算例，以讨论加载时间T0、孔隙水关于净应力的体积变化系数 以及土体分层性对非饱和土固结特性的影响。

城市地下管线采用承插式接口连接，管线结构的连续性与否取决于接口刚度。然而，现有研究大都假定地下管线为连续结构，忽略了接口转角对管线变形的影响。开展土工离心模型试验和两阶段有限元参数分析，系统研究了盾构下穿不同接口刚度管线的变形特性。发现非连续管线的薄弱接口导致其整体抗弯刚度明显小于连续管线，且接口转动导致非连续管线具有更好的与土体协调变形能力。因此，盾构下穿施工引起的非连续管线沉降明显大于连续管线。非连续管的结构连续性假定将大大低估管线沉降。盾构下穿引起的接口转角与地层损失呈线性关系。基于系统的有限元参数分析，建立了区分相对刚性、柔性管线的无量纲组。针对相对刚性、柔性的非连续管线，分别建立了管节相对长度、管?土相对刚度与接口转角的预测方法，利用离心模型试验结果验证了预测方法的合理性。新方法能有效地预测盾构穿越地下管线的接口转角，研究成果可应用于城市地下管网改造工程。

采用半解析法研究了移动荷载作用下多层非饱和铁路地基的共振问题，分析了土体饱和度对非饱和地基振动的影响。采用分层非饱和多孔介质描述铁路地基；基于非饱和多孔介质的控制方程，采用多重傅里叶变换和土层界面位移、应力条件，推导了频率?波速域分层非饱和地基总体刚度矩阵。结合包括钢轨、轨垫、轨枕和道砟的轨道模型，建立耦合轨道?地基半解析模型。研究了移动简谐荷载作用在轨道上的多层非饱和地基的动力响应；考虑饱和度对土体剪切模量的影响，讨论了地基首层土和夹层土饱和度变化对地面振动的影响。采用频散曲线研究了地基共振模态。研究发现土体饱和度变化对振动位移的影响规律，与地基分层、荷载速度和频率有关。随饱和度减小，地基共振模态出现频率增大，地基关键速度增大。

宽级配砾质土是由砾石料和黏土料按一定比例混合而成，其具有压缩性低、抗剪强度大等特点，目前常作为土石坝心墙料或路基填料而得到广泛应用。由于宽级配砾质土渗透系数较小，在常规三轴固结排水剪（CD）试验中固结排水较慢，导致其试验周期很长。为了提高固结排水的效率，可采用一种在试样中心加圆柱形砂芯的快速三轴CD试验方法。基于快速三轴CD试验，通过变化不同砂芯类别、不同砂芯直径、不同掺砾量等各因素，全面研究了各因素对宽级配砾质土快速三轴CD试验的影响。研究结果表明，快速三轴CD试验方法能够有效加快试样排水固结，从而加快整个试验进程；砂芯类别、砂芯直径、掺砾量等因素均对试验固结排水速度和剪切过程应力应变产生不同程度影响；砂芯直径越小，其试验成果与无砂芯试样越接近。

工程建设中当距离抽水井r = rb处水位基本没有变化或不受抽水影响时，或当此处存在止水帷幕时，含水层系统视为侧向有限延伸，rb为有限半径。为此，构建更加符合工程实际的侧向有限延伸的典型弱透水层?承压水层系统中非完整井非稳定流计算模型，同时考虑井径和井储效应的影响，应用Laplace变换和分离变量法得到了水位降深在拉氏空间下的解析解，并应用拉氏数值逆变换Stehfest法得到真实空间下的水位降深。新建立的解析解可以进一步退化为诸多已有解，并进一步将其与已知解和有限元数值解进行对比，验证了所得解的正确性和可靠性。基于新建解重点分析了侧向边界和井的完整性对承压水层水位降深的影响。结果表明：含水层系统的侧向有限边界仅对抽水后期的水位降深影响明显，含水层系统侧向无限延伸情况下的水位降深要大于情形1（在r = rb处为定水头边界）且明显小于情形2（在r = rb处为不透水边界）下的水位降深，rb越小，两者之间的误差越大；抽水井的完整性对整个抽水期间不同情形下的水位降深均有明显的影响，承压含水层顶板处的水位降深随着抽水井滤管的长度和埋深的增加而减小。

随钻跟管桩施工不能完全清除桩底岩土沉渣，从而影响桩基端承力。为揭示桩底沉渣对随钻跟管桩承载力的影响机制，开展了考虑桩底沉渣影响的随钻跟管桩竖向承载特性模型试验研究。试验结果表明：在密砂地层中，具有桩端水泥土扩大头的随钻跟管桩，其桩顶荷载?沉降曲线为缓降型，而模拟试验的其他管桩均为陡降型；桩底沉渣降低随钻跟管桩的极限承载力在22%以内，且其桩顶荷载主要由桩侧摩阻力承担，承担占比超过90%；与存在一定厚度沉渣的钻孔灌注桩相比，随钻跟管桩的桩底沉渣对降低承载力的影响相对较小；靠近桩端的轴力随着沉渣厚度的增加而减小，沉渣越厚，减少的幅度越明显；桩端水泥土扩大头施工可提高随钻跟管桩约37%的承载力，且桩端阻比均小于15%。现场原位测试（桩长为15.5 m，长径比为15.50）和室内模型试验（桩长为1 m，长径比为15.87）结果均表明：存在桩底沉渣时，随钻跟管桩是以发挥侧摩阻力为主的端承摩擦型桩。研究成果有助于进一步加深对随钻跟管桩承载性状的认识。

温度梯度、含水率梯度和压力梯度对非饱和土体中的水热迁移有重要影响。用非饱和黄土填筑室内模型，通入0.10 MPa高温水蒸气，研究蒸汽压梯度、温度梯度和含水率梯度耦合作用下的水热迁移规律。结果表明：在水蒸气迁移范围之内，土体的升温速率、增湿速率和增湿程度较大，以蒸汽传热和压力梯度驱动的水蒸气增湿为主。而在水蒸气未到位置，以温度梯度引起的热传导和含水率梯度与温度梯度耦合驱动的水热迁移为主；水蒸气迁移时，受土颗粒阻碍和蒸汽压消散的影响，土体升温速率、增湿速率和温度传导速率均随径向距离的增大而减小，且增湿速率小于升温速率；水蒸气增湿土体的最大含水率接近最优含水率，增湿效果较好，可有效提高一定范围内土体的压实性能；基于模型试验边界条件，确定含湿毛细多孔介质中二维热湿迁移方程的一组代数显式特解，并以20 cm测点为例，将温度和含水率实测值和计算值进行对比分析。研究结果可为非饱和黄土水?汽?热耦合传输规律和水蒸气增湿新技术研究提供理论支持。

重载铁路基床表层级配碎石渗透系数的快速、精确测定，对工程施工质量及工期的保证具有重要意义。基于达西定律研制适用于级配碎石的大型常水头渗透试验装置，并利用该装置考虑边壁效应进行一系列渗透试验，探究试样孔隙比e、不均匀系数Cu及特征粒径参数对其渗透系数K20的影响程度及规律。结合试验数据，以太沙基模型为基础，建立同时考虑孔隙及级配特征参数的渗透系数改进计算模型，并通过试验验证其精确性。结果表明：级配碎石的K20与其e2、平均粒径平方d502呈线性关系，与Cu呈负指数关系；d50对K20的影响权重最大，Cu次之，e最小；提出的计算模型具有较高准确性，可为新建重载铁路基床表层施工中填料渗透系数的快速测定及运营路基基床的渗透稳定性评估提供参考。

在求解非饱和态土中水分入渗问题时，水力函数是体积含水率或者吸力的函数，致使其控制方程呈现出强非线性的特征，进而使得其求解变得十分困难。基于水分在土体介质中流动耗时取极值路径的选择这一假定，引入时间泛函，基于变分法原理将水平入渗问题转化为泛函极值问题。通过求解Euler–Lagrange方程，结合边界条件，得到非线性瞬态水平入渗问题的显式解析解。结合Brooks-Corey型水力函数，显式地求解出该类型非饱和态土的体积含水率发展分布规律。通过计算4种不同类型土体的水平入渗规律，将求解结果与已有结果以及数值结果进行对比，验证了该方法的有效性。结果表明：体积含水率分布与位置距离和湿润峰距离比值呈幂函数关系，指数取决于土?水特征曲线的形状参数；初始条件与边界条件会对体积含水率分布造成不同程度的影响。

针对三峡库区“阶跃式”滑坡的变形特征，提出了一种新的滑坡位移预测方法。以白水河滑坡ZG118和XD-01监测点位移数据为例，采用基于软筛分停止准则的经验模态分解(SSSC-EMD)将累计位移?时间曲线和影响因子时间序列自适应地分解为多个固有模态函数(IMF)，并采用K均值(K-Means)聚类法对其进行聚类累加，得到有物理含义的位移分量（趋势性位移、周期性位移以及随机性位移）和影响因子分量（高频影响因子和低频影响因子）。使用最小二乘法对趋势性位移进行拟合预测；采用果蝇优化?最小二乘支持向量机(FOA-LSSVM)模型对周期性位移和随机性位移进行预测。将各位移分量预测值进行叠加处理，实现滑坡累计位移的预测。研究结果表明，所提出的(SSSC-EMD)-K-Means-(FOA-LSSVM)模型能够预测“阶跃式”滑坡的位移变化规律，且预测精度高于传统的支持向量机回归(SVR)、最小二乘支持向量机(LSSVM)模型；并通过改变训练集长度，进行单因素分析，发现其与预测精度之间呈正相关关系。

以岩土材料力学参数空间变异性的“点估计?有限元”分析方法为基础，结合节理分析时自身存在几何模型、网格划分等特性，扩展了该方法在节理分布空间变异性分析方面的适用性，明确了具体的研究步骤与方法。以某抽水蓄能水电站为例，通过分析节理空间变异性对围岩变形与塑性区的影响，验证了扩展后该方法的准确性和合理性。对工程案例开挖揭露的1 400余条节理进行概率统计，建立了节理空间变异性的有限元分析模型；采用扩展后的概率分析方法，研究了节理分布对地下洞室群围岩开挖稳定性的影响。研究结果表明：（1）对比概率分析得到的围岩变形概率分布与现场监测结果，发现剔除变形异常点后监测变形量值大部分位于得到的位移概率分布范围内，说明节理的空间变异性是导致监测变形波动的主要影响因素；（2）围岩变形概率分布的标准差能有效识别出围岩开挖变形受节理空间变异性的影响程度，对于所给出的案例依次为：机窝＞边墙＞顶拱；（3）围岩塑性区的概率分区能合理判断地下洞室群开挖时受节理影响较大的区域和范围，为工程施工的支护设计提供依据。

当盾构隧道平行侧穿建筑物时，大多关注建筑物的横向沉降规律，对其纵向沉降关注较少。为此，针对盾构隧道平行侧穿建筑物引发的空间变形开展研究。首先，对天津地铁6号线平行侧穿四座结构形式相近的砖混建筑的实测数据进行分析，得到建筑物基本变形模式；基于工程实测并考虑土体的小应变硬化特性建立三维有限元数值分析模型，研究了盾构侧穿引发的建筑物纵向挠曲、土体变形与应力变化规律，并分析了不同建筑平面长宽比的影响。结果表明，盾构隧道平行侧穿将诱发平面长宽比较大的建筑出现“下凹式”挠曲变形，纵墙中部沉降最大可为其角点沉降的2倍，平行侧穿并不能简化为平面应变问题进行分析。建筑物修建和盾构开挖将导致隧道上方土体经历较为复杂的应力变化过程，并可划分为6个阶段。沿建筑纵向基础中部的土体与边缘土体相比，其首先经历更大的压缩变形（建筑施工导致），在盾构穿越后又产生了更大的卸荷变形。当建筑平面长宽比小于2时，盾构开挖导致的纵向挠曲变形将显著减小。

以成昆铁路复线加筋土挡墙为工程依托进行现场原位试验，对比分析工后9个月内新型整体式面板、内嵌返包结构的模块式面板和模块式面板3类加筋土挡墙的结构特性。在此期间发生两次地震，监测了震后挡墙变形及土压力变化。结果表明：整体式面板加筋土挡墙整体稳定性最好，具有良好的抗震性能，格栅应变变化率、墙体压缩量和墙面水平位移均最小。整体式面板和模块式面板加筋土挡墙墙背土压力沿墙高近似呈M型，内嵌返包结构的模块式面板加筋土挡墙墙背土压力沿墙高近似呈倒S型。整体式面板加筋土挡墙的侧向土压力系数沿墙高变化较小，小于美国联邦公路局（FHWA）加筋设计指南计算值；内嵌返包结构的模块式面板加筋土挡墙侧向土压力系数沿墙高呈增大的趋势，挡墙中、下部小于FHWA计算值，上部接近或大于静止土压力系数；模块式面板加筋土挡墙侧向土压力系数大部分处在FHWA计算值与静止土压力系数之间。3类挡墙土工格栅应变沿墙高均呈非线性变化。挡墙墙体压缩量随着时间的增加逐渐增加，在工后前50 d时间内增加速率较快，随后增加速率减缓，进入雨季之后增加速率再次增大。发生地震后，3类加筋土挡墙均出现不同程度的墙背土压力减小、格栅应变增大、墙体压缩量增加和墙面板外移的现象。

数值模拟是露天矿山边坡稳定性分析的常用方法，边坡模型构建是数值模拟的基础，地形测量的准确性对于边坡稳定性分析精度具有十分重要的影响。近年来，低空低速小型无人机（UAV）被广泛应用于地质调查中，为复杂条件下地形快速测量提供了一种可行的方法。针对露天矿山边坡坡形更新快、地形复杂、人工测量困难的情况，采用搭载单镜头小型无人机，对露天矿山边坡进行低空摄影测量，并通过飞行试验，对比研究了无人机搭载镜头旋转角度对边坡测量精度的影响，最终得到高精度边坡数字高程模型（DEM），并通过Surfer、Rhino等软件转化为三维FLAC模型，实现从无人机测量到三维数值模型的快速建立。得到以下可供参考的研究结论：在目标边坡坡度小于45°时，无人机搭载镜头旋转角度不同对边坡测量精度有较大影响，当镜头旋转角度介于70°～90°时，获得的测量精度较高；与人工测量相比，无人机摄影测量获得的边坡坡形更符合真实边坡实际坡型；无人机航拍照片经过PhotoScan、Surfer、Rhino等软件的转化，能够快速建立三维数值模型，并获得更好的数值模拟计算结果，满足边坡稳定性分析需要。

为有效模拟裂缝性页岩储层中水力裂缝随机扩展过程，基于单元节点的拓扑数据结构，利用网格节点分裂方式，建立了一种基于有限元网格嵌入零厚度内聚力单元的水力裂缝随机扩展新方法。利用KGD模型解析解和2种室内试验，验证了新方法的准确性和有效性。同时，通过数值算例研究了水平地应力差和储层非均质性对水力裂缝随机扩展过程的影响。研究表明：（1）该方法弥补了ABAQUS平台内置的内聚力单元无法有效模拟水力裂缝随机扩展的不足；（2）在较高水平地应力差下页岩储层非均质性越强，与水力裂缝相交的高角度天然裂缝越容易开启。所建方法能准确地描述复杂水力裂缝的随机扩展行为，可为裂缝性页岩储层的数值模拟提供新手段。

常见的地质灾害如滑坡、泥石流、岩崩等通常都涉及不同形状的颗粒物质运动，这些形状不同的颗粒又多具有不同的尺寸和含量。基于典型的颗粒柱坍塌试验，首先根据试验方法确定了离散元模拟所需的各项参数，然后采用随机多面体方法生成了可控制长细比的大颗粒，利用离散元法就不同大颗粒含量下形态变化对二元颗粒柱坍塌特性的影响开展研究，研究结果表明：（1）利用离散元法可以较好地重现室内试验中小球和多面体组成的二元颗粒系统的颗粒柱坍塌过程；（2）在不同长细比的不规则大颗粒和小球组成的二元颗粒柱系统中，当大颗粒含量高于临界含量值20%时，二元颗粒柱坍塌持续的时间随非球形大颗粒长细比的增加而增加；（3）在不同长细比的不规则大颗粒和小球组成的二元颗粒柱中，当大颗粒含量高于临界含量值20%时，在相同百分比的大颗粒含量下，大颗粒长细比的增加会提高大颗粒平均配位数以及降低颗粒的运动能力，大颗粒间形成更强的互锁作用，降低了颗粒柱的整体流动性，使其最终堆积高度更高、最大跑出距离更短以及更小的归一化动能峰值。（4）在不同长细比的不规则大颗粒和小球组成的二元颗粒柱中，小颗粒可以较为明显降低大颗粒间摩擦及互锁作用，增加流动性，降低大骨料形态对坍塌过程的影响。

为了研究多年冻土表层的水热分布情况，在非饱和冻土的能量守恒方程和水分迁移的质量控制方程的基础上考虑冰水相变和水汽相变过程，并考虑水汽运移传热及温度势对水汽迁移的影响，建立了非饱和冻土的水?热?汽耦合模型。采用光滑粒子流体动力学（smoothed particle hydrodynamics，简称SPH）方法可方便地计算它们的演化过程。为此，在计算中先求解能量守恒方程的含冰量及气态水含量，再对未冻水含量和温度场进行求解，从而实现了温度场与水汽场的耦合。在此基础上，模拟计算了第1类热边界条件下半无限空间介质内非稳态温度场、体积含水率及水汽通量的分布情况，并将计算结果与未考虑耦合的解析解进行比较，结果显示水汽耦合的作用不容忽略。最后，针对处于季节性周期温度边界下路基的水热场的分布情况进行计算。研究表明，相比于水?热耦合模型，所建立的水?热?汽耦合模型得到的计算结果更为接近实际监测结果，可很好地揭示非饱和冻土中的水热汽迁移特征及其相变过程。