• Explosives and Blasting
• /
• v.39 no.2
• /
• pp.1-14
• /
• 2021

Recently, with the development of high-performance processing devices such as GPGPU, a three-dimensional dynamic analysis technique that can replace expensive rock material impact tests has been actively developed in the defense and aerospace fields. Experimentally observing or measuring fracture processes occurring in rocks subjected to high impact loads, such as blasting and earth penetration of small-diameter missiles, are difficult due to the inhomogeneity and opacity of rock materials. In this study, a three-dimensional dynamic fracture process analysis technique (3D-DFPA) was developed to simulate the fracture behavior of rocks due to impact. In order to improve the operation speed, an algorithm capable of GPGPU operation was developed for explicit analysis and contact element search. To verify the proposed dynamic fracture process analysis technique, the dynamic fracture toughness tests of the Straight Notched Disk Bending (SNDB) limestone samples were simulated and the propagation of the reflection and transmission of the stress waves at the rock-impact bar interfaces and the fracture process of the rock samples were compared. The dynamic load tests for the SNDB sample applied a Pulse Shape controlled Split Hopkinson presure bar (PS-SHPB) that can control the waveform of the incident stress wave, the stress state, and the fracture process of the rock models were analyzed with experimental results.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.67-78
• /
• 2012

A three-dimensional approximate computer-based method, YSPR (Yonsei Piled Raft), was developed for analysis of behavior of piled raft foundations. The raft was modeled as a flat shell element having 6 degrees of freedom at each node and the pile was modeled as a beam-column element. The behaviors of pile head and soil were controlled by using $6{\times}6$ stiffness matrix. To model the non-linear behavior, the soil-structure interaction between soil and pile was modeled by using nonlinear load-transfer curves (t-z, q-z and p-y curves). Comparison with previous model and FEM analysis showed that YSPR gave similar load-displacement behaviors. Comparison with field measurement also indicated that YSPR gave a reasonable result. It was concluded that YSPR could be effectively used in analysis and design of piled raft foundations.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.12
• /
• pp.59-69
• /
• 2018

In order to use the limit equilibrium theory, it is necessary to find a slip line under the ultimate failure condition. The strength reduction method using the Lagrangian finite element method defines the ultimate failure state at a time when the numerical solution cannot converge within the certain number of the iteration. When the coupled Eulerian-Lagrangian (CEL) method is used, however, such definition is inappropriate because the numerical solution of the CEL method can converge even under the ultimate failure condition. In this study, an objective condition designating the ultimate failure state in the finite element analysis adopting the CEL method was proposed. In the problem of the bearing capacity of the undrained soft ground subjected to the strip footing loading, we found that the rate of the plastic dissipated energy is highly sensitive at the load of the theoretical limit of the ultimate failure state.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.2
• /
• pp.83-91
• /
• 2018

The noise sources of a tip-jet driven rotor can be separated by rotor blade noise and jet noise. The rotor blade noise consists of thickness noise, loading noise, nonlinear quadrupole noise, and jet noise is divided into nozzle momentum noise and jet radiation noise. The flow analysis for the prediction of rotor blade noise is performed by CFD (Computational Fluid Dynamics) analysis, and the noise source of the rotor blade noise is identified by simultaneously applying the permeable and impermeable surface based FW-H (Ffowcs Williams-Hawkings) acoustic analogy. The nozzle momentum noise is obtained by permeable surface FW-H, and jet radiation noise is predicted by using empirical method for the fixed-wing jet. Both of jet noises use nozzle exit condition for noise analysis. The accuracy of the technique is verified based on the noise measurements of the tip-jet driven rotor, and the unique noise characteristics of the tip-jet driven rotor is confirmed by spectrum analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.4
• /
• pp.461-466
• /
• 2013

Advanced fiber-reinforced laminated composites are widely used in various fields of engineering to reduce weight. The material property of each ply is well known; specifically, it is known that ply is less reliable than metallic materials and very sensitive to the loading direction. Therefore, it is important to consider this uncertainty in the design of laminated composites. In this study, reliability analysis is conducted using COMSOL and MATLAB interactions for a laminated composite plate for the case in which the tip deflection is the design requirement and the material property is a random variable. Furthermore, the efficiency and accuracy of the approximation method is identified, and a probabilistic sensitivity analysis is conducted. As a result, we can prove the applicability of the advanced design method for the stabilizer of an underwater vehicle.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.1
• /
• pp.71-83
• /
• 2004

The purpose of this study is to present an analysis method for the prediction of the changes m ground conditions. To this end, three dimensional convergence displacements are analyzed in several ways to estimate the trend of displacement changes. Three-dimensional arching effect is occurred around the unsupported excavation surface including tunnel face when a tunnel is excavated in a stable rock mass. If the ground condition ahead of tunnel face changes or a weak zone exists, a diagnostic trend of displacement change is observed by the 3 dimensional measurement and numerical analysis. Therefore, the change of ground condition and the existence of a weak zone ahead of tunnel face can be predicted by monitoring 3-dimensional absolute displacements during excavation, and applying the methodology (the ratio of L/C, $C/C_o$, etc.) presented in this study.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.4C
• /
• pp.159-167
• /
• 2012

In this study, laboratory model pile-load test and finite element analysis were carried out to compare and analyze the strain behaviour around the model pile tip. In order to simulate the pile load, both the LCM(load control method)and DCM(displacement control method) were introduced to determine which one is appropriate for the FE simulation. In contrast to the previous simulation method, two interface elements around the model pile were used to consider the slip effect in the finite element analysis and its results were compared to the model test. Through this study it was found that the degree of non-associated flow was a dominant factor in terms of numerical solution convergence. In addition, an improved FE mesh was required to obtain the symmetric distribution of the maximum shear strain contour.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.21 no.3
• /
• pp.102-113
• /
• 2017

A new index for the buckling strength of non-structured plates is proposed. The external work or the deformation energy caused by the external loads or the boundary displacement controled by a load parameter is calculated along an equilibrium path of the member under consideration. If the second variation of the energy with respect to the parameter loses its positiveness, it defined as the limit of the stability. In contrast to the current method given in codes where the stability limit is evaluated by using only representative internal forces, the evaluation of the stability limit is always consistent even with the change of the distribution of the internal forces on the boundary. If the elasticity is concerned, the result from this proposed approach becomes identical to that from the classical methods.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.2
• /
• pp.86-97
• /
• 2005

Accurate estimation of the wave-induced pore water pressure in the seabed is key factor in studying the stability of the seabed in the vicinity of coastal structure. Most of the existing numerical models for wave structure seabed interaction have been linked through applying hybrid numerical technique which is analysis method separating the wave field and seabed regime. Therefore, it is necessary to develope a numerical model f3r simulating accurately wave$\cdot$structure$\cdot$ seabed interaction under wave loadings by the single domain approach for wave field and seabed regime together. In this study, direct numerical simulation is newly proposed. In this model, modeled fluid drag has been used to detect the hydraulic properties according to the varied geometrical shape inside the porous media by considering the turbulence resistance as well as laminar resistance. Contrary to hybrid numerical technique, direct numerical simulation avoids the explicit formulation of the boundary conditions at the fluid/porous media interface. A good agreement has been obtained by the comparison between existed experimental results by hydraulic model test and direct numerical simulation results far wave $\cdot$structure$\cdot$seabed interaction. Therefore, the newly proposed numerical model is a powerful tool for estimating the nonlinear dynamic responses among a structure, its seabed foundation and water waves.

• Computational Structural Engineering
• /
• v.4 no.1
• /
• pp.20-27
• /
• 1991

본 논문에서는 유체-구조 상호작용해석의 일종의 3차원 접수구조물의 진동해석을 효과적으로 수행하기 위한 해석방법을 제시하기 위하여 동적재해석기법을 검토하였다. 접수구조물의 유한구조 상호작용해석 결과는 구조진동의 관심 주파수역에서는 3차원 연성 부가수질량으로 표현되는 관성력으로 나타난다. 따라서 구조질량행렬에 부가수질량 행렬이 더해져서 전체 관성력으로 표현된다. 이 부가수질량을 추가질량으로 보고 재해석기법을 응용하는 방법을 수치실험을 통해 검증하였다. 이 때 재해석기법이 갖추어야 할 조건은 원구조의 질량과 거의 같은 정도의 질량이 추가되고 또한 완전 연성질량이 추가된 경우에도 정확한 해를 주어야 한다는 것이다. 이를 검증하기 위해 직접재해석기법과 섭동법을 이용한 재해석기법으로 4질량 스프링지지구조에 대한 해석을 수행한 결과 직접재해석기법의 응용이 적합함을 쉽게 입증할 수 있었다. 접수구조물의 예로는 3차원 잠수주상체에 대해 접수진동해석을 수행하였으며 그 결과 선체진동해석에 전통적으로 이용되고 있는 2차원 부가수질량과 3차원 수정계수를 사용한 기준차수법에서는 수지모드와 수평-비틔 연성모드와 같이 서로 독립적인 모드에 대해서는 따로 진동해석을 수행해 주어야 하는 단점이 발견되었다. 이 단점을 보완한 각 모드의 3차원 수정계수행렬을 이용한 재해석기법을 도입하여 모드에 상관없이 동시에 해를 구할 수 있었다. 그러나, 이 방법은 3차원 수정계수가 구해져 있는 경우에 한해서만 적용가능하며 실제 선체진동의 경우에는 10Hz 미만의 저차 주선체 진동에 한해서만 적용가능한 방법이다. 고차의 진도옴드에는 3차원 수정계수를 구할 수 없기 때문에 유체-구조 상호작용 해석결과로부터 얻은 3차원 연성 부가수질량을 이용하게 되며 이 때 이 행렬이 접수구조 표면의 전 자유도와 연성되어 있기 때문에 방대한 방정식을 푸어야 하지만 직접재해석기법을 적용함으로써 정확한 해를 구할 수 있었다. 또한 3차원 부가수질량을 이용한 직접재해석기법은 종래의 2차원 부가수질량과 3차원 수정계수를 이용한 방법에 비해 해석시간 면에서도 전혀 불리한 점이 없는 경제적 방법임이 밝혀졌다. 앞으로 Slamming 혹은 수중폭파 등의 충격하중에 의한 천이 구조응답 해석을 위한 효과적인 방법에 대해서도 연구결과를 발표할 계획이다.