Model reduction by moment matching with preservation of global stability for a class of nonlinear models