• 한국해안해양공학회지
• /
• 제17권3호
• /
• pp.166-177
• /
• 2005

본 연구에서는 정체수역에서 유입이론과 일부 형상계수의 조정을 통해 근역에서의 하$\cdot$폐수 혼합거동을 해석할 수 있는 근역제트적분모형을 개발하기위해 총 6개의 상미분 보존방정식에 6개의 미지수를 가지는 문제를 수치적으로 풀기 위한 4차의 Runge-Kutta기법을 사용하였다. 또한, LIF 시스템을 이용하여 검정과정을 통해서 단일수평부력제트의 수리실험을 수행하였다. 그리고 기존의 모형 CORMIX 1, WSJET,그리고 본 모형의 계산결과를 수리실험결과와 서로 비교하였다. VISJET모형에 의해 예측된 중심선 제적이 운동량과 부력이 지배적인 구간에서는 실험결과와 근접한 반면에 본 제트적분모형에 의해 예측된 결과는 천이영역에서 측정된 궤적과 잘 일치하였다. 중심선희석률에 있어서 운동량과 부력이 지배적인 구간에서 CORMIX1 모형의 결과와 잘 일치하는 반면에 초기영역과 천이 영역에서 본 모형의 결과와 대체로 잘 일치하는 경향을 보였다.

• 한국항공우주학회지
• /
• 제38권10호
• /
• pp.973-978
• /
• 2010

본 논문에서는 유격을 고려한 항공기 전방 착륙 장치의 쉬미 현상을 연구하였다. 쉬미는 항공기의 이착륙 시 랜딩기어가 주행도중 측방향과 조향방향으로 진동에 놓이는 현상이다. 이 현상은 스트럿의 낮은 강성, 랜딩기어 내부의 마찰과 유격, 휠의 불균형이나 마모된 부품 등으로 인해 발생하며, 항공기의 안정성을 저하시킨다. 유격은 비선형 요소이기 때문에 기술 함수로 선형화 하여 주파수 영역에서 안정성해석을 수행하였고, 4차 Runge - Kutta를 이용하여 시간영역에서 안정성해석을 수행하였다. 본 연구에서는 수치적인 해석법을 통해 쉬미현상의 선형 동특성과 비선형 동특성을 조사하였다. 유격을 고려한 비선형 수치 해석결과, 선형 임계속도보다 낮은 속도에서 제한주기진동이 발생하는 등, 유격으로 인해 쉬미 안정성이 저감되는 결과를 관찰하였다.

• 에너지공학
• /
• 제26권2호
• /
• pp.64-72
• /
• 2017

중동지역 환경에 적합한 70 W급 LED 조명등의 성능분석을 집중변수모델을 이용하여 수행하였다. LED 조명등은 발열기판, 히프파이프, 방열판으로 구성된다. LED 조명등을 4개의 물체로 구분하고, 각각에 대해 에너지평형을 적용하여, 4개의 연립 비선형미분방정식 형태의 집중변수모델을 수립하였다. 연립 방정식의 해는 Runge-Kutta법을 이용하여 구하였다. 집중변수모델의 대류열전달계수는 다차원해석을 통하여 구하였으며, 실험결과와 비교한 결과, 발열기판은 $1.5^{\circ}C$, 상부방열판에서 $1.8^{\circ}C$의 오차를 가지며, 상대오차는 약 0.6 %임을 확인하였다. 이 모델을 이용하여 대기온도가 $55^{\circ}C$인 정상 운전조건, 태양광만 주어질 때의 조건, 태양광이 주어진 상태에서의 비정상 운전조건, 상부방열판이 없는 경우 등에 대한 온도분포분석을 수행하였다.

• Journal of Astronomy and Space Sciences
• /
• 제2권1호
• /
• pp.19-33
• /
• 1985

3체 문제에 있어서 Newtonian 운동방정식의 특이접은 수치적분에 의한 특수해를 구하는 과정에서 정확도를 떨어뜨리고 computer 사용시간을 증가시킨다. 이러한 특이점은 Newtonian 운동방정식의 독립변수와 좌표축을 변환하는 정칙화(regularization)과정을 통해 제거할 수 있다.이 논문에서는 정칙화된 Newtonian 운동방정식을 $5^{th}$ -order Runge-Kutta 방법으로 특수해를 구하기 위해 BASIC 언어로 작성한 computer program 으로 적절한 초기위치와 초기속도를 가정하여 시간에 따른 3체의 위치와 속도를 계산하였다. 그 결과, 이러한 3중성계의 진화는 결국 분열되어 3체 중 1개는 쌍곡선 궤도를 그리면서 계를 탈출하고 나머지 2개는 연성계를 형성하게 되었다. 이는 연성계의 기원을 설명할 수 있는 하나의 방버이 되지 않을까 생각한다.

• 한국산학기술학회논문지
• /
• 제13권12호
• /
• pp.5689-5695
• /
• 2012

본 연구에서는 차량에 적용되는 시트 재질인 폴리우레탄 폼의 점탄성 특성을 고려하여 시트와 인체의 진동특성을 시험 및 수치해석 방법을 이용하여 분석하였다. 압축 시험을 통해 폴리우레탄 폼의 점탄성 특성인 비선형성과 준-정역학적 특성을 구하였다. 또한 컨벌루션 적분법 및 비선형 강성 모델을 이용하여 폴리우레탄 폼의 점탄성 특성을 수학적으로 모델링하였다. 시트의 승차감 기여도를 분석하기 위하여 시트의 동역학 모델과 ISO5982의 표준 인체 수직진동 모델을 이용하여 수직 진동모델을 구성하고 관련 운동방정식을 유도하였다. 비선형 운동방정식은 Runge-Kutta 적분법을 이용하여 수치해석 시뮬레이션을 수행하였다. 철도차량의 차체 바닥에서 측정한 진동가속도 입력에 대한 시트와 인체의 응답 특성을 분석하고 시트 설계 파라미터에 대한 승차감 지수 값들의 변화를 분석하여 시트 설계에 대한 방법론을 제시하고자 한다.

• 대한토목학회논문집
• /
• 제6권3호
• /
• pp.1-9
• /
• 1986

동하중(動荷重)을 받는 양단(兩端) 힌 지 포물선(抛物線) 아치의 연동방정식(連動方程式)을 Runge-Kutta 방법(方法)으로 수치해석(數値解析)하므로써 동적(動的) 임계하중(臨界荷重)을 구했다. 낮은 양단(兩端) 힌 지 포물선(抛物線) 아치에 step하중(荷重)과 impulse하중(荷重))이 작용(作用)하는 경우에 관해 Budiansky-Roth criterion을 적용하여 동적(動的) 임계하중(臨界荷重)을 정의(定義)하고, 이를 상관곡선(相關曲線)으로써 동적(動的) 안정영역(安定領域)을 제안하였다. 포물선(抛物線) 아치에 대하여 얻어진 결과를 정현(正弦) 아치의 경우와 비교(比較)하여 아치의 기하학적(幾何學的) 형상(形狀)이 동적(動的) 안정영역(安定領域)에 미치는 영향을 밝혔고, 동적(動的) 안정영역(安定領域)은 아치의 높이에 큰 영향을 받는다는 것을 밝혔다.

• 멤브레인
• /
• 제16권3호
• /
• pp.204-212
• /
• 2006

이산화탄소에 대한 가소화 안정성 및 이산화탄소/질소의 분리특성이 탁월한 폴리이서술폰(PES)중공사 분리막에 의한 이산화탄소 분리특성을 수치해석으로 알아보고자 하였다. 공정변수에 따른 이산화탄소 분리 거동을 예측하기 위하여 공급 기체와 투과기체가 같은 방향으로 흐르는 병류 흐름에 대한 분리막 공정 지배 방정식을 5차 Runge-Kutta-Verner 방법을 사용하여 Compaq Visual Fortran 6.6 소프트웨어를 이용 공정모사 프로그램을 개발하였다. 개발된 프로그램을 사용하여 수치해석을 수행한 결과, 이산화탄소 투과특성에 영향을 주는 가장 중요한 인자로서 공급 이산화탄소 분압, 투과측과 분리막 내부의 압력비 그리고 공급 기체가 분리막 내부에 머무르는 체류 시간임을 알 수 있었다.

• 한국항공우주학회지
• /
• 제44권4호
• /
• pp.281-289
• /
• 2016

시간에 따라 질량이 감소하는 낙하산 시스템의 궤적 및 낙하 시간 분석은 정밀한 투하가 요구되는 임무에 중요하므로 그 필요성이 더 커지고 있다. 본 연구에서는 질량 변동 물체인 조명탄을 투하하기 위한 십자형 낙하산 시스템의 동적 거동을 분석하는 연구를 수행하였다. 낙하산 시스템의 궤적을 분석하기 위해 유도된 상미분 형태의 운동방정식 시스템을 Runge-Kutta 수치기법을 적용하여 해석하였다. 그리고 동역학 방정식의 핵심적 입력정보인 십자형 낙하산과 조명탄의 항력 계수를 예측하기 위해 전산유체역학 해석을 수행하였다. 마지막으로 단순화된 대기교란 모델을 적용하여 풍향, 풍속에 따라 달라지는 낙하산 시스템 거동의 차이를 분석하였다.

• 융합정보논문지
• /
• 제10권7호
• /
• pp.131-137
• /
• 2020

가진 정현파 강제력에 의한 강제 진동의 특성이 연구되었고 또한, 주파수 영역에서 FRF의 수치 분석을 자세히 수행하였다. 이와 관련하여 강제 진동 모델에서 변위, 속도 및 가속도의 응답을 구하였다. 고유 주파수 주변의 실수부와 허수부의 FRF 특성도 각 경우에 따라 구하였다. 시간 영역에서의 강제 진동의 응답분석은 정현파 강제 진동의 특성을 식별할 수 있다. 변위, 속도 및 가속도 등의 응답을 얻기 위해 Runge-Kutta-Gill 방법의 수치해석 기법을 수행하여 강제력 주파수에 따른 응답을 얻었고 이 주파수는 응답에 큰 영향을 미치지 않았다. 또한, FRF는 강제 진동의 고유 특성을 나타내고 있으며 이러한 강제 진동 모델의 각 감쇠 조건에서 이러한 응답분석을 성공적으로 자세하게 얻을 수 있었다. 상이한 질량, 감쇠 및 강성에 대한 수치 분석 후, 정현파 강제력에 의한 강제 진동 응답 특성을 강제력의 진폭 및 주파수를 동시에 고려하여 분석되었다.

• Structural Engineering and Mechanics
• /
• 제12권6호
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.