• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.4
• /
• pp.267-274
• /
• 2024

Hydrogen embrittlement fracture poses a challenge in ensuring the structural integrity of materials exposed to hydrogen-rich environments. This study advances our comprehension of hydrogen-induced fracture through an integrated numerical modeling approach. In addition, it employs a ductile fracture model named the Gurson-cohesive model (GCM) and hydrogen diffusion analysis. GCM is employed as a fracture model that combines the Gurson model to illustrate the continuum damage evolution and the cohesive zone model to describe crack surface discontinuity and softening behavior. Moreover, porosity and stress triaxiality are considered as crack initiation criteria . A hydrogen diffusion analysis is also integrated with the GCM to account for hydrogen enhanced decohesion (HEDE) mechanisms and their subsequent impacts on crack initiation and propagation. This framework considers the influence of hydrogen on the softening behavior of the traction-separation relationship on the discontinuous crack surface. Parametric studies explore the sensitivity to diffusion properties and hydrogen-induced fracture properties. By combining numerical models of hydrogen diffusion and the ductile fracture model, this study provides an understanding of hydrogen-induced fracture and thereby contributes significantly to the ongoing efforts to design materials that are resilient to hydrogen embrittlement in practical engineering applications.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.4
• /
• pp.5-14
• /
• 2016

Numerical simulation of dynamic soil-pile-structure interaction embedded in a dry sand was carried out. 3D model of the dynamic centrifuge model tests was formulated in a time domain to consider nonlinear behavior of soil using the finite difference method program, FLAC3D. As a modeling methodology, Mohr-Coulomb criteria was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling (Kim et al., 2012) was used as boundary condition to reduce analysis time. Calibration process for numerical modeling results and test results was performed through the parametric study. Verification process was then performed by comparing numerical modeling results with another test results. Based on the calibration and validation procedure, it is identified that proposed modeling method can properly simulate dynamic behavior of soil-pile system in dry condition.

• Journal of the Korean Geotechnical Society
• /
• v.20 no.4
• /
• pp.89-102
• /
• 2004

The generated blasting pressure wave initiated under decoupled-charge condition is a function of peak blasting pressure, rise time, and wave-shape function. The peak blasting pressure and the rise time are also the function of explosive and rock properties. The probabilistic distributions of explosive and rock properties are derived from the results of their property tests. Since the probabilistic distributions of explosive and rock properties displayed a normal distribution, the peak blasting pressure and the rise time can also be regarded as a normal distribution. Parameter analysis and uncertainty analysis were performed to identify the most influential parameter that affects the peak blasting pressure and the rise time. Even though the explosive properties were found to be the most influential parameters on the peak blasting pressure and the rise time from the parameter analyses, the result of uncertainty analysis showed that rock properties constituted major uncertainties in estimating the peak blasting pressure and the rise time rather than explosive properties. Damage and overbreak of the remaining rock around the excavation line induced by blasting were evaluated by dynamic numerical analysis. A user-subroutine to estimate the rock damage was coded based on the continuum damage mechanics. This subroutine was linked to a commercial program called 'ABAQUS/Explicit'. The results of dynamic numerical analysis showed that the rock damages generated by the initiation of stopping hole were larger than those from the initiation of contour hole. Several methods to minimize those damages were proposed such as relocation of stopping hole, detailed subdivision of rock classification, and so on. It was found that fracture probability criteria and fractured zones could be distinctively identified by applying fuzzy-random probability.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.2
• /
• pp.303-312
• /
• 2008

In this paper, the seismic behavior of shear wal-frame systems is analyzed. The governing equations of the wall-frame systems with outrigger truss are formulated through the continuum approach and the whole structure is idealized as a shear-flexural cantileverwith rotational spring. The effect of shear deformation and flexural deformation of the wall-frame and outrigger trusses are considered and incorporated in the formulation of the wall-frame structures with and without outriggers are compared by using finite element analysis incorporated with the Newmark-${\beta}$ method. Numerical results are obtained and compared with the finite element package MIDAS. The proposed method is found to be simple and efficient, and provides reason ably accurate results in the early design stage of tall building structures.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.13-20
• /
• 2024

In this paper, a parallel analysis algorithm for Smoothed Particle Hydrodynamics (SPH), one of the numerical methods for fluidic materials, is introduced. SPH, which is a meshless method, can represent the behavior of a continuum using a particle-based approach, but it demands substantial computational resources. Therefore, parallel analysis algorithms are essential for SPH simulations. The domain decomposition algorithm, which divides the computational domain into partitions to be independently analyzed, is the most representative method among parallel analysis algorithms. In Discrete Element Method (DEM) and Molecular Dynamics (MD), the Cartesian coordinate-based domain decomposition method is popularly used because it offers advantages in quickly and conveniently accessing particle positions. However, in SPH, it is important to share particle information among partitioned domains because SPH particles are defined based on information from nearby particles within the smoothing length. Additionally, maintaining CPU load balance is crucial. In this study, a highly parallel efficient algorithm is proposed to dynamically minimize the size of orthogonal domain partitions to prevent excess CPU utilization. The efficiency of the proposed method was validated through numerical analysis models. The parallel efficiency of the proposed method is evaluated for up to 30 CPUs for fluidic models, achieving 90% parallel efficiency for up to 28 physical cores.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.9
• /
• pp.51-62
• /
• 2016

Parametric studies for various site conditions by using 3d numerical model were carried out in order to estimate dynamic behavior of soil-pile-structure system in dry soil deposits. Proposed model was analyzed in time domain using FLAC3D which is commercial finite difference code to properly simulate nonlinear response of soil under strong earthquake. Mohr-Coulomb criterion was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling was used as boundary condition to reduce analysis time. Also, initial shear modulus and yield depth were appropriately determined for accurate simulation of system's nonlinear behavior. Parametric study was performed by varying weight of superstructure, pile length, pile head fixity, soil relative density with proposed numerical model. From the results of parametric study, it is identified that inertial force induced by superstructure is dominant on dynamic behavior of soil-pile-structure system and effect of kinematic force induced by soil movement was relatively small. Difference in dynamic behavior according to the pile length and pile head fixity was also numerically investigated.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.4
• /
• pp.291-302
• /
• 2004

Blast-induced anisotropic rock damage around a blast-hole was analyzed by a using numerical method with user-defined subroutine based on continuum damage mechanics. Anisotropic blasting pressure was evaluated by applying anisotropic ruck characteristics to analytical solution which is a function of explosive and rock properties. Anisotropic rock damage was evaluated by applying the proposed anisotropic blasting pressure. Blast-induced isotropic rock damage was also analyzed. User-defined subroutines to solve anisotropic and isotropic damage model were coded. Initial rock damages in natural ruck were considered in anisotropic and isotropic damage models. Blasting pressure and elastic modulus of rock were major influential parameters from parametric analysis results of isotropic rock damage. From the results of anisotropic rock damage analysis, blasting pressure was the most influential parameter. Anisotropic rock damage area in horizontal direction was approximately 34% larger and about 12% smaller in vertical direction comparing with isotropic rock damage area. Isotropic rock damage area under fully coupled charge condition was around 30 times larger than that under decoupled charge condition. Blasting pressure under fully coupled charge condition was estimated to be more than 10 times larger than that of decoupled charge condition.

• Explosives and Blasting
• /
• v.24 no.2
• /
• pp.41-50
• /
• 2006

Excavation by explosives blasting necessarily involves noise and vibration, which is highly prone to face claims on the environmental and structural aspects from the neighbors. When the blasting carried out in the vicinity of a structure, the effect of blasting vibration on the stability of the structure should be carefully evaluated. In the conventional method of evaluation, an equation for blast vibration is obtained from test blasting which is later used to determine the amount of charge. This method, however, has limitations in use since it does not consider topography and change in ground conditions. In order to overcome the limitations, dynamic numerical analysis is recently used in continuum or discontinuous models, where the topography and the ground conditions can be exactly implemented. In the numerical analysis for tunnels and rock slopes, it is very uncommon to simulate multi-hole blasting. A single-hole blasting pressure is estimated and the equivalent overall pressure at the excavation face is used. This approach based on an ideal case usually does not consider the ground conditions. And this consequently results in errors in calculation. In this presentation of a case study, a new approach of using blast waves obtained in the test blast is proposed. The approach was carried out in order to improve the accuracy in calculating blasting pressure. The stability of a structure in the vicinity of a slope blasting was examined using the newly proposed method.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.9
• /
• pp.1-11
• /
• 2016

The effect of the migration of nanoparticles near the wall of a channel on the convective heat transfer in a laminar flow of $SiO_2$ nanoparticle suspensions (nanofluids) under constant wall heat flux boundary conditions was numerically and experimentally investigated in this study. The dynamic thermal conductivity of the aqueous $SiO_2$ nanofluids was measured using T-type thermocouples attached to the outer surface of a stainless steel circular tube (with a length of 1 m and diameter of 1.75 mm). The nanofluids used in this study were synthesized by dispersing $SiO_2$ spherical nanoparticles with a diameter of 24 nm in de-ionized water (DIW). The enhancement of the thermal conductivity of the nanofluids (e.g., an increase of up to 7.9 %) was demonstrated by comparing the temperature profiles in the flow of the nanofluids with that in the flow of the basefluids (i.e., DIW). However, this trend was not demonstrated in the computational analysis, because the numerical models were based on continuum assumptions and flow features involving nanoparticles in a stable colloidal solution. Thus, to explore the non-continuum effects, such as the modification of the morphology caused by nanoparticle-wall interactions on the heat exchanging surfaces (e.g., the isolated and dispersed precipitation of the nanoparticles), additional experiments were performed using DIW right after the measurements using the nanofluids.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.18 no.2
• /
• pp.129-142
• /
• 2016

The underground cavity known as one of the reasons of ground surface settlement is a discontinuous character. Therefore, it is limited to analyze with continuum analysis. In this research, The spherical underground cavity affecting the ground surface settlement is studied with Discrete Element Method. Ground properties, depth and diameter of the spherical underground cavity are chosen as factors of the spherical underground cavity and the effect of the each factor variations on the ground surface settlement is analyzed. Relative depth to the diameter of the spherical underground cavity is also studied. The result of the research suggests the basis of underground cavity collapse prediction and standard of support.