Deformable image registration is an important tool for image-guided radiotherapy. Physics-model-based deformable image registration using finite element analysis is one of the methods currently being investigated. The calculation accuracy of finite element analysis is dependent on given boundary conditions, which are usually based on the surface matching of the organ in two images. Such a surface matching, however, is hard to obtain from medical images. In this study, we developed a new boundary condition to circumvent the traditional difficulties. Finite element contact-impact analysis was employed to simulate the interaction between the organ of interest and the surrounding body. The displacement loading is not necessarily specified. The algorithm automatically deforms the organ model into the minimum internal energy state. The analysis was performed on CT images of the lung at two different breathing phases (exhalation and full inhalation). The result gave the displacement vector map inside the lung. Validation of the result showed satisfactory agreement in most parts of the lung. This approach is simple, operator independent and may provide improved accuracy of the prediction of organ deformation.