• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.6
• /
• pp.699-705
• /
• 2011

High-maneuver missile is a projectile which demands a strong momentum at short time. To produce a necessary thrust for the flight, the gas of high temperature and pressure is generated through explosive combustion of solid propellant, and a great thrust can be obtained by expanding this high temperature and pressure gas. Although the operating time of a rocket motor is less than a few seconds, a failure of part or ablation near the throat of nozzle may take place during the expansion of high temperature and pressure gas for great thrust. In other words, for the precise control of a missile an exact stress analysis considering both, the thermal stress caused by the heat transfer between combustion gas and wall, and the mechanical stress caused by the pressure change in the flow, should be considered first. In this connection, this study investigated the safety, as a point of view of stress and melting point of the material, of the pre-designed thrust generating structure which is subjected to high temperature and pressure as a function of motor operating time.

• 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.

• The Journal of the Acoustical Society of Korea
• /
• v.43 no.4
• /
• pp.383-390
• /
• 2024

When measuring the radiated noise of an underwater vehicle, the range information between the vehicle and the receiver is an important factor, but since Global Positioning System (GPS) is not available in underwater, an alternative method is needed. As an alternative, the range is measured by estimating the arrival time, arrival time difference, and arrival frequency difference using a separate acoustic signal. However, errors occur due to the channel environment, and these outliers become obstacles in continuously measuring range. In this paper, we propose a method to reduce errors by curve fitting with a function in the form of a V-curve as a post-processing to remove outliers that occurred in the process of measuring range information. Simulation, lake and sea trials were conducted to verify the performance of the proposed method. In the results of the lake trial, the range estimation error was reduced by about 85 % from the Root Mean Square Error (RMSE) point of view.

• 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 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.4D
• /
• pp.387-393
• /
• 2010

This paper studies static crack growth of asphalt pavement using the extended finite element method (XFEM). To consider nonlinear characteristics of asphalt concrete, a viscoelastic constitutive equation using the Maxwell chain is used. And a linear cohesive crack model is used to regularize the crack. Instead of constructing the viscoelastic constitutive law from the Prony approximation of compliance and retardation time measured experimentally, we use a smooth log-power function which optimally fits experimental data and is infinitely differentiable. The partial moduli of the Maxwell chain from the log-power function make analysis easy because they change more smoothly in a more stable way than the ordinary method such as the least square method. Using the developed method, we can simulates the static crack growth test results satisfactorily.

• Journal of the Korean Geotechnical Society
• /
• v.35 no.5
• /
• pp.5-19
• /
• 2019

Considering the natural phenomenon in which steep slopes ($65^{\circ}{\sim}85^{\circ}$) consisting of rock mass remain stable for decades, slopes steeper than 1:0.5 (the standard of slope angle for blast rock) may be applied in geotechnical conditions which are similar to those above at the design and initial construction stages. In the process of analysing the stability of a good to fair continuum rock slope that can be designed as a steep slope, a general method of estimating rock mass strength properties from design practice perspective was required. Practical and genealized engineering methods of determining the properties of a rock mass are important for a good continuum rock slope that can be designed as a steep slope. The Genealized Hoek-Brown (H-B) failure criterion and GSI (Geological Strength Index), which were revised and supplemented by Hoek et al. (2002), were assessed as rock mass characterization systems fully taking into account the effects of discontinuities, and were widely utilized as a method for calculating equivalent Mohr-Coulomb shear strength (balancing the areas) according to stress changes. The concept of calculating equivalent M-C shear strength according to the change of confining stress range was proposed, and on a slope, the equivalent shear strength changes sensitively with changes in the maximum confining stress (${{\sigma}^{\prime}}_{3max}$ or normal stress), making it difficult to use it in practical design. In this study, the method of estimating the strength properties (an iso-angle division method) that can be applied universally within the maximum confining stress range for a good to fair continuum rock mass slope is proposed by applying the H-B failure criterion. In order to assess the validity and applicability of the proposed method of estimating the shear strength (A), the rock slope, which is a study object, was selected as the type of rock (igneous, metamorphic, sedimentary) on the steep slope near the existing working design site. It is compared and analyzed with the equivalent M-C shear strength (balancing the areas) proposed by Hoek. The equivalent M-C shear strength of the balancing the areas method and iso-angle division method was estimated using the RocLab program (geotechnical properties calculation software based on the H-B failure criterion (2002)) by using the basic data of the laboratory rock triaxial compression test at the existing working design site and the face mapping of discontinuities on the rock slope of study area. The calculated equivalent M-C shear strength of the balancing the areas method was interlinked to show very large or small cohesion and internal friction angles (generally, greater than $45^{\circ}$). The equivalent M-C shear strength of the iso-angle division is in-between the equivalent M-C shear properties of the balancing the areas, and the internal friction angles show a range of $30^{\circ}$ to $42^{\circ}$. We compared and analyzed the shear strength (A) of the iso-angle division method at the study area with the shear strength (B) of the existing working design site with similar or the same grade RMR each other. The application of the proposed iso-angle division method was indirectly evaluated through the results of the stability analysis (limit equilibrium analysis and finite element analysis) applied with these the strength properties. The difference between A and B of the shear strength is about 10%. LEM results (in wet condition) showed that Fs (A) = 14.08~58.22 (average 32.9) and Fs (B) = 18.39~60.04 (average 32.2), which were similar in accordance with the same rock types. As a result of FEM, displacement (A) = 0.13~0.65 mm (average 0.27 mm) and displacement (B) = 0.14~1.07 mm (average 0.37 mm). Using the GSI and Hoek-Brown failure criterion, the significant result could be identified in the application evaluation. Therefore, the strength properties of rock mass estimated by the iso-angle division method could be applied with practical shear strength.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.16 no.4
• /
• pp.403-415
• /
• 2014

Recently, the temperature rise in the summer due to climate change, power usage is increasing rapidly. As a result, power generation facilities have been newly completed and the need for ultra-high-voltage transmission line for power transmission of electricity to the urban area has increased. The mechanized tunnelling method using a shield TBM have an advantage that it can minimize vibrations transmitted to the ground and ground subsidence as compared with the conventional tunnelling method. Despite the popularity of shield TBM for cable tunnel construction, study on the mechanical behavior of cable tunnel driven by shield TBM is insufficient. Thus, in this study, the effect of fractured zone ahead of tunnel face on the mechanical behavior of the shield TBM cable tunnel is investigated. In addition, it is intended to compare the behavior characteristics of the fractured zone with continuous model and applying the interface elements. Tunnelling with shield TBM is simulated using 3D FEM. According to the change of the direction and magnitude of the fractured zone, Sectional forces such as axial force, shear force and bending moment are monitored and vertical displacement at the ground surface is measured. Based on the stability analysis with the results obtained from the numerical analysis, it is possible to predict fractured zone ahead of the shield TBM and ensure the stability of the tunnel structure.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1994.10a
• /
• pp.116-120
• /
• 1994

최근 항공우주 및 생산자동화 분야의 급격한 발달에 따라 정밀도와 효율성을 향상시키기 위해 저중량, 고강도 구조물이 요구된다. 그러나 경량화 추세로 인해 수반되는 구조물의 유연성 증가로 외력에 대한 구조응답의 진폭이 커지고 구조물의 피로 수명이 단축되어 매우 위험한 상황에 이를 수 있다. 이런 바람직하지 않은 진동현상을 제어하기 위해 여러 제어이론을 응용한 진동억제시스템의 연구가 활발하며, 신소재인 압전재료의 개발로 새로운 방향이 제시되고 있다. 압전재료는 유연한 구조물에 부착되어 압전재료의 수축, 팽창 운동에 의해 발생된 에너지를 부착된 구조물에서의 제어력으로 사용되어, 진동 혹은 자세 및 형상 제어에 활용되고 있다. 압전재료에 대한 연구로는 Crawley, de Luis3가 보의 표면 혹은 내부에 압전세라믹을 부착하여 액튜에이터로 사용하는 경우 집중모멘트를 가하는 역할을 함을 밝혔고, Hanagud, obal은 압전재료를 센서와 액튜에이터로 사용해 복합재료 보에 대한 최적 진동제어 알고리즘을 개발, 그 성능에 대한 효과를 조사하였고 임의의 위치에 부착된 센서 및 액튜에이터를 고려한 복합재료 보의 운동방정식을 유한요소법을 이용 유도하였으며 변위율 피드백(rate feedback)과 모달피드백(modal feedback) 제어기를 적용하여 진동제어 효과를 고찰하였다. 그리고 Tomas, James, Hubbard는 압전필름을 액튜에이터로 사용해 복합재료 보에 Liapunov 제어기와 변위율 피드백 제어기를 사용하여 능동감쇠기를 설계하였고, Lee, Chaing, Sullivan은 압전필름을 센서와 액튜에이터로 사용해 평판에 변위율 피드백 제어기를 적용한 능동감쇠기를 설계하고 실험적으로 수행하였다. Base가 백색잡음가진을 지속적으로 받을 때 보끝의 움직임이 최소가 되도록 제어하고자 연구를 수행 중인 바 그 결과로소 본 논문에서는 적용시켰고 F-P-K 방정식을 이용해 확률영역으로 변환하여 LQR 제어기와 pole allocation 제어기를 시스템에 적용시켜 우수한 특성을 갖음을 제어 시뮬레이션의 결과를 통해 입증하였다.