new
def test0
tim0=clock
end
def test1
tim=(clock-tim0)/100.0
end
gen zone brick size 6 1 20 p1=3.0,0.0,0.0 ratio 0.9 1.0 0.97
gen zone brick size 20 1 20 p0=3.0,0.0,0.0 p1=20.0,0.0,0.0 &
ratio 1.08 1.0 0.97
ini z mul 0.5
model mohr
prop bul 2.e8 shea 1.e8 cohesion 1.e5
prop friction 0. dilation 0. tension 1.e10
fix x range x -.1 .1
fix x y z range z -.1 .1
fix x y z range x 19.9 20.1
fix y
fix x y z range x -.1 3.1 z 9.9 10.1
free x range x 2.9 3.1 z 9.9 10.1
ini zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1
def p_cons
array LoadPoints(50) ; Must increase if grid zones are increased
pdis1 = gp_near(0.,0.,10.)
pdis2 = gp_near(0.,1.,10.)
p_sol = (2. + pi)
setCohes = z_prop(zone_head,’cohesion’)
pnt = gp_head
n =0
xnext = 20.0
loop while pnt # null
if gp_zpos(pnt) > 9.9
if gp_xpos(pnt) < 3.1
n=n+1
LoadPoints(n) = pnt ; ... save addresses
endif
if gp_xpos(pnt) > 3.1 ; ... get 1st gp outside footing
xnext = min(xnext,gp_xpos(pnt))
endif
endif
pnt = gp_next(pnt)
endLoop
nLoads = n
EffectiveWidth = (3.0 + xnext) / 2.0
end
p_cons
def p_load
pload = 0.0
loop n (1,nLoads)
pnt = LoadPoints(n)
pload = pload + gp_zfunbal(pnt)
endLoop
actLoad = pload / (EffectiveWidth * setCohes)
c_disp = -(gp_zdisp(pdis1) + gp_zdisp(pdis2)) / 2.0
p_load = actLoad
end
def p_err
p_err = 100. * (actLoad - p_sol) / p_sol
end
hist n 50
hist p_load
hist p_sol
hist c_disp
hist unbal
test0
plo set pers off
cyc 25000
test1
print tim
print p_err
save pran3.sav
ret

建议你到技术邻提问，里面有专门的专家回答相关问题，提问免费，不亏哟~