# An FDM-Based Dynamic Zoning Method for Disturbed Rock Masses above a Longwall Mining Panel
## Appendix A  
```
loop whilep_z#null id=zone.id(p_z) w=zone.extra(p_z,11) if zone.model(p_z)#'null' then if zone.group(p_z)#50' then if zone.group(p_z)#2' then if zone.group(p_z)#1' then if zone.group(p_z)#0' then ail=zone.strain.inc.xx(p_z) ai2=zone.strain.inc.yy(p_z) ai3=zone.strain.inc.zz(p_z) ai4=zone.strain.inc.xy(p_z) ai5=zone.strain.inc.yz(p_z) ai6=zone.strain.inc.xz(p_z) dum2=zone.stress.prim.dir(p_z,z,as,ad) l=ad(l)/math.sgrdt(d)(1)*ad(d)(1)*ad(2)*ad(2)+ad(3)*ad(3) m=ad(2)/math.sgrdt(d(1)*ad(1)+ad(2)*ad(2)+ad(3)*ad(3)) n=ad(3)/math.sgrdt(ad(1)*ad(1)+ad(2)*ad(2)+ad(3)*ad(3)) el = -(! *ail+m*m*ail+2+n*n*ail3+2*l*m*ail+2*m*n*ail+2*m*n*ail5+2*l*n*ail6) if el > el_max(id) then e1_max(id) = e1 endif if zone.state(p_z,1)=0 then ee(id) = e1_max(id) endif
```  
**Algorithm A1** Algorithm Code of Dynamic Rock Parameter Adjustment. ```
if w>0 then if zone.group(p_z) # 'caving_zone' then if zone.group(p_z) # 'failure_shear' then if zone.group(p_z) # 'failure_tension' then if sigma_3<5 then \(\text{epr(id)=sigma_3*mp(w,5)+mp(w,9)}\) else if sigma_3<10 then \(\text{epr(id)=sigma_3*mp(w,6)+mp(w,10)}\) else if sigma_3<15 then \(\text{epr(id)=sigma_3*mp(w,7)+mp(w,11)}\) else if sigma_3<20 then \(\text{epr(id)=sigma_3*mp(w,8)+mp(w,12)}\) endif endif endif endif endif endif endif if zone.group(p_z)='fracture_zone_sr' then if ep(id)>10 then \(\text{zone.prop(p_z,;cohesion)=0}\) \(\text{zone.group(p_z)='fracture_zone_sr-0'}\) \(\text{if zone.pos.z(p_z)<0.034012*(zone.pos.y(p_z)-81500)-486.424}\) \(\text{zone.group(p_z)='caving_zone'}\) endif endif if zone.group(p_z)='fracture_zone_tr' then if ep(id)>1 then \(\text{zone.prop(p_z,;cohesion)=0}\) \(\text{zone.group(p_z)='fracture_zone_tr-0'}\) if zone.pos.z(p_z)<0.034012*(zone.pos.y(p_z)-81500)-486.424 \(\text{zone.group(p_z)='caving_zone'}\) endif endif endif
```
**Algorithm A1**_Cont_. ``` zone.prop(p_z,'tension')=0.1*2*zone.prop(p_z,'cohesion')*math.cos(zone.prop(p_z,'friction')*math.degrad)/(1-math.sin(zone.prop(p_z,'friction')*math.degrad)) :1/10*sigma_c zone.prop(p_z,'poisson)=(1-math.sin(zone.prop(p_z,'friction')*math.degrad))/2+0.01



## References

* (1) [PERSON]; [PERSON]; [PERSON]; [PERSON]; [PERSON] A zoning model for coal mining-induced strata movement based on microseismic monitoring. _Int. J. Rock Mech. Min._**2017**, 94, 123-138. [CrossRef]
* (2) [PERSON] _Luminal Mining_, 3 rd ed. ; CRC Press: Boca Raton, FL, USA, 2020; pp. 45-65. * (3) [PERSON]; [PERSON]; [PERSON] Evaluation of coal longwall caving characteristics using an innovative UDEC Trigon approach. _Comput. Geetech._**2014**, 55, 448-460. [CrossRef]
* (4) [PERSON]; [PERSON]; [PERSON]; [PERSON] The height of fractured water-conducting zone in undermined rock strata. _Eng. Geol._**2011**, 120, 32-39. [CrossRef]
* (5) [PERSON]; [PERSON]; [PERSON]; [PERSON] A RBFNN-based method for the prediction of the developed height of a water-conductive fractured zone for fully mechanized mining with sublevel caving. _Arab. J. Geosci._**2017**, _10_, 172. [CrossRef]
* (6) [PERSON]; [PERSON]; [PERSON]; [PERSON] Investigation of Water-Flow Fracture Zone Height in Fully Mechanized Cave Mining Beneath Thick Alluvium. _Geetech. Geol. Eng._**2017**, 35, 1745-1753. [CrossRef]
* (7) [PERSON]; [PERSON]; [PERSON]; [PERSON] Predicting the distribution of ground fissures and water-conducted fissures induced by coal mining: A case study. _SpringerPlus_**2016**, 5, 977. [CrossRef]
* (8) [PERSON]; [PERSON]; [PERSON] In situ explorations on zonal disintegration of roof strata in deep coal mines. _Int. J. Rock Mech. Min._**2012**, 49, 113-124. [CrossRef]
* (9) [PERSON]; [PERSON]; [PERSON] Mining-induced movement properties and fissure time-space evolution law in overlying strata. _Int. J. Min. Sci. Technol._**2012**, 22, 817-820. [CrossRef]
* (10) [PERSON]; [PERSON]; [PERSON]; [PERSON] Analysis of Overlying Strata Movement and Behaviors in Caving and Solid Backfilling Mixed Coal Mining. _Energies_**2017**, _10_, 1057.