模拟井筒程序
;模拟井筒程序new
;第一层(砂土)
gen zone radcyl p0 0 0 0 p1 50 0 0 p2 0 0 -25 p3 0 50 0 &
p8 5 0 0 p9 0 5 0 p10 5 0 -25 p11 0 5 -25 &
size 5 5 10 20ratio 1 1 1 1.1
gen zone reflect dd 270 dip 90 range z -25 0
group z1 range z -25 0
mo mo range group z1
prop bulk 1e8 shear 2e7 fric 38 coh8e3 rangegroup z1
;第二层(粘土质砂)
gen zone radcyl p0 0 0 -25 p1 50 0 -25 p2 0 0 -50 p3 0 50 -25 &
p8 5 0 -25 p9 0 5 -25 p10 5 0 -50 p11 0 5 -50 &
size 5 5 10 20ratio 1 1 1 1.1
gen zone reflect dd 270 dip 90 range z -50 -25
group z2 range z -50 -25
mo mo range group z2
prop bulk 6e7 shear 1.4e7 fric 30 coh 3e4 range group z2
;第三层(砂质粘土)
gen zone radcyl p0 0 0 -50 p1 50 0 -50 p2 0 0 -75 p3 0 50 -50 &
p8 5 0 -50 p9 0 5 -50 p10 5 0 -75 p11 0 5 -75 &
size 5 5 10 20ratio 1 1 1 1.1
gen zone reflect dd 270 dip 90 range z -75 -50
group z3 range z -75 -50
mo mo range group z3
prop bulk 4e7 shear 1e7 fric 22 coh4e4 range group z3
;第四层(粘土)
gen zone radcyl p0 0 0 -75 p1 50 0 -75 p2 0 0 -100 p3 0 50 -75 &
p8 5 0 -75 p9 0 5 -75 p10 5 0 -100 p11 0 5 -100 &
size 5 5 10 20ratio 1 1 1 1.1
gen zone reflect dd 270 dip 90 range z -100 -75
group z4 range z -100 -75
mo mo range group z4
prop bulk 2e7 shear 6e6 fric 15 coh6e4 range group z4
;基岩段(井筒共深120m)
gen zone radcyl p0 0 0 -100 p1 50 0 -100 p2 0 0 -120 p3 0 50 -100 &
p8 5 0 -100 p9 0 5 -100 p10 5 0 -120 p11 0 5 -120 &
size 5 4 10 20ratio 1 1 1 1.1
;井筒支撑(取30m的基岩)
gen zone radcyl p0 0 0 -120 p1 50 0 -120 p2 0 0 -150 p3 0 50 -120 &
p8 5 0 -120 p9 0 5 -120 p10 5 0 -150 p11 0 5 -150 &
size 5 6 10 20ratio 1 1 1 1.1
gen zone cshell p0 0 0 -120 p1 5 0 -120 p2 0 0 -150 p3 0 5 -120 &
p8 4 0 -120 p9 0 4 -120 p10 4 0 -150 p11 0 4 -150 &
size 1 6 10 4ratio 1 1 1 1 fill
gen zone reflect dd 270 dip 90 range z -150 -100
group z5 range z -150 -100
mo mo range group z5
prop bulk 3.2e9 shear 1.6e10 fric 43 coh 4e6 range group z5
;接触面的定义(在井壁外部及井端部)
interface 1 face range cylinder end1 0 0 0 end2 0 0 -25 radius 5.01
interface 2 face range cylinder end1 0 0 -25 end2 0 0 -50 radius 5.01
interface 3 face range cylinder end1 0 0 -50 end2 0 0 -75 radius 5.01
interface 4 face range cylinder end1 0 0 -75 end2 0 0 -100 radius 5.01
interface 5 face range cylinder end1 0 0 -100 end2 0 0 -120.1 radius 5.01 &
cylinder end1 0 0 -100 end2 0 0 -120.1 radius 4.99
interface 6 face range cylinder end1 0 0 -119.9 end2 0 0 -120.1 radius 5.01
;界面参数(不很确定)
interface 1 prop kn 3e8 ks 1e9 fric 20 coh 30000
interface 2 prop kn 1e8 ks 1e9 fric 20 coh 30000
interface 3 prop kn 5e7 ks 1e9 fric 20 coh 30000
interface 4 prop kn 6e7 ks 1e9 fric 20 coh 30000
interface 5 prop kn 4e8 ks 1e9 fric 20 coh 30000
interface 6 prop kn 4e8 ks 1e9 fric 20 coh 30000
;井壁(厚1m)
gen zone cshell p0 0 0 121 p1 5 0 121 p2 0 0 1 p3 0 5 121 &
p8 4 0 121 p9 0 4 121 p10 4 0 1 p11 0 4 1&
size 1 24 10 ratio 1 1 1
gen zone reflect dip 90 dd 270 range z 121 1
group jingbi range z 121 1
ini z add -121 range group jingbi
mo mo range group jingbi
prop bulk 1.4e10 shear 4.6e12 fric 38 coh 3.2e6 range group jingbi
plot surf
plot add interface red
save jingtong1.sav
;初始参数
set grav 0 0 -10
ini density 2500 range group jingbi
;围岩应力条件
;砂土层
ini szz 0 grad 0 0 2100 range z -25 0
ini sxx 0 grad 0 0 525range z -25 0
ini syy 0 grad 0 0 525range z -25 0
;粘土质砂
ini szz -52500 grad 0 0 2000 range z -50 -25
ini sxx -18900 grad 0 0 718range z -50 -25
ini syy -18900 grad 0 0 718range z -50 -25
;砂质粘土
ini szz -102500 grad 0 0 1950 range z -75 -50
ini sxx -52788grad 0 0 1004 range z -75 -50
ini syy -52788grad 0 0 1004 range z -75 -50
;粘土
ini szz -150000 grad 0 0 1900 range z -100 -75
ini sxx -100050 grad 0 0 1267 range z -100 -75
ini syy -100050 grad 0 0 1267 range z -100 -75
;基岩
ini szz -202000 grad 0 0 2600 range z -120 -100
ini sxx -67266grad 0 0 867range z -120 -100
ini syy -67266grad 0 0 867range z -120 -100
;边界条件(井筒轴对称,沿xoz面固定,下卧基岩视为稳定没有位移)
fix x range x 49.9 50.1
fix x range x -49.9 -50.1
fix y range y 49.9 50.1
fix y range y -0.1 .1
fix z range z -119.9 -120.1
hist unbal
solve ratio 5e-3
plot surf
plot add interface red
save jingtong2.sav
Re:模拟井筒程序
结论?Re:模拟井筒程序
岩石与土的参数我用的是冻结井筒的,过后才发现参数很大,井筒不会破坏Re:模拟井筒程序
鼓励一下,以后多和大家交流。 不错,顶顶顶顶顶 最好把 命令流 作为附件传以下。这样方便大家的学习哦 谢谢楼主,我们这的井筒是斜的,建模怎么建楼主有没有好办法 呵呵,用ansys建呗 学习了,谢谢楼主. 考下来好好学学!!!! #6,正解,:lol 好东西,正需要 学习了 很有帮助 谢谢楼主 学习学习 好东西啊,呵呵:):):)
页:
[1]