陈育民《FLAC3D基础与工程实例》全部命令流

《FLAC/FLAC 3D基础与工程实例》全部命令流 1 gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1.0 0.9 1.1 group brick_1 gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group brick_2 gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot sur 2-1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3-1 ;-------------------------------------工程信息 ;Project Record Tree export ;Title:Simple test ;---------------------------------计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81 history 999 unbalanced solve save state1.sav ;----------------------------------计算第二步 ;... STATE: STATE2 .... initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete 'null' history 1 xdisp i=4, j=11 solve save state2.sav ;--------------------------------绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis-A plot hold history 1 line 5-1 new ; =============================== ; 定义球体半径和半径方向上单元网格数 ; =============================== def parm rad=10.0 rad_size=5 end parm ; =============================== ; 建立八分之一球体外接立方体网格 ; =============================== gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 1 gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 2 gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 3 ; ================================== ; 利用FISH语言将内部立方体节点调整到球面 ; ================================== def make_sphere p_gp=gp_head loop while p_gp#null ; 获取节点点坐标值：P=(px,py,pz) px=gp_xpos(p_gp) py=gp_ypos(p_gp) pz=gp_zpos(p_gp) dist=sqrt(px*px+py*py+pz*pz) if dist>0 then ; 节点位置调整 maxp=max(px,max(py,pz)) k=(maxp/rad)*(rad/dist) gp_xpos(p_gp)=k*px gp_ypos(p_gp)=k*py gp_zpos(p_gp)=k*pz end_if p_gp=gp_next(p_gp) end_loop end make_sphere ; =============================== ; 利用镜像生成完整球体网格 ; =============================== gen zone ref gen zone ref dip 90 gen zone ref dip 90 dd 90 ; =============================== ; 显示球体网格 ; =============================== plot surf pl set back wh pl bl gr 5-3 n gen zon bri size 1 1 2 group soil ran z 1 1 group rock ran z 0 1 expgrid 1.flac3d 6-1 gen zon bri size 3 3 3 model mohr prop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 ini dens 2000 hist unbal set grav 10 solve elastic save 6-1.sav 6-2 rest 6-1.sav ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 solve save 6-2.sav 6-3 rest 6-1.sav ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3 hist id=3 gp zdis 1 1 2 hist id=4 gp xdis 1 1 3 hist id=5 gp xdis 1 1 3 hist id=6 zone szz 1 1 3 hist id=7 zone szz 1.5 1.5 2.5 hist id=8 zone sxz 1.5 1.5 2.5 solve save 6-3.sav 6-4 rest 6-3.sav set log on set logfile 6-2.log print zone stress print gp dis set log off 6-5 rest 6-1.sav ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30 plot con szz ou on magf 10 plot add hist 1 set movie avi step 1 file 6-5.avi movie start solve movie finish 7-1 n gen zon bri size 1 1 2 model elas prop bulk 3e7 shear 1e7 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 set grav 0 0 -10 solve plo con sz 7-2 n gen zon bri size 1 1 2 model mohr prop bulk 3e7 shear 1e7 c 1e10 f 15 tension 1e10 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 set grav 0 0 -10 solve prop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0 solve plo con sz 7-3 n gen zone brick size 1 1 2 model mohr prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 set grav 0 0 -10 solve elas plo con sz 7-4 new gen zone brick size 1 1 2 model mohr prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 ini szz -40e3 grad 0 0 20e3 ran z 0 2 ini syy -20e3 grad 0 0 10e3 ran z 0 2 ini sxx -20e3 grad 0 0 10e3 ran z 0 2 set grav 0 0 -10 solve plo con sz 7-5 n gen zon bri size 1 1 2 model m prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 ran z 0 1 ini dens 1500 ran z 1 2 ini szz -35e3 grad 0 0 20e3 ran z 0 1 ini syy -17.5e3 grad 0 0 10e3 ran z 0 1 ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1 ini szz -15e3 grad 0 0 15e3 ran z 1 2 ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2 ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2 ini pp 10e3 grad 0 0 -10e3 ran z 0 1 set grav 0 0 -10 solve plo con sz 7-6 n gen zon bri size 1 1 2 model m prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 ran z 0 2 ini szz -50e3 grad 0 0 20e3 ran z 0 1 ini syy -30e3 grad 0 0 10e3 ran z 0 1 ini sxx -30e3 grad 0 0 10e3 ran z 0 1 ini pp 30e3 grad 0 0 -10e3 ran z 0 2 app nstress -10e3 ran z 2 set grav 0 0 -10 solve plo con sz 7-7 new gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10 model elas pro bulk 10e10 she 10e10 ini den 2500 apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120 set grav 0 0 -10 step 30000 ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 plo cont szz 7-8 new gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10 model elas pro bulk 10e10 she 10e10 ini den 2500 ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 ini sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 ini syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 ini syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 ini szz -1e8 ran z -.1 .1 fix x y z ran z -.1 .1 set grav 0 0 -10 solve ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 plo cont szz 7-9 new gen zone brick p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group clay gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio 1.1 1 1 size 12 1 10 group clay gen zone brick p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio 1 1 0.8 size 8 1 4 group soil gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio 1.1 1 0.8 size 12 1 4 group soil gen zone brick p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 & p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group dam fix x y z ran z -49.9 -50.1 fix x ran x -.1 .1 fix x ran x 99.9 100.1 fix y model mohr ran z -50 0 model null ran z 0 5 prop bulk 7.8e6 shear 3.0e6 coh 10e10 tension 1e10 ran group soil ini dens 1500 ran group soil prop bulk 3.91e6 shear 1.5e6 coh 10e10 tension 1e10 ran group clay ini dens 1800 ran group clay set grav 0 0 -9.8 hist id=1 unbal solve prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil prop bulk 3.91e6 shear 1.5e6 coh 20e3 fric 20 ran group clay solve save elastic.sav ini xdis 0 ydis 0 zdis 0 ;将节点位移清零 ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 hist id=2 gp zdis 0 0 0 ;地基顶部中心点的沉降 hist id=3 gp zdis 27.5 0 0 ;记录路基坡脚处的沉降 hist id=4 gp xdis 27.5 0 0 ;记录路基坡脚处的水平位移 model elastic ran z 0 1 ; ;激活0 m ~ 1 m的单元 prop bulk 7.8e6 shear 3.0e6 ran z 0 1 ini dens 1500 ran z 0 1 solve ;按软件默认精度求解 save fill-1.sav model elastic ran z 1 2 prop bulk 7.8e6 shear 3.0e6 ran z 1 2 ini dens 1500 ran z 1 2 solve save fill-2.sav model elastic ran z 2 3 prop bulk 7.8e6 shear 3.0e6 ran z 2 3 ini dens 1500 ran z 2 3 solve save fill-3.sav model elastic ran z 3 4 prop bulk 7.8e6 shear 3.0e6 ran z 3 4 ini dens 1500 ran z 3 4 solve save fill-4.sav model elastic ran z 4 5 prop bulk 7.8e6 shear 3.0e6 ran z 4 5 ini dens 1500 ran z 4 5 solve save fill-5.sav pau ;plo bl gr ;gen zone brick p0 0 0 0 p1 100 0 0 p2 0 5 0 p3 0 0 5 size gen zone brick p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group clay gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio 1.1 1 1 size 12 1 10 group clay gen zone brick p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio 1 1 0.8 size 8 1 4 group soil gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio 1.1 1 0.8 size 12 1 4 group soil gen zone brick p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 & p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group dam 7-10 set log on ;打开log记录 set logfile 1.log ;设置记录文件名为：1.log restore fill-1.sav ;调用保存的文件 print gp dis range id 517 any id 533 any ;输出两个节点的变形值 restore fill-2.sav print gp dis range id 517 any id 533 any restore fill-3.sav print gp dis range id 517 any id 533 any restore fill-4.sav print gp dis range id 517 any id 533 any restore fill-5.sav print gp dis range id 517 any id 533 any set log off ;关闭log记录 8-1 def abc abc = 1 + 2 * 3 abcd = 1.0 / 2.0 end abc print fish 8-2 def abc if aa < 0 then abc = 0.0 else abc = 2.0 * aa endif end abc 8-3 def abc loop aa (1, 2.5) command print aa endcommand endloop end 8-4 new gen zon bri size 3 3 3 model elastic prop bu 3e7 sh 1e7 ini dens 2000 fix x y z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 10 solve ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 save 8-4.sav 8-5 rest 8-4.sav def E_modify p_z = zone_head d_k = 704 d_n = 0.38 d_pa = 101325.0 ;//大气压 loop while p_z # null sigma_3 = -1.0 * z_sig1(p_z) E_new = d_k * d_pa * (sigma_3 / d_pa) ^ d_n z_prop(p_z,'young') = E_new p_z = z_next(p_z) endloop end E_modify 8-6 rest 8-5.sav table 1 name load_settlement def add_load p_gp = gp_near(2,1,3) loop n (1,5) app_load = n * (-1000e3) file_name = '7-6_add_step' + string(n) + '.sav' command app nstress app_load ran z 2.9 3.1 x 1 2 y 1 2 solve save file_name endcommand xtable(1,n) = -1.0 * app_load ytable(1,n) = gp_zdisp(p_gp) endloop end add_load save 8-6.sav 8-7 rest 8-6.sav def find_max_disp p_gp = gp_head maxdisp_value = 0.0 maxdisp_gpid = 0 loop while p_gp # null disp_gp = sqrt(gp_xdisp(p_gp) ^ 2 + gp_ydisp(p_gp) ^ 2 + gp_zdisp(p_gp) ^ 2) if disp_gp > maxdisp_value maxdisp_value = disp_gp maxdisp_gpid = gp_id(p_gp) endif p_gp = gp_next(p_gp) endloop end find_max_disp print maxdisp_value maxdisp_gpid rest 8-6.sav config zextra 1 def get_sigma_dif p_z = zone_head loop while p_z # null sigma_dif = z_sig3(p_z) - z_sig1(p_z) z_extra(p_z,1) = sigma_dif p_z = z_next(p_z) endloop end get_sigma_dif plot con zextra 1 9-1 ;--------------------------------------------------- ; ; 移来移去法接触面的建立 ;--------------------------------------------------- n gen zone radcyl p0 (0,0,0) p1 (8,0,0) p2 (0,0,-5) p3 (0,8,0) & p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) & p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-5) p11 (0,.3,-5) & size 3 10 6 15 ratio 1 1 1 1.15 gen zone radcyl p0 (0,0,-5) p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5) & p4 (8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) & p8 (.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) & size 3 6 6 15 ratio 1 1 1 1.15 fill gen zone reflect dd 270 dip 90 group clay ; interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .31 & cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .29 not interface 2 face range cylinder end1 (0,0,-4.9) end2 (0,0,-5.1) radius .31 ; pause gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1) p3 (0,.3,6) & p4 (.3,0,1) p5 (0,.3,1) & size 3 10 6 gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) & p4 (.3,0,6) p5 (0,.3,6) & size 3 1 6 gen zone reflect dd 270 dip 90 range z 1 6.1 group pile range z 1 6.1 pause ini z add -6.0 range group pile save pile_geom.sav 9-2 ;导来导去法 ;--------------------------------------------- n gen zone brick size 3 3 3 group 2 range x 1 2 y 1 2 z 1 2 group 1 range gr 2 not save 1.sav del ran group 2 not interface 1 face range x 1 y 1 2 z 1 2 interface 1 face range x 2 y 1 2 z 1 2 interface 1 face range x 1 2 y 1 z 1 2 interface 1 face range x 1 2 y 1 2 z 1 interface 1 face range x 1 2 y 1 2 z 2 rest 1.sav del ran group 2 expgrid 1.fac3d impgrid 1.flac3d model ela fix x y z ran z 0 ini den 2000 set grav 0 0 -10 interface 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15 app nstr -200e3 ran x 0 1 y 1 2 z 3 solve 9-3 ;切割模型法 ;------------------------------ n gen zone brick size 3 3 3 group 1 range x 1 2 y 1 2 z 2 3 group 2 range group 1 not gen separate 1 int 1 wrap 1 2 int 1 maxedge 0.5 plo int red 9-4 ; Create Material Zones gen zone brick size 5 5 5 & p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) & p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5) gen zone brick size 5 5 5 p0 (0,0,5) edge 5.0 group Material ; Create Bin Zones gen zone brick size 1 5 5 & p0 (3,0,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) & p6 add (3,0,5) p7 add (3,6,5) gen zone brick size 1 5 5 & p0 (5,0,5) p1 add (1,0,0) p2 add (0,5,0) & p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) & p6 add (1,0,5) p7 add (1,6,5) gen zone brick size 5 1 5 & p0 (0,3,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) & p6 add (5,2,5) p7 add (6,3,5) gen zone brick size 5 1 5 & p0 (0,5,5) p1 add (5,0,0) p2 add (0,1,0) & p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) & p6 add (5,0,5) p7 add (6,1,5) group Bin range group Material not ; Create named range synonyms range name=Bin group Bin range name=Material group Material ; Assign models to groups model mohr range Material model elas range Bin gen separate Material interface 1 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 above interface 2 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 below int 1 maxedge 0.55 int 2 maxedge 0.55 ; Assign properties prop shear 1e8 bulk 2e8 fric 30 range Material prop shear 1e8 bulk 2e8 range Bin ini den 2000 int 1 prop ks 2e9 kn 2e9 fric 15 int 2 prop ks 2e9 kn 2e9 fric 15 ; Assign Boundary Conditions fix x range x -0.1 0.1 any x 5.9 6.1 any fix y range y -0.1 0.1 any y 5.9 6.1 any fix z range z -0.1 0.1 Bin ; 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 ) &