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
学习了,谢谢楼主
页:
[1]