• 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 Aerospace System Engineering
• /
• v.17 no.2
• /
• pp.1-8
• /
• 2023

Power transmission shafts in rotary wing aircraft use a hollow shaft to reduce weight. We can apply linear elastic fracture mechanics to predict crack propagation behavior. This paper predicted crack growth life of a hollow shaft with a circumferential through-type crack by finite element analysis. A 2D finite element model was created by applying a torsion and forming elements considering cracks. We defined the initial crack length and performed the finite element analysis by increasing the crack length to derive stress intensity factor at crack tips. We defined the length just prior to the stress intensity factor exceeding the fracture toughness as the crack limit length. We calculated the crack limit length using a handbook and numerically integrated the crack growth rate equation to derive growth life of each crack. The growth life of each crack was compared to verify the proposed finite element analysis method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.7
• /
• pp.713-721
• /
• 2012

The finite element simulations of fatigue crack growth are carried out. Using only the mechanical properties usually obtained from the tensile test as input data, we attempted to predict the fatigue crack growth behavior. The critical crack opening displacement is determined by monitoring the change in displacements at the node close to the crack tip. Crack growth is simulated by debonding the crack tip node. The exponent in the Paris law was determined and compared to the published exponent. Plotting with respect to the effective stress intensity factor range yielded more consistent results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.537-549
• /
• 1999

선형탄성 파괴해석은 균열을 갖는 변형률 경화재료의 파괴거동을 예측하는데 불충분하기 때문에 최근에는 균열 선단 부에서 대규모 소성 역을 갖는 균열 체에 적용할 수 있는 많은 파괴역학개념이 제안되고 있다. 따라서, 본 연구에서는 대규모항복 조건하의 연성파괴를 보이는 평판을 정확하게 해석할 수 있는 새로운 유한요소모델을 제시하고자 한다. 균열 선단 부의 응력 장을 정의하는데 가장 지배적인 파괴매개변수인 J-적분 값과 소성 역의 크기 및 형상을 J-적분법과 등가영역적분법을 통해 파괴거동을 설명할 수 있도록 증분소성이론에 기초를 둔 p-version 유한요소해석이 채택되었다. 제안된 유한요소모델에 의한 수치해석결과는 이론 해와 h-version 유한요소해석과 비교되었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• Magazine of the Korea Concrete Institute
• /
• v.9 no.5
• /
• pp.157-164
• /
• 1997

일반적으로 콘크리트와 철근간의 경계면을 나타내는 유한요소법에는 균열의 부근에서 발생하는 부착열화 현상을 고려하지 않고 있다. 이것은 균열 부근에서 과도한 부착을 초래하고 , 국소 변형과 균열의 진전에도 영향을 준다. 본 연구에서는 철근콘크리트 구조물의 균열부근에서 일어나는 부착거동의 변화를 고려한 비선형 부착응력-미끄럼 모델을 제안하였다. 철근과 콘크리트간의 경계면에는 링크요소를 이용하였고, 링크의 특성은 철근을 가로지르는 균열의 상태에 따라 변하도록 조정하였다. 균열의 형성상태를 정량화하고, 부착거동을 두 포락선 1) 균열로부터 충분히 떨어진 위치에서의 부착상태를 모델링한 외연포락선, 2)횡균열면에 있어서의 부착상태를 모델링한 내연포락선의 사이에 변이시키기 위하여 비국소적 손상도 개념을 도입하였다. 이 방법의 유효성을 알아보기 위하여 편재하중을 받는 T형 교각의 실험 및 해석결과를 제시하였다. 제안된 모델의 결과를 실험결과와 비교하여 본 모델의 유용성을 검증하였다.

• Magazine of the Korea Concrete Institute
• /
• v.9 no.5
• /
• pp.189-196
• /
• 1997

본 연구에서는 유한요소법을 이용한 콘크리트 인장균열의 발생 및 전파양상을 모형화하기 위하여 분산균열(smeared crack) 모델의 일종인 2차원 균질화된 균열(homogenized crack)모형을 제안하였다. 제안한 모형은 인장균열면을 따라 속도 불연속계(velocity discontinuity)를 도입하고 평형방정식 빛 적합방정식을 이용하여 인장균열을 포함한 콘크리트의 구성방정식을 유도할 수 있으며 인장균열이 소성연화거동을 하는 경우. 유한요소망내에서 객관성이 유지될 수 있음을 밝히기 위하여 1차원 영역내에서 엄밀해를 유도하였다. 제안한 모형을 이용한 1차원 또는 2차원 유한요소 해석은 기존의 노치를 포함한 콘크리트 시편에 대한 실험과 상응하는 결과를 보였을 뿐만 아니라 유한요소망의 객관성이 유지되고 있음을 밝혔다. 제안한 모형은 큰 어려움없이 3차원 영역으로 확대할 수 있을 것이며 이에 대한 추가적인 연구가 계속될 것이다.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.38-45
• /
• 1998

It is a tedious and excessive time consuming process to model the local area of crack tip part of structures in calculation of stress intensity factors by FEM. So, in this paper, the hybrid method of FEM and BEM approach was formulated to overcome this type of problems. The multi-domained BEM was adopted to simplify the modelling process of complex geometry and singularity characteristics of crack tip part and the ordinary FEM modelling was used in the rest part. The example calculations shows very good results compared with analytic solutions and other numerical method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.1
• /
• pp.145-153
• /
• 1996

The fracture process zone in concrete is a region ahead of a traction-free crack, in which two major mechanisms, microcracking and bridging, play important roles. The toughness due to bridging is dominant compared to toughness induced by microcracking, so that the bridging is dominani: mechanism governing the fracture process of concrete. Fracture mechanics does work for concrete provided that the fracture process zone is being considered, so that the development of model for the fracture process zone is most important to describe fracture phenomena in concrete. In this paper the bridging zone, which is a part of extended rnacrocrack with stresses transmitted by aggregates in concrete, is modelled by a Dugdale-Barenblatt type model with linear tension-softening curve. Two finite element techniques are shown for the analysis of progressive cracking in concrete based on the discrete crack approach: one with crack element, the other without crack element. The advantage of the technique with crack element is that it dees not need to update the mesh topology to follow the progressive cracking. Numerical results by the techniques are demonstrated.