FLAC3D声发射模拟代码

new set logfile myfilesnn.log set log on gen zone cylinder p0 0 0 0 p1 0.0125 0 0 p2 0 0.5 0 p3 0 0 0.0125 size 8 12 8 gen zone reflect norm 1,0,0 gen zone reflect norm 0,0,1 model mohr prop young 10e9 poisson 0.25 dens 2500 cohe 1e4 ten 1e5 fric 35 set gravity 0 0 -10 apply zvel -.01 range z 19.9 20.1 fix z range z -.1 .1 def ax_stress f_accum=0.0 pnt = gp_head loop while pnt#null if gp_zpos(pnt)<0.1 then f_accum=f_accum +gp_zfunbal(pnt) endif pnt=gp_next(pnt) endloop ax_stress=f_accum/900 end ax_stress def ax_strain str1=0 n=1.0 array arr(6) p_z=zone_head loop while p_z#null if z_zcen(p_z)<0.6 then oo=z_fsi(p_z,arr) str1=str1-arr(3) endif p_z=z_next(p_z) n=n+1.0 endloop ax_strain=str1/n end ax_strain table 1 0,0 0,0 def get_plast shearnow=1 tensionnow=2 shearpast=4 tensionpast=8 v1=0 v2=0 v4=0 v8=0 p_z=zone_head loop while p_z#null pid=z_id(p_z) if z_state(p_z,0)=shearnow then i=table_size(1) loop n(1,i) y=table(1,n) if pid#y then v1=v1+1 table(1,n)=pid endif endloop endif if z_state(p_z,0)=tensionnow then i=table_size(1) loop n(1,i) y=table(1,n) if pid#y then v2=v2+1 table(1,n)=pid endif endloop endif if z_state(p_z,0)=shearpast then i=table_size(1) loop n(1,i) y=table(1,n) if pid#y then v4=v4+1 table(1,n)=pid endif endloop endif if z_state(p_z,0)=tensionpast then i=table_size(1) loop n(1,i) y=table(1,n) if pid#y then v8=v8+1 table(1,n)=pid endif endloop endif p_z=z_next(p_z) endloop end get_plast hist n 1 ;plot zone colorby state plot block state set fishcall 0 get_plast his unbal his v1 his v2 his v4 his v8 his ax_stress his ax_strain step 430 set log off ;hist write 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 file xyz.txt