下面程序,点0,0,0 与点0,1,0 并不是同一个单元,为什么还可以表示其应力应变图,欢迎大家讨论
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 mohr
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
fix x y z range y -.1 .1
fix x y z range y 1.9 2.1
ini yvel 1e-7 range y -.1 .1
ini yvel -1e-7 range y 1.9 2.1
ini pp 1e5
hist gp ydisp 0,0,0 ;采样记录座标[0,0,0]处节点y方向位移
hist zone syy 0,1,0 ;采样记录座标[0,1,0]处单元体yy方向应力
hist zone syy 1,1,0 ;采样记录座标[1,1,0]处单元体yy方向应力
step 3000
plo his -2 vs 1