共享一个隧道计算的例子

shilw

首先给出参数
;问题描述 考虑圣维南原理,取周围岩土的尺寸维隧道尺寸的5~6倍,此处取为15m
;初衬 C50管片外径3.0m，内径2.7m，C50,E＝34.5GPa，V＝0.167,
;二衬 外径2.7m，内径2.3m,C30,E＝30GPa，V＝0.167（二衬先不考虑）
;粘质粉土，4.8m，φ＝26°，重度＝19，c＝11KPa，变形模量＝10MPa
;中细纱，1.7m，φ＝28°，重度＝20，c＝0KPa，变形模量＝13Mpa

;问题描述 考虑圣维南原理,取周围岩土的尺寸维隧道尺寸的5~6倍,此处取为15m
;初衬 C50管片外径3.0m，内径2.7m，C50,E＝34.5GPa，V＝0.167,
;二衬 外径2.7m，内径2.3m,C30,E＝30GPa，V＝0.167
;粘质粉土，4.8m，φ＝26°，重度＝19，c＝11KPa，变形模量＝10MPa
;中细纱，1.7m，φ＝28°，重度＝20，c＝0KPa，变形模量＝13Mpa
;--------------------------------前处理--------------------------------------
title '沉降计算'
plot set title text 'chenjiang'
plot set rotation 20 0 30
plot set center 0 0 5
plot set dist 300
plot set mag 0.8
plot add surface green
plot add axes
plot show

;----------------------------------建立模型---------------------------------
gen zone radcyl p0 0 0 0 p1 4.7 0 0 p2 0 50 0 p3 0 0 4.7 &
size 10 30 10 4 dim 3 3 3 3 group tunnel

gen zone brick p0 4.7 0 0 p1 18 0 0 p2 4.7 50 0 p3 4.7 0 4.7 &
size 10 30 5 group brick1

gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 50 0 p3 0 0 3 &
dim 2.7 2.7 2.7 2.7 size 1 30 10 6 group segment

gen zone cylinder p0 0 0 0 p1 2.7 0 0 p2 0 50 0 p3 0 0 2.7 &
size 6 30 10 group cy

gen zone reflect ori 0 0 0 norm 0 0 -1 range z 0 4.7

gen zone brick p0 0 0 4.7 p1 4.7 0 4.7 p2 0 50 4.7 p3 0 0 9.5 &
size 5 30 5 group brick2

gen zone brick p0 4.7 0 4.7 p1 18 0 4.7 p2 4.7 50 4.7 p3 4.7 0 9.5 &
size 10 30 5 group brick3

gen zone brick p0 0 0 -18 p1 4.7 0 -18 p2 0 50 -18 p3 0 0 -4.7 &
size 5 30 15 group brick4

gen zone brick p0 4.7 0 -18 p1 18 0 -18 p2 4.7 50 -18 p3 4.7 0 -4.7 &
size 10 30 15 group brick5

gen zon reflect norm 1 0 0 orig 0 0 0


;------------------------赋土参数-------------------------------------------
model mohr
;E=13Mpa,v=0.3
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 1e10 tension 1e10 range z -18 4.7
;E=10Mpa,v=0.25
prop density 1900 bulk 10e6 shear 4e6 fric 26 coh 1e10 tension 1e10 range z 4.7 9.5

;---------------------------初始地应力--------------------------------------
ini xdisp 0 ydisp 0 zdisp 0
ini szz -1.85e5 grad 0 0 1.9e4 range z 4.7 9.5             ;19×4.8+20×4.7＝185.2KPa
ini sxx -6.67e4 grad 0 0 6.84e3 range z 4.7 9.5          ;侧压力系数取为0.36 （66.67KPa）
ini syy -6.67e4 grad 0 0 6.84e3 range z 4.7 9.5
ini szz -1.85e5 grad 0 0 2.0e4 range z -18 4.7
ini sxx -6.67e4 grad 0 0 7.2e3 range z -18 4.7
ini syy -6.67e4 grad 0 0 7.2e3 range z -18 4.7

;-------------------------------边界条件-------------------------------------
apply nstress -1e4 range z 9.4 9.6
fix x range x -18.1 -17.9
fix x range x 17.9 18.1
fix y range y 49.9 50.1
fix y range y -0.1 0.1
fix z range z -18.1 -17.9
set grav 0 0 -10 small
;solve
;save small.sav

;-------------------------------开挖 求解-----------------------------------


model null range group cy

;---------------------------
ini xdisp 0 ydisp 0 zdisp 0
;---------------------------
model elas range group segment
prop density 2500 bulk 34.5e9 shear 14.78e9 range group segment

;E=13Mpa,v=0.3
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 0 range group brick1
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 0 range group tunnel
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 0 range group brick4
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 0 range group brick5

;E=10Mpa,v=0.25
prop density 2000 bulk 13e6 shear 5e6 fric 28 coh 0 range z 4.7 9.5




hist unbal
hist gp xdis 3 0 0
hist gp xdis -3 0 0
hist gp zdis 0 0 9.5
hist gp zdis 0 0 7.5
hist gp zdis 0 0 5.5
hist gp zdis 0 0 4.5
hist gp zdis 0 0 3
hist gp zdis 0 0 -3
set large
;solve
;save n.sav

大虾指正，多谢。
同时共享代码～

[ 本帖最后由 shilw 于 2006-9-10 16:32 编辑 ]

chunhuo420

请打包，否则没法复制。

wyb8254

对,打个包

shilw

已经修改了。高手斧正，多谢多谢

笑沧海

为何两侧沉降大，中间沉降小？开挖不会隆起？

zhlouv

你开始虽然写了求解初始平衡的语句，但是加了“  ；”  所以语句没有执行，这样即使你定义了初始应力，但是也会出现两侧沉降大，中间沉降小的错误情况

yxf4171

哈哈，很好的贴子，支持…………

gusiny

没法复制啊？？

hxbhwj

对,打个包

liuyongseu

对,打个包

liyun-hust

恩
没法复制呢

FredZhangFeng

是啊  复制下来的是乱码

dongwuu

实在太感谢了!

baby20023402

爽！正在找！

lq0719

多谢了，正在找~~~~~~~~~~~

天天见

顶一下

looloozheng

沉降云图有点问题好像，沉降规律应该是中间大，两侧小

boarer

请教一个问题，我挖的洞子怎么一计算沉降的时候就看不到洞子了啊？

mazhongyuan9049

楼主 想问你是怎摸考虑开挖后应力释放问题的?

csuroylin

因为在上表面增加了均布荷载，另外土体比较软，洞子参数大，因此形成两边沉降大，中间沉降小的趋势。如果没有均布荷载，就应该是中间大，两边小的趋势。
同一实验室的老师，看了一眼，就作出的判断。我看了一下命令流，果然有均布荷载，佩服。