jouleo 发表于 2018-3-14 17:28:45

flac3D主应力提取及记录

在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 发表于 2022-8-19 22:35:30

学习了,谢谢楼主
页: [1]
查看完整版本: flac3D主应力提取及记录