In this study, a comprehensive model for the numerical simulation of electrical and optical characteristics of multilayer organic light-emitting devices (OLED) is presented. The model consists of three parts: charge carrier transport model, exciton model, and a post-processing part for calculating luminance. The defect states in the organic layers are modeled and are introduced into the equations. For demonstrating the rule of the defect states in such devices, a multilayer OLED is considered consisting of a host–guest emissive layer sandwiched between electron and hole transport layers. The current density and luminance results of the analyzed device are compared with the experimental data. We have found that by considering the defect states, there is good accordance between the simulation results and the experimental data.