In this paper, contact modeling technique and dynamics analysis of piston and cylinder system are presented by using modal synthesis method. It is very important to select mode shapes representing a global or local behavior of a flexible body due to a specified loading condition. This paper proposes a technique to generate the static correction modes which are nicely representing a motion by a contact force between a piston and cylinder. First normal modes of piston and cylinder under a boundary condition are computed, and then static correction modes due to a contact force applied at contacted nodes are added to the normal modes. Also, this paper proposes an efficient dynamics analysis process while changing the shape of the piston and cylinder. In optimization process or design study, their geometric data can be changed a bit. The slight changes of their contact surfaces make a high variation of the magnitude of a contact force, and it can yield the different dynamic behavior of an engine system. But, since the variations of the normal and correction modes are very small, the re-computation of their normal and correction modes due to the change of contact surfaces can be useless. Until now, whenever their contact surfaces are changed at a design cycle, the modes have been recomputed. Thus, most engineers in industries have been spent many times in very tedious and inefficient design process. In this paper, the normal and correction modes from the basic geometry of the piston and cylinder are computed. If the geometry shape is changed, nodal positions of the original modal model are newly calculated from an interpolation method and changed geometry data. And then the updated nodes are used to compute a precise contact force. The proposed methods illustrated in this investigation have good agreement with results of a nodal synthesis technique and proved that it is very efficient design method.