模拟井筒程序

xiaogezhangwei

;模拟井筒程序
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 20  ratio 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 coh  8e3 range  group 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 20  ratio 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 20  ratio 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 coh  4e4 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 20  ratio 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 coh  6e4 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 20  ratio 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 20  ratio 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 4  ratio 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 525  range z -25 0
ini syy 0 grad 0 0 525  range z -25 0
;粘土质砂
ini szz -52500 grad 0 0 2000 range z -50 -25
ini sxx -18900 grad 0 0 718  range z -50 -25
ini syy -18900 grad 0 0 718  range z -50 -25
;砂质粘土
ini szz -102500 grad 0 0 1950 range z -75 -50
ini sxx -52788  grad 0 0 1004 range z -75 -50
ini syy -52788  grad 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 -67266  grad 0 0 867  range z -120 -100
ini syy -67266  grad 0 0 867  range 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

lakewater

结论？

xiaogezhangwei

岩石与土的参数我用的是冻结井筒的，过后才发现参数很大，井筒不会破坏

starsmoon

鼓励一下，以后多和大家交流。

wangjing0424

不错,顶顶顶顶顶

jfchai

最好把 命令流 作为附件传以下。这样方便大家的学习哦

wgy0207

谢谢楼主，我们这的井筒是斜的，建模怎么建楼主有没有好办法

lxrswf

呵呵,用ansys建呗

cumtsafe303

学习了,谢谢楼主.

bqshcumt2008

考下来好好学学！！！！

blue02

#6,正解，

yisihai

好东西,正需要

knight-hero

学习了 很有帮助 谢谢楼主

祸害

学习学习

lovemm

好东西啊，呵呵:):)