A probabilistic framework based on stochastic simulation is presented in this chapter for optimal design of supplemental dampers for multi - span bridge systems supported on abutments and intermediate piers through isolation bearings. The bridge model explicitly addresses nonlinear characteristics of the isolators and the dampers, the dynamic behavior of the abutments, and the effect of pounding between the neighboring spans to each other as well as to the abutments. A probabilistic framework is used to address the various sources of structural and excitation uncertainties and characterize the seismic risk for the bridge. Stochastic simulation is utilized for evaluating this seismic risk and performing the associated optimization when selecting the most favorable damper characteristics. An illustrative example is presented that considers the design of nonlinear viscous dampers for protection of a two-span bridge.