The current study extends the applications of the variational iteration
method for the analytical solution of fractional contact problems. The
problem involves Caputo sense while calculating the derivative of fractional
order, we apply the Penalty function technique to transform it into a system
of fractional boundary value problems coupled with a known obstacle. The
variational iteration method is employed to find the series solution of
fractional boundary value problem. For different values of fractional
parameters, residual errors of solutions are plotted to make sure the
convergence and accuracy of the solution. The reasonably accurate results
show that one of the highly effective and stable methods for the solution of
fractional boundary value problem is the method of variational iteration.