Terahertz Kerr nonlinearity analysis of a microribbon graphene array using the harmonic balance method