new
set logfile myfilesnn.log
set log on
gen zone cylinder p0 0 0 0 p1 0.0125 0 0 p2 0 0.5 0 p3 0 0 0.0125 size 8 12 8
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1
model mohr
prop young 10e9 poisson 0.25 dens 2500 cohe 1e4 ten 1e5 fric 35
set gravity 0 0 -10
apply zvel -.01 range z 19.9 20.1
fix z range z -.1 .1
def ax_stress
f_accum=0.0
pnt = gp_head
loop while pnt#null
if gp_zpos(pnt)<0.1 then
f_accum=f_accum +gp_zfunbal(pnt)
endif
pnt=gp_next(pnt)
endloop
ax_stress=f_accum/900
end
ax_stress
def ax_strain
str1=0
n=1.0
array arr(6)
p_z=zone_head
loop while p_z#null
if z_zcen(p_z)<0.6 then
oo=z_fsi(p_z,arr)
str1=str1-arr(3)
endif
p_z=z_next(p_z)
n=n+1.0
endloop
ax_strain=str1/n
end
ax_strain
table 1 0,0 0,0
def get_plast
shearnow=1
tensionnow=2
shearpast=4
tensionpast=8
v1=0
v2=0
v4=0
v8=0
p_z=zone_head
loop while p_z#null
pid=z_id(p_z)
if z_state(p_z,0)=shearnow then
i=table_size(1)
loop n(1,i)
y=table(1,n)
if pid#y then
v1=v1+1
table(1,n)=pid
endif
endloop
endif
if z_state(p_z,0)=tensionnow then
i=table_size(1)
loop n(1,i)
y=table(1,n)
if pid#y then
v2=v2+1
table(1,n)=pid
endif
endloop
endif
if z_state(p_z,0)=shearpast then
i=table_size(1)
loop n(1,i)
y=table(1,n)
if pid#y then
v4=v4+1
table(1,n)=pid
endif
endloop
endif
if z_state(p_z,0)=tensionpast then
i=table_size(1)
loop n(1,i)
y=table(1,n)
if pid#y then
v8=v8+1
table(1,n)=pid
endif
endloop
endif
p_z=z_next(p_z)
endloop
end
get_plast
hist n 1
;plot zone colorby state
plot block state
set fishcall 0 get_plast
his unbal
his v1
his v2
his v4
his v8
his ax_stress
his ax_strain
step 430
set log off
;hist write 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 file xyz.txt
- 1
- 2
前往页