An interference fit is a fastening between two parts which is achieved by friction after the parts are pushed together, rather than by any other means of fastening. A fit is depending on which part has its size controlled to determine the fit. To select an amount of interference see KS standard tables for class from light to heavy drive fits. The aim of this paper is to find the optimal gap under the interference fit with lowest contact pressure to the two cylinders. To accurately determine the amount of interference clearance, must consider the fluctuation of the response parameters such as contact pressure, displacements and stresses. The response of the cylinder system is modelled with the probabilistic finite element method using the Monte Carlo simulation. The probabilistic design is carried out using ANSYS probabilistic design software. And then the optimal design of the interference gap is sequentially solved to find the solution. The practical optimal design is proposed for the design of the interference fit system. The numerical results are obtained where the displacements and stresses treated as constraints.