Here, we present PySoftK, a modular and versatile code to model polymer structures with different topologies. We present the various modules that currently exist within PySoftK to build different polymer structures, and show how they can be combined in order to build highly complex polymer topologies in an automated way. We subsequently demonstrate the various functionalities that are found within the code to facilitate high-throughput molecular modelling calculations. PySoftK can also be used in the parameterization of dihedral terms in conjugated polymers, which are commonly poorly defined by standard classical forcefields. Finally, we briefly review the steps that have been taken when developing the code to ensure it can be successfully used in a broad range of applications.