Methodology for crack density evolution simulation in tension-tension fatigue and parameter determination in a model that relies on failure stress distribution (Weibull) in the 90 ply is presented. Cyclic loading is performed at one stress level to obtain detailed crack density dependence on the number of cycles. In addition, one data point (crack density at specified number of cycles) is necessary at a different stress level. Non-interactive crack density region is used to determine Weibull parameters. Then, crack density in a whole crack density range and for arbitrary stress level is predicted using the obtained Weibull parameters, the Monte Carlo method for failure stress distribution and a novel model for stress distribution between cracks. The predictions are in good agreement with test results. © 2022 Pakkam Gabriel et al.