To investigate the impact of plasticizer functional groups on starch plasticization, three distinct plasticizers were selected in this study: ethylene glycol (EG), ethylenediamine (EDA), and ethylenebisformamide (EBF). Three models of the plasticizer/starch system were constructed using molecular dynamics (MD) simulations, and the analysis encompassed the computation of mean square displacement (MSD), radial distribution function (RDF), and hydrogen bonding energy for each system. Additionally, the proportions of simulation were used to prepare thermoplastic starch films, which were subsequently subjected to examinations such as DSC, XRD, FT-IR, SEM, and mechanical property testing. Comparative analysis of the simulation data from the three systems and the properties of the manufactured thermoplastic starch (TPS) established that the diverse functional groups of plasticizers significantly influenced starch plasticization. In different plasticizer functional group types, it was observed that hydroxyl groups in EG and amino groups in EDA predominantly form hydrogen bonds with hydroxyl groups in starch molecular chain. In contrast, amide groups in EBF can establish hydrogen bonds not only with hydroxyl groups of starch but also with ether bonds on the starch main chain, thereby resulting in more effective starch plasticization.