A pilot plant of underground cold storage for food has been excavated as a R&D program. For the stability assessment of underground cold storage opeinengs in shallow jointed rocks, three kinds of stability problems were analyzed by numerical methods. For the analysis of unstability by rock block movements, DEM was used considering the statistical distribution of rock joints. Concerning thermally induced cracking, FDM was used with thermomechanical stress analysis. Finally, in order to evaluate the joint failure during the thawing process, BE algorithm was applied. Numerical examples applied for the pilot plant show that the possibility of unstable failure of opeings exists but can be avoided with proper rock reinforcements provided.