救助：三轴模拟求轴向力的办法

everillusion

三轴试验时轴向的是速度边界，求轴向应力的方法是下面这个函数
def ax_str
str = 0
pnt = gp_head
loop while pnt # null
if gp_ypos(pnt) < 0.1
str = str + gp_yfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
ax_str = str / pi ; cylinder radius = 1
end

我有两个地方不太明白，网格的定义是Y向的单元长度是4/30=0.13，那么gp_ypos(pnt) < 0.1是不是就是用来限制下面的求和时用到的节点都是端面上的点（也就是y=0）的呢
第二个问题，为什么求端面上的力的时候用的是不平衡力总和（yfunbal)，请高手们指教了，谢谢

everillusion

自己顶下，谢谢了