Analysis of neuromuscular fatigue finds various applications ranging from clinical studies to biomechanics. Surface electromyography (sEMG) signals are widely used for these studies due to its non-invasiveness. During cyclic dynamic contractions, these signals are nonstationary and cyclostationary. In recent years, several nonstationary methods have been employed for the muscle fatigue analysis. However, cyclostationary based approach is not well established for the assessment of muscle fatigue. In this work, cyclostationarity associated with the biceps brachii muscle fatigue progression is analyzed using sEMG signals and Spectral Correlation Density (SCD) functions. Signals are recorded from fifty healthy adult volunteers during dynamic contractions under a prescribed protocol. These signals are preprocessed and are divided into three segments, namely, non-fatigue, first muscle discomfort and fatigue zones. Then SCD is estimated using fast Fourier transform accumulation method. Further, Cyclic Frequency Spectral Density (CFSD) is calculated from the SCD spectrum. Two features, namely, cyclic frequency spectral area (CFSA) and cyclic frequency spectral entropy (CFSE) are proposed to study the progression of muscle fatigue. Additionally, degree of cyclostationarity (DCS) is computed to quantify the amount of cyclostationarity present in the signals. Results show that there is a progressive increase in cyclostationary during the progression of muscle fatigue. CFSA shows an increasing trend in muscle fatiguing contraction. However, CFSE shows a decreasing trend. It is observed that when the muscle progresses from non-fatigue to fatigue condition, the mean DCS of fifty subjects increases from 0.016 to 0.99. All the extracted features found to be distinct and statistically significant in the three zones of muscle contraction (p < 0.05). It appears that these SCD features could be useful in the automated analysis of sEMG signals for different neuromuscular conditions. © 2015, Springer Science+Business Media New York.