In this paper, an analytical and closed-form solution to the lumped parameter equivalent circuit model of organic solar cells is proposed to complete the simulations of the S-shaped I-V characteristics. Based on the model previously proposed by Mazhari, the set of terminal current and voltage equations describing the three diodes is solved and the effects from the model parameters are illustrated. Our solutions are verified by being compared with the least square method results and experimental data, respectively. Good agreements show that our solution calculation scheme is not only both accurate and efficient, but also valid in the whole operation regime of solar cells, especially for the S-shaped kink on the condition where the terminal voltage is larger than the open circuit voltage. Such an analytical solution can play an important role in the simulations for I-V characteristics of solar cells, fast extractions of the model parameters, and implements into practical photovoltaic device simulators.