Microsoft Word - ETASR_V13_N1_pp9843-9848


Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9843 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

The Effect of Rigid Inclusions on the Dynamic 
Response of Highway Embankment 

 

Faris Bouabdallah 

Department of Civil Engineering, Laboratory LMGHU, 20 August 1955 University of Skikda, Algeria 
f.bouabdallah@univ-skikda.dz 
(corresponding author) 
 
Kamel Goudjil 

Department of Civil Engineering, University of Mohamed-Cherif Messaadia, Algeria 
k_goudjil@yahoo.com 
 
Messast Salah 

Department of Civil Engineering, Laboratory LMGHU, 20 August 1955 University of Skikda, Algeria 
msalah2007msalah@gmail.com 
 

Received: 8 October 2022 | Revised: 28 October 2022 | Accepted: 6 November 

 

ABSTRACT 

Soft soils are widespread in many areas of the East-West highway of Algeria, covering large areas. 

Construction projects such as highway embankments, airfields, roads, and railways in such areas 

experience several problems due to their low strength and permeability. This study used a fully dynamic 

analysis to evaluate the dynamic response of a highway embankment. The numerical analysis was carried 

out using the FLAC 2D version 8.10 software. The main objective was to evaluate the effect of rigid 

inclusions on the dynamic response of highway embankments by monitoring settlement, shear-strain 

curves, and maximum soil accelerations and displacements at different depth points. The obtained results 

showed that the configuration of the rigid inclusions significantly influences the dynamic response of the 

model. For instance, the settlement at 0.9m is 16% higher than at 0.2m from rigid inclusion. 

Keywords-inclusions; dynamic response; flac; settlement; shear strain curve  

I. INTRODUCTION  

Several studies have been carried out to better understand 
the behavior of soft soils, due to the spread of soft soils, which 
cover several areas in Algeria, and develop disaster prevention 
measures, estimate their parameters, and study the effect of 
some materials such as rubber concrete and lime on the 
behavior of composite clayey soils [1-5]. Structures such as 
highway or railway embankments, tanks, walls, or slopes face 
several problems with such soft foundations. These problems 
relate to bearing capacity failures, intolerable total and 
differential settlements, large lateral pressures, and slope 
instabilities. Several studies showed that the principal cause of 
several landslides, which affected the East-West highway in 
Algeria, is the presence of clayey formations [6-7]. A variety of 
ground improvement techniques have been used to address the 
concerns on soft soil foundations, such as geosynthetic [8-11], 
geogrid [12-14], pile [15-16], PVD [17], soil nail [18-19], 
chemical injection and grouting [20], stone column [21-24], 
and rigid inclusions [25-27]. However, the suitable solution 
depends on several parameters, such as the thickness of the soft 

soil, the nature of the project to build, the requirements of the 
client, etc. 

Rigid inclusions are a ground improvement technique that 
has a significant advantage over deep foundations, as they can 
provide the needed settlement control and retain the shallow 
foundation support of the structure. Rigid inclusions, known by 
different names such as Controlled Modulus Columns (CMC), 
pile-supported earth platforms, jet grouting columns, or soil 
column reinforcement, are vertical elements across soil layers 
with low bearing capacity and/or high compressibility, 
extending down to a more resistant layer. Due to their higher 
stiffness compared to surrounding soil, they support a portion 
of the loads applied at the ground's surface. Therefore, the 
loads taken by the soft soil can be reduced to an acceptable 
level of soil-bearing capacity and settlement. Several studies 
investigated the effect of rigid inclusions to improve soft soils. 
Some studies used the Stress Reduction Ratio (SRR) to find the 
optimal design [28-29]. However, they did not consider several 
factors that can potentially influence the SRR, such as the 
characteristics of the subgrade soils and the influence of 
dynamic loadings. Other approaches adopted rigid inclusions 
techniques to improve soft soil intended to support a railroad 



Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9844 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

embankment. The rigid inclusions reduced the ground 
settlements by 15-20%, depending on the subgrade plasticity 
index [30]. Several numerical studies investigated the seismic 
behavior of highway embankments in improved soft soils. 
These studies used the linear elastic-perfectly plastic 
constitutive model with a Mohr-Coulomb failure criterion to 
represent the behavior of the soil and LTP [31-36]. Some 
studies used the constitutive models UBCYST and PM4silt to 
take into account the nonlinear behavior of soils during seismic 
motions [37-40]. 

This paper presents a fully dynamic analysis and evaluation 
of the seismic response of highway embankments over rigid 
inclusions to improve soft soil. The UBCHYST model [41] 
was adopted to consider the dynamic characteristics of the soil 
during seismic movements. The effect of ground motions of the 
CHI-CHI earthquake and two types of rigid inclusion 
configurations, floating and placed in hard clay, were studied. 

II. NUMERICAL MODELING 

A. The Geometry of the Numerical Model 

A 3-layer soil model was considered, where 10m thick soft 
soil was built over a 5m thick hard soil layer. The highway was 
constructed on an embankment with an 11.75m width and a 
slope angle of approximately 18°. The toe of the embankment 
was at an elevation of 15m. The embankment includes a load 
transfer platform 0.6m thick as part of the embankment 
materials. To simplify the problem, the soil layers were 
supposed to be horizontal, as shown in Figure 1. The water 
table was not set up in the model and drained conditions were 
considered. A finite difference mesh with a size of 0.25cm for 
each element was chosen to avoid wave distortion during 
seismic wave transmission through the soil [42]. The 
constitutive model UBCHYST can be used with low and high-
permeability soils, such as clayey, silty, and granular. In this 
hysteretic model, the shear modulus is a function of the stress 
ratio and varies throughout the loading cycle to give hysteretic 
stress-strain loops of varying amplitude and area damping 
throughout the earthquake excitation. A detailed description of 
all UBCHYST model parameters was given in [41]. 

 

 
Fig. 1.  Geometry of the numerical model 

TABLE I.  RIGID INCLUSION PARAMETERS  

Density 

(kg/m³)  

Shear modulus G 

(MPa) 

Young modulus E 

(MPa) 

25000 12.5 30 

Tables I and II present the parameters of the soil and 
structural elements. Soil characteristics were taken from [31, 
39, 43]. The calibration parameters of the constitutive 
UBCHYST model were determined based on laboratory 
experiments [39, 43]. Since the soft soil used in [31] was 
modeled as the Mhor-Coulomb constitutive model, a cyclic 
shear test was conducted to determine the calibration 
parameters of the UBCHYST model using the modulus 
reduction curves from [44]. 

TABLE II.  SOIL PARAMETERS OF THE UBCHYST MODEL 

Soil parameters Stiff clay Soft clay LTP Earth fill 

Bulk modulus K (MPa) 76 8.33 41.66 76 
Shear modulus G (MPa) 38 3.8 19 38 

Density (kg/m³),  2000 1600 2000 2000 
Friction angle Ф (deg) 30 25 25 35 

Cohesion C (kPa) 15.2 5 50 0 
Tensile T 0 0 0 0 

Dilation angle dil 0 0 0 0 
Stress rate factor Rf 1 0.8 0.98 1 

Stress rate exponent n 2 2 2.5 2 
First cycle factor Mod1 0.75 0.75 1 0.75 

Large strains exponent Rm 0.5 0.5 0.5 0.5 
Large strains factor dfac 0 0 0 0 

 

B. Rigid Vertical Elements  

Numerical analysis was conducted using 12 rigid inclusions 
per row in the soft soil, with 10m length and 0.6m diameter. 
The separation between elements in both directions was equal 
to 1.8m. The rigid inclusion characteristics were taken from 
[31].  

C. Boundary Conditions and Interfaces 

Artificial boundaries were used in the static analysis to 
depict the semi-infinite nature of the soil. The side boundaries 
were fixed in the horizontal direction while the bottom part was 
fixed in all directions. However, these boundary conditions 
were replaced in the dynamic analysis by quiet boundaries 
along the bottom of the model to minimize the effects of 
reflected waves at the bottom and free-field boundaries on the 
sides to avoid wave reflections[45]. To apply quiet boundary 
conditions along the same boundary as the dynamic input, the 
dynamic input must be applied as a shear stress boundary, 
because the effect of the quiet boundary will be nullified if the 
input is applied as an acceleration or velocity wave. To do this, 
the velocity record was converted into a shear stress boundary 
condition using 

�� � �����	 
 �� 
 
� � 
 ��    (1) 
where σs is the applied shear stress, ρ is the mass density of the 
material at the boundary, Cs is the speed of s-wave propagation 
through the medium at the boundary, vs is the input shear 
particle velocity, and the factor is generally equal to 2.  

Figures 2 and 3 show the time history acceleration and 
velocity of the CHI-CHI earthquake. Regarding soil-rigid 
inclusion interfaces, since the soil foundations of this study 
were clayey soils, the interfaces were assumed to have a zero-
friction angle and the same cohesive strength as the 
surrounding soil [46]. The shear and normal stiffness values 
were determined using (2), based on the FLAC manual [45]. 



Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9845 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

�� � �� � 10 �
���.� �
∆��� 

!   (2) 

where K and G are the bulk and shear modulus of the adjacent 
soil, respectively, and ∆Zmin is the smallest edge of the adjacent 
soil element. 

 

 
Fig. 2.  Input acceleration time history of the CHI-CHI earthquake in 1999. 

 
Fig. 3.  Input velocity time history of the CHI-CHI earthquake in 1999. 

III. RESULTS AND DISCUSSION 

The soil response in terms of settlements at the bottom of 
the embankment and along the depth, and shear stresses and 
strains were evaluated for both cases of Rigid Inclusions (RI), 
placed on hard clay and floating. All these results took into 
account only the influence of the dynamic loading on the soft 
soils improved with rigid inclusions. In the first step, an initial 
stress state was generated. Then, the vertical reinforcements 
were installed, and the last static calculation step considered the 
activation of the earth platform and the highway embankment. 
The required seismic boundary conditions, free field 
boundaries, and quiet boundaries were added for the dynamic 
analysis, and the dynamic analysis was executed by applying 
the corresponding horizontal wave using shear stress (2) at the 
base of the models.  

A.  Settlements at the Bottom of the Embankment 

Figure 4 shows the histories of soil settlements at the base 
of LTP for different points relative to the RI. It is noticeable 
that the settlement at 0.9m from RI is greater than at 0.2m. For 
example, the settlement of the point at 0.9m is 16% higher than 
at 0.2m from RI. The settlement decreases the closer the point 
is to the rigid inclusion. This can be explained by the effect of 
shear stiffness provided via the interaction of the rigid elements 
with the surrounding soil, which leads to reduce the settlement 
of the surrounding soil. 

B. Settlements along the Depth 

Figure 5 shows the histories of soil settlements at different 
depths along the centreline axis. It is clear that the maximum 

settlements (6.4cm) occur at the top of the LTP layer and are 
reduced with depth. The maximum settlement in the middle of 
soft clay (3.3cm) is 23.6% greater compared to those in its base 
(0.78cm). 

 

 
Fig. 4.  Settlement at the bottom of the embankment. 

 
Fig. 5.  Settlements along the depth. 

C. Shear Strains and stress curves 

Figures 6-8 show the hysteresis curve that relates shear 
stresses and shear strains obtained after dynamic loading in the 
soil.  

 

 
Fig. 6.  Shear stress and strain curve at the middle of hard clay. 

 
Fig. 7.  Shear stress and strain curve at the middle of soft clay. 

0 10 20 30 40 50 60 70 80 90
-3

-2

-1

0

1

2

3

Dynamic time (second)

In
p

ut
 a

cc
el

er
at

io
n

 (
m

/s
2
)

0 20 40 60 80 100
-0.4

-0.2

0

0.2

0.4

0.6

Dynamic time (second)

In
pu

t 
ve

lo
ci

ty
 m

/s
)

0 20 40 60 80 100
-0.06

-0.05

-0.04

-0.03

-0.02

-0.01

0

0.01

Dytime (second)

S
et

tl
em

en
t 

(m
)

 

 

0.9 m from R.I

0.4 m from R.I

0.2 m from R.I

0 10 20 30 40 50 60
-0.07

-0.06

-0.05

-0.04

-0.03

-0.02

-0.01

0

0.01

D y n amic t i me ( s econ d )

S
e
tt

le
m

e
n

t 
 (

m
)

 

 

at the base of model
at the top of hard clay
at the middle of soft clay
at 2m from the top of soft clay
at the top of LTP layer

-5 0 5 10

x 10
-4

-3

-2

-1

0

1

2

3

4
x 10

4

 

 

R.I Placed on hard clay

Max shear st ress = 30.2 kpa
Max shaer st rain= 9.75*10-4

-5 0 5 10 15

x 10
-3

-2

-1

0

1

2

3
x 10

4

Shear Strain

S
h
ea

r 
S

tr
es

s 
(k

p
a)

 

 

R.I placed on hard clay



Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9846 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

 

Fig. 8.  Shear stress and strain curve at the middle of LTP. 

 
Fig. 9.  Floating rigid inclusions. 

The values obtained were recorded in the middle of the soft 
clay and the LTP layer. Figures 3-5 show maximum shear 
strain at 9.75×10-4, 13×10-3, and 8.55×10-4, respectively. The 
distribution of the maximum shear strain is higher in the middle 
of the soft clay. At the LTP layer and the middle hard clay, the 
maximum shear strains are very close, but the shear stresses 
differ. These results are due to the higher value of shear 
stiffness of the hard clay and the LTP layer compared to the 
soft clay. Figure 9 shows 8m long floating rigid inclusions 
embedded in the soft soil. Shear stresses and strains at the same 
points were estimated to show the efficiency of rigid 
inclusions. Figures 10 and 11 show the hysteresis curve that 
relates shear stresses and strains obtained after dynamic loading 
in the soil for floating rigid inclusions. The obtained values 
were recorded at the same points for rigid inclusions placed on 
hard clay. The maximum shear strain value in the middle of 
soft clay was 15.2×10-3, having an increase of 2.2×10-3. Also, 
the maximum value of shear strain in the LTP layer reached 
9.63×10-4, having an increase of 1.08×10-4. 

 

 
Fig. 10.  Shear stress and strain curve at the middle of soft clay. 

 
Fig. 11.  Shear stress and strain curve at the middle of LTP. 

IV. CONCLUSION 

A fully dynamic numerical modeling analysis was 
conducted to evaluate the seismic response of a highway 
embankment constructed over reinforced soft soil. Two 
configuration types were considered to study the effect of rigid 
inclusions, placed on hard soil and floating. The numerical 
modeling results of all configurations were compared in terms 
of settlement, shear stress, and shear strain. Settlements at the 
base of the embankment were evaluated in three different 
positions relative to the position of the rigid inclusion. The 
settlement was important at 0.9m from the rigid inclusions 
while the minimum settlement was recorded at 0.2m. This can 
be explained by the shear stiffness provided via the interaction 
of the rigid elements with the surrounding soil, which reduces 
the settlement of the surrounding soil.The maximum value of 
shear strains was recorded in the middle of the soft clay 
compared to the LTP layer and the middle of the hard clay for 
all cases. The cohesion of the LTP layer and the hard clay 
increases the strength of the soil and reduces shear strain. In 
floating rigid inclusions, the maximum shear strain in the 
middle of soft clay was 15.2×10-3, which is an increase of 
2.2×10-3. The maximum value of shear strain in the LTP layer 
reached 9.63×10-4, which is an increase of 1.08×10-4. These 
results are due to the shear force that can develop along the 
rigid inclusion/soil interface, which is a function of the 
cohesive strength of the interface. Several studies investigated 
the effect of various parameters such as rigid inclusion 
diameter, spacing, seismic force intensity, etc. However, few 
studies investigated the influence of rigid inclusions using a 
fully dynamic analysis. This paper presented a dynamic study 
using the UBCHYST model to highlight the nonlinear dynamic 
behavior of soft soils. 

REFERENCES 
[1] S. M. Aissa Mamoune, "Contribution à la mesure, prévision et 

modélisation du comportement des sols expansifs," MSc Thesis, 
University of Tlemcen, Algeria, 2002. 

[2] M. Afes, "Contribution à la détermination des paramètres de gonflement 
des sols, et étude de l’argile de Mila (Algérie) traitée à la chaux," Ph.D. 
dissertation, Lyon, INSA, 1996. 

[3] D. Athmania, A. Benaissa, A. Hammadi, and M. Bouassida, "Clay and 
Marl Formation Susceptibility in Mila Province, Algeria," Geotechnical 
and Geological Engineering, vol. 28, no. 6, pp. 805–813, Nov. 2010, 
https://doi.org/10.1007/s10706-010-9341-5. 

[4] F. Klouche, E. Kara Mostefa, A. Sekkel, M. Maliki, and M. 
Bouguenaya, "Effect of waste concrete on the compressibility of clay 
from the Bouhanak region (Algeria)," Euro-Mediterranean Journal for 

-6 -4 -2 0 2 4 6 8 10

x 10
-4

-1.5

-1

-0.5

0

0.5

1

1.5
x 10

4

Shear Strain

S
he

ar
 S

tr
es

s 
(K

p
a)

 

 

R.I placed on hard clay

-5 0 5 10 15 20

x 10
-3

-2

-1

0

1

2

3
x 10

4

Shear St rain

S
h

ea
r 

S
tr

es
s 

(p
a)

 

 

Case of Float ing
Rigid Inclusions

-1 -0.5 0 0.5 1

x 10
-3

-3

-2

-1

0

1

2

3
x 10

4

Shear strain

S
he

ar
 s

tr
es

s 
(k

pa
)

 

 

Case of Floating
Rigid inclusion



Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9847 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

Environmental Integration, vol. 7, no. 2, pp. 191–200, Jun. 2022, 
https://doi.org/10.1007/s41207-022-00308-4. 

[5] H. Trouzine, M. Bekhiti, and A. Asroun, "Effects of scrap tyre rubber 
fibre on swelling behaviour of two clayey soils in Algeria," 
Geosynthetics International, vol. 19, no. 2, pp. 124–132, Apr. 2012, 
https://doi.org/10.1680/gein.2012.19.2.124. 

[6] B. E. Kechebour and A. Talah, "Durability of New Roads: Case Study of 
the Algerian East West Highway," International Journal of Structural 
and Civil Engineering Research, vol. 7, no. 1, pp. 87–91, Feb. 2018, 
https://doi.org/10.18178/ijscer.7.1.87-91. 

[7] A. Benaissa and M. A. Bellouche, "Propriétés géotechniques de 
quelques formations géologiques propices aux glissements de terrains 
dans l’agglomération de Constantine (Algérie)," Bulletin of Engineering 
Geology and the Environment, vol. 57, no. 3, pp. 301–310, Mar. 1999, 
https://doi.org/10.1007/s100640050049. 

[8] J. Han and M. A. Gabr, "Numerical Analysis of Geosynthetic-
Reinforced and Pile-Supported Earth Platforms over Soft Soil," Journal 
of Geotechnical and Geoenvironmental Engineering, vol. 128, no. 1, pp. 
44–53, Jan. 2002, https://doi.org/10.1061/(ASCE)1090-0241(2002)128: 
1(44). 

[9] P. Ariyarathne and D. S. Liyanapathirana, "Review of existing design 
methods for geosynthetic-reinforced pile-supported embankments," Soils 
and Foundations, vol. 55, no. 1, pp. 17–34, Feb. 2015, https://doi.org/ 
10.1016/j.sandf.2014.12.002. 

[10] T. Thyagaraj, Ed., "Ground Improvement Techniques and 
Geosynthetics: IGC 2016 Volume 2," in IGC 2016 Vol. 2, Chennai, 
India, 2019, vol. 14, https://doi.org/10.1007/978-981-13-0559-7. 

[11] M. Touahmia, "Performance of Geosynthetic-Reinforced Soils Under 
Static and Cyclic Loading," Engineering, Technology & Applied Science 
Research, vol. 7, no. 2, pp. 1523–1527, Apr. 2017, https://doi.org/ 
10.48084/etasr.1035. 

[12] M. S. S. Almeida, M. Ehrlich, A. P. Spotti, and M. E. S. Marques, 
"Embankment supported on piles with biaxial geogrids," Proceedings of 
the Institution of Civil Engineers - Geotechnical Engineering, vol. 160, 
no. 4, pp. 185–192, Oct. 2007, https://doi.org/10.1680/geng.2007.160. 
4.185. 

[13] H. L. Liu, C. W. W. Ng, and K. Fei, "Performance of a Geogrid-
Reinforced and Pile-Supported Highway Embankment over Soft Clay: 
Case Study," Journal of Geotechnical and Geoenvironmental 
Engineering, vol. 133, no. 12, pp. 1483–1493, Dec. 2007, https://doi.org/ 
10.1061/(ASCE)1090-0241(2007)133:12(1483). 

[14] B. Mazouz, T. Mansouri, M. Baazouzi, and K. Abbeche, "Assessing the 
Effect of Underground Void on Strip Footing Sitting on a Reinforced 
Sand Slope with Numerical Modeling," Engineering, Technology & 
Applied Science Research, vol. 12, no. 4, pp. 9005–9011, Aug. 2022, 
https://doi.org/10.48084/etasr.5131. 

[15] R. P. Chen, Z. Z. Xu, Y. M. Chen, D. S. Ling, and B. Zhu, "Field Tests 
on Pile-Supported Embankments over Soft Ground," Journal of 
Geotechnical and Geoenvironmental Engineering, vol. 136, no. 6, pp. 
777–785, Jun. 2010, https://doi.org/10.1061/(ASCE)GT.1943-5606. 
0000295. 

[16] C. Yun-min, C. Wei-ping, and C. Ren-peng, "An experimental 
investigation of soil arching within basal reinforced and unreinforced 
piled embankments," Geotextiles and Geomembranes, vol. 26, no. 2, pp. 
164–174, Apr. 2008, https://doi.org/10.1016/j.geotexmem.2007.05.004. 

[17] X. Jiang, Y. Jiang, C. Y. Wu, W. Q. Wang, J. Y. Geng, and Y. J. Qiu, 
"Numerical analysis for widening embankments over soft soils treated 
by PVD and DJM columns," International Journal of Pavement 
Engineering, vol. 21, no. 3, pp. 267–279, Feb. 2020, https://doi.org/ 
10.1080/10298436.2018.1461869. 

[18] D. A. Mangnejo, S. J. Oad, S. A. Kalhoro, S. Ahmed, F. H. Laghari, and 
Z. A. Siyal, "Numerical Analysis of Soil Slope Stabilization by Soil 
Nailing Technique," Engineering, Technology & Applied Science 
Research, vol. 9, no. 4, pp. 4469–4473, Aug. 2019, https://doi.org/ 
10.48084/etasr.2859. 

[19] W. R. Azzam and A. Basha, "Utilization of soil nailing technique to 
increase shear strength of cohesive soil and reduce settlement," Journal 
of Rock Mechanics and Geotechnical Engineering, vol. 9, no. 6, pp. 
1104–1111, Dec. 2017, https://doi.org/10.1016/j.jrmge.2017.05.009. 

[20] Sina Kazemain, "Review of soft soils stabilization by grouting and 
injection methods with different chemical binders," Scientific Research 
and Essays, vol. 7, no. 24, Jun. 2012, https://doi.org/10.5897/ 
SRE11.1186. 

[21] Z. Guetif, M. Bouassida, and J. M. Debats, "Improved soft clay 
characteristics due to stone column installation," Computers and 
Geotechnics, vol. 34, no. 2, pp. 104–111, Mar. 2007, https://doi.org/ 
10.1016/j.compgeo.2006.09.008. 

[22] J. Castro, M. Karstunen, and N. Sivasithamparam, "Influence of stone 
column installation on settlement reduction," Computers and 
Geotechnics, vol. 59, pp. 87–97, Jun. 2014, https://doi.org/10.1016/ 
j.compgeo.2014.03.003. 

[23] S. W. Abusharar and J. Han, "Two-dimensional deep-seated slope 
stability analysis of embankments over stone column-improved soft 
clay," Engineering Geology, vol. 120, no. 1, pp. 103–110, Jun. 2011, 
https://doi.org/10.1016/j.enggeo.2011.04.002. 

[24] S. Saxena and L. B. Roy, "Suitability Analysis of Stone Column 
Materials with PLAXIS," Engineering, Technology & Applied Science 
Research, vol. 12, no. 2, pp. 8421–8425, Apr. 2022, https://doi.org/ 
10.48084/etasr.4761. 

[25] J. A. B. Garcia, J. F. Rodríguez Rebolledo, D. V. dos Santos 
Mützenberg, B. Caicedo, and G. de Farias Neves Gitirana, 
"Experimental Investigation of a Load-Transfer Material for 
Foundations Reinforced by Rigid Inclusions," Journal of Geotechnical 
and Geoenvironmental Engineering, vol. 147, no. 10, Oct. 2021, Art. no. 
04021110, https://doi.org/10.1061/(ASCE)GT.1943-5606.0002649. 

[26] L. Briançon, D. Dias, and C. Simon, "Monitoring and numerical 
investigation of a rigid inclusions–reinforced industrial building," 
Canadian Geotechnical Journal, Mar. 2015, https://doi.org/10.1139/cgj-
2014-0262. 

[27] V. D. Tran, J. J. Richard, and T. Hoang, "Soft Soil Improvement Using 
Rigid Inclusions: Toward an Application for Transport Infrastructure 
Construction in Vietnam," in New Prospects in Geotechnical 
Engineering Aspects of Civil Infrastructures, HangZhou, China, Jul. 
2018, pp. 89–99, https://doi.org/10.1007/978-3-319-95771-5_8. 

[28] D. Mazzei, K. Kniss, F. Elsaid, and Y. Zhang, "Rigid Inclusions Ground 
Improvement for a New Energy Facility: Design, Construction, and Full-
Scale Embankment Load Testing and Results," pp. 101–114, Mar. 2019, 
https://doi.org/10.1061/9780784482117.009. 

[29] P. Burtin and J. Racinais, "Embankment on Soft Soil Reinforced by 
CMC Semi-Rigid Inclusions for the High-speed Railway SEA," 
Procedia Engineering, vol. 143, pp. 355–362, Jan. 2016, https://doi.org/ 
10.1016/j.proeng.2016.06.045. 

[30] D. Wang, M. Sánchez, and J.-L. Briaud, "Numerical study on the effect 
of rigid inclusions on existing railroads," International Journal for 
Numerical and Analytical Methods in Geomechanics, vol. 43, no. 18, pp. 
2772–2796, 2019, https://doi.org/10.1002/nag.3001. 

[31] G. A. L. Jimenez, "Static and Dynamic behaviour of pile supported 
structures in soft soil," Ph.D. dissertation, Université Grenoble Alpes, 
2019. 

[32] J. L. Rangel-Nunez, A. Gomez-Bernal, J. Aguirre-Gonzalez, E. Sordo-
Zabay, and E. Ibarra-Razo, "Dynamic Response of Soft Soil Deposits 
Improved with Rigid Inclusions," presented at the The 14 th World 
Conference on Earthquake Engineering, Beijing, China, Oct. 2008. 

[33] M. Á. Mánica Malcom, E. Ovando-Shelley, and E. Botero Jaramillo, 
"Numerical Study of the Seismic Behavior of Rigid Inclusions in Soft 
Mexico City Clay," Journal of Earthquake Engineering, vol. 20, no. 3, 
pp. 447–475, Apr. 2016, https://doi.org/10.1080/13632469.2015. 
1085462. 

[34] M. Houda, O. Jenck, and F. Emeriault, "Soft soil improvement by rigid 
inclusions under vertical cyclic loading: numerical back analysis," 
European Journal of Environmental and Civil Engineering, vol. 25, no. 
3, pp. 409–428, Feb. 2021, https://doi.org/10.1080/19648189.2018. 
1531268. 

[35] R. M. Patel, B. R. Jayalekshmi, and R. Shivashankar, "Effect of 
Reinforcement Width on Dynamic response of Basal Geosynthetic-
Reinforced Embankment," Transportation Infrastructure 
Geotechnology, vol. 9, no. 4, pp. 516–542, Aug. 2022, https://doi.org/ 
10.1007/s40515-021-00188-1. 



Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9843-9848 9848 
 

www.etasr.com Bouabdallah et al.: The Effect of Rigid Inclusions on the Dynamic Response of Highway Embankment 

 

[36] S. C. Edem, M. Padmavathi, and V. Padmavathi, "Dynamic Response 
Analysis of Highway Embankment with Different Fill Material 
Modifications," International Journal of Latest Technology in 
Engineering, Management & Applied Science, vol. 4, no. 10, 2015. 

[37] R. Hadidi, Y. Moriwaki, J. Barneich, R. Kirby, and M. Mooers, "Seismic 
Deformation Evaluation Of Lenihan Dam under 1989 Loma Prieta 
Earthquake," presented at the Tenth U.S. National Conference on 
Earthquake Engineering, Anchorage, Alaska, Jul. 2014. 

[38] F. R. Olaya and L. M. Cañabi, "Seismic Assessment of a Dam on a 
Clayey Foundation," in Proceedings of the 4th International Conference 
on Performance Based Design in Earthquake Geotechnical Engineering 
(Beijing 2022), Cham, 2022, pp. 1959–1967, 
https://doi.org/10.1007/978-3-031-11898-2_179. 

[39] A. Osouli and S. Zamiran, "The effect of backfill cohesion on seismic 
response of cantilever retaining walls using fully dynamic analysis," 
Computers and Geotechnics, vol. 89, pp. 143–152, Sep. 2017, 
https://doi.org/10.1016/j.compgeo.2017.04.007. 

[40] J. Macedo, P. Torres, L. Vergaray, S. Paihua, and C. Arnold, "Dynamic 
effective stress analysis of a centreline tailings dam under subduction 
earthquakes," Proceedings of the Institution of Civil Engineers - 
Geotechnical Engineering, vol. 175, no. 2, pp. 224–246, Apr. 2022, 
https://doi.org/10.1680/jgeen.21.00017a. 

[41] E. Naesgaard, P. Byrne, and A. Amini, "Hysteretic model for non-
liquefiable soils (UBCHYST5d)," UBCHYST5d Memo, 2011. 

[42] R. L. Kuhlemeyer and J. Lysmer, "Finite Element Method Accuracy for 
Wave Propagation Problems," Journal of the Soil Mechanics and 
Foundations Division, vol. 99, no. 5, pp. 421–427, May 1973, 
https://doi.org/10.1061/JSFEAQ.0001885. 

[43] G. Candia and N. Sitar, "Seismic Earth Pressures on Retaining 
Structures in Cohesive Soils," Berkeley, CA, USA, UCB GT 13-02, 
Aug. 2013. 

[44] M. B. Darendeli, "Development of a new family of normalized modulus 
reduction and material damping curves," Ph.D. dissertation, University 
of Texas, Austin, TX, USA, 2001. 

[45] "FLAC2D." Itasca Consulting Group, Inc., Minneapolis, MN, USA. 

[46] L. Hazzar, M. N. Hussien, and M. Karray, "Influence of vertical loads on 
lateral response of pile foundations in sands and clays," Journal of Rock 
Mechanics and Geotechnical Engineering, vol. 9, no. 2, pp. 291–304, 
Apr. 2017, https://doi.org/10.1016/j.jrmge.2016.09.002.