This paper presents a 3D-modeling numerical method using finite element method (FEM) to simulate the electromagnetic behavior of high-temperature superconductors (HTSC). The models are formulated by the magnetic field vector method (H-method). The resolving code was written by FROTRAN language. The electromagnetic properties of HTSC are described though Kim critical-state model. The magnetic fields and current distribution in the bulk HTSC in the applied non-uniform external magnetic fields generated by the permanent magnetic guideway (PMG) are obtained using the proposed method. The magnetic levitation forces by the interaction between the bulk HTSC and the PMG are calculated. In order to validate the method, measurement of the vertical force between a bulk YBaCuO(YBCO) and a PMG is obtained. The measurement and simulation results show good matching. This method could be used in the HTSC magnetic levitation transportation system optimization design.