刘兴荣，魏新平，陈豫津，王翔宇

（1.甘肃省科学院 地质自然灾害防治研究所，甘肃 兰州 730000；2.甘肃省地质环境监测院，甘肃 兰州 730050；3.中国石油西南管道贵阳输油气分公司，贵州 贵阳 550000）

2010年8月8日凌晨，舟曲县城北侧三眼峪沟和罗家峪沟同时暴发特大山洪泥石流，城区三分之一被淹，共造成1 435 人死亡，330 人失踪，直接经济损失超过10 亿元[1-2]，给舟曲县城居民生命财产造成了巨大损失，也给当地生产生活带来严重困难。灾害引起了党中央、国务院、中央军委及全国人民的高度关注，启动国家二级救灾应急响应，同时批复专项资金进行治理。其中，重力式拦挡坝是舟曲泥石流治理中的最主要工程之一。

周龙茂等[3]认为拦挡坝在泥石流治理中发挥着重要作用，但同时拦挡坝也是最易遭受破坏、失去防灾功能的泥石流防治构筑物，因此，在泥石流设计中，为了拦挡坝不被破坏，往往设计得很保守，王念秦等[4]提出这种设计容易造成两个极端现象：保守，造成资金浪费；冒进，防治工程失败。要做到既能保证拦挡坝安全，又能将投资最小化，就要求对泥石流拦挡坝抗冲击力验算方法提出新要求。传统的泥石流冲击力计算经验公式[5]只能通过大量试算表述结果，不能表述过程。而将三维有限元数值分析方法应用到泥石流拦挡坝稳定性验算中[6-7]，能将过程和结果同时呈现，即通过分步加载的方法，逐步呈现拦挡坝的位移情况和抗冲击力过程中的破损情况。

关于拦挡坝的数值模拟研究还比较少，本文将考虑损伤的混凝土本构模型与有限元计算方法相结合，对舟曲泥石流混凝土拦挡结构进行力学分析，最终确定了拦挡坝的抗冲击力合理区间，以期为泥石流治理工程的设计提供借鉴。……