输出接触面部分区域的应力

sinlk

如题，想输出接触面1纵坐标小于-0.4面上的stress，我写了下面的一个小程序，运行通过，只是没有输出22.txt这个文件，请问这是怎么回事呢？请大家帮帮忙，急！谢谢了
def get_sstress
  pnt=gp_head
  if gp_zpos(pnt)<-0.4 then
    command
      set log on
      set logfile 22.txt
      print int 1 stress
      set log off
    endcommand
    pnt=gp_next(pnt)
  endif
end
get_sstress