Bối cảnh và vấn đề nghiên cứu

Trong lịch sử phát triển của cơ học cổ điển kéo dài hơn 400 năm qua các công trình của Galileo Galilei, Robert Hooke, Charles-Augustin de Coulomb, James Clerk Maxwell, John von Neumann và Albert Einstein, cơ học vật rắn luôn giữ vai trò then chốt trong cả nghiên cứu lý thuyết lẫn kỹ thuật ứng dụng như cơ khí chế tạo và kết cấu công trình. Phương pháp tiếp cận truyền thống đối với các bài toán động lực học vật rắn là sử dụng công cụ giải tích thuần túy (đạo hàm, tích phân, giải phương trình vi phân) nhằm thiết lập nghiệm giải tích tường minh mô tả mối liên hệ giữa các đại lượng vật lý. Tuy nhiên, khi đưa vào các điều kiện vật lý thực tế như lực cản môi trường tỉ lệ với bình phương vận tốc, lực đẩy Archimedes hay tác động của gió biến thiên theo thời gian, việc tìm nghiệm giải tích dạng hàm sơ cấp trở nên bất khả thi hoặc chưa có lời giải chính xác trong các giáo trình cơ học chuẩn mực. Khi đó, phương pháp giải số trên máy tính bằng các ngôn ngữ lập trình như Fortran, Python, Matlab hoặc C trở thành công cụ tất yếu để mô phỏng và khảo sát các quy luật chuyển động phức tạp.

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Huệ, thực hiện năm 2022 tại Khoa Vật lý - Trường Đại học Sư phạm Thành phố Hồ Chí Minh dưới sự hướng dẫn của PGS. Phạm Nguyễn Thành Vinh, tập trung giải quyết bài toán mô phỏng chuyển động của các hệ vật rắn khi có sự tham gia của các yếu tố cản thực tế. Khóa luận xác lập các mục tiêu và nhiệm vụ cụ thể sau:

  1. Hệ thống hóa cơ sở lý thuyết động học và động lực học vật rắn quanh một trục cố định, cùng quy trình 4 bước chuẩn để giải một bài toán động lực học.
  2. Trình bày phương pháp giải số bài toán Cauchy đối với phương trình vi phân cấp hai thông qua giải thuật Runge-Kutta bậc 4 (RK4) và hiện thực hóa thuật toán bằng ngôn ngữ lập trình Fortran.
  3. Ứng dụng chương trình giải số để mô phỏng, khảo sát và so sánh hành vi động lực học của bốn hệ dao động: con lắc đơn với các biên độ góc khác nhau; con lắc vật lý thanh đồng chất; con lắc antenna dưới đáy đại dương chịu lực cản chất lỏng bậc hai và lực đẩy Archimedes; con lắc quả cầu đặc chịu lực cản không khí bậc hai cùng tác động của gió ngang (không gió, gió không đổi và gió điều hòa).

Đối tượng nghiên cứu của đề tài là các đại lượng động học và động lực học (tọa độ góc $\theta(t)$, vận tốc góc $\omega(t)$, gia tốc góc $\gamma(t)$, moment quán tính, moment lực) của các cơ hệ vật rắn quay quanh trục cố định. Phạm vi nghiên cứu giới hạn ở các mô hình dao động thuộc chương trình Vật lý đại cương, không xét đến lý thuyết cơ học giải tích Lagrange hay Hamilton, với thời gian khảo sát mô phỏng số từ 10 giây đến 60 giây.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của khóa luận dựa trên cơ học vật rắn cổ điển theo quan điểm động lực học Newton. Động học chuyển động quay được mô tả thông qua góc quay $\theta(t)$, vận tốc góc $\omega(t) = d\theta/dt$ và gia tốc góc $\gamma(t) = d\omega/dt = d^2\theta/dt^2$. Mối liên hệ quán tính được xác định qua moment quán tính $I_O = \sum m_i \rho_i^2$ (đối với hệ chất điểm) hoặc tích phân $I_O = \int \rho^2 dm$ (đối với vật rắn liên tục). Trường hợp trục quay không đi qua khối tâm, tác giả sử dụng định lý Steiner-Huygens (định lý trục song song): $I_O = I_{cm} + md^2$, trong đó $I_{cm}$ là moment quán tính qua khối tâm và $d$ là khoảng cách giữa hai trục song song. Động lực học quay được xác lập bằng phương trình cơ bản $\sum \vec{\tau} = I \vec{\gamma}$ và định luật bảo toàn/biến thiên moment động lượng $d\vec{L}/dt = \vec{\tau}$.

Về mặt toán học và phương pháp tính, khóa luận tiếp cận bài toán chuyển động quay dưới dạng bài toán giá trị ban đầu Cauchy đối với phương trình vi phân cấp hai:

$$y''(t) = f(t, y(t), y'(t))$$

với các điều kiện đầu cho trước $y(t_0) = y_0$ và $y'(t_0) = y'_0$. Tác giả thực hiện phân rã phương trình cấp hai thành hệ hai phương trình vi phân cấp một bằng cách đặt biến phụ: $y_1(t) = y(t)$ và $y_2(t) = y'(t)$, dẫn đến hệ:

$$\begin{cases} y_1'(t) = y_2(t) \ y_2'(t) = f(t, y_1(t), y_2(t)) \end{cases}$$

Miền thời gian $[a, b]$ được rời rạc hóa thành $n$ khoảng bằng nhau với bước thời gian $h = \Delta t = (b - a)/n$. Thuật toán giải số được sử dụng là phương pháp Runge-Kutta bậc 4 (RK4), được phát triển từ khai triển chuỗi Taylor nhằm đạt độ chính xác bậc 4 mà không đòi hỏi tính toán giải tích các đạo hàm cấp cao:

$$y_{i+1} = y_i + \frac{1}{6}(k_1 + 2k_2 + 2k_3 + k_4)h$$

trong đó $k_1, k_2, k_3, k_4$ là các hệ số gia góc xác định tại các điểm mốc trong bước chia $h$. Toàn bộ giải thuật RK4 được lập trình bằng ngôn ngữ Fortran, gọi hàm con xử lý đạo hàm DERIV với hai biến trạng thái temp(1) và temp(2) tương ứng với góc và vận tốc góc.

Nguồn dữ liệu tính toán bao gồm các tham số mô hình hóa được trích xuất từ tài liệu cơ học chuẩn (bài giảng của Richard Feynman tại CalTech năm 1961 về giải số chuyển động hành tinh), dữ liệu khảo sát thực nghiệm về thiết bị antenna địa vật lý dưới đáy đại dương và các thông số khí động học tiêu chuẩn của không khí và quả cầu kim loại.

Nội dung chính theo từng chương

Chương 1: Lý thuyết cơ học vật rắn

Chương 1 hệ thống hóa toàn bộ các kiến thức trọng tâm của cơ học vật rắn phục vụ cho việc thiết lập phương trình chuyển động. Khóa luận trình bày sự tương ứng chặt chẽ giữa chuyển động tịnh tiến và chuyển động quay quanh trục cố định, được tổng hợp trong bảng so sánh động học:

Đại lượng động học Chuyển động tịnh tiến Chuyển động quay
Vị trí / Tọa độ $x(t)$ $\theta(t)$
Vận tốc $v(t) = dx/dt$ $\omega(t) = d\theta/dt$
Gia tốc $a(t) = dv/dt = d^2x/dt^2$ $\gamma(t) = d\omega/dt = d^2\theta/dt^2$
Phương trình vi phân $d^2x/dt^2 = a(t, x, v)$ $d^2\theta/dt^2 = \gamma(t, \theta, \omega)$

Đối với các chuyển động đặc biệt có gia tốc góc không đổi hoặc vận tốc góc không đổi, phương trình chuyển động quay suy biến về các công thức tích phân trực tiếp. Khi vật rắn quay quanh trục bất kỳ, vận tốc dài và gia tốc dài của một điểm tại bán kính $\vec{r}$ được mô tả qua tích vector: $\vec{v} = \vec{\omega} \times \vec{r}$ và $\vec{a} = \vec{\gamma} \times \vec{r} + \vec{\omega} \times (\vec{\omega} \times \vec{r})$.

Về đặc trưng quán tính, tác giả trình bày công thức moment quán tính khối tâm của các vật thể đồng chất có dạng hình học đối xứng:

Dạng hình học vật thể đồng chất Moment quán tính khối tâm ($I_{cm}$)
Thanh thẳng mảnh chiều dài $L$ $\frac{1}{12} M L^2$
Vành dây / Trụ rỗng bán kính $R$ $M R^2$
Đĩa phẳng / Trụ đặc bán kính $R$ $\frac{1}{2} M R^2$
Quả cầu đặc bán kính $R$ $\frac{2}{5} M R^2$
Quả cầu rỗng bán kính $R$ $\frac{2}{3} M R^2$

Về động lực học, chương này phân tích moment của lực $\vec{\tau} = \vec{r} \times \vec{F}$, định luật hai Newton cho chuyển động quay quanh khối tâm và quanh trục cố định ngoài khối tâm, định nghĩa moment động lượng $\vec{L} = \vec{r} \times \vec{p}$ và định lý biến thiên moment động lượng. Cuối chương, tác giả đúc kết quy trình 4 bước chuẩn để giải quyết bài toán động lực học vật rắn:

  1. Nhận dạng: Xác định tọa độ góc $\theta(t)$, chọn hệ quy chiếu, xác định các điều kiện ban đầu $\theta(t_0)$ và $\omega(t_0)$.
  2. Mô hình: Xác định các lực tác dụng, tính moment của từng lực, áp dụng định luật hai Newton cho chuyển động quay để thiết lập phương trình vi phân $\gamma = \gamma(\theta, \omega, t)$.
  3. Giải: Giải phương trình vi phân bằng phương pháp giải tích hoặc giải số.
  4. Phân tích: Kiểm tra tính đúng đắn của nghiệm $\theta(t)$ và $\omega(t)$, đánh giá đặc tính chuyển động của cơ hệ.

Chương 2: Phương pháp tính toán

Chương 2 trình bày cơ sở toán học của việc giải số phương trình vi phân. Tác giả mở đầu bằng bài toán giá trị ban đầu Cauchy, chứng minh cách đưa một phương trình vi phân cấp hai bất kỳ về hệ hai phương trình vi phân cấp một thông qua việc đặt biến phụ $y_1(t) = y(t)$ và $y_2(t) = y'(t)$. Quá trình rời rạc hóa trên lưới thời gian với bước chia $h = \Delta t$ được giải thích làm cơ sở cho tính toán số.

Tiếp đó, tác giả phân tích nguồn gốc giải thuật Runge-Kutta từ chuỗi Taylor. Bằng cách ngắt chuỗi Taylor ở bậc nhất, ta thu được phương pháp Euler:

$$y_{i+1} = y_i + h f(t_i, y_i)$$

Phương pháp Euler có sai số cục bộ lớn nên được mở rộng lên phương pháp Runge-Kutta bậc 2 (RK2). Dựa trên khai triển Taylor hàm hai biến với phần dư Peano và đồng nhất các hệ số, khóa luận dẫn ra hệ phương trình điều kiện: $a_1 + a_2 = 1$, $a_2 p_1 = 1/2$, $a_2 q_{11} = 1/2$. Chọn $a_2 = 1/2$ cho phép xác định nghiệm $a_1 = 1/2$, $p_1 = q_{11} = 1$, từ đó thiết lập công thức RK2:

$$y_{i+1} = y_i + \frac{h}{2}(k_1 + k_2)$$

với $k_1 = f(t_i, y_i)$ và $k_2 = f(t_i + h, y_i + h k_1)$.

Mở rộng quy trình đồng nhất thức chuỗi Taylor lên bậc 4, tác giả đưa ra công thức giải thuật Runge-Kutta bậc 4 (RK4) hoàn chỉnh với 4 hệ số gia $k_1, k_2, k_3, k_4$. Ở phần hiện thực hóa thuật toán, chương trình tính toán được viết bằng ngôn ngữ Fortran sử dụng cấu trúc chương trình con rk4. Đối với phương trình cấp hai ($N = 2$), chương trình khai báo hàm DERIV chứa hai thành phần: hàm thứ nhất gán đạo hàm bậc nhất qua biến temp(2), hàm thứ hai mô tả phương trình vi phân gia tốc góc phụ thuộc vào tọa độ temp(1) và vận tốc temp(2). Đây là hạt nhân tính toán được áp dụng xuyên suốt cho các bài toán ở Chương 3.

Chương 3: Kết quả nghiên cứu

Chương 3 trình bày các kết quả mô phỏng và khảo sát số đối với 4 bài toán dao động cơ học vật rắn:

1. Bài toán con lắc đơn: Phương trình chuyển động quay chính xác của chất điểm khối lượng $m$, chiều dài dây $l$ là:

$$\frac{d^2\theta}{dt^2} = -\frac{g}{l}\sin\theta$$

Khi góc nhỏ ($\sin\theta \approx \theta$), phương trình có nghiệm giải tích dao động điều hòa: $\theta(t) = \theta_0 \cos(\omega t)$ với $\omega = \sqrt{g/l}$. Chọn $g/l = \pi^2$, nghiệm giải tích trở thành $\theta(t) = \theta_0 \cos(\pi t)$. Tác giả giải số trực tiếp phương trình phi tuyến bằng RK4 với khai báo DERIV = -3.14159265359**2 * sin(temp(1)) và bước thời gian tstep = 0.01 s. Khảo sát tại 4 biên độ góc ban đầu $\theta_0 = \pi/20, \pi/10, \pi/6, \pi/3$ cho thấy:

  • Tại $\theta_0 = \pi/20$ và $\pi/10$ (tương ứng $\le 18^\circ$), nghiệm giải tích trùng khớp với nghiệm giải số RK4.
  • Tại $\theta_0 = \pi/6$ và $\pi/3$, nghiệm giải tích gần đúng góc nhỏ lệch dần so với nghiệm số chính xác, độ lệch tăng dần theo thời gian khảo sát, minh chứng sự phá vỡ tính điều hòa dù chuyển động vẫn bảo toàn tính tuần hoàn.

2. Bài toán con lắc vật lý thanh đồng chất: Xét một thanh đồng chất khối lượng $m = 1\text{ kg}$, chiều dài $d = 0.25\text{ m}$, trục quay cố định tại một đầu thanh. Khối tâm nằm tại trung điểm ($l_{cm} = d/2 = 0.125\text{ m}$). Theo định lý Steiner-Huygens, moment quán tính quanh điểm treo đầu thanh là:

$$I_S = I_{cm} + m(d/2)^2 = \frac{1}{12}md^2 + \frac{1}{4}md^2 = \frac{1}{3}md^2 = \frac{1}{48}\text{ kg}\cdot\text{m}^2$$

Với gia tốc trọng trường $g = \pi^2$, tần số góc dao động nhỏ là $\omega_0 = \sqrt{m g l_{cm} / I_S} = \sqrt{6}\pi\text{ rad/s}$. Phương trình vi phân chính xác là $d^2\theta/dt^2 = -6\pi^2 \sin\theta \approx -59.2176 \sin\theta$. Chương trình RK4 mô phỏng với tstep = 0.01 s trong 10 giây đầu tiên tại các góc thả $\theta_0 = \pi/10$ và $\pi/6$. Đồ thị số xác nhận tính điều hòa của thanh bị suy giảm nhanh hơn theo thời gian khi biên độ góc tăng lên $\pi/6$.

3. Bài toán con lắc antenna dưới đáy đại dương: Mô phỏng thiết bị antenna dài $L = 5\text{ m}$ treo cố định dưới đáy biển dùng để đo tín hiệu điện từ dò tìm hydrocarbon. Thiết bị chịu tác động của trọng lực, lực đẩy Archimedes và lực cản chất lỏng tỉ lệ với bình phương vận tốc. Phương trình gia tốc góc có dạng:

$$\frac{d^2\theta}{dt^2} = -c_1 \sin\theta - c_2 \left|\frac{d\theta}{dt}\right|\frac{d\theta}{dt}$$

với các hệ số thực nghiệm $c_1$ (thể hiện gia tốc biểu kiến do hợp lực trọng lực và lực đẩy Archimedes) và $c_2$ (đặc trưng cho lực cản chất lỏng). Điều kiện ban đầu: $\theta_0 = 10^\circ = \pi/18\text{ rad}$, $\omega_0 = 0$. Kết quả giải số trong 60 giây mô tả chi tiết sự suy giảm biên độ dao động của tọa độ góc $\theta(t)$, vận tốc góc $\omega(t)$ và gia tốc góc $\gamma(t)$, khẳng định vai trò chi phối của số hạng phi tuyến $-c_2 |\omega|\omega$ trong quá trình tắt dần dưới nước.

4. Bài toán con lắc vật lý chịu lực cản không khí và tác dụng của gió: Hệ gồm quả cầu đặc khối lượng $M$, bán kính $R = 0.02\text{ m}$, khối lượng riêng $\rho_m = 7000\text{ kg/m}^3$, treo bằng sợi dây nhẹ dài $L = 1\text{ m}$. Lấy $g = 9.8\text{ m/s}^2$, mật độ không khí $\rho_a = 1.293\text{ kg/m}^3$, hệ số cản $D = \frac{1}{2}\rho_a \pi R^2$. Moment quán tính của hệ là $I = \frac{2}{5}MR^2 + ML^2$. Vận tốc tương đối giữa quả cầu và không khí là $\vec{v} - \vec{u}$ với gió thổi ngang $\vec{u} = u(t)\vec{i}$. Tác giả khảo sát 3 trường hợp:

  • Không có gió ($u = 0$): Với kích thích ban đầu nhỏ $\omega_{max} = 2\text{ rad/s}$ và $6\text{ rad/s}$, con lắc dao động tắt dần quanh vị trí cân bằng $\theta = 0$. Với kích thích mạnh $\omega_{max} = 10\text{ rad/s}$, cơ năng ban đầu đủ lớn khiến con lắc quay tròn 4 vòng liên tiếp trong 5 giây đầu (tọa độ góc vượt qua $\theta = 8\pi \approx 25.13\text{ rad}$), sau đó mới chuyển sang dao động tắt dần quanh vị trí cân bằng mới $\theta = 8\pi$.
  • Gió thổi ngang đều ($u_0 = 20\text{ m/s}$): Sự tắt dần diễn ra nhanh hơn rõ rệt so với khi không có gió do lực cản tương đối tăng mạnh. Vị trí cân bằng bị lệch khỏi các góc $2n\pi$, tương tự như con lắc đặt trong điện trường nằm ngang.
  • Gió thổi điều hòa ($u(t) = \frac{u_0}{2}[1 + \cos(2\pi t / T)]$): Khảo sát với $\omega_{max} = 8\text{ rad/s}$ tại hai chu kỳ gió $T = 5\text{ s}$ và $T = 10\text{ s}$. Kết quả cho thấy chu kỳ dao động của con lắc bị biến điệu rõ rệt; dao động không tắt dần hoàn toàn do được gió cấp năng lượng tuần hoàn. Khi $T = 10\text{ s}$, tại thời điểm $t = 15\text{ s}$ biên độ giảm gần về 0, nhưng đến $t = 20\text{ s}$ biên độ lại tăng vọt trở lại do gió đạt cực đại.

Kết quả và đóng góp

Khóa luận đã tổng hợp và lượng hóa thành công các tham số động lực học trong các bài toán khảo sát số:

Cơ hệ khảo sát Tham số đặc trưng Điều kiện ban đầu Phát hiện động lực học chính từ giải số
Con lắc đơn $g/l = \pi^2$, $h = 0.01\text{ s}$ $\omega_0 = 0$, $\theta_0 = \pi/20 \to \pi/3$ Nghiệm giải tích góc nhỏ khớp nghiệm số ở $\theta_0 \le \pi/10$; phân kỳ mạnh ở $\theta_0 \ge \pi/6$
Con lắc thanh thẳng $m = 1\text{ kg}$, $d = 0.25\text{ m}$, $I_S = 1/48\text{ kg}\cdot\text{m}^2$ $\omega_0 = 0$, $\theta_0 = \pi/10, \pi/6$ Khảo sát trong 10 s; tính điều hòa suy giảm nhanh ở $\theta_0 = \pi/6$ so với $\pi/10$
Antenna dưới biển $L = 5\text{ m}$, lực cản bậc hai $-c_2|\omega|\omega$ $\theta_0 = \pi/18$, $\omega_0 = 0$ Dao động tắt dần phi tuyến rõ rệt trên cả $\theta(t)$, $\omega(t)$, $\gamma(t)$ trong 60 s
Con lắc quả cầu (không gió) $R = 0.02\text{ m}$, $L = 1\text{ m}$, $\rho_m = 7000\text{ kg/m}^3$ $\theta_0 = 0$, $\omega_{max} = 2, 6, 10\text{ rad/s}$ $\omega_{max} = 10\text{ rad/s}$ làm quả cầu quay 4 vòng ($8\pi$) trong 5 s đầu trước khi tắt dần
Con lắc quả cầu (gió đều) $u_0 = 20\text{ m/s}$ $\theta_0 = 0$, $\omega_{max} = 2, 4, 6, 8, 10\text{ rad/s}$ Tắt dần nhanh hơn không gió; vị trí cân bằng dịch chuyển lệch khỏi $2n\pi$
Con lắc quả cầu (gió điều hòa) $u_0 = 20\text{ m/s}$, $T = 5\text{ s}, 10\text{ s}$ $\theta_0 = 0$, $\omega_{max} = 8\text{ rad/s}$ Dao động không tắt hẳn; biên độ phục hồi chu kỳ (tăng trở lại lúc $t = 20\text{ s}$ khi $T = 10\text{ s}$)

Đóng góp của khóa luận được tác giả tự xác định bao gồm:

  • Chứng minh tính hiệu quả và ưu việt của phương pháp tính toán giải số bằng thuật toán RK4 trên ngôn ngữ Fortran trong việc giải quyết các phương trình động lực học phi tuyến phức tạp mà phương pháp giải tích không thể xử lý.
  • Xây dựng được mô hình mô phỏng chuyển động trực quan cho các bài toán dao động có lực cản phi tuyến (lực cản tỉ lệ với bình phương vận tốc) và chịu ngoại lực phụ thuộc thời gian (gió thổi điều hòa), cung cấp bức tranh vật lý định lượng chuẩn xác về các đại lượng tọa độ, vận tốc và gia tốc góc.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận tự chỉ ra một số giới hạn và định hướng mở rộng:

  • Hạn chế: Nghiên cứu mới dừng lại ở các bài toán cơ học vật rắn phẳng, chuyển động quay một chiều quanh một trục cố định trong phạm vi chương trình Vật lý đại cương. Chưa mở rộng sang các phương trình vi phân bậc cao hơn hai hoặc các hệ nhiều bậc tự do sử dụng cơ học giải tích Lagrange/Hamilton.
  • Hướng nghiên cứu tiếp: Ứng dụng chương trình tính toán giải số vào đổi mới phương pháp giảng dạy học phần Vật lý đại cương nhằm trực quan hóa các hiện tượng cơ học phi tuyến. Mở rộng bài toán mô phỏng cho các hệ cơ học phức tạp hơn như con lắc chịu tác dụng của gió thổi theo phương bất kỳ trong không gian 3 chiều hoặc gió đổi hướng liên tục. Đồng thời, nâng cấp thuật toán để xử lý các hệ phương trình vi phân bậc cao với nhiều biến trạng thái tương tác.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo hữu ích cho:

  • Sinh viên ngành Sư phạm Vật lý và Vật lý học: Cung cấp tài liệu mẫu về việc kết hợp giữa lý thuyết cơ học cổ điển và phương pháp tính toán số, phục vụ học tập các học phần Cơ học đại cương và Vật lý tính toán.
  • Giảng viên và người làm công tác giảng dạy: Sử dụng các mô hình và đồ thị mô phỏng trong khóa luận làm học liệu trực quan giải thích sự khác biệt giữa gần đúng giải tích góc nhỏ và nghiệm thực tế, cũng như bản chất của dao động tắt dần phi tuyến và dao động có kích thích ngoại vi.
  • Phần đáng tham khảo nhất: Chương 2 (cách biến đổi hệ phương trình Cauchy và thuật toán RK4) cùng Chương 3 (kỹ thuật xây dựng hàm đạo hàm DERIV mô hình hóa lực cản phụ thuộc bình phương vận tốc tương đối $|\vec{v} - \vec{u}|(\vec{v} - \vec{u})$).

Câu hỏi thường gặp

1. Tại sao tác giả phải dùng phương pháp giải số RK4 mà không giải bằng phương pháp giải tích truyền thống?
Phương pháp giải tích chỉ giải được khi phương trình vi phân là tuyến tính hoặc sử dụng gần đúng góc nhỏ ($\sin\theta \approx \theta$). Khi đưa các yếu tố thực tế vào bài toán như lực cản môi trường tỉ lệ với bình phương vận tốc (số hạng $-c_2 |\omega|\omega$), lực cản tương đối của gió hay góc lệch ban đầu lớn ($\theta_0 = \pi/6, \pi/3$), phương trình vi phân trở nên phi tuyến mạnh và không tồn tại nghiệm giải tích dạng hàm sơ cấp trong các giáo trình cơ học hiện hành. Do đó, phương pháp giải số RK4 là bắt buộc để tìm nghiệm số chính xác.

2. Điểm khác biệt căn bản giữa dao động con lắc đơn ở góc nhỏ và góc lớn trong kết quả mô phỏng là gì?
Tại các góc nhỏ $\theta_0 = \pi/20$ và $\pi/10$ ($\le 18^\circ$), nghiệm số RK4 hoàn toàn trùng khớp với nghiệm giải tích dao động điều hòa $\theta(t) = \theta_0 \cos(\pi t)$. Ngược lại, tại các góc lớn $\theta_0 = \pi/6$ và $\pi/3$, nghiệm giải tích điều hòa bị lệch dần và phân kỳ so với nghiệm số theo thời gian, chứng tỏ hệ mất tính điều hòa nhưng vẫn duy trì tính chất tuần hoàn.

3. Hiện tượng đặc biệt nào xảy ra khi kích thích quả cầu với vận tốc góc ban đầu $\omega_{max} = 10\text{ rad/s}$ trong trường hợp không có gió?
Với $\omega_{max} = 10\text{ rad/s}$, cơ năng ban đầu của hệ rất lớn, khiến quả cầu không dao động ngay mà quay tròn 4 vòng liên tiếp trong 5 giây đầu tiên (tọa độ góc tăng vượt ngưỡng $\theta = 8\pi \approx 25.13\text{ rad}$). Sau khi lực cản không khí làm tiêu tán bớt năng lượng sau 5 giây đầu, con lắc mới bắt đầu chuyển sang chế độ dao động tắt dần quanh vị trí cân bằng mới là $\theta = 8\pi$.

4. Khi có gió thổi ngang điều hòa, tại sao dao động của con lắc không tắt dần hoàn toàn?
Gió thổi điều hòa theo quy luật $u(t) = \frac{u_0}{2}[1 + \cos(2\pi t / T)]$ cung cấp năng lượng tuần hoàn bù đắp cho phần năng lượng bị tiêu tán bởi lực cản không khí. Điển hình với chu kỳ $T = 10\text{ s}$, tại thời điểm $t = 15\text{ s}$ biên độ dao động đã giảm tiến sát về 0, nhưng đến thời điểm $t = 20\text{ s}$ biên độ lại tăng vọt trở lại tương ứng với lúc vận tốc gió thổi vào con lắc đạt cực đại.

5. Cấu trúc hàm DERIV trong chương trình Fortran của khóa luận nhận nhiệm vụ gì?
Hàm DERIV đóng vai trò là đầu vào quan trọng nhất của chương trình giải số, có nhiệm vụ biểu diễn vế phải của hệ hai phương trình vi phân cấp một. Thành phần thứ nhất xác định đạo hàm của tọa độ góc ($d\theta/dt$) tương ứng với biến vận tốc temp(2), thành phần thứ hai xác định gia tốc góc ($d\omega/dt$) biểu diễn tường minh theo góc temp(1), vận tốc temp(2) và biến thời gian $t$.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Huệ đã hệ thống hóa toàn diện cơ sở lý thuyết động lực học vật rắn và xây dựng thành công bộ chương trình tính toán giải số bằng giải thuật Runge-Kutta bậc 4 trên ngôn ngữ Fortran. Thông qua việc giải quyết bài toán con lắc đơn, con lắc vật lý thanh đồng chất, antenna dưới nước và con lắc quả cầu chịu lực cản của gió, công trình đã làm sáng tỏ các quy luật động lực học phi tuyến phức tạp vốn không thể tiếp cận bằng giải tích thuần túy. Kết quả nghiên cứu không chỉ khẳng định tính ưu việt của phương pháp giải số trong vật lý mà còn cung cấp nguồn học liệu mô phỏng có giá trị ứng dụng thiết thực cho công tác giảng dạy Vật lý đại cương.