In this paper a computationally efficient methodology to establish the seismic fragility curve for underground tunnel structures is developed by integrating the response acceleration method and the maximum likelihood approach with a large set of artificial ground motion time-histories. It is elaborated with an illustrative numerical example.