Numerical studies of thermal effects on the on-state characteristics of IGBTs have been carried out using two-dimensional simulation. A discretisation scheme suitable for electrothermal modelling of semiconductor devices has been derived. The coupled Poisson's equation and continuity equations together with the heat flow equation are solved self-consistently. The simulation results show that the on-resistance and hence current handling capability of IGBTs can be significantly degraded by their self-heating.