• 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차원 변절점요소를 이용한 복잡한 요소망 생성에 대한 예제로 대퇴골의 요소망 생성을 보였다.

• Tunnel and Underground Space
• /
• v.22 no.2
• /
• pp.93-105
• /
• 2012

This paper compares the two- and three-dimensional (2D and 3D) approaches for the numerical determination of the equivalent mechanical properties of fractured rock masses. Both orthogonally-fractured model and discrete fracture networks (DFN) were used for the geometry and 2D models were cut in various directions from 3D model to compare their mechanical properties. Geological data were loosely based on the data available from Sellafield, UK. Analytical method based on compliance tensor transformation was used for investigation in orthogonally fractured rock and numerical experiments were conducted on fractured rock mass with DFN geometry. It is shown that 2D approach always overestimates the elastic modulus of fractured rock masses by a factor of up to around two because fractures are assumed to be perpendicular to the model plane in 2D problems. Poisson ratios tend to have larger values in 2D analysis while there is opposite trend in some sections. The study quantitatively demonstrates the limitation of the 2D approach that uses the simplified model from true 3D geometry.

• 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.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 2004.11a
• /
• pp.297-302
• /
• 2004

• 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 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$.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.326-326
• /
• 2016

The fault can be defined, in a geological context, as a rupture plane showing a significant displacement generated in the case that the local tectonic stress exceeds a threshold of rupture along a particular plane in a rock mass. The hydrogeological properties of this fault can be varied with the spatial distribution and the connectivity of void spaces in a fault. When the formation of fault includes the process of the creation and the destruction of void spaces, a complex relation between the displacement along the fault and the variation of void spaces. In this study, the variation of flow with the geometrical characteristics of the fault is simulated and analyzed by using the three-dimensional discrete fracture network model. Three different geometrical characteristics of the faults are considered in this study: 1) simple hydraulic conductive plane, 2) damaged zone, and 3) relay structure of faults.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.247-253
• /
• 2007

Determination of crack depth in filed using the self-calibrating surface wane transmission measurement and the cutting frequency in the transmission function (TRF) is very difficult due to variations of the measurement conditions. In this study, it is proposed to use the measured full TRF as a feature for crack depth assessment. A Principal component analysis (PCA) is employed to generate a basis of the measured TRFs for various crack cases. The measured TRFs are represented by their projections onto the most significant principal components. Then artificial neural networks (NNs) using the PCA-compressed TRFs is applied to assess the crack in concrete. Experimental study is carried out for five different crack cases to investigate the effectiveness of the proposed method. Results reveal that the proposed method can be effectively used for the crack depth assessment of concrete structures.

• 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.

• Tunnel and Underground Space
• /
• v.11 no.1
• /
• pp.42-58
• /
• 2001

In order to protect the environment from waste disposal activities, the prediction of the flux and flow paths of the contaminants from underground facilities should be assessed as accurately as possible. Especially, the prediction of the pathways and travel times of the nuclides from high level radioactive wastes in a deep repository to biosphere is one of the primary tasks for assessing the ultimate safety and performance of the repository. Since the contaminants are mainly transported with groundwater along the discontinuities developed within rock mass, the characteristics of groundwater flow through discontinuities is important for the prediction of contaminant fates as well as safety assessment of a repository. In this study, the actual fracture network could be effectively generated based on in situ data by separating geometric parameter and hydraulic parameter. The calculated anisotropic hydraulic conductivity was applied to a 3D porous medium model to calculate the path flow and travel time of the large studied area with the consideration of the complex topology in the area. Using the model, the pathways and travel times for groundwater were analyzed. From this study, it was concluded that the suggested techniques and procedures for predicting the pathways and travel times of groundwater from underground facilities to biosphere is acceptable and those can be applied to the safety assessment of a repository for radioactive wastes.