彭文斌的书中例8-3中:
new
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1
model ss ;应变硬化/软化模型
; --- 材料参数 ---
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
prop ttab 1
table 1 0,2e5 2e-5,0
; --- 边界及初始条件 ---
fix x y z range y -.1 .1
fix x y z range y 1.9 2.1
ini yvel -1e-8 range y -.1 .1
ini yvel 1e-8 range y 1.9 2.1
def ax_str
str = 0
pnt = gp_head
loop while pnt # null
if gp_ypos(pnt) < 0.1 then ; Y轴方向<0.1坐标的所有节点的力
str = str + gp_yfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
ax_str = str / pi ; cylinder radius = 1,应力=力/面积
end
; --- 采样记录 ---
hist n 1
hist gp ydisp 0,0,0
hist ax_str
hist gp xdisp 1,1,0
; --- 试验 ---
step 1500
; --- 结果 ---
plot hist 2 vs -1 xlabel '轴向位移 ' ylabel '轴向应力';axial stress vs axial disp.
pause
plot hist 3 vs -1 xlabe '轴向位移 ' ylabel '径向位移'; circumferential disp. vs axial disp.
ret
请问最后plot hist 后面的编号 2,3,-1究竟是怎么对应的啊??有时候编号还是负数的?小弟是初学者,还望高人指点一下!!