The need to simulate complex and multi-physics systems is increasingly common in the industry, especially with the advent of mechatronic design. This coupled with economic pressure pushing the industry towards optimizing their products led to a strong increase in simulation time that technological advances can not compensate. An hybrid method analytical/finite element has been developed for the time domain simulation of gear transmissions involving contact non-linearities. Coupled with a reduction of finite element models, this approach allows fast resolution of the dynamics of a gearbox. Consequently, it is possible to integrate a gearbox as a part of a more complex mechatronic system. All physical phenomena involved in such a complex product are solved at the same time allowing strong coupling to be considered. The optimization of such a system is discussed with an example to conclude this study.