27.10.2014 Views

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

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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=0 jfric=35 jtens=20<br />

PROP jmat=3 jkn=700000 jks=99400 jcoh=0 jfric=20 jtens=20<br />

PROP jmat=4 jkn=700000 jks=99400 jcoh=50 jfric=50 jtens=50<br />

change mat=1 range 0,1564.64 0,1564.64<br />

;<br />

change jmat=1 range 0,1564.64 0,1564.64<br />

;<br />

change jmat=3 range 0,1564.64 246,1171 angle 15,30<br />

;<br />

change jmat=2 range -1,1565 303,1324 angle 140,160<br />

;<br />

change jmat=4 range 446,1118, 446,448 angle -1,1<br />

change jmat=4 range 446,1118 1117,1118 angle -1,1<br />

change jmat=4 range 446,448 446,1118 angle 89,91<br />

change jmat=4 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 />

PAUSE<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 />

191

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!