This paper pays attention to magnetic flux linkage optimization of a direct-driven surface-mounted permanent magnet synchronous generator (D-SPMSG). A new compact representation of the D-SPMSG nonlinear dynamic differential equations to reduce system parameters is established. Furthermore, the nonlinear dynamic characteristics of new D-SPMSG equations in the process of varying magnetic flux linkage are considered, which are illustrated by Lyapunov exponent spectrums, phase orbits, Poincaré maps, time waveforms and bifurcation diagrams, and the magnetic flux linkage stable region of D-SPMSG is acquired concurrently. Based on the above modeling and analyses, a novel method of magnetic flux linkage optimization is presented. In addition, a 2 MW D-SPMSG 2D/3D model is designed by ANSYS software according to the practical design requirements. Finally, five cases of D-SPMSG models with different magnetic flux linkages are simulated by using the finite element analysis (FEA) method. The nephograms of magnetic flux density are agreement with theoretical analysis, which both confirm the correctness and effectiveness of the proposed approach.