求解一段FISH函数

shpp14

n
title Structure_dynamic_analysis_lakewater
config dyn
sel pile id=1 beg 0 0 0 end 0 0 1
sel pile prop dens 2400 &
              Emod 1.0e10 Nu 0.3 XCArea 0.3 &
              XCJ 0.16375 XCIy 0.00625 XCIz 0.01575 &
              Per 2.8 &
              CS_sK 1.3e11 CS_sCoh 0.0 CS_sFric 10.0 &
              CS_nK 1.3e11 CS_nCoh 0.0 CS_nFric 0.0 CS_nGap off

def f1
  whilestepping
  f0=10000*sin(10*dytime)
  np = nd_head
  loop while np # null
    if nd_pos(np,1,3)=1
      nd_apply(np,1)=f0
    endif
    np = nd_next(np)
  endloop
end
sel node fix x y z xr yr zr ran id=1
sel set damp combined
这是FLAC/FLAC3D基础与工程实例中的一个例子，在桩顶部施加周期荷载作用，现在看不懂当中那段FISH函数，求指教，谢谢。