Poly(phosphazenes)(PZ) are versatile polymers due to their tunable properties that can be tailored for specific applications. Despite extensive experimental research, not all properties are tested, and the list of PZs studied via molecular simulations is limited. Further, a general procedure to generate and test PZ systems is lacking. We present an in situ polymerization procedure developed to make, test, and tune the thermo-mechanical properties of four PZs—poly(dichlorophosphazene)(PZ-DC), poly[bis(2,2,2-trifluoroethoxy)]phosphazene (PZ-TFE), poly(2,2,2-trifluoroethoxy-5,6-diazidohexanoxy) phosphazene (PZ-Azido), and poly(2,2,2-trifluoroethoxy-5,6-dinitratohexanoxy)phosphazene (PZ-Nitrato) via molecular dynamics simulations. The predicted thermo-mechanical properties (i.e., density and glass transition temperature) agreed with experimental values when a direct comparison of PZ systems was possible. This demonstrates the reproducibility and reliability of our procedure which will help understand the behaviour of PZs at the molecular scale.