本發(fā)明公開了巖石力學(xué)工程技術(shù)領(lǐng)域的一種考慮損傷滲流應(yīng)力耦合的數(shù)值模擬方法，旨在真實(shí)模擬巖體工程在施工運(yùn)營過程中所受應(yīng)力情況，并能充分考慮地質(zhì)環(huán)境等因素影響下的巖體變形行為。通用數(shù)值軟件內(nèi)置本構(gòu)模型很難滿足實(shí)際仿真的需要，因此如何將自定義的本構(gòu)模型應(yīng)用到數(shù)值計算已成為熱點(diǎn)和難點(diǎn)問題。所述方法包括如下步驟：巖石滲流應(yīng)力耦合損傷模型建立；建立損傷滲流應(yīng)力耦合流變增量本構(gòu)關(guān)系；廣義Kelvin損傷黏彈?塑性模型及其差分格式；D?P準(zhǔn)則二次開發(fā)。本發(fā)明解決了現(xiàn)有技術(shù)中無法考慮損傷滲流應(yīng)力耦合流變問題。