7.35 극좌표, 원통 좌표, 구면 좌표 변환

로봇공학에서 다루는 많은 물리 시스템은 직교 좌표계(Cartesian coordinate system)보다 극좌표(polar coordinates), 원통 좌표(cylindrical coordinates), 구면 좌표(spherical coordinates)로 표현할 때 수학적으로 더 간결하고 계산이 효율적이다. 회전 관절을 가진 매니퓰레이터의 작업 공간, 원통형 또는 구형 링크의 물성 계산 등에서 좌표 변환은 필수적이다.

1. 극좌표계 (Polar Coordinates)

1.1 정의와 변환 관계

2차원 평면에서 점 P의 위치를 원점으로부터의 거리 r과 양의 x축으로부터 반시계 방향으로 측정한 각도 \theta로 표현한 것이 극좌표이다.

직교 좌표 (x, y)와 극좌표 (r, \theta) 사이의 변환 관계는 다음과 같다.

x = r\cos\theta, \quad y = r\sin\theta

역변환은 다음과 같다.

r = \sqrt{x^2 + y^2}, \quad \theta = \text{atan2}(y, x)

여기서 \text{atan2}(y, x)는 4사분면을 올바르게 구분하는 역삼각 함수이다. r \geq 0이고 \theta \in [0, 2\pi) 또는 \theta \in (-\pi, \pi]로 제한한다.

1.2 극좌표의 면적 요소

극좌표에서 미소 면적 요소는 다음과 같다.

dA = r \, dr \, d\theta

이는 반경 방향 폭 dr과 원호 방향 폭 r \, d\theta로 구성되는 미소 직사각형의 넓이에 해당한다.

1.3 극좌표에서의 이중 적분

영역 D가 극좌표로 \alpha \leq \theta \leq \beta, h_1(\theta) \leq r \leq h_2(\theta)로 표현되면 다음과 같다.

\iint_D f(x, y) \, dA = \int_\alpha^\beta \int_{h_1(\theta)}^{h_2(\theta)} f(r\cos\theta, r\sin\theta) \, r \, dr \, d\theta

1.4 로봇공학 응용: 2자유도 평면 매니퓰레이터

링크 길이가 l_1, l_2인 2자유도 평면 매니퓰레이터의 작업 공간은 원환(annulus) 형태이다. 이때 작업 공간의 면적은 극좌표를 이용하면 쉽게 계산된다.

A = \int_0^{2\pi} \int_{\vert l_1 - l_2 \vert}^{l_1 + l_2} r \, dr \, d\theta = 2\pi \cdot \frac{1}{2}\left[(l_1 + l_2)^2 - (l_1 - l_2)^2\right] = 4\pi l_1 l_2

2. 원통 좌표계 (Cylindrical Coordinates)

2.1 정의와 변환 관계

3차원 공간에서 원통 좌표 (r, \theta, z)는 xy-평면에서의 극좌표 (r, \theta)에 높이 z를 추가한 것이다.

x = r\cos\theta, \quad y = r\sin\theta, \quad z = z

역변환은 다음과 같다.

r = \sqrt{x^2 + y^2}, \quad \theta = \text{atan2}(y, x), \quad z = z

매개변수의 범위는 r \geq 0, \theta \in [0, 2\pi), z \in (-\infty, \infty)이다.

2.2 원통 좌표의 체적 요소

원통 좌표에서 미소 체적 요소는 다음과 같다.

dV = r \, dr \, d\theta \, dz

이는 세 방향의 미소 길이 dr, r \, d\theta, dz의 곱으로부터 유도된다.

2.3 원통 좌표에서의 삼중 적분

영역 E가 원통 좌표로 표현될 때, 삼중 적분은 다음과 같이 변환된다.

\iiint_E f(x, y, z) \, dV = \iiint_{E'} f(r\cos\theta, r\sin\theta, z) \, r \, dr \, d\theta \, dz

여기서 E'는 (r, \theta, z) 공간에서의 대응 영역이다.

2.4 원통형 물체의 질량 계산 예제

밀도가 \rho(r, z) = \rho_0(1 - z/L)인 원통형 링크(반지름 R, 높이 L)의 질량을 계산한다.

m = \int_0^{2\pi} \int_0^R \int_0^L \rho_0\left(1 - \frac{z}{L}\right) r \, dz \, dr \, d\theta

= 2\pi \rho_0 \int_0^R r \, dr \int_0^L \left(1 - \frac{z}{L}\right) dz = 2\pi \rho_0 \cdot \frac{R^2}{2} \cdot \frac{L}{2} = \frac{\pi \rho_0 R^2 L}{2}

3. 구면 좌표계 (Spherical Coordinates)

3.1 정의와 변환 관계

3차원 공간에서 구면 좌표 (\rho, \theta, \phi)는 원점으로부터의 거리 \rho, 양의 z축으로부터의 천정각(polar angle) \phi, 그리고 xy-평면에서의 방위각(azimuthal angle) \theta로 위치를 표현한다.

x = \rho \sin\phi \cos\theta, \quad y = \rho \sin\phi \sin\theta, \quad z = \rho \cos\phi

역변환은 다음과 같다.

\rho = \sqrt{x^2 + y^2 + z^2}, \quad \phi = \arccos\left(\frac{z}{\rho}\right), \quad \theta = \text{atan2}(y, x)

매개변수의 범위는 \rho \geq 0, 0 \leq \phi \leq \pi, 0 \leq \theta < 2\pi이다.

3.2 구면 좌표의 체적 요소

구면 좌표에서 미소 체적 요소는 다음과 같다.

dV = \rho^2 \sin\phi \, d\rho \, d\phi \, d\theta

이 체적 요소의 유도 과정은 다음과 같다. 세 방향의 미소 길이는 각각 d\rho, \rho \, d\phi, \rho \sin\phi \, d\theta이며, 이들의 곱이 미소 체적이 된다.

3.3 구면 좌표에서의 삼중 적분

\iiint_E f(x, y, z) \, dV = \iiint_{E'} f(\rho\sin\phi\cos\theta, \rho\sin\phi\sin\theta, \rho\cos\phi) \, \rho^2 \sin\phi \, d\rho \, d\phi \, d\theta

3.4 구의 관성 모멘트 계산

균일 밀도 \rho_0인 반지름 R의 구에 대해 z축 둘레의 관성 모멘트를 구면 좌표로 계산한다. z축으로부터의 거리는 d_\perp = \rho \sin\phi이므로 다음과 같다.

I_z = \iiint_E \rho_0 d_\perp^2 \, dV = \rho_0 \int_0^{2\pi} \int_0^\pi \int_0^R (\rho\sin\phi)^2 \cdot \rho^2 \sin\phi \, d\rho \, d\phi \, d\theta

= \rho_0 \int_0^{2\pi} d\theta \int_0^\pi \sin^3\phi \, d\phi \int_0^R \rho^4 \, d\rho

각 적분을 계산하면 다음과 같다.

\int_0^{2\pi} d\theta = 2\pi, \quad \int_0^\pi \sin^3\phi \, d\phi = \frac{4}{3}, \quad \int_0^R \rho^4 \, d\rho = \frac{R^5}{5}

따라서 다음을 얻는다.

I_z = \rho_0 \cdot 2\pi \cdot \frac{4}{3} \cdot \frac{R^5}{5} = \frac{8\pi \rho_0 R^5}{15}

질량 m = \frac{4}{3}\pi \rho_0 R^3을 대입하면 다음과 같다.

I_z = \frac{2}{5}mR^2

4. 좌표계 비교와 선택 기준

5. 좌표계 간 변환 관계 요약

원통 좌표와 구면 좌표 사이의 관계는 다음과 같다.

r = \rho \sin\phi, \quad z = \rho \cos\phi, \quad \theta = \theta

\rho = \sqrt{r^2 + z^2}, \quad \phi = \arctan\left(\frac{r}{z}\right), \quad \theta = \theta

직교 좌표에서 원통 좌표로의 변환 행렬은 다음과 같다.

\begin{pmatrix} \hat{\mathbf{e}}_r \\ \hat{\mathbf{e}}_\theta \\ \hat{\mathbf{e}}_z \end{pmatrix} = \begin{pmatrix} \cos\theta & \sin\theta & 0 \\ -\sin\theta & \cos\theta & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} \hat{\mathbf{e}}_x \\ \hat{\mathbf{e}}_y \\ \hat{\mathbf{e}}_z \end{pmatrix}

직교 좌표에서 구면 좌표로의 변환 행렬은 다음과 같다.

\begin{pmatrix} \hat{\mathbf{e}}_\rho \\ \hat{\mathbf{e}}_\phi \\ \hat{\mathbf{e}}_\theta \end{pmatrix} = \begin{pmatrix} \sin\phi\cos\theta & \sin\phi\sin\theta & \cos\phi \\ \cos\phi\cos\theta & \cos\phi\sin\theta & -\sin\phi \\ -\sin\theta & \cos\theta & 0 \end{pmatrix} \begin{pmatrix} \hat{\mathbf{e}}_x \\ \hat{\mathbf{e}}_y \\ \hat{\mathbf{e}}_z \end{pmatrix}

이 변환 행렬들은 모두 직교 행렬(orthogonal matrix)이므로, 역변환은 전치 행렬로 구해진다.

6. 로봇공학에서의 좌표 변환 활용

로봇의 구조에 따라 적절한 좌표계를 선택하는 것이 중요하다.

• 원통형 로봇 (cylindrical robot): 작업 공간이 원통 형태이므로, 원통 좌표계를 사용하면 순기구학과 역기구학 계산이 간결해진다.
• 구형 로봇 (spherical robot): 작업 공간이 구의 일부이므로, 구면 좌표계가 자연스럽다.
• SCARA 로봇: xy-평면에서의 운동은 극좌표, z-방향은 직선 운동이므로 원통 좌표계와 유사한 구조를 가진다.

각 좌표계에서의 속도 벡터 표현도 중요하다. 극좌표에서 속도 벡터는 다음과 같다.

\mathbf{v} = \dot{r}\hat{\mathbf{e}}_r + r\dot{\theta}\hat{\mathbf{e}}_\theta

원통 좌표에서는 다음과 같다.

\mathbf{v} = \dot{r}\hat{\mathbf{e}}_r + r\dot{\theta}\hat{\mathbf{e}}_\theta + \dot{z}\hat{\mathbf{e}}_z

구면 좌표에서는 다음과 같다.

\mathbf{v} = \dot{\rho}\hat{\mathbf{e}}_\rho + \rho\dot{\phi}\hat{\mathbf{e}}_\phi + \rho\sin\phi\dot{\theta}\hat{\mathbf{e}}_\theta

참고 문헌

• Stewart, J. (2015). Calculus: Early Transcendentals. 8th ed., Cengage Learning.
• Kreyszig, E. (2011). Advanced Engineering Mathematics. 10th ed., Wiley.
• Craig, J. J. (2005). Introduction to Robotics: Mechanics and Control. 3rd ed., Pearson.
• Marsden, J. E., & Tromba, A. J. (2012). Vector Calculus. 6th ed., W. H. Freeman.

v 0.1