• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.09a
• /
• pp.1090-1094
• /
• 2009

Many problems in geotechnical engineering such as slop failure, debris flow, ground heaving due to embankment, and lateral flow caused by liquefaction are related to large deformation rather than small deformation. Traditional numerical methods such as finite element and finite difference methods have a difficulty to solve such large deformations because they use grids. A particle method was developed for fluid dynamics. The particle method can solve large deformation problems because it uses particles to discretize differential equations. It can also include soil constitutive model and thus solve soil behavior on various boundary conditions. In this study, a particle method, which is based on particles rather than grids, is introduced and used to simulate large deformation including soil failure. The developed method can be applied for various large deformation problems in geotechnical engineering because it incorporates soil constitutive models.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.5
• /
• pp.85-94
• /
• 2012

Recently, growing interests in frozen ground have stimulated us to advance fundamental theories and systematic researches on soil behavior under freezing conditions. Unlike the well-established soil mechanics theory, temperature variation and phase change of pore-water cause water migration to cold side, ground heaving, sharp increase in earth pressure, etc., which bring about serious problems in frozen geotechnical structures. Elasto-plastic mechanical constitutive model for frozen/unfrozen soil subjected to fully coupled THM phenomena is formulated based on a new stress variable that is continuous in frozen-unfrozen transitional regions. Numerical simulations are conducted to discuss numerical reliability and applicability of the developed constitutive model: one-dimensional heaving pressure, tri-axial compression test, and one-side freezing tests. The numerical results show that developed model can efficiently describe complex THM phenomena of frozen soil, and they can be utilized to analyze and design the geotechnical structures under freezing conditions, and predict their long-term behavior.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.10a
• /
• pp.1199-1203
• /
• 2008

It has been recognized that unsaturated soil behavior plays an importantrole in geomechanics. In the last decade several constitutive models have been proposed and used in the analysis. Many of them, however, are constructed in the frame work of rate independent model such as elasto-plastic one. Although rate dependency is an important characteristics of soil for both saturated and unsaturated soils, very few models have been developed taking account of rate dependency. In the present paper, we have developed an elasto-viscoplastic model considering an effect of suction based on the overstress-type viscoplasticity with soil structure degradation. In the model, we have adopted an averaged pore pressure composed of pore water pressure and air pressure to determine the effective stress.

• Geomechanics and Engineering
• /
• v.24 no.2
• /
• pp.193-203
• /
• 2021

The time-history finite element analysis is usually used to evaluate the seismic response of shallow foundations. However, the literature lacks studies on the influence of the soil constitutive model complexity on the seismic response of shallow foundations. This study, thus, aims to fill this gap by investigating the seismic response of shallow foundation resting on dry silica sand using the linear elastic (LE) model, elastic-perfectly-plastic (EPP) model, and hardening soil with small strain stiffness (HS small) model. These models have been used because it is intended to compare the results of a soil constitutive model that accurately captures the seismic response of the soil-structure interaction problems (which is the HS small model) with simpler models (the LE and EPP models) that are routinely used by practitioners in geotechnical designs. The results showed that the LE model produces a very small seismic settlement value which is approximately equal to zero. The EPP model predicts a seismic settlement higher than that produced using the HS small model for earthquakes with a peak ground acceleration (PGA) lower than 0.25 g for a relative density of 45% and 0.40 g for a relative density of 70%. However, the HS small model predicts a seismic settlement higher than the EPP model beyond the aforementioned PGA values with the difference between both models increases as the PGA rises. The results also showed that the LE and EPP models predict similar trend and magnitude of the acceleration-time relationship directly below the foundation, which was different than that predicted using the HS small model. The results reported in this paper provide a useful benchmark for future numerical studies on the response of shallow foundations subjected to seismic shake.

• Journal of the Korean Geosynthetics Society
• /
• v.10 no.4
• /
• pp.29-36
• /
• 2011

In this study, a series of the isotropic compression-expansion tests and the drained triaxial tests were performed on Pocheon granite soil with various the dry densities of $16.67kN/m^3$, $17.26kN/m^3$ and $17.65kN/m^3$. Using the tests results the characteristic of the parameters of Lade's single hardening constitutive model were investigated. The soil parameters such as kur and n related to elastic behavior, m and ${\eta}_1$ related to failure criterion, c and p related to hardening function and ${\psi}_2$ and ${\mu}$ related to plastic potential show in a positive linear relationship with the dry density. Since the soil parameters h and representing yield function do not change much to relative density and also are closely related to failure criterion, they can be replaced by failure criterion. We also observed that predicted values from the Lade's single hardening constitutive model were well consistent with the observed data.

• Journal of Navigation and Port Research
• /
• v.36 no.5
• /
• pp.395-400
• /
• 2012

The stress-strain relationship of soil is dependent on a number of factors such as soil type, density, stress level and stress path. Th accurate stress-stain relationship can be predict using a constitutive model incorporated all influencing factors. In this study, an isotropic compression-expansion test and a series of drained conventional triaxial tests with several stress paths were performed on Baekma river sand to investigate parameters characteristics of Lade's single work hardening model depending on the stress path.. Based on test results, the parameters of yield function (h, ${\alpha}$) are not much influenced by stress level and stress path, the these parameters do affect a little bit of stress-strain behavior. The parameters h and ${\alpha}$ are closely related to failure criterion ${\eta}_1$, they can be replaced by failure criterion parament. We also observed that predicted values from the Lade's single hardening constitutive model are well matched with the observed data.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2001.11a
• /
• pp.471-476
• /
• 2001

Even though structural performance evaluation techniques for reinforced concrete structures have been improved, there are still many problems in the evaluation of structural performance for underground structures which interacts with surrounding soils. Since experimental evaluation of underground RC structures considering the interaction with the surrounding soil medium is quite difficult to be simulated, the evaluation for underground RC structures using an analytical method can be applied very usefully, For underground structures interacted with surrounding soils, it is important to consider path-dependent RC constitutive model, soil constitutive model, and interface model between structure and soil, simultaneously. In this paper, an elastoplastic interface model which consider thickness was proposed and importance of interface model is discussed. The effects of stiffness of structures to entire underground RC system are investigated through numerical experiment for underground RC structure for different reinforcement ratios and thickness of interfaces.

• Geomechanics and Engineering
• /
• v.36 no.4
• /
• pp.355-366
• /
• 2024

The original elastoplastic Hardening Soil model is formulated actually partly under hexagonal pyramidal Mohr-Coulomb failure criterion, and can be only used in specific stress paths. It must be completely generalized under Mohr-Coulomb criterion before its usage in engineering practice. A set of generalized constitutive equations under this criterion, including shear and volumetric yield surfaces and hardening laws, is proposed for Hardening Soil model in principal stress space. On the other hand, a Mohr-Coulumb type yield surface in principal stress space comprises six corners and an apex that make singularity for the normal integration approach of constitutive equations. With respect to the isotropic nature of the material, a technique for processing these singularities by means of Koiter's rule, along with a transforming approach between both stress spaces for both stress tensor and consistent stiffness matrix based on spectral decomposition method, is introduced to provide such an approach for developing generalized Hardening Soil model in finite element analysis code ABAQUS. The implemented model is verified in comparison with the results after the original simulations of oedometer and triaxial tests by means of this model, for volumetric and shear hardenings respectively. Results from the simulation of oedometer test show similar shape of primary loading curve to the original one, while maximum vertical strain is a little overestimated for about 0.5% probably due to the selection of relationships for cap parameters. In simulation of triaxial test, the stress-strain and dilation curves are both in very good agreement with the original curves as well as test data.

• Journal of the Korean Geosynthetics Society
• /
• v.12 no.4
• /
• pp.87-98
• /
• 2013

This study was performed a series of the isotropic compression-expansion tests and the drained triaxial tests with various the relative densities 40%, 60%, 80% and 95% for bottom ash. Using the tests results the characteristic of the parameters of Lade's single hardening constitutive model were investigated. The soil parameters Kur and n representing elastic behavior are not much affected by the change of the relative density. The other parameters such as failure criterion(m, ${\eta}_1$), hardening function(c, p) and plastic potential(${\psi}_2$, ${\mu}$) are in a positive linear relationship with the relative density. Since the soil parameters h and ${\alpha}$ representing yield function do not change much to the change of relative density and also closely related to failure criterion, they can be replaced by failure criterion ${\eta}_1$. We also observed that predicted values from the Lade's single hardening constitutive model were well consistent with the observed data.