In this paper, we solve the three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind using the full multigrid (FMG) algorithm. We derive a finite difference approximations for the biharmonic equation on a 18 point compact stencil. The unknown solution and its second derivatives are carried as unknowns at grid points. In the multigrid methods, we use a fourth order interpolation to producing a new intermediate unknown functions values on a finer grid, and the full weighting restriction operators to calculating the residuals at coarse grid points. A set of test problems gives excellent results.