FLAC3D工程实例之MILES FIND_ADD CONT_EX CAVPRO CAV.doc

1、FLAC3D工程实例之MILES.dat（calculate and store Von Mises stress in zone extension 1）config zextra 1def mises; - calculate and store Von Mises stress in zone extension 1 -p_z = zone_headmax_mises = 0.0loop while p_z # nullmstr = (sxx(p_z) + syy(p_z) + szz(p_z) / 3.dsxx = sxx(p_z) - mstrdsyy = syy(p_z) - mstrdszz = szz(p_z) - mstrdsxy = sxy(p_z)dsxz = sxz(p_z)dsyz = syz(p_z)vmstr2 = 1.5 * (dsxx*dsxx + dsyy*dsyy + dszz*dszz)vmstr2 = vmstr2 + 3. * (dsxy*dsxy + dsxz*dsxz + dsyz*dsyz)if vmstr2 0.0 thenz_extra(p_z,1) = sqrt(vmstr2)elsez_extra(p_z,1) = 0.0endifmax_mises = max(max_mises,z_extra(p_z,1)p_z = z_next(p_z)end_loopendmisesplot cont zextra 1 alias Von Mises Stress averageprint max_mises2、FLAC3D工程实例之FIND_ADD.dat（find gridpoints at pile cap）; find gridpoints at pile cap;def find_addtop_head = nullgp_pnt = gp_headloop while gp_pnt # nullif gp_zpos(gp_pnt) 0.05 thennew = get_mem(2)mem(new) = top_headmem(new+1) = gp_pnttop_head = newendifgp_pnt = gp_next(gp_pnt)endloopendfind_add3、FLAC3D工程实例之CONT_EX.dat（完成当前开挖继续以下开挖）; finishes stepping for the current excavation step; and continues following excavation steps;def cont_excavrest_step = 3000*(step/3000+1)-stepif rest_step = 3000 thenrest_step = 0end_ifcommandstep rest_stepend_command;store displacements in tablesxtable(1,cut) = 3.0 * cutytable(1,cut) = zdisp(ipt_surf)xtable(2,cut) = 3.0 * cutytable(2,cut) = zdisp(ipt_crown)xtable(3,cut) = 3.0 * cutytable(3,cut) = zdisp(ipt_spring)commandsave m1.savend_commandif cut=10 thencommandsave m1_30.savend_commandend_ifexcavend4、FLAC3D工程实例之（剑桥模型介质的圆柱形空腔膨胀排水结果）;-; ; results from undrained cylindrical cavity expansion; in Cam-Clay medium;-res cav.sav; table 10 : pore pressure; table 11 : radial-component of effective stress (xx-stress); table 12 : axial-component of effective stress (yy-stress); table 13 : tangential-component of effective stress (zz-stress); rval : ln(r/a0)def b_tablepnt = zone_headloop while pnt # nullxval = xcen(pnt)yval = ycen(pnt)rval = 0.5*ln(xval*xval+yval*yval)c_pp = z_pp(pnt)table(10,rval) = c_pptable(11,rval) = -(sxx(pnt)*cc+syy(pnt)*ss+sxy(pnt)*sc) - c_pptable(12,rval) = -(sxx(pnt)*ss+syy(pnt)*cc-sxy(pnt)*sc) - c_pptable(13,rval) = -szz(pnt) - c_pppnt = z_next(pnt)end_loopendb_table; results -plot his 9 vs 4 ; Fig 9.3plot his -2 3 vs 4 ; Fig 9.4plot table 10 line 11 line 12 line 13 line ; Fig 9.5ret5、FLAC3D工程实例之CAV（剑桥模型介质的圆柱形空腔膨胀排水）;-; cav.dat; Undrained cylindrical cavity expansion in Cam-Clay medium;-titleUndrained cylindrical cavity expansion in Cam-Clay mediumconfig fluid; - model geometry -gen zone brick p0 1.0 0.0 1.0 p1 200.0 0.0 1.0 &p2 1.0 0.0 0.0 p4 200.0 0.0 0.0 &p3 0.9877 0.1564 1.0 p6 197.5377 31.2869 1.0 &p5 0.9877 0.1564 0.0 p7 197.5377 31.2869 0.0 &size 31 1 1 ratio 1.1 1 1def find_idpnt1 = gp_near(1.,0.,1.)pnt2 = gp_near(1.,0.,0.)pnt3 = gp_near(0.9877,0.1564,1.)pnt4 = gp_near(0.9877,0.1564,0.)id1 = gp_id(pnt1)id2 = gp_id(pnt2)id3 = gp_id(pnt3)id4 = gp_id(pnt4)pntz = z_near(1.,0.07,0.5)cc = cos(4.5*degrad)ss = sin(4.5*degrad)sc = 2.*cc*sscc = cc*ccss = ss*ssb_vel= 1.e-5b_vx = 0.9877*b_velb_vy = 0.1564*b_velendfind_idrange name out cyl end1 0 0 -10 end2 0 0 10 rad 197.5377 notrange name in cyl end1 0 0 -10 end2 0 0 10 rad 1.001; - model properties -mo cam-clayprop shear 74.prop mm 1.2 lambda 0.15 kappa 0.03 mp1 1.0 mv_l 2.3mo fl_iso; - boundary conditions -fix y z range y -.001 .001apply nvel 0.0 plane dip 90.0 dd -9.0 origin 0.0 0.0 0.0 &range in not y 0.07 200.0apply dvel 0.0 plane dip 90.0 dd -9.0 origin 0.0 0.0 0.0 &range in not y 0.07 200.0apply nstress -1.65 range outfix x y z range inini xvel b_vel range id id1 id1ini xvel b_vel range id id2 id2ini xvel b_vx yvel b_vy range id id3 id3ini xvel b_vx yvel b_vy range id id4 id4ini sxx -1.65 syy -1.65 szz -3.; model settings -set fluid biot onset fluid offset large; - fish functions -; . initial specific volume, tangent bulk modulus, porosity .def set_proppnt = zone_heads1 = -z_sxx(pnt)s2 = -z_syy(pnt)s3 = -z_szz(pnt)p0 = (s1 + s2 + s3) / 3.q0 = sqrt(s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1)*0.5)e0 = q0 / (z_prop(pnt,mm) * p0)pc = p0 * (1. + e0 * e0)p1 = z_prop(pnt,mp1)vc = z_prop(pnt,mv_l) - z_prop(pnt,lambda) * ln(pc/p1)v0 = vc - z_prop(pnt,kappa) * ln(p0/pc)b0 = v0 * p0 / z_prop(pnt,kappa)n0 = (v0 - 1.) / v0bb = 10. * b0bi = 100. * b0commandprop mpc pc bulk_bound bbini biot_mod biend_commandend; . numerical variables .def pathsrad = z_sxx(pntz)*cc+z_syy(pntz)*ss+z_sxy(pntz)*scp_fl = z_pp(pntz)c_tp = -(z_sxx(pntz)+z_syy(pntz)+z_szz(pntz) / 3.c_ep = z_prop(pntz,cp)c_q = z_prop(pntz,cq)c_cq = c_ep * z_prop(pntz,mm)c_v = z_prop(pntz,cv)c_a =