We present the whole basis of numerical method and useful formulae for general relativistic magnetohydrodynamic simulations in Kerr space-time.

서로 다른 기지의 성질을 갖는 재료들을 혼합하여 만든 합성재료의 새로운 물성치는 일반적으로 실험으로 규명하고 있다. 혼합하는 재료들의 체적비에 따라 실험으로 측정한 합성재료의 탄성계수와 포와송비는 그 합성재료로 만들어지는 구조물의 역학적 거동을 예측하는 해석적 모델의 기본자료로 사용된다. 합성재료 탄성물성치의 수치적 예측은 합성재료에 대한 유한요소 모델로 해석한 정적변위와 균질.등방성으로 가정한 모델을 해석한 정적변위와의 차이를 최소화하는 구속적 비선형 최적화기법을 사용하여 수행하였다. 유한요소 모델은 체적비에 따라 혼합물질을 분배하기 용이하도록 제안하였으며 구속조건 및 하중조건은 일축인장에 의한 거동을 예측하도록 설정하였다. 본 논문에서는 고체입자를 섞어 만든 합성재료의 탄성물성치를 예제를 통하여 수치적으로 예측하고 그 결과를 실험결과 및 이론식들과 비교.검토하였다.

본 논문은 기하학적 비선형과 재료적 비선형을 고려하여 강구조 부재 및 골조의 비선형 해석을 위한 프로그램을 개발하고 강재의 응력도 변형도 관계를 정식화하였다. 본 프로그램의 효율성을 검증하기 위하여 단조 하중을 받는 H형강 보와 각형강관 기둥, 그리고 반복 수평력을 받는 브레이스 강재 라멘에 대한 수치해석을 행하였다. 본 그로그램에 의해 얻어진 결과는 대체적으로 기존의 실험 및 해석결과와 일치하였다.

We have calculated the astronomical almanac 1994 and simulated the trajectory of a satellite orbit considering all perturbative forces with various initial conditions. In this work, Gauss Jackson multistep integration method has been used to calculate our basic equation of motion with high numerical accuracy. It has beer. found that our results agree well with the Astronomical Almanac Data distributed by JPL of NASA and the orbit simulations have been carried out with fast speed, stability and excellent round-off error accumulation, comparing with other numerical methods. In order to be carried out our works on almanac and orbit calculations easily by anyone who uses a personal computer, we have made a computer program on graphical user interface to provide various menus for detail works selected by a mouse.

This paper discusses the hydrodynamic characteristics of a catamaran at low speed. In this study, the Delft 372 catamaran model was selected as the target hull to analyze the hydrodynamic characteristics by using the RANS (Reynold-Averaged Navier-Stokes) numerical method. First, the turbulence study and mesh independent study were conducted to select the appropriate method for numerical calculation. The numerical method for the CFD (Computational Fluid Dynamic) calculation was verified by comparing the hydrodynamic force with that obtained experimentally at high speed condition and it rendered a good agreement. Second, the virtual captive model test for a catamaran at low speed was conducted using the verified method. The drift test with drift angle 0-180 degrees was performed and the resulting hydrodynamic forces were compared with the trends of other ship types. Also, the pure rotating test and yaw rotating test proposed by Takashina, (1986) were conducted. The Fourier coefficients obtained from the measured hydrodynamic force were compared with those of other ship types. Conversely, pure sway test and pure yaw test also were simulated to obtain added mass coefficients. By analyzing these results, the hydrodynamic coefficients of the catamaran at low speed were estimated. Finally, the maneuvering simulation in low speed conditions was performed by using the estimated hydrodynamic coefficients.

RC-SC 이종부재 접합부는 SC구조의 적용구간에 따라 RC-SC벽체 수평면 접합과 RC-SC벽체 연직면 접합 등 다양한 접합부 형태가 발생한다. 그러나 이종부재 접합부의 다양한 형태를 실증실험으로 모두 검증하기에는 비용과 시간 등 많은 제약사항이 존재한다. 선행연구에서 유한요소해석을 통해 다양한 이종부재 접합부 형태의 구조 건전성을 확인하고자 접합부 해석방법론을 개발 및 검증한바 있다. 본 연구에서는 RC-SC이종부재 접합부 해석 방법론을 개발하는 과정에서 선 수행한 요소별 해석결과, 제한조건 및 마찰계수별 해석결과 등 해석 시 요구되는 단계별 변수 해석결과를 분석/결정하여 RC-SC 이종부재 접합부 해석방법에 활용할 최적화된 변수를 선정하였다.

현재 전 세계적으로 화석연료의 고갈과 원자력에너지에 대한 불안감 등으로 인하여, 풍력, 조력, 조류, 지열, 태양광 등과 같은 신재생 에너지 개발을 위한 노력이 증가하고 있다. 신재생 에너지 중 발전효율이 높은 풍력에너지에 대한 연구가 활발히 진행되고 있으며, 해상풍력발전의 발전단가가 작게 산정되어 충분한 경쟁력을 확보하고 있다고 판단되며, 향후 많은 해상풍력발전 단지의 건설이 추진되고 있다. 해상에서는 육상과 비교하여 동일 풍속 대비 낮은 고도에서 양질의 풍력을 얻을 수 있기 때문에, 해상풍력발전의 단지건설이 증가하고 있으며, 발전용량의 증대를 위해 발전 터빈, Blase 및 지지구조물의 크기가 대형화되는 추세이다. 하지만, 해상풍력의 대형화에 따라, 강재로 제작된 타워 구조물의 세장비 또한 크게 증가하며, 지지구조의 좌굴 파괴 위험성도 증가하고 있다. 따라서 이의 해결을 위한 경제적이고 안전한 고강도 신형식 지지구조의 개발이 필요하다. 이에 따라 지지구조체에 대한 정확한 거동 모사를 할 수 있는 해석기법에 관한 연구가 활발히 진행중이다. 본 연구에서는 강재 풍력타워의 좌굴 문제 해결과 모듈형 풍력타워 구조체의 개발을 위해 이중관 충전콘크리트(DSCT ; Double-Skinned Concrete Fille Tube) 기둥의 비선형 거동을 모사하기 위하여 범용 유한요소 해석 프로그램을 이용한 해석을 수행하였으며, 실험 data와의 비교를 통하여 타당성을 검증하였다.

화재로 인한 건축물 구조부재의 파괴는 추가적인 경제적 손실과 인명피해를 가져오기 때문에 화재에 대한 구조부재의 안전성은 매우 중요하다. 화재에 노출된 철근콘크리트 구조부재의 경우, 부재의 열전도도에 따라 내부 온도가 상승하고 부재의 하중저항강도와 강성이 저하된다. 철근콘크리트 부재의 화재에 대한 구조성능 확보 여부는 일반적으로 표준화재실험을 통하여 입증한다. 콘크리트의 열전도도는 온도에 따라 변하기 때문에 20℃, 100℃, 200℃, 300℃, 400℃, 500℃, 600℃ 700℃, 800℃ 열적 정상상태의 콘크리트 블록을 이용하여 측정한다. 본 연구에서는 철근콘크리트 기둥의 화재에 대한 구조성능 확보 여부를 확인하기 위한 표준화재실험으로부터 얻어지는 열적 비정상상태의 일차원 온도분포 시간이력으로부터 콘크리트의 열전도도를 추출하는 방법을 제안하였다. 표준화재실험으로부터, 기둥의 표면에서 중심까지의 온도분포 시간이력을 얻었다. 콘크리트의 온도에 따른 열전도도 k(T)=α{(T/120)²-T/60+200}에 대하여 다양한 α값으로 온도분포 시간이력을 수치해석 하였고, 실험결과와 가장 오차가 작게 발생하는 α값을 결정하였다. 결정된 콘크리트의 온도에 따른 열전도도는 20℃, 100℃, 200℃, 300℃, 400℃, 500℃, 600℃ 700℃, 800℃ 열적 정상상태의 콘크리트 블록에서 측정된 콘크리트의 열전도도와 비교하였고, 제안된 수치해석에 의한 열전도도 산출방법이 유효함을 입증하였다.

토석류(debris flow)와 진흙흐름(mudflow)은 대부분 집중호우에 의해서 유발되며, 사면을 따라서 진흙흐름이나 토석류가 유하할 때, 강우의 직접 유출수 또는 선행강우에 의한 유출수가 존재하는 경우가 많다. 이 경우 토석류의 유동특성을 현저히 변화시킬 수 있는 사면 바닥에서의 강우 유출수의 연행(entrainment) 현상을 토석류 및 진흙흐름 발달과정에서 고려해야만 합리적이고 정확한 수치 모델링이 가능하다. 이 연구에서는 경계면 포착(interface-tracking)기법을 이용한 계산유체동력학 기법을 이용하여 진흙흐름과 토석류 흐름을 재현하기 위환 수치 모델링을 수행한다. 토석류 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축(interface reconstruction) 기법을 이용한다. 지배방정식은 비압축성(incompressible) 질량보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. USGS의 연구팀이 대형수로에서 수행한 실험을 통해서 확보한 토석류를 재현함으로서 모형의 적용성을 평가한다.

Linear and nonlinear models are available to simulate the effect of mine blasting. And the results are changed by selection of each model. Plaxis 8.0, FLAC 4.0 are used in this study for simulating blasting pressure. These programs offer linear model and nonlinear model. Linear model can properly simulate loss of impact wave in plastic failure region based on comparison of site vibration measurement and numerical analysis iteratively.

최근 물질의 특성과 구조를 해석하기 위해서 수치해석학적 모델과 컴퓨터 시뮬레이션이 많이 사용된다. 이러한 관점에서 물질의 미세구조를 해석하는 데 있어 분자동역학 해석법은 매우 유용한 방법이다. 이번 연구에서는 점토광물에 대한 확산계수 및 점착력과 같은 물리화학적 특성을 계산하기 위한 수치해석학적 방법을 소개한다. 이번연구에서 지질학적으로 심부에 위치하는 포화된 점토광물과 물을 포함한 점토광물에 대한 특성을 분자동역학을 이용해서 계산하고 균질화해석법을 활용하여 점착력과 같은 외부조건에 따라 결정되는 점토광물의 확산거동을 해석하였다. 그 결과 수치해석에 의한 해석결과 값과 기존의 실내 투수실험에 의한 결과 값이 매우 흡사한 결과를 보인다는 것을 알 수 있다. 이는 여러 가지 복합적인 조건하에서의 점토광물의 물리화학적 거동을 해석하는 데 수치해석학적 방법이 매우 유용하게 사용될 수 있다는 것을 의미한다.

In this study, specimens were produced to reach design standard strengths 40, 50 and 60MPa on age 28 days to evaluate the usability of high strength recycled aggregate concrete and the applicability of compressive strength by Impact Echo method and numerical analysis, and after 0, 30, 50 and 100% recycled coarse aggregates were substituted and mixed into each combination, totally 12 combinations were set. And, through compressive strength test, Impact Echo method, shock reverberation technique and numerical analysis on 1, 3, 7 and 28 days of age

항복강도 690MPa 이상의 초고강도 고성능강재 HSB800 후판이 최근 국내에서 개발되었으나, 플레이트 거더 등 주요 구조물에 최적 활용되기 위한 설계기술 기반은 아직 마련되지 못하고 있는 실정이다. 고강도 강재는 일반적인 강재에 비해 비교적 낮은 연신율과 큰 항복비를 갖고 있어 구조연성 성능에 대해 다소 불리한 것으로 예상되며, 현행 AASHTO LRFD 설계기준의 비탄성 설계과정에서도 공칭항복응력 485MPa를 초과하는 단면 설계에 대한 적용이 엄격히 제한되어 있다. 본 연구에서는 HSB800시편에서 표준시험을 통해 획득한 비탄성 응력-변형율 특성을 고려한 3차원 비선형 유한요소 해석 방안을 정립해 보고자 한다. 관련 규준의 표준화 연구를 위해 재료특성의 표준모형 정립이 선행되어야 한다. 기존 일반강에 적용하던 완전탄소성 거동 및 다중 선형 재료특성 곡선과 실제 재료특성을 적용한 비탄성 연성거동 해석결과에 대해 상호 비교연구를 수행하였다. 향후 실 구조물 규모의 실험적 검증을 통해 본 해석적 모형이 검증된다면, HSB800을 적용한 플레이트 거더의 비탄성 설계에 관한 광범위한 연구를 효과적으로 수행할 수 있을 것으로 판단된다.

적응적 메쉬세분화기법과 분할격자기법을 적용한 2차원 천수방정식모형을 활용하여 구조물을 고려한 극한 홍수 실험을 모의하였다. 본 연구에 사용된모형은 두 격자생성 기법을함께 사용함으로서 복잡한 경계를보다 적은 격자로 효율적으로 표현하는 것이 가능하며, 동적 적응 메쉬세분화기법을 사용하여 흐름이 빠르게 변하는 영역에서 정확도를 유지하면서도 효율적으로 계산하는 것이 가능하다. HLLC 리만근사해법과 MUSCL 기법을 적용하여 시공간상에서 2차정도를 유지하며,

The main objective of this research was to evaluate the effects of process variables which were forming ability, flow displacement, effective stress, effective strain, fluid vector and products defects on manufactured artificial lightweight aggregate made of both incinerated sewage sludge ash and clay by means of the numerical analysis of a rigid-plastic finite element method. CATIA (3D CAD program) was used for an extrusion metal mold design that was widely used in designing aircraft, automobile and metallic molds. A metal forming analysis program (ATES Co.) had a function of a rigid-plastic finite element method was used to analyze the program. The result of extrusion forming analysis indicated clearly that a shape retention of the manufactured artificial lightweight aggregate could be maintained by increasing the extrusion ratio (increasing compressive strength inside of extrusion die) and decreasing the die angle. The stress concentration of metal mold was increased by increasing an extrusion ratio, and it was higher in a junction of punch and materials, friction parts between a bottom of the punch and inside of a container, a place of die angle and a place of die of metal mold. Therefore, a heat treatment as well as a rounding treatment for stress distribution in the higher stress concentration regions were necessary to extend a lifetime of the metallic mold. A deformity of the products could have made from several factors which were a surface crack, a lack of the shape retention and a crack of inside of the products. Specially, the surface crack in the products was the most notably affected by the extrusion ratio.

Jindong bay at the northwestern part of Jinhae bay suffers from the occurrence of red tides in summer every year. In order to study the management methods of coastal environments, an ecological numerical model has been developed. The model experiments was forecasted that the load of nutrients from the land and field concentration will be cut down per 10% each. When we cut down 57.2% nitrogen load in the inner bay and 38.4% phosphorous load in the outer bay of bottom layer of the nutrients load from land and field concentration, the seawater quality standard levels up first grade. When we cut down 86.5% nutrients in the inner bay and 93.0% nutrients in the outer bay, the concentration of chlorophyll a decreases below 3.2(an individual concentration of phytoplankton : 10,000cell/ml), i.e. the red tides do not occur.

비보존성 오염물질의 종확산에 관한 수치모형을 개발하였다. 계산기법으로는 종확산 방정식을 이송, 감쇠 및 확산 방정식으로 분리하고, 이들 방정식을 1/3 시간 간격에 대하여 번갈아 계산하는 단계분리 유한차분기법을 사용하였다. 이송방정식에 대해서는 Holly-Preissmann 기법을, 감쇠방정식에 대해서는 해석적 방법을, 확산방정식에 대해서는 Crank-Nicholson 기법을 각각 사용하였다. 오염물질이 불균일 흐름 내로 연속적으로 유입되는 경우 및 균

The ocean wave is hydrodynamically investigated to get more reliable solution. To improve the computational accuracy more fine grids are used with relatively less computer storage on the free surface. One element of the free surface is discretized into more fine grids because the free-surface waves are much affected by the grid size in the finite difference scheme. Here the multi-grid method is applied to confirm the efficiency for the S103 ship model by solving the Navier-Stokes equation for the turbulent flows. According to the computational result approximately 30% can be improved in the free surface generation, Finally the limiting streamlines show numerical result is similar to the experiment by twin tuft.

모멘트법은 Lagrangian 벙법으로서 격자요소 내에서의 농도의 공간분포에 대한 0차, 1차, 2차 모멘트를 고려한고 각 모멘트의 보존성을 유지하면서 농도분포의 이송을 계산하는 방법이다. 따라 각 격자요소에서의 0차 모멘트, 즉 평균농도 뿐만 아니라 1차 및 2차 모멘트 값의 합리적인 초기 설정이 요구된다. 본 연구에서는 각 모멘트들의 초기값 설정방법을 검토하고, 기존 모멘트법의 Couuant 수에 대한 제약조건을 극복하기 위하여 모멘트법을 개선하였다