A novel numerical simulation method for predicting compressive properties of 3D multiaxial braided composites considering various defects