선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 본 논문에서는 에어포일의 단면 형상과 재료변화에 따른 고유진동수의 변화 경향을 파악하고자 CR-보를 이용한 에어포일 단면 형상을 가진 외팔보에 대하여 재료, spar 유무, 단일재료·복합재료, 에어포일 최대두께 변화에 따른 자유단의 끝단 진폭과 고유진동수를 비교하고 그에 대한 경향 및 차이점을 분석하였다.
- 본 연구에서는 에어포일의 단면형상과 재료변화에 따른 고유진동수의 변화 경향을 파악하고자 CR-보를 이용한 에어포일 단면 형상을 가진 외팔보에 대하여, 재료, spar 유무, 단일재료·복합재료, 에어포일 최대두께 변화에 따른 자유단의 끝단 진폭과 고유진동수를 비교하고, 그에 대한 경향 및 차이점을 분석하였다.
가설 설정
- 에어포일 단면 형상 변화에 따른 끝단 진폭 및 고유진동수 변화의 경향을 파악하기 위해 먼저 테이퍼가 없는 항공기 날개를 외팔보로 가정하고 시위의 길이를 1m 기준으로 날개 크기 비율을 조정하였다. 보에 작용하는 하중은 B-737 기준 25,000ft 고도로 순항하는 날개에 작용하는 양력(분포하중)으로 가정하였다. 또한 CR-보 해석에 필요한 단면의 영률(E), 단면적(A), 관성모멘트(I) 등 물성치는 VABS 프로그램을 통해 산출하였다.
제안 방법
- Fig. 2와 같이 해석 대상물에 대한 질량행렬, 강성행렬을 구성하여 고유치 해석을 수행하였다. 해석 대상물의 고유치 해석 결과는 Table 2와 같다.
- 이 장에서는 CR 기법을 이용한 외팔보 모델링과 상용프로그램 ANSYS와의 주파수 영역에서의 결과 비교를 통하여 CR 외팔보 모델의 신뢰성을 확보한다. 다음으로 에어포일 단면 형상, 재료 변화에 따른 에어포일 끝단 진폭과 고유 진동수 변화에 대한 경향과 차이점을 분석하였다.
- 보에 작용하는 하중은 B-737 기준 25,000ft 고도로 순항하는 날개에 작용하는 양력(분포하중)으로 가정하였다. 또한 CR-보 해석에 필요한 단면의 영률(E), 단면적(A), 관성모멘트(I) 등 물성치는 VABS 프로그램을 통해 산출하였다. Fig.
- 4는 sine sweep test 결과이다. 시작 가진주파수 3Hz부터 종단 가진주파수 60Hz까지 3의 배수로 가진주파수 변화를 주었으며, amplitude가 증가하는 특정구간에 대해 더 작은 크기의 가진주파수 변화를 주며 공진이 일어나는 지점을 찾아 고유진동수를 도출하였다.
- 6은 CR-보 해석을 위한 에어포일 단면 형상이다. 에어포일 단면 형상 변화에 따른 끝단 진폭 및 고유진동수 변화의 경향을 파악하기 위해 먼저 테이퍼가 없는 항공기 날개를 외팔보로 가정하고 시위의 길이를 1m 기준으로 날개 크기 비율을 조정하였다. 보에 작용하는 하중은 B-737 기준 25,000ft 고도로 순항하는 날개에 작용하는 양력(분포하중)으로 가정하였다.
- 에어포일 시위(1m)의 12~30% 범위에서 최대두께만 변화시킨 NACA0012, NACA0015, NACA0018, NACA 0024, NACA0030 형상에 기존 알루미늄 에어포일과 내부가 비어있고, 얇은 외벽 두께에 Spar를 가진 알루미늄 에어포일을 대상으로 끝단 진폭, 고유진동수로 나누어 결과를 비교하였다. Fig.
- 에어포일 시위(1m)의 12~30% 범위에서 최대두께만 변화시킨 NACA0012, NACA0015, NACA0018, NACA0024, NACA0030 형상에 spar를 가진 알루미늄 에어포일과 경량화된 등방성 재질의 Al 6061T6, Eglass 120과 이방성 재질의 Eglass 456, S2 Glass로 이루어진 복합재료 에어포일을 대상으로 끝단 진폭, 고유주파수로 나누어 결과를 비교하였다. Table 7과 Fig.
- 에어포일 시위(1m)의 12~30% 범위에서 최대두께만 변화시킨 NACA0012, NACA0015, NACA0018, NACA0024, NACA0030 형상에 각각 Isopink와 알루미늄 재료로 내부가 꽉 채워진(solid) 에어포일을 대상으로 CR-보 해석을 진행하여 끝단 진폭, 고유진동수로 나누어 결과를 비교하였다. Fig.
- 5는 ANSYS 모드 해석 결과이다. 해석 결과 1차 모드는 1차 bending에 일치하는 것을 확인할 수 있었으며, 이에 따라 1차 bending 모드에 대한 고유진동수를 비교하기로 결정하였다.
- 해석에 사용한 CR-보 해석 프로그램의 신뢰성을 확인하기 위해 특정 대상물에 대한 검증작업을 수행하였다. 검증작업은 시스템 인식, ANSYS 모드 해석, CR-보 고유치 해석을 통해 수행하였다.
데이터처리
- 보에 작용하는 하중은 B-737 기준 25,000ft 고도로 순항하는 날개에 작용하는 양력(분포하중)으로 가정하였다. 또한 CR-보 해석에 필요한 단면의 영률(E), 단면적(A), 관성모멘트(I) 등 물성치는 VABS 프로그램을 통해 산출하였다. Fig.
- 2는 상용프로그램 ANSYS를 통한 해석 대상물이다. CR-보에서 가정한 모델에 대한 결과를 상용프로그램인 ANSYS의 결과와 비교 분석하여 신뢰성을 확보하였다. Table 1은 비교 모델의 물성치와 외력의 크기를 나타낸다.
- 해석에 사용한 CR-보 해석 프로그램의 신뢰성을 확인하기 위해 특정 대상물에 대한 검증작업을 수행하였다. 검증작업은 시스템 인식, ANSYS 모드 해석, CR-보 고유치 해석을 통해 수행하였다.
- CR-보 해석에 필요한 단면의 영률, 단면적, 관성모멘트 등물성치는 VABS 프로그램을 통해 산출하였다. 계산결과의 신뢰성을 얻기 위해 상용프로그램 ANSYS와 CR-보 해석 결과의 끝단 진폭, 고유진동수를 비교하였고 오차(0.31%)의 크기가 매우 작으므로 신뢰성을 확보하였다고 판단하였다.
- 3은 시스템 인식과 ANSYS 모드 해석을 이용한 주파수 영역에서의 검증 흐름도이다. 시스템 인식은 sine 형태 가진 주파수 변화를 주어 CR-보를 통해 얻은 시간에 따른 변위 데이터를 MATLAB 신속 Fourier 변환을 이용하여 고유진동수를 도출하는 sine sweep 결과와 상용프로그램인 ANSYS 모드 해석을 통해 얻은 고유진동수 값과 비교하였다.
성능/효과
- 해석 대상물의 고유치 해석 결과는 Table 2와 같다. 1-3차 모드의 고유진동수의 값이 시스템 인식, ANSYS 모드 해석 결과의 값과 1% 미만의 오차로 일치함을 확인할 수 있었다.
- 위의 3가지 검증작업을 통해 얻은 해석 대상물의 고유 진동수 비교 결과는 Table 3과 같다. ANSYS 모드 해석, CR-보 해석, 고유치 해석을 통해 각각 얻은 고유진동수의 절대적 차이가 모두 1% 미만이므로 CR-보 해석의 신뢰성을 확보하였다고 판단하였다.
- 16은 단일재료(알루미늄)로 이루어진 spar를 가진 에어포일과 복합재료로 이루어진 spar를 가진 에어포일의 고유진동수 비교 결과이다. 그 결과 spar를 가진 알루미늄과 복합재료 에어포일 고유진동수는 평균 0.7%의 차이로 유사한 변화를 나타내었다. 그래프상의 전반적인 고유진동수의 증가는 에어포일의 최대두께 증가에 따른 2차 관성모멘트/단면적의 증가 때문이며, 단일재료로 이루어진 spar를 가진 에어포일 형상과 복합재료로 이루어진 spar를 가진 에어포일 형상의 2차 관성모멘트/단면적 크기가 유사한 값을 가지기 때문에 고유진동수의 변화가 유사한 경향을 나타낸다는 사실을 확인할 수 있었다.
- 13은 알루미늄으로 내부가 채워진 에어포일과 Spar를 고려한 에어포일의 고유진동수 비교 결과이다. 그 결과 알루미늄 대비 spar를 가진 알루미늄 고유진동수가 평균 2.74Hz(약 33.8%) 높게 나타났다. 이를 통해 최대두께가 증가할수록 2차 관성모멘트/단면적의 크기가 증가하여 전반적인 고유진동수의 증가 경향을 가져온다는 사실을 확인할 수 있었다.
- 7%의 차이로 유사한 변화를 나타내었다. 그래프상의 전반적인 고유진동수의 증가는 에어포일의 최대두께 증가에 따른 2차 관성모멘트/단면적의 증가 때문이며, 단일재료로 이루어진 spar를 가진 에어포일 형상과 복합재료로 이루어진 spar를 가진 에어포일 형상의 2차 관성모멘트/단면적 크기가 유사한 값을 가지기 때문에 고유진동수의 변화가 유사한 경향을 나타낸다는 사실을 확인할 수 있었다. Table 8은 복합재료에 따른 최대변위와 평균 고유진동수 비교 결과를 정리한 결과이다.
- 내부가 빈 얇은 두께의 외벽과 spar를 가진 에어포일은 내부가 재료로 꽉찬 에어포일 대비 끝단 최대진폭이 작고, 고유진동수는 평균 2.74Hz(약 33.8%) 높게 나타났다.
- 이를 통해 최대두께가 증가할수록 2차 관성모멘트/단면적의 크기가 증가하여 전반적인 고유진동수의 증가 경향을 가져온다는 사실을 확인할 수 있었다. 또한, spar 유무에 따른 고유진동수 차이가 발생한 원인으로는 재료의 영률과 밀도의 값이 일정할 때 고유진동수의 크기는 2차 관성모멘트/단면적의 값에 비례한다는 사실을 확인할 수 있었다. Table 6은 spar 유무에 따른 최대변위 및 고유진동수 비교 결과를 정리한 결과이다.
- 또한, spar를 가진 에어포일 단면 형상에 재료를 경량화된 복합재료로 구성된 에어포일의 고유주파수와 비교한 결과, 기존 spar를 가진 알루미늄 에어포일과 유사한 고유진동수 변화를 나타내었다.
- 5%) 높게 나타났다. 에어포일 최대두께 변화를 통해 에어포일의 2차 관성모멘트/단면적이 증가할수록 고유진동수가 증가한다는 사실을 확인할 수 있었고, 재료 비교를 통해 같은 단면적이라도 영률과 밀도 값과 고유진동수의 크기가 비례한다는 사실을 확인할 수 있었다.
- 이를 바탕으로 에어포일의 고유진동수 변화는 2차 관성 모멘트/단면적, 밀도, 영률의 변화와 종속적임을 확인할 수 있었다.
- Table 9는 각각의 해석에 대한 결과를 평균 고유진동수를 기준으로 요약한 결과이다. 이를 통하여 동일 단면적, 2차 관성모멘트 조건에서 고유진동수는 재료의 영률과 밀도에 비례하고, 동일 재료의 경우 고유진동수는 2차 관성모멘트/단면적에 비례한다는 사실을 확인할 수 있었다.
- 123m 크게 나타났다. 이를 통해 에어포일을 구성하는 재료의 부피가 작을수록 끝단 진폭이 크게 나타남을 확인할 수 있었다.
- 194m 더 크게 나타났다. 이를 통해 재료의 영률과 밀도가 작을수록 끝단 진폭이 크게 나타남을 확인할 수 있었다.
- 8%) 높게 나타났다. 이를 통해 최대두께가 증가할수록 2차 관성모멘트/단면적의 크기가 증가하여 전반적인 고유진동수의 증가 경향을 가져온다는 사실을 확인할 수 있었다. 또한, spar 유무에 따른 고유진동수 차이가 발생한 원인으로는 재료의 영률과 밀도의 값이 일정할 때 고유진동수의 크기는 2차 관성모멘트/단면적의 값에 비례한다는 사실을 확인할 수 있었다.
- 재료 변화에 따른 에어포일의 끝단 진폭과 고유진동수를 비교한 결과, 영률과 밀도가 큰 재료로 구성된 에어포일 끝단 최대진폭이 작고 고유진동수가 높게 나타났다.
- 10은 에어포일 재료 변화에 따른 고유진동수의 비교 결과이다. 최대두께가 증가할수록 알루미늄의 고유주파수가 더 큰 폭으로 증가하였으며, Isopink 대비 알루미늄의 고유진동수가 평균 3.57Hz(약 66.5%) 높게 나타났다. 에어포일 최대두께 변화를 통해 에어포일의 2차 관성모멘트/단면적이 증가할수록 고유진동수가 증가한다는 사실을 확인할 수 있었고, 재료 비교를 통해 같은 단면적이라도 영률과 밀도 값과 고유진동수의 크기가 비례한다는 사실을 확인할 수 있었다.
후속연구
- 이전 연구들은 받음각 변화, 재료 변화에 따른 에어포일의 고유진동수 연구가 일부 존재할 뿐, 에어포일의 형상 변화가 고유 진동수에 미치는 영향에 대한 연구는 미비하다. 따라서 실제 항공기의 초기 에어포일 설계 단계에서 단면 형상 변화에 따른 진동특성을 고려한다면, 보다 안정적인 설계 절차를 구축할 수 있을 것으로 예상된다.
- 본 논문에서 제시한 에어포일의 단면형상 변화와 재료 변화에 따른 고유진동수의 변화 경향은 초기 항공기 에어포일 설계 단계에서 안정성을 고려할 때, 근본적인 정보로 활용될 수 있을 것이다.
- 향후 복합재료의 적층각도 및 두께 변화를 고려하고, 1m로 가정한 시위의 크기를 실제 크기로 확대하여 비율을 조정하여 유동해석 결과를 함께 고려한다면 보다 정확한 에어포일 형상 설계 방향을 제시할 수 있을 것으로 예상된다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.