ansys施加正弦载荷的例子.doc

ansys施加正弦载荷的例子
/PREP7
ET,1,SHELL63
R,1,0.1

MP,EX,1,2e11
MP,NUXY,1,0.3
MP,DENS,1,7800

BLC4,,,10,2
AATT,1,1,1,0
AESIZE,ALL,0.8
MSHAPE,0,2D
MSHKEY,0
AMESH,ALL

FINISH
SAVE

/SOLU
ANTYPE,4
TRNOPT,FULL

LSEL,S,LOC,X,0
DL,ALL,,ALL,0
ALLSEL,ALL

FZ_KP=KP(10,2,0)
DT=0.1
*DO,J,1,40
   FK,FZ_KP,FZ,3000*SIN(2*3.14*DT*J)
   TIME,DT*J
   NSUBST,1,,,1
   KBC,0
   OUTRES,ALL,ALL
   SOLVE
*ENDDO

TIME,40*DT+2
NSUBST,20,,,1
KBC,0
SOLVE
FINISH

/POST26
N_FZ=NODE(10,2,0)
NSOL,2,N_FZ,U,Z,UZ
PLVAR,2

比如我要定义一个正弦载荷，如下：force=-1000+2000*sin(6.28*10*t),t是时间，分析的时间是0——5秒，并分成500份瞬态加载
*do,t,1,5,0.01
force=-1000+2000*sin(6.28*10*t)
f,node,fx,force
time,i
solve
*enddo

/SOLU
ANTYPE,4
TRNOPT,FULL
ALLSEL,ALL

DT=0.001
*DO,J,1,40
   F, NTOP,FZ,SIN(2*3.14*DT*J)
   TIME,DT*J
   NSUBST,1,,,1
   KBC,0
   OUTRES,ALL,ALL
   SOLVE
*ENDDO

TIME,4000*DT+2
NSUBST,20,,,1
KBC,0
SOLVE
FINISH

/SOLU
ANTYPE,4
TRNOPT,FULL
FZ_KP=KP(10,2,0)
DT=0.001
*DO,J,1,40
   FK,FZ_KP,FZ,3000*SIN(2*3.14*DT*J)
   TIME,DT*J
   NSUBST,1,,,1
   KBC,0
   OUTRES,ALL,ALL
   SOLVE
*ENDDO

TIME,4000*DT+2
NSUBST,20,,,1
KBC,0
SOLVE
FINISH

/POST26
N_FZ=NODE(10,2,0)
NSOL,2,N_FZ,U,Z,UZ
PLVAR,2