title of the thesis - Department of Geology - Queen's University
title of the thesis - Department of Geology - Queen's University
title of the thesis - Department of Geology - Queen's University
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
crack (838.2,737.616) (815.848,771.144)<br />
crack (815.848,771.144) (815.848,782.32)<br />
crack (815.848,782.32) (748.792,782.32)<br />
crack (748.792,782.32) (743.204,793.496)<br />
crack (743.204,793.496) (670.56,793.496)<br />
;==========================================================<br />
;===== define sub block ============================<br />
;==========================================================<br />
crack (447.04,447.04) (447.04,1117.6)<br />
crack (447.04,1117.6) (1117.6,1117.6)<br />
crack (1117.6,1117.6) (1117.6,447.04)<br />
crack (1117.6,447.04) (447.04,447.04)<br />
;==========================================================<br />
;===== define faults ==============================<br />
;==========================================================<br />
crack (0,251.92736) (1564.64,820.8772)<br />
crack (0,290.59632) (1564.64,859.54616)<br />
crack (0,355.2371429) (1564.64,1025.84504)<br />
crack (0,392.0061829) (1564.64,1062.61408)<br />
crack (0,437.4605714) (1564.64,1107.98864)<br />
crack (0,482.1629749) (1564.64,1158.05712)<br />
;===FW===<br />
crack (0,1309.295) (1564.64,307.925)<br />
;===Splays===<br />
crack (905.676,730.681) (864.557,605.485)<br />
crack (845.579,718.676) (787.063,578.045)<br />
crack (807.625,780.415) (746.192,713.585)<br />
;===1290/400E===<br />
crack (0,806.814) (670.243,806.814)<br />
crack (670.243,806.814) (922.184,789.484)<br />
crack (922.184,789.484) (1564.64,789.484)<br />
;<br />
;==========================================================<br />
;===== assign material properties =================<br />
;==========================================================<br />
PROP mat=1 dens=0.002965 k=24613 g=13856 fric=35 coh=10<br />
PROP jmat=1 jkn=700000 jks=99400 jcoh=5 jfric=35 jtens=20<br />
PROP jmat=2 jkn=700000 jks=99400 jcoh=50 jfric=50 jtens=50<br />
change mat=1 range 0,1564.64 0,1564.64<br />
change jmat=1 range 0,1564.64 0,1564.64<br />
;<br />
change jmat=2 range 446,1118, 446,448 angle -1,1<br />
change jmat=2 range 446,1118 1117,1118 angle -1,1<br />
change jmat=2 range 446,448 446,1118 angle 89,91<br />
change jmat=2 range 1117,1118 446,1118 angle 89,91<br />
;<br />
;==========================================================<br />
;===== elasto-plastic Mohr-Coulomb ================<br />
;==========================================================<br />
CHANGE cons 3 range 0,1564.64 0,1564.64<br />
;<br />
;==========================================================<br />
;===== null material to excavate ==================<br />
;==========================================================<br />
CHANGE cons 0 range 668,925 781,853<br />
CHANGE cons 0 range 812,925 726,786<br />
PAUSE<br />
;==========================================================<br />
;===== Discretize =================================<br />
;==========================================================<br />
gen edge 20.0<br />
;<br />
197