flac3d-三轴模拟命令流

nimpgrid1groupsoilrangegrouprocknotplotblockgroup;assignpropertiesmodelmohrproden2200bulk4e7she1.5e7co8e4fric28ten9e3dil10rangegroupsoilproden2900bulk2.5e10she1.2e10co4e6fric45ten2e6dil15rangegrouprock;boundaryconditionsfixxyzrangez-0.10.1fixxyzrangez7.98.1fixyappnstress-1e5rangex-0.10.1appnstress-1e5rangex4.95.1pause;*************inizvel-2.5e-5rangez7.98.1;-*************defaxi_stressf_accum=0.0pnt=gp_headloopwhilepnt#nullifgp_zpos(pnt)0.15thenf_accum=f_accum+gp_zfunbal(pnt)endifpnt=gp_next(pnt)endloopaxi_stress=f_accum/0.5enddeflat_stressg_accum=0.0pnt=gp_headloopwhilepnt#nullifgp_xpos(pnt)0.1theng_accum=g_accum+gp_xfunbal(pnt)endifpnt=gp_next(pnt)endlooplat_stress=g_accum/0.8endhistn1histaxi_stresshistlat_stresshistgpzdis008histgpxdis004histgpxdis054histunbal;**************plothist-1vs-3;axialstressvsaxialdisplace;**************step15000save1.savhistwrite1fileaxi_stress.dathistwrite2filelat_stress.dathistwrite3fileaxi_dis.dathistwrite4filelat_dis1.dathistwrite5filelat_dis2.dathistwrite6fileunbalance.dat