5、100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick &p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 60 0 60 &p4 100 2 40 p5 60 2 60 p6 100 0 60 p7 100 2 60 &size 30 1 10plot surf;初始应力场生成model elapro bulk 1e9 shear 2e9 dens 2500fix yfix x range x -.1 .1fix x range x 99.9 100.1fix
6、x y z range z -.1 .1set gravity 0 0 -10set mech ratio 1e-3solve ini xvel 0 yvel 0 zvel 0ini xdis 0 ydis 0 zdis 0 model mohrpro bulk 8.3e7 shear 3.8e7 coh 4.2e4 dila 20 ten 0.8e6 fric 17 dens 2500 save dizhen1.sav restore dizhen1.sav;动力计算部分def iniwavefr
7、eq =5.0endIniwave;定义动荷载（以速度形式）def wavevv=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)if dytime<0.2 wave=dytime/0.2*vv else if dytime<2.0 wave=vv else if dytime <=3.0 wave=(3.0-dytime)/1.0*vv endif e
8、ndifendif if dytime>3.0 wave=0.0endifend set dyn onfree x y z ran z -.1 .1 ;去除静力边界apply nquiet squiet dquiet ran z -.1 .1 ;在底部施加静态边界apply xvel 1.0 hist wave range z -.1 .1 ;施加动荷载apply ff;施加自由场边界 pause khist id =1 gp xvel 15 0 0hist id =2 gp xvel 15 0 10hist i
9、d =3 gp xvel 15 0 20hist id =4 gp xvel 15 0 30hist id=5 wavehist id=6 dytimeplot create hhhplot showset dyn time=0set dyn damp local .314 ;局部阻尼solve age 4plot hist 5 vs 6 pause Kplot 1 2 3 4 vs 6 Fig.1水平位移云图（最大水平位移53cm） Fig.2 剪切速率云图及速度矢量图 Fig.3地震波速度时程图
