flac3D主应力提取及记录

jouleo

在flac3D的后处理中，有时候需要对单元应力提取并查看，除了可以采用set log on来实现外，可以用下面提供的简单的fish函数来实现：
def zonesig
  array arr(20000) brr(20000) nrr(20000)
  znt=zone_head
  ii=1
  loop while znt # null
    arr(ii)=string(z_sig1(znt))  ；提取第一应力
    brr(ii)=string(z_sig3(znt))  ；提取第三应力
    nrr(ii)=string(z_id(znt))    ；提取单元编号
    znt=z_next(znt)
    ii=ii+1
  endloop
jj=ii-1
end   
zonesig
def recordsig
  status=open('sig1.txt',1,1)
  status=write(arr,jj)
  status=close
  status=open('sig3.txt',1,1)
  status=write(brr,jj)
  status=close
  status=open('zid.txt',1,1)
  status=write(nrr,jj)
  status=close
end
recordsig

772559732

学习了，谢谢楼主