Estudiamos las propiedades termodinámicas del modelo generalizado de Curie-Weiss multispecies no convexo, donde las interacciones entre diferentes tipos de partículas (formando las especies) están codificadas en una matriz genérica. Para espines con una distribución previa genérica, calculamos el límite termodinámico de la función generadora para los momentos de la medida de Boltzmann-Gibbs utilizando técnicas de interpolación simples. Para espines de Ising, analizamos además las fluctuaciones de la magnetización en el límite termodinámico bajo la medida de Boltzmann-Gibbs. Se muestra que un teorema del límite central (CLT) se cumple para un vector de magnetizaciones de especies escalado y centrado, que converge a una distribución normal multivariante centrada o no centrada, dependiendo de la tasa de convergencia de los tamaños relativos de las especies.