본 논문에서는 전산점근해석기법을 사용하여 복합재료 보에 대한 경계층 해를 계산하고, ANSYS 결과와 비교 검증하였다. 경계층 해는 내부해와 순수 경계층 효과의 합으로 표현되기 때문에, 내부 및 경계층에 대한 수학적으로 엄밀한 정식화를 요구한다. 전산점근 해석기법은 수학적으로 매우 강력한 기법으로, 이러한 문제에 유용하다. 그러나 경계층과 내부 해들의 연결을 시키기 쉽지 않은데, 본 연구에서는 가상일의 원리를 통해 생브낭의 원리와 내부 및 경계층 문제를 체계적으로 분리하였다. 경계층 해는 팝코비치-패들 고유 벡터를 계산하여, 실수부와 허수부 벡터들의 선형 조합으로 표현하고, 내부 해의 워핑 함수들을 보상할 수 있도록 최소오차 자승법을 적용하였다. 계산된 해들은 2차원 유한요소 해석 결과와 비교하여 정성적일 뿐만 아니라 정량적으로도 잘 일치하는 결과를 얻었다.

The radwaste repository consists of a multi-barrier, including natural and engineered barriers. The repository’s long-term safety is ensured by using the isolation and delay functions of the multi-barrier. Among them, natural barriers are difficult to artificially improve and have a long time scale. Therefore, in order to evaluate its performance, site characteristics should be investigated for a sufficient period using various analytical methods. Natural barriers are classified into lithological and structural characteristics and investigated. Structural factors such as fractures, faults, and joints are very important in a natural barrier because they can serve as a flow path for groundwater in performance evaluation. Considering the condition that the radioactive waste repository should be located in the deep part, the drill core is an important subject that can identify deep geological properties that could not be confirmed near the surface. However, in many previous studies, a unified method has not been used to define the boundaries of structural factors. Therefore, it is necessary to derive a method suitable for site characteristics by applying and comparing the boundary definition criteria of various structural factors to boreholes. This study utilized the 1,000 m deep AH-3 and DB-2 boreholes and the 500 m deep AH-1 and YS- 1 boreholes drilled around the KURT (KAERI Underground Research Tunnel) site. Methods applied to define the brittle structure boundary include comparing background levels of fracture and fracture density, excluding sections outside the zone of influence of deformation, and confining the zone to areas of concentrated deformation. All of these methods are analyzed along scanlines from the brittle structure. Deriving a site-specific method will contribute to reducing the uncertainties that may arise when analyzing the long-term evolution of brittle structures within natural barriers.

Considering the non-linear behavior of structure and soil when evaluating a nuclear power plant's seismic safety under a beyond-design basis earthquake is essential. In order to obtain the nonlinear response of a nuclear power plant structure, a time-domain SSI analysis method that considers the nonlinearity of soil and structure and the nonlinear Soil-Structure Interaction (SSI) effect is necessary. The Boundary Reaction Method (BRM) is a time-domain SSI analysis method. The BRM can be applied effectively with a Perfectly Matched Layer (PML), which is an effective energy absorbing boundary condition. The BRM has a characteristic that the magnitude of the response in far-field soil increases as the boundary interface of the effective seismic load moves outward. In addition, the PML has poor absorption performance of low-frequency waves. For this reason, the accuracy of the low-frequency response may be degraded when analyzing the combination of the BRM and the PML. In this study, the accuracy of the analysis response was improved by adjusting the PML input parameters to improve this problem. The accuracy of the response was evaluated by using the analysis response using KIESSI-3D, a frequency domain SSI analysis program, as a reference solution. As a result of the analysis applying the optimal PML parameter, the average error rate of the acceleration response spectrum for 9 degrees of freedom of the structure was 3.40%, which was highly similar to the reference result. In addition, time-domain nonlinear SSI analysis was performed with the soil's nonlinearity to show this study's applicability. As a result of nonlinear SSI analysis, plastic deformation was concentrated in the soil around the foundation. The analysis results found that the analysis method combining BRM and PML can be effectively applied to the seismic response analysis of nuclear power plant structures.

해양경계획정에 대한 국가들의 시도는 분쟁당사국 간의 합의, 국제재판소의 사법적 판결로써 일반적인 흐름을 형성해왔다. 이것은 법의 유권적 해석에 관 한 정형화된 절차 또는 기관이 확정되어 있는 국내법과는 달리 국제법에 있어서는 ICJ, ITLOS 또는 중재법원의 결정들이 국제법의 해석에 관한 한 가장 권 위적이고 결정적인 것으로서 받아들여지고 있기 때문이다. 다만 국제재판소의 재판의 준칙이 되어야 할 유엔해양법협약(UNCLOS)의 관련 규정은 해양경계획정의 목적이 ‘형평한 해결’에 있다는 것만 밝히고 있을 뿐 그 방법이나 고려해야할 관련사정 등에 대해서는 구체적인 기준을 제시하지 않고 있기 때문에, 해양경계획정문제는 아직까지 관련당사국들의 상이한 주장에 맡겨진 채 쉽게 해결에 이르지 못하고 있다. 더욱이 해양경계획정의 각 사건은 독특한 혹은 독자적인 것으로서 ‘선례구속의 원칙(Doctrine of stare decisis)’이 인정되지 않기 때문에 이들 판례로부터 일반적인 법원칙을 도출해 내는 것은 매우 어려운 일이었다. 그럼에도 불구하고 국제재판소가 자신의 선 판례를 존중하는 입장을 취하면서 후속 사건의 판단기준으로 인용하고 있다는 점은 주목 할 만하다. 이러한 관점에서 그동안 축적되어온 국제 판례상의 해양경계획정원칙과 방법, 관련사정 등을 종합적으로 검토함으로써 해양경계획정문제에 적용할 수 있는 일반적 규칙을 도출하고 이러한 규칙에 따라 경계획정을 권고하는 방안을 검토해 나가는 것은 실익이 있다고 판단된다. 다만 유엔해양법협약의 개정이 현실적으로 매우 어려운 점을 감안하여 국제 판례의 리스테이트먼트 작업을 통한 법적 정리 및 국제적 권고 수준의 가이드라인 제정방식을 제안하고자 한다. 무엇보다도 오늘날 해양권익의 확대에 따라 국가 간의 명확하고 신속한 경계획정이 우선적으로 요청된다는 점을 상기하면서, 국가들은 구체적인 해양경계획 정방법을 고안하는데 협력하여야 한다.

In this study, we propose a new scheme of nonlinear analysis for Incheon International Airport Terminal-2 which was opened on January of 2018 for the Olympic Winter Games of PyeongChang in South Korea. The terminal was built by a single layered irregular space frame. It has hard problems for nonlinear analysis geometrically, because of a limitation of personal computer's ability by the number of rigid joints in the roof. Therefore we attempt easier approach to be chosen a center part of the roof instead of the whole structure, and to substitute the other boundary parts as springs. The scheme shows some merits for saving memory and calculation time and so on.

Boundary element solution method is introduced for radiation heat transfer problem with objects inside a 2-D enclosure, where shadow zone exists. Surfaces are assumed as diffuse gray in a transparent medium. Boundary integral and boundary element equations were derived from radiation transfer equation, and their theory is reviewed. Also the process of solution methods to implement the boundary element method is analysed and explained with consideration of shadow effects. BEM solution results are compared with two test problems and are found to be good agreements with the both analytic and ANSYS Fluent numerical solutions. Therefore the current BEM analysis for radiation heat transfer problem can be considered as verified, and their efficacy with engineering applicability is established as a result.

1차원 등엔트로피 모델과 통합된 경계층 적분법은 초음속 노즐의 설계과정에서 내열재 표면의 열전달을 예측하는데 효과 적으로 사용되고 있지만 노즐 목과 같이 2차원 효과와 경계층과 노즐 코어유동의 상호작용이 발생하는 지점에서는 경계층 외부유동 해석의 부정확성으로 해석의 정확도가 감소한다. 따라서 본 연구에서는 경계층 적분법을 이용한 열전달 예측의 정 확도를 향상시키기 위해 CFD를 이용하여 2차원 효과와 노즐 코어유동의 상호작용이 고려된 경계층 외부유동 조건을 도출 하고 이를 경계조건으로 하는 해석기법을 개발하였다. 오일러 모델과 SST k-ω 모델을 CFD로 해석하여 경계조건으로 적용 했으며 계산방법을 검증하기 위해 선행문헌의 실험노즐에 대해 해석을 수행하였다. 계산 결과 CFD를 통해 경계층 외부유동 조건을 도출한 해석에서 노즐 열전달의 정확도가 향상되는 것을 확인하였으며 특히 노즐 목 후방과 팽창부에서의 차이가 크 게 나타났다. SST k-ω모델로 도출된 계산결과는 1차원 등엔트로피 모델과 비교 시 팽창부에서 실험결과와의 오차가 16% 감소하였다. 본 연구에서 개발된 해석기법은 향후 로켓노즐의 내열설계에 유용하게 사용될 것으로 평가된다.

This paper presents a detailed procedure for a nonlinear soil-structure interaction of a seismically isolated NPP(Nuclear Power Plant) structure using the boundary reaction method (BRM). The BRM offers a two-step method as follows: (1) the calculation of boundary reaction forces in the frequency domain on an interface of linear and nonlinear regions, (2) solving the wave radiation problem subjected to the boundary reaction forces in the time domain. For the purpose of calculating the boundary reaction forces at the base of the isolator, the KIESSI-3D program is employed in this study to solve soil-foundation interaction problem subjected to vertically incident seismic waves. Wave radiation analysis is also employed, in which the nonlinear structure and the linear soil region are modeled by finite elements and energy absorbing elements on the outer model boundary using a general purpose nonlinear FE program. In this study, the MIDAS/Civil program is employed for modeling the wave radiation problem. In order to absorb the outgoing elastic waves to the unbounded soil region, spring and viscous-damper elements are used at the outer FE boundary. The BRM technique utilizing KIESSI-3D and MIDAS/Civil programs is verified using a linear soil-structure analysis problem. Finally the method is applied to nonlinear seismic analysis of a base-isolated NPP structure. The results show that BRM can effectively be applied to nonlinear soil-structure interaction problems.

이 논문에서는 경계반력법을 이용한 비선형 지반-구조물 상호작용 해석을 위해 LS-DYNA나 MIDAS/Civil 등의 유한요소해석 프로그램과 연계하는 방법을 제시하였다. 경계반력법 적용시 유한요소프로그램에서 구조물과 지반은 선형 또는 비선형 유한요소를 이용하여 모델링하였다. 유한요소의 해석모델 외부의 무한영역으로 전달되는 탄성파를 최대한 흡수하기위해 유한요소 모델의 외측에 LS-DYNA의 경우에는 PML(Perfectly Matched Layer) 요소를, MIDAS/Civil의 경우에는 점성감쇠-스프링 요소를 적용하였다. 비선형 유한요소는 구조물영역에만 적용되는 것으로 가정하였다. 이 연구에서는 입사지진파에 의한 경계반력은 KIESSI-3D 프로그램을 이용하여 계산하였다. 선형 지반-구조물 상호작용 문제에 대해 일반적인 KIESSI-3D의 해석결과와 BRM해석결과를 비교하여 제시된 방법의 효율성을 제시하였다. 또한 수치적 비교를 통해 비선형 구조에 대해 보수적인 응답을 보이는 선형 SSI문제에 대하여 얻은 경계반력이 비선형 지반-구조물 상호작용해석에 효과적으로 적용 가능함을 알 수 있었다.

본 논문은 1차원 자유경계문제 해석의 정확도 향상을 위해 이동최소제곱 차분법을 이용하여 이동경계의 위상변화를 implicit하게 추적하는 기법을 제시한다. 기존의 이동최소제곱 차분법은 이동경계의 위치를 explicit하게 진전시켜 반복계산 은 필요없지만 해의 정확도 감소를 피할 수 없었다. 그러나 본 연구에서 제시한 implicit 기법은 전체 계방정식이 비선형 시 스템이 되어 반복계산 과정이 필요하지만, 실제로 수치예제를 통해 검증해 본 결과 계산량의 큰 증가없이 해석의 정확도를 획기적으로 향상시켰다. 이동하는 미분불연속 특이성을 갖는 융해(melting)문제를 수치계산한 결과, implicit 이동최소제곱 차분법을 통해 2차정확도를 얻을 수 있음을 보였다.

논문은 2차원 선형탄성 직접 경계요소법에서 저매개변수 요소를 사용할 때 Kernel의 적분방법에 대하여 논의하였 다. 일반적으로 등매개변수 요소의 경우 형상함수로 통칭되는 해의 기저함수와 요소의 적분을 위해 사용되는 사상함수를 동일하게 사용한다. 그러나 본 논문에서는 사상함수의 차수를 낮게 취하여 순수기저절점을 도입하고 그때 직접 경계요소의 Kernel을 적분하기 위한 방법이 모색되었다. 일반적으로 경계요소법의 적분 Kernel의 경우 Log수치적분과 코쉬주치 (Cauchy principal value) 등을 통해 해결하는데, 본 논문에서는 대수적 조작을 통해 적분값의 정확도를 높일 수 있도록 새 로운 수식을 유도하였다. 본 연구에서 저매개변수 기반의 직접 경계요소에 대한 강건성과 정확도를 검증하기 위해 2차원 타원형 편미분방정식으로 표현되는 평면응력과 평면변형문제에 대해 적용하였다. 적용 예제로는 단순연결영역(simple connected region)의 대표적 문제인 캔틸레버보와 다중연결영역(multiple connected region)의 대표적인 문제인 개구부가 있는 사각평면에 대해 각각 수치해석을 수행한 결과 대폭적인 자유도의 감소에 비해 정확도 측면에는 기존의 방법과 차이 가 없음을 볼 수 있었다. 본 논문에서 제시된 방법은 기저함수 고차화 저매개변수 직접 경계요소법(subparametric high order boundary element)과 이에 기초를 둔 저매개변수 고차 이중경계요소법(subparametric high order dual boundary element)의 초석이 될 수 있을 것이다.

본 연구에서는 뒤틀림을 받고 있는 정사각형 단면의 중공단면 보를 최소자승법과 경계요소법을 이용하여 수치 해석하고 구조물을 해석하였다. 임계하중은 하중을 점차적으로 증가하여 구조물이 파괴가 발생하여 안정성을 상실하는 상태에서 가 장 작은 하중을 의미한다. 뒤틀림을 받고 있는 beam은 일반 구조물에서 많이 발생하는 현상이며, 구조물의 안정성에 크게 영향을 미치고 있다. 최소자승법과 경계요소법은 복잡한 구조물에서도 물론, 다양한 경계조건을 포함하는 문제에 이르기까 지 구조물의 안정성을 검사하는데 효과적인 수치해석 방법이다. 특히 뒤틀림의 문제에서는 단순성 및 일반성에 기인하여 매우 적합한 해석방법이다. 본 연구에서는 뒤틀림을 받고 있는 중공단면 보의 해석해를 유도하여 최소자승법으로 수치 해 석하고 또한 경계요소법을 적용하여 빔의 안정성을 비교 검토하였다. 개발한 컴퓨터 프로그램의 타당성을 증명하기 위하여 삼각형, 사각형 그리고 타원형 단면에 대하여 각각 해석하여 해석해와 비교․검토하였다.

본 논문은 로터 시스템의 디스크 회전운동을 표현하는데 있어 운동방정식을 통합하는 과정에서 기존 연구자들이 채택한 오일러 각 사용법이 일관성이 없음을 지적하였다. 기존 연구자들은 오일러 각 순서가 달라서 속도와 운동에너지도 달리 산정하였음은 물론, 운동방정식은 오직 선형 시스템만 취급해 왔다 이러한 오일러 각 사용법의 단점을 극복하기 위하여 회전운동을 더욱 단순하게 매개화할 수 있는 4원법(quaternion)과 구 좌표계를 적용하여 비선형 시스템을 도출하였다. 이를 바탕으로 수치해석을 통하여 기존 방법과 비교하여 제안한 방법의 신뢰성과 우수성을 보였다.

본 연구는 2차원 및 3차원 동적 탄소성 응력 해석을 위한 특수 적분해 경계요소법의 공식 개발을 제시한다 정적 탄성에 대한 기본식이 일반해를 구하는데 이용되었으며, 전체형상함수 개념을 이용하여, 변위율과 traction rate의 특수 적분해를 구함으로써 지배 방정식의 가속도 부분을 근사화시켰다. 시간 적분을 위하여 Houbolt 시적분 방법을 이용하였으며, Newton-Raphson 알고리즘을 이용하여 수치 연산을 행하였다. 제시된 공식에 따른 예제 해석을 통하여 그 방법의 유효성과 정확성을 설명하였다.

The boundary element alternating method is applied to two dimensional and axisymmetric racture mechanics problems. Different crack depths are considered with particular emphasis on cracks. Several methods of calculating the stress intensity factor are presented and applied to several practical examples to demonstrate their accuracy.

본 연구는 계면경계에서 특이성을 갖는 이종재료 열전달문제를 효율적으로 해석할 수 있는 이동최소제곱 유한차분법을 제시한다 이동최소제곱 유한차분법은 격자망(grid)없이 절점만으로 이동최소제곱법을 이용하여 Taylor 다항식을 구성하고 차분식을 만들어 미분방정식을 직접 푼다. 초평면함수 개념에 근거한 쐐기함수를 이동최소제곱 센스(sense)로 근사식에 매입하여 쐐기거동과 미분 점프에 따른 계면경계 특성을 효과적으로 묘사하고 고속으로 미분을 근사하는 이동최소제곱 유한차분법의 강점을 발휘하도록 했다. 서로 다른 열전달계수를 갖는 이종재료 열전도문제 해석을 통해 이동최소제곱 유한차분법이 계면경계문제에서도 뛰어난 계산효율성과 해의 정확성을 확보할 수 있음을 보였다.

대기의 밀도 성층화는 대기 과학적 측면에서 다양한 현상을 유발하는 인자로 인식되어 오고 있다. 대기의 밀도 성층화와 지형의 상호 작용에 의한 현상은 비단 대기 과학적인 관심일 뿐 만 아니라 산악주위를 비행하는 항공기의 갑작스런 부침에도 영향을 미치며 간접적으로는 대기의 온도 역전현상으로 인해 도시 지역에 발생하는 대기 오염물질의 갈색 구름 (brown cloud)의 형성에도 영향을 미친다. 이와 같이 밀도 성층화는 건축, 토목, 환경, 해양, 기계와 관련된 유동에서 매우 중요한 역할을 하는 유동유발 인자이다. 하지만 밀도 성층이 관련된 유동의 해석은 국내에서는 극히 미미한 실정이고 외국에서는 실험을 중심으로 비교적 많은 연구가 이루어지고 있다. 본 연구에서는 3차원 지형 주위의 밀도 성층 유동 현상을 이해하기 위하여 전산 유체역학을 통한 해석을 수행하였다. 수치 기법으로 가상 경계법을 사용하여 반구(hemisphere)와 3차원 축대칭 언덕(Gaussian hill) 주위 성층유동에 대한 해석을 수행하여 Re 및 Fr 와 장애물의 형상 변화가 유동장에 주는 영향에 대해 중점을 두었다. 한 개의 장애물뿐 만 아니라 2개의 언덕이 주유동 방향으로 배열되어 있을 때 그 거리에 따른 상호작용이 연구되었으며 층류 영역에서 연구가 진행되었다. 연구 결과 Fr<0.4일 때 장애물 후류에서 와흘림이 발생하며 성층이 강해질수록 Re의 효과는 약해지고 와흘림의 성격은 일률화 된다. 특히 Fr=0.2일 때 St의 값이 형상에 조차 의존하지 않고 약 0.19로 수렴하는 결과를 낳았다. 상부로 갈수록 수평단면적이 급격히 감소하는 3차원 언덕의 형상적 특징은 Re를 감소시키는 효과를 발생시킨다. 이런 효과는 비성층 유동에서 잘 나타나고 있으며 같은 Re 임에도 불구하고 3차원 언덕 주위의 유동 현상과 반구 주위의 유동현상이 서로 같지 않은 것을 알 수 있다. 뿐만 아니라 3차원 언덕의 형상적 특징은 강성층에서 높이에 따른 와구조의 위상 변화를 유발한다. 반구의 경우 장애물의 상하부에서 거의 동시에 와류가 활성화되는 반면 3차원 언덕의 경우에는 언덕의 하부보다 상부에서 먼저 와류가 형성됨을 확인하였다.

본 논문에서는 다양한 문제들의 형상 설계 민감도 해석에 대한 효율적인 경계기반 기법을 제시하였다 우선 문제에서 정의된 일반적인 함수들에 대한 연속체 형태의 식에 근거하여, 경계 적분 형태의 해석적 민감도 식을 유도하였다. 이 식은 다양한 형상 설계 문제들의 경사를 계산하는데 편리하게 사용할 수 있다. 그리고 경계법은 형상 변분 벡터가 전체 도메인이 아닌 경계에서만 요구된다는 장점이 있는데, 여기서 경계 형상 변분은 형상 함수의 복잡한 해석적 미분 대신 형상을 미소 증분시킴으로써 편리하게 계산할 수 있다. 제시한 방법의 효율성을 보이기 위해 포텐셜 유동 문제와 필렛(fillet)에서의 응력 집중 문제에 이를 적용하였다.