[求助]FLAC3D单轴模拟问题（圆柱壁和底端固定），看代码有什么问题吗？？ [复制链接]

$De14  
new F<H[-k*t/  
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 RrB)u?  
gen zone reflect norm 1,0,0 J8IdQ:4^l  
gen zone reflect norm 0,0,1 HGIPz{/5U  
P5-1z&9O  
model ss 
{S[+hUl  
prop bulk 1.19e6 shear 1.1e6 dens 1500 |[qq $  
prop coh 2.72e10 fric 44 ten 2e5 ctab 1 ftab 2 %c:v70*h=  
table 1 0 100e3 0.05 50e3 0.1 25e3 1 25e3 Z1Y/2MVSb  
table 2 0 23 0.05 20 0.1 18 1 18 OI/m_x
x@j  
nIc:<w]  
fix x y z range y -.1 .1 qM}Uk3N0  
fix x y z range y 1.9 2.1 X)6}<A  
fix x y z range cylinder end1 0,0,0 end2 0,2,0 radius 1.1 cyl end1 0,0,0 end2 0,2,0 radius 0.9 not ;r<(n3"F  
ini yvel 1e-7 range y -.1 .1 jT/}5\  
Zcst$Aro  
hist gp ydisp 0,0,0 ; ;J[ed>v;3  
hist zone syy 0,0,0 ;  =ie8{j2:  
hist zone syy 0,1,0 ; /q[5-96c  
step 1000 (OS -v~{r@  
d^b(Uo=$  
plot hist -2 vs 1 /6S% h-#\  
hist write 1 2 file 单压.txt

，，，，，，，，，，，，，，，，，，，，