The engine speed fluctuation at idle operation mainly comes from cyclic variation of combustion in SI engine. In the present study, engineering model that is representing the cyclic variation of combustion was proposed for the sub-model of the engine cycle simulation. From the observed behaviors of the mass burn rates, probability density functions for the parameters of Wiebe function were defined. The mass burn rate of each cycle is obtained by Monte Cralo perturbation method with the probability function. The simulation results shows that trends of cylinder pressure variation and imep distribution follow up with those of experimental results at idle condition.