In this study, a stochastic finite element model is proposed with a view to consider the uncertainty of physical properties of discontinuous rock mass in the analysis of structural behavior on underground caverns. In so doing, the LHS(Latin Hypercube sampling) technique has been applied to make up weak points of the Crude Monte Carlo technique. Concerning the effect of discontinuities, a joint finite element model is used that is known to be superior in explaining faults, cleavage, things of that nature. To reflect the uncertainty of material properties, the variables such as the the elastic modulus, the poisson's ratio, the joint shear stiffness, and the joint normal stiffness have been used, all of which can be applicable through normal distribution, log-normal distribution, and rectangulary uniform distribution. The validity of the newly developed computer program has been confirmed in terms of verification examples. And, the applicability of the program has been tested in terms of the analysis of the circular cavern in discontinuous rock mass.