The paper presents the numerical modeling of a Y-shaped three-pole radial magnetic bearing based on two-dimensional (2D) and three-dimensional (3D) magnetic field computation with nonlinear model of the material. The used numerical method is the Finite Element Method (FEM). The nonlinear system of equations according to the nonlinear characteristics of ferromagnetic material can be handled by the Newton-Raphson technique and by the fixed-point method.