The purpose of the study is to find development of contact area, contact pressure and friction forces occurring at joints or connection areas inbetween structural members or mechanical parts. The problem has a pair of difficulties intrinsically; a constraint of displacement due to contact, and presence of work term by nonconservative friction force in the variational principle of the problem. Because of these difficulties, the variational principle remains in the form of inequality. It is resolved by penalty method and perturbation method making the inequality to an equality which is proper for computational purposes. A contact problem without friction is solved to find contact area and contact pressure, which are to be used as data for the analysis of the friction problem using perturbed variational principle. For numerical experiments, a Hertz problem, a rigid punch problem, and the latter one with friction effects are solved using $Q_2$-finite elements.