• 제목/요약/키워드: 4차 Runge-Kutta 방법

검색결과 30건 처리시간 0.026초

GLONASS 위성 가시성 분석을 위한 방송궤도력 기반 궤도 예측 (Orbit Prediction using Broadcast Ephemeris for GLONASS Satellite Visibility Analysis)

  • 김혜인;박관동
    • Journal of Astronomy and Space Sciences
    • /
    • 제26권2호
    • /
    • pp.199-210
    • /
    • 2009
  • 다양한 위성항법시스템이 개발 중이지만 현재 측위에 사용 가능한 것은 GPS와 GLONASS 뿐이다. 이 연구에서는 GLONASS의 궤도력 중에서 방송궤도력을 이용하여, 위성의 운동을 나타내는 미분방정식을 4차 Runge-Kutta 방법으로 수치적분하여 위성궤도를 예측하고, 그 정 확도를 평가하였다. 생성한 예측지도는 정밀궤도력과 비교하여 정확도를 검증하였는데, 1일간의 예측궤도와 7일간의 예측제도의 3차원 최대오차는 각각 17.4km, 40.1km로 나타났으며, RMS 오차는 각각 14.3km, 15.7km로 나타났다. 또한 예측제도를 이용하여 산출한 가시위성의 개수와 실제관측 결과를 비교하였다. 그 결과, 관측지점의 주변 건물에 의한 영향으로 발생하는 차이를 제외하고 결과가 일치하는 것을 확인하였다.

유한요소법과 운동방정식의 결합에 의한 운동체의 특성 해석 (Analysis of Moving Body Characteristics by Coupling Finite Element Method and Motion Equation)

  • 김영선;이준호;이기식;이복용
    • 한국자기학회지
    • /
    • 제6권5호
    • /
    • pp.281-286
    • /
    • 1996
  • 전기 기계의 운동 특성을 정확히 해석하기 위해서는 해석모델의 전자장방정식과 운동방정식이 함께 고려되어야 한다. 본 논문에서는 자장내에서 두 개의 스프링에 의해 진동 운동을 하는 도체를 해석함에 있어 순차적인 결합방법을 채책하였다. 한 시점에서 속도가 주어지면 전자장해석 및 유도전류 계산은 유한요소법을 이용하였고, 그때마다 전자력을 게산하였다. 또한 전자장방정식이 풀린다음 4차 runge-kutta 미분방정식 해법을 이용하여 운동방정식과 결합하였다. 위의 과정을 순차적으로 계속 반복함으로써 시간에 대한 위치, 속도, 유도전류 및 전자력을 알아내어 운동체의 운동특성을 알아보았다. 계산된 결과는 에너지보존 법칙에 적용하여 제안한 순차적인 결합방법의 유용성을 검증하였다.

  • PDF

선체의 회전을 고려한 선박과 교량의 간이충돌모델 (A Simplified Bridge-vessel Collision Model Considering with the Rotational Motions of the Vessel)

  • 이계희
    • 대한토목학회논문집
    • /
    • 제31권2A호
    • /
    • pp.43-49
    • /
    • 2011
  • 본 연구에서는 선교나 부선등에 실린 장비들과 같이 선박의 회전중심으로부터 이격된 선박구조와 교량의 상부구조가 충돌하는 경우의 충돌거동을 해석할 수 있도록 선박의 회전을 고려한 간이충돌모델을 개발하였다. 충돌모델은 선박과 교량의 질량과 강도 및 비선형성을 나타내도록 하였고, 모델의 비선형 운동방정식을 4차 Runge-Kutta방법을 이용하여 수치적으로 해석하였다. 현재 국내에서 운용되고 있는 부선의 제원을 선폭을 기준으로 표준화하고 작성된 충돌모델을 이용하여 각종 충돌조건이 교량의 충돌하중이 미치는 영향에 대한 매개변수연구를 수행하였다.

주기적인 Cellular 유동장내에 놓여진 비대칭형 입자의 운동에 관한 연구 (Motions of Asymmetric Particles in A Periodic Cellular Flowfield)

  • 김윤제
    • 유변학
    • /
    • 제8권2호
    • /
    • pp.139-146
    • /
    • 1996
  • 2차원 비점성 Euler 방정식의 Stuart 해에 의하여 주어지는 주기적인 유동장, 소위 말하여 Stuart와동장내에 놓여있는 아주 작은 비대칭성 입자의 다양한 부력 변수에 따른 침 강운동을 고찰하엿다. 본 연구에서 사용한 수치방법은 매질의 준전적인 Stokes 방정식과 입 자의 평형방정식을 연계하였다. 다양한 초기 배향각 또는 형상비에 따른 궤적이나 각변화를 예측하기 위하여 힘과 토오크 관계식을 입자의 운동방정식에 적용하였으며 4차 Runge-Kutta 방정식을 이용하여 입자의 운동을 규명하였다.

  • PDF

2계 선형상미방 경계치문제의 퍼지시스템 해법 (Solution of the boundary value problem for the second order ordinary differential equations by a fuzzy system)

  • 문병수;정종은;황인구;김정수
    • 한국지능시스템학회:학술대회논문집
    • /
    • 한국퍼지및지능시스템학회 2002년도 춘계학술대회 및 임시총회
    • /
    • pp.289-292
    • /
    • 2002
  • 2계 선형 상미분방정식의 경계치 문제는 보통 해를 구하고자 하는 구간의 양 끝점에서 도함수의 값을 임의로 선정한 후 각 점에서 초기치 문제의 해를 구한 다음 적절한 1차 결합을 이용하여 구하게 된다. 이 경우 초기값과 도함수 값을 사용한 반복연산이 수반되며 따라서 오차의 누적이 불가피 하게 된다. 이 논문에서는 이같은 오차의 누적을 피할 뿐 아니라 3차 Spline 함수를 사용함으로써 오차가 O( $h^2$)인 해를 구하는 방법에 대하여 기술한다 두 개의 경계조건과 근사값을 구하고자 하는 점에서의 함수 값을 "If x is $B_{i}$, then f is $C_{i}$"와 같은 Fuzzy Rule들로 변형하고 주어진 미분방정식을 상수 $C_{i}$들의 관계식으로 변형하여 해를 구하였다. 산출된 결과로부터의 보간 연산은 Fuzzy System사용에 의하여 대체되었다. 이상의 방법으로 산출한 해의 근사오차가 O( $h^2$).임을 증명하였으며 3개의 예제에 대한 계산결과를 4계 Runge-Kutta 방법에 의한 해와 비교하여 기술하였다였다였다였다

  • PDF

고속철 속도변화에 대한 PSC박스 교량의 경간길이 별 동적해석 (A Dynamic Analysis of PSC Box Bridge Varying Span Lengths for Increased Speeds of KTX)

  • 오순택;이동준;심영우;윤준관
    • 한국구조물진단유지관리공학회 논문집
    • /
    • 제15권4호
    • /
    • pp.204-211
    • /
    • 2011
  • 고속철 교량의 동적응답을 보다 정밀하게 해석하기 위한 동적해석방법을 개발하였다. 차후 증가될 초고속(450km/h)을 포함하여 고속 주행하는 KTX 동력차에 의한 교량의 동적거동을 면밀한 속도변수분석과 정밀한 해석을 위한 고속철, 교량 그리고 궤도구조물의 상호작용을 포함한 수치모델을 구성하였다. 네 가지 40~25미터 단순지간의 PSC 박스교를 3차원 유한 프레임요소 모델로 개발하였다. 스펙트럼밀도함수로 산출된 궤도불규칙값과 궤도간 상이한 거리차이를 수치모델화 하였다. 고속철차량은 (KTX) 38자유도로 구성하였다. 38자유도 모델은 3방향 변위와 상응하는 회전각을 고려하였다. 동적증폭계수는 다양한 불규칙 궤도, 켐버, 주행속도, 자갈도상과 같은 주행조건에 의해 결정된다. 이와 같은 동적증폭계수를 해석하기위한 Newmark-${\beta}$ 기법과 Runge-Kutta기법을 적용하여 고속철 속도별과 경간별로 면밀하게 비교 분석하였다.

단섬유보강 플라스틱 사풀성형제품의 섬유방향성에 따른 기계적성질 예 (Study for Predicting Mechanical Properties Depending on Fiber Orientation in Injection Molded Short-Fiber-Reinforced Plastic Parts)

  • 권태헌
    • 유변학
    • /
    • 제8권2호
    • /
    • pp.78-91
    • /
    • 1996
  • 단섬유 보강 플라스틱 복합재료에 사출성형에서 섬유배향은 금형 충전 공정 중의 유동장에 의해 결정되고 섬유의 배향 상태는 역으로 유동장에 영향을 미친다. 단섬유에 의 한 추가적인 응력을 포함하는 Dinh과 Armstrong의 이방성 구성방정식을 충전유동과 섬유 배향의 연계해석에 도입하였다. 충전유동의 해석은 새로운 압력 지배방정식과 에너지 방정 식을 유한요소법과 유한차분법을 이용하여 풀고 동시에 2차 배향 텐서의 변화방정식을 4차 Runge-kutta 방법으로 풀었다. 섬유의 배향상태를 구한 후에 일방향성 복합재료의 Halpin-Tsai 식과 배향 평균모델을 도입하여 사풀성형품의 이방성 기계적 성질이 예측되었 다. 직사각형 캐비티에서 수치해석결과를 실험결과와 비교하였다. 섬유배향과 유동과의 상호 연계작용을 특히 게이트 근처에서 섬유배향에 영향을 미치며 수치해석 결과는 벽면 근처에 서 유동방향으로 배향하는 shell층을 과대 예측함을 알수 있었는데 이는 배향 텐서 변화 방 정식의 최종근사에서 기인하는 오차로 판단된다. 수정된 복합최종 근사를 바탕으로 예측된 이방성 기계적 성질이 기존의 복합최종 근사에 기초한 예측보다 실험 결과에 정량적으로 보 다 잘 일치하였다.

  • PDF

마찰을 고려한 노즈 랜딩기어의 비선형 쉬미 해석 (Non-linear Shimmy Analysis of a Nose Landing Gear with Friction)

  • 이미선;배재성;황재혁
    • 한국항공우주학회지
    • /
    • 제39권7호
    • /
    • pp.605-611
    • /
    • 2011
  • 쉬미는 항공기의 이착륙 시 랜딩기어가 주행도중 측방향 및 조향방향의 진동이 발생하는 현상이다. 쉬미 현상은 스트럿의 낮은 강성, 랜딩기어 내부의 유격, 휠의 불균형이나 마모된 부품 등으로 인해 발생하며, 항공기의 안정성을 저하시킨다. 본 연구는 소형항공기의 쉬미 안정성 검토를 위해 수행되었다. 수치해석을 위하여 소형항공기의 전방 착륙장치를 선형시스템으로 모델링하고 상태방정식을 수립하였다. 근궤적 기법을 이용한 주파수 영역 해석과 4차 Runge-Kutta 방법을 이용한 시간영역 해석을 통해 쉬미 현상을 예측하였고 주요 변수의 설계범위를 검토하였다. 현 착륙장치는 와셔의 압축력을 이용하여 조향 방향 마찰을 가함으로써 쉬미현상을 저감하는 기법을 채택하고 있으므로 마찰을 기술함수를 이용하여 선형화시키고 상태방정식에 적용하여 해석을 수행함으로써 쉬미의 발생이 저감되는 결과를 확인하였다.

비선형 2자유도계의 비감쇠 강제진동 연구 (A Sudy on the Undamped Forced Vibration of Nonlinear Two-Degree-of-Freedom Systems)

  • 박철희;박선재;윤영석
    • 대한기계학회논문집
    • /
    • 제12권2호
    • /
    • pp.193-199
    • /
    • 1988
  • 본 연구의 목적은 SCPA방법을 이용한 비선형 2자유도 비감쇠계의 해석을 통하여 응답곡선을 구하고, 그 응답곡선의 분기현상을 규명함에 있다. 결과의 비교를 위하여 4차의 Runge-Kutta방법을 이용한 수치실험을 수행하였다.

단섬유보강 플라스틱재료의 사출성형 충전공정 수치해석 (Numerical Simulation of Injection Molding Filling Process of Short-Fiber-Reinforced Thermoplastics)

  • 권태헌
    • 유변학
    • /
    • 제5권2호
    • /
    • pp.109-124
    • /
    • 1993
  • 본 연구의 목표는 임의의 3차원 사출성형 금형 공간내에서 단섬유 강화 플라스틱의 충전 공정에서의 과도기적 섬유방향성을 예측하는 수치해석 프로그램의 개발에 있다. Hele-Shaw 방정식에 단섬유에 의해서 추가된 응력을 고려한 Dinh-Armstrong의 모델을 도 입함으로써 새로운 충전과정의 압력 지배 방정식이 유도되었다. 새로운 압력지배 방정식은 단섬유에 의한 응력 때문에 몇 개의 새로운 항들을 포함하고 있다. 충전 과정의 해석은 새 로운 압력지배방정식과 에너지 방정식을 유한효소법과 유한 차분법을 이용하여 풀고 동시에 배향텐서(roientation tensor)의 변화 방정식을 4차 Runge-Kutta 방법을 이용하여 풀었다. 단섬유 배향 텐서를 텐서의 변환 법칙을 이용하여 임의의 3차원 금형 공간내의 모든유한요 소의 중심에서 두께방향의 모든유한 차분 격자를 따라 계산하였다. 이러한 방법으로 임의의 3차원 사출성형 금형 공간내에서 비등온 충전유동과 과도기적 3차원 섬유배향상태를 서로의 상호작용을 고려하여 수치 모사하여 다양한 유동 형태에 따른 단섬유 배향 상태의 변화에 대하여 알아보고자 한다.

  • PDF