Focused on the sensitivity to climate change and the special mechanical characteristics of undisturbed expansive soil, an elasto-plastic damage constitutive model was proposed based on the mechanics of unsaturated soil and the mechanics of damage. Undisturbed expansive soil was considered as a compound of non-damaged part and damaged part. The behavior of the non-damaged part was described using non-linear constitutive model of unsaturated soil. The property of the damaged part was described using a damage evolution equation and two yield surfaces, i.e., loading yield (LY) and shear yield (SY). Furthermore, a consolidation model for unsaturated undisturbed expansive soil was established and a FEM program named UESEPDC was designed. Numerical analysis on solid-liquid-gas tri-phases and multi-field couple problem was conducted for four stages and fields of stress, displacement, pore water pressure, pore air pressure, water content,suction, and the damage region as well as plastic region in an expansive soil slope were obtained.