A method to simulate anomalous fermion-boson interactions in Wtb vertex is presented with a minimal set of simulated samples of single top quark events at the LHC energies. In proposed method additional auxiliary vector fields corresponding to the Standard Model gauge fields are implemented for each type of the anomalous vertex structure. The method allows to simulate a manifestation of anomalous gauge couplings in two approaches used in experimental analyses either keeping only the linear order contributions in the anomalous couplings or keeping all contributions in numerators and denominators as appeared in matrix elements. All the single top quark production mechanisms, $t$-, $s$- and associative tW-channels are considered. The correctness of the proposed method is demonstrated.