The high-density behavior of the stellar matter composed of nucleons and leptons under β-equilibrium and charge neutrality conditions is studied with the Skyrme parametrizations shown to be consistent (CSkP) with the nuclear matter, pure neutron matter, symmetry energy and its derivatives in a set of 11 constraints [Dutra et al., Phys. Rev. C 85, 035201 (2012)]. The predictions of these parametrizations on the tidal deformabilities related to the GW170817 event are also examined. The results points out to a correlation between the Love numbers and tidal deformabilities with the respective radii of the binary neutron stars system (BNSS). We also find that those CSkP supporting massive neutron stars (M ∼ 2M⊙) predict radii of the BNSS in full agreement with recent data from LIGO and Virgo Collaboration (LVC) given by R1 = R2 = 11.9 +1.4 −1.4 km. A correlation between dimensionless tidal deformability and radius of the canonical star is found, namely, Λ1.4 ≈ 5.87 × 10 −6 R 7.19 1.4 , with results for the CSkP compatible with the recent range of Λ1.4 = 190 +390 −120 from LVC. Finally, an analysis of the Λ1 ×Λ2 graph shows that the CSkP compatible with the recent bounds obtained by LVC, namely, GSkI, Ska35s20, MSL0 and NRAPR, can also support massive stars (M ∼ 2M⊙), and predict a range of 11.83 km R1.4 12.11 km for the canonical star radius.