Numerical analysis is important for the design, construction and maintenance of large caverns. The rock mass contains generally discontinuities such as faults, joints and fissures. The mechanical behavior and geometric characteristics of these discontinuities would have a significant impact on the stability of the caverns. In this research the Distinct Element Method(DEM) was used to analyze the structural stability of the large cavern. The Barton-Bandis Joint Model (B-B J.M) was used as a constitutive model for the joint. In addition, two different cases 1) analysis with a support system and 2) analysis with no support system, were analyzed to optimize a support system and to investigate reinforcing effects of a support system. The most significant parameters of in-situ stress, JRC of in-situ natural joints, and spatial distribution characteristics of discontinuities were acquired through field investigation. Displacement (horizontal, joint shear), maximum joint opening, maximum and minimum principal stresses, range of relaxed zone, rockbolt axial forces and shotcrete stresses were calculated at each excavation stage. As a result of analysis the calculated values proved to be under the allowable value Rockbolts also proved to be an efficient support measure to control joint shear displacement which had significant effects on extending the relaxed zone. As a consequence, the structural stability of the cavern was assured with an appropriate support system.