In an iron mine where the non-sill-pillared sublevel caving system is adopted, simulations are conducted, using the three dimensional finite elementary method, with regard to pre-stoping and stoping conditions. In performing these simulations, actual stope geometry, mechanical parameters and support patterns are used as conditions for the mathematical model. Through simulation, stress distribution in drift roof, pillar sidewalls and linkage drift are systematically analyzed. Based on the stress characteristics of the stope structure, corresponding support measures are proposed, with great success in field support tests.