To understand the bonding mechanism at the heterogeneous interface of hybrid thermoplastic composites which is difficult to obtain through experimental investigations, a series of molecular dynamic (MD) simulations were conducted in this paper.