• Structural Engineering and Mechanics
• /
• v.91 no.1
• /
• pp.103-121
• /
• 2024

The smoothed particle hydrodynamics (SPH) method is a numerical technique used in dynamic analysis to simulate the fluid-like behavior of materials under extreme conditions, such as those encountered in explosions or high velocity impacts. In SPH, fluid or solid materials are discretized into particles. These particles interact with each other based on certain smoothing kernels, allowing the simulation of fluid flows and predict the response of solid materials to shock waves, like deformation, cracking or failure. One of the main advantages of SPH is its ability to simulate these phenomena without a fixed grid, making it particularly suitable for analyzing complex geometries. In this study, the structural damage to a masonry arch bridge subjected to blast loading was investigated. A high-fidelity micro-model was created and the explosives were modeled using the SPH approach. The Johnson-Holmquist II damage model and the Mohr-Coulomb material model were considered to evaluate the masonry and backfill properties. Consistent with the principles of the JH-II model, the authors developed a VUMAT code. The explosive charges (50 kg, 168 kg, 425 kg and 1000 kg) were placed in close proximity to the deck and pier of a bridge. The results showed that the 50 kg charges, which could have been placed near the pier by a terrorist, had only a limited effect on the piers. Instead, this charge caused a vertical displacement of the deck due to the confinement effect. Conversely, a 1000 kg TNT charge placed 100 cm above the deck caused significant damage to the bridge.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.3
• /
• pp.53-64
• /
• 2003

Many transmission towers, high-rise buildings and bridges are constructed near steep slopes and are supported by large-diameter piles. These structures may be subjected to large lateral loads, such as violent winds and earthquakes. Widely used types of foundations for these structures are pier foundations, which have large-diameters with high stiffness. The behavior of a pier foundation subjected to lateral loads is similar to that of a short rigid pile because both elements seem to fail by rotation developing passive resistance on opposite faces above and below the rotation point, unlike the behavior of a long flexible pile. This paper describes the results of several numerical studies performed with a three-dimensional finite element method (FEM) of model tests of a laterally loaded short pile located near slopes, respectively. In this paper, the results of model tests of single piles and pile groups subjected to lateral loading, in homogeneous sand with 30$^{\circ}$ slopes and horizontal ground were analyzed by the 3-D FE analyses. The pile was assumed to be linearly elastic. The sand was assumed to have non-associative characteristics, following the MC-DP model. The failure criterion is governed by the Mohr-Coulomb equation and the plastic potential is given by the Drucker-Prager equation. The main purpose of this paper is the validation of the 3-D elasto-plastic FEM by comparisons with the experimental data.

• Geomechanics and Engineering
• /
• v.23 no.3
• /
• pp.245-259
• /
• 2020

Tunnels have become an indispensable part of metro cities. Blast resistance design of tunnel has attracted the attention of researchers due to numerous implosion event. Present paper deals with the non-linear finite element analysis of rock tunnel having shear zone subjected to internal blast loading. Abaqus Explicit schemes in finite element has been used for the simulation of internal blast event. Structural discontinuity i.e., shear zone has been assumed passing the tunnel cross-section in the vertical direction and consist of Highly Weathered Granite medium surrounding the tunnel. Mohr-Coulomb constitutive material model has been considered for modelling the Highly Weathered Granite and the shear zone material. Concrete Damage Plasticity (CDP), Johnson-Cook (J-C), Jones-Wilkins-Lee (JWL) equation of state models are used for concrete, steel reinforcement and Trinitrotoluene (TNT) simulation respectively. The Coupled-Eulerian-Lagrangian (CEL) method of modelling for TNT explosive and air inside the tunnel has been adopted in this study. The CEL method incorporates the large deformations for which the traditional finite element analysis cannot be used. Shear zone orientations of 0°, 15°, 30°, 45°, 60°, 75° and 90°, with respect to the tunnel axis are considered to see their effect. It has been concluded that 60° orientation of shear zone presents the most critical situation.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.10
• /
• pp.17-28
• /
• 2018

Important parameters for the stability checks of cantilever wall are the active earth pressure and the weight of soil above the heel of the base slab. If the heel length is so long enough that the shear zone bounded by the failure plane is not obstructed by the stem of the wall, the Rankine active condition is assumed to exist along the vertical plane which is located at the edge of the heel of the base slab. Then the Rankine active earth pressure equations may be theoretically used to calculate the lateral pressure on the vertical plane. However, in case of the cantilever wall with a short heel, the application of Rankine theory is not only theoretically incorrect but also makes the lateral earth pressure larger than the actual pressure and results in uneconomical design. In this study, for the cantilever wall with a short heel the limit analysis method is used to investigate the mechanism of development of the active earth pressure and then the magnitude and location of the resultants of the pressure and the weight of the soil above the heel are determined. The calculated results are compared with the existing methods for the stability check. In case of the cantilever wall with a short heel, the results by the Mohr circle method and Teng's method show max. 3.7% and 32% larger than those of the limit analysis method respectively.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.11
• /
• pp.33-40
• /
• 2015

Volcanic rocks in Jeju Island indicate the differences in geological and mechanical characteristics from region to region, and have vesicular structure caused by various environmental factors. In this study, triaxial compressive strength tests were conducted for intact rocks sampled in northeastern onshore and offshore, southeastern offshore and northwestern offshore of Jeju Island. The estimated cohesion and internal friction angle from the results of triaxial compression tests were compared and analyzed with absorption, a parameter representing the vesicular properties of basalts in Jeju Island. As a result, it was found that the relationship between cohesion and absorption could be classified clearly, considering two different linear relationships in bulk specific gravity and absorption. As the absorption increases, the cohesion decreases exponentially. In addition, the internal friction angle decreases almost linearly with increasing in the absorption, regardless of the relationships in bulk specific gravity and absorption.

• Geomechanics and Engineering
• /
• v.11 no.3
• /
• pp.439-455
• /
• 2016

This paper presented solutions of displacement and stress for a circular opening which is reinforced with grouted rock bolt. It satisfies the Mohr-Coulomb (M-C) or generalized Hoek-Brown (H-B) failure criterion, and exhibits elastic-brittle-plastic or strain-softening behavior. The numerical stepwise produce for strain-softening rock mass reinforced with grouted rock bolt was developed with non-associative flow rules and two segments piecewise linear functions related to a principle strain-dependent plastic parameter, to model the transition from peak to residual strength. Three models of the interaction mechanism between grouted rock bolt and surrounding rock proposed by Fahimifar and Soroush (2005) were adopted. Based on the axial symmetrical plane strain assumption, the theoretical solution of the displacement and stress were proposed for a circular tunnel excavated in elastic-brittle-plastic and strain-softening rock mass compatible with M-C or generalized H-B failure criterion, which is reinforced with grouted rock bolt. It showed that Fahimifar and Soroush's (2005) solution is a special case of the proposed solution for n = 0.5. Further, the proposed method is validated through example comparison calculated by MATLAB programming. Meanwhile, some particular examples for M-C or generalized H-B failure criterion have been conducted, and parametric studies were carried out to highlight the influence of different parameters (e.g., the very good, average and very poor rock mass). The results showed that, stress field in plastic region of surrounding rock with considering the supporting effectiveness of the grouted rock bolt is more than that without considering the effectiveness of the grouted rock bolt, and the convergence and plastic radius are reduced.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.381-395
• /
• 2012

Dynamic response of Pile Supported Structures is highly depended on Soil Pile Structure Interaction. In this paper, by comparison of experimental and numerical dynamic responses of a prototype jacket offshore platform for both hinge based and pile supported boundary conditions, effect of soil-pile-structure interaction on dynamic characteristics of this platform is studied. Jacket and deck of a prototype platform is installed on a hinge-based case first and then platform is installed on eight skirt piles embedded on continuum monolayer sand. Dynamic characteristics of platform in term of natural frequencies, mode shapes and modal damping are compared for both cases. Effects of adding and removing vertical bracing members in top bay of jacket on dynamic characteristics of platform for both boundary conditions are also studied. Numerical simulation of responses for the studied platform is also performed for both mentioned cases using capability of ABAQUS and SACS software. The 3D model using ABAQUS software is created using solid elements for soil and beam elements for jacket, deck and pile members. Mohr-Coulomb failure criterion and pile-soil interface element are used for considering nonlinear pile soil structure interaction. Simplified modeling of soil-pile-structure interaction effect is also studied using SACS software. It is observed that dynamic characteristics of the system changes significantly due to soil-pile-structure interaction. Meanwhile, both of complex and simplified (ABAQUS and SACS, respectively) models can predict this effect accurately for such platforms subjected to dynamic loading in small range of deformation.

• Geomechanics and Engineering
• /
• v.3 no.4
• /
• pp.291-321
• /
• 2011

The paper describes a simple numerical FLAC model that was developed to simulate the dynamic response of two instrumented reduced-scale model reinforced soil walls constructed on a 1-g shaking table. The models were 1 m high by 1.4 m wide by 2.4 m long and were constructed with a uniform size sand backfill, a polymeric geogrid reinforcement material with appropriately scaled stiffness, and a structural full-height rigid panel facing. The wall toe was constructed to simulate a perfectly hinged toe (i.e. toe allowed to rotate only) in one model and an idealized sliding toe (i.e. toe allowed to rotate and slide horizontally) in the other. Physical and numerical models were subjected to the same stepped amplitude sinusoidal base acceleration record. The material properties of the component materials (e.g. backfill and reinforcement) were determined from independent laboratory testing (reinforcement) and by back-fitting results of a numerical FLAC model for direct shear box testing to the corresponding physical test results. A simple elastic-plastic model with Mohr-Coulomb failure criterion for the sand was judged to give satisfactory agreement with measured wall results. The numerical results are also compared to closed-form solutions for reinforcement loads. In most cases predicted and closed-form solutions fall within the accuracy of measured loads based on ${\pm}1$ standard deviation applied to physical measurements. The paper summarizes important lessons learned and implications to the seismic design and performance of geosynthetic reinforced soil walls.

• Geomechanics and Engineering
• /
• v.7 no.3
• /
• pp.263-277
• /
• 2014

In the conventional design of retaining structures in a seismic zone, seismic inertia forces are commonly assumed to act upwards and towards the wall facing to cause a maximum active thrust or act upwards and towards the backfill to cause a minimum passive resistance. However, under certain circumstances this design approach might underestimate the dynamic active thrust or overestimate the dynamic passive resistance acting on a rigid retaining structure. In this study, a new analytical method for dynamic active and passive forces in c-${\phi}$ soils with an infinite slope was proposed based on the Rankine earth pressure theory and the Mohr-Coulomb yield criterion, to investigate the influence of seismic inertia force directions on the total active and passive forces. Four combinations of seismic acceleration with both vertical (upwards or downwards) and horizontal (towards the wall or backfill) directions, were considered. A series of dimensionless dynamic active and passive force charts were developed to evaluate the key influence factors, such as backfill inclination ${\beta}$, dimensionless cohesion $c/{\gamma}H$, friction angle ${\phi}$, horizontal and vertical seismic coefficients, $k _h$ and $k_v$. A comparative study shows that a combination of downward and towards-the-wall seismic inertia forces causes a maximum active thrust while a combination of upward and towards-the-wall seismic inertia forces causes a minimum passive resistance. This finding is recommended for use in the design of retaining structures in a seismic zone.

• Earthquakes and Structures
• /
• v.18 no.2
• /
• pp.201-213
• /
• 2020

Seismic damages that occurred by the effects of epicenter distance of the earthquake are one of the most important problems for the earthquake engineering. In this study, it is aimed to examine the nonlinear seismic behaviors of concrete gravity (CG) dams considering various epicenter distances. For this purpose, Boyabat CG dam that is one of the biggest concrete gravity dams in Turkey is selected as a numerical application. FLAC3D software based on finite difference method is used for modelling and analyzing of the dam. Drucker-Prager nonlinear material model is used for the concrete body and Mohr-Coulomb nonlinear material model is taken into account for the foundation. Special interface elements are used between dam body and foundation to represent interaction condition. Free-field and quiet non-reflecting boundary conditions are utilized for the main surfaces of 3D model. Total 5 various epicenter distances of 1989 Loma Prieta earthquake are considered in 3D earthquake analyses and these distances are 5 km, 11 km, 24 km, 85 km and 93 km, respectively. According to 3D seismic results, x-y-z displacements, principal stresses and shear strain failures of the dam are evaluated in detail. It is clearly seen from this study that the nonlinear seismic behaviors of the CG dams change depending to epicenter distance of the earthquake. Thus, it is clearly recommended in this study that when a CG dam is modelled or analyzed, distance of the earthquake fault to the dam should be strongly examined in detail. Otherwise, earthquake damages can be occurred in the concrete dam body by the effects of seismic loads.