• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.47-50
• /
• 2009

본 논문에서는 불균질 재료의 균열진전을 해석하기 위한 방법으로 변절점 유한요소를 이용한 멀티스케일 기법을 제시하였다. 효율적인 해석을 위하여 서로 다른 스케일의 요소망을 적용하여 전체 모델의 자유도를 감소시킨다. 균열선단과 비교적 멀리 떨어져 있는 영역은 균질화 기법을 도입하여 불균질 재료에 대한 등가물성을 갖는 성긴 요소망으로 대체하고, 균열선단 주변의 요소망은 재료의 기하학적 특성과 불균질성을 반영하도록 조밀하게 구성한다. 한편 균열선단에 존재하는 응력 특이성을 표현하기 위하여 균열선단을 포함한 요소를 더욱 조밀한 요소망으로 분할하여 구성한다. 여기에서 서로 다른 스케일의 요소망 경계에는 변절점 유한요소를 적용함으로써 경계에서의 절점 연결조건과 적합성을 만족시킬 수 있다. 제시한 멀티스케일 기법을 수치예제에 적용함으로써 정확성과 효율성을 검증하였으며, 특히 불균질 성분이 균열진전에 미치는 영향을 경계조건과 T-응력의 관점에서 분석하였다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.51-54
• /
• 2009

본 논문에서는 MLS기반 유한요소에 대한 개략적인 현재 개발상황과 향후 예상할 수 있는 응용분야에 대한 제안을 하였다. 이동최소제곱근사를 이용하여 형상함수를 생성하는 MLS기반 유한요소는, 요소의 경계에서 기존 유한요소의 성질-크로네커 델타 조건-을 가지면서도 기존 요소가 갖지 못했던 임의의 절점추가가 자유롭다는 장점이 있어 다양한 변절점요소로의 개발이 이루어져왔다. 선형 또는 이차형상함수를 갖는 2차원 변절점요소 뿐 아니라, 균열선단과 균열면을 포함하고 있는 2차원 균열요소와 3차원에서의 제한적인 변절점요소 등이 개발되어 다양한 불연속성 문제에 적용 가능함이 입증되었다. 이러한 MLS기반 유한요소는 향후 2차원 변절점 3각요소, 2차원 삼각균열요소, 변절점 쉘요소, 균열 쉘요소, 마칭큐브알고리즘에 적합한 3차원 변절점요소로의 개발이 가능할 것으로 예상되며, 본 논문에서는 3차원 변절점요소를 이용한 복잡한 요소망 생성에 대한 예제로 대퇴골의 요소망 생성을 보였다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.52-55
• /
• 2010

본 논문에서는 일반유한요소법(Generalized Finite Element Method)를 이용하여 응력확대계수를 계산하는 방법을 소개한다. 기존의 유한요소법을 사용하여 응력확대계수를 계산하기위해서는 J-integral 방법 등을 이용한 후처리 과정이 필수적으로 요구된다. 뿐만 아니라 균열선단 근방에서의 응력을 기술하기 위해서는 세밀한 요소망(mesh)이 요구된다. 후처리 과정과 균열선단 근방에서의 요소망은 수치적 오류를 발생시키고 이는 정확한 응력확대계수를 얻는데 어려움을 준다. 일반유한요소법은 근사함수를 요소망의 영향 없이 추가해서 사용할 수 있는 장점을 가지고 있지만, 활용성 측면에서 기존의 유한요소법보다 복잡하여 실용성이 떨어진다. 본 논문에서는 일반유한요소법의 장점을 충분히 살려 균열선단근방에서는 응력을 모델링하여 근사함수로 사용하고 균열선단에서 거리가 먼 곳은 기존의 유한요소를 써서 계산을 하였다. 특별한 후처리 과정(Post processing) 없이 비교적 정확한 응력확대계수를 손쉽게 얻을 수 있다. 일반유한요소법을 이용한 제시된 방법론이 타당함을 수치 예제를 통하여 확인하였다.

• Journal of Korea Water Resources Association
• /
• v.52 no.2
• /
• pp.149-161
• /
• 2019

In this study, the effect of the fluid flow characteristics on the length distribution of the fracture elements composing the fracture network is analyzed numerically using the 3D fracture crack network model. The truncated power-law distribution is applied to generate the length distribution of the fracture elements and the simulations of fluid flow are carried out with the exponent ${\beta}_l$ from 1.0 to 6.0. As a result of simulations, when the exponent ${\beta}_l$ increases, the length distribution of the fracture elements gradually decreases, and the connectivity between the fracture elements affecting the permeability of the fracture network becomes weak. When we analyzed the distributions of flow rate calculated at each fracture element with the exponent ${\beta}_l$, the mean flow rate at ${\beta}_l=1.0$ was estimated to be about 447 times larger than that at ${\beta}_l=6.0$ and for the flow calculated at the outflow boundary of the fracture network, the case of ${\beta}_l=1.0$ was estimated to be 6,440 times larger than that of ${\beta}_l=6.0$.

• Journal of the Korea Concrete Institute
• /
• v.13 no.6
• /
• pp.637-644
• /
• 2001

A regularization method based on the Duvaut-Lions viscoplastic scheme for plastic-damage and continuum damage models, which provides mesh-independent and well-posed solutions in nonlinear earthquake analysis of concrete dams, is presented. A plastic-damage model regularized using the proposed rate-dependent viscosity method and its original rate-independent version are used for the earthquake damage analysis of a concrete dam to analyze the effect of the regualarization and mesh. The computational analysis shows that the regularized plastic-damage model gives well-posed solutions regardless mesh size and arrangement, while the rate-independent counterpart produces mesh-dependent ill-posed results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.335-341
• /
• 2009

In this paper, review of developed MLS-based finite elements and a proposal for their applications are described. The shape functions and their derivatives of MLS-based finite elements are constructed using Moving-Least Square approximation. In MLS-based finite element, using the adequate influence domain of weight function used in MLS approximation, kronecker delta condition could be satisfied at the element boundary. Moreover, because of the characteristics of MLS approximation, we could easily add extra nodes at an arbitrary position in MLS-based finite elements. For these reasons, until now, several variable-node elements(2D variable element for linear case and quadratic case and 3D variable-node elements) and finite crack elements are developed using MLS-based finite elements concept. MLS-based finite elements could be extended to 2D variable-node triangle element, 2D finite crack triangle element, variable-node shell element, finite crack shell element, and 3D polyhedron element. In this paper, we showed the feasibility of 3D polyhedron element at the case of femur meshing.

• Tunnel and Underground Space
• /
• v.21 no.2
• /
• pp.128-137
• /
• 2011

The shut-in pressure calculated in common hydrofracturing test for vertical borehole equals generally to the minimum horizontal principal stress, so it should be considered as an essential parameter for determining the in-situ stress regime around the rock mass. It shows usually an ambiguous value in pressure-time history curves, however, because of the relationship between the behavior of hydraulic fractures and the condition of remote stress regime. In this study, a series of numerical analyses have been carried out to compare several methods for determining the shut-in pressure during hydrofracturing. The hydraulic-mechanical coupling has been applied to numerical analysis for simulating the fracture propagation by hydraulic pressure, and the different discontinuity geometry has been considered in numerical models to examine the effect of numerical element shape on fracture propagation pattern. From the numerical simulations with the four different discontinuity geometries, it was revealed that the shut-in pressure obtained from graphical methods rather than statistical method was relatively small. Consequently a care should be taken in selecting a method for determining the shut-in pressure when a stress anomaly around borehole and a fracture propagation with complicate mechanism are considered.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.293-300
• /
• 2012

In this study, the delamination analysis has been implemented to investigate the initiation and propagation of crack in composite laminates composed of orthotropic materials. A simple modeling was achieved by moving nodal technique without re-meshing work when crack propagation occurred. This paper aims at achieving two specific objectives. The first is to suggest a very simple modeling scheme compared with those applied to conventional h-FEM based models. To verify the performance of the proposed model, analysis of double cantilever beams with composite materials was implemented and then the results were compared with reference values in literatures. The second one is to investigate the behavior of interior delamination problems using the proposed model. To complete these objectives, the full-discrete-layer model based on Lobatto shape functions was considered and energy release rates were calculated using three-dimensional VCCT(virtual crack closure technique) based on linear elastic fracture mechanics.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.129-137
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.378-386
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.