For multileaf foil bearings,the interactions between adjacent foils, foils and gas film are key to predict the bearing performance. However, most of the analytical models were based on line contact rather than area contact between adjacent foils and the latter may occur and play an important role in load capacity. In this study, the curved beam model is adopted to build up the elastohydrodynamic model considering the area contact between adjacent foils, which is solved iteratively coupling the Reynolds equation by the finite element method. An algorithm is developed to judge and determine the contact condition of adjacent foils. The simulation results show that the area contact between adjacent foils does occur and the contact area expands with the increasing of eccentricity ratios and rotational speeds, which make the film thickness distribution smooth and are beneficial to the increase of load capacity. The area contact occurs near the trailing end of the foils first, and then gradually expands. The later formed contact area bears greater load capacity instead of the early formed contact area, and the load capacity is mainly provided by the gas film on the foils with more area contact.