• 한국전산구조공학회논문집
• /
• 제23권5호
• /
• pp.475-482
• /
• 2010

본 논문에서는 대규모 자유도 시스템의 병렬처리를 위하여 2단계로 이루어진 영역분할법(Domain Decomposition Method) 기반의 병렬 알고리즘을 제안하였다. 분할된 영역의 내부 및 외부 경계를 상위영역문제로 정의하고 국부영역문제는 변위 경계조건이 모두 주어지는 분할영역에서의 Dirichlet 문제로 구성한다. 상위영역에서는 전체 상위영역에 대한 강성 행렬의 어셈블이 필요없는 반복법을 통하여 변위를 구하고, 이를 바탕으로 국부영역에서 Multi-Frontal Sparse Solver (MFSS)를 이용하여 변위를 계산한다. 상위영역문제의 연산에서 프로세서 간의 데이터 교환을 최소화하여 계산효율을 유지하며, 동시에 해석 가능한 자유도를 증대시키는 병렬 PCG(Preconditioned Conjugate Gradient)법 기반의 알고리즘을 개발하였다. 제안된 알고리즘을 적용하여 수치해석을 수행한 결과, 프로세서 수가 증가할수록 계산성능의 손실없이 해석 가능한 자유도가 비례하여 증가하는 선형 확장성을 관찰할 수 있었으며, 대규모 자유도 문제에 효과적으로 사용 가능함을 확인하였다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.259-264
• /
• 2007

In this paper, to analyze the initial valued non-homogeneous elastic half space by the scaled boundary analysis, the infinite element approach was introduced. The free surface of the initial valued non-homogeneous elastic half space was mode1ed as a circumferential direction of boundary scaled boundary coordinate. The infinite element was used to represent the infinite length of the free surface. The initial value of material property(elastic modulus) was considered by the combination of the position of the sealing center and the power function of the radial direction. By use of the mapping type infinite element, the consistent e1ements formulation could be available. The performance and the feasibility of proposed approach are examined by two numerical examples.

• 수산해양기술연구
• /
• 제18권1호
• /
• pp.47-53
• /
• 1982

본 논문은 비균일대칭성 열Flux인 수직사각 Duct내의 층류조합대류 열전달 효과를 해석하기 위하여 그 유동의 특성 지배 방정식 및 비균일 열Flux의 경계조건을 무차원화 시켜 이를 Galerkin's 방법에 의해 유한요소식으로 정식화하고 이에 대하여 R 하(a) 수 및 압력구배 변수에 대해서 Duct 내의 온도분포, 속도분포 및 Nusselt 수의 값을 계산하였고 온도분포를 열 Flux가 일정 및 없는 경우와 비교하였으며 또 닥트내의 열전달 특성을 R 하(a) 수, 응력구배변수 및 Corner에 따른 변호경향을 조사하였다. 그 결과 1. 본 해석의 경계벽 온도분포 계산치와 유효자료들과의 비교에서 열 Flux가 일정 또는 없는 경우는 그 값이 일치하였다. 2. 닥트내의 온도분포와 Nusselt수의 값은 R 하(a) 수 및 압력구배 변수에 비례하여 증감하였다. 3. Nusselt수는 Corner에서 유속지연에 의한 온도분포의 특성 때문에 그 값이 감소하였으며 최대치는 0.7부근이었다

• 자원환경지질
• /
• 제36권6호
• /
• pp.491-499
• /
• 2003

오키나와트러프는 인장력에 의해 형성된 열개 분지이다. 한반도 남부해역에 위치하는 북부 오키나와트러프의 북서 주변부에서 취득한 다중 채널 탄성파 자료 해석에 의하면 오키나와트러프의 지하 지질 구조는 전형적인 열개 분지에서 나타나는 경사진 단층 블록, 단층 블록들 사이에 발달한 반지구, 그리고 반지구의 경계부에 발달한 점완단층으로 특징 지워진다. 오키나와트러프에 퇴적된 퇴적층은 퇴적 시기에 따라 선열개, 동시열개 및 후열개 퇴적층으로 나뉘며 이들은 서로 다른 탄성파 반사 특성을 나타낸다. 선열개 퇴적층은 경사진 단층 블록의 경사 방향에 평행하게 발달한다. 동시열개 퇴적층은 주로 경사진 단층 블록 사이에 발달한 반지구에 퇴적되며 경사진 단층 블록을 향하여 반사면들이 경사를 이루며 발산(divergence)하는 특징을 나타낸다. 이는 동시열개 퇴적층이 단층 블록의 기울어짐(tilting)과 동시에 퇴적이 이루어졌음을 시사한다. 후열개 퇴적층은 반사면의 연속성이 양호하며 평행하게 발달한다. 열개 분지에서 반지구의 폭과 분지 경계 단층의 낙차 변위율은 밀접한 관계가 있다. 낙차 변위율은 열개 작용이 가장 활발한 시기에 최대이며 반지구의 폭은 분지 경계 단층의 낙차 변위율에 비례한다.

• 한국전산구조공학회논문집
• /
• 제33권3호
• /
• pp.201-207
• /
• 2020

본 연구에서는 다공성 보와 논로컬 매개변수 사이의 관계에 대한 유한요소해석을 수행한다. 논로컬 매개변수는 다공성 보의 결함을 표현하는 변수들로 정의하여, 하중조건 및 경계조건에 대한 수치모사를 통해 계산한다. 다공성 보와는 반대 개념의 결함을 가지는 보에 대한 해석도 수행하였다. 이러한 보들의 거동은 논로컬 매개변수의 항으로 표현하였으며, 이 매개변수는 구멍의 지름의 제곱 그리고 원기둥 지름의 세제곱에 비례하는 것을 확인하였다. 특히 작은 원기둥을 가지는 보에 축 하중을 가하는 경우, 예상과는 다르게 3차원 유한요소 해석 결과와 2차원 평면응력 해석 결과는 다름을 알 수 있었다.

• 대한조선학회논문집
• /
• 제28권1호
• /
• pp.38-52
• /
• 1991

자유표면하에서 움직이는 임의의 형상의 3차원 물체로 인한 비선형 유체력을 준 Lagrangian 방법을 사용하여 해석하였다. 경계치 문제는 경계 적분 방법(Boundary Integral Method)을 이용하여 해결하였으며, 물체와 자유 표면의 형상은 곡면 Panel로 표현하였다. 이들 표면은 hi-cubic B-spline 방법을 사용하여 유한개의 작은 표면 요소로 나뉘어지게 되며, 또한 \phi와 (equation omitted) 표면 요소상에서 bi-linear하다고 가정한다. 특이점에 의한 유기 포텐시얼의 계산시 1/R에 비례하는 부분은 제거하고 해석적으로 처리하였다. 물체로부터 멀리 떨어진 곳에서의 유체 유통은 좌표계의 원점에 위치한 Dipole로 표현하였으며, Time Stepping시 Runge-Kutta의 4차 방법을 사용하였다. 3차원인 경우에 대한 적분 방정식과 포텐시얼의 시간 미분간에 대한 경계 조건이 유도되었으며, 이러한 식들을 사용하여 자유표면의 형상과 물체의 운동을 동시에 계산하였다. 대진폭을 가지고 규칙적으로 진동하는 물체에 작용하는 힘과 이때의 자유 표면 형상을 계산하고 기 발표된 자료와 비교하여 보았으며, 자유표면 근처에서 운동하는 물체에 작용하는 비선형 효과를 관찰하였다.

• 수산해양기술연구
• /
• 제38권2호
• /
• pp.110-118
• /
• 2002

본 논문에서는 무요소 이론을 정식화하였고 이를 이용한 1차원 및 2차원 EFG 프로그램을 Visual Basic과 C언어를 이용하여 작성해 보았다. 그리고 각각의 EFG 수치해석의 예를 작성된 프로그램을 이용하여 해를 구하였다. 해석결과는 다른 문헌의 결과와 일치하였으며 해석결과에서 나타나듯이 무요소 해의 정도는 영향영역의 비례축소인자 dmax와 가중함수의 종류, 절점 배치형태에 의해 좌우된다는 사실을 알 수 있었다. 특히 1, 2차원 EFG 해석결과에서 가장 최적의 해를 보이며 정해(exact solution)에 가장 근접한 조건은 dmax = 2 이고 가중함수가 3차 Spline형일 때로 나타났으며 유한요소법과 마찬가지로 절점의 수가 많을 수록 그리고 절점을 균일하게 배치할수록 높은 정도를 나타내는 것을 알 수 있었다. 특히 2차원의 경우 3차 Spline형 이외의 다른 가중함수를 사용할 경우에 상당히 큰 오차를 나타내는 점은 1차원 EFG 해석의 결과와는 다른 점이었지만 그 외 대부분 같은 결과를 나타내었다. 1차원에서 절점을 임의로 배치한 경우는 비교적 균일하게 배치한 경우가 해에 근접하는 형태를 나타내었으며 절점 간격이 상대적으로 적은 곳에서 큰 오차를 나타내었다. 그리고 절점을 임의로 선택할 때 변위가 모두 ‘0’의 값을 가지는 경우를 볼 수 있는데, 이는 화면상의 좌표계산에서 생긴 미소한 오차가 절점들에 의해 반복됨으로서 발생하는 것으로 보인다. 또한 탄성계수 값이 클 경우 dmax 에서 계산이 제대로 수행되지 못하는 경우가 있는데, 이는 수치가 double형의 크기를 초과하기 때문인 것으로 보인다. 결과에서 나타나듯이 무요소법에서 적당한 가중함수와 비례축소 인자를 사용하면 정해에 가까운 우수한 해를 얻을 수 있다는 것을 알 수 있다. 비록 프로그래밍 과정이나 이론의 정식화가 유한요소법에 비해 상당히 어려운 점은 있으나 무요소법은 요소의 정보를 필요치 않으므로 사용자 입장에서는 매우 편리할 것이다. 앞으로 경계조건을 효과적으로 만족시키는 문제를 해결하고 효과적인 알고리즘이 개발된다면 실용적으로 유한요소법을 대신할 수 있는 좋은 대안이 될 수 있을 것이라 생각된다.

• 한국전산구조공학회논문집
• /
• 제29권6호
• /
• pp.583-590
• /
• 2016

본 논문은 강형식 기반의 MLS 차분법에 Rayleigh 감쇠효과를 적용한 동적균열진전 해석기법을 제시한다. Rayleigh 감쇠 효과가 반영된 동적 평형방정식과 구성방정식을 도출하고, MLS 미분근사식을 이용하여 지배방정식들을 이산화하였다. 평형방정식뿐만 아니라 구성방정식에서도 감쇠효과를 적절하게 고려하여 기존의 무요소 강정식화 기법에서 고려하지 못했던 비례감쇠 알고리즘을 구현하였다. 시간관련 항을 포함한 동적 평형방정식은 중앙차분법(central difference method)을 이용하여 시간적분 하였고, 속도에 대한 차분식을 lagging시켜 이산화 방정식을 간소화시켰다. 균열의 기하학적 특성은 표면력 '0'인 자연경계 조건을 균열면에 놓인 절점들에 부과하여 묘사하였으며, 균열성장으로 인해 해석단계마다 변하는 절점의 생성 및 이동 효과를 계방정식 구성에 반영하였다. 단일균열과 다중균열을 갖는 수치예제를 통해서 제안된 수치기법의 정확성을 검증하였으며, 비례감쇠 효과의 고려가 동적균열진전 해석결과에 미치는 영향을 보였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.548-552
• /
• 2005

새만금호의 완공이후 내부 개발이 완료된 시점에서 만경강 및 동진강 홍수량 유입이 담수호의 수위 및 유속등 수리특성에 어떠한 변화를 가져올 것인가를 평가하기 위한 수치실험이 내부 관리수위(DL=-0.5m, DL=-1.0m, DL=-1.5m)별로 실시되었다. 홍수파의 전파에 따른 동적거동을 이상적으로 모의하기 위해서 본 연구에서는 2차원 수심적분의 모형인 ADCIRC 모형을 새만금 하구호에 적용하였다. 상류의 유입 경계조건으로는 50년 빈도의 홍수수문곡선도를 사용하여 만경강과 동진강의 유입량이 30분 간격으로 변화되는 홍수량을 이용하였다. 하류의 경계조건으로는 신시 수문이 외해의 조위와 내부 관리수위와 연동되어 외조위가 낮아질 때만 수문이 열리는 조건을 설정하는 동적인 경계조건이 설정되었다. 홍수 계산은 4일에 걸쳐 실시되어 새만금호 내부에서 수위와 유속의 변화를 관찰하였다. 만경수계내 홍수파의 동적 거동은 시간이 경과함에 따라 홍수파가 하류로 전파되어 수위와 유속 변화가 크게 나타난다. 만경강 유입부에서는 홍수수문곡선의 시간적 변화가 민감하게 작용하는 반면, 홍수파 전파가 하류로 진행하면서 유입부에서 높아졌던 수위가 새만금호 내부에서는 점차 낮아지면서 안정화 되는 것을 보여주고 있다. 그러나 수문이 위치한 지점에서는 유입 홍수에 의한 영향은 거의 나타나지 않고, 신시수문의 개방에 따라 동적으로 변동되는 것이 지배적이다. 유속 모의 결과 상류 수로부에서는 $1.8m/s\~0.5m/s$로 나타나지만 새만금호에 진입하면서 $0.5m/s\~0.2m/s$의 낮은 유속으로 변한다. 신시 갑문 인근에서는 수문을 닫을 때 일시적인 충격파가 내부로 전파되면서 와동현상을 보이는 등 유속과 수위 변화에 있어 동적으로 복잡한 거동을 보이는 것으로 해석 된다.한 영향을 미치는 인자에 대한 이론적인 분석을 수행하고, 배수갑문 개방에 의한 수질개선효과를 최대화하기 위한 환경관리 방안 제시에 중점을 두어 수행하였다.ncy), 환경성(environmental feasibility) 등을 정성적으로(qualitatively) 파악하여 실현가능한 대안을 선정하였다. 이렇게 선정된 대안들은 중유역별로 검토하여 효과가 있을 것으로 판단되는 대안들을 제시하는 예비타당성(Prefeasibility) 계획을 수립하였다. 이렇게 제시된 계획은 향후 과학적인 분석(세부평가방법)을 통해 대안을 평가하고 구체적인 타당성(feasibility) 계획을 수립하는데 토대가 될 것이다.{0.11R(mm)}(r^2=0.69)$로 나타났다. 이는 토양의 투수특성에 따라 강우량 증가에 비례하여 점증하는 침투수와 구분되는 현상이었다. 경사와 토양이 같은 조건에서 나지의 경우 역시 $Ro_{B10}(mm)=20.3e^{0.08R(mm)(r^2=0.84)$로 지수적으로 증가하는 경향을 나타내었다. 유거수량은 토성별로 양토를 1.0으로 기준할 때 사양토가 0.86으로 가장 작았고, 식양토 1.09, 식토 1.15로 평가되어 침투수에 비해 토성별 차이가 크게 나타났다. 이는 토성이 세립질일 수록 유거수의 저항이 작기 때문으로 생각된다. 경사에 따라서는 경사도가 증가할수록 증가하였으며 $10\% 경사일 때를 기준으로 $Ro(mm)=Ro_{10}{\times}0.797{\times}e^{-0.021s(\%)}$로 나타났다.천성 승모판 폐쇄 부전등을 초래하는 심각한 선천성 심질환이다. 그러나 진단 즉시 직접 좌관상동맥-대동맥 이식술로 수술적 교정을 해줌으로써 좋은 성적을 기대할 수 있음을 보여주었다.특히 교사들이 중요하게 인식

• 한국음향학회지
• /
• 제18권6호
• /
• pp.31-37
• /
• 1999

본 논문은 근거리 음장 측정 자료로부터 원거리 음장 예측을 위한 기술에 대한 것이다. 음원의 음장 분포 특성은 원거리에서 측정된 자료의 해석으로 이루어지는 것이 일반적 방법이나, 음향수조 또는 무향실과 같은 제한된 공간에서는 근거리 영역에서 측정이 이루어지는 경우가 발생한다. 따라서 근거리 영역에서의 측정으로부터 원거리 음장이 예측되어야 한다. 이 경우 음원을 둘러싼 근거리 음장의 측정점수는 원거리 음장 예측치의 정확도와 자료 처리의 계산량과 상관된다. 기존 연구 결과는 최적측정점수는 음원의 kL에 비례하고 음원의 기하학적 형태 또는 지향특성에 따라 kL의 의존성이 다르게 나타난다고 되어 있으나 정확한 기준이 없다. 따라서 본 논문에서는 최적측정지점수에 대한 기준을 유도하기 위해 Helmholtz 적분식과 Green 함수를 근간으로 한 원거리 음장 예측 기술인 경계배치법(Boundary Collocation Method)을 분석하여 최적측정점수는 kL이 증가함에 따라 0.54kL로 수렴한다는 결과를 얻었다. 기존의 연구 결과 보다 최적측정점수를 1/2 정도로 줄였다.