We present a theoretical model to analyse the propagation of a Gaussian laser beam through double-sided nonlinear media. This model is based on the Huygens-Fresnel diffraction integral method. This theoretical model is not only consistent with the cascade structure model for a small nonlinear phase-shift but also can be used for a large nonlinear phase-shift. It has been verified that it is suitable to characterize the double-sided nonlinear media compared with the cascade structure model. A good agreement between the experimental data and the results from the theoretical model is obtained. It will be useful for the design of multi-sided nonlinear materials.