항공우주/동역학

강체의 운동방정식 - 2

깊은대학 2022. 2. 5. 16:48

관성좌표계의 원점 O 에 대한 파티클 시스템의 총 각운동량 HO 를 다음과 같이 정의한 바 있다.

 

(1)HO=j=1nrj×mjvj

 

여기서 vj 는 파티클 j 의 속도로서 vj=idrjdt 이다.

 

 

임의의 점 A 에 대한 파티클 시스템의 총 각운동량 HA 는 다음과 같이 정의한다.

 

(2)HA=j=1nrj/A×mjvj

 

여기서 rj/A 는 점 A 에서 파티클 j 까지의 위치벡터다. 주의할 점은 속도는 여전히 vj=idrjdt 로서 관성좌표계의 원점 O 에 대한 파티클 j의 위치벡터 rj 를 관성좌표계에서 미분한 속도라는 것이다.

 

 

A 에 대한 파티클 시스템의 총 각운동량 HA 를 다음과 같이 정의하기도 한다.

 

(3)HA=j=1nrj/A×mjvj/A

 

여기서 vj/A 는 점 A 에 대한 파티클 j 의 상대속도로서 vj/A=idrj/Adt 이다. 일반적으로 HAHA 이지만 점 A 가 시스템의 질량중심 G 라면 두 각운동량은 같다.

 

(4)HG=HG

 

먼저 이를 증명해 보기로 하자. 식 (2)로부터 질량중심G에 대한 파티클 시스템의 총 각운동량 HG 는 다음과 같다.

 

(5)HG=j=1nrj/G×mjidrjdt=j=1nrj/G×mjiddt(rG+rj/G)=j=1nrj/G×mjidrGdt+j=1nrj/G×mjidrj/Gdt

 

위 식에서 점 G 는 질량중심이므로,

 

(6)j=1nrj/G×mjidrGdt=(j=1nmjrj/G)×idrGdt=0

 

이 된다. 따라서 식 (5)는 식 (3)의 정의에 의해서

 

(7)HG=j=1nrj/G×mjidrj/Gdt=HG

 

가 된다.

 

 

이제 임의의 점 A 에 대한 시스템의 총 각운동량 식 (2)를 미분해 보자.

 

(8)idHAdt=j=1nidrj/Adt×mjidrjdt+j=1nrj/A×mjid2rjdt2

 

여기서 위 그림에 나와 있듯이 다음과 같은 위치벡터 간의 관계식을 이용하면,

 

(9)rj/A=rG/A+rj/G ,     rj=rG+rj/G

 

식 (8)은 다음과 같이 된다.

 

(10)idHAdt=j=1n(idrG/Adt+idrj/Gdt)×mj(idrGdt+idrj/Gdt)+j=1nrj/A×Fj=midrG/Adt×vG+j=1nMjA

 

여기서 질량중심의 정의 j=1nmjrj/G=0 을 이용했으며, vG 는 질량중심의 속도로서 vG=idrGdt, m 은 총 질량으로서 m=j=1nmj 이고, MjA 는 점 A 에 대한 외력 Fj 에 의한 모멘트이다.

 

 

만약 점 A 가 질량중심이라면 식 (10)에서 rG/A=0 이므로 다음 식이 성립한다.

 

(11)idHGdt=j=1nMjG

 

또한 식 (4)에 의하면 위 식은 다음과 같다.

 

(12)idHGdt=j=1nMjG

 

한편 점 A 가 관성좌표계의 원점 O 라면 식 (10)에서 rG/A=rG 이므로 idrG/Adt×vG=vG×vG=0 이 되어서 다음 식이 성립한다.

 

(13)idHOdt=j=1nMjO

 

따라서 관성좌표계에서 고정된 점 O 또는 움직이는 질량중심점 G 에 대한 총 각운동량의 변화율은 점 O 또는 G에 관한 외력에 의한 모멘트와 같다는 것을 알 수 있다. 이와 같이 질량중심점에 관해서 모멘트 식을 전개하면 고정된 점에 관해서 식을 전개한 것처럼 식이 간단해진다.

파티클 시스템의 총 운동에너지는 각 파티클이 갖는 운동에너지의 합과 같다.

 

(14)T=12j=1nmjidrjdtidrjdt

 

식 (9)를 이용하면 위 식을 다음과 같이 전개할 수 있다.

 

(15)T=12j=1nmj(idrGdt+idrj/Gdt)(idrGdt+idrj/Gdt)=12mvGvG+12j=1nmjidrj/Gdtidrj/Gdt

 

식 (15)에 의하면, 파티클 시스템의 운동에너지는 질량중심에 집중된 시스템의 전체 질량을 고려하여 얻은 운동에너지에 질량중심에 대한 파티클의 상대적인 운동에너지를 더한 것과 같다.

이로서 강체 운동방정식을 구하기 위한 준비를 모두 마쳤다.