受BESIII实验在三体介子$D$衰变方面最新进展的启发，我们系统地分析了$D_{(s)} \to P_1 (V \to) P_2 P_3$衰变过程，其中$V$表示矢量共振态（$\rho, K^*, \omega, \phi$），而$P_{1,2,3}$为轻赝标介子（$\pi, K, \eta^{(\prime)}$）。通过因子化辅助拓扑振幅（FAT）方法，我们计算了中间过程$D_{(s)} \to P_1 V$，并采用相对论布雷特-维格纳分布来描述随后的$V \to P_2 P_3$强衰变。通过全面包含所有相关的共振态（$\rho, K^*, \omega, \phi$），我们计算了这些衰变模式的分支比以及$D_{(s)} \to P_1 (\omega \to) KK$过程中的布雷特-维格纳尾效应。我们的框架全面结合了可因子化和不可因子化的贡献，在由于粲介子质量尺度限制而使传统方法面临挑战的非微扰区域显著提高了理论预测精度。FAT方法的结果与实验数据吻合良好，表明其在捕捉不可因子化贡献方面具有较高的精确度。我们对尚未观测到的衰变模式的预测，特别是分支比在$10^{-4}$至$10^{-3}$范围内的模式，预计将在未来的BESIII和LHCb高精度实验中得到验证。

Motivated by recent experimental advances in three-body hadronic $D$ decays from BESIII, we present a systematic analysis of $D_{(s)} \to P_1 (V \to) P_2 P_3 $ decay processes, where $V$ denotes vector resonances ($\rho, K^*, \omega$, $\phi$) and $P_{1,2,3}$ are light pseudoscalar mesons ($\pi,K, \eta^{(\prime)}$). Using the factorization-assisted topological-amplitude (FAT) approach we calculate the intermediate subprocesses $D_{(s)} \to P_1 V$, incorporating relativistic Breit-Wigner distributions to model the subsequent $V \to P_2 P_3$ strong decays. By comprehensively including all relevant resonances ($\rho, K^*, \omega, \phi$), we calculate branching fractions for these decay modes as well as the Breit-Wigner-tail effects in $D_{(s)} \to P_1 (\omega \to) KK$ processes. Our framework comprehensively incorporates both factorizable and nonfactorizable contributions, significantly improving theoretical predictions in the nonperturbative regime where conventional methods face challenges due to the limited mass scale of charm mesons. The FAT approach yields results in good agreement with experimental data, demonstrating its effectiveness in capturing nonfactorizable contributions with improved precision. Our predictions for yet-unobserved decay modes, particularly those with branching fractions in the order of $10^{-4}$ to $10^{-3}$, are expected to be tested in future high-precision experiments at BESIII and LHCb.