FLAC3D工程实例之MILESFINDADDCONTEXCAVPROCAV.doc

1、FLAC3D工程实例之MILES.dat（calculate and store Von Mises stress in zone extension 1） config zextra 1 def mises ; --- calculate and store Von Mises stress in zone extension 1 --- p_z = zone_head max_mises = 0.0 loop while p_z # null mstr = (sxx(p_z) + syy(p_z) + szz(p_z)) / 3. dsxx = sxx(p_z) - mstr dsyy = syy(p_z) - mstr dszz = szz(p_z) - mstr dsxy = 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 then z_extra(p_z,1) = sqrt(vmstr2) else z_extra(p_z,1) = 0.0 endif max_mises = max(max_mises,z_extra(p_z,1)) p_z = z_next(p_z) end_loop end mises plot cont zextra 1 alias Von Mises Stress average print max_mises 2、FLAC3D工程实例之FIND_ADD.dat（find gridpoints at pile cap） ; find gridpoints at pile cap ; def find_add top_head = null gp_pnt = gp_head loop while gp_pnt # null if gp_zpos(gp_pnt) 0.05 then new = get_mem(2) mem(new) = top_head mem(new+1) = gp_pnt top_head = new endif gp_pnt = gp_next(gp_pnt) endloop end find_add 3、FLAC3D工程实例之CONT_EX.dat（完成当前开挖继续以下开挖） ; finishes stepping for the current excavation step ; and continues following excavation steps ; def cont_excav rest_step = 3000*(step/3000+1)-step if rest_step = 3000 then rest_step = 0 end_if command step rest_step end_command ;store displacements in tables xtable(1,cut) = 3.0 * cut ytable(1,cut) = zdisp(ipt_surf) xtable(2,cut) = 3.0 * cut ytable(2,cut) = zdisp(ipt_crown) xtable(3,cut) = 3.0 * cut ytable(3,cut) = zdisp(ipt_spring) command save m1.sav end_command if cut=10 then command save m1_30.sav end_command end_if excav end 4、FLAC3D工程实例之cav.pro（剑桥模型介质的圆柱形空腔膨胀排水结果） ;----------------------------------------------------------- ; cav.pro ; 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 eff