This paper presents a numerical simulation of hydrodynamic journal bearing lubrication by using finite element method to solve Reynolds equation in static load condition. Reynolds boundary condition is applied in order to compute oil film pressure distribution. Then determine the equilibrium position of the housing bearing by solving the equilibrium equation of the charge using Newton-Raphson method. The equilibrium positions are simulated in different values of the journal speed and the applied load. The results show that at the different sections of bearing along the axial direction, the starting disruption positions are different, at the middle section, the pressure reaches the maximum value and decreases toward two ends of bearing. The more load applied, the housing bearing center position get farther from the journal center. The higher journal rotation speed makes the balance point closer to the journal center.