*dim,frequence,,26,,
*dim,amplitude,,26,3,
frequence(1)=10,20,25,31.5,40,50,63,80,100,125
frequence(11)=160,200,250,315,400,500,630,800,1000
frequence(20)=1250,1600,2000,2500,3150,4000,5000
F,node(0.075,0.176,5.5),FY,-170000,
*do,jj,1,26 !!!循环表示
!!以下部分用于进行谐响应计算
!!!!!!!!!!
/SOL !*谐响应分析
ANTYPE,3
HROPT,FULL
HROUT,ON
LUMPM,0
EQSLV,FRONT,0,
PSTRES,0
HARFRQ,frequence(jj),frequence(jj),
NSUBST,1,
KBC,1
OUTRES,ALL,ALL,
solve
fini
/post26
!!!!!!!!!!!!!!!!!!!
nsol,2,node(0.075,0.176,5.5),u,y,UY_2
vget,amplitude(jj,1),1
vget,amplitude(jj,2),2,,0 !!获得计算结果的实部
vget,amplitude(jj,3),2,,1 !!获得计算结果的虚部
fini
*enddo
aaa
[ 本帖最后由 austin 于 2008-10-15 22:41 编辑 ] |