利用耗散粒子动力学 (dissipative particle dynamics, DPD) 方法对 C$_{m}$H$_{2m+1}$-PO$^{4-}$-(CH$_{2}) _{2}$-N$^{+}$(CH$_{3})_{2}$-C$_{n}$H$_{2n+1}$ (简记为 C$_{m}$-P-N-C$_{n}$, 其中 $m$, $n$=9, 9; 9, 12; 9, 15; 9, 18; 12, 12; 12, 15; 12, 18; 15, 15; 15, 18; 18, 18) 型系列两性杂双子表面活性剂在水溶液中的自组装行为开展了模拟研究, 观察到了球状、柱状、平面网格、层状、蜂窝状, 以及一维孔道、二维孔道和三维孔道等多种自组装结构的形成. 所有模拟体系均随表面活性剂浓度的逐步升高而呈现出“球状—柱状—平面网格—三维孔道—层状—二维孔道—一维孔道”的结构转化规律. 当 $m\ne n$ 时, 疏水链长度对于自组装行为有显著影响. 此外, C$_{m}$-P-N-C$_{n}$ 分子内部各功能基团与水分子间的相互作用强弱和亲近 程度可通过径向分布函数来完全体现. 研究结果可为两性杂双子表面活性剂自组装行为的进一步研究及其实际应用 提供新的见解.