In this work, we study the S, P, and D wave resonance contributions to three-body decays B ðsÞ → ½D ðÃÞ ;D ðÃÞ K þ K − by employing the perturbative QCD (PQCD) approach, where the kaon-kaon invariant mass spectra are dominated by the f 0 ð980Þ, f 0 ð1370Þ, ϕð1020Þ, ϕð1680Þ, f 2 ð1270Þ, f 0 2 ð1525Þ, f 2 ð1750Þ, and f 2 ð1950Þ resonances. The KK S-wave component f 0 ð980Þ is modeled with the Flatté formalism, while other resonances are described by the relativistic Breit-Wigner line shape. The corresponding decay channels are studied by constructing the kaon-kaon distribution amplitude Φ KK , which captures important final state interactions between the kaon pair in the resonant region. We found that the PQCD predictions for the branching ratios for most considered decays agree with currently available data within errors. The associated polarization fractions of those vector-vector and vector-tensor decay modes are also predicted, which are expected to be tested in the near future experiments. The invariant mass spectra for the corresponding resonances in the B ðsÞ → ½D ðÃÞ ;D ðÃÞ K þ K − decays are well established, which can be confronted with the precise data from the LHCb and Belle II experiments.