분야
hwp(a) 주어진 데이터는 응답의 절대값과 위상이다. x축을 주파수, y축을 축의 길이로 하고, z축은 각각 응답의 크기, 위상, 실수부, 허수부 값을 가지는 그래프를 1, 3, 5, 7번 노드에 대해 그려보고 고유 진동수를 찾으시오. 이 때 어느 그래프를 사용하는 것이 고유 진동수를 찾는데 유리하며, 그 이유는 무엇인가
(b) 위의 그림에 도시된 축 시스템의 경우 이론적으로 데이터가 Maxwell의 상반정리를 만족하게 된다. 즉, 1의 노드에 가속도계를 장착하고 3번에서 임팩트 해머를 때리는 것이나 3번에 가속도계를 장착하고 1번에서 임팩트 해머를 때리는 것의 주파수 응답은 동일하다. 그러나 실험과정에서 이런 상반 정리가 잘 적용되지 않는다. 이 이유를 설명하시오.
(c) 시스템을 수학적 해석을 위해 모델링을 했다면 그 모델링이 제대로 된 것인지 검증하는 것이 중요하다. 일반적으로 시스템의 해석방법에 대해서 불변 값인 고유 진동수를 통한 검증을 시행하게 된다. 그 중에서 가장 일반적은 방법이 바로 임펄스 테스트이다. 임펄스 테스트의 장․단점에 대해 조사하고 이 이외에 다른 어떤 검증 방법이 있는지 설명하시오
임펄스 테스트의 장,단점
이 외의 다른 검증 방법
(d) 임펄스 테스트를 통해 얻은 신호는 가속도계를 이용하여 측정할 수 있다. 이렇게 측정된 데이터를 적분하여 속도 및 변위 데이터를 얻을 수 있다. 그러나 이렇게 구한 데이터는 시간 도메인의 데이터이므로 고유값을 구하기 위해서는 주파수 도메인으로 변형시킬 필요가 있다. 이는 퓨리에 급수 혹은 변환을 통해 구할 수 있다. 연속적인 함수의 형태로 얻는 응답의 경우 일반적인 식에 대입하여 구할 수 있지만 이산적으로 얻어진 데이터의 경우 이산 푸리에 변환이 필요하다. 다음의 이산 푸리에 식을 이용하여 응답 스펙트럼의 크기 및 위상을 구해보라.
2. 제프콧 로터 시스템
(a) 가장 간단한 가정은 양 옆에 지지되는 베어링 성분에서 강성은 무한대이고 댐핑은 존재하지 않는 경우이다. 질량은 가운데의 디스크에 집중되어 있으며 디스크는 회전의 자유도가 구속되어 있다. 사용되는 변수는 다음과 같다.
(b) (a)의 경우에는 디스크의 병진(translation) 운동만이 가능했다. 베어링의 강성과 댐핑이 유한한 경우에는 베어링 부분에서의 병진 운동도 고려해야 한다. 디스크의 회전의 자유도가 구속되어 있다는 가정에서 베어링의 강성과 댐핑이 각각라고 할 때 (a)의 변수들을 참고해서 운동 방정식을 구하시오.
(c) 디스크의 회전의 자유도가 있는 보다 실제에 가까운 모델의 경우 와 함께 을 참고한다.
(d) 축 강성과 베어링 강성 비의 변화에 따른 고유 진동수의 변화를 논하고 그래프로 그려보시오.
(e) e=15, kx=3, ky=6, m=10, t=1 일 때, 회전속도(Ω)에 따른 기하학 중심과 무게 중심의 응답을 그래프로 그려 보아라. 이 때 고유진동수가 두 군데로 갈라지는 이유와 두 고유 진동수 사이에서 백워드 훨이 나옴을 확인하라.
(f) 축의 회전이 첫 번째 고유 진동수보다 작을 범위일 때 축의 가진 방향(불평형 질량의 위치)과 같은 방향으로 축의 응답 반응이 나타난다. 즉 불평형 질량이 위쪽에 놓여있을 때 가진이 위쪽으로 작용할 때 축은 위족으로 휘어있다. 그리고 두 번째 고유 진동수보다 큰 범위에서는 축의 가진 방향과 반대 방향으로 축의 응답이 나타나게 된다. 그리고 두 고유 진동수 사이에서 축이 회전할 때는 백워드 훨이 나타나게 되는데 이 같은 현상의 이유를 수식을 써서 설명하라. (단, 댐핑은 없다고 가정한다.)
3. 자기베어링 로터 시스템
(a) m=5kg, , , a=0.5m, b=-0.2m, k=200N/mm 라고 할 때 M, D, K 매트릭스를 만들고, state space에서의 방정식 형태로 바꾸어 A, B, C 매트릭스를 구하시오.
(b) 고유 진동수와 고유 벡터를 구하고 고유 벡터로 이 시스템의 모드를 알아볼 수 있는지 논의해보고, 제시된 모델의 한계를 기술하시오.
고유진동수와 고유 벡터
고유 벡터와 시스템의 모드
제시된 모델의 한계
(c) 강성의 변화에 따른 고유 진동수 그래프를 그리고 일반적인 회전체 시스템과 달리 진동수가 강성 증가에 따라서 계속 증가하는 이유를 서술하시오.
(d) 회전속도 Ω를 바꿔 가면서 A 매트릭스의 고유진동수를 구하고, 이를 위험속도선도(campbell Diagram)로 나타내시오. (x축 : 회전속도, y축 : 고유 진동수) 또, 이 고유 진동수 변화의 물리적 의미를 자세히 설명하시오. 또, 여기서 위험속도를 찾는 방법을 기술하시오.
(e) 스펙트럼 해석은 회전체 시스템을 진단하기 위한 가장 보편적은 신호처리 방법이다. 주파수 성분과 외에 상응하는 크기와 위상이 기계의 오류에 따라 변하기 때문이다. 기존의 스펙트럼 해석은 주파수 성분의 표현을 체계적으로 할 수 있으나, 로터의 진동을 실수의 양으로 표한하기 때문에 주파수 스펙트럼은 방향성-전방향 회전 또는 후방향 회전-과 같은 중요한 궤적 정보를 잃는다. 반면 복소 신호의 방향성 스펙트럼 해석은 방향성에 대한 정보, 평면운동에 대한 형상 정보를 가지고 있어 회전체를 진단하는 강력한 해석법이다. 이 같은 방향성 스펙트럼에 대해 조사하시오
4. Magnetic Bearing
(a) 자기베어링은 다음과 같은 스펙을 가졌다. 이때 구동력 F는 어떻게 선형화 되는가?
(b) 3번 문제의 모델 중 f에 (a)에서 구한 선형화된 식을 넣고 동역학적 모델을 만들어 고유 진동수를 구하여 시스템이 불안정함을 보이시오.
5. 동기 및 비동기 가진
(a) 동기 가진과 비동기 가진이 일어나는 이유에 대해 간단하게 설명하고 이와 같은 예를 서술하시오.
(b) 학부실험 내용 중 자기베어링을 이용해 축에 10Hz(600rpm과 동등)로 조화 가진을 넣고 축을 300, 1200, 1800, 2400rpm으로 회전시키면서 축의 거동을 살펴보았다. 이 때 축의 궤적을 MATLAB을 통해 구현해보시오. 여기서 중요한 요소는 주파수이며 크기는 큰 영향이 없다.
(c) 비동기 가진의 주파수가 600rpm이고 축의 회전 속도가 600rpm일 때 맥놀이(beat) 현상이 일어나게 된다. 이를 시간 축에서 구현해 보시오.
*Reference*
따라서 위의 식으로부터 주파수 응답은 가 된다. 이렇게 구해진 주파수 응답을 Euler Equation에 의해 다음과 같이 실수부분과 허수 부분으로 나눌 수 있게 된다.
여기서 주파수 응답의 실수부는 가 되고, 주파수 응답의 허수부는 가 된다. 이렇게 구해진 실수부와 허수부의 응답을 앞에서 그래프로 표현하였다. 그 결과 실수부보다는 허수부에서 고유 진동수 부분이 더 확연하게 드러나는 것을 알 수 있었다. 그 이유를 알아보도록 하자.
허수부가 실수부보다 고유주파수를 파악하기 쉬운 이유
다음 식에서 Phase의 Delay인 에 대해서 생각해보자.
위의 식을 정리하면 아래와 같이 에 대해서 정리할 수 있다.
여기에 다음의 3가지 식을 대입하면,
, ,
아래와 같이 에 대해서 정리할 수 있다.
여기에서 감쇠비 는 1보다 매우 작은 값이므로 이라고 하고 그래프를 그리면 다음과 같은 그림이 나온다.
위의 그래프를 보면 Ratio of Frequency(주파수비, 진동수비)가 1에 가까운 경우에 Phase(위상)이 90이라고 할 수 있고, 이외의 경우에는 0(또는 180)으로 간다고 할 수 있다. 즉, 가진주파수가 고유주파수와 일치하거나 비슷해지는 경우에 위상이 90도에 가까워지고 가진주파수가 고유주파수와 어느 정도 차이가 날수록 위상이 0도나 180도 쪽으로 간다. 이 때 가진주파수가 고유주파수와 거의 일치하는 경우는 값이 1에 가까워지게 된다. 여기서 주파수 응답의 실수부는 이고, 주파수 응답의 허수부는 이므로 고유 주파수 부근에서 허수 그래프는 최대값을 가지게 된다. 이에 비해서 실수부는 고유 주파수 부근까지는 응답이 계속 증가하다가 고유 주파수에서는 응답이 갑자기 0이 되어야 한다. 위의 그래프에서도 같은 사실을 확인할 수 있다. 따라서 허수부에서 응답의 최대값을 찾아서 고유주파수를 구하는 방식이 더 간단하다고 할 수 있다.
데이터 처리 코드(Java)
import java.io.BufferedReader;
import java.io.BufferedWriter;
import java.io.FileReader;
import java.io.FileWriter;
import java.io.PrintWriter;
class Dataprocess
{
public static void main(String[] args) throws Exception
{
BufferedReader in11 = new BufferedReader(new FileReader("m71"));
BufferedReader in12 = new BufferedReader(new FileReader("m72"));
BufferedReader in13 = new BufferedReader(new FileReader("m73"));
BufferedReader in14 = new BufferedReader(new FileReader("m74"));
BufferedReader in15 = new BufferedReader(new FileReader("m75"));
BufferedReader in16 = new BufferedReader(new FileReader("m76"));
BufferedReader in17 = new BufferedReader(new FileReader("m77"));
BufferedReader in21 = new BufferedReader(new FileReader("p71"));
BufferedReader in22 = new BufferedReader(new FileReader("p72"));
BufferedReader in23 = new BufferedReader(new FileReader("p73"));
BufferedReader in24 = new BufferedReader(new FileReader("p74"));
BufferedReader in25 = new BufferedReader(new FileReader("p75"));
BufferedReader in26 = new BufferedReader(new FileReader("p76"));
BufferedReader in27 = new BufferedReader(new FileReader("p77"));
PrintWriter out = new PrintWriter(new BufferedWriter(new FileWriter("imaginary7.txt")));
String node1 = null;
String node2 = null;
String node3 = null;
String node4 = null;
String node5 = null;
String node6 = null;
String node7 = null;
String node11 = null;
String node21 = null;
String node31 = null;
String node41 = null;
String node51 = null;
String node61 = null;
String node71 = null;
double no1 = 0;
double no2 = 0;
double no3 = 0;
double no4 = 0;
double no5 = 0;
double no6 = 0;
double no7 = 0;
double no11 = 0;
double no21 = 0;
double no31 = 0;
double no41 = 0;
double no51 = 0;
double no61 = 0;
double no71 = 0;
while((node1=in11.readLine())!= null)
{
node2 = in12.readLine();
node3 = in13.readLine();
node4 = in14.readLine();
node5 = in15.readLine();
node6 = in16.readLine();
node7 = in17.readLine();
node11 = in21.readLine();
node21 = in22.readLine();
node31 = in23.readLine();
node41 = in24.readLine();
node51 = in25.readLine();
node61 = in26.readLine();
node71 = in27.readLine();
no1 = Double.parseDouble(node1);
no2 = Double.parseDouble(node2);
no3 = Double.parseDouble(node3);
no4 = Double.parseDouble(node4);
no5 = Double.parseDouble(node5);
no6 = Double.parseDouble(node6);
no7 = Double.parseDouble(node7);
no11 = Double.parseDouble(node11);
no21 = Double.parseDouble(node21);
no31 = Double.parseDouble(node31);
no41 = Double.parseDouble(node41);
no51 = Double.parseDouble(node51);
no61 = Double.parseDouble(node61);
no71 = Double.parseDouble(node71);
no1 = no1*Math.sin(no11*(Math.PI/180));
no2 = no2*Math.sin(no21*(Math.PI/180));
no3 = no3*Math.sin(no31*(Math.PI/180));
no4 = no4*Math.sin(no41*(Math.PI/180));
no5 = no5*Math.sin(no51*(Math.PI/180));
no6 = no6*Math.sin(no61*(Math.PI/180
2. http://hydyn.hanyang.ac.kr/cyberclass/vibration/lecturenote/chap7a/chap7a.pdf
3. 조치영• 이종원, 방향 스펙트럼을 이용한 회전체의 비대칭성 및 비등방성 진단, 한국소음진동공학회지v.3, no.3, 한국소음진동공학회, 1993
