本人在模拟直剪试验的过程时想记录整个剪切面的剪力变化情况。
看了一些资料，发现孙书伟编著的《FLAC3D在岩土工程中的应用》一书中155页有一段记录接触面剪力的fish命令。但我不明白这串命令是如何做到记录剪力的。
希望大神能帮帮我！命令如下：
def ax_str
str = 0
sdis1=0
n1=0
i_pnt = i_head
node_pnt = i_node_head(i_pnt)
loop while node_pnt # null
n1=n1+1
str = str + in_sstr(node_pnt,1)*in_area(node_pnt)
sdis1=sdis1+in_sdisp(node_pnt,1)
node_pnt = in_next(node_pnt)
endloop
ax_str = str/0.085
sdis1=sdis1/n1
end
hist ns 1
hist ax_str sdis1
