Abstract: In this paper the no-interpenetration condition between the two structures is treated as a penalty term and then the cons-trained variational problem of the contact system will be transformed into a penalty-optimum one in which the substruction finite element analysis is applied The advantage of the method is that the iterative procedure will be only related to the displacements of the nodes and be irrelative to the nodal forces therefore the calculator time is reduced relatively.