• 한국지반신소재학회논문집
• /
• 제10권2호
• /
• pp.1-12
• /
• 2011

조립재료의 입자 결합 및 파쇄 형태가 전단거동 특히 잔류 전단거동 특성에 미치는 영향을 분석하기 위하여, 개별요소법(DEM, discrete element method)에 기초를 둔 프로그램인 PFC(Particle Flow Code)를 이용하여 링 전단시험을 수치해석적으로 모델링하였다. 본 연구에서는 PFC내의 clump 모델 및 cluster 모델을 이용하여 두 개의 비파쇄모델 그리고 두 개의 파쇄모델을 포함한 총 네 개의 모델을 제시하였다. Lobo-Guerrero and Vallejo(2005)가 제안한 Lobo-crushing 모델의 적합성을 검토하였다. 또한 링 전단시험 모델링의 결과 분석을 통하여 직접전단시험 모델링 결과와 비교하였다. 연구 결과, 잔류 전단거동 분석을 위해서는 링 전단시험의 모델링이 필수적임을 알 수 있었다. 또한 잔류 전단강도 분석을 위해서는 Lobo-crushing 모델이 부적합함을 알 수 있었다. 따라서 본 연구에서 제시한 수치해석 모델은 향후 입자 파쇄를 포함한 조립재료의 잔류 전단강도 특성 연구에 다양하게 적용될 수 있다고 판단된다.

• 한국지반공학회논문집
• /
• 제25권10호
• /
• pp.41-53
• /
• 2009

조립재료의 입자 형상이 입자 파쇄 전개 및 전단 거동 특성에 미치는 영향을 분석하기 위하여, 개별요소법(DEM, discrete element method)을 이용하여 직접전단시험을 수치해석적으로 모델링하였다. PFC(Particle Flow Code)내의 clump 모델 및 cluster 모델을 이용하여 6가지 형상의 입자를 생성하여 이를 원형입자의 직접 전단거동과 비교 분석함으로써 입자형상의 영향을 연구하였다. 연구결과, PFC에 의해 모델링된 직접 전단모델의 수치해석 결과는 실내 실험결과와 잘 일치하였으며, 따라서 본 연구 결과의 타당성을 입증하였다. 입자 형상 관점에서 모나고 거친 입자의 내부마찰각이 상대적으로 둥글고 매끄러운 입자에 비해 큰 값을 나타냈으며, 입자 파쇄 또한 많이 발생하는 것을 확인하였다. 이때 입자파쇄는 전단면근처에 집중되며 전단대를 형성하였다. 따라서 본 연구에서 제시한 수치해석 모델은 향후입자 파쇄를 포함한 조립재료의 전단강도 특성 연구에 다양하게 적용될 수 있다고 판단된다.

• 한국터널지하공간학회 논문집
• /
• 제12권1호
• /
• pp.31-41
• /
• 2010

본 연구에서는 실험 및 수치해석적인 접근방법을 통하여 대심도 과지압 구간에서 발생하는 암반의 스폴링(spalling) 및 팽창모드에 대한 모델링 기법을 연구하였다. 이에 대한 실험적 접근 방법으로서 축 방향 압축을 받는 직사각형 인공암석보(beam)에 4점 휨 시험을 결합한 축방향 압축 휨 시험을 수행하여 대심도 지하공간의 응력모드와 유사한 조건 하에서의 암석의 균열 팽창 및 스폴링 과정을 고찰하였다. 또한, 수치해석적 접근방법으로서 기존의 연속체 해석으로는 모사하기 힘든 암석의 균열과정 및 팽창특성을 개별 입자해석 프로그램인 PFC2D를 이용하여 모델링 하였다. 본 연구 결과 휨 실험에서 구한 팽창시점은 스폴링에 요구되는 응력수준을 평가하는데 중요한 지표가 됨을 알 수 있었으며, 또한 수치해석 결과도 유사한 결과를 모사할 수 있음을 보여주었다.

• 터널과지하공간
• /
• 제33권1호
• /
• pp.57-69
• /
• 2023

본 연구에서는 2차원 수리-역학적연계 개별요소모델링(DEM)을 사용하여 네델란드 그로닝엔(Groningen) 천연가스전 저류층의 유발지진을 모사하였다. 수치해석 코드는 ITASCA社의 상용프로그램인 PFC2D (Particle Flow Code 2D)를 사용하였으며 본 수치해석 연구에 적용하기 위해 수리-역학적 연계 모델 외 1) 비균질 저류층 압력분포 초기화, 2) 비선형 압력-시간이력 경계조건, 3) 국소 응력 분포 계산 등의 개별모듈을 추가개발, 적용하였다. 그로닝엔 가스전에 분포하는 복잡한 단층 형상을 포함하는 40 × 50 km² 크기의 2차원 모델을 생성하였고, 1960년부터 2020년까지 약 60년 동안의 가스생산, 즉 압 력저하로 인한 단층의 파괴거동을 모사하였다. 유발지진의 시공간적 발생을 수치해석모델로 재현하였고 그 발생 메커니즘을 규명하였다. 또한 저류층 압축으로부터 지표에서의 지반침하의 분포를 예측하였고 그로닝엔에서의 실측자료 사이에 유사성을 확인할 수 있었다. 이를 통해 본 연구에서 소개한 2차원 수리-역학적연계 개별요소모델링(DEM)의 복잡한 지질조건과 수리-역학적 연계 프로세스에 의한 단층거동을 구현할 수 있는 툴(tool)로서의 활용성을 확인하였다.

• 한국인터넷방송통신학회논문지
• /
• 제16권3호
• /
• pp.203-211
• /
• 2016

유전자 코드는 바이오 정보 처리에 키 포인트로 인체의 유기적인 조직체이다. 현대 과학에서는 유전자 코드 분자구조의 신비스러운 특성을 체계적으로 설명하고 이해하는데 연구가 집중되고 있다. 본 논문에서는 유전자 시스템을 대칭적으로 해석하는데 중점을 두었고, Jacket 행렬로 무잡음 RNA 유전자 코드를 가장 단순하게 해석했다. 이유는 Jacket 행렬과 RNA는 그 역행렬이 Element (Block)-wise Inverse로 그 역(Inverse)도 자신이란 점과 대칭적 성질, 그리고 Kronecker곱을 갖기 때문이다. 제안된 방법이 유전자 RNA 코돈(Codon : 괘(卦))의 견지에서 Jacket 행렬의 분해를 통해 간단하고 명료함을 보인다.

• Structural Engineering and Mechanics
• /
• 제13권2호
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• 한국도로학회논문집
• /
• 제16권2호
• /
• pp.43-52
• /
• 2014

PURPOSES : To solve problems in current compaction control DCPT(Dynamic Cone Penetrometer Test), highly correlated with various testing methods, simple, and economic is being applied. However, it、s hard to utilize DCPT results due to the few numerical analyses for DCPT have been performed and the lack of data accumulation. Therefore, this study tried to verify the validation of numerical modeling for DCPT by comparing and analyzing the results of numerical analyses with field tests. METHODS: The ground elastic modulus and PR(Penetration Rate) value were estimated by using PFC(Particle Flow Code) 3D program based on the discrete element method. Those values were compared and analyzed with the result of field tests. Also, back analysis was conducted to describe ground elastic modulus of field tests. RESULTS : Relative errors of PR value between the numerical analyses and field tests were calculated to be comparatively low. Also, the relationship between elastic modulus and PR value turned out to be similar. CONCLUSIONS : Numerical modeling of DCPT is considered to be suitable for describing field tests by carrying out numerical analysis using PFC 3D program.

• Structural Engineering and Mechanics
• /
• 제68권5호
• /
• pp.537-547
• /
• 2018

This research presents the effect of anisotropy of the hollow disc mode under Brazilian test using PFC3D. The Brazilian tensile strength test was performed on the hollow disc specimens containing the bedding layers and then these specimens were numerically modeled by using the two dimensional discrete element code (PFC3D) to calibrate this computer code for the simulation of the cracks propagation and cracks coalescence in the anisotropic bedded rocks. The thickness of each layer within the specimens varied as 5 mm, 10 mm and 20 mm and the layers angles were changed as $0^{\circ}$, $25^{\circ}$, $50^{\circ}$, $75^{\circ}$ and $90^{\circ}$. The diameter of internal hole was taken as 15 mm and the loading rate during the testing process kept as 0.016 mm/s. It has been shown that for layers angles below $25^{\circ}$ the tensile cracks produce in between the layers and extend toward the model boundary till interact and break the specimen. The failure process of the specimen may enhance as the layer angle increases so that the Brazilian tensile strength reaches to its minimum value when the bedding layers is between $50^{\circ}$ and $75^{\circ}$ but its value reaches to maximum at a layer angle of $90^{\circ}$. The number of tensile cracks decreases as the layers thickness increases and with increasing the layers angle, less layer mobilize in the failure process.

• 터널과지하공간
• /
• 제14권4호
• /
• pp.248-260
• /
• 2004

응력 증가에 의한 취성 암석의 손상은 미세균열의 개시로부터 시작하여 각 개별 균열들의 전파 및 결합에 의해 거시적인 파괴면을 발생시킨다. 전통적으로 암반의 손상 및 파괴현상을 설명하기 위해 거시적인 파괴 기준이나 탄소성 모델과 같은 연속체적인 접근법이 주류를 이루어왔다. 하지만 개별적인 균열들의 개시와 전락 과정을 명시적으로 고려할 수 있다면 현상론적인 관점에서 보다 실제에 가까운 암석 손상 및 파괴 과정을 재현할 수 있을 것이다. 본 연구에서는 암석의 균열 진전 모델링을 위해 개발된 경계요소 코드인 FRACOD를 이용하여 암석의 손상 및 파괴 과정을 모사한 결과를 제시한다. 수치일축압축시험을 통해 개발된 모델의 적정성을 검증하고 암반의 치수효과를 고려한 현실적인 암석 파괴 과정을 재현하였다. 또한 이러한 접근법의 적용 사례로서, 실제 굴착이 진행중인 심부 수갱 암반 주변에서 심도와 암반 특성에 따라 균열 진전과 이에 따른 암반 손상의 범위를 예측한 결과를 제시하였다. 이 접근법은 취성도가 큰 암반에서 발생하는 안정성 문제에 대한 공학적인 해법을 찾는데 기여를 할 수 있을 것으로 기대된다.

• 화약ㆍ발파
• /
• 제41권3호
• /
• pp.62-72
• /
• 2023

1970년대 산업화에 따른 플랜트 구조물이 노후화됨에 따라 구조적 기능을 상실해 발파 해체 공법을 활용한 해체 철거 수요가 증가하고 있다. 발파 해체공법은 기계식 해체공법에 비해 해체 공기가 짧아 환경공해 발생 노출 시간을 최소화할 수 있지만 잘못된 발파 설계 및 시공 계획에 따른 붕괴거동의 실패는 안전성에 매우 큰 위험을 유발한다. 따라서 붕괴거동 모사를 통해 최적의 발파 해체 조건과 이에 따른 영향을 고려하는 것이 중요하다. 본 연구에서는 Finite element method (FEM)와 Discrete element method (DEM)의 장점을 활용해 구축된 3-D Combined finite discrete element method (FDEM) 코드 기반 3-D DFPA 를 적용해 (구)서천화력발전소의 연돌 구조물에 대한 해체 모사를 수행하였으며 실제 구조물의 연돌 구조물 발파해체의 붕괴거동과 비교 분석하였다. 수치 모사 결과, 실제 구조물과 붕괴 거동 및 붕괴 완료 시간이 동일하게 나타났다. 또한, 발파구간 개구부 상부에 위치한 후드부의 크기가 연돌의 붕괴거동에 미치는 영향을 분석하기 위해 후드부의 면적을 조정하여 해체 모사를 수행하고 균열 발생 양상 및 z-방향 변위 곡선을 통해 비교 분석하였다. 분석 결과, 후드부의 면적이 증가함에 따라 하중을 지지하는 면적이 줄어들고 그에 따른 균열 발생 증가 및 전도 시간 감소를 확인하였다.