Effects Of Geosynthetic Reinforcement Spacing On The Behavior Of Mechanically Stabilized Earth Walls
Appendix. Typical Input Data File
title
Effects of Geosynthetic Reinforcement Spacing on the Behavior of MSEW
;
define set_parameters
;
nstep='t5';
home_dir='c:\chris\';
fmov=string(home_dir)+'flac'+string(nstep)+'.mov'
command
movie file @fmov
end_command
;
; WALL DIMENSIONS AND GRID DENSITY PARAMETERS
hw=15
k1=0.6
lr=k1*hw
lra=10
lbc=5
nsl=2
;
; MATERIAL PROPERTY PARAMETERS
fri_g=45
ni_g=0.3
emod_gi=60000.0
smod_gi=emod_gi/(2.0*(1.0+ni_g))
bmod_gi=emod_gi/(3.0*(1.0-2.0*ni_g))
;
fri_s=45
ni_s=0.3
emod_si=60000.0
smod_si=emod_si/(2.0*(1.0+ni_s))
bmod_si=emod_si/(3.0*(1.0-2.0*ni_s))
;
fri_b=45
ni_b=0.3
emod_bi=60000.0
smod_bi=emod_bi/(2.0*(1.0+ni_b))
bmod_bi=emod_bi/(3.0*(1.0-2.0*ni_b))
;
emod_cb=600000.0
ni_cb=0.15
smod_cb=emod_cb/(2.0*(1.0+ni_cb))
bmod_cb=emod_cb/(3.0*(1.0-2.0*ni_cb))
;
mod_cab=1e6
stif_gr=smod_si
stre_gr=100
nstif_if=bmod_cb
sstif_if=smod_si
;
nl=hw/0.2
tnsl=nl*nsl
ah=0.2/nsl
nz=lr/ah
nz1=lra/ah
;
nzx=nz+nz1+43
nzy=(nsl+2)*nl+20
ngx = nz + 1
ngy = nzy + 1
;
c_x_1 = 5.2 + lr
c_x_2 = 10.2 + lr + lra
c_x_3 = 5.2 + lr + lra
c_y_1 = 5 + hv
;
z_i_1 = 22 +(lr+lra)/0.1
z_i_2 = z_i_1 + 1
z_i_3 = z_i_1 + 20
z_i_4 = nz + 22
z_i_5 = z_i_4 + 1
z_i_6 = nz + nz1 + 24
z_i_7 = nz + nz1 + 43
z_i_8 = nz + nz1 + 22
z_i_9 = nz + nz1 + 23
z_j_1 = 19+(nsl+2)*nl
z_j_2 = 21+tnsl
z_j_3 = 21+hw/0.1
;
g_i_1 = z_i_1+1
g_i_2 = z_i_3+1
g_i_3 = z_i_5
g_i_4 = z_i_6
g_i_5 = z_i_6 + 1
g_i_6 = z_i_9
g_j_1 = 20+(nsl+2)*nl
g_j_2 = 22+tnsl
g_j_3 = 22+hw/0.1
cx1=10.5 + lr + lra
end; "set_paremeters"
;
set_parameters
movie on
;
define e_modulus_ground
;
loop ig (1,z_i_3)
loop jg (1,20)
;
ag=abs(sxx(ig,jg))
agg=sqrt(ag)
ex_1(ig,jg)=20720000*agg
ex_2(ig,jg)=ex_1(ig,jg)/(2.0*(1.0+ni_g))
ex_3(ig,jg)=ex_1(ig,jg)/(3.0*(1.0-2.0*ni_g))
smod_g=ex_2(ig,jg)
bmod_g=ex_3(ig,jg)
;
command
prop bulk=bmod_g shear=smod_g i=ig j=jg
end_command
endloop
endloop
end; "e_modulus_ground"
;
define e_modulus_resoil
;
loop ir (23,z_i_8)
loop jr (22,wtzj)
;
ar=abs(sxx(ir,jr) )
arr=sqrt(ar)
ex_1(ir,jr)=2072*arr
ex_2(ir,jr)=ex_1(ir,jr)/(2.0*(1.0+ni_s))
ex_3(ir,jr)=ex_1(ir,jr)/(3.0*(1.0-2.0*ni_s))
smod_s=ex_3(ir,jr)
bmod_s=ex_2(ir,jr)
;
command
prop bulk=bmod_s shear=smod_s i=ir j=jr
end_command
end_loop
end_loop
end; "e_modulus_resoil"
;
define e_modulus_backfill
;
loop ib (z_i_6,z_i_7)
loop jb (22,ba_tzj)
;
ab=abs(min(sxx(ib,jb),szz(ib,jb)))
abb=sqrt(ab)
ex_1(ib,jb)=2072*abb
ex_2(ib,jb)=ex_1(ib,jb)/(2.0*(1.0+ni_b))
ex_3(ib,jb)=ex_1(ib,jb)/(3.0*(1.0-2.0*ni_b))
smod_b=ex_2(ib,jb)
bmod_b=ex_3(ib,jb)
;
command
prop bulk=bmod_b shear=smod_b i=ib j=jb
end_command
end_loop
end_loop
end; "e_modulus_backfill"
;
window -0.25 cx1 -0.25 cx1
;
; GROUND: GRID GENERATION & MODELLING
grid ngx, ngy
model mohr i=1,z_i_3 j=1,20
prop dens=2.201 bulk=bmod_gi shear=smod_gi fric=fri_g dil=15 coh=1000 &
i=1,z_i_3 j=1,20
gen 0,0 0,5 5,5 5,0 rat 0.91878121,0.91878121 i=1,21 j=1,21
gen same same 5.2,5 5.2,0 rat 1, 0.981878121 i=21,23 j=1,21
gen same same c_x_3,5 c_x_3,0 rat 1, 0.91878121 i=23,g_i_1 j=1,21
gen same same c_x_2,5 c_x_2,0 rat 1.0883984,0.91878121 i=g_i_1,g_i_2 j=1,21
;
; GROUND: BOUNDARY CONDITIONS
fix x i=1
fix y j=1
fix x i=g_i_2 j=1,21
;
; GROUND: INITIAL STRESSES
set gravity 9.81
ini sxx -33 var 0,33 j=1,20
ini syy -100 var 0,100 j=1,20
ini szz -33 var 0,33 j=1,20
e_modulus_ground
;
; HISTORIES
his 1 u
his 2 ydis i=38 j=20
his 3 ex_2 i=38 j=20
his 4 ex_3 i=38 j=20
;
define ground_equilibrium
;
grfname=string(home_dir)+'gr'+string(nstep)+'.sav'
prfname=string(home_dir)+'gr'+string(nstep)+'.pri'
;
command
prop tens 1e5 coh 1e5
solve srat 0.01
prop tens 0 coh 0
prop coh=1000 i=1,z_i_2 j=1,20
solve srat 0.01
end_command
;
e_modulus_ground
end; "ground_equilibrium"
;
ground_equilibrium
;
ini xdis 0 j=1,21
ini ydis 0 j=1,21
ini x 5.0 y 5.0 i=21 j=21
ini x 5.2 y 5.0 i=23 j=21
ini x c_x_3 y 5.0 i=g_i_1 j=21
;
save @grfname
;set log on
;set log @prfname
;print xdisp i=25
;set log off
;
define one_layer_parameters
;
; parameters for block generation
bbgj=22+(layer-1)*(nz+2)
btgj=bbgj+nsl
bbzj=bbgj
btzj=btgj-1
bfgi=17
bfzi=bfgi
;
y1=5+(layer-1)*0.2
y2=y1+0.2
;
; parameters for wall generation
wbzj=22+(layer-1)*nsl
wtzj=wbzj+nsl-1
wbgj=wbzj
wtgj=wtzj+1
;
; parameters for backfill generation
ba_bzj=22+(layer-1)*2
ba_tzj=ba_bzj+2
ba_bgj=ba_bzj
ba_tgj=ba_tzj+1
;
end; "one_layer_parameters"
;
define first_layer_generation
;
command
; concrete block generation
m e i=bfzi,21 j=bbzj,btzj
gen 5,y1 5,y2 5.2,y2 5.2,y1 i=bfgi,22 j=bbgj,btgj
prop d=2.199 b=bmod_cb s=smod_cb i=bfzi,21 j=bbzj,btzj
;
; a layer of wall generation
m m i=23,z_i_8 j=wbzj,wtzj; wall
gen 5.2,y1 5.2,y2 c_x_3,y2 c_x_3,y1 i=23,g_i_6 j=wbgj,wtgj
prop d=2.2 bulk=bmod_si shear=smod_si fric=fri_s dil=15 &
i=23,z_i_8 j=wbzj,wtzj
;
; a layer of backfill generation
m m i=z_i_6,z_i_7 j=ba_bzj,ba_tzj; backfill
gen c_x_3,y1 c_x_3,y2 c_x_2,y2 c_x_2,y1 rat 1.0883984,1 &
i=g_i_4,g_i_5 j=ba_bgj,ba_tgj
prop dens=2.2 bulk=bmod_bi shear=smod_bi fric=fri_b dil=15 &
i=z_i_6,z_i_7 j=ba_bzj,ba_tzj
;
att as fr g_i_6,wbgj to g_i_6,wtgj bs fr g_i_4,ba_bgj to g_i_4,ba_tgj
fix x i=g_i_5 j=ba_bgj,ba_tgj
;
end_command
;
end; "first_layer_generation"
;
define one_layer_generation
;
command
; concrete block generation
m e i=bfzi,21 j=bbzj,btzj
gen x1,y11 x6,y6 x7,y7 x2,y22 i=bfgi,22 j=bbgj,btgj
prop d=2.199 b=bmod_cb s=smod_cb i=bfzi,21 j=bbzj,btzj
;
; a layer of wall generation
m m i=23,z_i_8 j=wbzj,wtzj; wall
gen x3,y3 x8,y8 x9,y9 x4,y4 i=23,g_i_6 j=wbgj,wtgj
prop d=2.2 b=bmod_si s=smod_si fric=fri_s dil=15 i=23,z_i_8 j=wbzj,wtzj
;
; a layer of backfill generation
m m i=z_i_6,z_i_7 j=ba_bzj,ba_tzj
gen x5,y5 x10,y10 c_x_2,y2 c_x_2,y1 rat 1.0883984,1 &
i=g_i_4,g_i_5 j=ba_bgj,ba_tgj
prop d=2.2 b=bmod_bi s=smod_bi fric=fri_b dil=15 &
i=z_i_6,z_i_7 j=ba_bzj,ba_tzj
att as fr g_i_6,wbgj to g_i_6,wtgj bs fr g_i_4,ba_bgj to g_i_4,ba_tgj
fix x i=g_i_5 j=ba_bgj,ba_tgj
end_command
;
end; "one_layer_generation"
;
define top_coord_upgrade
;
x1=x(bfgi,btgj)
y11=y(bfgi,btgj)
x2=x(22,btgj)
y22=y(22,btgj)
x3=x(23,wtgj)
y3=y(23,wtgj)
x4=x(g_i_6,wtgj)
y4=y(g_i_6,wtgj)
x5=x(g_i_4,ba_tgj)
y5=y(g_i_4,ba_tgj)
;
end; "top_coord_upgrade"
;
define define_top_coord
;
x6=x1-y22+y11
y6=y11+x2-x1
x7=x2-y22+y11
y7=y22+x2-x1+1
x8=x7
y8=y7
x9=x4
y9=y7
x10=x4
y10=y7
;
end; "define_top_coord"
;
; BEAM & CABLE PROPERTIES
struct prop=2 area=0.002 peri=2 e=mod_cab yield=200 d=2.1
struct prop=2 kbond=stif_gr sbond=stre_gr sfriction=35
struct prop=3 area=0.002 e=mod_cab i=6.67e-10 d=2.1
;
define sos
;
fname=string(home_dir)+'l'+string(nstep)+'_'+string(layer)+'.sav'
fnameend=string(home_dir)+'end'+string(nstep)+'.sav'
prfname=string(home_dir)+'l'+string(nstep)+'_'+string(layer)+'.pri'
;
command
solve srat 0.01
end_command
;
command
save @fname
;set log on
;set log @prfname
;print xdisp i=25
;set log off
end_command
;
e_modulus_ground
e_modulus_resoil
e_modulus_backfill
;
end; "sos" = solve & save
;
define wall_construction
;
h_lbc=lbc*0.2
nl1_f=(hw-0.2)/h_lbc
nl1=int(nl1_f)
h_top=hv-0.2-nl1*h_lbc
if h_top=0 then
tnc=nl1
else
tnc=nl1+1
endif
nrl=1
;
layer=1
;
one_layer_parameters
first_layer_generation
;
command
att as fr 23,21 to g_i_1,21 bs fr 23,22 to g_i_6,22
att as fr g_i_1,21 to g_i_2,21 bs fr g_i_6,22 to g_i_5,22
int 1 as fr 19,21 to 24,21 bs long fr bfgi,btgj to 22,btgj
int 1 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int 2 as fr 22,bbgj to 22,btgj bs fr 23,wbgj to 23,wtgj
int 2 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
sos
top_coord_upgrade
b1e_no=6
;
command
struct beam begin x1,y11 end x2,y2 seg=5 prop=3
struct node range 1 6 pin
struct cable begin no b1e_no end c_x_1,y2 seg=nz prop=2
int 3 as fr no b1e_no to no 1 bs fr 22,btgj to bfgi,btgj
int 3 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
tot_nu_int=3
;
loop layer (2,nl)
;
top=1+nrl*lbc
bottom=top-lbc+1
;
if top=bottom then
;
one_layer_parameters
define_top_coord
one_layer_generation
;
bbf_no=(nz+6)*(layer-1)+1
bbb_no=bbf_no+5
btf_no=(nz+6)*(layer-2)+1
btb_no=btf_no+5
b_gj=bbgj-2
b_gjb=bbgj-nsl-2
f_1=(layer-2)*5+4
f_2=f_1+1
f_3=f_1+2
f_4=f_1+3
f_5=f_1+4
tot_nu_int=f_5
;
command
int f_1 as fr no bbf_no to no bbb_no bs fr bfgi,bbgj to 22,bbgj
int f_1 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_2 as fr 22,bbgj to 22,btgj bs fr 23,wbgj to 23,wtgj
int f_2 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_3 as fr bfgi,bbgj to 22,bbgj bs fr bfgi,b_gj to 22,b_gj
int f_3 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_4 as fr 23,wbgj to 23,wtgj bs fr 22,b_gjb to 22,b_gj
int f_4 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
sos
top_coord_upgrade
;
command
stru beam begin x6,y6 end x7,y7 seg 5 prop 3
stru cable begin no btb_no end c_x_1,y9 seg nz prop 2
int f_5 as fr no btb_no to no btf_no bs fr 22,btgj to bfgi,btgj
int f_5 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
nrl=nrl+1
;
else
;
if layer=bottom then
;
one_layer_parameters
define_top_coord
one_layer_generation
;
bbf_no=(nz+6)*(nrl-2)+1
bbb_no=bbf_no+5
b_gj=bbgj-2
b_gjb=bbgj-nsl-2
f_6=tot_nu_int+1
f_7=f_6+1
f_8=f_6+2
f_9=f_6+3
tot_nu_int=f_9
;
command
int f_6 as fr no bbf_no to no bbb_no bs fr bfgi,bbgj to 22,bbgj
int f_6 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_7 as fr 22,bbgj to 22,btgj bs fr 23,wbgj to 23,wtgj
int f_7 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_8 as fr bfgi,bbgj to 22,bbgj bs fr bfgi,b_gj to 22,b_gj
int f_8 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_9 as fr 23,wbgj to 23,wtgj bs fr 22,b_gjb to 22,b_gj
int f_9 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
sos
top_coord_upgrade
;
else
;
if layer<top then
;
one_layer_parameters
define_top_coord
one_layer_generation
;
b_gj=bbgj-2
b_gjb=bbgj-nsl-2
f_10=tot_nu_int+1
f_11=f_10+1
f_12=f_10+2
tot_nu_int=f_12
;
command
int f_10 as fr bfgi,b_gj to 22,b_gj bs fr bfgi,bbgj to 22,bbgj
int f_10 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_11 as fr 22,bbgj to 22,btgj bs fr 23,wbgj to 23,wtgj
int f_11 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_12 as fr 23,wbgj to 23,wtgj bs fr 22,b_gjb to 22,b_gj
int f_12 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
sos
top_coord_upgrade
;
else
;
one_layer_parameters
define_top_coord
one_layer_generation
b_gj=bbgj-2
b_gjb=bbgj-nsl-2
tbf_no=(nz+6)*nrl+1
tbb_no=tbf_no+5
f_13=tot_nu_int+1
f_14=f_13+1
f_15=f_13+2
f_16=f_13+3
tot_nu_int=f_16
;
command
int f_13 as fr bfgi,b_gj to 22,b_gj bs fr bfgi,bbgj to 22,bbgj
int f_13 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_14 as fr 22,bbgj to 22,btgj bs fr 23,wbgj to 23,wtgj
int f_14 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
int f_15 as fr 23,wbgj to 23,wtgj bs fr 22,b_gjb to 22,b_gj
int f_15 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
sos
top_coord_upgrade
;
command
stru beam begin x6,y6 end x7,y7 seg 5 prop 3
stru cable begin no tbb_no end c_x_1,y9 seg nz prop 2
int f_16 as fr no tbb_no to no tbf_no bs fr 22,btgj to bfgi,btgj
int f_16 f 30 kn nstif_if ks sstif_if coh 0; tbond 30
end_command
;
nrl=nrl+1
;
endif
endif
endif
endloop
end; "wall_construction"
;
wall_construction
solve srat 0.01
save @fnameend
;
; = END OF FILE =
|