1)是正弦波的地震波反应- /solu
- *dim,baseace1,array,600 ! 定义正弦波的“基础加速度”, 500时间点
- *do,i,1,600
- baseace1(i)=0.1*sin(i*2) !
- *enddo
- NSUBST,1, , ,1 !1个子步
- OUTRES,ALL,1 !输出每个子步的结果
- ANTYPE,TRANS
- CSYS,1
- NSEL,S,LOC,X,13000,13000
- CSYS,0
- *do,i,1,600
- f,all,fx,3e4*baseace1(i)
- f,all,fz,3e4*baseace1(i)
- TIME,i*0.02
- solve
- *enddo
复制代码 2)做z向行波分析- /nerr,,99999999
- fini
- /solu
- *dim,baseace1,array,300 ! 定义正弦波的“基础加速度”, 500时间点
- *dim,baseace2,array,300
- *do,i,1,300
- baseace1(i)=0.1*sin(i*0.2)
- baseace2(i)=0.1*sin(i*0.2-31)
- *enddo
- ANTYPE,TRANS !时程分析
- CSYS,1
- NSEL,S,LOC,X,13000,13000
- CSYS,0
- *do,i,1,300
- f,all,fx,3e4*baseace1(i)
- f,all,fz,3e4*baseace1(i) ! F=m*a
- TIME,i*0.02
- *enddo
- NSUBST,300, , ,1 !500个子步
- OUTRES,ALL,1 !输出每个子步的结果
- LSWRITE
- !在Z方向加行波效应
- ANTYPE,TRANS
- CSYS,1
- NSEL,S,LOC,X,13000,13000
- CSYS,0
- *do,j,1,300
- f,all,fz,3e4*baseace2(j)
- f,all,fx,0
- TIME,(j+301)*0.02
- *enddo
- NSUBST,300, , ,1 !500个子步
- OUTRES,ALL,1 !输出每个子步的结果
- LSWRITE
- LSSOLVE,1,2
复制代码 |