Hi @everyone, @itascan
Can you please have a look at the program below. Do you see anything missing with the geogrids or anything else that I am not getting change in the output with and without addition of geogrid layers?
model restore ‘initial’
zone gridpoint ini disp 0 0 0
zone gridpoint initialize fluid-modulus 2e6
zone gridpoint ini fluid-tension -1e20
zone history name ‘S_pg_1=0’ displacement-z position (0.5,0.5,-0.1)
zone history name ‘S_Pg_2=-2’ displacement-z position (0.5,0.5,-2)
zone history name ‘S_Pg_3=-8’ displacement-z position (0.5,0.5,-8)
; mechanical properties for fill
zone cmodel assign mohr-coulomb range group ‘L1’
zone property density 2100 young 70e6 poiss 0.3 …
friction 36 cohesion 61000 dilation 3 range group ‘L1’
; set saturation to 0 for nodes in the fill
zone gridpoint ini saturation 0 range group ‘L1’ …
pos-z 0.1 0.51
; create piles - put top node in fill
structure pile create by-line (0.0975, 0.0975, 0.05) (0.0975, 0.0975, -8.4) segments=36
structure pile property young=1.8e10 poisson=0.2 cross-sectional-area=0.038025 …
torsion-constant=0.0 moi-y=0 moi-z=0 perimeter=0.78 …
coupling-stiffness-shear=6.03e9 …
coupling-cohesion-shear=35000 …
coupling-friction-shear=0 …
coupling-stiffness-normal=0 …
coupling-cohesion-normal=0 …
coupling-friction-normal=0 coupling-gap-normal=off
; =======================================================
;; Assigning interface
structure pile property coupling-stiffness-shear=6.03e9 …
coupling-cohesion-shear=45000 coupling-friction-shear 0 range position-z -4 0.05;
; Include end-bearing effect
structure link delete range position-z -8.4
structure link create target zone group ‘End’ range position-z -8.4
structure link attach x=normal-yield y=free z=free …
range group ‘End’
structure link attach rotation-x=free rotation-y=free rotation-z=free …
range group ‘End’
structure link property x area=1 stiffness=1.8e11 yield-compression=44691 …
range group ‘End’
;
; Create two geogrid layers, and assign properties
structure geogrid create by-quadrilateral (0,0,0.2) (1,0,0.2) …
(1,1,0.2) (0,1,0.2) size (1,1)
structure geogrid create by-quadrilateral (0,0,0.4) (1,0,0.4) …
(1,1,0.4) (0,1,0.4) size (1,1)
structure geogrid cmodel assign elastic
structure geogrid property young 1.04e8 poisson 0.3 thickness=5e-3 …
coupling-stiffness=2.3e6 coupling-cohesion=0.0 …
coupling-friction=25 ; 25 coupling friction!
model fluid flow off;
model solve ratio 1e-3