麻烦哪位大神解释一下下面这段代码是如何求轴向应力的？

佚名

def ax_str
        str = 0
        pnt = gp_head   ;把网格点指针的首个地址（gp_head）赋予变量pnt
        loop while pnt # null
                if gp_ypos(pnt) < 0.1 then   ;当变量pnt所对应的y轴方向坐标小于0.1
                        str = str + gp_yfunbal(pnt)    ;变量str为上次循环中的值与变量pnt所对应的y轴方向上的不平衡力的和
                endif
                pnt = gp_next(pnt)    ; 把网格点指针的下一个地址（gp_next）赋予变量pnt
        endloop
        ax_str = str / pi ; cylinder radius = 1，应力=力/面积
end