An analytical method is presented in this paper for simulating piston secondary dynamics and piston-bore contact for an asymmetric half piston model including elastohydrodynamic (EHD) lubrication at the bore-skirt interface. A piston EHD analysis is used based on a finite-difference formulation. The oil film is discretized using a two-dimensional mesh. For improved computational efficiency without loss of accuracy, the Reynolds’ equation is solved using a perturbation approach which utilizes an “influence zone” concept, and a successive over-relaxation solver. The analysis includes several important physical attributes such as bore distortion effects due to mechanical and thermal deformation, inertia loading and piston barrelity and ovality. A Newmark-Beta time integration scheme combined with a Newton-Raphson linearization, calculates the piston secondary motion. The EHD behavior is coupled with a contact algorithm at the piston-bore interface as well as the piston secondary dynamics. The effects of some piston design parameters such as piston pin offset, clearance and location of cylinder peak pressure are investigated.