地震波

对于地震波的输入，可以把荷载记录做成文件，利用apdl的读取功能读入倒数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。有用到的朋友自己作一下把。
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句：
/SOLU
ANTYPE,trans
！求解其自己选了
TM_START=0.01
TM_END=15.00
TM_INCR=0.01
*DO,TM,TM_START,TM_END,TM_INCR
TIME,tm
alpha,
BETAD,
ACEL,acex(tm),acey(tm),acez(tm)
SOLVE
*ENDDO
fini
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句：
/SOLU
ANTYPE,trans
！求解其自己选了
TM_START=0.01
TM_END=15.00
TM_INCR=0.01
*DO,TM,TM_START,TM_END,TM_INCR
TIME,tm
alpha,
BETAD,
ACEL,acex(tm),acey(tm),acez(tm)
SOLVE
*ENDDO
fini
*************************************************************
转贴]计算地震反映谱分析的命令流
FINISH
/solu
ANTYPE,MODAL
MODOPT,LANB,30
LUMPM,ON
SOLVE
FINISH

/SOLU
ANTYPE,SPECTR
SPOPT,SPRS,30,YES
SVTYPE,2,.107*9.8!加速度反应谱,0.107=Ci*Cz*Kh查规范
SED,1,0,0
!I类场地土,0.05 IS THE DAMPING COEFFICIENT 查规范
!FREQ,0.16667,0.2,0.2222,0.25,0.2857,0.3333,0.4,0.5,0.66667,
!FREQ,0.7143,0.83333,0.9091,1,1.25,1.66667,2.5,3.33333,10,
!SV,0.05,0.3409,0.3879,0.4114,0.4349,0.4585,0.4820,0.5055,0.5290,0.5525
!SV,0.05,0.5879,0.6754,0.7305,0.7959,0.9729,1.209,1.815,2.352,2.352
!II类场地土
FREQ,0.1,0.2,0.4,0.5,0.

66667,0.83333,0.9091,1,1.1111
FREQ,1.25,2,2.5,3.333,4,5,10,20,100
SV,0.05,0.3,0.3,0.3,0.35055,0.4647,0.5783,0.6298,0.69145,0.76666
SV,0.05,0.8605,1.3639,1.6972,1.9345,2.25,2.25,2.25,2.25,1.725,1.305
SOLVE
FINISH

/SOLU
ANTYPE,MODAL
EXPASS,ON
MXPAND,30,,,YES,0.01
SOLVE
FINISH

/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SRSS,0.01,DISP
OUTRES,ALL,ALL
SOLVE
FINISH
/post1
/input,spec,mcom
lcwrite,11
lcase,11
prrsol,
PRNLD, , ,
ETABLE,FXI,SMISC,1
ETABLE,MZI,SMISC,2
ETABLE,MYI,SMISC,3
ETABLE,MXI,SMISC,4
ETABLE,FYI,SMISC,5
ETABLE,FZI,SMISC,6
ETABLE,FXJ,SMISC,14
ETABLE,MZJ,SMISC,15
ETABLE,MYJ,SMISC,16
ETABLE,MXJ,SMISC,17
ETABLE,FYJ,SMISC,18
ETABLE,FZJ,SMISC,19
finish