土石坝FLAC3D模拟算例|FLAC/FLAC3D - 领先的岩土技术社区，岩土领域的媒体、社区与应用平台！

离线haobo

工程师

发帖: 202

土币: 30

威望: 580

原创币: 0

只看楼主倒序阅读使用道具楼主发表于: 2011-04-16

我建的土石坝模型，坝高20.5m，FLAC3D进行稳定性分析，命令流如下：

new
;1 几何建模============================
Generate zone brick &
p0=(0,0,0)&
p1=(55,0,0)&
p2=(0,1,0)&
p3=(0,0,20.5)&
p4=(55,1,0)&
p5=(0,1,20.5)&
p6=(3.5,0,20.5)&
p7=(3.5,1,20.5)&
    size 25 1 10 ratio 1 1 1;矩形
plot gp surface
gen zone brick &
p0 0 0 -55 p1 add 165 0 0 p2 add 0 1 0 p3 add 0 0 55 &
size 75 1 25 ratio 1 1 1 ;矩形
attach face range z -0.1 0.1 ;联接网格
;镜像网格
Generate zone reflect normal(1,0,0) origin=(0,0,0)
group soildam range z 0 20.5
group groundsill1 range z -17 0
group groundsill2 range z -55 -17

;2边界条件===========================================
fix x range x 164.9 165.1      ;固定模型右侧边界的X方向速度
fix x range x -164.9 -165.1    ;固定模型左侧边界的X方向速度
fix y range y -0.1 0.1         ;固定模型前侧边界的Y方向速度
fix y range y 0.9 1.1         ;固定模型后侧边界的Y方向速度
fix x y z range z -54.9 -55.1   ;固定模型底部边界的三个方向速度

;3指定地基土为Mohr Coulomb模型======================
model null range z 0 20.5          ;对坝体设置空单元
model moh range z -55 0           ;对地基土网格进行赋值
pro bulk 12.5e6 shear 2.7e6 coh 10e10 ten 1e10 range group groundsill1
ini density 1850 range group groundsill1
pro bulk 4.67e9 shear 2.8e9 coh 10e10 ten 1e10 range group groundsill2
ini density 2740 range group groundsill2

;4初始应力状态=======================================
set grav 0 0 -9.81
hist id=1 unbal            ;设置不平衡力监测
solve                    ;第一次求解

pro bulk 12.5e6 shear 2.7e6 coh 2.2e4 fric 20 range group groundsill1  ;赋真实参数
pro bulk 4.67e9 shear 2.8e9 coh 3.47e6 fric 32 range group groundsill2
solve                        ;第二次求解
save 地基初始应力.sav         ;保存文件

“主要问题是下边这步！！！”

;5激活坝体实际模型
ini xdisp 0 ydisp 0 zdisp 0           ;将节点位移清零
ini xvel 0 yvel 0 zvel 0              ;将节点速度清零
model moh range group soildam     ;对坝体激活实际模型
pro coh 50e3 tens 13.6e3 fric 26 dila 0 bulk 20e6 she 9.2e6  range group soildam
set grav 0 0 -9.81
ini density 2650
solve
save土石坝.sav

上述5这一步命令是“表示一次填筑坝体”吗？如果是，为什么显示的Z方向位移云图是下面的：“最大位移在坝体中部”？
一次筑坝不是位移在坝顶最大吗？
请教大家，共同探讨！谢谢！QQ：1176890449

Z方向位移云图