To simulate the joint network around mining tunnels and to conduct the stability analysis of tunnels in jointed rock mass, a numerical model for jointed rock mass containing particle flow block model, original joint network and secondary joint network is established. As the orientation and area of internal secondary joints are unmeasurable, a method for calculating the secondary joint orientation and area under unloading is put forward. The unloading direction is calculated according to the change of stress...