In this work, we report rovibrational energy levels for four isotopologues of methylene (CH2, CHD, CD2, and 13CH2) in their ground triplet electronic state (3B1) from variational calculation up to ~10,000 cm−1 and using a new accurate ab initio potential energy surface (PES). Triplet methylene exhibits a large‐amplitude bending vibration and can reach a quasilinear configuration due to its low barrier (~2000 cm−1). To construct the ab initio PES, the Dunning's augmented correlation‐consistent core‐valence orbital basis sets were employed up to the sextuple‐ζ quality [aug‐cc‐pCVXZ, X = T, Q, 5, and 6] combined with the single‐ and double‐excitation unrestricted coupled cluster approach with a perturbative treatment of triple excitations [RHF‐UCCSD(T)]. We have shown that the accuracy of the ab initio energies is further improved by including the corrections due to the scalar relativistic effects, DBOC and high‐order electronic correlations. For the first time, all the available experimental rovibrational transitions were reproduced with errors less than 0.12 cm−1, without any empirical corrections. Unlike more “traditional” nonlinear triatomic molecules, we have shown that even the energies of the ground vibrational state (000) with rather small rotational quantum numbers are strongly affected by the very pronounced rovibrational resonance interactions. Accordingly, the polyad structure of the vibrational levels of CH2 and CD2 was analyzed and discussed. The comparison between the energy levels obtained from the effective Watson A‐reduced Hamiltonian, from the generating‐function approach and from a variational calculation was given.