超长超深隧道施工及泄水沉降模拟

T-MACLeo

模拟的具体情况：四层水平土层共计810m深，长7.3Km（由于要看差异沉降所以建3D模型并且考虑隧道的周围土体不同引起的差异沉降不同，所以希望建全长的情况）、纵剖面斜度6度、内径7m，厚0.4m、每环宽1.5m的超长隧道，由于埋深大导致水压力、隧底每隔5环设置一个泄水孔；求开挖引起的沉降，隧顶应力变化曲线，浅埋、中埋、深埋典型断面以及整个隧道的应力云图，再求泄水引起的固结沉降
问题：
1.假设一环管片一个单元大小的话需要的单元格数=7300*810*20（对称模型宽度）=118260000，需要1千万有多的单元，FLAC3D中超限吗？
2.水平地层，倾斜隧道的模型能直接在FLAC3D中建出来吗？
3.泄水孔怎么模拟？

T-MACLeo

对了，感兴趣的童鞋，可以加我建的QQ群：102289697（渗流讨论），进群请冒泡，大家一起分享资料、程序和经验，一起解决难题。。。先谢过了。。。

T-MACLeo

2# T-MACLeo

T-MACLeo

3# T-MACLeo
网速太慢了哈。。。回了条空白的！！！！不好意思
还有一个重要问题：有没有谁知道一个叫STEAD(Subway Tunneling-induced ground-Enviroment-damage Assessment and conrtol Design system)基于 FLAC3D的二次开发计算软件，如果有的话可否发我一份？邮箱：liyanshenxian@sohu.com，先谢一个

T-MACLeo

补充问题。。。我想建得太多了也不好算，就算长度为30m的（20环管片）
首先，模型太高了（不能简化为力加上去，因为要考虑潜水的泄水问题，当然，如果仍可简化的话，希望大家能给意见），这样模型就细长了，可以改变它的显示比例吗？比如纵向缩小10倍。
其次，也是最关键的，模型里面建出来的斜隧道很不成形啊。。。轴歪了两个半圆就跟着歪，最后不成隧道
我会把问题的图片发到我的空间，这里把命令流给大家，大家运行一下就知道我的问题的意思了
n
config fluid
set fl off

;建立几何模型
plo blo gro
gen zon radcyl p0 0 0 60 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 20 p4 a 20 30 3.153 p5 a 0 30 20 p6 a 20 0 20 p7 a 20 30 20 dim  3.5 3.5 3.5 3.5 size 1 20 8 16 ratio 1 1 1 1.2 group 56层 fill group 隧道
pau
gen zon radcyl p0 0 0 60 p1 a 0 0 -20 p2 a 0 30 3.153 p3 a  20 0 0 p4 a 0 30 -20 p5 a 20 30 3.15 p6 a 20 0 -20 p7 a 20 30 -20 dim  3.5 3.5 3.5 3.5 size 1 20 8 16 ratio 1 1 1 1.2 group 56层 fill group 隧道
pau

T-MACLeo

昨天没上上网，所以没更新得这么快，有朋友建议用ANAsys导入，可我不会ANAsys啊!!!,高手给个招儿，可以直接在FLAC里面建么？即使建成了，在加Liner衬砌单元的时候会不会也遇到同样的问题呢？

T-MACLeo

倾斜地层的问题解决了，可初始平衡就出错咯，我这里的网貌似只有早上才是好的！！现把命令流奉上：这命令流最后是不能算的，出来警告“zero volume tet in 3249”大家有没有知道是什么问题的，献计献策啊！！！谢谢咯。。。

n
config fluid
set fl off

;建立几何模型
plo blo gro
gen zon cshell p0 0 0 60 p1 a 4 0 0 p2 a 0 30 3.153 p3 a 0 0 4 p4 a 4 30 3.153 p5 a 0 30 7.153 p8 a 3.5 0 0 p9 a 0 0 3.5 p10 a 3.5 30 3.153 p11 a 0 30 6.653 dim 3.5 3.5 3.5 3.5 size 1 20 8 1 group 注浆层 fill group 隧洞
gen zon radcyl p0 0 0 60 p1 a 20 0 0 p2 a 0 30 3.153 p3 a 0 0 20 p4 a 20 30 3.153 p5 a 0 30 20 p6 a 20 0 20 p7 a 20 30 20 p8 a 4 0 0 p9 a 0 0 4 p10 a 4 30 3.153 p11 a 0 30 7.153 dim 4 4 4 4 size 1 20 8 16 ratio 1 1 1 1.2 group 56层
gen zon cshell p0 0 0 60 p1 a 0 0 -4 p2 a 0 30 3.153 p3 a 4 0 0 p4 a 0 30 -0.847 p5 a 4 30 3.15 p8 a 0 0 -3.5 p9 a 3.5 0 0 p10 a 0 30 -0.347 p11 a 3.5 30 3.153 dim  3.5 3.5 3.5 3.5 size 1 20 8 1 ratio 1 1 1 1.2 group 注浆层 fill group 隧洞
gen zon radcyl p0 0 0 60 p1 a 0 0 -20 p2 a 0 30 3.153 p3 a  20 0 0 p4 a 0 30 -20 p5 a 20 30 3.135 p6 a 20 0 -20 p7 a 20 30 -20 p8 a 0 0 -4 p9 a 4 0 0 p10 a 0 30 -0.847 p11 a 4 30 3.153 dim 4 4 4 4 size 1 20 8 16 ratio 1 1 1 1.2 group 56层
gen zon bri p0 0 0 0 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 40 size 20 20 4 ratio 1 1 0.6 group 56层
gen zon bri p0 0 0 80 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 80 size 20 20 6 ratio 1 1 1.6 group 56层
gen zon bri p0 0 0 160 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 194 size 20 20 1 group 34层
gen zon bri p0 0 0 354 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 441 size 20 20 1 group 2层
gen zon bri p0 0 0 795 p1 a 20 0 0 p2 a 0 30 0 p3 a 0 0 15 size 20 20 1 group 1层
;镜像处理（考虑到水的边界关系）
gen zon ref ori 0 0 0 norm -1 0 0
;由于网格划分不同使用面连接
attach face
;加位移边界
;plo gpfix red sk 显示约束
fix x range x -20.1 -19.9
fix x  range x 19.9 20.1
fix y range y 29.9 30.1
fix y range y -0.1 0.1
fix z range z -0.1 0.1
;为模型赋本构参数
;设置Mohr_coulomb模型
model mohr
def derive
s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
s_mod2=E_mod2/(2.0*(1.0+p_ratio2))
b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))
s_mod3=E_mod3/(2.0*(1.0+p_ratio3))
b_mod3=E_mod3/(3.0*(1.0-2.0*p_ratio3))
s_mod4=E_mod4/(2.0*(1.0+p_ratio4))
b_mod4=E_mod4/(3.0*(1.0-2.0*p_ratio4))
end
set E_mod1=22e6 p_ratio1=0.3 E_mod2=26.3e6 p_ratio2=0.22 E_mod3=24e6  p_ratio3=0.3 E_mod4=40e6 p_ratio4=0.21
derive
prop bulk b_mod1 shear s_mod1 cohe 0 tens 1e10 fric 35 ran group  1层
prop bulk b_mod2 shear s_mod2 cohe 27.2e6 tens 1e10 fric 27.8 ran group 2层
prop bulk b_mod3 shear s_mod3 cohe 5e6 tens 1e10 fric 35 ran group 34层
prop bulk b_mod4 shear s_mod4 cohe 40e6 tens 1e10 fric 34.7 ran group 56层
prop bulk b_mod4 shear s_mod4 cohe 40e6 tens 1e10 fric 34.7 ran group 注浆层
prop bulk b_mod4 shear s_mod4 cohe 40e6 tens 1e10 fric 34.7 ran group 隧洞
;设置初始密度
ini dens=1400 range group 1层
ini dens=1700 range group 2层
ini dens=1600 range group 34层
ini dens=1700 range group 56层
ini dens=1700 range group 注浆层
ini dens=1700 range group 隧洞
;设置流体模型
model fl_iso ;各向同性渗流模型
prop perm 9.26e-9 porosity 0.5 ran group 1层
prop perm 8.1e-9 porosity 0.5 ran group 2层
prop perm 10.42e-9 porosity 0.5 ran group 34层
prop perm 8.1e-9 porosity 0.5 ran group 56层
prop perm 8.1e-9 porosity 0.5 ran group 注浆层
prop perm 8.1e-9 porosity 0.5 ran group 隧洞
ini fdens 1000 ;流体密度
ini fmod 2e9 ;流体模量
ini sat 1.0 ran z 0 809 ;流体饱和度
set grav 0 0 -10 ;必需在water table前设置
water dens 1000 ;在water table 命令之前用于计算孔压
water table ori 0 0 809 norm 0 0 1 ;设置水位面
;(ini pp 809e4 grad 0 0 -10000 不用水位面时用此命令设置初始水压）
;初始土压平衡
ini szz -0.285e6 grad 0 0 19e3 ran z 795 810
ini szz -9.987e6 grad 0 0 22e3 ran z 354 795
ini szz -14.082e6 grad 0 0 21e3 ran z 160 354
ini szz -17.062e6 grad 0 0 22e3 ran z 0 160
ini syy -0.1425e6 grad 0 0 19e3 ran z 795 810
ini syy -4.9935e6 grad 0 0 22e3 ran z 354 795
ini syy -7.041e6 grad 0 0 21e3 ran z 160 354
ini syy -8.531e6 grad 0 0 22e3 ran z 0 160
ini sxx -0.1425e6 grad 0 0 19e3 ran z 795 810
ini sxx -4.9935e6 grad 0 0 22e3 ran z 354 795
ini sxx -7.041e6 grad 0 0 21e3 ran z 160 354
ini sxx -8.531e6 grad 0 0 22e3 ran z 0 160
set large ;开启大位移模式
;记录最大不平衡力
hist unbal
sav 初始平衡前模型.sav
set mech on fl off
solve
sav 初始平衡.sav
ini xdis=0 ydis=0 zdis=0                               ;位移归零
sav 初始平衡归零.sav

T-MACLeo

上一个不能算的（已将长度减小，已解决倾斜隧道问题），出来警告“zero volume tet in 3249”的问题解决：原因是埋深太大，E值太小，加了三个数量级就可以了
现在的问题是：如何去设置一个泄水孔？

T-MACLeo

还有一个问题，昨天在编程的时候用到一个CASEOF的表达式，你们看看这问题出在哪儿？提示说CASEOF 是不存在的编码
CASEOF n
CASE 1
save 第一环开挖完成.sav
CASE 2
save 第二环开挖完成.sav
ENDCASE
我又试了，IF也不行，是不是我的是绿色版的没有这些功能的原因呢？