We study the elastic time-harmonic wave scattering problems on unbounded domains with boundaries composed of finite collections of disjoints finite open arcs (or cracks) in two dimensions. Specifically, we present a fast spectral Galerkin method for solving the associated weakly-and hyper-singular boundary integral equations (BIEs) arising from Dirichlet and Neumann boundary conditions, respectively. Discretization bases of the resulting BIEs employ weighted Chebyshev polynomials that capture the solutions' edge behavior. We show that these bases guarantee exponential convergence in the polynomial degree when assuming analyticity of sources and arcs geometries. Numerical examples demonstrate the accuracy and robustness of the proposed method with respect to number of arcs and wavenumber.