It is known that the growth and development of the mosquito are greatly affected by the change of the meteorological factors. In particular, temperature and precipitation are closely related to the life cycle of the mosquito, and their effects have different characteristics for each species of mosquito. Therefore, to develop a mosquito activity index based on mosquito density, it is essential to develop a prediction model based on weather data. In this study, we developed a functional formula that can estimate the change of mosquito density according to the change of meteorological factors using the mosquito classification data of Incheon region collected from 2011 to 2017. Also, using the data of the digital mosquito monitoring system (DMS) from April to October 2018, mosquito activity index according to characteristics of space in city was developed. In order to reflect the temporal characteristics of the mosquito life-cycle, we used data of temperature and precipitation prior to 1-2 weeks, and used land cover data to reflect the spatial characteristics of mosquito density. Density of Culex pipiens collected in the Incheon area were gradually increased when the average temperature increased two weeks ago after adjusting the precipitation. However, when the average temperature reached 22°C, the density of Culex pipiens marked a peak, and above the 22°C, the density was decreased. The predicted mosquito activity index calculated by applying the machine learning method to the DMS data of the Incheon area is designed to calculate from 1 to 10 grades. The accuracy of the mosquito activity index was 87% when the one grade error was allowed.

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.

PURPOSES: Implementation and verification of the simple linear cohesive viscoelastic contact model that can be used to simulate dynamic behavior of sticky aggregates. METHODS: The differential equations were derived and the initial conditions were determined to simulate a free falling ball with a sticky surface from a ground. To describe this behavior, a combination of linear contact model and a cohesive contact model was used. The general solution for the differential equation was used to verify the implemented linear cohesive viscoelastic API model in the DEM. Sensitivity analysis was also performed using the derived analytical solutions for several combinations of damping coefficients and cohesive coefficients. RESULTS : The numerical solution obtained using the DEM showed good agreement with the analytical solution for two extreme conditions. It was observed that the linear cohesive model can be successfully implemented with a linear spring in the DEM API for dynamic analysis of the aggregates. CONCLUSIONS: It can be concluded that the derived closed form solutions are applicable for the analysis of the rebounding behavior of sticky particles, and for verification of the implemented API model in the DEM. The assumption of underdamped condition for the viscous behavior of the particles seems to be reasonable. Several factors have to be additionally identified in order to develop an enhanced contact model for an asphalt mixture.

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

동역학의 새로운 변분이론인 확장 해밀턴 이론은 수학물리학을 비롯한 공학에 있어 초기치-경계치 문제해석에 광범위하게 적용될수 있는 기반을 제공하는 것으로 본 논문에서는 이 이론을 기반으로 선형탄성 단자유도계에 적용한 새로운 수치해석법을 제안하였다. 곧, 변분이론의 특성을 감안해, 전체 time-step에 대한 수치해를 한번에 산정하는 해석법을 제안하였고, 주요 예제를 통해 이 해석법의 특성을 살펴보았다. 에너지 보존 시스템의 경우(비감쇠 시스템에 외력이 작용치 않는 경우), time-step에 관계없이 에너지와 모멘텀이 보존되는 symplecticity property를 가지고 있음을 확인할 수 있었고, 감쇠 시스템인 경우, time-step이 점점 작아질수록 정확한 해에 빠르게 수렴하는 것을 확인하였다.

논문은 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)의 초석이 될 수 있을 것이다.

막 구조물은 연성의 막에 초기 장력을 주고 외관의 강성을 늘림으로써 외부하중에 안정된 형태를 유지하는 구조물로 두께를 얇게 하여 대공간 구조에 많이 채택된다. 이러한 막 구조는 자유로운 곡선을 표현할 수 있는 특성이 있어, 구조적 형태의 선정은 매우 중요하다. 이에 본 논문에서는 넙스를 기저함수로 하는 비정형 곡면으로 형상을 표현하고, 최적의 곡면 형상 탐색을 위한 대변형 결과값 도출을 위해 기하학적 비선형을 고려한 유한요소해석법을 제안하였다. 또한, 형상 탐색 결과로 나타난 곡면의 형상 근사화의 최소화를 위해 유한 요소망으로 표현된 최종 형상을 다시 넙스로 구현하는 인터페이스 기법을 제안하여, 비정형 막 구조물의 최적 곡면을 표현하였다.

  Ultra-high voltage transformer industry has characteristic of small quantity batch production system by other order processing unlike general mass production systems. In this industry, observance of time deadline is very important in market competitive

무요소기법이 공통적으로 내재하고 있는 수치적분의 부정확성을 해결하기 위해, 페트로프-갤러킨 자연요소법이라 불리는 향상된 자연요소법을 제안한다. 제안된 방법은 라플라스 기저함수를 시도 형상함수로 사용하는 반면, 시험 형상함수로서 델라우니 삼각형이 지지영역이 되는 함수를 새롭게 정의한다. 이러한 접근은 통상적인 적분영역과 적분함수 지지영역간의 불일치를 제거하게 하며, 이는 적용이 편리할 뿐만 아니라 수치적분의 정확성을 보장한다 본 논문에서는 2차윈 선형 탄성의 대표적인 검증문제를 통하여 제안된 방법의 타당성을 검증한다. 비교를 위해 기존의 부브노프-갤러킨 자연요소법과 일정 변형률 유한요소법을 이용한 해석을 동시에 수행한다. 조각 시험과 수렴율 평가를 통해 제안된 기법의 우수성을 확인할 수 있다.

절연 특성이 기존의 SiO2 보다 우수한 500 두께의 SiN4층을 두 단결정 실리콘사이의 절연막질로 채택하고 직접접합시켜 직경 10cm의 Si(100) /500 -Si3N4/Si (100) 기판쌍을 제조하였다. p-type (100) 실리콘기판을 친수성, 소수성을 갖도록 습식방법으로 세척한 두 그룹의 시편들을 준비하였다. 기판전면에 LPCVD로 500 Å 두께의 Si3N4∥Si(100) 기판을 성장시키고 실리론 기판과 고청정상태에서 가접시킨 후, 선형열원의 이동속도를 0.1mm/s로 고정시키고 선형 입열량을 400~1125w 범위에서 변화시키면서 직접접합을 실시하였다. 접합된 기판은 적외선 카메라로 계면 접합면적을 확인하고 razor blade creek opening 측정법으로 세정 방법에 따른 각 기판쌍 그들의 접합강도를 확인하였다. 접합강도가 측정된 기판쌍은 high resolution transmission electron microscopy (HRTEM )을 사용하여 수직단면 미세구조를 조사하였다. 입열량의 증가에 따라 두 그를 모두 접합율은 큰 유의차 없이 765% 정도로, 소수성 처리가 된 기판쌍의 접합강도는 1577mJ/m2가지 선형적으로 증가하였으나, 친수성 처리가 된 기판쌍은 주어진 실험 범위에서 입열량의 증가에 따라 큰 변화 없이 2000mj/m2이상의 접합 강도를 보였다 친수성 처리가 된 기판쌍의 수직단면 미세구조를 고분해능 투과전자현미경으로 각인한 결과 모든 시편의 실리콘과 Si3N4사이에 25 Å 정도의 SiO2 자연산화막이 존재하여 중간충 역할을 함으로서 기판접합강도를 향상시키는 것으로 판단되었다.

실온에서 직접 접합된 실리콘 기판의 접합강도를 향상기키기 위하여 기존의 고온 로내 열처리법을 대체할 수 있는 선형 열처리법을 개발하였다. 한 개의 열원과 타원형 반사경으로 구성된 선형 열처리법은 접합면의 간격이 열처리 온도의 증가와 더불어 감소하는 특성과 온도 증가와 더불어 접합면에 생성된는 기체상의 밀도가 증가하는 현상을 응용하여 접합면의 기체상을 밀도차이를 이용하여 기판 외부로 방출시키는 방법으로 SimidmidSi 기판쌍 및 SimidmidSiO2/Si 기판쌍의 직접 접합에 적용하여 보았다. IR camera와 HRTEM으로 직접 관찰한 접합면은 실온에서 접합면에 침투한 외부 불순물에 의한 비접합 영역을 제외하고는 자제 생성된 기체상에 의한 비접합 영역은 나타나지 않았고 매우 깨끗한 접합계면을 나타내었다. 접합된 기판쌍을 Crack opening법과 인장시험법을 적용하여 접합 강도를 측정하였다. 접합 강도는 열처리 온도의 증가와 더불어 점차로 증가하였고 두 측정방법 모두 동일한 경향성을 나타내었다.

본 연구에서는 지하철구조물, 터널구조물 및 제방 등과 같은 2차원 지반-구조계의 지진응답해석을 위한 주파수영역 동적해석법을 제시하였다. 제시한 방법에서는 구조물과 구조물 주변 근역지반은 유한요소를 이용하고 원역지반은 주파수종속 동적무한요소를 이용하여 모형화하였다. 지진입력은 입력지진파를 수직으로 입사하는 P-파와 SV-파로 가정하여 자유장응답을 구하였으며 외부고정경계법을 적용하여 등가지진하중을 산정하였다. 본 연구기법의 검증을 위하여, 층상 자유장지반과 균질 반무한지반에 매입된 원동형 공동에 대하여 지진응답을 수행하였다. 이들을 다른 기법에 의한 해와 비교한 결과, 본 연구의 기법이 매우 정확함을 알 수 있었다. 마지막으로 지하철 역사의 지진응답해석 예제를 제시하여 본 연구의 적용성을 보였다.

실제로 구조시스템들의 최적설계는 설계변수가 연속값이 아닌 이산값을 요하는 경우가 대부분이다. 본 논문은 이산형 설계변수를 갖는 비대칭 복합 적층평판에 대해 선형 근사화방법을 이용한 이산최적설계를 수행하였으며, 이 방법이 매우 효율적임을 보였다. 대상 문제는 축력, 전단력, 그리고 휨과 비틀림 모멘트의 평면 내하중들(in-plane loads)의 다중하중조건을 받는 것으로 고려하였으며, 복합 적층평판을 구성하는 플라이들에 대한 최대변형률 규준을 설계 제약조건으로 부과하였다. 이산 최적화를 위한 초기 접근방법으로 단 한번의 연속변수 최적화 과정이 FDM(Feasible Direction Method)을 이용하여 수행되었으며, 차후 이산 및 연속변수를 포함하는 비선형 이산최적화문제를 SLDP(Sequential Linear Discrete Programming)방법에 의해 선형 근사화된 혼합정수계획문제로 형성하여 풀었다. 수치예에서 6개의 플라이로 구성된 비대칭 복합 적층평판을 대상으로 회전식 적층배열([(90-.theta.)/-(60+.theta.)/-.theta./-(45+.theta.)/(45-.theta.)]s)에 따른 이산최적해를 구하였다. 효율성 입증을 위해 똑같은 문제를 비선형 분기한계법을 이용하여 풀었으며, 그 결과를 비교 분석하였다.

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

본 연구에서는 Biot의 선형압밀이론에 근거한 2차원 압밀문제의 근사해를 구하기 위한 경계요소법을 제시한다. 먼저 선형 압밀문제의 기초미분방정식의 시간의존성을 제거하기 위하여 시간에 대한 Laplace변환을 적용시키고, 변환공간에서의 미분방정식을 대상으로 정식화를 한다. 변환공간에서의 변위와 간극수압에 대한 경계적분방정식계를 유도하고, 변환공간에서의 연성문제에 대한 기본해를 구체적으로 보인다. 변환공간에서의 해를 실공간의 해로 변환하기 위하여 Hosono의 수직 Laplace역변환법을 적용하였으며, 해석예로서 2차원 반무한 지반의 국소재하에 의한 압밀문제를 해석예로 선택하였고, 암밀해와 비교하여 제안해법의 적용성 및 타당성을 보였다.