冰川碎屑冰冻体的冰水赋存状态对温度的变化非常敏感。近年受到气候变暖影响，由碎屑冰冻升温融化而引发的冰崩灾害被广泛报道。为了研究温升条件下碎屑冰冻体的抗剪强度劣化行为，提出了正融碎屑冰冻体的有限-离散元数值模型（F-DEM）。模型将碎屑冰冻体概化为冻土、碎石实体单元以及粘聚力单元，将根据实验所得的碎屑冰冻体抗剪强度升温劣化描述为粘聚力单元的粘结强度的升温劣化过程。首先，在冻土内预置使用“牵引-分离”准则的粘聚力单元来代表孔隙冰。然后，在ABAQUS中使用用户自定义子程序VUSDFLD，通过温度场变量对粘聚力单元的粘结强度劣化规律进行控制。将数值结果与实验进行对比，发现峰值抗剪强度、变形模式、破坏模式等宏观响应特征和强度劣化规律与试验结果吻合较好。对融冰率、碎石率和载荷水平等条件对碎屑冰冻体剪切力学行为的影响进行了研究，结果表明当融冰率小于等于2%，剪切破坏模式为粗糙的“锯齿状”；随着融冰率增加（>2%），剪切破裂面逐渐向较为平滑的“圆弧状”转变。在冻土-碎石界面强度差异较大容易形成应力集中，导致裂缝沿该界面形成渐进式贯通。含石率的增加会导致碎屑冰冻体抗剪强度的降低。随融冰率增加，碎屑冰冻体抗剪强度对含石率的敏感性降低。对于维持着较高剪切负荷的碎屑冰冻体来说，环境温度的轻微上升可能导致剪切应变的突变。恒载下的变形可分为三个阶段：初始阶段，发育阶段以及快速变形阶段。在初始阶段，剪切应变先增加后稳定；在发育阶段，应变出现拐点；在快速变形阶段剪切应变快速增加。拐点的温度随着初始剪切应力的增加而降低，两者近似为线性关系。较高的应力水平下，碎屑冰冻体的剪切应变对温度的改变十分敏感。未来应就模型简化方法、参数取值规律、冰水相变和尺寸效应等问题进行进一步研究，以适应实际冰川案例中剪切强度劣化行为的模拟。

1 基于改进F-DEM的强度升温劣化行为模拟
1.1 F-DEM模型说明
1.2 强度劣化控制方程
1.3 数值验证
2 数值结果
2.1 常应变率模型
2.2 恒定载荷模型
3 讨论
4 结论