发布时间 : 星期二 文章FLAC3D命令流更新完毕开始阅读
plo cre pile plo current pile
plo set back black fore white mag 0.8 plo add sel geo id on sca .04 magf 1e3 plo add sel fapp lgreen magf 1e3
plo add sel pile mom my lblue lred magf 1e3 axe yel set movie avi step 100 file pile.avi movie start sol age 1 movie finish 11-1 new
conf dyn ;打开动力计算功能 gen zone brick size 10 5 10 mod elas
mod null range x=0,5 z=5,10 ;删除部分网格 fix z range x=-.1 .1 z=.1 10.1 ;设置静力边界条件 fix z range x=9.9,10.1 z=.1 10.1 fix y range y=-.1 .1 fix y range y=4.9 5.1
prop bulk 2e8 shear 1e8 ;设置土体参数 prop bulk 4e9 shear 2e9 range x=5,6 z=5,10 ;设置墙体参数（土体参数的20倍）
ini dens 2000 ;设置密度
def setup ;动荷载中的变量赋值 freq = 1.0
omega = 2.0 * pi * freq old_time = clock end
setup ;执行变量赋值 def wave ;定义动荷载函数 wave = sin(omega * dytime) ;定义动荷载变量 end
apply xvel = 1 hist wave range z=-.1 .1 ;施加动荷载 apply zvel = 0 range z=-.1 .1 hist gp xvel 5,2,0 hist gp xvel 5,2,10 hist gp zvel 5,2,10 hist dytime
def tim ;估算程序运行的时间 tim = 0.01 * (clock - old_time) end set dyn multi on ;设置动态多步 solve age 1.0
print tim ;输出计算时间
print dyn ;输出动力计算相关信息 save mult1.sav 11-2 new
config dyn
gen zone brick size 1,1,50 model elas
prop shear 1e7 bulk 2e7 ini dens 1000 def setup
omega = 2.0 * pi * freq pulse = 1.0 / freq end
set freq=4.0 setup def wave
if dytime > pulse wave = 0.0 else
wave = 0.5 * (1.0 - cos(omega * dytime)) endif end
range name bottom z=-.1 .1
fix z range z=.5 55 ;将上部网格都施加数值向约束 apply dquiet squiet range bottom
apply sxz -2e5 hist wave syz 0.0 szz 0.0 range bottom ;-2e5的系数来源于 的值
apply nvel 0 plane norm 0,0,1 range bottom hist gp xvel 0,0,0 hist gp xvel 0,0,25 hist gp xvel 0,0,50 hist dytime hist wave
plot create hhh
plot add hist 1 2 3 vs 4 plot show solve age 2
11-3 new
;第一步：静力计算阶段 config dyn set dyn off
gen zone brick size 6 3 2
gen zone brick size 2 3 2 p0 0 0 2 gen zone brick size 2 3 2 p0 4 0 2 gen zone wedge size 1 3 2 p0 2 0 2
gen zone wedge size 1 3 2 p0 4 3 2 p1 3 3 2 p2 4 0 2 p3 4 3 4 &
p4 3 0 2 p5 4 0 4
model elastic
prop bulk 66667 shear 40000 ini dens 0.0025 set grav 0 0 -10
fix x range x -0.01 0.01 fix x range x 5.99 6.01 fix y range y -0.01 0.01 fix y range y 2.99 3.01 fix z range z -0.1 0.1 hist unbal solve
save 11-3_1.sav
;第二步：动力计算阶段 set dyn on def iniwave
per = 0.01 end iniwave def wave
wave = 0.5 * (1.0 - cos (2*pi*dytime/per)) end
free x y z ran z -0.1 0.1 ;去掉模型底部原有的静力条件
apply nquiet squiet dquiet ran z -0.1 0.1 ;静态边界条件
apply dstress 1.0 hist wave ran z -0.1 0.1 ;加动力荷载 apply ff ;施加自由场边界条件 group ff_corner
group ff_side ran x 0 6 group ff_side ran y 0 3
group main_grid ran x 0 6 y 0 3
set dyn time = 0 ;设置动力计算从0s开始 hist reset ;清空已有的历史信息 hist unbal
hist dytime ; 主体网格 hist gp xvel 2 1 0
hist gp xvel 2 1 5.0 ; 柱体网格 hist gp xvel -1 -1 0
hist gp xvel -1 -1 5 ;平行于y方向的二维自由场网格 hist gp xvel -1 0 0 hist gp xvel -1 0 5.0
;平行于x方向的二维自由场网格 hist gp xvel 2 -1 0 hist gp xvel 2 -1 5.0 solve age 0.015 save 11-3_2.sav 11-4 conf dy
gen zone brick size 3,3,3 model elas
prop bulk 1e8 shear 0.3e8 ini dens 1000 fix z range z -.1,.1
set dyn=on, grav 0 0 -10, hist_rep=1 hist gp zdisp 3.0,1.5,3.0 hist dytime plot create hh plot add his 1 vs 2 save 11-4_1.sav
;保存文件，在后续计算中需调用该文件 cyc 150 11-5
;（1）质量分量和刚度分量共同作用 rest 11-4_1.sav
set dyn damp rayleigh 1 22.8 solve age=0.2 title
vertical displacement versus time (mass & stiffness damping) plot show pause
;（2）只有质量分量 rest 11-4_1.sav
set dyn damp rayleigh 2 22.8 mass solve age=0.08 title
vertical displacement versus time (mass damping only) plot show pause
;（3）只有刚度分量 rest 11-4_1.sav
set dyn damp rayleigh 2 22.8 stiffness solve age=0.08 title
vertical displacement versus time (stiffness damping only) plot show 11-6
rest 11-4_1.sav
set dyn damp rayleigh 0.05 22.8 set hist_rep=5 solve age=0.5 title
vertical displacement versus time (5% Rayleigh damping) plot show pause
rest 11-4_1.sav
set dyn damp local 0.1571 ; = pi * 0.05 set hist_rep=5 solve age=0.5 title
vertical displacement versus time (5% Local damping) plot show 11-7
;振动台试验的例子 new
config dynamic fluid def model_dim h_R = 0
h_R1 = h_R + 1.0 end
model_dim gen zon bri p0 0 0 -10 p1 30 0 -10 p2 0 1 -10 p3 0 0 0 p4 30 1 -10 p5 0 1 0 p6 30 0 h_R p7 30 1 h_R size 30 1 10 group sand
gen zon bri p0 0 0 0 p1 30 0 h_R p2 0 1 0 p3 0 0 1 p4 30 1 h_R p5 0 1 1 p6 30 0 h_R1 p7 30 1 h_R1 size 30 1 1 group top
;gen zon bri p0 0 0 -.5 p1 3 0 -.5 p2 0 1 -.5 p3 0 0 0 p4 3 1 -.5 p5 0 1 0 p6 3 0 0 p7 3 1 0 size 30 1 10 model elastic
prop bulk=3e7 shear=1e7 fric=35 ini dens 2000 model fl_iso
prop poro 0.5 perm 1e-8 ini fmod 2e8 fdens 1000
ini pp 0 grad 0 0 -10e3 ran z 0 -10.0 fix z ran z -9.9 -10.1 fix x ran x -.1 .1 fix x ran x 29.9 30.1 fix y
set grav 10
set fluid off dyn off ini fmod 0
set mech rat 1e-6 solve
def ini_conf _k0 = 1.0
pnt = zone_head loop while pnt # null
val = _k0 * z_szz(pnt) + (_k0-1.) * z_pp(pnt) z_sxx(pnt) = val z_syy(pnt) = val pnt=z_next(pnt) end_loop end ini_conf solve
save 11-7.sav 11-8
rest 11-7.sav
set dyn on fluid on ini fmod 2e8 set fluid pcut on
model finn ran gro sand
prop bulk=3e7 shear=1e7 fric=35 ran gro sand endloop ini dens 2000 ran gro sand end
prop ff_latency=50 ff_switch = 0 ff_c1=0.8 solve_ages ff_c2=0.79 ff_c3=0.45 ff_c4=0.73 ran gro sand ;扭剪试验结果
def setup freq=5.0 ampl=2 omega = 2.0 * pi * freq end setup def sine_wave vv = 9.36e-2 * sin(omega*dytime) if dytime < 2.0 sine_wave = dytime / 2.0 * vv else if dytime < 20.0 sine_wave = vv else if dytime <= 30. sine_wave = (30.0 - dytime) / 10.0 * vv endif endif endif if dytime > 30.0 sine_wave = 0.0 endif end free x apply xvel = 1.0 hist sine_wave ran z -9.9 -10.1 apply xvel = 1.0 hist sine_wave ran x -.1 .1 apply xvel = 1.0 hist sine_wave ran x 29.9 30.1 set dyn damp local .314 call ppr.dat set hist_rep 100 set large ;set dyn dt 3e-4 set mech rat 1e-20 def solve_ages loop n(1,39) save_file = '11-8_'+string(n)+'s.sav' command sol age n ;save save_file endcommand save save_file
hist write 20 21 22 23 24 30 31 32 33 34 vs 2 file 10-8_Outfile_pp.txt hist write 40 41 42 43 44 vs 2 file 10-8_Outfile_xdis.txt 12-1 gen zon bri size 3 3 3 ini pp 30e3 grad 0 0 -10e3 plot con pp outline on ;outline的作用是在云图中显示网格的轮廓 12-2 gen zon bri size 3 3 3 water density 1000 set grav 10 water table ori 0 0 3 norm 0 0 1 12-3 gen zon bri size 3 3 3 ini pp 30e3 grad 0 0 -10e3 ran x -.1 .1 fix pp ran x -.1 .1 plot con pp outline on： 显示节点上的孔压云图 plot bcon pp outline on： 显示单元中心点处平均孔压块图 12-4 gen zon bri size 3 3 3 app pp 30e3 grad 0 0 -10e3 ran x -.1 .1 plot con pp outline on： 显示节点上的孔压云图 plot bcon pp outline on： 显示单元中心点处平均孔压块图 12-5 config fluid gen zone brick size 20 1 10 model mohr prop bulk 8.33e6 shear 3.85e6 fric 15 coh 10e3 tens 1e10; fix x range x -.1 .1 fix x range x 19.9 20.1 fix x y z range z -.1 .1