flac3d地震作用分析命令流

newconfigdynsetdynoffgenzonebrickp0000p110000p2020p30040size50110genzonebrick&p040040p1100040p240240p360060&p4100240p560260p6100060p7100260&size30110plotsurf;初始应力场生成modelelaprobulk1e9shear2e9dens2500fixyfixxrangex-.1.1fixxrangex99.9100.1fixxyzrangez-.1.1setgravity00-10setmechratio1e-3solveinixvel0yvel0zvel0inixdis0ydis0zdis0modelmohrprobulk8.3e7shear3.8e7coh4.2e4dila20ten0.8e6fric17dens2500savedizhen1.savrestoredizhen1.sav;动力计算部分definiwavefreq=5.0endIniwave;定义动荷载（以速度形式）defwavevv=9.36e-2*sin(2.0*pi*freq*dytime);wave=0.5*(1.0-cos(2*pi*dytime/per));wave=sin(2*pi*dytime/per)ifdytime0.2wave=dytime/0.2*vvelseifdytime2.0wave=vvelseifdytime=3.0wave=(3.0-dytime)/1.0*vvendifendifendififdytime3.0wave=0.0endifendsetdynonfreexyzranz-.1.1;去除静力边界applynquietsquietdquietranz-.1.1;在底部施加静态边界applyxvel1.0histwaverangez-.1.1;施加动荷载applyff;施加自由场边界pausekhistid=1gpxvel1500histid=2gpxvel15010histid=3gpxvel15020histid=4gpxvel15030histid=5wavehistid=6dytimeplotcreatehhhplotshowsetdyntime=0setdyndamplocal.314;局部阻尼solveage4plothist5vs6pauseKplothist1234vs6