Validation of Using Modulus of Subgrade Reaction To Consider The Soil Structure Interaction (Farouk 2015)
Validation of Using Modulus of Subgrade Reaction To Consider The Soil Structure Interaction (Farouk 2015)
1
Founder and Coordinator, Soil-Structure Interaction Group in Egypt (SSIGE), Non-
Profit Organization, Nasr City, Cairo, Egypt. E-mail: [email protected]
2
Chairman of EHE Consulting Group, Nasr City, Cairo, Egypt; and Ajman, UAE. E-
mail: [email protected]
Abstract
This paper presents the effect of the soil structure interaction on the modulus of
subgrade reaction. Many of the structural designers represent the footings as a hinged
support in the structural models, while the other designers represent the footings as a
group of springs with spring coefficient K (kN/m). Designers assumed this modulus
by different equations that are not considering the soil structure interaction. This
paper showed that the soil structure interaction has a significant effect on the modulus
of subgrade reaction. The subgrade reaction modulus is not uniformly distributed
under the footings. Using the soil structure interaction in the calculation of the normal
forces in columns results in a relatively change. Plane two-bay frames and two-bay by
two-bay frames were studied using geotechnical and structural, 2D and 3D finite
element programs. It is investigated for the plane frames that the normal force in the
inner wall is decreased and it is increased for the outer walls. For the three
dimensional frames, the inner column’s normal force is decreased and the edge
column’s normal force is increased with the same normal force in the corner walls.
For the complex structural systems, the distribution cannot be estimated, so that the
use of geotechnical finite element programs can consider the soil, foundation, and
superstructure interaction, which is highly recommended for the practical and
researcher engineers.
INTRODUCTION
The development of modern cities within large urban areas with limited
surface space, has led to an in-crease in the rate of construction of high rise buildings
with complex architectural design. The foundation of such complex buildings
presents a geotechnical challenge where the soil structure interaction plays an
important role to achieve the most economical design that satisfies all safety and
serviceability requirements. The cooperation between both geotechnical and structural
engineers is necessary to reach a successful design. The subgrade reaction modulus
"ks" can be considered as an appropriate interface between the geotechnical and
© ASCE
AEI 2015
AEI 2015 639
structural engineers. It is well known that the subgrade reaction modulus is not a soil
constant, but it depends on many factors such as the dimensions of footings, soil
conditions, load levels, and superstructure rigidity. There are different ways to
calculate the modulus of the subgrade reaction. These ways are mainly based on the
Winkler model.
as a bed of evenly spaced, independent, linear springs, as shown in Figure (1). It also
assumes that each spring deforms in response to the vertical stress applied directly to
that spring, and does not transmit any shear stress to the adjacent springs, however, in
real soils the displacement distribution is continuous. The deflection under a load can
occur beyond the edge of the slab and the deflection diminished at some finite
distance.
The elastic constant of the springs K (kN/m) is taken equal to the area of
coefficient of subgrade of the soil ks (kN/m2/m). The springs assumed to resist both
compression and tension. The intensity of soil reaction and the vertical deflection of a
raft or foundation beams at every point are related by the following equation:
Pi = K * yi …………………. (1)
Where:
Pi = intensity of soil reaction at point (i)
K = spring’s coefficient
yi = vertical deflection at point (i)
Tension
© ASCE
AEI 2015
AEI 2015 640
Winkler model. Pasternak considered the shear interactions between the spring
elements of Winkler foundation by connecting the ends of the springs with a beam or
plate consisting of incompressible vertical elements, which are deformed only by
transverse shear, as shown in Figure (3).
Shear Layer
Vlasov Model (1966). There is no value can be adopted for ks from the
previous modification models. Vlasov and Leont'ev (1966), Figure (4), tried to use a
new mathematical approach to solve the above problem. They developed a two-
parameter model for plates on elastic foundations using a variational method. This
method accounted for the neglected shear strain within the soil continuum. Vlasov
and Leont'ev introduced another parameter, γ, to characterize the vertical deformation
profile within the soil continuum. Vlasov and Leont'ev assumed a displacement
pattern of the soil medium by using an arbitrary value of, γ. The model is provided for
the effects of the neglected shear strain energy in the soil and the subsequent shear
forces on the plate edges are because of the soil displacement.
© ASCE
AEI 2015
AEI 2015 641
uniformly distributed loads, they found that the value of γ is also dependent on the
ratio between the depth of the soil and the length of the beam. They called their
model a modified Vlasov model. All these parameters depended on the material
properties of the soil, soil structure, geometry of the structure, depth of the stratum,
and the distribution of the loading on the plate. They used the finite-difference
method for the calculation. Straughan (1990) applied the modified Vlasov model to
the analysis of rectangular slabs on elastic foundations. He used the finite-difference
Downloaded from ascelibrary.org by University of Saskatchewan on 03/19/15. Copyright ASCE. For personal use only; all rights reserved.
technique to solve the field equation for several loading cases. Recent works by
Vallabhan and Das (1987, 1988, and 1989) strengthened the Vlasov postulation for
beams on elastic foundations, but stopped short of developing computational
techniques for plates.
Due to the mathematical complexity, the improved soil models are less
attractive than the Winkler model. Some practical designers whom consider the soil
structure interaction in the structural design are still using the Winkler model to
represent the soil. The Winkler model does not transmit any shear stress outside the
footing. This was the main disadvantage for this model as stated in many references.
In our opinion, there are some other factors that affect the modulus of subgrade
reaction. Some of these factors were discussed before by the authors in other
researches such as: the foundation width, the foundation thickness, the type of soil. It
was proved that, the subgrade reaction modulus (ks) is not a soil constant and not
uniformly distributed under the footing. The distribution of the modulus is mainly
based on the foundation thickness and width, i.e. foundation rigidity. If the footing is
flexible, the subgrade is concentrated under the columns and low at edges. If the
footing is rigid, the modulus is concentrated near edges and low under the columns.
Therefore, using of an average value of the subgrade affects the resulted bending
moment on the footings. Sometimes, the practical engineers use the plate loading test
to find the average value for the subgrade modulus and they correct the result by the
width factor only. It was proved by authors that, the plate materials and dimensions
have a significant effect on the plate rigidity factor that affects the test results. The
authors preferred to use a new modification factor for the subgrade reaction to
consider the foundation rigidity. It can be done by using the numerical analysis in the
calculation of the subgrade modulus with soil models that can consider the soil
plasticity. The numerical analysis can be performed using a geotechnical finite
element program that can model the footing with the underlying soil.
© ASCE
AEI 2015
AEI 2015 642
Two frame-soil models are studied. Model one is modeled using the 2D
geotechnical and structural finite element programs, while Model two is modeled
using the 3D geotechnical and structural finite element programs.
The frame has bay span of 10.00 meters and 7.00 meters height. The inner footing is
4.00 meters breadth and 1.00 meter thickness, and the outer footings are 2.00 meters
breadth and 1.00-meter thickness. The foundation level is assumed at the ground
surface with no embedment depth to neglect the effect of the embedment depth on the
results. The outer and inner walls are 0.45 meters thickness. Rotations of the walls are
permitted at foundation level. The slab has thickness of 0.60 meters with applied
distributed load of 40.0 kPa. Self-weights of walls and slabs are neglected. The soil is
sand with Young’s modulus (Es) of 20.0 MPa, Poisson’s ratio is 0.3, and the soil’s
unit weight is 18.0 kN/m3. The frame is modelled using the geotechnical finite
element program PLAXIS 2D-AE with modelling the soil by the linear elastic model,
as shown in Figure (5-a). The frame also modelled using the structural finite element
program SAP2000 by representing the footings as hinged supports, as shown in
Figure (5-b). These two representations are to show the effect of the soil structure
interaction on the resulted normal forces in the frame walls. The normal forces in
walls influence the resulted bending moment on frame slabs. It also affects the design
of the foundations.
The normal forces (NF) in outer and inner walls are determined from PLAXIS
and SAP2000 models. Table (1) illustrates values of normal forces (NF) for the
corner and inner walls resulted from the two programs. Difference ratio (∆) is
depicted for the normal forces, (∆ = 100 – (NFSAP /NFPLAXIS *100)). It is the
difference between considering and neglecting the soil structure interaction in the
analysis of the superstructure members.
It is investigated that the soil structure interaction leads to increase the normal
forces in the outer walls and decrease that in the inner wall. Therefore, the resulted
settlement under the outer walls should be increased and the inner settlement should
be decreased. Hence, the differential settlement is decreased without any additional
members; i.e. tie beams that connect the footings. We can say that, the soil structure
interaction affects the superstructure and the geotechnical analyses. The distributions
of the subgrade modulus are shown in Figure (6- a, b).
© ASCE
AEI 2015
AEI 2015 643
affects also the distribution and the values of the modulus of subgrade reactions (ks)
under the inner and outer footings.
The use of plate loading test results a uniformly distributed and unique value
for the modulus of subgrade reactions under the footings. The value that is resulted
from the test is then corrected by the width and rigidity factors. The factors are based
on the footings widths and thicknesses. The footings dimensions are calculated using
the applied forces from the walls. The applied forces in the walls mainly depend on
Downloaded from ascelibrary.org by University of Saskatchewan on 03/19/15. Copyright ASCE. For personal use only; all rights reserved.
the soil structure interaction as showed before. The soil structure interaction cannot be
correctly performed, because the subgrade modulus that is resulted from the plate
loading test is not corrected until the superstructure and the foundations are designed.
This circulation is the main factor that makes some practical designers to neglect the
soil structure interaction. The circulation can be stopped by using a modification
factor for the plate loading test by using a finite element geotechnical program that
can model the superstructures, foundations, and the underlying soils.
(a) PLAXIS 2D
© ASCE
AEI 2015
AEI 2015 644
Table 1. Values of normal forces in walls (NF) and the percentage of difference
(∆ %).
A new advanced design procedure is suggested and could be used for the
special projects to consider the superstructure rigidity in the geotechnical analysis.
The procedure is as the following:
1- The geotechnical engineers should present a preliminary modulus of
subgrade using the plate loading tests after the modifications.
2- The structural engineers should finish the preliminary design based on
the ks suggested in the geotechnical report.
3- The preliminary structural design should be returned to the
geotechnical engineers.
© ASCE
AEI 2015
AEI 2015 645
………… (2)
Where H is the half width of the building in the plane of deformation and EI is
the bending rigidity of the structure. EI can be calculated using the following
Equation.
………. (3)
© ASCE
AEI 2015
AEI 2015 646
Figure (8) shows the relation between the rigidity factor and the resulted average
spring’s coefficients.
Downloaded from ascelibrary.org by University of Saskatchewan on 03/19/15. Copyright ASCE. For personal use only; all rights reserved.
It can be seen that, the superstructure has a significant effect on the subgrade
modulus under the inner and corner footings and has no effect on the modulus under
external footings. In addition, the applied load levels affect the resulted values of the
average subgrade modulus. The footings dimensions are the same for the different
load levels and for the different structural rigidities. Nevertheless, the subgrade
© ASCE
AEI 2015
AEI 2015 647
modulus is not constant under the footings at the same load level. This is the result of
the superstructure rigidity effect.
The resulted normal forces (NF) in the columns from PLAXIS and SAFE are
collected and presented in Figure (10) for the load level of 10.00 kPa. It can be seen
that the inner column is significantly affected. The corner columns also are affected
and the edge columns approximately are not affected.
© ASCE
AEI 2015
AEI 2015 648
Figure (11). Validation of the suggested procedure for the cooperation between
the geotechnical and structural engineers
The resulted normal forces in walls from considering soil structure interaction
by SAFE and PLAXIS are closed to other. There are small differences because of
using an average value for the subgrade modulus in the structural program.
CONCLUSIONS
© ASCE
AEI 2015
AEI 2015 649
RECOMMENDATIONS
ACKNOWLEDGEMENTS
Praise and glory are due to Allah whom grants us all our knowledge and
success in our life. We would like to express sincere appreciation to the team of EHE-
Consulting Group on their help and encouragement. Furthermore, the researchers are
grateful to the Soil-Structure Interaction Group in Egypt (SSIGE) for their continuous
support to the design process.
REFERENCES
© ASCE
AEI 2015
AEI 2015 650
© ASCE
AEI 2015