Interlaminar damage model for polymer matrix composites