节点不平衡力监测

zsq

1.ad1 = gp_near(0,2,0)，0 2 0这点已经精确，最近意义是什么？2.监测四个角点意义是什么，请大家帮忙看看
new
gen zone brick size 1 2 1
model mohr
prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10
fix x y z range y -0.1 0.1
apply yvel -1e-5 range y 1.9 2.1
plot set rotation 0 0 45
plot block group
def get_ad
ad1 = gp_near(0,2,0)
ad2 = gp_near(1,2,0)
ad3 = gp_near(0,2,1)
ad4 = gp_near(1,2,1)
end
get_ad
def load
load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)
end
hist load
hist gp ydis 0,2,0
step 1000
plot his 1 vs -2