![[기계항공실험] 동역학제어실험-1](http://images.happyhaksul.com/thumb/370/369370-0001.gif)
![[기계항공실험] 동역학제어실험-2](http://images.happyhaksul.com/thumb/370/369370-0002.gif)
![[기계항공실험] 동역학제어실험-3](http://images.happyhaksul.com/thumb/370/369370-0003.gif)
![[기계항공실험] 동역학제어실험-4](http://images.happyhaksul.com/thumb/370/369370-0004.gif)
![[기계항공실험] 동역학제어실험-5](http://images.happyhaksul.com/thumb/370/369370-0005.gif)
![[기계항공실험] 동역학제어실험-6](http://images.happyhaksul.com/thumb/370/369370-0006.gif)
![[기계항공실험] 동역학제어실험-7](http://images.happyhaksul.com/thumb/370/369370-0007.gif)
![[기계항공실험] 동역학제어실험-8](http://images.happyhaksul.com/thumb/370/369370-0008.gif)
![[기계항공실험] 동역학제어실험-9](http://images.happyhaksul.com/thumb/370/369370-0009.gif)
![[기계항공실험] 동역학제어실험-10](http://images.happyhaksul.com/thumb/370/369370-0010.gif)
![[기계항공실험] 동역학제어실험-11](http://images.happyhaksul.com/thumb/370/369370-0011.gif)
![[기계항공실험] 동역학제어실험-12](http://images.happyhaksul.com/thumb/370/369370-0012.gif)
![[기계항공실험] 동역학제어실험-13](http://images.happyhaksul.com/thumb/370/369370-0013.gif)
![[기계항공실험] 동역학제어실험-14](http://images.happyhaksul.com/thumb/370/369370-0014.gif)
![[기계항공실험] 동역학제어실험-15](http://images.happyhaksul.com/thumb/370/369370-0015.gif)
![[기계항공실험] 동역학제어실험-16](http://images.happyhaksul.com/thumb/370/369370-0016.gif)
![[기계항공실험] 동역학제어실험-17](http://images.happyhaksul.com/thumb/370/369370-0017.gif)
![[기계항공실험] 동역학제어실험-18](http://images.happyhaksul.com/thumb/370/369370-0018.gif)
![[기계항공실험] 동역학제어실험-19](http://images.happyhaksul.com/thumb/370/369370-0019.gif)
![[기계항공실험] 동역학제어실험-20](http://images.happyhaksul.com/thumb/370/369370-0020.gif)
분야
hwp1-(a) 주어진 데이터를 이용하여 x 축을 주파수, y 축을 축의 길이, z 축을 응답의 크기(magnitude)로 하여 3 차원의 그래프를 1 번, 3 번, 5 번, 그리고 7 번 노드에서 그리시오. 이 경우 복소수로 주어진 데이터는 응답의 절대값과 위상이다. 이를 실수와 허수 값으로 각각 나타내어 그래프를 그려보고, 고유 진동수를 찾으시오.
① Magnitude를 이용한 3차원 그래프
② 실수와 허수부를 이용한 3차원 그래프
③ 그래프 비교 및 고유 진동수
1-(b) 위의 그림에 도시된 축 시스템의 경우 이론적으로 데이터가 Maxwell의 상반정리를 만족하게 된다. 즉, 1의 노드에 가속도계를 장착하고 3번에서 임팩트 해머를 때리는 것이나 3번에 가속도계를 장착하고 1번에서 임팩트 해머를 때리는 것의 주파수 응답은 동일하다. 그러나, 실험과정에서 이런 상반 정리가 잘 적용되지 않는다. 이 이유를 설명하시오.
1-(c) 시스템을 수학적 해석을 위해 모델링을 했다면 그 모델링이 제대로 된 것인지 검증하는 것이 중요하다. 일반적으로 시스템의 해석방법에 대해서 불변값인 고유 진동수를 통한 검증을 시행하게 된다. 그 중에서 가장 일반적은 방법이 바로 임펄스 테스트이다. 임펄스 테스트의 장․단점에 대해 조사하고 이 이외에 다른 어떤 검증 방법이 있는지 설명하시오.
1-(d) 임펄스 테스트를 통해 얻은 신호는 가속도계를 이용하여 측정할 수 있다. 이렇게 측정된 데이터를 적분하여 속도 및 변위 데이터를 얻을 수 있다. 그러나 이렇게 구한 데이터는 시간 도메인의 데이터이므로 고유값을 구하기 위해서는 주파수 도메인으로 변형시킬 필요가 있다. 이는 퓨리에 급수 혹은 변환을 통해 구할 수 있다. 연속적인 함수의 형태로 얻는 응답의 경우 일반적인 식에 대입하여 구할 수 있지만 이산적으로 얻어진 데이터의 경우 이산 푸리에 변환이 필요하다. 다음의 이산 푸리에 식을 이용하여 응답 스펙트럼의 크기 및 위상을 구해보라.
MATLAB CODE of Pro.1
2. 제프콧 로터 시스템
2-(a) 가장 간단한 가정은 양 옆에 지지되는 베어링 성분에서 강성은 무한대이고 댐핑은 존재하지 않는 경우이다. 질량은 가운데의 디스크에 집중되어 있으며 디스크는 회전의 자유도가 구속되어 있다. 사용되는 변수는 다음과 같다.
2-(b) (a)의 경우에는 디스크의 병진(translation) 운동만이 가능했다. 베어링의 강성과 댐핑이 유한한 경우에는 베어링 부분에서의 병진 운동도 고려해야 한다. 디스크의 회전의 자유도가 구속되어 있다는 가정에서 베어링의 강성과 댐핑이 각각 라고 할 때 (a)의 변수들을 참고해서 운동 방정식을 구하시오.
2-c) 디스크의 회전의 자유도가 있는 보다 실제에 가까운 모델의 경우 와 함께 을 참고한다.
2-(d) 축 강성과 베어링 강성 비의 변화에 따른 고유 진동수의 변화를 논하고 그래프로 그려보아라.
제프콧 로터의 가능한 궤적
MATLAB CODE of Pro.2-d
2-(e) e=15, =3, =6, m=10, t=1 일 때, 회전속도(Ω)에 따른 기하학 중심과 무게 중심의 응답을 그래프로 그려보아라. 이 때 고유진동수가 두 군데로 갈라지는 이유와 두 고유 진동수 사이에서 백워드 훨이 나옴을 확인하라.
① 고유진동수가 두 군데로 갈라지는 이유
②백워드 휠 확인
MATLAB CODE of Pro.2-e
2-(f) 축의 회전이 첫 번째 고유 진동수보다 작을 범위일 때 축의 가진 방향(불평형 질량의 위치)과 같은 방향으로 축의 응답 반응이 나타난다. 즉 불평형 질량이 위쪽에 놓여있으며 가진이 위쪽으로 작용할 때 축은 위쪽으로 휘어있다. 그리고 두 번째 고유 진동수보다 큰 범위에서는 축의 가진 방향과 반대 방향으로 축의 응답이 나타나게 된다. 그리고 두 고유 진동수 사이에서 축이 회전할 때는 백워드 훨이 나타나게 되는데 이 같은 현상의 이유를 수식을 써서 설명하라.(단, 댐핑은 없다고 가정한다.)
①백워드 휠 이유
②축의 응답 방향 이유
3. 자기베어링 로터 시스템
3-(a) m=5kg, , , a=0.5m, b=-0.2m, k=200N/mm 라고 할 때 M, D, K 매트릭스를 만들고, state space에서의 방정식 형태로 바꾸어 A, B, C 매트릭스를 구하시오.
3-(b) 고유 진동수와 고유 벡터를 구하고 고유 벡터로 이 시스템의 모드를 알아볼 수 있는지 논의해보고, 제시된 모델의 한계를 기술하시오.
1) 고유진동수
2) 고유벡터
3)고유벡터가 직선으로 나오는 이유와 제시된 모델의 그 한계
3-(c) 강성의 변화에 따른 고유 진동수 그래프를 그리고 일반적인 회전체 시스템과 달리 진동수가 강성 증가에 따라서 계속 증가하는 이유를 서술하시오.
3-(d) 회전속도 Ω를 바꿔 가면서 A 매트릭스의 고유진동수를 구하고, 이를 위험속도선도(campbell Diagram)로 나타내시오. (x축 : 회전속도, y축 : 고유 진동수) 또, 이 고유 진동수 변화의 물리적 의미를 자세히 설명하시오. 또, 여기서 위험속도를 찾는 방법을 기술하시오.
3-(e) 스펙트럼 해석은 회전체 시스템을 진단하기 위한 가장 보편적은 신호처리 방법이다. 주파수 성분과 외에 상응하는 크기와 위상이 기계의 오류에 따라 변하기 때문이다. 기존의 스펙트럼 해석은 주파수 성분의 표현을 체계적으로 할 수 있으나, 로터의 진동을 실수의 양으로 표한하기 때문에 주파수 스펙트럼은 방향성-전방향 회전 또는 후방향 회전-과 같은 중요한 궤적 정보를 잃는다. 반면 복소 신호의 방향성 스펙트럼 해석은 방향성에 대한 정보, 평면운동에 대한 형상 정보를 가지고 있어 회전체를 진단하는 강력한 해석법이다. 이 같은 방향성 스펙트럼에 대해 조사하시오.
4. Magnetic Bearing
4-(a) 자기베어링은 다음과 같은 스펙을 가졌다. 이때 구동력 F는 어떻게 선형화 되는가?
4-(b) 3번 문제의 모델 중 f에 (a)에서 구한 선형화된 식을 넣고 동역학적 모델을 만들어 고유 진동수를 구하여 시스템이 불안정함을 보이시오.
(단, 자이로 효과는 무시한다. Hint; 자이로 효과를 무시하면 x축의 성분들과 y축의 성분들이 각각 독립적으로 계산 되어진다. 그러므로 x나 y 두 성분 중 한가지 성분만 해석해도 무리가 없다.)
5. 동기 및 비동기 가진
5-(a) 동기 가진과 비동기 가진이 일어나는 이유에 대해 간단하게 설명하고 이와 같은 예를 서술하시오.
5-(b) 학부실험 내용 중 자기베어링을 이용해 축에 10Hz(600rpm과 동등)로 조화 가진을 넣고 축을 300, 1200, 1800, 2400rpm으로 회전시키면서 축의 거동을 살펴보았다. 이 때 축의 궤적을 MATLAB을 통해 구현해보시오. 여기서 중요한 요소는 주파수이며 크기는 큰 영향이 없다.
MATLAB CODE of Pro.5-b
5-(c) 비동기 가진의 주파수가 600rpm이고 축의 회전 속도가 600rpm일 때 맥놀이(beat) 현상이 일어나게 된다. 이를 시간 축에서 구현해 보시오.
MATLAB CODE of Pro.5-c
FEM 이나 전달 행렬 방법으로 모델링 한 회전축을 실험적으로 검증하기 위해서 임펄스 테스트를 수행한다. 베어링이 없는 축을 매달고 매단 방향과 직각으로 가속도계를 축에 붙이고 축을 임팩트 해머로 때려 그 둘 사이의 주파수 응답을 구한다. 이러한 주파수 응답을 통해서 회전축의 고유진동수와 모드 형상을 측정할 수 있고 이론적으로 구성된 모델과 비교해서 모델의 신뢰성을 검증한다. 이 과제와 함께 첨부되는 데이터는 이번 실험을 통하여 그림 1 의 축을 임펄스 테스트를 통해서 얻은 것이다. 아래 그림과 같이 8 개의 지점을 정한 후 i 번째 위치에 가속도계를 두고 j 번째 위치에서 임팩트 해머로 때리고 FFT analyzer 를 이용하여 그 전달함수를 얻는다. 전달함수 데이터는 magnitude, phase 가 각각 Mij , Pij (degree)로 저장되었으며 총 64 번 (8x8)의 실험을 통해 각각 64 개의 데이터가 주어진다. 주파수데이터는 frequency(Hz)로 주어진다. 가속도계 위치는 sensor로 주어진다.
1-(a) 주어진 데이터를 이용하여 x 축을 주파수, y 축을 축의 길이, z 축을 응답의 크기(magnitude)로 하여 3 차원의 그래프를 1 번, 3 번, 5 번, 그리고 7 번 노드에서 그리시오. 이 경우 복소수로 주어진 데이터는 응답의 절대값과 위상이다. 이를 실수와 허수 값으로 각각 나타내어 그래프를 그려보고, 고유 진동수를 찾으시오.
① Magnitude를 이용한 3차원 그래프
Data로부터 magnitude 에 대해 3차원 그래프를 그려보았다. 그리고 고유진동수를 찾아내기 위해 3차원 그래프를 frequency와 응답의 2차원 그래프로 변환하였고 고유진동수를 가장 잘 나타내어 주는 그래프 하나만 따로 그려보았다.
고유진동수(Hz): 167, 451, 718
고유진동수(Hz) : 171, 453, 714
고유진동수(Hz) : 172, 439, 720
고유진동수(Hz) : 171, 438, 715
