Based on the equations of electrodynamics and the concept of a critical state for hard superconductors of the 2nd kind, numerical simulation of the magnetic properties of axisymmetric superconducting samples, in particular, granules, is performed for a number of models of the dependence of the critical current density on the magnetic field induction. The magnetic moment loops are calculated directly by integrating the integral equation for the current density over time. The phenomena of the peak effect and the asymmetry of the magnetization hysteresis loop are also considered using the indicated equation. Various versions of the functions used in the literature were used as peak functions. In addition to the hysteresis loop of the magnetic moment, the magnetic field induction at the center of axisymmetric samples, and the total penetration field, the profiles of the critical current density Jc(B) and the equilibrium magnetic moment for spherical granules were obtained. The method used for calculating the magnetic moment of superconductors makes it possible to take into account the equilibrium and nonequilibrium regions of the magnetization of the samples independently Keywords: magnetic moment, axisymmetric granules, type II superconductor, critical state, peak effect, equilibrium magnetic moment.