A Bayesian model selection approach is presented for selecting the most robust model class among a set of physics-based models for fatigue damage prognostics in composites. Three families of micro-damage mechanics models (shear-lag, variational and crack-opening displacement) are chosen to capture the relationship between stiffness reduction and matrix-cracks density in composite fiber-reinforced polymers (CFRP). First, the candidate models are parameterized by global sensitivity analysis, and then, they are ranked through probabilities that measure the extent of agreement of their predictions with observed SHM data, while avoiding the extremes of over-fitting or under-fitting. These probabilities are computed based on the evidence provided by the data and the modeler’s choice of prior probability for each model class. A case study is presented using multi-scale fatigue damage data from a cross-ply…