FLAC3D命令流

; Monitor histories hist unbal hist gp zdisp (6,6,10) hist gp zdisp (0,0,10) hist gp zdisp (0,0,0) ; Settings set large set grav 0,0,-10 ; Cycling step 4000 save bin.sav 9-5 ; ;------------------------------------------------------------------ rest pile_geom.sav model mohr range group clay prop bulk 8.333e7 shear 3.846e7 coh 30000 fric 0 range group clay model elas range group pile prop bulk 8.333e7 shear 3.846e7 range group pile interface 1 prop kn 1e8 ks 1e8 fric 20 coh 30000 interface 2 prop kn 1e8 ks 1e8 fric 20 coh 30000 ; ini dens 1230 range group clay ini dens 1230 range group pile model null range z -0.1 0.15 ; fix z range z -8.1 -7.9 fix x range x -8.1 -7.9 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 7.9 8.1 set grav 0 0 -10 ini szz 0. grad 0 0 12300. range z -5.5 0. ini szz 17600 grad 0 0 15500 range z -8 -5.5 ini sxx 0. grad 0 0 5271.4 range z -5.5 0. ini sxx 7542.86 grad 0 0 6642.86 range z -8 -5.5 ini sxx add 31428.6 grad 0 0 5714.3 range z -8 -5.5 ini syy 0. grad 0 0 5271.4 range z -5.5 0. ini syy 7542.86 grad 0 0 6642.86 range z -8 -5.5 ini syy add 31428.6 grad 0 0 5714.3 range z -8 -5.5 ; water density 1000 water table origin 0,0,-5.5 normal 0 0 -1 ini dens 1550 range z -8 -5.5 hist unbal ; solve rat 1.e-6 save pile0.sav ; model elas range group pile prop bulk 13.9e9 shear 10.4e9 range group pile ini dens 2500 range group pile call find_add.fis solve rat 1.e-6 save pile1.sav rest pile1.sav ;调用保存文件 ini state 0 ini xdis 0.0 ydis 0.0 zdis 0.0 ;位移清零

apply szz -0.4e6 range z 0.05 0.15 group pile ;桩顶加第一级荷载 solve save app0.4.sav print gp disp range id 1 ;输出第一级荷载下的桩顶位移，假定桩顶中心的id号为1 apply szz -0.6e6 range z 0.05 0.15 group pile ;桩顶加第二级荷载 solve save app0.6.sav print gp disp range id 1 ;输出第二级荷载下的桩顶位移 ???????????????? ;依次加载，直到桩破坏 ;-------------------------------------------------------------------- ;速度加载法 rest pile1.sav ini state 0 ini xdis 0 ydis 0 zdis 0 def zs_top ;检测桩顶竖向荷载 ad = top_head zftot = 0.0 loop while ad # null gp_pnt = mem(ad+1) zf = gp_zfunbal(gp_pnt)

zftot = zftot + zf ad = mem(ad) endloop

zs_top = zftot / 0.1414 end

fix z range z 0.05 .15 group pile ;固定桩顶速度，用速度来确定位移 def ramp

while_stepping if step < ncut then

udapp = float(step) * udmax / float(ncut) else

udapp = udmax endif

ad = top_head

loop while ad # null gp_pnt = mem(ad+1) gp_zvel(gp_pnt) = udapp ad = mem(ad) endloop end

hist gp zdis 0,0,0 hist gp zvel 0,0,0 hist zs_top

hist zone szz 0,0,-.1 set mech damp comb

set udmax = -1e-8 ncut 30000 step 225000 save pile2.sav ;

;--------------------------------------------------------------------

;位移控制法 def solve_steps loop n (1,21)

save_file = string(n) + '-step.sav' command step 40000 save save_file

pri zone stress ran id 2381 a id 2361 a id 2341 a ;输出桩顶网格单元的应力 endcommand endloop end

solve_steps

;----------------------------------------------------------------------------- 10-1

sel beam beg 0 0 0 end 2 0 0 nseg 2 sel beam beg 2 0 0 end 4 0 -1 nseg 3

sel beam id=2 beg 4 0 -1 end 5 0 -2 nseg 2 plot sel geo id on nod on scale 0.04 plot ad ax 10-2

sel node id=1 0 0 0 sel node id=2 2 0 0 sel node id=3 4 0 -1 sel node id=4 5 0 -2

sel beamsel id=1 cid=1 node 1 2 sel beamsel id=1 cid=2 node 2 3 sel beamsel id=1 cid=3 node 3 4 plot sel geo id on nod on scale 0.04 plot ad ax 10-3

def set_vals

ptA = 25.0 * sin( 40.0*degrad ) ptB = 25.0 * cos( 40.0*degrad ) end set_vals

gen zone cylinder p0=( 0.0, 0.0, 0.0 ) &

p1=( ptA, 0.0, ptB ) & p2=( 0.0, 25.0, 0.0 ) & p3=( 0.0, 0.0, 25.0 ) & p4=( ptA, 25.0, ptB ) & p5=( 0.0, 25.0, 25.0 ) & size=(1, 2, 2)

sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &

end2=(0.0,25.0,0.0) radius=24.5 not plot blo gro

plot ad sel geom black black cid on scale=0.03 plot ad ax pau

delete ; delete all zones

sel node init zpos add -25.0 10-4 new

gen zone brick size 6 8 8 model mohr

prop bulk 1e8 shear 0.3e8 fric 35 prop coh 1e10 tens 1e10 set grav 0 0 -9.81 ini dens 1000

fix x range x -0.1 0.1 fix x range x 5.9 6.1 fix y range y -0.1 0.1 fix y range y 7.9 8.1 fix z range z -0.1 0.1 hist n 5 hist unbal

set mech force 50 solve

save beam-brace0.sav ;

prop coh 1e3 tens 1e3

model null range x 2 4 y 2 6 z 5 10 set large

ini xdis 0 ydis 0 zdis 0

sel beam begin=( 2, 4, 8) end=( 4, 4, 8) nseg=2 sel beam prop emod=2.0e11 nu=0.30

sel beam prop XCArea=6e-3 XCIz=200e-6 XCIy=200e-6 XCJ=0.0 hist gp zdisp 4 4 8 solve

save beam-brace1.sav ;

plot create GravV

plot set plane dip 90 dd 0 origin 3 4 0 plot set rot 15 0 20

plot set center 2.5 4.2 4.0 plot set cap size 25

plot add cont disp plane behind shade on plot add sel beam force fx

plot add sel geom black black node=off shrinkfac=0.0 plot add axes plot show 10-5

;非全长锚固、预紧力锚杆（锚索）模拟

;方法1、通过删除－建立link链接来模拟托盘

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 10e5 gr_k 2e7 range cid 11,17 sel delete link range id 1 sel link id=100 1 target zone

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

step 2000 sav 10-5.sav 10-6

;非全长锚固、预紧力锚杆（锚索）模拟 ;方法2、通过设置极大锚固剂参数模拟托盘

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 10e5 gr_k 2e7 range cid 11,17

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 10e8 gr_k 2e10 range cid 1,1

sel cable id=1 pretension 60e3 range cid 1,10 step 2000 sav 10-6.sav 10-7

;非全长锚固、预紧力锚杆（锚索）模拟

;方法3：借助别的结构单元（如liner单元）来模拟托盘

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6 sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 10e5 gr_k 2e7 range cid 11,17 sel liner range y=-.1, .1 x=-1,1 z=-1,1

sel liner PROP iso=( 25e9, 0.15) thick=0.1 ; concrete sel liner PROP cs_nk=8e8 cs_sk=8e8 &

cs_ncut=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 18

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

step 2000 sav 10-7.sav 10-8 n

title Structure_dynamic_analysis_lakewater config dyn

sel pile id=1 beg 0 0 0 end 0 0 1 sel pile prop dens 2400 &

Emod 1.0e10 Nu 0.3 XCArea 0.3 &

XCJ 0.16375 XCIy 0.00625 XCIz 0.01575 & Per 2.8 &

CS_sK 1.3e11 CS_sCoh 0.0 CS_sFric 10.0 &

CS_nK 1.3e11 CS_nCoh 0.0 CS_nFric 0.0 CS_nGap off

def f1

whilestepping

f0=10000*sin(10*dytime) np = nd_head

loop while np # null if nd_pos(np,1,3)=1 nd_apply(np,1)=f0 endif

np = nd_next(np) endloop end

sel node fix x y z xr yr zr ran id=1 sel set damp combined