This paper proposes an improved disc brake system optimization method for squeal instability reduction using slip-dependent eigenvalue results. Although complex eigenvalue analysis is widely used for minimizing brake squeal instability, conventional optimization approaches still have the limitation of not being able to reflect slip rate-varying squeal instability characteristics. While relative angular velocity between the pad and disc declines due to braking, disc brake system instability gradually increases up to a specific peak velocity point and decreases until the vehicle stops, which means a maximum instability point exists during the braking process. Therefore, instability optimization should target the prevention of a maximum value during a braking scenario. The proposed optimization formulation is conducted considering maximum instability during full braking. To obtain braking time profiles, a model-based design method is employed and utilized instead of full finite element transient dynamic analysis to reduce computational cost. Kriging surrogate modeling is also used for solving the optimization problem and better express the multi-variable squeal problem. The proposed optimal design method produces minimal squeal instability during the full vehicle braking time range. The effectiveness of the proposed disc brake optimal design is demonstrated via acceleration power value comparison of the structure acceleration with that derived by conventional optimization approach.