We have developed a microscopic model based on molecular dynamics (MD) to describe strongly interacting quark plasma. The MD simulation method has successfully been applied to a variety of many-particle systems, e.g. one-component plasmas (OCP), binary ionic mixtures and complex (dusty) plasmas. The traditional MD method makes use of point-like particles and a well defined inter-particle potential which can principally be any kind of pair-potential. Here the non-Abelian nature of the quark-quark interaction is taken into account by introducing a continuously fluctuating interaction for q-q pairs. Recent results show an interesting resonance-like behavior of the heating rate inside the strongly interacting plasma, when the applied interaction-fluctuation frequency is near the nominal (density and coupling constant dependent) plasma frequency.