初始应力的平衡
为什么我这个隧道列子初始应力不正常啊New
;园柱体外环绕放射状网格
Generate zone radcylinder size 8 25 10 5 ratio 1 1 1 1 &
p0=(36,0,0) p1=(54,0,0) p2=(36,100,0) p3=(36,0,18) &
dimension 8,8,8,8 fill group tunnel
;平行六面体外环绕放射状网格
Generate zone radtunnel size 5 25 5 5 ratio 1 1 1 1 &
p0=(36,0,0) p1=(36,0,-12) p2=(36,100,0) p3=(54,0,0) &
dimension 2,8,2,8 fill group tunnel
Generate zone reflect normal (-1,0,0) origin=(36,0,0)
Plot set rotation 20 0 30
Plot set center auto
Plot set magnificatio 1.5
Plot show
Generate zone radcylinder size 4 25 10 5 ratio 1 1 1 1 &
p0=(0,0,0) p1=(18,0,0) p2=(0,100,0) p3=(0,0,18) &
dimension 4,4,4,4 fill group tunnel
Generate zone radcylinder size 4 25 10 5 ratio 1 1 1 1 &
p0=(0,0,0) p1=(0,0,-12) p2=(0,100,0) p3=(18,0,0) &
dimension 4,4,4,4 fill group tunnel
gen zon brick p0=(18 0 18) p1=(54 0 18) p2=(18 100 18) &
p3=(18 0 50) size 10 25 5 &
ratio 1,1,1.2
gen zon brick p0=(18 0 -40) p1=(54 0 -40) p2=(18 100 -40) &
p3=(18 0 -12) size 10 25 3 &
ratio 1,1,0.8
gen zon brick p0=(0 0 18) p1=(18 0 18) p2=(0 100 18) &
p3=(0 0 50) size 5 25 5 &
ratio 1,1,1.2
gen zon brick p0=(0 0 -40) p1=(18 0 -40) p2=(0 100 -40) &
p3=(0 0 -12) size 5 25 3 &
ratio 1,1,0.8
gen zon brick p0=(54 0 0) p1=(100 0 0) p2=(54 100 0) &
p3=(54 0 18) size 5 25 5 &
ratio 1.2,1,1
gen zon brick p0=(54 0 18) p1=(100 0 18) p2=(54 100 18) &
p3=(54 0 50) size 5 25 5 &
ratio 1.2,1,1.2
gen zon brick p0=(54 0 -40) p1=(100 0 -40) p2=(54 100 -40) &
p3=(54 0 -12) size 5 25 3 &
ratio 1.2,1,0.8
gen zon brick p0=(54 0 -12) p1=(100 0 -12) p2=(54 100 -12) &
p3=(54 0 0) size 5 25 5 &
ratio 1.2,1,1
Generate zone reflect normal (-1,0,0) origin=(0,0,0)
Plot add surface yellow range z -40 25
Plot add surface red range z 25 35
Plot add surface green range z 35 45
Plot add surface blue range z 45 50
group section1 range y 0 5 group tunnel
group section2 range y 5 10 group tunnel
group section3 range y 10 15 group tunnel
group section4 range y 15 20 group tunnel
group section5 range y 20 25 group tunnel
group section6 range y 25 30 group tunnel
group section7 range y 35 40 group tunnel
group section8 range y 40 45 group tunnel
group section9 range y 45 50 group tunnel
group section10 range y 50 55 group tunnel
group section11 range y 55 60 group tunnel
group section12 range y 60 65 group tunnel
group section13 range y 65 70 group tunnel
group section14 range y 70 75 group tunnel
group section15 range y 75 80 group tunnel
;设置边界条件
fix x range x -100.1 -99.9
fix x range x 99.9 100.1
fix z range z -40.1 -39.9
fix y range y -.1 .1
fix y rang y 99.9 100.1
;参数
mo mohr
pro bulk 6.3e7 she 2.6e7 range z -40 25
pro bulk 3.3e7 she 1.3e7 range z 25 35
pro bulk 2.3e7 she 0.8e7 range z 25 35
pro bulk 1e7 she 0.2e7 range z 35 45
pro bulk 6e6 she 2e6 range z 45 50
;设置初始应力状态
set grav 0 0 -10
ini density 2000
solve
SAV 1.SAV
mo mohr
pro bulk 6.3e7 she 2.6e7 fric 70 coh 1e5 ten 1e5 range z -40 25
pro bulk 3.3e7 she 1.3e7 fric 50 coh 4e3 ten 1e3 range z 25 35
pro bulk 2.3e7 she 0.8e7 fric 30 coh 2e4 ten 8e3 range z 25 35
pro bulk 1e7 she 0.2e7 fric 25 coh 8e3 ten 5e3 range z 35 45
pro bulk 6e6 she 2e6 fric 20 coh 4e3 ten 0 range z 45 50
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
mod null range group section1
solve 网格剖分图
file:///C:/Users/hp/AppData/Local/Temp/LV3T@I6%7D3V_N3U8B03Q1 看看更改强度参数的弹塑性初始应力求解法。
现将粘聚力和抗拉强度设置为大值。
页:
[1]