Physical Phenomena’s located around us are primarily nonlinear in nature and their solutions are of highest significance for scientists and engineers. In order to have a better representation of these physical models, fractional calculus is used. Fractional order oscillation equations are included among these nonlinear phenomena’s. To tackle with the nonlinearity arising, in these phenomena’s we recommend a new method. In the proposed method, Picard’s iteration is used to convert the nonlinear fractional order oscillation equation into a fractional order recurrence relation and then Legendre wavelets method is applied on the converted problem. In order to check the efficiency and accuracy of the suggested modification, we have considered three problems namely: fractional order force-free Duffing–van der Pol oscillator, forced Duffing–van der Pol oscillator and higher order fractional Duffing equations. The obtained results are compared with the results obtained via other techniques.