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.
;==========================================================<br />
PROP mat=1 dens=0.002965 k=24613 g=13856 fric=35 coh=10<br />
PROP jmat=1 jkn=700000 jks=99400 jcoh=0 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 1 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 />
;==========================================================<br />
;===== Discretize =================================<br />
;==========================================================<br />
gen edge 20.0<br />
;<br />
;==========================================================<br />
;===== set initial stress conditions ==============<br />
;==========================================================<br />
bound stress -102.70,0,0 range -1,1 0,1564.64<br />
bound stress -102.70,0,0 range 1563.64,1565.64 0,1564.64<br />
bound stress 0,0,-73.43 range 0,1564.64 -1,1<br />
bound stress 0,0,-73.43 range 0,1564.64 1563.64,1565.64<br />
insitu stress -102.70, 0, -73.43<br />
INITIAL sxx -102.70<br />
INITIAL syy -73.43<br />
INITIAL szz -83.43<br />
;==========================================================<br />
;===== Set boundary velocity to zero ==============<br />
;==========================================================<br />
BOUNDARY yvel=0.0 range 0,1564.64 -1,1<br />
BOUNDARY yvel=0.0 range 0,1564.64 1563,1565<br />
BOUNDARY xvel=0.0 range -1,1 0,1564.64<br />
BOUNDARY xvel=0.0 range 1563,1565 0,1564.64<br />
;==========================================================<br />
;===== time steps =================================<br />
;==========================================================<br />
hist unbal<br />
hist xvel(50,50)<br />
hist yvel(50,50)<br />
hist xdisp(50,50)<br />
hist ydisp(50,50)<br />
SOLVE step 1000<br />
set plot po color<br />
;<br />
def ratios3s1<br />
bi=block_head<br />
loop while bi#0<br />
zi=b_zone(bi)<br />
loop while zi#0<br />
xval= abs(z_sxx(zi))<br />
204