Based on an efficient numerical integration, called LEBEDEV quadrature, a fast ellipsoid shape model for asteroid is presented to simulate the variation of observed brightness, i.e. lightcurves. The fast method can simultaneously derive the shape model with other physical parameters, such as rotational period, orientation of spin axis in an efficient way. Furthermore, in order to simulate the asymmetric surface of asteroids, a more general shape model, called ‘CELLINOID’, is applied to simulate the asymmetric lightcurves with different extreme values. The shape model consists of eight octants of ellipsoids having different semi-axes, with the constraint that adjacent octants must have two equal semi-axes in common. Simply extended from ellipsoid shape model, CELLINOID model can derive the shape and other physical parameters more accurately than ellipsoid model. Additionally, a convex method is introduced with the acceleration based on Lebedev quadrature. Employing these lightcurves, the shape model and physical parameters of asteroids can be refined.