Damages in composite materials are mostly governed by interfacial debonding. In order to take this into account, we employ an analytical model considering partial debonding along interface, and formulate a finite element. The model utilizes a concept where a moduli of an inclusion is replaced by an equivalent one of debonded inclusions. By several examples, eligibility of the model is shown. As an application, excavation analysis of a tunnel in jointed rock mass is carried out, and the results are compared with other observations.