3DEC算例1.doc

下载文档

35
0
约4.4千字
约 7页
2017-08-10 发布于重庆
举报
版权申诉
保障服务

3DEC算例1.doc

1、本文档共7页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

3DEC算例1

new poly brick 0,20 0,20 0,20 pl bl pl reset pl set dip 70 dd 210 gen edge 4 prop mat=1 dens 2500 bulk 5e9 g 3e9 phi 35 change cons 2 gravity 0 0 -10 bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range x 0 bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range x 20 bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range y 0 bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range y 20 bound stress 0,0,-5e6 0,0,0 range z 20 bound zvel=0 range z 0 insitu stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 step 500 plot block fill off stress colorby maximum reverse ret ！！！！！！！！！！！！！ new poly brick -1.5 1.5 -1.5 1.5 -1.5 1.5 pl bl color by material pl reset pl set dip 70 dd 200 tunnel radial region 1 a (-0.30 -1.5 0)(-0.3 -1.5 0.4)(-0.25 -1.5 0.47) (-0.15 -1.5 0.52)( 0 -1.5 0.55)( 0.15 -1.5 0.52) ( 0.25 -1.5 0.47)( 0.3 -1.5 0.4)( 0.3 -1.5 0) b (-0.30 1.5 0)(-0.3 1.5 0.4)(-0.25 1.5 0.47) (-0.15 1.5 0.52)( 0 1.5 0.55)( 0.15 1.5 0.52) ( 0.25 1.5 0.47)( 0.3 1.5 0.4)( 0.3 1.5 0) remove region 1 ret ！！！！！！！！！！！！ new poly brick -1.5 1.5 -1.5 1.5 -1.5 1.5 pl bl colorby material pl reset pl set dip 70 dd 200 tunnel radial region 1 a (-0.30 -1.5 0)(-0.3 -1.5 0.4)(-0.25 -1.5 0.47) (-0.15 -1.5 0.52)( 0 -1.5 0.55)( 0.15 -1.5 0.52) ( 0.25 -1.5 0.47)( 0.3 -1.5 0.4)( 0.3 -1.5 0) b (-0.30 1.5 0)(-0.3 1.5 0.4)(-0.25 1.5 0.47) (-0.15 1.5 0.52)( 0 1.5 0.55)( 0.15 1.5 0.52) ( 0.25 1.5 0.47)( 0.3 1.5 0.4)( 0.3 1.5 0) remove region 1 ret ！！！！！！！ new poly brick 0 10 0 10 0 10 plot block axes pl reset plot set dip 70 dd 150 jset dip 0 dd 180 origin 0 0 5 gen quad ndiv 4 4 4 change cons 1 prop mat=1 dens 1000 bulk 8e9 g 5e9 prop jmat=1 jkn 1e10 jks 1e10 gravity 0 0 -10 bound stress 0 0 -1e5 0 0 0 range x 0 10 y 0 10 z -.1 .1 bound xvel=0 range x -.1 .1 bound xvel=0 range x 9.9 10.0 bound y