flac3d模拟锚杆与初衬支护的开挖出问题，大家帮忙呀
本帖最后由 zf5566 于 2010-5-22 22:59 编辑new
title
zhi hu
;建立坐标网格
gen zone brick size 60 5 26 p0 0 0 0 p1 60 0 0p2 0 5 0 p3 0 0 26
;选取土体本构模型
model mohr
;设定土体属性
prop bulk=2.594e7,shear=5.772e6,cohesion=4.5e4,friction=10.0 range z=22.0,26.0 ;土层1
prop bulk=4.390e7,shear=3.02e6,cohesion=5.51e7,friction=21.5 range z=18.0,22.0 ;土层2
prop bulk=2.594e7,shear=5.772e6,cohesion=4.5e4,friction=12.4 range z=14.0,18.0 ;土层3
prop bulk=4.390e7,shear=3.02e6,cohesion=5.51e7,friction=35.0 range z=0.00,14.0 ;土层4
;土体天然密度
initial density=1800 range z=22.0,26.0
initial density=1910 range z=18.0,22.0
initial density=2020 range z=14.0,18.0
initial density=2000 range z=0.0,14.0
;设定应变边界条件
fix z range z=-0.1,0.1
fix x range x=-0.1,0.1
fix x range x=59.9,60.1
fix y range y=-0.1,0.1
fix y range y=4.9,5.1
set gravity 0 0 -9.81
plot con xdisp ou on magf 20
plot set rot 30 0 30
plot set center 20 6 13
plot set mag 1.8
hist id 1 unbal
plot add hist 1
set movie avi step 10 file xin2.avi;记录
movie start
solve elas ratio=10e-5 ;分阶段塑性求解法的初始应力计算
save 011.sav
plot con szz ou on
pause;暂停计算，查看应力，应该有szz=ρgz ，
initial zd 0 yd 0 xd 0 zv 0 yv 0 xv 0
;=====================================================================
; 开挖基坑第一层6米
;=====================================================================
model null range x 0 20 z 20 26
plot con ssi ou on vel red ;x方向应力及速度矢量图
plot add hist 1
solve
save 012.sav
pause;暂停计算，查看应力和位移分布
plot con xdisp ou on magf 20
step 50
plot con sxx ou on
step 50
;=====================================================================
; 设置预应力锚杆
;=====================================================================
sel cable id 1 beg 20 2 22 end 35 2 17 nseg 10 ;自由段
sel cable id 1 beg 35 2 17 end 47 2 13 nseg 6;锚固段
sel cable id 1 prop emod 200e9 xcarea 3.8e-4 ytens 1.414e5 gr_k 1 gr_coh 1 gr_fric 1 gr_per 0.314 range cid 1 10
sel cable id 1 prop emod 200e9 xcarea 3.8e-4 ytens 1.414e5 gr_k 7e6 gr_coh 1e2 gr_fric 25 gr_per 0.314 range cid 11 16
sel liner range x 19.9 20.1 y 1 3 z 21 23 ;设置初衬
sel liner prop iso=(25e9,0.15) thick 0.1 cs_nk 8e8 cs_sk 8e8 cs_nk 0.0 cs_scoh 0.0 cs_scohres 0.0 cs_sfric 0.0
sel delete link range id 1
sel link id 100 1 target node tgt_num 20 ;建立锚杆与初衬的联系
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel cable id 1 pretension 60e3 range cid 1 10 ;在自由段施加预应力
plot block group white
plot add sel geom black red cid on scale 0.02
plot sel cable stress sk
pause;暂停计算,查看施加的预应力情况
plot add con ssi ou on vel red
solve ratio=10e-5
;=====================================================================
;开挖基坑第二层6米，共计开挖12米
;=====================================================================
model null range x 0 20 z 14 20
plo sel liner coupling stress normal
plo add sel liner coupling stress shear
pause;暂停计算
plot con ssi ou on vel green
plot add sel cable stress
plot add hist 1
solve ratio=10e-5运行到此提示错误，下面的截图
save 013.sav
pause;暂停计算，查看应力和位移分布
plot con ssi ou on vel red ;剪切应变增量及速度矢量图
step 50
plot sk dis red scale 0.07
step 50
plot sel cable stress sk
step 50
movie finish 大家帮帮我吧 弄了一晚上了 也没解决 solve后面语句去掉试试 解决了 找到问题啦 呵呵 你也不说一下是怎么解决的 让大家都学习一下 这个错误提示指的是材料赋值的错误，我检查了下锚杆材料的赋值，就找到错误了，大家以后注意就是了，是我马虎了 你的锚杆拉反了 到底是什么错误啊，楼主说的明白一些吧 楼主能不能把错误说一下呀，好让大家学习学习 是哪里错了，请版主指明，学习中，谢谢！ 应该是除零了吧
页:
[1]