 Open Access
 Total Downloads : 0
 Authors : Er. Prabhat Kumar Jha
 Paper ID : IJERTV7IS080084
 Volume & Issue : Volume 07, Issue 08 (August – 2018)
 Published (First Online): 05012019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Simulation of Triaxial Compressional Behaviour of Sand based on State Concept
Simulation of Triaxial Compressional Behaviour of Sand based on State Concept
Er.Prabhat Kumar Jha1 Department of Roads Nepal
Abstract An elasticplastic constitutive model, which account for the combined influence of density and effective stress in a unique way, is developed for the triaxial compressional behaviour of sand under monotonic loading. The research addresses many important issues relevant to sandy soil behaviour and represents a unite experimental and theoretical effort to develop a constitutive model for sandy soils comprising balanced features of accuracy, simplicity and versatility. Particular attention in both experimental and theoretical considerations is given to the effects of density and effective stress as well as their combined influence on the sandy soil behaviour.
A distinct feature of the constitutive model of the research is that a single set of values of the material parameter is used to model the behaviour of the given sand. The parameters of this relation are expressed as function of an index parameter established in the framework of the steady state of deformation concept. The parameter employed in the relation is the State Index, Is which characterizes the sandy soil behaviour by accounting for the combined effects of the density and effective stresses related to a given initial fabric.
The elasticplastic constitutive model is developed in the framework of the incremental theory of plasticity. The framework is developed by the modification of the torsional simple shear test framework, in stress space definition, yielding criterion, stressdilatancy relationship, plastic modulus, and loading index. The model is defined in a stress space that enables to account failure surface which incorporates the effects of density and effective stress.
The accuracy and effectiveness of the elasticplastic constitutive model is assessed through a comparison of the measured and predicted behaviour of the Bagmati sand as well as the Toyoura sand in monotonic triaxial test. The model correctly predicts the stressstrain characteristics observed in the undrained triaxial compression tests on medium dense, loose, very loose sand, providing gradual change in these characteristics from those typical for dense sand, with continuous increase of the shear stress with straining, to those typical for loose sand, with drop in the shear stress associated with developmenat of the large strains including the extreme cases with zero residual strength.
Keywords Constitutive Model; Sandy Soil Behaviour; State Index; Bagmati Sand

INTRODUCTION
The need to solve increasingly sophisticated soil engineering and soilstructure interaction problems has created a demand for better experimental devices and techniques as well as better analytical tools, in order to have more clear understanding and representation of soil behavior. In the recent years, considerable experimental and theoretical efforts have been made to research and model the soil behaviour under various loading conditions, including cyclic loads as caused by earthquakes, traffic and sea wave loads. On the
other hand, intensive research has also been done on its monotonic behavior, providing abundant valuable information on the essence of soil behaviour from very small strains with nearly elastic behaviour up to large deformations and flow of the material.
The rapid development in the recent years of advanced computation procedures for nonlinear stress analysis, such as the finite element method, especially emphasized the need for appropriate constitutive models, undoubtedly a crucial element in the numerical analysis. The development of more advanced constitutive models for soils has also been stimulated and thoroughly supported by the increasing supply of high quality experimental evidence and a continuous widening of our knowledge of soil behavior. A considerable challenge is posed, however, by the fact that soil behaviour is remarkably complex, involving various complicated processes and phenomena both under monotonic and cyclic loading.
Constitutive modeling of soils confronts two opposing requirements: on the one hand, the model should be able to describe the extreme complexity of soil behavior, and on the other hand, it should be comprehensible and simple enough to assure easy implementation in numerical analysis and application to practical problems. Whether complexity or simplicity will prevail is the first important decision to be made in the process of the development of a constitutive model. Yet, this judgment and the choice of the relevant aspects of the soil behaviour that should be modeled is probably the most critical for the accuracy, effectiveness and applicability of the model.
It is general trend in elasticplastic sand models, however, that neither density nor the combined influence of density and normal stress is considered as a variable. In fact, a sand with different densities is regarded as a different material, and the effects of the normal stress are considered to be independent of the density. Roscoe and Poorooshasb (1963) were the first to emphasize the importance of the combined influence of density and normal stress on soil behavior. On the basis of both theoretical and experimental considerations, Roscoe and Poorooshasb introduced the critical state of the soil as a reference for quantification of the effects of the initial state on soil behavior. Here, the initial state refers to the void ratio and the mean normal stress of the soil prior to shearing. In a similar vein, Been and Jefferies (1985), Bolton (1986), Ishihara (1993) and Verdugo (1992) have reported extensive studies on the combined influence of the initial density and normal stress on sand behaviour. A distinctive feature of these studies is that particular indices and parameters serving as measures for the effects of the initial state on sand behaviour
have been proposed. The proposed relative Dilatancy index Ir (Bolton, 1986), the state parameter (Been and Jefferies, 1985) and the state index Is (Ishihara, 1993; Verdugo, 1992) integrate the effects of the density and normal stress into a single parameter and as such, they are different from the commonly used parameter, the relative density Dr .Recently, several attempts have been made to use the state parameter and the state index Is for sand modelling (Jefferies, 1993; Cubrinivski, 1993; Wood et al., 1994). Jefferies (1993) incorporated the state parameter as a hardening parameter in the Cam clay model to describe monotonic behaviour of sand and Wood et al. (1994) used the same parameter for modelling drained monotonic stressstrain behaviour associated with strain softening. Ishihara and Cubrinovski ( 1998 ) have simulated the undrained behaviour of sand using their modified hyperbolic relation along with the stress Dilatancy relation based on the energy approach of Roscoe et al 1963, and Gutierrez ( 1993 ) flow rule, in plane strain case.
The research deals with the sandy soils and presents an integrated experimental and theoretical investigation on their behaviour under monotonic loading conditions. The research contains all the necessary steps in evaluation of a constitutive model for sandy soils, including experimental research where the experimental basis for the model is produced, then development of a stressstrain relation, formulation of an elasticplastic constitutive model , verification of the model. The present research model follows the framework of the state concept, in which a relative initial state with respect to some eference states of sand (Critical or Steady state; Quasisteady state) is used for characterization of sand behavior.
Constitutive Modelling
Without aiming to review present constitutive models for sandy soils at this point, a brief discussion on the general trends in the modeling and basic requirements that a constitutive model has to fulfill will be given. To start with, it is important to recognize the existence and the need of various types of models since there are considerable difference as regards the purpose of the model, the problem that it considers and the required level of accuracy and simplicity among others. The development of the advanced constitutive models has already been under way for several decades, and these models have passed well beyond the speculative or qualitative description of soil behavior. During this period, a great number of constitutive models have been proposed involving various theories and concepts. Those concepts are either of a micromechanical or phenomenological nature. Even tough the principal feature of soils is their particulate structure and thus the ultimate adequate representation of the true nature of soils lies in the understanding of the microscopic behavior, the complexity of these processes and the diversity of the soil field conditions will definitely restrict the quantitative predictions of micromechanical models and especially their applicability in the foreseeable future. Therefore, it is quite understandable that the majority of the present advanced constitutive models rely on phenomenological or macroscopic observations of soil behavior. Only, these will be of interest in the present work.
Recent development in the constitutive modeling of the sandy soils and characterized by usage of a variety of approaches and theories including elasticity, viscosity and plasticity, as well as combinations of these theories: viscoelasticity, visco plasticity and elastoplasticity. This trend in the modeling is simply due to the fact that soils exhibit behaviour related to these theories. However, as pointed out earlier, it is crucial to select and model only the important aspects of the sandy soil behaviour that a constitutive model should cover, since that is the only approach which will contribute to the simplicity, comprehensibility, and applicability. That this has not been always done in frequent constitutive models is clearly indicated by the comment of Scott (1988) on the constitutive models presented at the workshop on Constitutive Equations for Granular Non Cohesive Soils (Saada and Bianchini
,1988).
Clearly, we cannot go on like this, with the ultimate goal of achieving a thermoviscoelastic hypoplastic multiple critical bounding capped yielding endochronic state surface models. At some stage, physical intuition has to reenter the problem.
It has to be recognized , on the other hand, that plasticity models are already proved to be capable of capturing the important features of soil behaviour such as dilatancy, nonlinear and irrecoverable deformations, dependency on stress path and inherent and induced anisotropy among others. In fact, plasticity theory has been the cornerstone of the constitutive modeling of soils for more than two decades. The balance between the simplicity and accuracy of a constitutive model seems to be more or less achievable within the framework of the incremental theory of plasticity.
The requirements that a particular model has to fulfill are strongly dependent on the purpose of the model and in particular on the problem that it aims to solve. However, in general, the model has to satisfy theoretical, experimental and computational requirements. In other words, besides the requirements for theoretical rigors, it is essential to model the soil behaviour on a physical relevant basis in order to attain a realistic representation of the soil behaviour and versatile model. No less important is the fact that the number of the model parameters should be reasonable and that their experimental evaluation can be easily assessed in conventional tests. Computational efficiency and memory requirements are certainly features that can seriously influence the applicability of a constitutive model.
Sandy Soil Behavior
The behaviour of the granular materials can normally be accurately characterized for the practical application from the results of the drained tests, since in many problem, the material can usually be considered fully drained. However, if the rate of loading is very high, such as those resulting from explosions or earthquakes, or the permeability of the cohesionless soil is relatively low, as in fine sands or silts, undrained conditions can exists. If the material possesses low relative density, pore pressure increments, may result from these undrained conditions. This in turn may lead to instability, liquefaction, and failure earth dam and submarine. Slope failure has been attributed to the development of
undrained condition in loose, fine sands and silts. These catastrophic failures illustrate the importance of understanding the undrained behaviour of sand.
Undrained monotonic loading of saturated cohesionless materials has been primarily used to investigate the mechanics of the soil instability and liquefaction. The requirement for soil instability and liquefaction are that a saturated, cohesionless soil exhibits significant contractive volume change tendencies during the undrained shearing with increase in the positive excess pore pressure such that continuously decreasing effective stress results as shearing progress. This can be achieved by shearing very loose soil at low as well as high initial confining pressure, but it can also be realized by shearing the dense specimen at high initial confining pressure. When sand is subjected to a monotonic or cyclic shear load it exhibits very complex behaviour which is governed by the coupling between the shear and the volumetric strains. The link between the distortional and volumetric strains is certainly the most distinctive feature of the deformation characteristics of granular materials. In accordance with this feature the behaviour of sand is profoundly affected by the initial density and fabric of packing, as well as by the initial normal and shear stresses. In other words, given sand behaves either in contractive or dilative manner and shows different cyclic strength or resistance to liquefaction depending on the initial
Most probably, the worldwide use of the triaxial test in geotechnical engineering is because of some of the following advantages of this laboratory test (Baldi et al.1986): relative simplicity of drainage control and measurement of pore pressure, ability to apply principal stresses in known directions, ease of measuring axial and volumetric strains, use of solid cylindrical specimens which can be conveniently obtained from standard tube samples or easily trimmed from block samples, and versatility of the equipment which may be used for a variety of determinations besides triaxial strength and stiffness such as for instance, consolidation and permeability parameters, wave velocity, dynamic parameters and so forth. On the other hand, the conventional triaxial test has some limitations. Perhaps, the two most important are: the principal stresses can not be rotated continuously, but only by a sudden rotation of 900, and it is not possible to have a truly generalized stress state, but at least, two principal stresses must to be identical. Consequently, the state of stresses is always axially symmetric.
State of Stress Conventional TriAxial Shearing Mode Tests
1 0 0
2
0 0
packing and stress conditions.
Triaxial Shearing
0 0
1,1
3
0 0
Without any doubt, the triaxial test on cylindrical specimen is presently the most widely used procedure to evaluate several important properties of the soil behavior such as: strength, dilatancy, stressstran relationship and anisotropy among others. The triaxial test has been used for research as well as
i, j
0
0
1 0
0
0
0 0
0
for design projects and a very significant part of the present knowledge about soil behavior comes from triaxial test results.
In the conventional triaxial test, the sample is cylindrical in shape, laterally surrounded by a rubber membrane and vertically enclosed by rigid end platens through which the axial loading is applied. This configuration immediately creates specific boundary conditions that, as in any other laboratory test device, are different from the field conditions. This shortcoming has been accepted due to the fact that the
0
0
1,1
0
i, j
0
2 0
0 3
0 0
0
2,2
0 3,3
field conditions cannot be exactly known nor can they be completed simulated by any laboratory test. In this sense, laboratory measurements have always to be seen an approximation of what may actually occur in the field.

RESEARCH METHODOLOGY

StressStrain Representation
A major simplification applied to the constitutive relation to be developed is the limitation of the formulation to two dimensions. The cylindrical specimen of the conventional triaxial test is in axisymmetric stress state ( fig.1 ). It is considered that the assumption of a two dimensional condition will avoid complicating the formulation significantly and yet on the other hand, will be appropriate for idealization of many insitu loading condition within reasonable computational efforts and required accuracy.

Incremental Formulation
1 2 1
The model is formulated in the XYp space (fig. 2), where the
p 1 X Y
(7)
stress variables are defined as; ( equ.1 , 2 by K,H, Roscoe,
3 3
1
1.5
1
q 1 X Y
(8)
3 6
For the convention triaxial test with axisymmetric stress state,
dq
(9)
2 d
3 11

d
33
Since , in case of undrained shearing , d v is zero, which suggests
d q = d11 , i.e. axial strain is similar to the shear strain. The failure surface, f ( X,Y ) = 0 can be expressed as
A.N. Schofield, 1961 )
2 2 0
(10)
Fig. 1 State of stress, strain in triaxial compression test
1
f X 2 Y 2

rf

To start , it is assumed that the total strain increment
dij is
decomposed into elastic and plastic components d e and d p by a
simple superposition ,
ij ij
ij
d
ij
dij = d e + p
(11)
The elastic and plastic strain increment components are calculated separately, by an elastic stressstrain law of the form,
d ij
E d
e
ijkl kl
(12)
Fig.2 Typical shape of Failure surface in XYp space ( K,H, Roscoe,
A.N. Schofield, 1961 )
1
Where,
e
2G
X 1 2 3
3
(1 )
Eijkl
1 2
ij
kl
Ge ik
jl
il
jk
(13)
Y 2
1 3
(2)
is Poissons Ratio,
3
2
ij
= Kronecker Delta ( ij =0 if i
j , and = 1 for i
j ),
X 2 Y 2 p2 3
2
(3)
Ge is Elastic Shear Modulus, dependence on the epplane and
3
empirically defined as ( Iwasaki and Tatsuoka, 1977 )
2
2.17 2 n
r p
3
2
(4)
G A p
e p
(14)
3
e a 1 e
pa
r p
3
2 2
(5)
Where, A is dimensionless parameter, pa
is atmospheric pressure
f 3
p
max
1104 Kpa, p is the mean normal stress, n is the dimensionless
material parameter determining the variation of Ge with p. Thus, for
EN . q . max
triaxial case,
, (6)
G 2 3
q
max EN . p
d11 e d11 1 2
, d 22 0,
d 33 0,
d12 0
(15)
The incremental plastic relations are given by


Material Parameters Determination
A distinct feature of the constitutive model of the present
ij
d p
L g
ij
(16)
study is that a single set of values of the material parameter is used to model the behaviour of the given sand. That is whatever is the relative densities of the sand and the confining
Where, direction of plastic strain increment is specified as the normal to the plastic potential, g whereas the magnitude of plastic strain increment is given by the loading index , L.
Finally, the general expression for the plastic strain increment can be written as,
stress, the values of the material parameters remain the same. Thus, it is a prerequisite to use a determination procedure and to evaluate material parameters which will able resonablely accurately overall representation of a given sand. In general, except for the elastic parameters, the model parameters are dependent on the fabric of packing. And, except the mean
p 1
f
g
normal pressure of the reference lines, all the parameters are dimensionless. The model requires four groups of parameters.
d ij
HP ij
ij
d ij
(17)

Elastic parameters.

State lines (reference states).
The total strain increment can be calculated as

Stress strain parameters.
d ij
= E 1 e
ijkl d kl +
1
H
f

d
ij
g
(18)

Dilatancy parameters.
P ij
ij

Elastic parameters.
There is no specific requirement for the evaluation of the elastic parameters, and insistence on determination of these parameters at
Since, our experiment is strain controlled, so its necessary to model a strain controlled simulation.
very small strain levels is not necessary since modeling of small strain behaviour is not attempted herein. Moreover, the Poissons ratio and the exponent, n may be assigned fixed value of 0.2 0.3 and
0.5 0.6 respectively. The shear constant, A is back calculated from
1
f
g
known value of elastic shear modulus, Ge ; equation 14.
d ij Eijkl d kl
HP
ij



d ij
ij
(19)
The empirical relation for Ge in undrained case is given by (Hicher,
The calculation of the plastic strain increment components in
the framework of the incremental theory of plasticity is based
PierreYves, 1994 );
on definition of the yield function, plastic potential function, and hardening function. The stressdilatancy relationship for the case of triaxial shearing, given by Dafalias (1999) is used
Ge
450 p'
e
1 u
1
(22)
in the formulation.
The loading index, L is
Where, u is the Poissons ratio in undrained condition (=0.5) for perfectly saturated soil.

Reference states.
L 1 XC
(20)
The state lines, the URlines and QSSline, will be evaluated using
And,
2 HP
r f
results of several monotonic undrained strain controlled tests on loose samples having different void ratio with theonly requirement being that in the course of the undrained shearing, the samples will exhibit drop in the shear stress or QSS. Assembling the quasi steady states
d p
2
from each test in the ep diagram, QSS line will be drawn. If during
H E

f E E
r p
the test of undisturbed soil samples, the requirement for drop in the
p N
0.01 N
N ,min 1
r f
(21)
shear stress often may not be satisfied, and in such case it will be necessary to replace the QSSline with SSline as a reference line.
Again, if the zero steady state cannot be directly observed in the undrained tests, the position of the URline ( threshold void ratio )
Where, HP
is plastic modulus, and l is loading surface.
will be reasonably approximated in the ep plane by extending either
Along the loading surface, plastic modulus is constant. The loading surface which is moving in the stress space, ij , according to the load increments, is used to evaluate the current magnitude of the plastic modulus as well as to trace the
the QSS/SS line to the point of the intersection with the void ratio
axis. This point of intersection then denotes the void ratio eU .

Stressstrain parameters
The stress strain parameters of the model consists of linear
loading intensity. It has to notice that the
GN is the initial
coefficients of the expressions,
modulus whereas, HP is the tangential plastic modulus.
max
a1 b1Is
, Is 0
EN ,max a2 b2 Is
EN ,min a3 b3Is
, Is 0
, Is 0
(23)
The coefficients, a1,b1, a2 ,b2 , a3,b3 will be determined using results from a series of drained pconstant tests on samples

Dilatancy parameters
For the evaluation of dilatancy parameters, the critical state
with distinct relative initial states or samples having different
stress ratio,
M q
will be first determined from the
values of the state index.
p
cs
The peak stress ratio, max
of each test will directly read off
aforementioned monotonic undrained tests or alternately it
may be approximated by the peak stress ratio of the drained
from the measured normalized stress strain curve and the
tests, for
Is 1, or the final slope or slope at large strains of
initial plastic moduli,
EN ,max and
EN ,min
will be determined
the curve, normalized plastic shear work versus plastic shear
either as the best fit small strain and large strain values
respectively, in trial and error simulations of the experimental data or by establishing an experimental relationship between
p
strain from drained tests. The parameter o is the initial slope
or slope at very small strains of the curve, normalized plastic shear work versus plastic shear strain from drained tests. The
EN and
d q ( Cubrinovski and Ishihara, 1998 ). For each
dilatancy parameter, Sc
which controls the pore pressure
test, the state index
Is will be calculated using the reference
development and has an especially large influence on the
line and the initial ep state of the sample, and the values of
max , EN ,max and EN ,min as defined above will be plotted
cyclic strength. The parameter will be determined by
mathematical corroboration due to absence of cyclic triaxial equipment in the laboratory, or can be determined as the value
verses
Is ; these correlation with
Is will be approximated
of the plastic strain at which average value of is attained.
with linear relationship and the
For a larger value of Sc
reflects smaller at a given strain
coefficients, a1, b1, a2 , b2 , a3, b3
these data arrangements.
Computational Flow Chart
will be determined from
or less contractive behavior. The material parameters for the Bagmati sand will be used to model the monotonic behaviour of samples with various combinations of initial relative densities and confining pressure.
2.4 Physical Parameters of the Bagmati Sand (Bagmati River within Kathmandu Valley)
Maximum Void Ratio, emax 0.946 Minimum Void Ratio, emin 0.610 Effective Size, D50 0.6 mm Coeffi. of Uniformity, CU 2.65 Coeffi of Curvature, CC 1.12
% Fines less than 1 % Specific Gravity, Gs 2.68
Source: Jha,2004


RESULTS AND EVALUATION

Material Parameters of the model for Bagmati Sand

Elastic Parameters
Since, small strain behaviour is not attempted; there is no specific requirement for the evaluation of the elastic parameters. The Poissons ratio, and the exponent, n are fixed to be the value of 0.25 and 0.55 respectively. And, the shear constant, A is computed by the back calculation of equ.14 and 22, and found to be 245.31.

Reference Line

3.1.4 Dilatancy Parameters
The normalized shear work has been computed for the Triaxial case for confining pressures 100Kpa, 200 Kpa, 400 Kpa, and 500 Kpa, and plotted for corresponding plastic shear strain up to 27 %. The plot has justified the relation that the normalized work for various stress paths are approximately proportional to the shear strengths along those stress path ( fig.7 ).
700
600
500
EN,max = 50 Is + 150.0
Fig. 3 Characteristic Curves of Bagmati Sand, ( Author, Karn,2005, Jha,2004 )
The characteristic curves of Bagmati Sand is shown in fig. 3, representing the loosest isotropic consolidation line ( IC Line
), the steady state lines ( SS Line ) and the upper reference line ( UR Line ). The SS line is best fitted from various set
EN,max , ( Kpa )
400
300
200
100
0
0 1 2 3 4 5 6 7 8 9
State Index, Is
of void ratios and effective mean pressures at steady state of large number of drained and few undrained tests. These curves has been extended up to the effective mean pressure of 500 Kpa. The upper reference line is horizontal up to 50 Kpa, and
Fig. 5 Relationship between stressstrain parameters and state index for Bagmati Sand : Normalized Initial Elastic Modulus, EN,max (Source:Author, Karn,2005, Jha,2004 )
EN,min = 13.20 Is + 68.57
after 50 Kpa, the IC line serves as the UR line. The IC –
line almost follows parallel to the SSline. But, the UR line for Toyoura sand is horizontal up to 400 Kpa. These lines are used for the computation of the state index.
3.1.3 Stress – Strain Parameters
The peak stress ratios and the corresponding tests state index are plotted . Though points are little scattered, it justifies linear relationship between these. The slope of the line ( b1 ) is 0.0205, and the intercept ( a1 ) is 1.39, ( fig. 4).
(q/p)max = 0.0205 Is + 1.39
250
200
EN,min , ( Kpa )
150
100
50
0
0 1 2 3 4 5 6 7 8 9 10 11 12
State Index, Is
1.70
1.65
1.60
Fig. 6 Relationship between stressstrain parameters and state index for Bagmati Sand : Normalized Initial Elastic Modulus, EN,min (Source:Author, Karn,2005, Jha,2004 )
1.55
(q/p)max
1.50
1.45
Normalized Plastic Shear Work, p
0.6
0.5
100 200 400 500
1.40
1.35
0.4
1.30 0.3
1.25
0 1 2 3 4 5 6 7 8 9
State index, Is
0.2
Fig. 4 Relationship between stressstrain parameters and state index for Bagmati Sand : Peak Stress Ratio ( Source:Author, Karn,2005, Jha,2004 )
Similarly, the normalized initial elastic modulus ( maximum and minimum ) and the corresponding tests state index are
0.1
0
0 0.05 0.1 0.15 0.2 0.25 0.3
Plastic Shear Strain, p
plotted ( fig. 5, and 6 ) . Though points are little scattered, it justifies linear relationship between these. The slope of the line representing the relationship between EN,max and state index ( b2) is 50, and the intercept ( a2 ) is 150. And, the slope of the line representing the relationship between EN,min and state index ( b3 ) is 13.20, and the intercept ( a3 ) is 68.57.
Fig. 7 Normalized shear work versus plastic shear strain relationship measured
in monotonic drained Triaxial compression test on dry sample of Bagmati Sand with various initial ep state ( e = 0.682 0.863 ) ( Source :Author, Karn,2005, Jha,2004 )
p
0.09
0.08
0.07
Normalized Plastic Shear Work,
0.06
0.05
0.04
100 200 400 500
450
400
Deviatoric Stress, q ( Kpa )
350
300
250
200
150
Simulation Sc — 0.0409 Simulation Sc — 0.0125 Simulation Sc — 0.0090 Simulation Sc — 0.0445 Lab
0.03
0.02
Âµ
0.01
0
100
50
0
0 50 100 150 200 250 300 350
Effective Mean Pressure , p' ( Kpa )
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
Plastic Shear Strain, p
Fig 8 Normalized shear work versus plastic shear strain relationship measured in monotonic drained Triaxial compression test on dry sample of Bagmati Sand enlarge view for small strain range (Source :Author, Karn,2005, Jha,2004 )
Fig. 10 Mathematical Corroboration for Sc value for Bagmati Sand at confining pressure 200 Kpa
350
The dilatancy parameter , o
has been determined by the initial
Deviatoric Stress, q ( Kpa )
300
250
Lab Simulation TSP
slope or say, slope at very small strains of the normalized plastic
shear work versus plastic shear strain curve ( fig. 8 ), from drained tests, and found to be 1.057, which is much higher than that of Toyoura sand ( = 0.22 ). The another dilatancy parameters, the critical state stress ratio, M has been first approximated by the peak stress ratio of the drained tests, for Is 1, and found to be 1.42. Again, it has been determined by the value of the final slope of the curve of normalized plastic shear work versus plastic shear strain from drained tests ( fig.7 ), and found to be the same, 1.42. This value is also much
200
150
100
50
0
0 50 100 150 200 250 300
Effective Mean Pressure, p' ( Kpa )
higher than that of Toyoura sand ( = 0.607 ).
400
Simulation Sc — 0.0409
Fig. 11 Comparison of model simulation with undrained Triaxial test results
on Bagmati Sand having void ratio of 0.779 and state index 4.96 at confining pressure 100 Kpa
350
Deviatoric Stress, q ( Kpa )
300
250
200
Simulation Sc — 0.0125 Simulation Sc — 0.0090 Simulation Sc — 0.0445 Lab
350
300
Deviatoric Stress, q ( Kpa )
250
Lab Simulation
150
100
50
0
0 50 100 150 200 250 300
Effective Mean Pressure, p' ( Kpa )
200
150
100
50
Fig. 9 Mathematical Corroboration for Sc value for Bagmati Sand at confining pressure 100 Kpa
And, the final dilatancy parameter, Sc which controls the pore
pressure development has been determined to be 0.0410 as the value of the plastic strain at which average value of is
0
0 2 4 6 8 10 12 14
Strain ( %)
Fig. 12 Comparison of model simulation with undrained Triaxial test results on Bagmati Sand having void ratio of 0.779 and state index 4.96 at confining pressure 100 Kpa
40
attained ( fig. 7 ). Also, the mathematical corroboration for the parameter value determination has been performed with the trial Sc values of 0.009, 0.0125, 0.0409,and 0.0445 ( fig. 9 and 10 ), the predicted effective stress path ( ESP ) corresponding to the value 0.0409 seems to be closer to that of the laboratory results for both the confining pressures, 100 Kpa and 200 Kpa. But, finally the value from laboratory ( = 0.0410 ) is accepted for model simulation.
30
Pore Pressure, u ( Kpa )
20
10
0
10
20
30
40
50
60
Lab Simulation
0 2 4 6 8 10 12 14 16
strain ( % )
Fig. 13 Comparison of model simulation with undrained Triaxial test results on Bagmati Sand having void ratio of 0.779 and state index 4.96 at confining pressure 100 Kpa
400
350
43.8% 20% 10% TSP
350
300
57.1% 20%
10% TSP
Deviatoric Stress, q ( Kpa )
300
250
200
150
Deviatoric Stress, q ( Kpa )
250
200
150
100
50
0
0 50 100 150 200 250 300
Mean Effective Pressure, p' ( Kpa )
100
50
0
0 50 100 150 200 250 300
Mean Effective Pressure, p' ( Kpa )
Fig. 14 Model simulations with undrained Triaxial test results on Bagmati Sand having relative density 10%, 20%, 43.8% at confining pressure 100 Kpa
Fig. 18 Model simulations with undrained Triaxial test results on Bagmati Sand having relative density 10%, 20%, 57.1% at confining pressure 200 Kpa
400
Deviatoric Stress, q ( Kpa )
350
300
250
200
150
100
50
0
Lab Simulation TSP
0 50 100 150 200 250 300
Mean Effective Pressure, p' ( Kpa )
1000
Deviatoric Stess, q ( Kpa )
900
800
700
600
500
400
300
200
100
0
Lab Simulation TSP
0 100 200 300 400 500 600 700
Effective Mean Pressure, p' ( Kpa )
Fig. 15 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.754 and state index 5.16 at confining pressure 200 Kpa
Fig. 19 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.802 and state index 1.54 at confining pressure 400 Kpa
450
400
Deviatoric Stress, q ( Kpa )
350
300
250
200
150
100
50
0
Lab Simulation
0 2 4 6 8 10 12 14
Strain ( % )
1000
900
Deviatoric Stress, q ( Kpa )
800
700
600
500
400
300
200
100
0
lab Simulation
0 5 10 15 20 25
Strain ( % )
Fig. 16 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.754 and state index 5.16 at confining pressure 200 Kpa
Fi. 20 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.802 and state index 1.54 at confining pressure 400 Kpa
80
Lab Simulation
70
Pore Pressure, u ( Kpa )
60
50
40
30
20
180
160
Pore Pressure, u ( Kpa )
140
120
100
80
60
40
Lab Simulation
10
0
0 2 4 6 8 10 12 14
strain ( % )
20
0
0 5 10 15 20 25
Strain ( % )
Fig. 17 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.754 and state index 5.16 at confining pressure 200 Kpa
Fig. 21 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.802 and state index 1.54 at confining pressure 400 Kpa
Deviatoric Stress, q ( Kpa )
Deviatoric Stress, q ( Kpa )
1000
900
42.7%
20%
10%
TSP
800
350
700
300
57.1%
10%
600
500
250
20% TSP
400
300
200
100
0
0 100 200 300 400 500 600 700
Mean Effective Pressure, p' ( Kpa )
200
150
100
50
0
0 50 100 150 200 250 300
Mean Effective Pressure, p' ( Kpa )
Fig. 22 Model simulations with undrained Triaxial test results on Bagmati
Sand having relative density 10%, 20%, 42.7% at confining pressure 400 Kpa
1000
Fig. 26 Model simulations with undrained Triaxial test results on Bagmati Sand having relative density 10%, 20%, 57.1% at confining pressure 200 Kpa
900
Deviatoric Stress, q ( Kpa )
800
700
600
500
400
300
200
100
0
L100
S100
L200
S200
L400
S400
Kf – Line S
Kf – Line – L
0 100 200 300 400 500 600 700
Effective Mean Pressure, p' ( Kpa )
3.2 Performance of Model
For the verification and validation of the model, the predictions of the model are compared to the laboratory undrained monotonic loading of Bagmati sand, and also for Toyoura sand. For Bagmati sand, the experimental data is taken from Jha (2004), and for Toyoura sand, the experimental data is taken Mitsutoshi Yoshimine ( 1996 1999 ). The data from the these tests comprise various aspects of the undrained behaviour of the sand and particularly emphasize the effects of density on stressstrain behaviour and effective stress path along the pore pressure development.
Fig. 23 Comparison of model simulations with undrained Triaxial test results on Bagmati Sand having void ratio of 0.779, 0.754, 0.802 and state index 4.96, 5.16, 1.54 at confining pressure 100, 200, 400 Kpa respectively
350
Toyoura Sand – RD. 10% – p' 100 Kpa
3.2.1 Simulation for Bagmati Sand
The triaxial compressional behaviour of the Bagmati sand has been predicted for the three isotropic consolidation pressures, 100 Kpa, 200 Kpa, 400 Kpa, and the relative densities 42.7 %, 57.1 %, 43.8 %, respectively. The simulation is also done for relative densities 10 %, 20 %. The simulation is done in strain controlled manner at a strain
Deviatoric Stress, q ( Kpa )
300
250
200
150
100
50
0
Bagmati Sand – RD. 10% – p' 100 Kpa
Toyoura Sand – RD. 20% – p' 100 Kpa " Bagmati Sand – RD. 20% – p' 100 Kpa "
0 50 100 150 200 250
Mean effective Pressure, p' ( Kpa )
rater of 0.01 mm.
Fig. 11 and fig. 12 shows the comparison between predicted and measured ( i.e. laboratory ) behaviour of isotropic consolidated the sand with relative density 42.7 % and consolidation pressure of 100 Kpa. Apparently, good agreement between the predicted and measured behaviour is obtained for both effective stress path ( ESP ) and stressstrain curve. Initially contractive response associated with development of pore pressure and movement of the ESP toward right, approximately uniform
Fig. 24 Model simulations with undrained Triaxial test results on Toyoura Sand having relative densities 22.9 % and 65.5 % at confining pressure 100 respectively.
350
Deviatoric Stress, q ( Kpa )
300
250
200
Toyoura Sand – RD. 10% – p' 100 Kpa Bagmati Sand – RD. 10% – p' 100 Kpa Toyoura Sand – RD. 20% – p' 100 Kpa " Bagmati Sand – RD. 20% – p' 100 Kpa "
150
100
50
0
0 50 100 150 200 250
Mean effective Pressure, p' ( Kpa )
Fig. 25 Comparison of simulated results of Bagmati Sand and Toyoura Sand for relative density 10%, and 20% at confining pressure 100 Kpa respectively
All the material parameters are summarized in table 1.
slope, until phase transformation state is reached , and then dilative response with further movement of the ESP toward right with lower uniform slope than the previous one. Up to 6
% of strain, there is development of positive pore pressure and only after 6 %, there is development of negative pore pressure, even though, the Bagmati sand dilates before 6 % of strain. (fig. 12, and 13). This is due to the fact that negative pore pressure development due to shearing , is lesser than , due to increase in effective mean pressure ( i.e. mechanical phenomenon ). For the relative densities 10 % and 20 %, the ESPs have been predicted, and even at 10 % of the relative density, the Bagmati sand does not exhibit complete contractive behaviour at the confining pressure, 100 Kpa (fig. 14). At 10 %, initially contractive response associated with
development of pore pressure and movement of the ESP toward left, until phase transformation state is reached , and then dilative response with movement of the ESP toward right.
Similarly, fig. 15 and fig. 16 shows the comparison between predicted and measured behaviour of isotropic consolidated sand with relative density 57.1 % and consolidation pressure of 200 Kpa. It may seen that good agreement between the predicted and measured behaviour is obtained for both effective stress path and stressstrain curve. Initially contractive response associated with development of pore pressure and movement of the ESP toward right concaving to the left, until phase transformation state is reached , and then dilative response with further movement of the ESP toward
right with approximately uniform slope. The pore pressure development has been shown in fig. 17. Again, for the relative densities 10 % and 20 %, the ESPs have been predicted, and even at 10 % of the relative density, the Bagmati sand does not exhibit complete contractive behaviour at the confining pressure, 200 Kpa ( fig. 18). At 10 % and 20 % , initially contractive response associated with development of pore pressure and movement of the ESP toward left, until phase transformation state is reached , and then dilative response with movement of the ESP toward right.
Figure 19 and fig. 20 also show the comparison between predicted and measured behaviour of isotropic consolidated sand with relative density 42.7 % and consolidation pressure of 400 Kpa. Good agreement between the predicted and measured behaviour is obtained for both effective stress path and stressstrain curve. Initially contractive response associated with development of pore pressure and movement of the ESP toward right concaving to the left, until phase transformation state is reached , and then dilative response with further movement of the ESP toward right with approximately uniform slope. The pore pressure development has been shown in fig. 21. Again, for the relative densities 10
% and 20 %, the ESPs have been predicted, and even at 10 % of the relative density, the Bagmati sand does not exhibit complete contractive behaviour at the confining pressure, 400 Kpa ( fig. 22). At 10 % and 20 % , initially contractive response associated with development of pore pressure and movement of the ESP toward left, until phase transformation state is reached , and then dilative response with movement of the ESP toward right
The failure envelops, represented by Kf line, are drawn for both the predicted as well as measured case, shown in fig. 23. These simulation and laboratory failure envelops gives the M,
i.e. qp stress at steady state, with corresponding values 1.42,
concaving to the right, until phase transformation state is reached , and then dilative response with further movement of the ESP toward right with approximately uniform slope. At relative density 65.5 %, the pore pressure development due to the shearing phenomenon, is much lesser than that of due to the mechnical phenomenon.
3.2.3 Comparison between Bagmati Sand and Toyoura Sand The behaviour of the Bagmati sand and Toyoura sand has been compared, and it is remarked that at 10 % of relative density, the toyoura sand behaves fully contractive, where as the Bagmati sand shows contractive behaviour initially and then dilates ( fig.25 ) at 100 Kpa. But, for relative density 20 %, both behave similar behaviour with different magnitudes. This pictorial presentation reflects that for same relative density and confining pressure, Toyoura sand is more susceptible to liquefaction and flow deformation.
The model correctly predicts the stressstrain characteristics observed in the undrained triaxial compression tests on medium dense, loose, very loose sand, providing gradual change in these characteristics from those typical for dense sand, with continuous increase of the shear stress with straining, to those typical for loose sand, with drop in the shear stress associated with development at of the large strains including the extreme cases with zero residual strength.


CONCLUSION AND RECOMMENDATION
The research is an attempt to represent the triaxial undrained monotonic behaviour on a physically relevant basis and provide a constitutive model comprising balanced features of accuracy, simplicity, and versatility. The theoretical part of the study is focused on the development of a proper modeling concept in order to account for the combined effects of density and mean effective stress i.e. state index as state parameter, on the behaviour of cohesionless soils. The influence of various
and 1.50. And,
M 6 sin ss
3 sin ss
Finally, the angle of inter
factors on the monotonic behaviour of sandy soils is accounted for in an elasticplastic constitutive model as a close result of
particle friction during steady state is calculated to be 35.03,
and 36.86 respectively. Thus, the model gives little overestimated value of the interparticle friction during steady state.
3.2.2 Simulation for Toyoura Sand
Also, the triaxial compressional behaviour of the Toyoura sand has been predicted for the three isotropic consolidation pressures, 100 Kpa and the relative densities 22.9 %, 65.5 % respectively ( fig. 24). The simulation is done at strain rate of
0.01 mm.
Fig. 24 shows the comparison between predicted and measured behaviour. The model also set good agreement between the predicted and measured behaviour for the sand, for both effective stress path ( ESP ) and stressstrain curve. At
22.9 %, initially contractive response with movement of the ESP toward left concaving to the leftdown, until phase transformation state is reached , and then dilative response with further movement of the ESP toward right with approximately uniform slope. Similarly, at 65.5 %, initially contractive response with movement of the ESP toward
this work. The framework is developed by the modification of the torsional simple shear test framework, in stress space definition, yielding criterion, stressdilatancy relationship, plastic modulus, and loading index.

Characteristic Curves of Bagmati Sand
The three characteristic curves of Bagmati sand; the loosest isotropic consolidation line, the steady state lines and the upper reference line, are achieved for the mean effective pressure up to 500 Kpa. With these curves reference, the flow and nonflow case can be assessed for the given initial state ( e, p ).

Undrained Monotonic behaviour of Bagmati Sand
For these three set of relative densities and isotropic consolidation pressure, the pore pressures development due to the dilation, are not so dominant over the mechanical phenomenon that it generates net negative pore pressure, except for 42.7 % relative density and 100 Kpa after 6 % of strain.
Though, its very much difficult to prepare a sample of 10 % relative density and below, the behaviour of the sand is predicted, and even at this relative density and confining pressure up to 400 Kpa, the Bagmati sand does not exhibit complete contractive behaviour. It also conclude that the sand is less susceptable to liquefaction and flow deformation.
The failure envelops, Kf line, of Bagmati sand has slope 1.42:1, from the lab result, and the angle of interparticle friction during steady state is and 35.03 .

Performance and Application of the Model
The model has ability to simulate undrained response from that typical of dense to that typical of loose sand including contractive flow or limited flow, for wide range of effective stresses. To check the effectiveness of the proposed model, simulated results are compared to the measured behaviour of sand in triaxial undrained monotonic tests. For the purpose, strain controlled version of the model are developed. Besides, the physical relevance of the model, it is important to emphasize the practical benefits arising from the fact that a single set of parameters characterizes the behaviour over a wide range of densities and stresses. Moreover, it is proved that the model is computationally efficient and has negligible computer memory requirements. The model can be applied for the various geotechnical studies like liquefaction, shear strength evaluation, residual strength evaluation, etc.
Key features of the proposed constitutive model can be summarized as follows :

The characterization of sand behaviour is based on the state concept via state index.

An energy based stressdilatancy relation is employed.

The stressstraindilatancy model is incorporated into the elastoplastic framework of sand modeling in order to establish a mathematical formulation.

The plasticity formulation includes the assumptions for continuous yielding and the dependence of the plastic strain increment direction on the stress increment direction.

The integration of the state concept and the elasto plasticity into a single framework leads to the unique capability of the model to capture the prominent features of sand behaviour over a wide range of density and stress states.

Since, the research is limited for the model simulation of triaxial undrained monotonic behaviour, it is recommended to extend the model framework for the triaxial undrained cyclic behaviour of sand.


REFERENCES