13.4 질량 분포와 밀도 함수

1. 질량 분포의 개요

질량 분포(mass distribution)는 강체 내부의 질량이 공간적으로 어떻게 분포되어 있는지를 나타내는 양이며, 강체의 동역학적 성질(질량 중심, 관성 텐서, 운동 에너지 등)을 결정하는 가장 기본적인 정보이다. 질량 분포는 일반적으로 밀도 함수(density function)로 기술된다. 본 절에서는 밀도 함수의 정의, 종류, 그리고 강체 역학에서의 활용을 자세히 다룬다.

2. 밀도 함수의 정의

2.1 부피 밀도

3차원 강체의 부피 밀도(volume density) $\rho(\mathbf{r})$ 은 위치 $\mathbf{r}$ 에서의 단위 부피당 질량이다.

$\rho(\mathbf{r}) = \lim_{\Delta V \to 0}\frac{\Delta m}{\Delta V}$

여기서 $\Delta m$ 은 부피 $\Delta V$ 내의 질량이다.

2.2 단위

부피 밀도의 단위는 kg/m³ (SI 단위계)이다.

2.3 밀도와 질량의 관계

전체 질량은 밀도의 부피 적분이다.

$M = \int_B\rho(\mathbf{r})\,dV$

여기서 $B$ 는 강체의 영역이다.

3. 다른 종류의 밀도

3.1 면 밀도

얇은 평판이나 표면의 면 밀도(surface density 또는 areal density) $\sigma(\mathbf{r})$ 은 단위 면적당 질량이다.

$\sigma(\mathbf{r}) = \frac{dm}{dA}$

단위는 kg/m²이다.

3.2 선 밀도

가는 봉이나 줄의 선 밀도(linear density) $\lambda(\mathbf{r})$ 은 단위 길이당 질량이다.

$\lambda(\mathbf{r}) = \frac{dm}{dL}$

단위는 kg/m이다.

4. 균일 밀도와 비균일 밀도

4.1 균일 밀도

균일 밀도(uniform density)는 강체 내부의 모든 점에서 밀도가 같은 경우이다.

$\rho(\mathbf{r}) = \rho_0 = \text{const}$

이 경우 질량 분포가 단순하며, 적분 계산이 쉽다.

4.2 비균일 밀도

비균일 밀도(non-uniform density)는 위치에 따라 밀도가 변하는 경우이다. 다양한 함수 형태가 가능하다.

선형 변화: $\rho(x) = \rho_0 + ax$
지수 변화: $\rho(x) = \rho_0 e^{-x/L}$
다항식: $\rho(\mathbf{r}) = a + bx + cy + dz + \ldots$

4.3 응용

비균일 밀도는 합성 재료, 다층 구조, 유체 충전 컨테이너 등에서 발생한다.

5. 점 질량의 디랙 델타 함수 표현

5.1 디랙 델타 함수

디랙 델타 함수(Dirac delta function) $\delta(\mathbf{r})$ 은 다음의 성질을 가진 분포이다.

$\int_{\mathbb{R}^3}\delta(\mathbf{r})\,dV = 1$

$\int_{\mathbb{R}^3}f(\mathbf{r})\delta(\mathbf{r} - \mathbf{r}_0)\,dV = f(\mathbf{r}_0)$

5.2 점 질량의 표현

위치 $\mathbf{r}_i$ 에 있는 질량 $m_i$ 의 점 질량은 다음과 같이 밀도 함수로 표현될 수 있다.

$\rho(\mathbf{r}) = m_i\delta(\mathbf{r} - \mathbf{r}_i)$

5.3 다중 점 질량

여러 점 질량의 시스템은 다음과 같이 표현된다.

$\rho(\mathbf{r}) = \sum_{i=1}^{N}m_i\delta(\mathbf{r} - \mathbf{r}_i)$

5.4 통합

이러한 표현으로 점 질량과 연속체를 통합적으로 다룰 수 있다.

6. 밀도와 다른 물리량

6.1 질량 중심

질량 중심은 밀도의 1차 모멘트이다.

$\mathbf{r}_c = \frac{1}{M}\int_B\rho(\mathbf{r})\mathbf{r}\,dV$

6.2 관성 텐서

관성 텐서는 밀도의 2차 모멘트와 관련된다.

$I_{ij} = \int_B\rho(\mathbf{r})(r^2\delta_{ij} - r_ir_j)\,dV$

6.3 전체 질량

$M = \int_B\rho(\mathbf{r})\,dV$

7. 균일 밀도 강체의 질량과 부피

7.1 단순 도형

균일 밀도 강체의 질량은 단순히 밀도와 부피의 곱이다.

$M = \rho_0 V$

7.1.1 직육면체

부피: $V = abc$ , 질량: $M = \rho_0 abc$

7.1.2 구

부피: $V = \frac{4}{3}\pi R^3$ , 질량: $M = \frac{4}{3}\pi R^3\rho_0$

7.1.3 원기둥

부피: $V = \pi R^2 h$ , 질량: $M = \pi R^2 h\rho_0$

7.1.4 원뿔

부피: $V = \frac{1}{3}\pi R^2 h$ , 질량: $M = \frac{1}{3}\pi R^2 h\rho_0$

8. 비균일 밀도의 예시

8.1 지구

지구의 밀도는 균일하지 않다. 핵에서 약 13,000 kg/m³, 표면 근처에서 약 3,000 kg/m³이다. 평균 밀도는 약 5,515 kg/m³이다.

8.2 합성 재료

합성 재료(섬유 강화 플라스틱, 적층판 등)는 비균일 밀도를 가진다.

8.3 다공성 재료

다공성 재료는 빈 공간이 있어 평균 밀도가 낮다. 빈 공간의 분포에 따라 밀도가 변한다.

9. 밀도의 좌표계 변환

9.1 밀도의 변환

밀도 함수는 좌표계 변환에 대해 다음과 같이 변환된다.

$\rho'(\mathbf{r}') = \rho(\mathbf{r})$

여기서 $\mathbf{r}' = \mathbf{R}\mathbf{r} + \mathbf{t}$ 이다. 밀도는 스칼라이므로 좌표 변환에 대해 불변이다.

9.2 본체 좌표계

본체 좌표계에서 밀도는 시간에 따라 변하지 않는다. 강체의 정의에 따라 질량 분포가 고정되어 있기 때문이다.

10. 적분 계산

10.1 직교 좌표계

직교 좌표계에서 부피 적분은 다음과 같이 표현된다.

$\int_B f(\mathbf{r})\rho(\mathbf{r})\,dV = \int_{x_1}^{x_2}\int_{y_1}^{y_2}\int_{z_1}^{z_2}f(x, y, z)\rho(x, y, z)\,dz\,dy\,dx$

10.2 원통 좌표계

원기둥 대칭 강체에서는 원통 좌표계가 편리하다.

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

10.3 구 좌표계

구 대칭 강체에서는 구 좌표계가 편리하다.

$dV = r^2\sin\theta\,dr\,d\theta\,d\phi$

10.4 좌표계의 선택

응용에 따라 적절한 좌표계를 선택하면 적분 계산이 크게 단순화된다.

11. 수치 적분

11.1 메쉬 기반 계산

복잡한 형상의 강체에서는 수치적 적분이 필요하다. CAD 모델의 메쉬 데이터로부터 부피, 질량, 질량 중심 등을 계산한다.

11.2 Monte Carlo 방법

Monte Carlo 방법은 무작위 점을 사용하여 적분을 근사한다. 매우 복잡한 형상에 적합하다.

11.3 가우스 적분

가우스 적분(Gaussian quadrature)은 특정 점들에서의 함수값을 가중합으로 적분을 근사한다. 효율적이고 정확하다.

12. 밀도와 동역학

12.1 운동 에너지

강체의 운동 에너지는 밀도 적분으로 표현된다.

$T = \frac{1}{2}\int_B\rho(\mathbf{r})\lVert\mathbf{v}(\mathbf{r})\rVert^2\,dV$

여기서 $\mathbf{v}(\mathbf{r})$ 은 위치 $\mathbf{r}$ 에서의 속도이다.

12.2 각운동량

각운동량도 밀도 적분으로 표현된다.

$\mathbf{L} = \int_B\rho(\mathbf{r})(\mathbf{r} \times \mathbf{v}(\mathbf{r}))\,dV$

12.3 강체의 운동 방정식

이러한 적분 표현이 강체의 운동 방정식의 유도에 사용된다.

13. 매니퓰레이터 링크의 밀도

13.1 모델링

매니퓰레이터의 각 링크는 일반적으로 균일 밀도로 모델링된다. 이는 단순화이지만 합리적인 근사이다.

13.2 정확한 모델

CAD 데이터를 사용하면 더 정확한 비균일 밀도 모델이 가능하다. 복잡한 내부 구조(모터, 와이어, 강화 리브 등)를 반영할 수 있다.

13.3 매개변수 식별

링크의 밀도 분포는 일반적으로 알려지지 않으며, 매개변수 식별 실험을 통해 추정된다. 질량, 질량 중심, 관성 텐서가 식별의 대상이다.

14. 응용 예시: 균일 막대

길이 $L$ , 단면적 $A$ 인 균일 막대를 고려하자. 선 밀도는 다음과 같다.

$\lambda = \rho A$

총 질량은 $M = \rho A L = \lambda L$ 이다.

질량 중심은 막대의 중간이다.

15. 응용 예시: 비균일 막대

선 밀도가 $\lambda(x) = \lambda_0(1 + x/L)$ 인 막대를 고려하자.

총 질량은

$M = \int_0^L\lambda_0\!\left(1 + \frac{x}{L}\right)dx = \lambda_0 L\!\left(1 + \frac{1}{2}\right) = \frac{3}{2}\lambda_0 L$

질량 중심은

$x_c = \frac{1}{M}\int_0^L x\lambda_0\!\left(1 + \frac{x}{L}\right)dx = \frac{\lambda_0}{M}\!\left[\frac{x^2}{2} + \frac{x^3}{3L}\right]_0^L = \frac{\lambda_0}{M} \cdot \frac{5L^2}{6}$

$M$ 을 대입하면 $x_c = 5L/9$ 이다. 이는 막대의 중간보다 무거운 끝쪽으로 치우쳐 있다.

16. 응용 예시: 균일 직사각형 판

가로 $a$ , 세로 $b$ , 두께 $t$ 인 균일 직사각형 판의 부피는 $V = abt$ 이고, 질량은 $M = \rho V$ 이다.

면 밀도(두께를 따라 적분)는 $\sigma = \rho t$ 이다.

17. 응용 예시: 구의 비균일 밀도

반지름 $R$ 인 구의 밀도가 중심에서 표면으로 갈수록 줄어든다고 하자.

$\rho(r) = \rho_0\!\left(1 - \frac{r}{R}\right)$

총 질량은

$M = \int_0^R\rho(r)4\pi r^2\,dr = 4\pi\rho_0\int_0^R\!\left(r^2 - \frac{r^3}{R}\right)dr$

$= 4\pi\rho_0\!\left[\frac{R^3}{3} - \frac{R^3}{4}\right] = 4\pi\rho_0 \cdot \frac{R^3}{12} = \frac{\pi\rho_0 R^3}{3}$

18. 응용 예시: 매니퓰레이터 링크

매니퓰레이터의 한 링크가 균일 밀도 원기둥으로 근사된다. 길이 $L$ , 반지름 $R$ , 밀도 $\rho$ 일 때

질량: $M = \rho\pi R^2 L$
질량 중심: 축의 중간
관성 모멘트: 별도의 절에서 다룬다

19. 응용 예시: 복합 부품

매니퓰레이터의 손목과 같은 복합 부품은 여러 균일 밀도 부품으로 분해하여 분석한다. 각 부품의 질량과 질량 중심을 계산한 후 합성한다.

20. 응용 예시: 드론

드론의 질량 분포는 매우 비균일하다. 중심에 무거운 배터리가 있고, 가장자리에 모터가 있다. 정확한 질량 분포 모델이 비행 안정성 분석에 필요하다.

21. 응용 예시: 인간형 로봇

인간형 로봇의 각 부위(머리, 몸통, 팔, 다리)는 다른 질량과 밀도를 가진다. 전체 시스템의 동역학 분석에는 이러한 분포가 필수이다.

22. 밀도 측정

22.1 직접 측정

부피와 질량을 직접 측정하여 평균 밀도를 계산할 수 있다.

$\bar\rho = \frac{M}{V}$

22.2 비균일 밀도의 측정

비균일 밀도는 측정하기 어렵다. CT 스캔, X선 단층 촬영 등이 사용된다.

22.3 시뮬레이션 검증

시뮬레이션 결과와 실제 측정의 비교를 통해 밀도 모델을 검증한다.

23. 밀도 함수의 일반화

23.1 시변 밀도

일반적으로 강체의 밀도는 시간에 따라 변하지 않지만, 변형체나 유체는 시변 밀도를 가질 수 있다.

$\rho(\mathbf{r}, t)$

23.2 응용

연속체 역학과 유체역학에서 시변 밀도가 중요하다.

23.3 강체의 한계

강체 가정 하에서는 밀도가 본체 좌표계에서 시간에 무관하다. 이는 강체 모형의 본질적 단순화이다.

24. 학습 권장사항

다양한 종류의 밀도(부피, 면, 선)를 이해한다.
균일 밀도와 비균일 밀도의 차이를 인식한다.
좌표계 선택의 중요성을 학습한다.
적분 계산을 연습한다.
디랙 델타 함수의 활용을 익힌다.

25. 본 절의 의의

본 절은 강체의 질량 분포와 밀도 함수를 자세히 다루었다. 밀도 함수는 강체의 모든 동역학적 성질을 결정하는 기본 정보이며, 후속 절에서 다룰 관성 텐서와 운동 방정식의 토대이다. 다양한 강체의 정확한 모델링과 분석에서 밀도의 정확한 표현이 필수이다.

26. 참고 문헌

Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
Symon, K. R. (1971). Mechanics (3rd ed.). Addison-Wesley.
Beer, F. P., Johnston Jr., E. R., Mazurek, D. F., & Cornwell, P. J. (2013). Vector Mechanics for Engineers: Statics and Dynamics (10th ed.). McGraw-Hill.
Hibbeler, R. C. (2016). Engineering Mechanics: Dynamics (14th ed.). Pearson.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.

version: 1.0