[分享]地震波的输入和求解
!对于地震波的输入，可以把荷载记录做成文件，利用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
在地震输入中怎样考虑行波效应呢？ 我正在用ANSYS做抗震分析，有用！ 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 感谢阿！ 我也要做地震分析 谢谢拉！ 谢谢！ 如何获得位移地震波 有没有位移地震波啊 我也正在考虑行波效应，谢谢了！ 谢谢大侠的代码！！！！！！！！！！！！！
页:
[1]