!ע⣺λƱһҪ
!ԪĿΪ50100200ʱõĿλһģԵԪĿλӦûжӰ졣

/config,nres,2000
/PREP7                          !ģ
ET,1,BEAM3                      !Ԫ
MP,EX,1,207E9                    !ģ
MP,DENS,1,2000                 !ܶ
MP,NUXY,1,0.3
R,1,0.1,0.0001/12,0.1        !Ծࡢ
NELEM=100
LBEAM=32                      !
deltL=LBEAM/NELEM
f=1000
V=120                          !ƶٶ
NNODE=NELEM+1
deltt=deltL/V*3.6
*do,i,1,NNODE,1
N,i,(i-1)*deltL                  !ڵΪ0.3
*enddo
*do,i,1,NELEM
E,i,i+1
*enddo
finish


/solu
d,1,all
d,NNODE,all
antype,trans 

*do,ii,1,nnode+8,1
outres,all,all 					!мӲн
time,ii*deltt               		!ÿزʱ
nsubst,5                    	!ÿزӲ
*if,ii,lt,8,then
  fdele,all,all               		!ɾеļ
  nsel,s,loc,x,(ii-1)*deltl       		!ѡýڵ
  f,all,FY,-f				 	!ʩӼк
  allsel                        	!ѡ
  solve

*elseif,ii,lt,nnode
  fdele,all,all
  nsel,s,loc,x,(ii-1)*deltl
  nsel,a,loc,x,(ii-7-1)*deltl
  f,all,fy,-f
  allsel
  solve

*else
  fdele,all,all
  nsel,s,loc,x,(ii-7-1)*deltl
  f,all,fy,-f
  allsel
  solve

*endif
*enddo
finish

/post26
NSOL,2,nelem/2+1,U,Y,UY_2       !ڶΪUY_2,ֵΪ101ŽڵYλ
xvar,1                   !ʱΪ
plvar,2                  !λͼ
deriv,3,2,1,,,,,1        !λƶʱ󵼣õٶȣΪ3
plvar,3                  !ٶͼ
fini
