A trapped field magnet is one application for melt-processed bulk high-Tc superconductors (HTSs). Evaluation of the maximum stress during magnetization is important from the viewpoint of the destruction of bulk HTS. Stresses in a bulk HTS are numerically evaluated during the magnetization process by field cooling. Shielding current distributions are obtained through macroscopic numerical simulation using Maxwell equations and the critical-state model. In the study of maximum stress, 3D axisymmetirc numerical solutions with open and fixed boundary conditions are compared with 1D analytical solutions. Analytical solutions are different from numerical solutions. It is shown that the evaluation of experimental results on the destruction of a bulk HTS will need numerical analysis.