/SOLU
ANTYPE,4 !定义分析类型
TRNOPT,FULL !瞬态分析
LUMPM,0
NSUBST,20,0,0 !定义加载步的子步长
!ALPHAD,0.452709 !定义阻尼矩阵
!BETAD,0.004295 !定义刚度矩阵
n=750
*dim,time,array,n,1,1
*dim,acc,array,n,1,1
*dim,acc1,array,3000,1,1
*CREATE,ansuitmp
*VREAD,acc1,'aaaa','txt',' ',1, , , , !El centro地震波文件aaaa.txt,单位cm/s^2
*END
/INPUT,ansuitmp
*do,j,1,n
time(j,1,1)=j*0.02
acc(j,1,1)=acc1(0+j*1,1,1)*0.5*0.01*392/341.7 !注意,要转换成0.4g峰值,并且单位为m/s^2
*enddo
allsel
*do,i,1,n,1
acel,acc(i),0,0
time,i*0.02
solve
*enddo
FINISH
这是我用的地震波加载命令流,麻烦您给看看 有什么不对的地方吗?至于地震波文件就是一列数值,非常感谢!急! |