大家帮忙看一下

zg15848112854

config flu dynamic
set flu off dyn off
;
def setup
   aa=1.0
   bb=8.0
   rad=0.165
   piledi=-6
   modeldi=-20
end
@setup
gen z radcyl p0 0 0 0 p1 @aa 0 0 p2 0 0 @piledi p3 0 @aa 0 p4 @aa 0 @piledi p5 0 @aa @piledi ...
             p6 @aa @aa 0 p7 @aa @aa @piledi p8 @rad 0 0 p9 0 @rad 0 p10 @rad 0 @piledi ...
             p11 0 @rad @piledi size 5 6 6 3 r 1 1 1 1.14  

gen z radcyl p0 0 0 @piledi p1 @aa 0 @piledi p2 0 0 @modeldi p3 0 @aa @piledi p4 @aa 0 @modeldi p5 0 @aa @modeldi ...
             p6 @aa @aa @piledi p7 @aa @aa @modeldi p8 @rad 0 @piledi p9 0 @rad @piledi p10 @rad 0 @modeldi ...
             p11 0 @rad @modeldi size 5 14 6 3 r 1 1 1 1 fill


gen z radtun p0 0 0 0 p1 @bb 0 0 p2 0 0 @modeldi p3 0 @bb 0 p4 @bb 0 @modeldi p5 0 @bb @modeldi ...
             p6 @bb @bb 0 p7 @bb @bb @modeldi p8 @aa 0 0 p9 0 @aa 0 p10 @aa 0 @modeldi p11 0 @aa @modeldi ...
             p12 @aa @aa 0 p13 @aa @aa @modeldi size 5 20 5 3 r 1 1 1 1.4
gro clay ran z 0 -15
gro sand ran z -15 -20
gen z refl dip 90 dd 90
gen z refl nor 0 -1 0 ori 0 0 0

;?¨???￡¤??
def int_pra
    pos_s=-5.9
    pos_x=-6.1
    rad_=0.265
    radn_=0.065
end
@int_pra
int 1 face ran cylin end1 0 0 0 end2 0 0 @pos_x radius @rad_ ...
               cylin end1 0 0 0 end2 0 0 @pos_x radius  @radn_ not
int 2 face ran cylin end1 0 0 @pos_s end2 0 0 @pos_x radius @rad_  



;??§¨￠?é?á?

gen z cylin p0 0 0 7.15 p1 @rad 0 7.15 p2 0 0 1 p3 0 @rad 7.15 p4 @rad 0 1 p5 0 @rad 1 ...
            size 5 15 6 r 1 1 1 gro AA
gen z refl dip 90 dd 90  ran gro AA
gen z refl nor 0 -1 0 ori 0 0 0 ran gro AA
gro pile1 ran z 1 7.15

ini z add -7 ran gro pile1


attach face

;??§¨￠?é3D¨??§
gen z radcyl p0 0 0 0.15 p1 1 0 0.15 p2 0 0 0 p3 0 1 0.15 p4 1 0 0 ...
             p5 0 1 0 p6 1 1 0.15 p7 1 1 0 p8 @rad 0 0.15 p9 0 @rad 0.15 ...
             p10 @rad 0 0 p11 0 @rad 0 size 5 2 6 3  r 1 1 1 1.14 gro cap1
gen z reflect dip 90 dd 0 ran gro cap1
gen z reflect normal -1 0 0 ori 0 0 0.15 ran gro cap1
int 1 face ran cylin end1 0 0 0.15 end2 0 0 0 radius @rad_ ...
               cylin end1 0 0 0.15 end2 0 0 0 radius  @radn_ not
int 2 face ran cylin end1 0 0 0.1 end2 0 0 -0.1 radius @rad_  

int 3 face ran x -1 1 y -1 1 z 0.15
gen z b p0 0 0 0.3 p1 1 0 0.3 p2 0 0 0.15 p3 0 1 0.3 p4 1 0 0.15 ...
        p5 0 1 0.15 p6 1 1 0.3 p7 1 1 0.15 size 5 2 5 r 1 1 1  gro cap2
gen z refl dip 90 dd 90 ran  gro cap2
gen z refl nor 0 -1 0 ori 0 0 0.3 ran gro cap2
gro cap ran gro cap1
gro cap ran gro cap2


;material property

model mohr ran gro clay
model mohr ran gro sand
model elas ran gro pile1
model elas ran gro cap
prop bulk 4.1667e7 shear 3e6 coh 3000 fric 35 range group clay
prop bulk 4.1667e7 shear 3e6 coh 3000 fric 45 range group sand
prop bulk 3e8 shear 3e8 range group pile1
prop bulk 3e8 shear 3e8 range group cap
               
int 1 prop kn 1e8 ks 1e8 coh 3000 fri 28
int 2 prop kn 1e8 ks 1e8 coh 3000 fri 20
int 3 prop kn 1e8 ks 1e8 coh 3000 fri 0

ini dens 2200 ran gro clay
ini dens 1800 ran gro sand
ini dens 2800 ran gro pile1
ini dens 2800 ran gro cap

;k0=u/(1-u)
ini szz 0 grad 0 0 17000 ran z 0 -20
ini sxx 0 grad 0 0 7293 ran z 0 -20
ini syy 0 grad 0 0 7293 ran z 0 -20
;
fix x ran x -8
fix x ran x 8
fix y ran y -8
fix y ran y 8
fix z ran z -20
water density 1000
set grav 0 0 -9.8
water table ori 0 0 0 nor 0 0 -1

set mech ratio 1e-5
hist id 1 unbal

solve



;
ini state 0
ini xdis 0 ydis 0 zdis 0 xv 0 yv 0 zv 0
model fl_iso
set flu on
pro perm 1e-13 poros 0.5 biot_c 1.0 ran gro clay
pro perm 1e-17 poros 0.5 biot_c 1.0 ran gro sand
ini pp 1.0e4 grad 0 0 -1.0e4
ini fmod 2e9 fdens 1000
set flu pcut on
fix pp 0 range z -0.1 0.1
fix pp 0 ran x -8
fix pp 0 ran  x 8
fix pp 0 ran y -8
fix pp 0 ran y 8

;
ini state 0
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
set dyn on flu off
app nquiet squiet ran z 0.3
table 1 read dynload.txt
app szz -2e5 hist table 1 ran z 0.3
table 2 read dynload1.txt
app syz -2e5 hist table 2 ran z 0.3   
ini pp 30e3 grad 0 0 -10e3 ran z -20
app nquiet dquiet squiet ran z -20
set dyn damp rayleigh 1 22.8
;
ini state 0
ini xdis 0 ydis 0 zdis 0 xv 0 yv 0 zv 0

config thermal
set thermal on
model th_iso
set thermal off
pro conduct 1.8 spec_heat 1000  ran gro clay
pro conduct 1.21 spec_heat 1406  ran gro sand
pro conduct 1.74 spec_heat 920  ran gro pile1
pro conduct 1.74 spec_heat 920  ran gro cap


def th_add
   loop n(15,30)
   ini_th=5+n
     command
        ini  t @ini_th ran gro pile1
     endcommand
   endloop
end
@th_add
ini t 5 ran gro clay
ini t 5 ran gro sand
app vsource 30 ran gro pile1
app vsource 30 interior ran gro pile1

solve