Advanced
Vibration and Stability Analysis of Composite Spinning Shafts
Vibration and Stability Analysis of Composite Spinning Shafts
Transactions of the Korean Society for Noise and Vibration Engineering. 2015. Jul, 25(7): 510-517
Copyright © 2015, The Korean Society for Noise and Vibration Engineering
  • Received : June 25, 2015
  • Accepted : July 09, 2015
  • Published : July 20, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
정석 서
Member, Department of Mechanical Engineering, Chungnam National University
창기 안
Member, Department of Mechanical Engineering, Chungnam National University
상윤 박
Member, Department of Mechanical Engineering, Chungnam National University
오섭 송
Corresponding Author ; Member, Department of Mechanical Engineering, Chungnam National University E-mail :songos@cnu.ac.kr

Abstract
The free vibration and stability analysis of a spinning composite shaft modelled as a thin-walled closed beam is performed for several design parameters, such as ply angle, aspect ratio, and spin speed. The governing equations of spinning shafts based on the Timoshenko beam theory are derived via Hamilton’s variational principle. Coriolis acceleration and anisotropy of constituent materials are incorporated in the derivation. The equations of motion are then transformed to the standard form of an eigenvalue problem for free vibration and stability analysis. Analytical results both for uniform circular cylindrical shaft and rectangular cross-section shaft are obtained by using extended Galerkin method, and the results are compared with those from FEM ANSYS analysis for a verification.
Keywords
1. 서 론
섬유강화 복합재료(fiber-reinforced composites)는 높은 비강도(specific strength) 및 비강성(specific stiffness), 내부식성, 높은 피로 수명, 낮은 열팽창률 등의 장점 때문에 토목, 건축, 기계, 항공 우주 산업 등에서 광범위하게 사용되고 있다. 또한 최근 구조물의 자중 절감 및 구조적 효율성을 높이기 위해 박판 보(thin-walled beam)의 사용이 항공, 우주, 기계, 토목, 건축분야에서 증가되고 있다. 특히, 복합재료 회전축은 항공 추진시스템 및 헬리콥터 구동장치, 산업기계, 가스터빈, 터보 발전기, 정밀기계 등의 동력전달장치로서 사용이 기대되고 있으며, 장비의 신뢰도는 회전축의 안정성에 크게 의존한다 (1) . 이러한 회전축의 진동특성을 향상시키고, 시스템의 정상작동에 위험이 되는 불안정성의 발생을 제거하거나 지연시켜야 하는 필요성이 제기된다.
Bert and Kim은 복합재료의 이방성에 의한 복합재료 회전축의 굽힘-비틀림 연성효과를 고려한 안정성 문제를 연구하였다 (2) . Lancaster는 감쇠를 무시한 회전체의 고유치와 고유벡터가 공액복소수 형태로 이루어진다는 것을 밝혀냈다 (3) . 복잡한 형상의 회전체의 유한요소 해석이 Linkins에 의해서 수행되었으며 (4) , 이를 이용하여 Patel과 Seltzer가 일정한 각속도로 회전하는 구조물의 탄성변형만이 고려된 고유값 문제를 풀 수 있는 컴퓨터 프로그램을 제시하였다 (5) .
이 연구에서 축을 중심으로 일정 각속도로 회전하는 원형 및 사각단면 복합재료 회전축의 Timoshenko보 이론에 기초한 회전체동역학 구조 모델을 정립하고 이를 ANSYS 유한요소 해석을 통하여 비교 검증하여 이 연구에서 제시한 구조 모델이 회전체 진동 설계에 사용될 수 있음을 보였다. 또한 섬유각과 사각단면의 종횡비 등의 변수들이 회전축의 위험속도에 미치는 영향을 도출하여 그 값들을 비교 검증하였다.
2. 박판 보의 운동학
- 2.1 변위장
원형단면 및 사각단면 박판 보(thin-walled beam)에서 임의의 한 점 A의 변위장을 6개의 변위변수(kenematic variables)로 표시하면 다음과 같다( Fig. 1 참조).
PPT Slide
Lager Image
Coordinates of composite thin-walled beam, (X,Y,Z): global coordinates, (n,s,z): local(surface) coordinates
PPT Slide
Lager Image
여기서 u 0 , v 0 , w 0 는 각각 x , y , z 축 방향의 병진변위를 의미하며, θ x , θ y , ϕ x , y , z 축 중심의 회전변위를 나타낸다. F w , na 는 각각 1, 2차 와핑 함수를 나타내며 단면이 원형인 경우에 와핑 함수는 영이 된다.
- 2.2 속도, 가속도 벡터
변형된 보의 임의의 한 점 M ( x , y , z ) 의 위치벡터, 속도벡터 그리고 가속도 벡터는 다음과 같다.
PPT Slide
Lager Image
여기에서
PPT Slide
Lager Image
는 각각 x , y , z 축 방향의 단위벡터이며,
PPT Slide
Lager Image
이다.
PPT Slide
Lager Image
Configurations of composite circular and rectangular cross-section shafts
3. 지배방정식 및 경계조건
- 3.1 해밀턴의 변분원리
아래 해밀턴 변분원리를 이용하여 원통 및 사각단면 보의 운동방정식과 경계조건을 유도한다.
PPT Slide
Lager Image
여기서 δU δK 는 각각 탄성에너지와 운동에너지의 변분이다.
구조물에 저장되는 탄성변형에너지는, 박판 보의 단면은 면내 변형을 하지 않는다는 가정 하에 ( xx , nn , γsn =0) 다음과 같이 유도된다.
PPT Slide
Lager Image
여기서 N zz , N sz , N ns 은 복합재료 적층 구조물의 stress resultants이고 L zz 는 stress couple이다.
박판 보의 운동에너지는 다음과 같다.
PPT Slide
Lager Image
- 3.2 운동방정식 및 경계조건
식 (4),(5)를 식 (3)에 대입하여 정리하면 다음 운동방정식과 경계조건을 동시에 유도할 수 있다 (6) .
종방향 및 횡방향 연성운동
PPT Slide
Lager Image
Clamped-free 경계조건
PPT Slide
Lager Image
PPT Slide
Lager Image
축방향 및 비틀림 연성운동
PPT Slide
Lager Image
Clamped-free 경계조건
PPT Slide
Lager Image
PPT Slide
Lager Image
식 (5)~(9)에서 a ij 는 global stiffness, b i 는 질량 또는 관성 모멘트를 나타낸다.
4. 고유치 문제
자유진동의 동기운동(synchronous motion)을 가정하면 변위 요소는 다음과 같이 표현된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 식 (6)의 gyroscopic항 때문에 λ는 복소수 형태의 고유치가 되며 trial function u j , v j , x j , y j 는 구조물의 경계조건을 최대한 만족시킬 수 있는 admissible 함수이다.
확장된 Galerkin 방법(EGM : extended Galerkin method)을 적용하기위해 식 (11), (12)를 해밀턴 원리식 (2)에 대입한 후에 z -좌표에 대하여 적분을 수행하면 다음과 같은 행렬식을 얻는다.
PPT Slide
Lager Image
B T = [ a 1 , a 2 a N , b 1 , b 2 , ⋯ b N , c 1 , c 2 , ⋯ c N , d 1 , d 2 , ⋯ d N ]
특성방정식은 다음과 같다.
PPT Slide
Lager Image
여기에서 [M]과 [K]는 질량행렬, 강성행렬로서 실수 대칭행렬(real symmetric matrix)이며, [G]는 Gyroscopic 행렬로서 실수 역대칭 행렬(real skew symmetric matrix)이다.
식 (12)는 다음과 같은 고유치문제로 변형시켜 그 해를 구할 수 있다.
PPT Slide
Lager Image
여기에서 Z행렬의 고유치를 구하면 항상 N쌍의 공액복소수(complex conjugate)로 나타나며, 고유벡터 또한 N쌍의 공액복소수로 이루어진다.
5. 유한요소 해석
이 연구의 수치해석에 사용된 복합재료의 물성치 및 크기는 Table 1 과 같다.
Material properties of composite material
PPT Slide
Lager Image
Material properties of composite material
Elements of ANSYS model
PPT Slide
Lager Image
Elements of ANSYS model
Figs. 3 , 4 는 원형단면 및 사각단면 회전축의 ANSYS 모델을 나타낸다. ANSYS 프로그램을 사용하여 섬유각 θ =0, 30, 45, 60, 90인 5가지 경우에 대하여 1차 고유진동수 변화를 고찰하였다. Shell281 요소의 자유도는 UX, UY, UZ, ROTX, ROTY, ROTZ 6개이며 고정단에서는 6개의 자유도를 모두 0으로 고정하였고 자유단에서는 아무 구속도 하지않았다.
PPT Slide
Lager Image
Circular composite shaft-ANSYS model
PPT Slide
Lager Image
Rectangular composite shaft-ANSYS model
6. 결과 및 고찰
Tables 3 , 4 는 6개의 섬유각에 대하여 식 (5)~(10)의 고유치 문제(EVP, beam theory)를 풀어 얻은 결과와 ANSYS 해석결과를 비교한 것이다.
Comparison of fundamental frequencies from EVP(beam theory) and FEM for different ply angles(lateral-vertical bending coupled motion)
PPT Slide
Lager Image
Comparison of fundamental frequencies from EVP(beam theory) and FEM for different ply angles(lateral-vertical bending coupled motion)
Comparison of fundamental frequencies from EVP and FEM for different ply angles(extension-twist coupled motion)
PPT Slide
Lager Image
Comparison of fundamental frequencies from EVP and FEM for different ply angles(extension-twist coupled motion)
Table 3 Ω = 0 일 때의 종방향 및 횡방향 연성 1차 고유진동수를 비교한 표이고 Table 4 Ω = 0 일 때 축방향 및 비틀림 연성운동의 1차 고유진동수를 비교한 표이다.
종방향 및 횡방향 연성운동의 고유진동수는 섬유각 θ =0°, 30°에서는 오차가 5 % 이하이지만 섬유각이 증가함에 따라 그 오차가 커짐을 알 수 있다. 이는 보 이론에서는 굽힘 운동 시 면내 단면은 변형하지 않는다는 조건인 xx , nn , γsn =0으로 가정하지만 ANSYS에서는 면내 단면의 변형이 허용되어 보이론의 경우보다 고유진동수가 낮아지는 것으로 판단된다. 축방향-비틀림 연성 고유진동수는 유한요소 해석과 보 이론의 차이가 모든 섬유각의 경우에 대하여 1.1 % 이하 인 것을 알 수 있다.
PPT Slide
Lager Image
1st mode shapes of circular and rectangular cross-section shafts(lateral-vertical bending coupled motion)
PPT Slide
Lager Image
Fundamental mode shapes of circular and rectangular cross-section shafts(extension-twist coupled motion)
Figs. 7 , 8 은 원형단면 회전축의 경우에 섬유각 변화에 따른 종-횡방향 연성 고유진동수의 변화를 나타낸다 (6) . 가로축은 무차원 회전속도( Ω / ω 0 ) 이고 세로축은 무차원 고유진동수( ω / ω 0 ) 이다. Figs. 7 , 8 에서 섬유각 마다 회전속도 변화에 따라 고유진동수가 영(0)이 되는 위험속도(critical speed)가 존재함을 알 수 있다. 위험속도에서는 고유치 λ = σ + 의 실수값이 양(+)이 되어 시간에 따라 운동이 지수적으로 증가하는 발산현상(divergence)이 관찰된다. 위험속도는 섬유각이 클수록 증가하는 것을 알 수 있다.
PPT Slide
Lager Image
Variation of lateral-vertical coupled bending natural frequencies of circular shaft with spin speed for different fiber angles(ANSYS)
PPT Slide
Lager Image
Counterpart of Fig. 7(EVP, beam theory)
또한 Figs. 7 , 8 에서 Ω = 0 일 때 하나였던 고유진동수가 회전속도 증가에 따라 forward와 backward 두 개의 진동수로 분기되는 현상을 확인할 수 있다. Fig. 9 는 원형 단면 회전축의 회전속도 변화에 따른 extension-twist 연성 진동수의 변화이다. Figs. 7 , 8 과 달리 회전속도 증가에 따라 고유진동수는 지속적으로 감소한다.
PPT Slide
Lager Image
Variation of extension-twist coupled natural frequencies of circular shafts with spin speed for different fiber angles
Figs. 10 , 11 은 섬유각 θ =0˚, 30˚, 45˚, 60˚, 90˚인 사각단면 회전축에 대하여 회전속도 증가에 따른 고유진동수의 변화를 보여주는 그래프이고 원형단면 회전축과 마찬가지로 섬유각이 증가할수록 고유진동수가 증가하고 위험속도도 증가함을 알 수 있다 (7) .
PPT Slide
Lager Image
Lateral-vertical bending natural frequencies vs. spin speed for a rectangular shaft(ANSYS)
PPT Slide
Lager Image
Counterpart of Fig. 10(EVP, beam theory)
Figs. 12 , 13 은 회전축 단면의 종횡비에 따른 고유진동수 변화를 보여주는 그래프이다. 종횡비 R은 가로의 길이가 고정일 때 가로에 대한 세로의 비를 의미한다. R=0.5일 때 위험속도 영역이 가장 낮았으며 R이 커질수록 위험속도 영역이 증가함을 확인할 수 있다 (8) . 이는 단면 비대칭성에 의하여 불안정영역이 한 점이 아닌 구간으로 나타난다 (9) .
PPT Slide
Lager Image
Variation of natural frequencies of rectangular cross-section shaft with spin speed for different aspect ratios(ANSYS)
PPT Slide
Lager Image
Counterpart of Fig. 12(EVP, beam theory)
Fig. 14 는 섬유각 θ =0°이고 R=1, 0.5, 1.5, 2인 사각단면 회전축에 대하여 ANSYS와 EVP(beam theory)에 의한 결과를 비교한 것으로서 두 결과가 매우 잘 일치함을 보여 주고 있다. R이 1이 아닌 경우에는 위험 속도가 한 점이 아닌 영역으로 나타남을 알 수 있다. 또한 고유진동수의 분기(forward and backward)를 확인할 수 있다.
PPT Slide
Lager Image
Comparison of the results from ANSYS and EVP(beam theory) for a rectangular cross-section shaft
7. 결 론
이 연구를 통하여 다음과 같은 결론을 얻었다.
(1) 종방향 및 횡방향 굽힘 연성 운동에서 회전축 단면의 종횡비가 커질수록 고유진동수가 증가하고 불안정 영역이 증가함을 알 수 있다.
(2) 섬유각에 따른 영향은 종방향 및 횡방향 굽힘 연성 운동에서는 섬유각이 0°에서 90°로 증가할수록 고유진동수가 증가하고 불안정영역 또한 더 큰 회전 속도 방향으로 이동하였으나 축방향 및 비틀림 연성 운동에서는 별다른 변화를 보이지 않았다.
(3) 이 연구에서 제시된 보 이론에 기초한 고유치 문제(EVP)는 ANSYS 유한요소 모델 결과와 비교하여 그 변화 경향이 유사함을 확인할 수 있었다. 따라서 이 연구에서 제시된 회전축 모델을 통하여 다양한 parametric study가 가능하여 복합재료 회전축 기초설계의 도구로 사용될 수 있음을 보였다.
기 호 설 명
c.b : 사각단면 박판 보의 가로, 세로 길이 h : 박판 보의 두께 L,R,R : 원형단면 박판 보의 길이와 반경 (n,s,z) : 반경방향, 원주방향, 축방향 좌표 Ω : 회전속도 rad/s θ : 섬유각 θx,θy : X, Y축에 관한 회전 변위 u0,v0 : X, Y축 방향의 병진 변위
Acknowledgements
이 연구는 교육부와 한국연구재단의 지역혁신인력양성사업으로 수행된 연구결과임(No.2012026181).
BIO
Jungseok Seo received his B.S. degree from the Department of Mecahnical Engineering at Chungnam National University in 2012. He is currently the master’s course in Chungnam National University His research interests encompass vibration and control of composite structures and rotor dynamics.
Ohseop Song received his B.S. degree from the Department of Mechanical Design Engineering at Seoul National University in 1978 and M.S. degree from Mechanical Engineering at New Jersey Institute of Technology. His Ph.D. degree was awarded by Virginia Polytechnic Institute and State University. He is currently a professor in the Department of Mechanical Engineering at Chungnam National University, Daejeon, Korea. His research interests encompass vibration and control of composite structures, instability of nanopipes, and rotor dynamics.
References
Bauchau O. 1981 Design, Manufacturing and Testing of High Speed Rotating Graphite-epoxy Shaft, Ph.D. Thesis Massachusetts Institute of Technology, Department of Aeronautical and Astronautics
Bert C. W. , Kim C. D. 1995 Whirling of Composite–material Drive-shaft Including Bending-twisting Coupling and Transverse Shear Deformation Journal of Vibration and Acoustics 117 (1) 17 - 21
Lancaster P. 1972 Lambda-matrices and Vibration Systems Pergamon Press Oxford, England 7.7 Sec
Linkins P. W. 1972 Finite Element Appendage Analysis International Journal of Solids and Structures 8 (5) 709 - 731    DOI : 10.1016/0020-7683(72)90038-8
Patel J. S. , Seltzer S. M. 1972 Complex Eigenvalue Analysis of Rotating Structures in NASTRAN ; User’s Experiences 197 - 234
Kim S.-K. 2010 Dynamic Response Analysis of Composite H-type Cross-section Beams Transactions of the Korean Society for Noise and Vibration Engineering 20 (6) 583 - 592    DOI : 10.5050/KSNVE.2010.20.6.583
Librescu L. , Song O. S. 2006 Thin-walled Composite Beams : Theory and Application Springer
Yoon H. W. 2004 Vibration and Stability of Composite Thin-walled Spinning Shaft Proceedings of the KSNVE Annual Autumn Conference 1083 - 1088
Lee C.-W. 1993 Vibration Analysis of Rotor Kluwer Academic Publishers