Skip to contentUnited States Department of Transportation - Federal Highway AdministrationSearch FHWAFeedback

Geotechnical Engineering

 

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 =

 

<<Previous TOC Next>>
Updated: 04/07/2011

FHWA
United States Department of Transportation - Federal Highway Administration