T-1.txt

节选段落一：
n
restore sanzhou
def modulus
jianqie=young/(2*(1+v))
tiji=young/(3*(1-2*v))
end

set @young=4035e6
set @v=0.2575
@modulus

set grav 0 0 -9.81

model mohr
prop density 2652 shear @jianqie bulk @tiji coh 5.52e6 fric 41.3 range group 2

fix x y z range y -0.001 0.001
fix x y z range


节选段落二：
)
endloop
axi_stress = f_accum / (3.141592653*0.025*0.025)
end
def lat_stress
g_accum=0.0
pnt = gp_head
loop while pnt # null
if gp_xpos(pnt) <0.00125 then
g_accum = g_accum + gp_xfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
lat_stress = g_accum / (3.141592653*0.025*2*0.1)
end
hist n 1


节选段落三：
hist @axi_stress
hist @lat_stress
hist gp ydis 0 0.1 0
hist gp ydis -0.025 0.1 0
hist gp xdisp -0.025 0.05 0
hist gp xdisp 0.025 0.05 0

plot hist -1 vs -4
;plot hist -2 vs -5step 1500
step 10000