This study describes the methodology employed to construct a seismic fragility function based on a pre-existing numerical model tailored for underground stations. Employing a dynamic numerical model, a comprehensive analysis encompassing 110 distinct cases was conducted, each varying in soil depth and classification. Seismic waves, conforming to the standard design spectrum, were utilized within these numerical analyses. The formulation of the fragility function within the constructed model follows a structured approach, segmented by damage indices and severity levels. This systematic breakdown serves to outline the fundamental framework for establishing the fragility function, providing insights into its development process. Subsequently, the derived fragility function underwent a rigorous comparative analysis against established seismic fragility functions from prior studies. This comparative assessment serves as a critical evaluation tool, allowing for an appraisal of the suitability and robustness of the newly developed fragility function in relation to existing benchmarks.