With the growing demand for reliability and efficiency in heavy machinery, this study establishes a bearing simulation model that accounts for elastic deformation using the finite difference method and finite element method. Fortran programs were employed to conduct numerical simulations on large-diameter, heavily-loaded journal bearings, yielding performance parameters such as bearing pressure field, temperature field, and oil film thickness. By comparing these results with simulations that do not consider elastic deformation, it was found that the maximum oil film pressure, maximum oil film temperature, and maximum temperature rise are all reduced. The minimum oil film thickness shows minimal variation due to changes in the bearing's equilibrium position.