On the basis of continuum damage mechanics theory, the effective stress theory of damaged porous media was presented and its modified selfcontained seepage differential equation was deduced. The traditional DruckerPrager yield criterion subsequently was modified, and a coupled seepage and damage finite element computer program compiled that was combined with the modified Biot consolidation equation. A simple roadbed foundation model then was analyzed and calculated using the computer program we compiled, and the typical node settlement, pore pressure, and damage variables were discussed and compared with those from the analytic solution. The results show that the damage finite element program we compiled is reasonable and effective.