Motion stops after physical time 1s

[Xingyuan]

I am simulating a plate, which is moving forced by a following force  P=P_max*sin(2*pi*frequency*time). I found it strange that no matter how large or small the frequency is  the plate does not move until physical time 1second and the calculation is very very slow after then. The motion of the middle point of the plate e.g. with force's frequency=5Hz can be seen in the attached figure. What is the problem? 

Code: [Select]

FEAP * * moving plate
  0 0 0 3 3 8

PARAmeter
  nx=16
  ny=2
  nz=1
  x1=18.e-3
  y1=4.e-3
  z1=0.0
  x2=48.e-3
  y2=4.5e-3
  z2=2.e-3

BLOCK
  CARTesian nx ny nz 0 0 1 10
    1 x1 y1 z1
    2 x2 y1 z1
    3 x2 y2 z1
    4 x1 y2 z1
    5 x1 y1 z2
    6 x2 y1 z2
    7 x2 y2 z2
    8 x1 y2 z2

BLOCK
  CARTesian nx nz 0 0 2 10
    1 x1 y2 z1
    2 x2 y2 z1
    3 x2 y2 z2
    4 x1 y2 z2

EBOUnd
 1 x1 1 1 1
 1 x2 1 1 1
 3 z1 0 0 1
 3 z2 0 0 1

MATErial,1
  SOLId
  finite mixed
  enhanced
  ELAStic STVK 169.e8 0.3
  DENSity data       8000 !rho

MATErial,2 
  PRESsure
  LOAD 1 1
  follower

END

TIE

batch
  prop,,1
  dt,,0.005
  tplot
  loop,,1000
    time
    loop,newton,300
      utan,line,1        ! Unsymmetric solver -- needed for follower forces
    next
!    plot mesh
!    plot,load,,1
!    plot bound
    plot deformed
    plot,cont,2
  next
END
1 1 0 1 0 0 1.e5 2*3.1415*10
disp,,2,(x1+x2)/2,y1,(z1+z2)/2

inte
stop

[Prof. R.L. Taylor]

There is a bug in the 3-d follower loading (presld.f in the element/frame directory).  Try the attached file to see if it fixes your problem.

You stop your prop load at time = 1 so after that nothing changes.

[Xingyuan]

Thank you Prof. Taylor.
I have tried the attached file. However, during the compiling it reports that the variable 'gdeffl' and 'pstyp' do not have a type.  Maybe my version of FEAP is too old. It is version 7.5. So how can I modify the routine(s) in order to solve the problem in presld.f?

There is a bug in the 3-d follower loading (presld.f in the element/frame directory).  Try the attached file to see if it fixes your problem.

You stop your prop load at time = 1 so after that nothing changes.

[Prof. R.L. Taylor]

We do not have any ver 7.5 that runs on our machines.  To many changes dealing with pointers, etc., since introduction of 64 bit machines.  You will have to debug yourself.