• 제목/요약/키워드: Implicit numerical integration

검색결과 119건 처리시간 0.023초

전투기 형상의 외부장착물이 꼬리날개에 미치는 영향에 대한 수치적 연구 (Numerical Investigation of The Effect of External Stores on Tail Wing Surfaces of a Generic Fighter Aircraft)

  • 김민재;권오준;김지홍
    • 한국항공우주학회지
    • /
    • 제36권3호
    • /
    • pp.211-219
    • /
    • 2008
  • 본 연구에서는 삼차원 비정렬 비점성 유동 해석 코드를 이용하여 전투기 형상의 외부 장착물이 꼬리 날개에 미치는 공력 간섭효과에 대한 연구를 수행하였다. 수치적 기법으로는 격자점 중심(vertex-centered)에 기초한 유한체적법과 시간 적분을 위한 내재적인 point Gauss-Seidel 반복 계산법을 사용하였다. 해석 코드의 검증을 위해 전투기 형상에 대한 정상 유동 계산을 수행하였으며 결과를 실험 데이터와 비교하였다. 외부 장착물의 후류(wake)를 정확히 포착하기 위해 예상되는 후류 영역에 대해 국부적인 격자 조밀화를 수행하였으며 천음속 영역에서의 비정상 유동 해석을 통해 외부 장착물에서 발생하는 후류가 수평꼬리날개에 미치는 간섭효과를 확인하였다. 공력 간섭효과를 감소시키기 위한 대안으로는 뒷전 플랩 꺽임각(trailing edge flap deflection)을 고려하였으며 이에 대한 정량적인 감소효과를 제시하였다.

비정렬격자계를 이용한 헬리콥터 덮개 꼬리 로터의 제자리 비행 유동 해석 (Detailed Flow Analysis of Helicopter Shrouded Tail Rotor in Hover Using an Unstructured Mesh Flow Solver)

  • 이희동;권오준;강희정;주진
    • 한국항공우주학회지
    • /
    • 제31권5호
    • /
    • pp.1-9
    • /
    • 2003
  • 본 연구에서는 비정렬격자계를 사용하여 덮개 꼬리 로터의 제자리 비행에 대한 압축성, 비정성 유동을 해석하였다. 유동계산을 위한 수치적 기법으로는 셀 중심에 기초한 유한체적법과 내재적인 시간적분법을 사용하였다. 계산은 로터의 한 블레이드에 대해 수행되었으며, 블레이드와 블레이드 사이에는 주기적 경계조건을 설정하였다. 덮개가 없는 로터 형상에 대한 성능은 실험 결과와 잘 일치함을 보였다. 덮개를 포함하는 로터 형상에 대한 계산은 비교된 실험 형상의 불확실성을 고려하여 추력이 일치하는 피치각을 가지는 경우에 대해 수행하였으며, 자세한 유동은 실험결과와 잘 일치함을 확인할 수 있었다. 그 결과로부터 본 방법이 블레이드 끝단간극을 포함하는 복잡한 3차원 덮개 꼬리 로터 형상 해석에 매우 유용하게 사용될 수 있음을 알 수 있었다.

불균일 단면을 갖는 저장력 예인케이블의 동적해석을 위한 수치해석적 연구 (A Study on Numerical Simulation for Dynamic Analysis of Towed Low-Tension Cable with Nonuniform Characteristics)

  • 정동호
    • 한국전산구조공학회논문집
    • /
    • 제16권1호
    • /
    • pp.69-76
    • /
    • 2003
  • 최근 들어 심해역에 대한 개발과 합섬섬유 재질 케이블의 발달로 인하여 저장력 케이블의 사용이 증가되었다. 저장력 케이블은 장력에 의한 복원력이 작기 때문에 대변위가 발생하게 되며, 따라서 기하학적 비선형이 강하게 나타나게 된다. 또한 해양환경에서는 유체 비선형도 작용하게 된다. 본 연구에서는 수치해석적 방법을 통하여 불균일하게 구성된 예인되는 저장력 케이블의 3차원 동적거동 해석을 수행한다. 수치해석에서는 유체 및 기하학적 비선형과 굽힘강성이 고려되며, 유한차분법(음해법)을 적용하여 풀이된다. 비선형 해를 구하기 위해서 뉴톤-랍슨 방법을 사용한다. 대형 행렬을 풀이하기 위하여 불록삼중대각행렬 풀이법이 적용되는데, 이 방법은 일반적인 행렬 풀이법인 가우스-조르단 방법에 비하여 계산시간을 상당히 줄일 수 있었다. 선배열 음탐 케이블에 대한 다양한 예제해석을 수행하였으며, 해석결과는 미국 우즈홀 해양연구소에서 개발된 프로그램 결과와 잘 일치하였다.

지진격리된 원전배관의 지진취약도 분석 (Seismic Fragility Analysis of Base Isolated NPP Piping Systems)

  • 전법규;최형석;함대기;김남식
    • 한국지진공학회논문집
    • /
    • 제19권1호
    • /
    • pp.29-36
    • /
    • 2015
  • Base isolation is considered as a seismic protective system in the design of next generation Nuclear Power Plants (NPPs). If seismic isolation devices are installed in nuclear power plants then the safety under a seismic load of the power plant may be improved. However, with respect to some equipment, seismic risk may increase because displacement may become greater than before the installation of a seismic isolation device. Therefore, it is estimated to be necessary to select equipment in which the seismic risk increases due to an increase in the displacement by the installation of a seismic isolation device, and to perform research on the seismic performance of each piece of equipment. In this study, modified NRC-BNL benchmark models were used for seismic analysis. The numerical models include representations of isolation devices. In order to validate the numerical piping system model and to define the failure mode, a quasi-static loading test was conducted on the piping components before the analysis procedures. The fragility analysis was performed by using the results of the inelastic seismic response analysis. Inelastic seismic response analysis was carried out by using the shell finite element model of a piping system considering internal pressure. The implicit method was used for the direct integration time history analysis. In addition, the collapse load point was used for the failure mode for the fragility analysis.

유한차분법을 이용한 저장력 예인케이블의 비선형 동적해석 (Nonlinear Dynamic Analysis on Low-Tension Towed Cable by Finite Difference Method)

  • 박한일;정동호
    • 대한조선학회논문집
    • /
    • 제39권1호
    • /
    • pp.28-37
    • /
    • 2002
  • 본 연구에서는 저장력 예인케이블의 비선형 동적거동을 수치적으로 해석하였다. 고장력 케이블해석에서는 흔히 무시되는 굽힘강성의 효과가 저장력 케이블에서는 중요한 역할을 하므로 본 연구에서는 이를 고려하였다. 또한 저장력 케이블에서는 대변위가 발생하기 쉬우므로 기하학적인 비선형 및 유체 비선형 효과가 크므로 이를 고려하였다. 저장력 예인케이블에 대한 3차원 비선형 운동방정식을 수립하고 유한차분법을 적용하여 이산화 시켰다. 시간적분에 있어서 안정적인 해를 얻을 수 있는 음해법(implicit method)을 적용하였으며 비선형 해를 구하기 위하여 Newton-Raphson 반복법을 사용하였다. 케이블과 같이 양단경계조건을 갖고 대각선 주변 성분만 있는 행렬식을 계산하는 경우에는 Gauss-Jordan 방법 등과 같이 일반적인 방법 보다 블록삼중대각행렬 풀이법이 계산시간을 상당히 줄일 수 있음을 알 수 있었다. 몇 가지 예제해석을 수행하였으며 실해역 실험결과에 의해 이미 검증되어 있는 케이블 해석프로그램인 WHOI Cable 프로그램의 해석결과와 비교 검토한 결과 서로 잘 일치함을 알 수 있었다.

수평가진에 의한 계류라인의 장력응답 해석 (Analysis on Tension Response of Mooring Line by Lateral Excitation)

  • 정동호;김현주;문덕수;박한일;최학선
    • 한국해양환경ㆍ에너지학회지
    • /
    • 제7권4호
    • /
    • pp.185-191
    • /
    • 2004
  • 계측주라인은 부체를 일정 위치에 유지시키기 위한 구조요소로 이용된다 본 연구에서는 부소파제의 계류라인으로 사용되는 체인에 대한 동적거동 특성을 파악하기 위해서 체인에 대한 구조해석을 수행하였다. 먼저 체인에 대한 3차원 동적방정식을 유도하는데, 장력이 0부근인 영역에서 복원력을 표현하기 위하여 굽힘강성 성분을 포함시켰다. 수치적인 방법 통하여 3차원 동적 방정식을 해석하였다. 유한차분법을 적용하였는데, 조건에 상관없이 안정적인 음해법과 함께 비선형 해법인 뉴톤-랍슨 반복법을 사용하였다. 본 연구를 통하여 가진 주기에 따른 계류라인의 위치에 따른 장력응답 경향을 파악할 수 있었으며, 본 연구의 해석결과는 참고문헌의 실험결과와 그 경향이 잘 일치하였다. 본 연구결과는 향후 부소파제를 위한 계류라인 설계에 잘 이용될 수 있을 것이다.

  • PDF

그물망과 대량입자의 멀티 스케일 접촉해석 (Multi-Scale Contact Analysis Between Net and Numerous Particles)

  • 전철웅;손정현
    • 대한기계학회논문집A
    • /
    • 제38권1호
    • /
    • pp.17-23
    • /
    • 2014
  • 그래픽 처리장치(GPU)는 병렬적인 정보를 포함하는 문제를 해결하는데 이상적이다. 본 연구에서는 GPU 는 입자동역학과 함께 다물체 동역학 시뮬레이션을 효율적으로 수행하기 위해 사용되었다. 수치계산을 위해서 HHT 암시적 적분 알고리즘이 사용되었다. 입자들 사이의 접촉을 판별하기 위해서 공간 분할 알고리즘과 입자 거동 해석법으로 이산 요소법(DEM)이 사용되었다. 개발된 다물체 동역학 프로그램은 해는 ADAMS 프로그램의 결과와 비교 검증하였다. CPU 기반의 순차해석 프로그램과 GPU 기반 병렬 프로그램은 입자의 수에 따른 수치계산 효율성을 알아보기 위해 서로 비교되었으며, 입자의 수가 많아질수록 계산시간은 단축되었다. 본 예제에서 입자의 수가 1,300 개일 때, 순차 해석 프로그램보다 병렬 프로그램이 약 5 배 가량 빠른 계산 속도를 보였다.

폐쇄형 풍동 시험부내의 원형 실린더 유동에 대한 비정상 벽면효과 연구 (Unsteady Wall Interference Effect on Flows around a Circular Cylinder in Closed Test-Section Wind Tunnels)

  • 강승희;권오준;홍승규
    • 한국항공우주학회지
    • /
    • 제33권7호
    • /
    • pp.1-8
    • /
    • 2005
  • 풍동 시험부 비정상 벽면효과에 대한 연구를 위해 폐쇄형 시험부내의 원형 실린더 주위 유동장에 대한 수치적 연구를 수행하였다. 수치기법은 Roe의 flux-difference-splitting을 사용한 격자점 중심 유한체적법과 이중시간 전진 기법을 사용하는 내재적 시간적분법을 사용하였다. 계산 결과 폐쇄형 시험부에는 실린더 주위 비정상 유동장에 압력구배를 강화시켜 실린더의 양력 및 항력의 진폭을 크게 하고, 실린더 뒷전에서의 기저압력을 작게 하여 항력을 증가시키는 벽면효과가 있음을 확인하였다. 또한, 이러한 시험부 벽면은 실린더 와류 shedding 주파수를 커지게 하는 효과가 있다. 시험부 벽면에서의 압력은 벽면효과가 포함된 shedding 주파수를 기본으로 하는 고조파 현상을 보인다.

평판형 전개판의 3차원 운동 모델링 (Modeling of flat otter boards motion in three dimensional space)

  • 최무열;이춘우;이건호
    • 수산해양기술연구
    • /
    • 제43권1호
    • /
    • pp.49-61
    • /
    • 2007
  • Otter boards in the trawl are the one of essential equipments for the net mouth to be spread to the horizontal direction. Its performance should be considered in the light of the spreading force to the drag and the stability of towing in the water. Up to the present, studies of the otter boards have focused mainly on the drag and lift force, but not on the stability of otter boards movement in 3 dimensional space. In this study, the otter board is regarded as a rigid body, which has six degrees of freedom motion in three dimensional coordinate system. The forces acting on the otter boards are the underwater weight, the resistance of drag and spread forces and the tension on the warps and otter pendants. The equations of forces were derived and substituted into the governing equations of 6 degrees of freedom motion, then the second order of differential equations to the otter boards were established. For the stable numerical integration of this system, Backward Euler one of implicit methods was used. From the results of the numerical calculation, graphic simulation was carried out. The simulations were conducted for 3 types of otter boards having same area with different aspect ratio(${\lambda}=0.5,\;1.0,\;1.5$). The tested gear was mid-water trawl and the towing speed was 4k't. The length of warp was 350m and all conditions were same to each otter board. The results of this study are like this; First, the otter boards of ${\lambda}=1.0$ showed the longest spread distance, and the ${\lambda}=0.5$ showed the shorted spread distance. Second, the otter boards of ${\lambda}=1.0$ and 1.5 showed the upright at the towing speed of 4k't, but the one of ${\lambda}=0.5$ heeled outside. Third, the yawing angles of three otter boards were similar after 100 seconds with the small oscillation. Fourth, it was revealed that the net height and width are affected by the characteristics of otter boards such as the lift coefficient.