贡献我计算巷道的流变命令流，请教高手看一下为啥最后曲线出不来？哪里的原因？
请教高手看一下我计算的巷道流变，为啥最后曲线出不来？哪里的原因？是不是加载出现问题了？？？？同时希望大家将看完的程序改正后可以尊重劳动果实，互相学习！！
rest aaa.sav
config fluid
def constants
c_perm= 1e-10
c_biotc = 1.
c_biotm = 4.e9* 100.
c_bulk= 5.e8
c_shear = 2.e8
comod = c_bulk + 4. * c_shear / 3.
storage = 1. / c_biotm + c_biotc * c_biotc / comod
cv = c_perm / storage
c_bulku = c_bulk + (c_biotc^2) * c_biotm
c_biotma= 1. / storage
hh = 27.5
bt = cv / (hh * hh)
pi2 = pi * .5
py = 1e5
sig0 = -24000*500
origr1= sig0
p0 = py * c_biotc / (comod * storage)
uy0 =py * hh / comod
csig = c_biotc * 2. * c_shear / comod
end
constants
def point
pnt = gp_near(0.,14.0,0.)
ypnt =y_near(17.5,14.5,0.35)
yy = (hh - gp_ypos(pnt))/ hh
yyy = (hh - y_ycen(ypnt)) / hh
end
point
; --- mechanical model ---
model mo
set fluid biot on
prop bu c_bulk sh c_shear fric 30 coh 0.5e6 ten 1.0e6dil 34 range group 1 anygroup 2 any
fix x z
fix y range y -11.9 -12.1
ini density 2400
set grav 0 -100
ini syy 0.
apply syy sig0 range y 15.4 15.6
apply sxx sig0 range y 15.4 15.6
apply szz sig0 range y 15.4 15.6
ini syyorigr1 grad 0 2.4e40
ini sxxorigr1 grad 0 2.4e40
ini szzorigr1 grad 0 2.4e40
;install bolt
sel cable id1 begin-2.5000000-1.31471500.3500000 end -5.3000000-1.31471500.3500000 nseg1
sel cable id2 begin2.5000000-1.31471500.3500000 end 5.3000000-1.31471500.3500000 nseg1
sel cable id3 begin-2.5000000-0.61471500.3500000 end -5.3000000-0.6147150 0.3500000 nseg1
sel cable id4 begin2.5000000-0.61471500.3500000 end 5.3000000-0.61471500.3500000 nseg1
sel cable id5 begin-2.45685400.4624600 0.3500000 end -5.1638641.1780680.3500000 nseg1
sel cable id 6begin2.45685400.4624600 0.3500000 end 5.1638641.1780680.3500000 nseg1
sel cable id 7begin-2.10000001.3564660 0.3500000 end -4.46789602.85081100.3500000 nseg1
sel cable id 8begin2.10000001.3564660 0.3500000 end 4.46789602.85081100.3500000 nseg1
sel cable id 9begin-0.70000002.4000000 0.3500000 end -1.9592420 4.90086200.3500000 nseg1
sel cable id 10 begin0.70000002.4000000 0.3500000 end 1.9592420 4.90086200.3500000 nseg1
sel cable id3 begin-1.40000002.071232 0.3500000 end -4.98195905.55972200.3500000 nseg 7
sel cable id3 begin-4.98195905.55972200.3500000end -6.41474306.9551190 0.3500000 nseg 11
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+10 gr_c .12E+09 range cid 1,1
sel cable pro emod .200E+12xcarea.18E-03gr_per1.000 yten .16E+06 gr_k 1.000 gr_c 1.000 range cid 2,7
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+07 gr_c .12E+06 range cid 8,18
sel cable id4 begin1.40000002.071232 0.3500000 end 4.98195905.55972200.3500000 nseg 7
sel cable id4 begin4.98195905.55972200.3500000end 6.41474306.9551190 0.3500000 nseg 11
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+10 gr_c .12E+09 range cid 1,1
sel cable pro emod .200E+12xcarea.18E-03gr_per1.000 yten .16E+06 gr_k 1.000 gr_c 1.000 range cid 2,7
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+07 gr_c .12E+06 range cid 8,18
sel cable id5 begin0.00000002.5000000 0.3500000 end 0.00000007.5000000 0.3500000 nseg 7
sel cable id5 begin0.00000007.5000000 0.3500000end 0.00000009.50000000.3500000 nseg 11
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+10 gr_c .12E+09 range cid 1,1
sel cable pro emod .200E+12xcarea.18E-03gr_per1.000 yten .16E+06 gr_k 1.000 gr_c 1.000 range cid 2,7
sel cable pro emod .200E+12xcarea.18E-03yten .16E+06 gr_k .18E+07 gr_c .12E+06 range cid 8,18
; --- fluid flow model ---
model fl_iso
prop perm c_perm biot_c c_biotc
ini biot_mod c_biotm
ini pp 0
; --- fish functions ---
def pp10
pp10= gp_pp(pnt) / py
ft = fltime
c_syy = y_syy(ypnt)/ sig0
c_esyy= (y_syy(ypnt) + y_pp(ypnt)) / sig0
c_uy= gp_ydisp(pnt) / uy0
end
def ppsol
tt = bt * fltime
val = 0.0
valy = 0.0
valu = 0.0
loop m (0,20)
mm = pi2 * (2. * m + 1.)
c_e = exp(-mm*mm*tt)/mm
val = val + sin(mm*yy)*c_e
valu = valu + cos(mm*yy)*c_e/mm
end_loop
ppsol= val*2.*p0/py
ydsol= valu*2.*c_biotc*p0/py + yy - 1.
end
; --- first establish undrained response ---
set fluid off mech on
prop bu c_bulku
ini biot_mod 0.0
solve
ini pp p0
save ucons_und.sav
; --- drained response ---
fix pp 0 range y 15.4 15.6
def my_solve
loop ii (1,10)
c_age = 2e5*ii
command
prop bu c_bulk
ini biot_mod c_biotma
set fluid on mech off
solve age c_age
set fluid off mech on
ini biot_mod 0.0
solve ratio 2e-5
end_command
ytable(1,ii) = pp10
xtable(1,ii) = ft
ytable(2,ii) = ppsol
xtable(2,ii) = ft
ytable(3,ii) = c_uz
xtable(3,ii) = ft
ytable(4,ii) = ydsol
xtable(4,ii) = ft
ytable(5,ii) = c_syy
xtable(5,ii) = ft
ytable(6,ii) = c_esyy
xtable(6,ii) = ft
end_loop
end
my_solve
save ucons.sav
plot create uz
plot set title text
Vertical displacement versus time at column mid-height
plot add table 3 4 line
plot show 这是结果图！没有看到曲线 可以把表格print出来，看看有没有数据。 竖向的位移是有的啊，应力值也有，还请指教。 难道和X轴重合了，这个问题好像出现过
页:
[1]