从多孔介质的角度,首先建立岩石孔隙度与体积应变之间的关系,在此基础上定义材料的渗透系数、孔隙度、损伤变量和饱和度的基于孔隙度变化的动态演化关系式;然后基于多孔介质有效应力原理,提出非饱和渗流应力耦合的弹塑性损伤本构模型;最后将所建模型通过子程序导入大型非线性有限元分析软件ABAQUS,对锦屏水电站引水隧洞开挖过程中的围岩变形、应力、塑性区、孔隙水压力、损伤及孔隙度演化进行计算与分析,研究结果对引水隧洞的开挖与支护设计具有一定的参考价值和借鉴意义。数值计算结果表明,所建模型合理,可反映围岩饱和–非饱和–再饱和过程中变形、损伤演化的机制。