In this work, a novel method for damage evolution analysis of concrete under uniaxial compression is proposed based on the multi-feature acoustic emission (AE) and the Gaussian mixture model (GMM) clustering. The hierarchical clustering algorithm is adopted to select optimal AE parameters, while multiple features of these parameters are generated through the principal component analysis (PCA). Then the concrete damage signals are separated by using the GMM clustering with multiple features. Meanwhile, the waveform signals associated with concrete damage in each cluster are validated using the Fast Fourier transform (FFT) and the continuous wavelet transform (CWT). Finally, the damage evolution process of concrete under uniaxial compression is divided based on the variation of AE characteristics. The results show that the frequency ranges of signals for micro-cracks, mixed cracks and friction of concrete under uniaxial compression are 20–30 kHz, 20–65 kHz, and 100–120 kHz, respectively. The damage process of concrete under uniaxial compression is divided into three stages according to the inflection points of the cumulative energy curve, namely the initial compaction stage, the crack formation stage and the failure stage.