A numerical approach for limit analysis of orthotropic composite laminates