In this paper, we systematically calculate two-body strong decays of newly observed $D_J(3000)$ and $D_{sJ}(3040)$ with 2P$(1^+)$ and 2P$(1^{+prime})$ assignments in an instantaneous approximation of the Bethe-Salpeter equation method. Our results show that both resonances can be explained as the 2P$(1^{+})$ with broad width via $^3P_1$ and $^1P_1$ mixing in $D$ and $D_s$ families. For $D_J(3000)$, the total width is 229.6 MeV in our calculation, close to the upper limit of experimental data, and the dominant decay channels are $D_2^*pi$, $D^*pi$, and $D^*(2600)pi$. For $D_{sJ}(3040)$, the total width is 157.4 MeV in our calculation, close to the lower limit of experimental data, and the dominant channels are $D^*K$ and $D^*K^*$. These results are consistent with observed channels in experiments. Given the very little information that has been obtained from experiments and the large error bars of the total decay widths, we recommend the detection of dominant channels in our calculation.