This paper establishes a general dynamic two-dimensional nonlinear time-varying (TDNLTV) vibration model with multi-drive shaft for shield tunnel boring machine (TBM) cutterhead driving system. Multipleinput and multiple-output (MIMO) state-space vibration model of cutterhead driving system is also presented. TDNLTV vibration model is simulated, and physical parameters' effects on the vibration response are detailed investigated. A preliminary approach is proposed to reduce vibration or nonlinear vibration behaviors. Dynamic TDNLTV vibration model provides a basic foundation of fault detection and diagnosis (FDD) and health monitoring for cutterhead driving system. i-th induction motor rotor viscous damped after equivalent coupling-i J z,i inertia of the i-th coupling between motor-i and reducer-i b z,i viscous damped of the i-th coupling between motor-i and reducer-i J r,ii -th induction motor rotor inertia b r,i i-th induction motor rotor viscous damped J m large gear inertia b m large gear viscous damped r m large gear radius θ p,i angular displacement of i-th active pinion θ i angular displacement of i-th motor rotor θ m large gear angular displacement ω m large gear angular speed J c,