Prediction of composite fatigue life is not a straightforward matter, depending on various failure modes and their interactions. In this paper, a methodology is presented to predict fatigue life and residual strength of composite materials based on Phenomenological Model(non-linear fatigue damage model). It is assumed that the residual strength is a monotonically decreasing function of the number of loading cycles and applied fatigue stress ratio and the model parameters(strength degradation parameter and fatigue shape parameter) are assumed as function of fatigue life. Then S-N curve is used to extract model parameters that are required to characterize the stress levels comprising a randomly-ordered load spectrum. Different stress ratios (${sigma}_{min}/{;}{sigma}_{max}$) are handled with Goodman correction approach(fatigue envelope) and the residual strength after an arbitrary load cycles is represented by two parameter weibull functions.