Among various ceramic matrix composites developed, self-healing ceramics have been studied as new functional materials. Self-healing occurs in such materials by high-temperature oxidation triggered by a micro-crack initiation on the surface, and the strength of the material autonomously recovers to its robust state since the micro-crack is re-bonded. To facilitate the use of self-healing ceramics in machines and equipment, a novel numerical simulation method based on finite element analysis (FEA) needs to be applied. In this study, we applied a previously proposed constitutive model to a series of self-healing and damage processes. In the constitutive model, the damage process is formulated on the basis of fracture mechanics, while the self-healing process is formulated on the basis of empirical oxidation kinetics. The FEA model implemented the constitutive model to simulate a series of experiments of the alumina/15 vol% SiC composites. The self-healing process was targeted to a prescribed damage by Vickers indentation. Thereafter, the self-healing behavior was quantitatively compared with that observed in the experiment. The results suggest that the proposed FEA approach can be applied to the analysis of ceramic matrix composites with self-healing properties.