Abstract
A Galerkin's finite element approach based on weighted-residual formulation is presented to find approximate solutions to obstacle, unilateral and contact second-order boundary-value problems. The approach utilizes a piece-wise linear approximations utilizing linear Langrange polynomials. Numerical studies have shown the superior accuracy and lesser computational cost of this scheme in comparison to collocation, finite-difference and spline methods.