문제집에 예뻐보이는 그림을 하나 골라봤습니다.
언제나 그렇듯이 m=k=1로 합니다.
g=10, M=2, l=4입니다.
이번에도 처음에 불길한 느낌은 있었다만, 상황의 선택에 거창하게 실패했습니다.
하지만 발빼지 않고 어떻게든 식 3개를 연립해서 각 변수 (x, r, theta)의 이계도함수에 대한 식을 구해보았습니다.
그런데 마지막 식을 다 세우고 코드 창을 열어보는 순간 떠올렸습니다.
아... 퍼텐셜 에너지에 탄성에너지를 빼먹었구나...
그래서 의욕이 훅 떨어져 이번 문제는 그만 풀도록 하겠습니다.
용수철 진자는 다음에 한번 더 도전하도록 하지요.
(용수철 진자는 2차원도 꽤 까다롭긴 한데, 차원이 더 올라갈수록 더욱 어려워지는 문제입니다)
한편 위의 3계연립미분방정식을 풀이하면 탄성력이 작용하지 않는 '평면 위에 미끄러지는 물체에 달린 진자'가 됩니다.
'[硏 Ⅰ] 연구하다 - 교양 > 수치해석' 카테고리의 다른 글
벽에 용수철로 연결된 그릇 속, 두개의 용수철로 이어진 물체 (0) | 2021.07.09 |
---|---|
스프링으로 연결되어 서로 평행한 직선 위에서 움직이는 일련의 질점 (5) | 2021.07.07 |
[미완] 천장에 실 두개로 매달린 평행 선형 진자에 매달린 두 개의 진자 (1) | 2021.07.06 |
두개의 수평 용수철에 매달린 질점의 운동 (0) | 2021.07.06 |
python - n체 문제(n-body problem) (0) | 2021.06.23 |
댓글