The simulation signals generated by the bearing dynamics model have a big gap with the actual signals, which limits their efficacy in bearing fault diagnosis. Therefore, it is valuable to build an accurate digital twin model of faulty rolling bearing . Firstly, a multi-degree-of-freedom bearing fault dynamics model is constructed in the virtual space for generating the vibration responses of bearing parts. Then considering that the frequency spectrum contains more characteristic information than the time-domain signal, a frequency-domain bi-directional long short-term memory (Bi-LSTM) cycle generative adversarial network (CycleGAN) named FBC-GAN is proposed to construct the frequency-domain coupling mapping relationship between the multipart vibration responses and the measured signals. In the proposed network, Bi-LSTM is used for enhancing the feature extraction ability. Meantime, a new spectrum-constraint loss is proposed to ensure the frequency-domain mapping. Next, the simulated fault bearing signals close to the actual signals are generated by FBC-GAN and Fourier transform . Finally, the results of two experiments show the superiority of the proposed method over other advanced data augmentation methods in bearing fault diagnosis with the imbalanced samples.