def setup
IO_READ = 0
IO_WRITE = 1
IO_FISH = 0
IO_ASCII = 1
FlacFile='jdstress.dat'
end
setup
def jdl
status = open(FlacFile, IO_WRITE, IO_ASCII)
array rface(1)
array zhig(1)
array lface(1)
array face(1)
array rkou(1)
array ckou(1)
;--------------右侧面---------------
loop i (1,16)
loop j (1,36)
z=1675.88-(i-1)*20
x=xtable(5,j)
y=ytable(5,j)
p_z=z_near(x,y,z)
sx=z_sxx(p_z)
sy=z_syy(p_z)
sz=z_szz(p_z)
xy=z_sxy(p_z)
xz=z_sxz(p_z)
yz=z_syz(p_z)
rface(1)=string(sx)+' '+string(sy)+' '+string(sz)+' '+string(xy)+' '+string(xz)+' '+string(yz)
if j=1 then
zhig(1)='右侧面高度'+string(z)+'\n'
status=write(zhig,1)
status=write(rface,1)
else
status=write(rface,1)
endif
endloop
endloop
;------------左侧面---------------
loop i (1,16)
loop j (1,36)
z=1675.88-(i-1)*20
x=xtable(6,j)
y=ytable(6,j)
p_z=z_near(x,y,z)
sx=z_sxx(p_z)
sy=z_syy(p_z)
sz=z_szz(p_z)
xy=z_sig1(p_z)
xz=z_sig2(p_z)
yz=z_sig3(p_z)
lface(1)=string(sx)+' '+string(sy)+' '+string(sz)+' '+string(xy)+' '+string(xz)+' '+string(yz)
if j=1 then
zhig(1)='左侧面高度'+string(z)+'\n'
status=write(zhig,1)
status=write(lface,1)
else
status=write(lface,1)
endif
endloop
endloop
;--------------上表面--------------------
loop i (1,16)
loop j (1,36)
z=1675.88
x=xtable(7,j+(i-1)*36)
y=ytable(7,j+(i-1)*36)
p_z=z_near(x,y,z)
sx=z_sxx(p_z)
sy=z_syy(p_z)
sz=z_szz(p_z)
xy=z_sig1(p_z)
xz=z_sig2(p_z)
yz=z_sig3(p_z)
face(1)=string(sx)+' '+string(sy)+' '+string(sz)+' '+string(xy)+' '+string(xz)+' '+string(yz)
if j=1 then
zhig(1)='上表面第'+string(i)+'行'+'\n'
status=write(zhig,1)
status=write(face,1)
else
status=write(face,1)
endif
endloop
endloop
;--------------入洞口-----------------------
loop i (1,16)
loop j (1,16)
z=1675.88-(i-1)*20
x=xtable(7,j*36)
y=ytable(7,j*36)
p_z=z_near(x,y,z)
sx=z_sxx(p_z)
sy=z_syy(p_z)
sz=z_szz(p_z)
xy=z_sig1(p_z)
xz=z_sig2(p_z)
yz=z_sig3(p_z)
rkou(1)=string(sx)+' '+string(sy)+' '+string(sz)+' '+string(xy)+' '+string(xz)+' '+string(yz)
if j=1 then
zhig(1)='入口高度'+string(z)+'\n'
status=write(zhig,1)
status=write(rkou,1)
else
status=write(rkou,1)
endif
endloop
endloop
;--------------出洞口-------------------
loop i (1,16)
loop j (1,16)
z=1675.88-(i-1)*20
x=xtable(7,(j-1)*36+1)
y=ytable(7,(j-1)*36+1)
p_z=z_near(x,y,z)
sx=z_sxx(p_z)
sy=z_syy(p_z)
sz=z_szz(p_z)
xy=z_sig1(p_z)
xz=z_sig2(p_z)
yz=z_sig3(p_z)
ckou(1)=string(sx)+' '+string(sy)+' '+string(sz)+' '+string(xy)+' '+string(xz)+' '+string(yz)
if j=1 then
zhig(1)='出口高度'+string(z)+'\n'
status=write(zhig,1)
status=write(ckou,1)
else
status=write(ckou,1)
endif
endloop
endloop
status=colse
end
jdl
