为研究具有负泊松比效应（negative Poisson’s ratio，NPR）的恒阻大变形锚索在地震作用下的吸能特性与动力响应规律，利用振动台开展含断层隧道使用普通锚索支护以及使用NPR锚索支护的振动台物理模型试验，对比分析地震作用下处于不同围岩条件中以及使用不同锚索进行支护时隧道结构的动力响应规律。研究结果表明：在同一地震波作用下，使用NPR锚索进行支护的隧道结构破坏程度与使用普通锚索进行支护的隧道结构相比有明显的改善，表明NPR锚索在地震作用下具有良好的吸能特性，能够为隧道结构提供有效支护；通过对比分析普通锚索与NPR锚索的锚索轴力时程曲线可知，在地震波作用下，普通锚索会发生破断现象，轴力突降至零，而NPR锚索的高恒阻力能够持续提供有效支护；分析隧道结构各测点的加速度、应变、土压力等动力响应可知，隧道断层区域在地震作用下的动力响应更为明显，隧道结构在地震作用下呈现出拱肩→拱腰→拱顶的破坏规律，这一研究成果为类似隧道工程的支护对策提供了理论参考。

受降雨作用，球状风化花岗岩类土质边坡的土-岩差异风化界面极易演化为优势渗流通道而发生渗流潜蚀，进而加速该类边坡的变形失稳，然而当前有关其渗流潜蚀作用特征、细颗粒迁移规律等的研究仍鲜见开展。基于多孔介质非饱和渗流理论，综合考虑细颗粒运移、潜蚀启动响应与非饱和渗流的耦合关系，提出一种可准确描述土-岩界面渗流潜蚀过程的数值计算框架。采用有限元方法，构建优势流作用下非饱和花岗岩残积土的渗流潜蚀模型，并以均质土柱的渗流潜蚀过程为参考，系统研究3种典型土-岩界面埋藏状态下的优势流潜蚀特性。结果表明：球状风化花岗岩类土质边坡的土-岩界面与基质渗透性存在高度差异性，湿润锋形成向下凹陷的渗透漏斗，且随着降雨的持续，湿润锋的凹陷程度愈发明显；细颗粒流失程度与土-岩界面的埋藏状态相关，其中下填土体工况的优势流潜蚀最为显著，其界面处甚至出现超孔隙水压力，最不利于该类边坡的稳定性。研究成果可为降雨条件下球状风化花岗岩类土质边坡稳定性的准确评价提供科学依据。

为了更准确地分析桩顶在承受水平动荷载H(t)和扭转振动T(t)组合荷载作用下的受力变形特性，基于Pasternak模型，考虑多向荷载间的相互作用，建立层状地基中H(t)-T(t)受荷桩简化分析模型。利用有限杆单元法推导出桩身内力位移数值解，并与已有相关理论解、模型试验和有限元模拟结果进行对比验证。参数分析表明：（1）与Winkler模型计算结果相比，考虑土体剪切效应后，桩顶水平位移和桩身最大弯矩分别减小了11.8%和10.5%。（2）增大外荷载无量纲频率可减小桩身水平位移和弯矩，但同时会减小H(t)-T(t)耦合刚度。（3）在层状地基中，表层土对桩身内力位移的影响最大，且存在临界影响深度，表层硬土的临界影响厚度是表层软土临界影响厚度的3.5～6.5倍。（4）采用杆单元模型可减少单元划分数量和计算时间，有效提高计算效率。

为深入探索断层黏滑失稳的阶段特性以及声-电多源联合前兆特征，对含不同预制倾角断层的正长花岗岩开展了侧压分别为5、15、25 MPa的双向剪切摩擦试验，观察和分析了断层黏滑失稳过程的力学行为以及声发射、电荷感应信号的时频响应特征。结果表明：（1）侧压和断层倾角改变了断层的黏滑特征，进而影响了感应电荷量和声发射事件震级。侧压越大，间黏滑和黏滑失稳阶段的电荷量占比以及电荷累积速度越大。相较于56º断层，45º断层黏滑失稳时基于振幅的统计指标b值更小，大尺度破裂占比更大。（2）声-电信号幅值的分形特征可以反映断层面黏滑错动释放能量的差异。声-电信号频域的阶段性响应特点可有效辨识实验室尺度下断层黏滑失稳的演化过程，可以将声-电信号主频幅值急剧提高作为断层黏滑失稳前兆特征之一。临近断层黏滑失稳时声发射和电荷感应参量（能量、分形维数、频谱幅值）均呈现明显的“突增”特性。（3）声-电信号能量突增与黏滑应力降之间具有较好的对应性，但并不完全同步。断层黏滑失稳时，声发射会先于电荷感应产生高幅值信号。仅采用单一声发射或电荷感应信号来预警有可能造成信息漏报、误报，需结合声发射和电荷感应信号特征以实现声-电二元信号的相互补充与验证。

为探明微震信号能量在采空区“三带”结构中的衰减规律，拟开展采空区覆岩相似模型试验，采集人工激发微震波经由采空区结构传播的微震信号，通过变分模态分解（variational mode decomposition，VMD）处理微震信号，获取各频率下模态分量。针对采空区微震信号在VMD下各模态分量中心频率与能量之间的关系展开分析。根据中心频率法确定微震信号最佳模态数量，并计算微震信号欠分解状态、最佳分解状态、过分解状态下各分量能量；对各震源下信号最佳分解状态时各模态分量能量与中心频率分布关系进行拟合，分析在“三带”结构中，微震信号不同传播状态下各结构层对信号能量影响作用。研究结果表明：（1）在VMD过程中，人工激发震动信号有效模态数量在6～11范围内，微震信号能量随模态数量变化明显。（2）采用幂函数可实现对微震信号模态能量与频率关系的拟合，且拟合状态良好（决定系数大于0.9），其中低频模态分量包含能量占信号总能量近50%；采用高斯函数可以拟合震源各分量能量在频域上的分布表现，拟合状态较好，且表现出高斯单峰特征。（3）微震信号穿越采空区“三带”结构，微震信号能量随震源位置与传感器距离增加而减小，同时信号能量随震源位置到达传感器穿越岩层数量增加而减小，信号能量在经由垮落带时，能量变化明显，相较于裂隙带和弯曲下沉带，垮落带对信号能量衰减作用明显。

锚托板支护结构因为其安全稳定、经济和安装简便等优点在填方加固边坡中得到推广应用。然而目前在锚托板支护结构设计过程中，锚托板尺寸仍没有成熟的计算方法，设计过于保守。利用对数螺旋线滑移面结合极限平衡法计算得到非饱和土的主动土压力。根据土体内部平衡求出每层锚托板所需承担的拉力，得到了锚托板尺寸的计算方法，根据滑移面的位置得到了锚托板自由段拉杆的最小设计长度。将主动土压力计算结果与现有的经典文献进行了对比分析，发现主动土压力系数变化趋势与现有结论一致。分析了降雨工况和地震工况对锚托板尺寸的影响，结果表明，在降雨强度较小时对锚托板锚固尺寸的影响可以忽略，但在同时考虑强降雨和地震时，会显著增加所需的锚托板锚固尺寸，最大增大幅度可达136.71%。分析结果可为锚托板锚固尺寸的设计提供参考。

采用极限上限分析定理构建了路堤失稳的三维旋转破坏机构，继而根据Bishop非饱和土抗剪强度理论，考虑非饱路堤内部基质吸力的空间分布及其随地下水位的变化，建立了路堤破坏土体外力功率与内部能量耗散功率的能量守恒方程，并利用遗传算法编写了非饱和路堤最小上限解的高效搜索算法。通过将非饱和土路堤基底破坏模式退化为边坡破坏模式，并与现有非饱和土边坡稳定性计算结果对比，验证了所提上限解的正确性和遗传搜寻算法的准确性。进一步地，对路堤三维稳定性的关键影响因素展开了系统分析，研究了路堤填土孔径分布、基质吸力、进气值倒数、路堤倾角以及有效内摩擦角等因素对路堤三维稳定性的影响规律。研究表明，非饱和土质路堤的稳定性不仅取决于路堤填土性质，而且依赖于影响土体吸力大小与分布的填土孔径参数和进气值等因素。地下水位升降引起的基质吸力变化对路堤稳定性存在显著影响。研究结果为路堤稳定性精细化分析提供了重要的理论依据。

假设瞬态荷载沿隧道方向按指数衰减，建立非饱和介质中圆柱形衬砌隧道受到内源爆炸荷载作用时衬砌隧道及周围非饱和介质计算模型。基于多孔介质混合物理论和连续介质力学理论，推导爆炸荷载作用下衬砌结构及周围非饱和土的控制方程，运用Laplace及Fourier变换求得隧道衬砌及周围土体在频域-波域内的解析解答，研究内源爆炸作用下饱和度对衬砌结构动力响应的影响。计算结果表明：饱和度对衬砌内表面径向位移影响较大，在轴向上，随着饱和度由0.15增至0.98，爆源处的径向位移值升高约28%；衬砌内表面径向应力随时间变化受饱和度的影响不明显；衬砌内表面环向应力幅值随饱和度的增大而增大，且达到幅值所用的时间逐渐延迟。在轴向上，随着饱和度由0.15增至0.98，爆源处环向应力升高约29%。

流态性滑坡的大变形运动过程受到诸多因素的影响，其中一个重要因素便是岩土体的强度参数。由于现场测试和室内试验的局限性，强度参数往往呈现出明显的空间变异性，同时在不同方向上具有不同的波动范围，即随机场波动范围具有各向异性。针对黏聚力随机场波动范围的各向异性对流态性滑坡滑动距离的影响问题，引入基于乔列斯基分解的中心点法实现了各向异性随机场的离散，采用摩尔-库仑破坏准则和非牛顿流体模型相结合的光滑粒子流体动力学（smoothed particle hydrodynamtcs, SPH）分析方法模拟了流态性滑坡的滑动距离，在蒙特卡洛模拟框架上提出了流态性滑坡的随机分析方法。然后，通过羊宝地滑坡和水平地层模型的模拟，验证了流态性滑坡确定性分析方法和随机场离散方法的适用性。最后，根据汶川地震中王家岩滑坡的地形资料构建了概念性的滑坡分析算例，讨论了黏聚力随机场中各向异性波动范围影响下滑坡运动过程的变化，分析了滑动距离的概率分布规律。结果表明：竖直波动范围的增加会增大滑坡滑动距离的分布范围，使其表现出更强的离散性；在黏聚力参数符合对数正态分布的前提下，滑动距离的分布也符合对数正态分布，说明流态性滑坡滑动距离的概率分布与强度参数的不确定性具有一定的联系。

拓展传统极限平衡方法在复杂条件下边坡稳定性分析中的适用性一直是岩土工程界的研究热点。众所周知，边坡岩土体剪切失效一般服从于非线性强度准则，而传统极限平衡方法仅适用于线性强度准则下边坡稳定性分析。此外，预应力锚索作为一种主动防护措施，可极大地即时提高边坡稳定性而被广泛应用于中、大型边坡工程加固。然而，在应力松弛、蠕变效应以及腐蚀环境影响下预应力锚索加固边坡稳定性时效特征往往被忽视。由此，基于Morgenstern-Price（M-P）法，嵌入非线性Mohr-Coulomb（M-C）强度准则，并融入应力松弛与蠕变效应下锚索预应力损失模型以及腐蚀环境下锚索性能劣化模型，从而，建立非线性强度准则下预应力锚索加固边坡安全系数时程表达式。进一步，为了解决非线性强度准则下隐式嵌套公式难以求解的问题，归零化条间法向力初始值，并应用迭代循环求解策略，在引入循环终止允许误差判别条件后，可实现边坡稳定性理论真实结果的逼近求解。经过后续算例对比分析，验证了所提方法的可行性和合理性，与此同时，还研究了不同预应力损失速率和不同锈蚀环境下预应力锚索加固边坡稳定性的时效变化规律。

对含不同强度粗颗粒的滑带土进行了可视化中型直剪试验，通过竖直置入铝丝观察剪切后铝丝变形情况获取剪切破坏带的空间信息，并建立了剪切带空间曲面方程。通过粒子图像测速（particle image velocimetry，PIV）技术对可视面处二维剪切带信息进行提取，并与剪切试验得到的空间曲面方程的边界外推值进行对比，证明了空间曲面方程与PIV技术能够描述剪切带特征。然后通过PIV技术对不同剪切总位移与某一剪切位移下短时间内相对位移规律进行了分析。最后分析了预制损伤与粗颗粒强度对剪切带特征的影响。研究结果表明，剪切带形状为两端薄中间厚的带状，可采用高斯曲面方程对剪切带轮廓进行拟合。剪切带在发展过程中经历压密、自由损伤、损伤沿固定路径发展与剪切带贯通4个阶段，损伤使得剪切带发展各阶段前移，粗颗粒强度会对局部剪切带变形产生较大影响而对整体剪切带的变形影响较小。该研究成果对于揭示滑坡剪切带形成与演变规律从而设计合理的边坡防治方案具有重要意义。

砂土中夹层的性状会影响饱和砂土孔压发展，从而影响砂土层变形。为研究夹层位置、厚度和种类等不同夹层状态下砂土液化过程中的孔压变化规律，设计了冲击荷载作用下层状砂液化试验，建立了含夹层饱和砂土理论模型，并将试验结果与理论分析进行对比。结果表明：含夹层饱和砂土的孔压发展呈现3个阶段，即快速上升、快速消散、缓慢消散阶段。高渗透性夹层高度越高，其下方土层孔压快速消散时长越短，越快趋于稳定值，但消散总时长无明显影响；低渗透性夹层高度或厚度的增大，均会使夹层上方孔压快速消散阶段速率加快，孔压消散稳定阶段延长，孔压消散总时长随之线性增长；同时，孔隙水会在低渗透性夹层下方形成水膜，夹层高度或厚度的增加均会使水膜持续时间增长，但水膜形态主要受夹层厚度影响。试验结果与理论分析较为一致，说明了试验的可靠性。

为深入研究干湿-冻融循环对原状黄土力学特性的影响及细观损伤演化规律，通过不同次数干湿-冻融循环条件下的固结排水三轴剪切试验（consolidation drainage triaxial shear test, CD）和核磁共振试验，从宏-细观角度分析其应力-应变曲线及强度指标变化规律和细观孔隙的损伤变化规律。在此基础之上，假定黄土微元强度分布服从复合函数，建立了黄土的损伤统计本构模型，并验证其适用性。研究结果表明：土体的应力-应变曲线表现为应变软化，并且随着循环次数的增大软化程度逐渐减弱。偏应力峰值随循环次数减小并逐渐趋于稳定，在循环2次时其衰减程度最大，不同围压下分别衰减了17.6%、23.2%、24.5%和18.1%。土体内胶结块在循环作用下发生破损，使得内部孔隙面积逐渐增大，主要为小孔隙向中大孔隙的转变，随着循环次数的增大，内部结构逐渐趋于稳定。

场地基本周期是场地效应模型的重要变量，是场地分类的重要指标。已有基本周期的计算假定基岩地震波为垂直入射，常用方法包括利用垂直入射情形下场地内SH波竖向一维传播模型建立的模拟与经验方法，其中逐层单自由度法这一经验方法因计算简易、结果与模拟方法相当而被日本抗震设计规范所采纳。然而垂直入射假定大多与实际不符，为此给出可考虑基岩地震波斜入射的场地基本周期计算方法。首先利用Snell定律将基岩斜入射地震波作用下场地内的SH波二维传播表征为一维竖向视传播建立了斜入射情形场地内SH波传播模型。其次利用这一模型给出了斜入射情形基本周期计算的模拟和修正逐层单自由度法并应用于KiK-net台站场地，结果表明基本周期随入射角的增加而变小，就统计所得平均入射角57°而言，基本周期将缩减11%；修正逐层单自由度法可准确考虑斜入射的影响且计算简易度和原方法相当。最后，讨论了基岩面选取对基本周期计算的影响，结果表明，若依据土体首次达到给定剪切波速以确定基岩面，当剪切波速大于700 m/s时，不同基岩面所对应的基本周期具有高度线性相关性。

弄清温度作用下海相软土的力学变化对海底管道的建设与长期运行有着重要作用。针对广东汕尾海相软土，进行了一系列的基本物理力学性质试验，并通过X射线衍射仪、X射线荧光光谱仪和扫描电镜仪对其矿物成分、元素组成和微观结构进行了分析。通过温控三轴仪对其开展了等向固结与偏压固结条件下的固结不排水室内试验，试验结果表明：当轴向应变较小时，环境温度和固结应力比K_c对饱和广东软土的不排水剪切特性影响表现为环境温度越高或固结应力比越小，其割线模量越大，且其割线模量与固结应力比呈反比，与环境温度呈幂函数关系。轴向应变增大时，其峰值强度受环境温度和固结应力比共同影响，而其峰值孔压受环境温度影响较小。试样初始干密度会对等向固结土（K_c=1.00）的不排水剪切特征产生影响：试样初始干密度越大，峰值孔压越小，峰值强度越大，且温度升高时其峰值强度的下降幅度越大。固结应力比会影响不同环境温度下饱和广东软土的峰值强度与有效应力路径：当K_c =1.00、0.67时，环境温度上升对饱和广东软土起软化作用；当Kc =0.50、0.40时，环境温度上升则对其起硬化作用。

强夯法是一种加固某些土体的施工操作较为简单的方法，合理分析其夯击作用在土体中产生的动应力是工程实践中的关键环节之一。针对成都平原地区的典型土-石混合填土，通过现场夯击试验，测得了在4 000 kN·m的夯击能作用下土体中不同深度处的竖向夯击动应力，鉴于既有算法与试验值的差异，从两个方面确定了适于土-石混合填土夯击动应力的简单计算方法。一方面针对较为复杂的既有理论公式问题，采用敏感性分析与回归分析法，建立了以夯锤质量、夯锤落距、夯锤半径、土体密度、土体动力剪切模量、土体泊松比、土体阻尼比、夯锤入土的速度损失率等8个参数表征夯击点冲击应力的线性回归方程；另一方面采用竖向动应力传递指数对理想弹性静力学理论公式进行修正，得到了基于拟静力法的土体中夯击动应力计算表达式，进而可分析确定强夯有效加固深度。试验结果显示，土-石混合填土的竖向动应力传递指数约为1.673 57，土-石混合填土地层的强夯有效加固深度约为8.0 m。

为获取位于雅鲁藏布江断裂带东段加查县新建水电站工程场区的地应力状态，探讨雅鲁藏布江断裂带东段活动性，采用自主研发的新一代水压致裂地应力测试技术，实测获取了工程场区的地应力状态，并在收集雅鲁藏布江断裂带沿线地应力数据基础上，结合摩尔-库仑破裂准则，对其断层活动性进行了探讨。结果表明：（1）在工程厂区内122.75～418.75 m深度，实测最大水平主应力SH介于6.07～37.62 MPa之间，最小水平主应力Sh介于3.13～20.33 MPa之间，最大水平主应力优势方向为北北东向；（2）主应力关系总体表现为S_H＞S_h＞S_v（S_v为垂向主应力），应力结构有利于逆断层的孕育和活动；（3）实测及收集的大部分数据应力莫尔圆与摩擦系数为0.6的破坏临界线相交，表明整体上雅鲁藏布江断裂带东段发生滑动失稳的风险较高，其中西侧贡嘎至朗县一带风险高于东侧林芝地区。

岩碴图像分析是隧道掘进机（tunnel boring mechine，TBM）智能掘进中实时破岩效率评价和掘进参数优化的重要手段。为了研究掘进参数与岩碴图像特征的关系，依托青岛地铁6号线创石区间TBM施工段，在一段较完整花岗岩地层中开展了分步推力掘进试验，并通过搭载的岩碴图像分析系统对岩碴图像进行采集和分析。结果表明，岩碴图像的中位粒径d50、最大粒径dmax和粗糙度指数CI与推力、扭矩和贯入度之间均表现出正相关，与现场贯入指数FPI和工程比能SE之间均呈现出良好的负相关。当推力超过破岩临界值后，推力的提升将显著提高岩碴的d50和CI，而只有当推力进一步增大到一定水平后dmax才显著增大。此外，随着推力的提升，大岩片的二维形状逐渐呈现出长条形，其平均长短轴比也将增加，但是推力超过一定水平后，大岩片长短轴将保持不变。上述研究结果表明，岩碴粒径和形状参数均反映出破岩效率，可作为实时TBM掘进参数优化的评价指标。

为探究全断面隧道岩石掘进机（tunnel boring machine，TBM）有压输水隧洞内张钢圈-管片-围岩组合结构的受荷状态，以榕江―关埠引水工程中Ⅴ类围岩TBM有压输水隧洞段为研究对象，建立内张钢圈-管片-围岩组合结构三维精细化有限元模型，研究内水压力、围岩类型对组合结构力学特性的影响。研究结果表明：采用内张钢圈加固隧洞结构能有效控制隧洞变形，减小管片拉应力及受拉区范围，提高组合结构的承载能力；内水压力作用下管片压应力、竖向变形、接缝张开度、连接螺栓应力、内张钢圈应力、锚杆应力均较无内水压力作用下有所减小，但管片拉应力和接缝错台分别增大了19.68%、39.25%，其主要原因为内水压力作用致使隧洞结构整体呈现向外膨胀，在充水运营期间应加强内水压力作用下的安全监测；在外部水土压力和内水压力共同作用下，围岩类型的改变使组合结构中连接螺栓应力、管片接缝错台、锚杆应力分别提高37.11%、15.29%、14.75%，其主要原因为围岩岩性越差承担荷载能力越弱，在围岩类型的过渡区域应进行重点监测；此外，在外部水土压力和内水压力共同作用下，组合结构中围岩、管片、内张钢圈、锚杆的荷载分担率分别为21.38%、43.08%、24.01%和11.53%，内张钢圈荷载分担率较无内水压力作用提高了34.06%，内水压力作用提高了内张钢圈分担荷载的效果；组合结构中围岩的分担率随围岩类型（Ⅲ类、Ⅳ类、Ⅴ类）的增加而减小，相同荷载下Ⅲ类围岩的荷载分担率较Ⅴ类围岩提高了16.96%。该研究成果可为类似隧洞工程的衬砌设计及后期加固措施提供理论参考。

不同于常规隧道衬砌结构，无中导洞连拱隧道先行洞衬砌施作稳定后还要受到邻近后行洞开挖的二次扰动，后行洞开挖对先行洞衬砌受力的影响不可忽视。依托杨家寨公路隧道工程，基于考虑了后行洞开挖影响的隧道围岩压力计算方法，结合现场监测位移数据采用荷载-结构法进行位移反分析，提出了一种确定浅埋无中导洞连拱隧道先行洞衬砌承担围岩荷载的方法，最后和衬砌背后围岩荷载现场监测值进行对比验证。研究结果表明：后行洞开挖对先行洞衬砌承担围岩荷载的影响较为显著，监测点增长幅度均值高达142%；先行洞衬砌受力呈现“拱腰较小，拱顶及拱肩较大”的分布规律，且靠近后行洞一侧拱肩位置受力最大；通过现场监测位移反分析得到的先行洞衬砌受力结果和现场监测结果较为一致，论证了该方法的合理性和实用性，研究成果可为无中导洞连拱隧道结构设计提供参考。

随着城市地下空间桩基与隧道立体交叉分布体系的日渐复杂，隧道渗流侵蚀引发的桩基稳定性问题越发突出。然而，目前关于隧道破损引发渗流侵蚀方面的研究工作主要聚焦于侵蚀发生机制方面，其与邻近桩基的相互作用尚欠关注。因此，在物理试验揭示渗流侵蚀机制的基础上，通过建立土体临界细粒含量表达式，量化表征了地层侵蚀形态；随后，采用随机定量删除细颗粒的方法，基于离散单元法DEM建立了多区域侵蚀地层模型，进而分析了不同桩基位置、沉桩方式、桩基荷载和基础类型条件下桩基结构的侵蚀响应特征。结果表明，桩基受侵蚀后均产生不同程度的下沉位移，并在不平衡力的驱动下向侵蚀区倾斜；桩端阻力受损后随桩体下沉逐渐上升，侧摩阻力损失则基本无法恢复。此外，不同沉桩方式与桩基荷载下，桩体的位移与阻力变化模式基本一致，主要是变化量的不同；相比单桩，群桩具有更好的抗侵蚀能力。

高位远程滑坡快速解体后，往往形成碎屑流，会沿程动力侵蚀坡体表层第四纪松散物质，改变滑坡滑移动力学特征，导致滑坡体积显著放大，灾害效应强烈，已成为近年来防灾减灾研究的重点。首先，针对高位堆积体远程滑坡动力侵蚀过程量化预测问题，基于深度积分水动力计算理论与固液两相界面剪应力侵蚀模型架构，构建了高位堆积体远程滑坡动力侵蚀过程量化预测的连续介质数值模型，自主进行了数值实现，编制了具有二阶计算精度的软件程序DisasterFlow V1.0。然后，通过经典的溃坝水流算例验证，结果表明所构建的数值模型的数值求解格式正确并具有较好的总变差不增特性（total variation diminishing, TVD）。再后，通过对柏杨林高位堆积体远程滑坡案例动力侵蚀过程预测可知，所构建的数值模型与传统雪橇模型对该滑坡滑移速度预测结果基本一致，但相较于传统雪橇模型，所构建的数值模型对滑动路径地形的适应能力强，滑移速度预测值振荡性较小。最后，利用所构建的数值模型在时间和空间维度上对柏杨林高位堆积体远程滑坡的动力侵蚀过程进行量化表达，结果表明所构建的数值模型能够对该滑坡动力侵蚀过程中各历时的滑移深度、速度及侵蚀堆积深度进行量化，这是传统雪橇模型所不具备的。所做研究可充实高位远程滑坡动力侵蚀过程量化预测理论与技术体系。

目前针对盾构隧道开挖面稳定性的研究还较少考虑土体分层特性和地震作用的耦合影响，因此构建了一种在分层土体中考虑地震作用的开挖面三维对数螺旋破坏模型进行研究。首先，将地震引起的动态响应通过拟静力法简化为水平和竖直两个方向上的惯性力作用；其次，在均质土体中的三维对数螺旋破坏机制的基础上，将其改进为适用于分层土体的三维对数螺旋破坏机制；再次，根据上限定理，在虚功率方程中引入地震惯性力所做功率，得到考虑土体分层特性和地震作用条件下的盾构隧道开挖面支护力的上限解；最后，将上限理论解与三维数值模拟结果和既有模型试验结果进行对比，得到了较好的一致性。此外，针对水平地震加速度系数和地层厚度对关键物理特征进行了影响因素分析。结果表明，当比例系数z＞0时，极限支护力随水平地震加速度系数的增大显著增大；当z＜0时，极限支护力随水平地震加速度系数增大而增大的趋势减弱。当水平地震加速度系数k_h=0时，即在无地震作用情况下，归一化极限支护力不随z的变化而变化；在上硬下软土层中，下部土层厚度比的增大会引起极限支护力的增大，在上软下硬土层中，下部土层厚度比的增大会引起极限支护力的减小。

为了研究初始孔隙对气孔杏仁状玄武岩力学行为的影响规律，提出了一种考虑计算孔隙率、体积参数、形状参数、角度参数和构造参数的气孔杏仁状玄武岩随机三维模型构建方法，基于建立的随机模型开展了五因素五水平的单轴压缩正交数值模拟，分析不同模型的应力-应变曲线、破坏模式和破坏特征。结果表明：计算孔隙率、角度参数和构造参数对于岩石的单轴抗压强度和弹性模量有较大影响，岩石的单轴抗压强度和弹性模量均与计算孔隙率、角度参数呈负相关，与构造参数呈正相关，与其余因素相关性较弱；随着模型孔隙率的增加，气孔杏仁状玄武岩应力-应变曲线中塑性应变占比随之增加，破坏模式由单剪切破坏面、多剪切破坏面转为局部粉碎状破坏，岩石由脆性破坏逐渐转为延性破坏，5%和12.5%可近似作为强脆性―脆性―延性转换的临界孔隙率；建立的损伤统计本构方程能够根据初始孔隙特征较好地预测气孔杏仁状玄武岩的力学行为，对于玄武岩地层中多孔岩石力学性质的确定具有重要意义。

在降雨和水位升降的联合作用下，库岸边坡易发生滑塌现象。采用理想弹塑性Mohr-Coulomb准则，在干湿循环作用下，分析库岸边坡的稳定性，难以反映边坡土的超固结消散等复杂力学特性。以江西省新干航电枢纽库岸边坡为背景，针对弱超固结非饱和红黏土，开展非饱和直剪试验，采用反正切函数，描述吸力应力与基质吸力的数学物理关系，提出红黏土的超固结非饱和统一硬化（unified hardening, UH）模型。然后，在UH模型的基础上，采用C++语言，开发FLAC3D的应用程序；利用GeoStudio软件的SEEP/W模块，计算降雨入渗的非饱和渗流场；编写接口程序，导入FLAC3D，开展库岸边坡的稳定性计算。库岸边坡的水平位移明显大于经典的非饱和理想弹塑性模型计算结果，结果表明，修正UH模型能够合理地预测超固结非饱和红黏土库岸边坡的稳定性，为同类库岸边坡工程稳定性分析提供了参考。

粒子冲击辅助破岩技术凭借其快速、高效等优点，对硬质岩体有着较好的破岩效果。通过粒子冲击试验和离散元模拟相结合的方法研究单、双粒子冲击速度、双粒子间距等因素对高强度花岗岩表面、三维及剖面形貌特征的影响，探寻冲击坑深度、冲击坑体积及冲击坑表面面积随冲击参数的变化规律，统计粒子冲击破岩裂纹的分布规律，并且从能量吸收率的角度评价双粒子迟滞冲击破岩的效能。结果表明：冲击坑深度与冲击速度呈正相关；随着粒子间距的增大，冲击坑由相交逐渐相离，形貌随之变化，冲击坑体积随之减小，而冲击坑表面面积增大；通过模拟发现，裂纹主要分布在斜长石与钾长石的晶界处，以拉伸破坏为主；当选用5 mm直径的钢粒子破碎强度200 MPa左右的极坚硬花岗岩时，双粒子迟滞冲击破岩的能量吸收率曲线随粒子冲击速度增大趋于平缓，在双粒子迟滞冲击破岩的粒子间距为8～10 mm且冲击速度400 m/s左右时能达到最佳的冲击辅助破岩效果。

网格分布形式对三维极限分析下限有限元法的计算精度影响较大，为获取精确下限解，通常需对破坏区域进行密集网格划分，从而极易导致计算规模过大，求解效率低下。针对上述问题，提出一种基于单元应力的三维极限分析下限有限元“后验”网格自适应加密策略。首先，构建基于Mohr-Coulomb（M-C）准则和半定规划技术的三维极限分析下限有限元模型，避免了屈服准则的近似处理。其次，引入基于M-C准则的网格自适应加密策略，通过判断各个单元应力接近屈服的程度，确定加密点的坐标。然后，将加密点与原先节点组合构成新的点集，并重新划分网格，建立新的下限有限元计算模型。最后，利用所提方法研究隧道稳定性问题，表明利用所提网格自适应加密策略能够以较少单元精确模拟破坏区域应力分布，从而获取高精度下限解。