A receptance based method for the calculation of nonlinear normal modes of large ordered structures with distributed localized nonlinearities. (December 2022)