- 积分
- 0
- 注册时间
- 2016-1-8
- 仿真币
-
- 最后登录
- 1970-1-1
|
new
config fluid dynamic
set dyn off
gen zone brick size 2 10 2 p0 (0,0,0) p1 (.3,0 ,0) p2 (0,0,-10) p3 (0,.3,0) group soil1
gen zone brick size 2 10 5 p0 (0 ,0.3,0) p1 (0.3,0.3,0) p2 (0,0.3,-10) p3 (0,6,0) p4 (0.3,0.3,-10) &
p5 (0,6,-10) p6 (2, 6, 0) p7 (2 ,6 ,-10) group soil2
gen zone brick size 2 10 5 p0 (0.3,0.3,0) p1 (0.3,0,0) p2 (0.3,0.3,-10) p3 (6,2,0) p4 (0.3,0,-10) &
p5 (6,2,-10) p6 (6, 0, 0) p7 (6 ,0 ,-10) group soil2
gen zone brick size 5 10 5 p0 (0.3,0.3,0) p1 (6,2,0) p2 (0.3,0.3,-10) p3 (2,6,0) p4 (6,2,-10) &
p5 (2,6,-10) p6 (8, 8, 0) p7 (8 ,8 , -10) group soil2
gen zone brick size 2 10 2 p0 (0,6,0) p1 (2,6,0) p2 (0,6,-10) p3 (0,10,0) p4 (2,6,-10) &
p5 (0,10,-10) p6 (2, 10, 0) p7 (2 ,10 , -10) group soil2
gen zone brick size 2 10 6 p0 (6,2,0) p1 (6,0,0) p2 (6,2,-10) p3 (20,2,0) p4 (6,0,-10) &
p5 (20,2,-10) p6 (20, 0, 0) p7 (20 ,0 , -10) group soil2
gen zone brick size 5 10 2 p0 (2,6,0) p1 (8,8,0) p2 (2,6,-10) p3 (2,10,0) p4 (8,8,-10) &
p5 (2,10,-10) p6 (8, 10, 0) p7 (8,10 , -10) group soil2
gen zone brick size 5 10 6 p0 (8,8,0) p1 (6,2,0) p2 (8,8,-10) p3 (20,8,0) p4 (6,2,-10) &
p5 (20,8,-10) p6 (20, 2, 0) p7 (20 ,2 , -10) group soil2
gen zone brick size 6 10 2 p0 (8,8,0) p1 (20,8,0) p2 (8,8,-10) p3 (8,10,0) p4 (20,8,-10) &
p5 (8,10,-10) p6 (20, 10, 0) p7 (20 ,10 , -10) group soil2
gen zone reflec dip 90 dd 0
gen zone reflec dip 90 dd 90
;边界条件
fix x y z ran z -9.9 -10.1
fix x ran x -19.9 -20.1
fix x ran x 19.9 20.1
fix y ran y -9.9 -10.1
fix y ran y 9.9 10.1
; 属性
model mo range z 0 -10.1
model fl_iso ran z 0 -10.1
ini fdens 1000
set grav 10.0
ini pp 0 grad 0 0 -10e3 ran z 0 -10.1
hist unbal
set mech on fluid off
ini fmod 0
prop dens 1040 por 0.66 bu 14.71e6 sh 5.640e6 f 25 c 10e3 perm 1e-13 range z 0 -1
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7 f 30 c 0 perm 1e-8 range z -1 -9
prop dens 1100 por 0.60 bu 14.71e6 sh 5.640e6 f 25 c 15e3 perm 1e-13 range z -9 -10.1
ini szz 0 grad 0 0 17e3 ran z 0 -1 ;soila
ini syy 0 grad 0 0 13.5e3 ran z 0 -1
ini sxx 0 grad 0 0 13.5e3 ran z 0 -1
ini szz 2e3 grad 0 0 19e3 ran z -1 -9 ;soilb
ini syy 1e3 grad 0 0 14.5e3 ran z -1 -9
ini sxx 1e3 grad 0 0 14.5e3 ran z -1 -9
ini szz -16e3 grad 0 0 17e3 ran z -9 -10.1 ;soilc
ini syy -8e3 grad 0 0 13.5e3 ran z -9 -10.1
ini sxx -8e3 grad 0 0 13.5e3 ran z -9 -10.1
solve
save 初始平衡.sav
rest 初始平衡.sav
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
set mech on fluid off
del range group soil1 z -10 0
;生成桩
gen zone brick size 4 13 4 p0 (-.3,-.3,15) p1 (.3,-.3 ,15) p2 (-.3,-.3,2) p3 (-.3,.3,15) group pile2
interface 1 face range x -0.29 -.31 y -.29 -.15 z -10 0
interface 1 face range x -0.29 -.31 y .15 .29 z -10 0
interface 1 face range x 0.29 .31 y -.29 -.15 z -10 0
interface 1 face range x 0.29 .31 y .15 .29 z -10 0
interface 1 face range y -0.29 -.31 x -.29 .29 z -10 0
interface 1 face range y 0.29 .31 x -.29 .29 z -10 0
interface 2 face range x -0.29 -.31 y -.15 .15 z -10 0
interface 2 face range x 0.29 .31 y -.15 .15 z -10 0
;interface 1 face range x -.29 .29 y -.29 .29 z -10.1 -9.9
int 1 prop kn 1e8 ks 1e8 f 20 c 0 range z -1 -9
int 1 prop kn 1e8 ks 1e8 f 15 c 7e3 range z 0 -1
int 1 prop kn 1e8 ks 1e8 f 15 c 10e3 range z -9 -10.1
int 2 prop kn 1e8 ks 1e8 f 20 c 0 range z -1 -9
int 2 prop kn 1e8 ks 1e8 f 15 c 7e3 range z 0 -1
int 2 prop kn 1e8 ks 1e8 f 15 c 10e3 range z -9 -10.1
interface 1 perm = off
interface 2 perm = off
interface 1 effe = off
interface 2 effe = off
ini z add -12 ran gro pile2
fix z ran z -10.1 -9.9
;分组
group soil range group pile2 not
group soila range group soil z 0 -1
group soilb range group soil z -1 -9
group soilc range group soil z -9 -10.1
;桩属性
group pile1 range group pile2 z 2 3
group pile range group pile2 z 2 -10
model elastic ran gro pile
prop bu 1.67e10 sh 7.69e9 ran gro pile
ini dens 2400 ran gro pile ;z 0 2
;桩顶重物
;属性
model elastic ran gro pile1
prop bu 1.67e10 sh 7.69e9 ran gro pile1
ini dens 2400 ran gro pile1
;土属性
model mo range group soila any group soilb any group soilc any
prop dens 1040 por 0.66 bu 14.71e6 sh 5.640e6 f 25 c 10e3 range group soila
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7 f 30 c 0 range group soilb
prop dens 1100 por 0.60 bu 14.71e6 sh 5.640e6 f 25 c 15e3 range group soilc
;pau
;水属性
model fl_iso range group soila any group soilb any group soilc any ;group toushui any
model fl_nul range group pile a group pile1 a
prop perm 1e-13 range group soila
prop perm 1e-8 range group soilb
prop perm 1e-13 range group soilc
ini fdens 1000 ran gro soila any group soilb any group soilc any
ini pp 0 grad 0 0 -10e3 ran z 0 -10.1 gro soila any group soilb any group soilc any
set mech on fluid off
ini fmod 0 ran gro soila any group soilb any group soilc any
ini sat = 0 range group pile any group pile1 any
;pau
def add_top
loop n (1,10)
ini_add = n * 2400
save_file = string(n) + 'G.sav'
command
ini dens ini_add ran gro pile1
set mech rat 1e-5
solve
save save_file
endcommand
endloop
end
add_top
set mech rat 1e-6
solve
save 加桩后初始应力.sav
rest 加桩后初始应力.sav
ini state 0
set dyn on mech on fluid on
set dyn mult on
set fluid time 0
set mech rat 1e-100
set fluid pcut on ;土体不产生负孔压
model finn ran group soilb
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7 f 30 c 0 perm 1e-8 ran group soilb
prop ff_latency=50 ran group soilb
prop ff_switch = 0 ran group soilb
prop ff_c1=0.8 ff_c2=0.79 ran group soilb
prop ff_c3=0.45 ff_c4=0.73 ran group soilb
ini fmod 2e8 ran gro soila any group soilb any group soilc any
;阻尼
ini damp local 0.1571 range group pile a group pile1 a
ini damp local 0.3142 range group soila a group soilc a
ini damp local 0.4713 range group soilb
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
free x range z -10.1 -9.9
;底部速度施加
def setup
freq=5.0
ampl=2
omega = 2.0 * pi * freq
end
setup
def sine_wave
vv = ampl/omega*cos(omega*dytime)
if dytime < 2.0
sine_wave = dytime / 2.0 * vv
else
if dytime < 15.0
sine_wave = vv
else
if dytime <= 17.0
sine_wave = (17.0 - dytime) / 2.0 * vv
endif
endif
endif
if dytime > 17.0
sine_wave = 0.0
endif
end
apply xvel=1.0 hist sine_wave range z -9.9 -10.1
apply ff
call ppr.dat
;历史变量监测
his reset
hist id 1 unbal
hist id 2 dytime
hist id 3 sine_wave
;加速度(底部输入加速度时程)
hist id 111 gp xacc 0 0 -10
hist id 112 gp xacc 0.3 0 -10
hist id 113 gp xacc 1.4 0 -10
hist id 114 gp xacc 2.5 0 -10
hist id 115 gp xacc 3.6 0 -10
hist id 116 gp xacc 4.7 0 -10
hist id 117 gp xacc 5.9 0 -10
hist id 118 gp xacc 8 0 -10
hist id 119 gp xacc 12 0 -10
hist id 120 gp xacc 16 0 -10
def get_ppr_zone
whilestepping
p_z43 = z_near(1 , 0.1, -5.5)
ppr_431 = z_extra(p_z43,2)
ppr_432 = z_extra(p_z43,3)
ppr_433 = z_extra(p_z43,6)
ppr_434 = z_extra(p_z43,7)
p_z44 = z_near(2, 0.1, -5.5)
ppr_441 = z_extra(p_z44,2)
ppr_442 = z_extra(p_z44,3)
ppr_443 = z_extra(p_z44,6)
ppr_444 = z_extra(p_z44,7)
p_z46 = z_near(4, 0.1, -5.5)
ppr_461 = z_extra(p_z46,2)
ppr_462 = z_extra(p_z46,3)
ppr_463 = z_extra(p_z46,6)
ppr_464 = z_extra(p_z46,7)
p_z49 = z_near(12, 0.1, -5.5)
ppr_491 = z_extra(p_z49,2)
ppr_492 = z_extra(p_z49,3)
ppr_493 = z_extra(p_z49,6)
ppr_494 = z_extra(p_z49,7)
end
;监测 超孔压比、超孔压、有效应力
his id 639 ppr_431 ;(1 0.1 -5.5)
his id 640 ppr_432
his id 641 ppr_433
his id 642 ppr_434
his id 643 ppr_441 ;(2 0.1 -5.5)
his id 644 ppr_442
his id 645 ppr_443
his id 646 ppr_444
his id 651 ppr_461 ;(4 0.1 -5.5)
his id 652 ppr_462
his id 653 ppr_463
his id 654 ppr_464
his id 663 ppr_491 ;(12 0.1 -5.5)
his id 664 ppr_492
his id 665 ppr_493
his id 666 ppr_494
set hist_rep 300
set dyn dt 4e-5
def solve_ages
loop n (1,20)
time_n = n
save_file = 'Dyn-' + string(n) + 'x1s.sav'
command
sol age time_n
save save_file
endcommand
endloop
end
solve_ages |
|