• Journal of the Korean GEO-environmental Society
• /
• 제12권7호
• /
• pp.57-65
• /
• 2011

Numerical simulations by three-dimensional Particle Flow Code($PFC^{3D}$, Itasca) considering distinct element method (DEM) were carried out for prediction of triaxial compression test with sand material. The effect of scale conditions for numerical model and distinct material on final prediction results was analyzed by numerical models under various scale conditions, and following observations were made from the numerical experiments. It is very useful to model the initial material condition without any porosity conversion from 2-D to 3-D DEM. Numerical experiments have shown that in all cases considered, 3D distinct element modeling could provide good agreement on stress-strain behavior, volume change and strength properties with laboratory testing results. It was important thing to assess reasonable scale ratio of numerical model and distinct elements for saving calculation time and securing calculation efficiency under condition with accuracy and appropriateness as numerical laboratory. As results of DEM simulations under various scale conditions, most of results show that shear strength properties as cohesion and internal friction angle are similar in condition of $D_{mod}/D_{gmax}$ < 10. It shows that 3-D distinct element method could be used as efficient tool to assess strength properties by numerical laboratory technique.

• Proceedings of the Korea Water Resources Association Conference
• /
• 한국수자원학회 2020년도 학술발표회
• /
• pp.305-305
• /
• 2020

2019년 10월에 발생한 태풍 미탁은 강원도 삼척시 근덕면 궁촌리 기준 시간당 최대 129mm/hr, 누적강수량 488.5mm를 기록하며 폭우와 강풍, 홍수, 토석류 및 산사태로 사망 13명, 실종 2명, 이재민 910세대 1,442명의 인명피해와 공공시설 1,835건, 사유시설 3,700건 등의 재산 피해가 발생하였다. 한편 Hyper KANAKO 모형은 일본에서 개발된 토석류 수치 모델로 오픈소스인 QGIS와 연동되어 도구의 형태로 사용할 수 있어 사용자에게 편리함을 제공하고 토석류 해석 시 침식 및 퇴적작용의 해석이 가능하다. 1차원의 지형(계류)과 2차원의 지형(퇴적부)의 지반설정과 사방댐의 종류(Closed, Slit, Landslide) 및 개수, 위치, 높이 등을 쉽게 조절할 수 있다는 장점이 있다. 또한, DEM을 기반으로 변환 프로그램을 사용하여 구축한 자료인 LP(Laser Profiler) 데이터를 바탕으로 3차원의 지형을 볼 수 있어 2차원의 지형보다 사용자의 직관적인 해석이 가능하다. 본 연구에서는 QGIS와 Hyper KANAKO를 이용하여 태풍 미탁으로 토석류 피해가 발생한 강원도 삼척시 신남지역을 대상으로 DEM과 LP 데이터를 구축하고 수치모델링을 실시하여 현장조사 결과와 비교분석 하였다.

• Journal of the Korean Geotechnical Society
• /
• 제34권8호
• /
• pp.15-26
• /
• 2018

Recently, as the seismic performance based design methods have been introduced, dynamic numerical analyses need to be performed to evaluate the actual performance of structures under earthquakes. The verification of the numerical modeling is the most important for the performance based design. Therefore, 2-dimensional numerical analyses were performed to simulate the seismic behavior of a pile-supported structure, to provide the proper numerical modeling and to determine of input parameters. A dynamic centrifuge test of a pile group in dry loose sand was simulated to verify the applicability of the numerical model. The numerical modeling was carefully made to reflect the actual condition of the centrifuge test including dynamic soil properties, soil-pile interaction, boundary condition, the modeling of the group pile and structure and so on. The predicted behavior of the numerical analyses successfully simulated the acceleration variation in ground, the moment and displacement of the pile, and the displacement and acceleration of the structure. Therefore, the adopted numerical modeling and the input parameters can be used to evaluate the seismic performance of pile groups.

• Geophysics and Geophysical Exploration
• /
• 제2권1호
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed, when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• The Journal of Engineering Geology
• /
• 제6권3호
• /
• pp.119-130
• /
• 1996

The theory and numerical technique using boundary elements method (BEM) are given to solve 2-dimensional resistivity problems. Potential distributions from homogeneous resistivity model and layered model are calculated by using BEM for a point source of current injection. The potential distributions of BEM are compared with those of finite difference method (FDM) and finite elements method (FEM). Among the three numerical methods to solve 2-dimensional resistivity problem, it is proved that BEM is more efficient tool than FDM and FEM in consideration of computing storage and time as weU as the accuracy of solutions.

• Geophysics and Geophysical Exploration
• /
• 제16권1호
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

• Proceedings of the Korea Water Resources Association Conference
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.419-423
• /
• 2005

2차원 하천흐름. 유사이동, 오염확산 해석을 위해서는 유한요소법(FEM) 등을 활용한 수치해석 모델이 사용되며, 이 때 모델링을 위해서 모의영역은 2차원의 요소망으로 구성하게 된다. 기존 국내 연구에서는 이러한 2차원 요소망 생성 및 수치해석을 위해 SMS 및 CCHE2D와 같은 외국의 상용 프로그램을 이용하였으나, "수자원의 지속적 확보기술개발 사업"의 일환으로 수행중인 "RAMS(River Analysis and Modeling System) 개발" 과제를 통해 순수 국산 2차원 수리해석 소프트웨어가 개발 중에 있다. 본 연구에서는 RAMS를 비롯한 기타 수리해석 모델에서 사용 가능한 2차원 범용 요소망 생성 프로그램을 개발하고자 다양한 요소망 생성 기법 및 국외 상용 2차원 수리해석 모델 프로그램에서의 메쉬 생성방법 등을 분석하여 프로토타입 메쉬 생성기를 개발하였다. 현재 개발 중인 요소망 생성 프로그램은 요소망 생성을 위한 기본 기능 외에도 생성된 요소망 데이터를 RAMS에 포함되어 있는 다양한 2차원 유한요소 모형의 입력자료 형태로 변환함으로써 수리해석에 적용가능하다. 이러한 범용 2차원 요소망 생성 프로그램의 개발은 안정적인 수치해석의 기반을 제공하고 다양한 하천흐름, 유사, 수질 해석 모델과 연계함으로써 하천의 수리학적인 거동을 보다 정확히 모의할 수 있을 것으로 기대된다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제29권3호
• /
• pp.285-292
• /
• 2016

In this study, FE-BE direct coupling methods of 1D and 2D problems are considered for the pontoon-type floating structure and the difference of the modeling dimensions is investigated for the hydroelastic response. The modeling dimensions are defined as the 1D problem consisting 1D beam-2D fluid coupling and the 2D problem consisting 2D plate-3D fluid coupling with zero-draft assumption. For case studies, hydroelastic responses of the 1D Problem are compared to those of the 2D Problem for a wide range of aspect ratio and regular waves. It is shown that the effects of the elastic behavior are increased by decreasing the incident wavelength, whereas the effects of the rigid behavior are increased by increasing the incident wavelength. In 2D problem, the incident wave angle can be considered, and slightly more accurate results can be obtained, but the computational efficiency is lower. On the other hand, in 1D problem with plate-strip condition, the incident wave angle cannot be considered, but when the aspect ratio is large, the overall responses can be analyzed through a simplified model, and the computational efficiency can be improved.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 한국복합재료학회 1999년도 추계학술발표대회 논문집
• /
• pp.107-110
• /
• 1999

본 연구에서는 진공압추출법에 의해 제조되는 알루미나(AI$_2$O$_3$) 단섬유예비성형체의 제조공정시 발생하는 섬유들의 배향분포의 특징을 2차원 수치모사를 통해 밝혔다. 수치해석은 섬유의 배향분포를 2차원으로 가정한 후 제조공정시 섬유들이 적층되는 현상을 섬유의 낙하, 회전, 미끄러짐, 이동의 단계의 나누어 모델링하였다. 해석결과 섬유의 배향과 부피분율은 섬유의 종횡비(길이/지름)에 의해 크게 영향을 받으며, 배향 분포는 정규분포와 유사한 것으로 밝혀졌다.

• Economic and Environmental Geology
• /
• 제27권2호
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.