Abstract: In this paper, solutions of three-dimensional lubrication problem of bearing are presented, The basic equations are converted into boundray integral equations using Green's function me thod. Spatial triangular boundary elements and the isopara metric elements with six nodes are used to discretize the boundary of the solution domain. The integral equations are solved numerically. The velocity and pressure fields and the surface stress are determined. When we are compared with analytic solutions, this method has higher accuracy. Some satisfactory computation results for other typical problems are indicate the affectiveness of the method for numerical analysis of these kinds of problems.