Cold mechanical vibration of air-cooled brake disc is mainly caused by the non-uniform geometry, including disk thickness variation (DTV) and side face run-out (SFR). A mathematical model consisting of the non-uniform geometry is constructed, which is imported into the dynamic brake model established by Lagrange equation as initial parameters. The braking moment fluctuation during braking process is obtained by MATLAB software. The experimental verification is realized based on Link3900 NVH test platform, where result shows that the calculated brake pressure is quite consistent with the test value. By single variable method, the key influence factor that affected braking moment fluctuation is analyzed, including DTV, SFR, friction coefficient, contact stiffness coefficient and damping coefficient. The conclusion shows that reducing these parameters within a certain range can significantly reduce the amplitude of cold vibration, except for the contact damping coefficient.