In this article, in order to investigate the nonlinear forced vibrations in multi-physics fields, a coupled nonlinear modeling for composite cylindrical shells is developed by the improved Donnell nonlinear shell theory and Maxwell static electricity/magnetism equations. The nonlinear governing equations are obtained by Hamilton principle, and then the partial differential equations are transformed into sixty-degree-of-freedom nonlinear ordinary differential equations through multi-mode Galerkin technique. Afterwards, numerical simulations are conducted by the pseudo-arclength continuation algorithm. The comparisons and validations show that the present procedure is capable of solving coupled multi-modes nonlinear partial differential equations. Finally, the effects of the external temperature change, magnetic potential, electric potential, excitation amplitude on nonlinear vibration response of composite cylindrical shells are evaluated.