질량(선밀도: $\lambda$)이 있는 줄이 매달린 공($M$)을 위로 $v_0$ 쏘아 올린다. 최고점의 높이는 공의 질량에 (줄의 질량이 없는 경우는 최고점은 공의 질량에 무관하다)

1. 무관하다.

2. 의존한다.

더보기

떠있는 공+줄 부분에 작용하는 힘:  공의 무게와 줄의 무게, 그리고 바닥의 정지해 있는 줄이 뜨기 위해서는 impulsive force가 필요하다. 만약 질량 $dm$이 속도 $0\rightarrow v$로 변하면 운동량 변화가 $dp = dm(v-0)$이므로 필요한 충격력은 $f = v \frac{dm}{dt}$(위쪽)이다. 따라서 공중에 떠 있는 부분이 받는 반작용(장력) 충격력은 $-f $이다. 

공의 높이가 $y$일 때 공중에 떠 있는 부분에 뉴턴의 운동방정식을 쓰면(위쪽:+)

$$ (M+ \lambda y) \frac{dv}{dt}=  \sum F = -(M + \lambda y) g -f $$ 

운동량이 $P= (M+ \lambda y)v$이므로 운동 방정식은 다시 

$$ \frac{dP}{dt} = -(M+ \lambda y ) g$$

로 쓸 수 있다. (Note: 이는 이미 공중에 있는 부분과 추가되는 $dm$ 부분을 하나의 계로 보므로 충격력은 내력이 되어서 합력에 나타나지 않는다)

 

$P$를 높이 $y$의 함수로 볼 수 있으므로 $\frac{dP}{dt} = \frac{dP}{dy}\frac{dy}{dt} =   \frac{dP}{dy} v = \frac{P}{M + \lambda y}\frac{ dP}{dy}$

$$ P \frac{dP}{dy} = -(M+ \lambda y)^2 g$$

적분을 해서,

$$\frac{1}{2}( P^2 - M^2 v_0^2 ) = \frac{g}{3\lambda} (M^3 - (M+\lambda y)^3)$$ 

따라서 최고점의 높이는 ($P=0$)

$$  y_\text{max} = \frac{M}{\lambda} \left( \sqrt[3]{ 1 + \frac{3\lambda v_0^2}{2Mg}}-1\right)$$

따라서 초기속도뿐만 아니라 공의 질량에도 의존하게 된다. 줄의 질량을 무시할 수 있는 경우인 $\lambda\rightarrow 0$이면

$$y_\text{max} = \frac{M}{\lambda}\left(   1 + \frac{\lambda v_0^2}{2Mg} +... - 1\right) \rightarrow \frac{v_0^2}{2g}$$ 이어서 잘 알고 있는 결과로 수렴한다.

 

그리고 impulsive force 때문에 역학적 에너지는 보존이 되지 않음을 체크할 수 있다.

참고: https://kipl.tistory.com/119

 

떨어지는 쇠사슬이 바닥에 주는 힘

사슬(길이=$L$, 질량=$M$)이 수직으로 바닥으로 떨어질 때 바닥이 받는 힘은? 우선 떨어지는 부분은 사슬고리 사이의 마찰 등을 무시하면 오직 중력에 의해서 자유낙하한다. 사슬이 바닥에 닿을 때

kipl.tistory.com

 

728x90

'Physics > 역학' 카테고리의 다른 글

포탄이 날아가는 각도는?  (0) 2022.01.24
개는 얼마나 가속할 수 있을까?  (0) 2022.01.21
등시진자  (0) 2022.01.14
체인이 떨어지는 가속도는?  (0) 2022.01.12
더 멀리 날아갈까?  (0) 2022.01.12
Posted by helloktk
,

표준형 타원

$$ \frac{x^2}{a^2} + \frac{y^2}{b^2} = 1$$

의 두 주축에 대한 central moment (moment of inertia)는 각각

$$ \mu_{20} = \frac{\pi}{4} a^3 b, \quad \mu_{02} = \frac{\pi}{4} ab^3$$ 

으로 주어짐은 쉽게 계산할 수 있다. 그런데 타원의 면적인 0차 central moment가 $\mu_{00}= \pi ab$이므로 normalized central moment (위치의 분산을 의미한다)는 각각

$$\tilde\mu_{20} = \frac{\mu_{20}}{\mu_{00}} = \frac{1}{4}a^2, \quad  \tilde\mu_{02} = \frac{\mu_{02}}{\mu_{00}} = \frac{1}{4}b^2$$

따라서 장축과 단축의 반지름은 주축에 대한 2차 normalized central moment을 구하면 얻을 수 있다.

기울어진 타원의 경우는 $\tilde\mu_{pq}$가 픽셀의 분포를 알려주므로 이를 이용한 covariance 행렬의 고유값을 구하면 장축과 단축의 반지름을 알 수 있고, 기울어진 정도도 알 수 있다. 공변행렬

$$ \Sigma = \left( \begin{array}{cc} \tilde\mu_{20} & \tilde\mu_{11} \\ \tilde\mu_{11} & \tilde\mu_{02} \end{array} \right)$$

의 두 고윳값은

$$ \lambda_1 = \frac{1}{2}\left( \tilde\mu_{20}+ \tilde\mu_{02} + \sqrt{ (\tilde\mu_{20} - \tilde\mu_{02})^2 + 4 \tilde\mu_{11}^2} \right),$$

$$ \lambda_2 = \frac{1}{2}\left( \tilde\mu_{20}+ \tilde\mu_{02} - \sqrt{ (\tilde\mu_{20} - \tilde\mu_{02})^2 + 4 \tilde\mu_{11}^2} \right).$$

따라서 기울어진 타원의 장축과 단축의 반지름, 그리고 회전각은

$$ a = 2\sqrt{ \lambda_1}, \quad b = 2\sqrt{ \lambda_2}, \quad \tan(2\theta)=\frac{2 \tilde\mu_{11}}{ \tilde\mu_{20}-\tilde\mu_{02}   }$$

로 주어지므로 회전된 좌표계 $(u, v)$에서 타원의 방정식은

$$ \frac{ u^2 }{4 \lambda_1} + \frac{v^2 }{4\lambda_2} = 1$$

 

로 쓰인다. 

다시 원 좌표계 $(x, y)$로 돌아가기 위해서 (물론 질량중심이 원점인 좌표계이다)

$$ u =  \cos \theta x+  \sin \theta y, \quad v = -\sin \theta x + \cos \theta y$$

를 대입하면 타원을 2차식 형태

$$ A x^2 + 2B xy + C y^2 = 1$$

로 쓸 수 있는데, 그 계수는

$$ A= \frac{1}{8}\left( \frac{1}{\lambda_1} + \frac{1}{\lambda_2} +\Big( \frac{1}{\lambda_1} - \frac{1}{\lambda_2} \Big) \cos (2\theta) \right)=\frac{\tilde\mu_{02}}{ 4(\tilde\mu_{20}\tilde\mu_{02}-\tilde\mu_{11}^2 )} $$

$$ C= \frac{1}{8}\left( \frac{1}{\lambda_1} + \frac{1}{\lambda_2} -\Big( \frac{1}{\lambda_1} - \frac{1}{\lambda_2} \Big) \cos (2\theta) \right)=\frac{\tilde\mu_{20}}{ 4(\tilde\mu_{20}\tilde\mu_{02}-\tilde\mu_{11}^2 )} $$

$$ B= \frac{1}{8}\left( \frac{1}{\lambda_1} - \frac{1}{\lambda_2} \right)\sin (2\theta)=\frac{-\tilde\mu_{11}}{ 4(\tilde\mu_{20}\tilde\mu_{02}-\tilde\mu_{11}^2 )} $$

로 주어진다. 주어진 object의 normalized central moment를 구하면 타원 fitting에 대한 정보를 구할 수 있음을 보였다.

 

 

 
 
728x90

'Image Recognition > Fundamental' 카테고리의 다른 글

Least Squares Fitting of Ellipses  (0) 2022.01.27
Circle Fitting: Pratt  (0) 2022.01.20
Best-fit Ellipse  (0) 2022.01.16
Image Moment  (0) 2021.12.04
Orientation 추정  (0) 2021.11.30
Posted by helloktk
,

영상에 담고 있는 object을 간단히 근사를 할 때 타원으로 많이 기술한다(e.g: head tracking). 타원으로 기술하면 장축의 방향으로 object의 기울어진 방향을, 장축과 단축의 길이로 object의 크기를 가늠할 수 있다. object의 픽셀 분포에서 형상에 대한 정보는 2차 moment를 계산해서 얻을 수 있다. 이는 타원이 2차 곡선이기 때문에 가능하다. 그리고 질량중심을 원점으로 잡으면 2차 central moment를 계산해야 한다. 그런데 통계적인 의미를 부여하기 위해서는 central moment를 object의 픽셀로 나눈 normalized central moment로 구성한 covariance matrix를 사용하면 된다.

$$\Sigma = \left(\begin{array}{cc} \tilde\mu_{20} &  \tilde\mu_{11} \\ \tilde\mu_{11} & \tilde\mu_{02} \end{array} \right), \quad\quad\tilde\mu_{pq} \equiv \frac{\mu_{pq}}{\mu_{00}},  \quad (p+q=2)$$ 

(Note: $\mu_{00}$로 정규화를 하지 않더라도 문제는 없다).

$\Sigma$는 영상에서 object pixel의 $x$ 뱡향 분산($\mu_{20}$: $x$축에 대한 회전관성), $y$ 방향 분산($\mu_{02}$: $y$축에 대한 회전관성), $x$-$y$의 correlation을 나타낸다. $\Sigma$가 대칭행렬이므로 두 개의 음이 아닌 고윳값을 가진다. 

\begin{gather} \lambda_1 = \frac{\tilde\mu_{20} + \tilde\mu_{02}}{2} + \frac{ \sqrt{(\tilde\mu_{20} - \tilde\mu_{02})^2 + 4\tilde\mu_{11}^2 } }{2}, \\  \lambda_2 = \frac{\tilde\mu_{20} + \tilde\mu_{02}}{2} - \frac{ \sqrt{(\tilde\mu_{20} - \tilde\mu_{02})^2 + 4\tilde\mu_{11}^2 } }{2}\end{gather}

큰 고유값에 해당하는 고유벡터의 방향이 타원의 장축 방향에 해당하고 (픽셀 변동이 심하므로) , 작은 고윳값의 고유벡터 방향은 단축 방향이다. 그리고 고윳값은 각각 장축과 단축의 반지름의 제곱에 비례한다($\tilde\mu_{pq}$는 단위가 거리 제곱이다). object의 orientation인 타원의 장축 방향은 

$$ \theta = \frac{1}{2} \tan^{-1} \Big( \frac{   2\tilde\mu_{11} }{\tilde\mu_{20} - \tilde\mu_{02}}\Big)= \frac{1}{2} \tan^{-1} \Big( \frac{   2\mu_{11} }{\mu_{20} - \mu_{02}}\Big) $$

로 계산된다. (https://kipl.tistory.com/58)

 

타원의 orientation 각도를 구했으므로 두 주축을 나타내는 단위벡터는

$$\text{major axis: }(\cos \theta, \sin \theta), \quad \text{minor axis: }(-\sin\theta, \cos \theta)$$

로 쓸 수 있다. 그리고 이 두 축에 대한 object의 회전관성은 정의에 의해서 다음 식으로 구할 수 있다: $(\bar{x}, \bar{y})=\text{center of mass}$

\begin{align}\text{major axis:} ~I_\text{min} &= \sum_{(x, y )\in \text{object}} | - (x-\bar{x})\sin \theta + (y - \bar{y})\cos \theta |^2 \\ &=\frac{\mu_{20}+\mu_{02}}{2} -\frac{\mu_{20} - \mu_{02}}{2} \cos (2\theta)-  \mu_{11} \sin(2\theta)\end{align}

\begin{align}\text{minor axis:} ~I_\text{max} &= \sum_{(x, y)\in \text{object}} | (x-\bar{x}) \cos \theta + (y - \bar{y}) \sin \theta|^2 \\&= \frac{\mu_{20}+\mu_{02}}{2}+\frac{\mu_{20}-\mu_{02}}{2}\cos (2\theta)+ \mu_{11}\sin(2\theta)\end{align}

(note: object의 orientation 각은 $I_\text{min}$을 최소화시키는 값이다)

 

표준 타원의 장축 반지름이 $a$고 단축 반지름이 $b$일 때 ($x^2/a^2 + y^2/b^2 =1$) 2차 central moment(회전관성)는 간단한 계산에 의해서

$$\mu_{20}^\text{(e)} = \frac{\pi}{4} a^3 b, \quad \mu_{02}^\text{(e)} = \frac{\pi}{4} ab^3,\quad \mu_{11}^\text{(e)}= 0$$

으로 구해짐을 알 수 있다.

 

주어진 타원이 object을 잘 피팅하려면 타원의 두 축에 대한 회전관성이 각각 object의 주축에 대한  회전관성과 같은 값을 가져야 할 것이다:

$$ I_\text{min} = \mu_{20}^\text{(e)}, \quad I_\text{max} =\mu_{02}^\text{(e)}$$

이 두 식을 풀면 타원의 장축과 단축의 반지름을 구할 수 있다. 

$$a = \Big( \frac{4}{\pi}\Big)^{1/4} \Big( \frac{I_\text{max}^3}{I_\text{min}}\Big)^{1/8}, \quad b = \Big( \frac{4}{\pi}\Big)^{1/4} \Big( \frac{I_\text{min}^3}{I_\text{max}}\Big)^{1/8}$$

 

보통 object을 타원 피팅할 때 윤곽선 정보를 이용하는데, 이 방법은 윤곽선을 추출할 필요가 없어서 편리하다. 단 내부에 빈 곳이 있는 object의 경우 회전관성을 감소시키므로 좋은 결과를 기대할 수 없다. 영상이 다수의 object를 담고 있을 때는 connected component labeling을 한 후 각각의 component에 대해서 fitting을 수행하면 된다.

void getEllipse(CRaster& raster) {
    const double four_pi = 1.0 / atan(1.0);
    CSize sz = raster.GetSize();
    double xsum = 0, ysum = 0;
    double x2sum = 0, y2sum = 0, xysum = 0;
    int count = 0;
    for (int y = 0; y < sz.cy; y++) {
        BYTE *p = (BYTE *)raster.GetLinePtr(y);
        for (int x = 0; x < sz.cx; x++) {
            if (*p++) {
                count++;
                xsum += x;      ysum += y;
                x2sum += x * x; y2sum += y * y; 
                xysum += x * y;
            }
        }
    }
    if (!count) return;
    double xm = xsum / count;
    double ym = ysum / count;
    double m20 = x2sum - xm * xm * count;
    double m02 = y2sum - ym * ym * count;
    double m11 = xysum - xm * ym * count;
    double theta2 = atan2(2 * m11, m20 - m02);
    double ct = cos(theta2), st = sin(theta2);
    double Imin = 0.5 * (m20 + m02) - 0.5 * (m20 - m02) * ct - m11 * st;
    double Imax = 0.5 * (m20 + m02) + 0.5 * (m20 - m02) * ct + m11 * st;
    double major = pow(four_pi, 0.25) * pow(Imax * Imax * Imax / Imin, 0.125);
    double minor = pow(four_pi, 0.25) * pow(Imin * Imin * Imin / Imax, 0.125);
    drawEllipse(raster, mx, my, major, minor, theta2 / 2);
}
 
 
728x90

'Image Recognition > Fundamental' 카테고리의 다른 글

Circle Fitting: Pratt  (0) 2022.01.20
Best-fit Ellipse 2  (0) 2022.01.18
Image Moment  (0) 2021.12.04
Orientation 추정  (0) 2021.11.30
Poisson Image Editing  (0) 2021.08.03
Posted by helloktk
,

등시진자

Physics/역학 2022. 1. 14. 19:25

단순진자는 원호 위에서 반복운동을 한다. 수직에 대해 벌어진 각이 $\theta$일 때 각에 대한 운동 방정식은

$$ \ddot \theta = - \frac{g}{L} \sin \theta.$$

진폭이 작은 경우 ($|\theta| \ll 1$) 윗 식은 용수철 진자의 운동인 단순 조화 운동이 되고 주기는 진폭에 무관하게 일정한 값을 가지게 된다:

$$ \ddot \theta \approx - \frac{g}{L} \theta    \quad \Leftrightarrow  \quad \ddot{x} = - \omega^2 x.$$

그러나 단순진자의 진폭이 일정 이상 커져 작은 각 근사에서 벗어나면 주기는 진폭에 따라 달라짐이 잘 알려져 있다.

 

진자가 원호가 아닌 다른 곡선 위를 움직일 때 주기가 진폭에 무관하게 주어질 수 있는지 알아보자. 이 경우는 각보다는 용수철 진자처럼 평형점에서 움직인 거리를 이용해서 운동을 기술하는 것이 더 편리하다. 단순진자의 경우 평형점에서 잰 원호의 거리를 $s$라면 $s=L \theta$로 표현되고 작은 각 근사에서 운동 방정식은 
$$ \ddot {s} =-\omega^2  s.$$

이제 진자가 움직이는 곡선에 어떤 제약이 들어오는지 살펴보자. 먼저 움직인 거리에 대한 운동 방정식은 단순조화운동식과 같아야 하므로 위의 형태는 변하지 않아야 한다. 진자가 움직이는 곡선이 $y(x)$로 표현되면 평형점($x=0$) 에서 움직인 거리($x <0$이면 움직인 거리의 음수)는

$$s = \int_0^x \sqrt{ 1+ \Big(\frac{dy}{dx}\Big)^2} dx$$

로 표현된다. $s$에 대한 단순 조화 진동은 위치 에너지가 $U(s) = \frac{1}{2} m\omega^2 s^2$인 경우에 해당하는데 진자는 중력의 영향을 받으므로 이 위치 에너지는 결국 중력 위치 에너지 표현되어야 한다: 

$$ \frac{1}{2} m\omega^2 s^2 = mgy.$$

양변을 $x$에 대해 미분한 후 정리하면

$$ \frac{dy}{dx} = \sqrt{ \frac{2\omega^2 y}{  g - 2\omega^2 y}}.$$

이 방정식을 풀기 위해서 새로운 매개변수 $\psi$를 도입하는데, $dy/dx$가 접선의 기울기이므로 $dy/dx=\tan (\psi/2)$로 놓으면(각은 나중을 위해서 2배 한 것임), 

$$   \sin^2(\psi/2) = \frac{2\omega^2 }{g} y  \quad \Longrightarrow \quad y = \frac{g}{4\omega^2} (1 - \cos\psi).$$

$dx/d\psi  = (dy/d\psi) / (dy/dx)$을 계산해서 적분하면 ($x(\psi=0)=c$)

$$  \frac{dx}{d\psi} = \frac{g}{4\omega^2} (1+\cos \psi)\quad \Longrightarrow \quad x= \frac{g}{4\omega^2} (c+\psi + \sin \psi).$$

이제 $\psi \rightarrow \pi +\psi, c=-\pi $로 변수 치환을 하면(이 경우 $dy/dx= -\cot(\psi/2)$) 진자가 움직여야 하는 곡선이 우리가 잘 알고 있는 cycloid임을 알 수 있다.

$$\begin{matrix} x= a ( \psi -\sin \psi )  \\ y=a(1 + \cos \psi)\end{matrix} ~~~~\left(a \equiv \frac{g}{ 4\omega^2}\right).$$

이 곡선은 반지름 $a$인 원이 $y=2a$인 수평선을 따라 구를 때, 원이 처음 $y=2a$와 접하는 점이 그리는 자취를 나타내고, $\psi$는 원의 중심과 그 점이 잇는 선분이 수직과 이루는 각이다.

Cycloid 모양을 결정하는 $a$가 정해지면 진자의 각진동수 $\omega=\sqrt{g/4a}$를 알 수 있고 주기는 $$T=\frac{2\pi}{\omega }= 4\pi \sqrt{\frac{a}{g}}$$

로 주어진다. 곡선이 주어졌으므로 처음 $y=h$에서 출발할 때 구체적으로 주기를 확인해 보자. 역학적 에너지가 보존되므로 

\begin{align} E &=\frac{1}{2} m \Big( \frac{ds}{dt}\Big)^2 +\frac{1}{2} m \omega^2 s^2 \\&= \frac{1}{2} \Big( \frac{ds}{dt}\Big)^2 + mgy = mgh\\ \Longrightarrow~~ dt &= \pm \frac{ds}{\sqrt{2g(h-y)}} \end{align} 주기는 $y=h$에서 출발해서 바닥에 도달하는데 걸리는 시간의 4배이므로 \begin{align} T &= 4 \int_0^h \frac{ds}{\sqrt{2g(h-y)}} \\ &= \frac{4}{\sqrt{2g}} \int_{\pi/2}^{\psi_0} \frac{-4a d \cos (\psi/2)}{\sqrt{2a[\cos^2(\psi_0/2) - \cos^2(\psi/2)]}} \\ &=4\pi \sqrt{\frac{a}{g}}\end{align} 즉, 주기는 출발 높이에 무관하게 주기가 일정함을 알 수 있다.

 

공을 cycloid 모양으로 생긴 골짜기에 굴리면 등시운동을 하지만, 그럼 등시진자는 어떻게 만들수 있을까? 이 문제도 역시 cycloid로 해결이 된다. 위에서 구한 cycloid를 y방향으로 $-2a$ 만큼 평행이동시킨 모양을 고려하자. $$ \begin{matrix} x = a (\psi - \sin \psi) \\ y = a( \cos \psi - 1)\end{matrix}$$이 식으로 표현된 cycloid 모양의 천정을 만든 후(그림의 실선), 원점(꼭대기)에 길이 $L$인 줄을 고정시키고 끝에는 무거운 추를 매단다. 추을 진동을 시키면 줄의 일부는 cycloid 모양의 천정을 따라 접하고 나머지 부분은 직선의 형태로 된다.

줄과 cycloid가 접하는 끝지점을 $(x, y)$라 할 때 접하는 부분의 줄의 길이는 $$ \ell = \int_0^\psi ds = a \int_0^\psi \sqrt{ (1-\cos \psi')^2 + (\sin \psi')^2} d\psi' = 4a [1- \cos (\psi/2)]$$로 주어진다. $(x, y)$ 이후의 줄은 접선의 방향으로 나간다. 접선의 기울기를 $\frac{dy}{dx}=\tan \phi$로 놓으면 $\phi=\psi/2 - \pi/2$이고, 직선 부분의 길이가 $L- \ell$이므로 추의 위치는 $$ \begin{matrix} X = x + (L- \ell)\cos \phi = (L - 4a) \cos \phi + a(\psi + \sin \psi) \\ Y = y+(L-\ell) \sin \phi =(L-4a)\sin \phi - a(3+ \cos \psi)  \end{matrix}$$로 주어진다. 줄의 길이를 $L=4a$로 선택하면 추의 위치 $(X, Y)$도 (평행이동된) cycloid(그림의 점선) 상에서 움직임을 알 수 있다. 따라서 이렇게 만들어진 추의 주기는 등시성을 갖는다. $$ \text{추의 위치:}~~\left\{ \begin{array}{l} X= a (\psi + \sin \psi)  \\ Y = -a( 3 + \cos \psi)\end{array}\right. $$

다시 이 진자의 주기를 구체적으로 확인하면, 

$$ds=  \sqrt{dX^2 + dY^2} = \pm a\sqrt{2(1+\cos \psi)} {d\psi} $$

이고, 역학적 에너지가 보존되므로($\Delta y = a (\cos\psi - \cos \psi_0)$)  

\begin{align} \tau &= 4 \int_{0}^{\psi_0} \frac{ds}{\sqrt{2g \Delta y} }  \\  &= 4\sqrt{\frac{a}{g}} \int_{0}^{\psi_0} \sqrt{\frac{1+\cos \psi}{\cos \psi - \cos\psi_0}} d\psi  \\ &= 4\pi \sqrt{ \frac{a}{g}} = 2\pi \sqrt{ \frac{L}{g}}\end{align}이다. 즉, 줄의 길이가 $L=4a$인 단진자가 작은 진동을 할 때의 주기와 동일하고, cycloid에서 미끄러지면서 운동하는 물체의 진자와도 같음을 알 수 있다.

728x90

'Physics > 역학' 카테고리의 다른 글

개는 얼마나 가속할 수 있을까?  (0) 2022.01.21
공은 얼마나 높이 올라갈까?  (0) 2022.01.19
체인이 떨어지는 가속도는?  (0) 2022.01.12
더 멀리 날아갈까?  (0) 2022.01.12
먼저 바닥에 닿는 다트는?  (0) 2022.01.12
Posted by helloktk
,

Catenary

Mathematics 2022. 1. 14. 10:18

체인이나 줄을 느슨한 상태로 양끝을 고정시킬 때 모양은 포물선처럼 보이지만 실제로는 그렇지 않고 현수선(catenary)라고 불리는 곡선이다. 양끝을 고정시킨 줄을 보자.

늘어진 줄에는 자신의 무게를 지탱하기 위해서 장력이 걸린다. 그런데 중력이 수직방향으로 걸리므로 장력은 줄의 위치에 따라 달라져야 한다.  수평방향은 움직임이 없으므로 장력의 수평 성분은 모두 같아야 하는데, 줄의 가장 아래로 처진 부분의 접선방향이 수평이므로 이 지점에서 장력($T_0$)와 같아야 한다. 줄의 선밀도가 $\lambda$이고, 가장 아래 지점을 기준으로 곡선의 길이를 $s$라고 하자. 현수선의 수평 위치를 $x$, 수직 위치를 $y$로 하면 $y$는 $x$의 함수로 생각할 수 있고, 가장 아랫부분($x=0$으로 잡자)에서 잰 줄의 길이는

$$s =\int\sqrt{dx^2 + dy^2} =  \int_{0}^x \sqrt{1 + (dy/dx)^2} dx$$

로 쓸 수 있다.

그러면 $(x, y)$에서 줄의 장력을 $T$, 접선이 이루는 각도를 $\psi$라면, 힘의 평형 조건에서

\begin{gather} T \cos \psi = T_0 \\ T \sin \psi = \lambda s g . \end{gather}

따라서 접선의 기울기는

$$ \frac{dy}{dx} = \tan \psi = \frac{\lambda gs }{T_0} = \frac{s}{a}= \frac{1}{a} \int_0^{x} { \sqrt{ 1+\Big( \frac{dy}{dx}\Big)^2}} dx, \quad a \equiv \frac{T_0}{\lambda g}$$

이 식은 $y(x)$에 대한 미분-적분 방정식 형태이므로 미분 방정식으로 바꾸기 위해서 한번 더 미분을 하면 

$$ \frac{d^2 y}{dx^2}  = \frac{1}{a} \sqrt{1 + \Big(\frac{dy}{dx}\Big)^2 }.$$

한번 적분하면 ($\frac{dy}{dx} (x=0)=0$)

$$ \frac{dy}{dx}  =\sinh (x/a).$$ 

다시 적분하면 

$$ y = a \cosh( x/a) + C$$

을 얻는 데, 원점을 이동해서 현수선의 가장 아랫부분이 ($0,a$)가 되도록 조정하면 $C=0$이 된다. 현수선은 선밀도와 꼭지점에서의 장력 비 $a$로 모양이 결정된다. 이 값은 한 지점에서 꼭짓점까지 수평거리($x$)와 기울기를 측정하면 결정할 수 있다: 특징 1번. 그리고 $|x| \ll a$일 때

$$y=  a\left[ 1 + \frac{1}{2}\Big(\frac{x}{a}\Big)^2 + \frac{1}{24}\Big( \frac{x}{a}\Big)^4+...\right] $$

이므로 꼭짓점 근방에서는 포물선으로 근사가 가능하다.

몇 가지 특징:

1. 접선의 기울기:  $\tan \psi =\frac{dy}{dx}= \sinh(x/a)$. 

2. 곡선의 길이: $s = a \tan\psi = a \sinh( x/a)$.

3. $y^2 = a^2+ s^2 ~~\rightarrow ~~ y = a\cosh(x/a) = a\sec(\psi)$.

4. 장력: $T = \lambda g \sqrt{a^2 + s^2} =\lambda g \cosh (x/a)= \lambda g y$. 

5. 현수선의 길이($s$)를 매개변수로 선택하면,

$$\frac{dx}{ds} = \frac{1}{ds/dx}= \frac{1}{\sqrt{1 + (dy/dx)^2} }= \frac{a}{\sqrt{a^2 + s^2}},$$

$$ \frac{dy}{ds} = \frac{dy/dx}{ds/dx}=\frac{s}{\sqrt{a^2+s^2}},$$

따라서 $(dx/ds)^2 + (dy/ds)^2 = 1$임을 알 수 있다.

6. $\frac{dx}{d s} = \frac{1}{a}   \cos^2(\psi)\frac{dx}{d\psi} = \cos(\psi) ~~\rightarrow~~ \frac{dx}{d\psi} = a \sec{\psi} $  이므로 $x= a \ln ( \sec(\psi) + \tan(\psi))$

Note: 미분방정식의 유도를 local version으로 바꾸자. 줄의 장력의 $x$의 함수로 볼 수 있고, 그림에서 힘의 평형을 적용하면, 우선 수평방향에 대해서

$$ \sum F_x = T(x+dx) \cos[ \psi(x+dx)] - T(x) \cos [\psi(x)]=0$$ 

이 식은 장력의 수평성분은 어디서나 같음을 의미하므로 이 값을 꼭지점($\psi=0$)의 값 $T_0$로 고정하자.

수직방향에 대해서는

$$ \sum F_y = T(x+dx) \sin[ \psi(x+dx)] - T(x) \sin[\psi(x)] -dm g=0$$

$T(x+dx)$와 $T(x)$을 소거하고 정리하면

$$  \tan  [\psi(x+dx)] -\tan[\psi(x)] = \frac{dm g}{T_0} = \frac{\lambda g}{T_0} \sqrt{dx^2+dy^2}.$$ 

그런데 $dy/dx = \tan(\psi)$이므로 

 $$ \tan[\psi(x+dx)]- \tan[\psi(x)] = dx \frac{d}{dx} \tan(\psi) =    \frac{d^2y}{dx^2 }dx$$

정리되어 현수선 방정식을 얻을 수 있다.

    

728x90
Posted by helloktk
,