System restoration models are usually adopted for seismic resilience analysis of bridges. However, no analytical process can be found in the existing literature to develop the restoration models for bridges. This paper proposed a new Monte Carlo-based method to derive the generalized system restoration models for ductility highway bridges. In the proposed method, a large number of random samples for ductility bridges were generated by considering uncertainty of structural parameters. The IDA was adopted for producing the damaged bridge samples at different damage states. The repair time of each bridge sample was estimated to generate the actual restoration curve. The Monte Carlo simulations were then adopted to estimate the expected mean functionality curves, which were used to derive the generalized system restoration models by using mathematical functions. Finally, seismic resilience analysis based on the derived generalized system restoration models was conducted and compared with the traditional method to illustrate the effectiveness of the proposed method. It is concluded that the proposed Monte Carlo-based method is an efficient and reliable method for developing the restoration models for ductility highway bridges.