We investigate the CP-averaged branching ratios, the polarization fractions, the relative phases, and the CPviolating asymmetries of the penguin-dominated B → φK1(1270) and φK1(1400) decays in the perturbative QCD(pQCD) approach, where K1(1270) and K1(1400) are believed to be the mixtures of two distinct types of axial-vector K1A( 3 P1) and K1B( 1 P1) states with different behavior, however, their mixing angle θK 1 is still a hot and controversial topic presently. By numerical evaluations with two different mixing angles θK 1 ∼ 33• and 58• and phenomenological analysis, we find that: (a) the pQCD predictions for the branching ratio, the longitudinal polarization fraction and the direct CP violation of B ± → φK1(1270) ± decay with the smaller angle 33• are in good agreement with the currently available data; (b) though the central values significantly exceed the available upper limit, both pQCD predictions of Br(B ± → φK1(1400) ± ) with two different mixing angles are consistent with that obtained in QCD factorization and with the preliminary data in 2σ errors. These results and other relevant predictions for the considered decays will be further tested by the LHCb and the forthcoming Super-B experiments; (c) the weak annihilation contributions can play an important role in B → φK1(1270) and φK1(1400) decays; (d) these pQCD predictions combined with the future precision measurements can examine the reliability of the factorization approach employed here, but also explore the complicated QCD dynamics and mixing angle θK 1 of the axial-vector K1(1270) and K1(1400) system.