• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1992.10a
• /
• pp.90-96
• /
• 1992

경계면 요소를 이용하여 철근콘크리트 구조물의 접촉면 문제를 유한 요소법으로 해석하는 기법에 대하여 연구한다. 본 연구에서는 경계면 요소의 수치해석의 이론과정을 전개하고, 실험 관찰된 부착 시험체에 적용하여 이형철근과 콘크리트 부착기구의 접촉면을 해석한다. 경계면은 특별한 연결요소를 이용하여 재현하며 Mohr-Coulomb의 마찰 이론을 응응한다. 해석의 주요점으로 하중상태에 따라 변화되는 경계면의 접촉상태, 즉 고정(stick), 미끄러짐(slide), 분리(separation)를 묘사하여 경계면 재료의 비선형 거동을 관찰한다. 부착모델의 해석결과는 실험실의 결과와 대체로 일치되며 따라서 철근콘크리트 접촉면의 응력해석을 위해 경계면 요소가 활용될 수 있음을 보여준다.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2006.03a
• /
• pp.600-606
• /
• 2006

실제 지반은 경계가 없는 무한상태로 존재하기 때문에 지반구조물의 동적거동을 유한요소법을 이용하여 해석할 시 모델의 영역을 성립하는 것은 특별한 고려가 필요하다. 유한요소법에서의 동적해석은 파동의 전달을 포함하기 때문에 모델의 경계에서 인공적인 경계조건이 필요하다. 인공적인 경계 조건은 유한요소내의 지반상태를 무한상태로 변형시킬 수 있어야 하며, 경계에 도달하는 응력 파동을 모델내로 반사시키지 않고 흡수 할 수 있어야 한다. 본 논문에서는 간단한 점 탄성 반무한 불연속 요소를 이용하여 지반구조물의 동적해석을 수행하는 방법을 보여준다. 반무한 요소의 실행은 OpenSees라는 유한요소 해석프로그램을 이용하여 수행되었으며, 예를 통하여 불연속 요소가 경계에 도달하는 응력 파동을 충분히 흡수하여 유한요소 모델을 반무한 상태로 전환 시킬 수 있다는 것을 보여준다.

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

For the propagation of elastic waves in unbounded domains, absorbing boundary conditions at the fictitious numerical boundaries have been proposed. In this paper we focus on both first and second order paraxial boundary conditions(PBCs) in the framework of variational approximations which are based on paraxial approximations of the scalar and elastic wave equations. We propose a penalty function method for the treatment of PBCs and apply these into finite element analysis. The numerical verification of the efficiency is carried out through comparing PBCs with Lysmer-Kuhlemeyer's boundary conditions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.512-515
• /
• 2011

터널 해석은 주로 지표침하와 터널 라이닝 내 단면력 산정에 초점이 맞춰지며 이는 시공단계를 고려한 3차원 수치해석 모델을 이용해 결정할 수 있다. 수치해석 시 shield는 응력 제어, shell element로 모델링하는 방법 등으로 모사될 수 있다. 한편 변위 제어를 통한 쉴드 모사는 shield를 적절한 경계조건으로 처리함으로서, 다른 shield 모사 방법에 비해 모델링 작업을 간소화하고 해석의 효율성을 높일 수 있다. 본 연구에서는 변위 제어에 의한 shield 모사를 위한 적정 경계조건을 제안한다. 이를 위해 시공단계가 고려된 유한요소해석을 사용하여, 쉴드 및 굴착면에서의 경계조건 변화와 이에 따른 지표침하 관측 수행하였다. 제안된 shield 변위 제어로부터 얻어진 해석결과를 이론적인 해와 비교함으로서, 제시된 shield 모델링 방법의 적정성과 지반 거동 변화를 평가하고자 한다. 해석 결과는 지반 모델의 지표침하를 기준으로 관찰되었으며, 변위제어에 의한 결과와 요소에 의한 모델링 결과가 유사하게 얻어짐을 보여준다. 또한 변위제어의 쉴드 모사에서 회전 구속보다 변위 구속 조건에 지배적으로 영향을 받음을 확인하였다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.19-24
• /
• 1995

단순한 형상의 소음기는 평면파이론에 의해 비교적 간단하게 음향성능을 해석적으로 구할 수 있다. 그러나 소음기의 형상이 복잡해지거나 해석하고자 하는 주파수의 범위가 평면파의 차단주파수 이상이 될 경우 소음기 내부의 음장이 평면파에서 벗어나게 되어 평면파 이론에 의한 해석은 실제와 상당한 오차가 발생하게 되므로 음장에 대한 3차원 해석이 필요하다. 이론적으로 3차원 문제를 해석할 수 있는 경우는 형상이 극히 단순한 경우에 국한되므로 유한요소법(FEM), 경계요소법(BEM)과 같은 수치해석적인 방법이 이용되고 있다. 경계요소법은 적분 커넬(kernel)의 특이성(singularity) 문제가 있지만 대상 영역의 경계면만을 이산화함으로써 모델링에 소요되는 시간과 노력을 절약할 수 있으므로 음향문제 해석에 있어서 효율적인 방법이라고 할 수 있다. 본 연구의 목적은 3차원 경계요소법 프로그램을 개발하고 평면파이론에 의한 해석이 어려운 여러가지 형태의 소음기에 대한 음향성능을 예측하고 실험으로 검증하는것이다. 특히, 단일영역으로 해석이 불가능한 다공형 소음기에 영역분할법을 적용하여 계산하고 결과를 검토하였다.

• Journal of the Korean Geosynthetics Society
• /
• v.20 no.3
• /
• pp.11-20
• /
• 2021

In this paper, the approximately coupled method of finite element method and boundary element method to obtain efficient and accurate analysis results is proposed for a two-dimensional elasto-static problem with a geometrically abruptly changing part. As the finite element of a two-dimensional problem, three-node and four-node plane stress element is applied, and as the boundary element of a two-dimensional problem, three-node boundary element is applied. In the modeling stage, firstly, an entire analysis target object is modeled as finite elements, and then a geometrically abruptly changing part is modeled as boundary elements. The boundary element is defined using the nodes defined for modeling finite elements. In the analysis stage, finite element analysis is firstly performed on a entire analysis target object, and boundary element analysis is automatically performed afterwards. As for the boundary conditions at boundary element analysis, displacement conditions and stress conditions, which are the results of finite element analysis, are applied. As a numerical example, the analysis results for a two-dimensional elasto-static problem, a plate with a crack, are presented and investigated.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.127-134
• /
• 1996

This paper is about a progressive fracture analysis of concrete by boundary element method. From both displacement boundary integral equation and traction boundary integral equation of solids with cracks, a boundary integral equation for crack problem is derived. For the analysis of progressive fracture of concrete, fracture process zone is modelled based on Dugdale-Barenblatt model with linear tension-softening curve. By using the boundary element modeling, the progressive fractures of concrete beam and compact-tension specimens with various loading conditions are analyzed and compared with experiments. The analysis results show that the technique in this paper can predict the maximum strength and the nonlinear behavior of concrete including post-peak behavior.

• 전기의세계
• /
• v.39 no.3
• /
• pp.47-54
• /
• 1990

유한요소법과 경계요소법의 합성으로 전자계 해석을 하는 기법은 각 방법의 장점을 수용하여 경계가 없는 무한영역의 전자장을 분석하는 기법으로서 어떤 복잡하고 어려운 기하학적 구조의 문제도, 비선형이나 비균질성 재질의 문제도 쉽게 formulation이 가능하여 용이하게 해석할 수 있지만 전체 System matrix방정식이 비대칭이며 부분적인 full matrix를 형성하여 계산시간이 길어 진다는 단점도 있다. 적용예에서 보여 준 것과 같이 합성요소법은 그 해가 실제에 근사한 값을 가질수 있다고 생각되며, 계산시간을 단축시키기 위하여 직접법이나 반복법을 사용한 새로운 해법들이 도입되고 있다. 최근에는 system전체 node의 순서를 고려한 NDRA(Nested Dissection Reordering Algorithm)이 도입되고 있고, System matrix자체를 유한 요소법의 형태로 유지시키며 풀수 있는 방법으로 알려진 Absorbin 경계조건을 사용하여 전자파에 대한 해석을 하고 있다. 유한 및 경계요소 합성법은 초고압 옥외용 전력기기의 전자장 해석과 설계, 레이다나 안테나 등의 전자파 해석문제, 초전도 응용, 전력기기의 전자장해석과 설계, 우주공간에서의 전력전송문제 등을 쉽게 model화하여 적용할 수 있을 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.398-402
• /
• 2007

본 연구에서는 삼각 요소망을 구성하는데 사용되는 기존의 Delaunay Triangulation 기법 외에 Advanced Front Technique 알고리즘을 적용하여 하천의 경계가 중요시되는 해석 모델에 적합한 삼각 요소망을 생성하고자 한다. Advanced Front Technique 알고리즘은 외부 경계를 따라 삼각 요소를 형성하기 시작하여 전체 해석 영역으로 요소를 채워가는 방법이다. 이 방법은 위에서 언급한 Delaunay Triangulation 기법의 단점인 해석 영역의 경계에서 요소의 형질이 좋지 못한 문제를 극복할 수 있을 뿐만 아니라, 미리 절점을 생성하지 않으므로 메모리를 절약하는 효과도 얻게 된다. Advanced Front Technique 알고리즘은 외부 경계의 노드(i)와 노드(i+1)을 Front로 하고, 이등변삼각형을 이루기 위해 해석 영역 내부의 적절한 위치에 새로운 절점(C)을 생성한다. 이렇게 Front와 새로운 절점(C)으로 이루어진 이등변 삼각형 요소가 생성된 후에, Front는 다음 Front로 전진하면서 점차 해석 영역의 내부로 삼각형 요소를 채우게 된다. 이와 같은 방법에 의해서 해석영역의 경계에서는 비교적 이등변 삼각형에 가까운 삼각 요소가 생성된다.