Over the past two decades, extensive work has been conducted on the dynamic effect of joint clearances in multibody mechanical systems. In contrast, little work has been devoted to optimizing the performance of these systems. In this study, the analysis of revolute joint clearance is formulated in terms of a Hertzian-based contact force model. For illustration, the classical slider-crank mechanism with a revolute clearance joint at the piston pin is presented and a simulation model is developed using the analysis/design software MSC.ADAMS. The clearance is modeled as a pin-in-a-hole surface-to-surface dry contact, with an appropriate contact force model between the joint and bearing surfaces. Different simulations are performed to demonstrate the influence of the joint clearance size and the input crank speed on the dynamic behavior of the system with the joint clearance. In the modeling and simulation of the experimental setup and in the followed parametric study with a slightly revised system, both the Hertzian normal contact force model and a Coulomb-type friction force model were utilized. The kinetic coefficient of friction was chosen as constant throughout the study. An innovative design-of-experiment (DOE)-based method for optimizing the performance of a mechanical system with the revolute joint clearance for different ranges of design parameters is then proposed. Based on the simulation model results from sample points, which are selected by a Latin hypercube sampling (LHS) method, a polynomial function Kriging meta-model is established instead of the actual simulation model. The reason for the development and use of the meta-model is to bypass computationally intensive simulations of a computer model for different design parameter values in place of a more efficient and cost-effective mathematical model. Finally, numerical results obtained from two application examples with different design parameters, including the joint clearance size, crank speed, and contact stiffness, are presented for the further analysis of the dynamics of the revolute clearance joint in a mechanical system. This allows for predicting the influence of design parameter changes, in order to minimize contact forces, accelerations, and power requirements due to the existence of joint clearance.