Advanced
Vibration Analysis of Rotating Pre-twisted Inward Beams with a Concentrated Mass
Vibration Analysis of Rotating Pre-twisted Inward Beams with a Concentrated Mass
Transactions of the Korean Society for Noise and Vibration Engineering. 2015. Jun, 25(6): 384-390
Copyright © 2015, The Korean Society for Noise and Vibration Engineering
  • Received : November 27, 2014
  • Accepted : June 09, 2015
  • Published : June 20, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
건 호 이
Member, School of Mechanical Engineering, Hanyang University
홍 희 유
Corresponding Author ; Member, School of Mechanical Engineering, Hanyang University E-mail :hhyoo@hanyang.ac.kr

Abstract
The vibration analysis of rotating inward beams considering the pre-twisted is presented based on Euler-Bernoulli beam theory. The frequency equations, are calculated using hybrid deformation variable modeling along with the Rayleigh-Ritz assumed mode methods. In this study, resulting system of ordinary differential equations shows the effects of angular speed, and Young's modulus ratio. It is believed that the results will be a reference with which other researchers and commercial FE analysis program, ANSYS can compare their result.
Keywords
1. 서 론
회전중심방향 자유단 외팔보의 구조물은 터보홴, 믹싱 터보홴 엔진에 활용되고 있다. 원심력을 받는 이러한 구조물의 진동에 관한 연구는 현재까지 활발히 수행 되어왔다 (1~3) . 구조물의 회전으로 인해 외팔보의 강성이 감소하여 진동특성 또한 변하게 된다. 따라서 회전으로 인한 진동특성 변화를 정확히 예측하기 위한 연구의 중요성은 강조되고 있어 많은 연구가 수행되고 있다.
회전중심방향 자유단 외팔보 구조물의 진동특성에 관한 연구를 설명하기 앞서 일반적인 회전 구조물의 연구는 1920년대 초부터 시작되었다. 초기의 연구들에서는 현재까지도 그 간편성 때문에 널리 쓰이고 있는 Southwell (4) 방정식이라 불리는 고유진동수 예측 식이 소개되고 발전 개선되었다. 1970년도 이후에는 컴퓨터를 이용한 수치해석적 방법을 통하여 고유진동수 및 모드형상을 구하는 연구결과들이 등장하게 되었으며 또한 참고문헌과 같이 진동방정식을 해석적으로 직접 풀이하는 방법도 등장하게 되었다 (5~7) .
이 논문에서는 일반적인 회전 구조물의 연구와 다르게 회전중심축방향 자유단 외팔보의 진동특성을 주제로 연구를 수행하였으며, 이에 관련된 논문으로 Mostaghel, N. et al (8) 의 연구가 있다. 그 후 Fox, C. H. J. et al (9) 는 Mostaghel, N. et al (8) 의 연구 결과를 발전시켜 초기 장착각을 고려하여 회전중심축방향 자유단 외팔보의 진동특성 변화에 관한 연구를 수행했다. 그러나 이 연구들에서는 회전중심축방향 자유단 외팔보에서 흔히 나타나는 집중질량과 초기 비틀림을 갖는 구조형태를 고려하지 않았다. 반면 이 논문에서는 초기 비틀림을 갖는 회전 중심축 방향 자유단 외팔보가 임의의 위치에 집중질량을 갖는 경우의 진동특성 연구를 수행하는데 목적을 두었다. 이 연구는 참고문헌 (10) 에서 제시된 일반 회전 구조물에 집중질량과 초기 비틀림을 갖는 회전 외팔보의 모델링에 근거 하였으며 여기에 참고문헌 (9) 에 제시된 회전중심방향 자유단 외팔보의 모델링을 종합하여 외팔보의 길이에 대한 축반경비, 회전각속도, 비틀림각, 집중질량의 크기 및 위치 등에 의한 진동특성 변화를 종합적으로 조사하였다는 점에 고유성을 갖는다. 운동방정식을 얻기 위해 복합 변형 모델링과 Kane의 방법을 이용하여 유도하였고 상용 프로그램인 ANSYS를 통해 얻은 결과와 유도된 방정식을 통해 얻은 결과를 비교하여 제안된 모델의 정확성을 검증하였다.
이 논문의 구성으로는 2장에서 동적 모델링에 의한 운동방정식을 유도하며 3장에서는 유도된 운동방정식에 근거하여 정확성을 검증하기 위해 회전중심 방향 자유단 외팔보의 진동해석결과를 결과를 참고문헌 (11) 과 비교하였다. 기존의 결과에서 비틀림을 추가하여 얻어낸 해석 결과를 상용 프로그램 ANSYS를 통해 비교하여 검증하였다. 마지막 장에서는 이 논문의 결론을 수록하였다.
2. 운동방정식
- 2.1 운동방정식
Fig. 1 은 반경 r 의 디스크에 고정되어 평면에서 â3 방향으로 ω 만큼 회전 운동하는 외팔보를 보여준다. 디스크 축에 부착된 서로 수직인 단위벡터들이 나타나 있으며, L 은 빔의 길이를 이고 d 는 집중질량의 위치이다.
PPT Slide
Lager Image
Configuration of a pre-twisted inward beam with a concentrated mass fixed to a rotating rigid ring
Fig. 2 는 외팔보의 회전에 의한 변형 시 형상을 나타내며 탄성변위의 의미와 고정단에 고정되어있는 외팔보의 기준 좌표계를 나타낸다. x 는 외팔보와 디스크가 고정점 O 로부터 변형 전 임의 질점까지 거리를, s 는 외팔보 중립축의 임의 점에서 보가 인장된 길이를 나타낸다. P 0 는 변형 후 x + s 의 길이를 갖는 P의 위치로 이동된다. 모드방법을 이용해 상미분 운동방정식으로 구하기 위해 s , u 2 , u 3 를 다음과 같이 근사화 한다.
PPT Slide
Lager Image
Inward beam with a concentrated mass attached to a rotating rigid ring
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, ϕ 1j , ϕ 2j , ϕ 3j s , u 2 , u 3 를 위한 가상모드이다. q 1j , q 2j , q 6j 는 일반좌표 그리고 μ 1 , μ 2 , μ 3 는 각각 사용된 일반좌표의 개수이다. 반경 r 인 강체축 A 의 각속도와 외팔보의 고정점 O 의 속도는 다음과 같이 사용한다.
PPT Slide
Lager Image
PPT Slide
Lager Image
운동방정식을 유도하기 위해 먼저 P 점의 속도
PPT Slide
Lager Image
는 다음과 같이 구할 수 있다.
PPT Slide
Lager Image
식 (6)에서 u 1 은 근사화를 위해 최종적으로 s u 2 로 나타내야 하며, 이를 위하여 다음의 근사화된 관계식이 유도되어(참고문헌 (12) 참조) 사용된다.
PPT Slide
Lager Image
속도
PPT Slide
Lager Image
를 이용하면 시스템의 운동방정식은 다음 식으로부터 유도될 수 있다
PPT Slide
Lager Image
여기서, L ρ 는 블레이드의 길이와 단위길이당 질량을 나타내며, U 는 탄성에너지로서 다음과 같이 나타낼 수 있다.
PPT Slide
Lager Image
여기서, E 는 블레이드의 영의계수, A 는 단면적이다. 임의 단면에서 블레이드의 2차면적 모멘트와 프로덕트를 나타내는 I 2 , I 3 , I 23 를 주면적 모멘트
PPT Slide
Lager Image
PPT Slide
Lager Image
로 나타내면
PPT Slide
Lager Image
여기서, θ 는 블레이드의 임의의 단면에서 비틀림각을 나타내며, 따라서 블레이드의 고정단에서 θ = 0 , 블레이드의 자유단에서 θ = θ 0 이다. 이상의 식들을 이용하여 운동방정식을 구하면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
블레이드의 임의 지점에 집중질량이 추가로 붙어 있은 것을 모델링하기 위해 블레이드의 단위길이당 질량을 나타내며 집중질량의 임의의 지점 x = d 에서 충격함수(Dirac’s delta function)로 가정하고, 방정식에서 ρ 와 관련된 항들을 다시 구하면 다음과 같다.
PPT Slide
Lager Image
여기서, m은 임의의 지점에 위치한 집중질량의 크기를 나타낸다. 블레이드의 인장방향 고유진동수와 굽힘방향 고유진동수는 그 차이가 크므로 이들의 연성에 의한 효과는 무시될 수 있으므로 이 논문에서는 그 효과를 무시하고 수치해석을 하였다. 연성 효과를 무시할 때 집중질량을 고려한 항들의 식 (11~13)에 대입하여 운동방정식을 구하면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
- 2.2 무차원 운동방정식 유도
유도된 운동방정식, 식 (16)과 (17)을 무차원화 하기 위해 다음과 같이 시간과 공간 그리고 일반좌표 변수들을 무차원화 한다.
PPT Slide
Lager Image
식 (18)에서 T 는 다음과 같이 정의되는 변수이다.
PPT Slide
Lager Image
무차원 변수로 정의된 식 (18)를 이용하여 식 (10)을 다음과 같이 나타낼 수 있다.
PPT Slide
Lager Image
이상의 관계식들을 이용하여 식 (16)과 (17)을 무차원화된 방정식들로 만들면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
φ aj ϕ aj 와 같은 함수값을 가지나 ξ 의 함수이며
PPT Slide
Lager Image
ϑ aj τ 에 관한 2차 미분을 나타낸다. 식 (21)와 (22)에서
PPT Slide
Lager Image
PPT Slide
Lager Image
로 구성되어 있는 column 벡터를
PPT Slide
Lager Image
라 한다. 이 벡터를 시간에 대한 조화함수로 가정하면 다음과 같다.
PPT Slide
Lager Image
여기서 ω 는 고유진동수이며
PPT Slide
Lager Image
는 모드형상을 나타내는 모드벡터이다. 식 (24)를 이용하여 식 (21)와 (22)은 다음과 같이 표현한다.
PPT Slide
Lager Image
여기서 M K 는 대칭 정방행렬들로 그 요소들은 다음과 같이 표시된다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
3. 수치해석
이 장에서는 앞에서 유도된 무차원 진동방정식을 사용하여 수치 해석을 수행하였다. 수치 해석을 위해 굽힘방향 변위의 가상모드(assumed modes)들로는 회전중심방향 자유단 외팔보의 굽힘 진동모드들을 사용하였으며 10개의 모드를 이용하였다.
Fig. 3 은 회전중심방향 자유단 외팔보의 유도된 운동방정식의 정확성과 검증하기 위해 회전경우에 대해서 기존결과 참고문헌 (11) 과 비교하였다. 최대 오차율 0.3 % 오차 결과 비교는 이 연구에서 제시된 모델링 방법을 검증하고 신뢰성을 보여준다고 판단된다.
PPT Slide
Lager Image
Comparison of natural frequencies obtained by proposed EOM and reference
Fig. 4 는 집중질량과 초기 비틀림을 갖는 회전중심방향 자유단 외팔보가 집중질량 비의 변화에 따른 고유진동수를 나타냈다. 수치해석 결과는 ANSYS와 비교하여 검증하였다. ANSYS의 유한요소 모델은 beam 188 요소를 사용하였고 집중질량은 mass 21요소를 사용하였다. 집중질량 비는 고유진동수를 0에서 2사이에 급격히 변화시키며 그 이후는 완만한 변화를 보여주고 있다.
PPT Slide
Lager Image
Concentrated mass ratio effect on the natural frequencies( β =1, γ = 6, κ = 0.25, θ0 = 90, δ = 1)
Fig. 5 에서는 집중질량의 임의의 위치가 변할 때의 고유진동수를 나타낸다. 수치해석 결과는 ANSYS와 비교하여 검증하였다. 무차원 변수 δ 가 1로 가까워질수록 자유단으로 가까워지며 1이었을 때는 자유단 끝에서의 집중질량의 위치를 나타낸다. 무차원 집중 질량의 위치가 0.3정도까지는 고유진동수의 작은 변화만 초래하나 그 이후는 급격한 변화를 보여주고 있다.
PPT Slide
Lager Image
Concentrated mass location effect on the natural frequencies( α =1, γ = 6, κ = 0.25, θ0 = 90, δ = 1)
4. 결 론
이 연구에서는 임의의 지점에 집중질량과 초기 비틀림을 갖는 회전중심방향 자유단 외팔보의 진동특성을 예측하기 위한 모델링 방법을 제시하였으며 보다 일반적인 결과를 얻기 위해 무차원 운동 방정식의 형태로 전환하였고 이에 근거한 수치해석을 수행하였다. 그 해석결과의 정확성은 수치 해석 결과를 참고문헌과 ANSYS의 무차원 고유진동수를 비교하여 규명하였다. 수치해석 결과는 회전 각속도, 집중질량의 위치 비, 집중질량의 무게, 초기 비틀림, 빔의 길이와 단면적의 비의 변화를 예측하였다. 그리고 이 연구에서 사용된 모델링 방법은 터보홴 또는 가정용 풍력 발전기의 회전익 분야에 사용되고 있으므로 논문에서 제시된 방법이 위의 조건들을 고려한 외팔보의 구조설계 시에 유용하게 사용될 수 있으리라 판단된다.
A part of this paper was presented at the KSNVE 2014 Annual Autumn Conference
Acknowledgements
이 연구는 2011년도 산업통상자원부 재원으로 한국에너지기술평가원(KETEP)의 지원을 받아 수행한 연구 과제입니다(No. 20111510100050).
BIO
Hong Hee Yoo received his B.S. and M.S. degrees in the Department of Mechanical Design in Seoul National University in 1980 and 1982. He received his Ph.D. degree in the Department of Mechanical Engineering and Applied Mechanics in the University of Michigan at Ann Arbor in 1989. He is a professor in the Department of Mechanical Engineering at Hanyang University, Seoul, Korea. His research interests include multi-body dynamics, structural vibration, and statistical uncertainty analysis in mechanics.
Gun Ho Lee received his B.S. degrees in the Department of Mechanical Engineering in Kyungpook National University, Daegu, Korea in 2013. He is working as a M.S. candidate in the Department of Mechanical Engineering in Hanyang University, Seoul, Korea. His research interests include structural vibration.
References
Fox C. H. J. 1985 The Free Vibration of Compact Rotating Radial Cantilevers Journal of Sound and Vibration 98 (3) 325 - 336    DOI : 10.1016/0022-460X(85)90279-2
Naguleswaran S. 1994 Lateral Vibration of a Centrifugally Tensioned Uniform Euler-Bernoulli Beam Journal of Sound Vibration 176 (5) 613 - 624    DOI : 10.1006/jsvi.1994.1402
Turhan O. , Bulut G. 2005 On Nonlinear Vibrations of a Rotating Beam Journal of Sound and Vibration 322 (1-2) 314 - 335
Southwell R. , Gough F. 1921 The Free Transverse Vibration of Airscrew Blades
Schihansl M. 1958 Bending Frequency of a Rotating Cantilever Beam Journal of Applied Mechanics 25 28 - 30
Yoo H. H. , Ryan R. R. , Scott R. A. 1995 Dynamics of Flexible Beams Undergoing Overall Motions Journal of Sound and Vibration 181 (2) 261 - 278    DOI : 10.1006/jsvi.1995.0139
Lee J. H. , Yoo H. H. 2010 Vibration Analysis of Rotating Blades with the Cross Section Taper Considering the Pre-twist Angle and the Setting Angle Transactions of the Korean Society for Noise and Vibration Engineering 20 (1) 10 - 21    DOI : 10.5050/KSNVE.2010.20.1.010
Mostaghel N. , Tadjbakhsh I. 1973 Buckling of Rotating Rods and Plates International Journal of Mechanical Sciences 15 (6) 429 - 434    DOI : 10.1016/0020-7403(73)90026-X
Fox. C. H. J. 1979 The Natural Frequencies of a Thin Rotating Cantilever with Offset Root Journal of Sound Vibration 65 (2) 151 - 158    DOI : 10.1016/0022-460X(79)90509-1
Yoo H. , Shin S. 1998 Vibration Analysis of Rotating Cantilever Beams Journal of Sound and Vibration 212 (5) 807 - 828    DOI : 10.1006/jsvi.1997.1469
Naguleswaran S. 1997 Out-of-plane Vibration of a Uniform Euler-Bernoulli Beam Attached to the Inside of a Rotating Rim Journal of Sound Vibration 200 (1) 63 - 81    DOI : 10.1006/jsvi.1996.0667
Kane T. R. , Levinson D. A. 1985 Dynamics: Theory and Applications McGraw-Hill Book Co. New York