Abstract
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.