Multi phase-field modeling of anisotropic crack propagation in 3D fiber-reinforced composites based on an adaptive isogeometric meshfree collocation method