십자성군의 비밀의 방

식(2)

눈이 아픈 식이다...
cs[y1(t)-y2(t)], k[y1(t)-y2(t)]
용수철의 복원력의 방향은 위쪽이다. 아래방향의 힘을 +로 하고있으니 분명 저 힘은 -일 것이다.
그리고 이 식은 y2(t)의 관계식에서도 볼 수 있다. 즉, y2(t)를 y1(t)에 대한 관계식으로 바꾸어 식(2)를 y1(t)과 f만의 관계식으로 바꿀 수 있다는 것이다.

나는 바보같이 m(s^2)y2(t)를 통째로 옮겨버려서 헤맸었다.
또는 용수철 k와 c를 무시하고 mg의 힘이 질량체 M에 영향을 미치는 식으로 하다가 헤맸었다.

-균형상태에서의 변화를 생각하자. y1(t)과 y2(t)는 균형상태에서의 변위로 생각하고(좌표와는 다르다)
-mg와 Mg는 평형상태에서 상쇄된 상태로 생각.

y2(t)를 y1(t)의 관계식으로

해답

y2(t)를 y1(t)에 대한 식에 대입하자.

Y(s)는 y1(t)였었다. f와 함께 전달함수로 하면서 s에 관한 식으로 변경.

결과적으로 위의 해답식이 도출된다.

(1)네이버 지식인 출처

 

(2)Maxwell material(출처, 위키피디아)

(3)Kelvin(캘빈)모델

(4)Maxwell material(2번째)

설명들어간다. 일단 (3)그림의 캘빈모델을 예로 들어보겠다.
그림은 귀찮기에 말로 설명한다.

작용하는 힘
중력 : 아래방향, mh
용수철 : 위쪽방향, mg+k[y(t)-u(t)]
댐퍼가 미치는 힘 : 위쪽방향, b(d/dt){y(t)-u(t)}

몇가지 설명이 잘 안되있는것이 있다. 정의를 하겠다.
일단 입력을 u(t)라고 하고, 출력을 y(t)라고 한다.
b는 댐퍼의 점성마찰계수, k는 용수철상수이다.
u(t)는 천정의 좌표이다. 저 그림은 다른곳에서 가져왔기에 표시되어 있지 않지만. 입력변수를 천정으로 하였다. 즉, 천정이 위아래로 움직일 수 있다는 뜻이겠지.
y(t)는 질량m의 물체의 좌표이다.(단, 천정기준이 아니다. 공간기준으로 생각하자)

일단 기본조건을 잡자.
u(t)=0일때, 그리고 질량m에 의한 중력 mg(아래방향)와, 댐퍼와 용수철의 힘(위쪽방향)이 균형을 갖출때를 기본으로 삼는다.

여기서 주의할것은 용수철에 대한 위의 식이다.
댐퍼쪽은 별 문제가 없다고 생각하지만...
용수철의 경우는 다시 한번 잘 생각하자.
일단 나는 기준을, 질량 m에 의한 아래방향 힘과 용수철 및 댐퍼에 의한 힘이 균형을 갖춘상태로 잡았다. 즉, 그 상태가 0인 상태라는 것이다. 그 상태에서 y(t)=0 이며 u(t)=0 이다.

착각하면 안될것은 내가 좌표라고 했다고 단순하게 천정에서 잰 용수철의 길이 자체를 y(t)로 생각하고 k에 곱해서는 안된다는 것이다.
공학도라면 알고 있듯이, kx(t)에서 x(t)란 kx(t)=0인 상태의 용수철 길이를 기준으로 삼고있다. x(t)란 그 상태를 0으로 잡은 후의 '변위'이다. 이를 착각하지 말자.
용수철 식에서 mg가 들어간 이유는 '질량과 균형'인 상태이기 때문이다. 균형이란 상태는 변화가 없는 상태이다. 당연히 y(t)와 u(t)변화가 없으니 댐퍼의 힘은 없는 상태로 하고, 그 상태에서 질량m의 물체는 계속해서 mg로 아래방향 힘을 가하며, 이에 대해 변위가 없다는 것은 용수철이 물체에 대해서 위쪽 방향으로 mg의 힘을 가하고 있는 상태라는 것이다. 따라서 기본 mg를 갖춘 상태에서 용수철의 변위에 따라 용수철쪽의 힘이 변하는 것이다.


알다시피 '계'에서 모든 힘은 평형인 상태이다.
물체가(이제부터 M이라 하겠다.), 좌표가 변할경우, 즉 M이 아래방향으로 움직이는 식으로 y(t)가 변할경우를 생각하자.

M이 아래로 움직인다는 것은 M에 대해 mg말고도 '여분'의 힘이 작용한다고 생각하자.
아래로 내려가는 물체 M의 힘은
m[(d/dt)^2]y(t)라고 할 수 있다.
물체 mg에 주어지는 힘은
mg-{mg+k[y(t)-u(t)]}-b(d/dt){y(t)-u(t)}이다.

이에 따른 결과가 m[(d/dt)^2]y(t) 인 것이다.

(d/dt)를 제어에서 미분연산자 s로 표현한다.

예쁘게 바꿔 표현하면
m(s^2)y(t)=-k[y(t)-u(t)]-bs[y(t)-u(t)]
가 되겠다.

Y(s)/U(s)= [k+bs]/[ms^2+k+bs]
되겠다.
물론 여기서 쭉~ 더 풀어서 시스템 방정식을 구해도 되지만 넘어가자...

참고로 u(t)를 0으로 하면 캘빈 방정식이 나오는데.
m(s^2)y(t)=F(t)= -ky(t)-b[dy(t)/dt]가 된다.
위키피디아에서 나오는 식과는 조금 다른데, 아마 방향때문이겠지.

===================================================================
그럼 Maxwell model로 넘어가자. 그림(4)가 되겠다.

문제를 다시한번 상기한다.

문제.
다음 기계시스템의 전달함수를 구하라. 단, 힘f가 입력, 변위 y가 출력이다.
또한, 물체는 수평면상에서 운동한다.

(4)그림에는 안적혀 잇지만 질량 m물체 M밑에 아래화살표로 f가 되어있다.
y가 물체의 좌표이고
댐퍼의 점성마찰계수는 c라 한다.
용수철 상수는 k이다.

내가 잘못했던것부터 차근차근 가겠다.

일단 기본부터가 잘못됐었다.

처음에 생각했던것이... 아주 추상적인 생각이었다.

1)=> 물체 M과 용수철을 한덩어리로 본다. 용수철이 늘어난다고 해도 용수철과 물체를 한덩어리로 보면 댐퍼에 mg가 가해지는것 아니겠는가? 댐퍼와 물체 사이에서 일어나는 갖가지 힘은 결국 그 '계'안에서 상쇄될것이다.
------> 해봤지만 정답과는 아주 거리가 멀었다. 추상적인 생각이었기에 뭘 잘못했는지도 알 수 없었다.

참고로 추상적인 생각만 6시간동안 했다... 책부터 뒤져볼껄 그랬다.

2)=>댐퍼의 변위를 x로 생각한다. 용수철의 변위는 {y(t)-x(t)}가 될것이다.
m(s^2)y=-c[d/dt]x(t)-k[y(t)-x(t)]

----> 쓸데없는 변수를 만들었다. 식이 맞고틀리고는 둘째치고 결과낼때 x(t)를 못없애겠드라...

이 외에도 여러가지 생각해봤지만, 추상적으로 생각해봤자 무슨 결과가 나오겠는가...
결국 인터넷과 책들을 뒤져보았다.

일단 위키피디아에서 나오는 Maxwell material을 참조해보자.

Maxwell material
1. σ_Total = σ_D = σ_S

와 같은 식이 된다. y(t)에 대한 운동방정식이라 하겠다.

뭔가, 캘빈식과 같지 않은가? 분명 캘빈모델과는 다른 모델일 터인데 말이다. 나중에 교수님께 여쭤봐야겠다.

분명 이 식과 관련시켜 모델의 작용을 상상해보면 분명 틀린말은 아닌것 같다.
물론, 이를 확인해보지 않아 증명못하겠지만...

같은 모델일텐데 Maxwell material과는 식이 전혀 다르다.
일단 전제부터가 다르다.

Maxwell material에서는 댐퍼와 스프링에 걸리는 응력이 각각 같다고 나오는데.
여기에서는 각각에 걸리는 힘은 다르고 각각의 힘을 합친것이 힘의 총합이다.

'응력'과 '힘'의 개념을 잘못잡고있는 것일까?

그런데, 책을 읽다가 발견한것... 스프링-댐퍼-질량모델은 RLC모델과 같은 형태를 취한다고 볼 수 있다.

스프링-댐퍼-질량모델과 같은것을 전기회로인 RLC의 형태로 바꾸는것을 '전기적 상사'라고 한다.
--------------------------------------------------------------------------------------------

위와같은 RLC회로로 바꾸는것이 가능하다.

RLC회로에 대해서는 직접 알아보기를 권하지만 일단 간단하게 설명하자면.

각각 R, L, C를 저항과 같이 표현하면

이와같이 표현이 가능하다.

스프링-댐퍼-질량모델과 비교하자면

(1/C)<=>k
L<=>m
R<=>c
q<=>y
E<=>f

여기서 주의,  (1/C)<=>k 와, E<=>f 임을 명심하자
f는 질량 m에 주어지는 임의의 힘이다. 즉, 위에서 논의한 F~F3를 제외하고 임의로 힘이 주어졌을 경우를 이야기 한다.

최종적으로 위의 세 식중 제일 마지막 식이 스프링-댐퍼-질량 모델의 식과 유사한 형태를 띄고 있음을 알 수 있다. 물론 기본적인 식은 제일 위의 첫번째 식이고 이 역시 스프링-댐퍼-질량 모델의 식과 같다.
실제로 쓰이는 식은 첫번째 식이다.

여기서 또 주의. 과연 위의 RLC회로는 스프링-댐퍼-질량의 어떤 형태와 같은것일까?
Kelvin모델일까? Maxwell모델일까?
문제의 해결은 이것에서부터 시작한다.
그럼 문제를 해결하도록 하겠다.
-------------------------------------------------------------------------------------------
일단 결론부터 말하자면 위의 RLC회로는 Kelvin식이라고 할 수 있겠다.
단순하게 이치를 따지면서 생각해 보자.
위의 RLC회로는 '직렬'이다.
RLC회로의 '저항의'계산식은 '직렬식'이다.

그럼, 여기에서 발생할 수 있는 오류.
용수철의 직렬연결=RLC회로의 '저항의'직렬연결.

이는 잘못된 생각이다.

단순하게 생각해보자.
물체를 몇개의 용수철에 연결하여 천장에 매단다고 했을 때.
.용수철을 직렬로 연결할 경우, 이를 하나의 물체로 생각할 수 있다. 이 하나의 물체에 하나의 힘, 물체의 무게가 걸릴것이다.
.용수철을 병렬로 연결할 경우, 이를테면 5개의 용수철에 병렬로 연결하여 천장에 매달경우. 단순한 모델에서는 각각의 용수철에 물체 무게의 1/5씩 걸릴것이다.

전기회로에서는 저항을 직렬로 연결하여야 각각의 저항이 전력E(v)를 나눠가질것이다.
같은크기의 전류 I가 지나지만, 각각의 저항에 걸리는 전압은 다르다.

어쨌든, 위의 두 경우를 보면 용수철의 병렬연결이 전기회로의 저항의 직렬연결과 같다는걸 알 수 있다.

이에서 내가 수시간동안 바보같은 행동을 하고있었다는걸 알았다. 이를 금방 깨닫고 문제를 해결한 다른이들과 비교하면 나는 바보가 아닐까?

어쨌든 이에 의거하여 식을 다시 세워보자.
일단 [d/dt]=s라고 정의 해두겠다.

문제의 시스템은, 용수철과 댐퍼가 '직렬'로 연결되어 있지만, 전기회로와 비교해 보면 '병렬'연결이다.