Giải một số bài toán cơ học vật rắn một cách hiệu quả

Cơ học vật rắn là một lĩnh vực quan trọng trong vật lý, nghiên cứu về chuyển động và biến… 55 trang · Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh · 2022.

• Tác giả: Nguyễn Thị Huệ
• Người hướng dẫn: PGS. Phạm Nguyễn Thành Vinh
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Sư Phạm Vật Lý
• Đề tài: Giải Số Một Số Bài Toán Cơ Học Vật Rắn
• Loại tài liệu: khóa luận tốt nghiệp
• Năm xuất bản: 2022
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• MỤC LỤC
• DANH SÁCH BẢNG
• DANH SÁCH HÌNH VẼ
• KÍ HIỆU VIẾT TẮT
• GIỚI THIỆU
• 1. CHƯƠNG 1: LÝ THUYẾT CƠ HỌC VẬT RẮN
• 1.1. Chuyển động quay
• 1.2. Động lực học vật rắn
• 2. CHƯƠNG 2: PHƯƠNG PHÁP TÍNH TOÁN
• 2.1. Bài toán Cauchy
• 2.2. Phương pháp Runge-Kutta
• 2.2.1. Phương pháp Runge-Kutta bậc 2
• 2.2.2. Phương pháp Runge-Kutta bậc 4
• 2.3. Chương trình giải số
• 3. CHƯƠNG 3: KẾT QUẢ NGHIÊN CỨU
• 3.1. Con lắc đơn
• 3.2. Con lắc vật lý
• 3.3. Con lắc vật lý đặt trong môi trường chất lỏng
• 3.4. Con lắc vật lý dưới tác dụng của gió
• KẾT LUẬN VÀ ĐỊNH HƯỚNG
• PHỤ LỤC A: CHƯƠNG TRÌNH GIẢI SỐ RK4
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Giới thiệu về giải quyết bài toán cơ học vật rắn hiệu quả

Cơ học vật rắn là một lĩnh vực quan trọng trong vật lý, nghiên cứu về chuyển động và biến dạng của vật rắn dưới tác động của lực. Việc giải quyết bài toán cơ học vật rắn không chỉ có ý nghĩa lý thuyết mà còn có ứng dụng thực tiễn trong nhiều lĩnh vực như cơ khí, xây dựng và kỹ thuật. Khóa luận này sẽ trình bày các phương pháp giải số hiệu quả nhằm phân tích các bài toán cơ học vật rắn.

1.1. Tổng quan về cơ học vật rắn

Cơ học vật rắn nghiên cứu các hiện tượng vật lý liên quan đến chuyển động và lực tác động lên vật rắn. Các khái niệm cơ bản như chuyển động tịnh tiến và quay sẽ được trình bày.

1.2. Tầm quan trọng của việc giải bài toán cơ học

Giải bài toán cơ học vật rắn giúp hiểu rõ hơn về các hiện tượng vật lý trong thực tế, từ đó ứng dụng vào thiết kế và cải tiến các sản phẩm kỹ thuật.

II. Những thách thức trong giải quyết bài toán cơ học vật rắn

Mặc dù có nhiều phương pháp giải, nhưng việc giải quyết bài toán cơ học vật rắn vẫn gặp nhiều thách thức. Các yếu tố như điều kiện biên phức tạp, lực cản và các yếu tố môi trường có thể làm cho việc giải bài toán trở nên khó khăn.

2.1. Vấn đề với điều kiện biên phức tạp

Điều kiện biên phức tạp có thể làm cho các phương trình vi phân trở nên khó giải. Việc xác định các điều kiện đầu và biên là rất quan trọng trong quá trình giải.

2.2. Ảnh hưởng của lực cản và môi trường

Lực cản từ môi trường như không khí hoặc nước có thể ảnh hưởng đến chuyển động của vật rắn, làm cho việc tính toán trở nên phức tạp hơn.

III. Phương pháp giải bài toán cơ học vật rắn hiệu quả

Để giải quyết các bài toán cơ học vật rắn, nhiều phương pháp đã được phát triển. Trong đó, phương pháp giải số là một trong những phương pháp hiệu quả nhất, cho phép tính toán các đại lượng vật lý một cách chính xác.

3.1. Phương pháp Runge Kutta

Phương pháp Runge-Kutta là một trong những phương pháp giải số phổ biến, cho phép giải các phương trình vi phân một cách chính xác và hiệu quả.

3.2. Ứng dụng của phương pháp giải số

Phương pháp giải số có thể được áp dụng để mô phỏng chuyển động của các vật rắn trong nhiều tình huống khác nhau, từ đó giúp phân tích và dự đoán hành vi của chúng.

IV. Kết quả nghiên cứu từ việc giải bài toán cơ học vật rắn

Kết quả nghiên cứu cho thấy rằng việc áp dụng phương pháp giải số đã mang lại những hiểu biết sâu sắc về chuyển động của vật rắn. Các mô hình được xây dựng đã cho thấy sự chính xác trong việc dự đoán hành vi của các hệ thống vật lý.

4.1. Kết quả từ mô hình con lắc vật lý

Mô hình con lắc vật lý đã cho thấy sự tương đồng giữa kết quả giải số và giải tích, từ đó khẳng định tính chính xác của phương pháp giải số.

4.2. Ứng dụng trong môi trường thực tế

Kết quả nghiên cứu cũng chỉ ra rằng các mô hình có thể được áp dụng trong các tình huống thực tế, như trong thiết kế các cấu trúc kỹ thuật.

V. Kết luận và định hướng tương lai trong nghiên cứu cơ học vật rắn

Nghiên cứu về cơ học vật rắn vẫn còn nhiều tiềm năng phát triển. Việc áp dụng các phương pháp giải số sẽ tiếp tục được mở rộng để giải quyết các bài toán phức tạp hơn trong tương lai.

5.1. Tương lai của nghiên cứu cơ học vật rắn

Nghiên cứu sẽ tiếp tục tập trung vào việc phát triển các phương pháp giải số mới, nhằm nâng cao độ chính xác và hiệu quả trong việc giải quyết các bài toán cơ học.

5.2. Ứng dụng công nghệ mới trong nghiên cứu

Công nghệ mới như trí tuệ nhân tạo và máy học có thể được tích hợp vào nghiên cứu cơ học vật rắn, mở ra nhiều hướng đi mới cho các nghiên cứu trong tương lai.

Nội dung chính

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

Trong phân môn cơ học cổ điển, cơ học vật rắn đóng vai trò nền tảng khi nghiên cứu trạng thái chuyển động và biến dạng của các vật thể dưới tác dụng của ngoại lực, sự biến đổi nhiệt độ cùng các yếu tố tương tác cơ học. Phạm vi ứng dụng của lý thuyết này trải rộng từ nghiên cứu vật lý cơ bản đến các phân ngành kỹ thuật ứng dụng như cơ khí chế tạo, công trình xây dựng và kỹ thuật đại dương.

Phương pháp truyền thống để giải quyết các bài toán cơ học vật rắn chủ yếu dựa vào công cụ giải tích toán học (đạo hàm, tích phân, giải phương trình vi phân vi tích phân). Lối tiếp cận này cho phép thiết lập các công thức tường minh mô tả mối liên hệ giữa các đại lượng động học và động lực học. Tuy nhiên, khi đưa các điều kiện thực tế vào mô hình dao động—chẳng hạn lực cản môi trường phụ thuộc phi tuyến vào vận tốc (tỉ lệ với bình phương vận tốc) hoặc ngoại lực biến đổi phức tạp theo thời gian như luồng gió thổi—việc giải phương trình vi phân bằng phương pháp giải tích giải được nghiệm chính xác giải tích trở nên bất khả thi hoặc gặp nhiều khó khăn. Trong bối cảnh đó, phương pháp giải số trên máy tính bằng các ngôn ngữ lập trình khoa học (Fortran, Python, Matlab, C) trở thành giải pháp hữu hiệu nhằm xác định tọa độ, vận tốc và gia tốc của cơ hệ thông qua các bảng số liệu và đồ thị mô phỏng.

Mục tiêu và nhiệm vụ nghiên cứu của khóa luận được xác định cụ thể như sau:

1. Hệ thống hóa các cơ sở lý thuyết động học và động lực học vật rắn, tập trung vào chuyển động quay quanh trục cố định và quy trình 4 bước giải bài toán động lực học vật rắn.
2. Trình bày cơ sở toán học của bài toán giá trị ban đầu (bài toán Cauchy), phương pháp Runge-Kutta bậc 2 (RK2) và bậc 4 (RK4), xây dựng thuật toán giải hệ phương trình vi phân cấp một tương đương.
3. Ứng dụng ngôn ngữ lập trình Fortran và phương pháp RK4 để giải số và phân tích các bài toán dao động cơ học vật rắn từ cơ bản đến phức tạp có xét yếu tố cản và ngoại lực thực tế.

Đối tượng và phạm vi nghiên cứu:

• Đối tượng nghiên cứu: Quy luật chuyển động quay và dao động của vật rắn quanh trục cố định, bao gồm con lắc đơn, con lắc vật lý dạng thanh thẳng, mô hình thiết bị antenna đo tín hiệu điện từ dưới đáy đại dương và con lắc vật lý chịu lực cản không khí cùng luồng gió thổi.
• Phạm vi không gian: Khảo sát chuyển động trong mặt phẳng hai chiều (2D), hệ tọa độ cực và tọa độ Descartes.
• Phạm vi thời gian và điều kiện: Mô phỏng dao động trong các khoảng thời gian hữu hạn từ 10 giây, 20 giây đến 60 giây đầu tiên kể từ thời điểm kích thích ban đầu.

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

Khóa luận sử dụng hệ thống lý thuyết cơ học vật rắn cổ điển kết hợp lý thuyết giải tích số cho bài toán Cauchy. Cơ cấu chuyển động quay của vật rắn quanh trục cố định được xây dựng dựa trên sự tương đồng với chuyển động tịnh tiến, thông qua các đại lượng góc, moment quán tính và định luật hai Newton cho chuyển động quay.

Chuyển động tịnh tiến	Chuyển động quay quanh trục cố định	Mối liên hệ tiếp tuyến
Tọa độ vị trí: $x(t)$	Tọa độ góc: $\theta(t)$	$s = R\theta$
Vận tốc: $v(t) = dx/dt$	Vận tốc góc: $\omega(t) = d\theta/dt$	$v = R\omega$
Gia tốc: $a(t) = dv/dt = d^2x/dt^2$	Gia tốc góc: $\gamma(t) = d\omega/dt = d^2\theta/dt^2$	$a_t = R\gamma$
Khối lượng quán tính: $m$	Moment quán tính: $I = \sum m_i \rho_i^2$	$I_O = I_{cm} + md^2$ (Steiner-Huygens)
Định luật 2 Newton: $\sum \vec{F} = m\vec{a}$	Phương trình quay: $\sum \tau_z = I_z \gamma = I_z \frac{d^2\theta}{dt^2}$	$\vec{\tau} = \vec{r} \times \vec{F}$

Để xác định moment quán tính khi thay đổi vị trí trục quay song song với trục đi qua khối tâm, tác giả sử dụng định lý Steiner-Huygens: $$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.

Vật thể đồng chất	Vị trí trục quay qua khối tâm	Moment quán tính ($I_{cm}$)
Thanh thẳng chiều dài $L$	Trục vuông góc với thanh tại trung điểm	$\frac{1}{12}ML^2$
Vòng dây hoặc hình trụ rỗng bán kính $R$	Trục đối xứng hình học	$MR^2$
Đĩa tròn hoặc hình trụ đặc bán kính $R$	Trục đối xứng hình học	$\frac{1}{2}MR^2$
Quả cầu đặc bán kính $R$	Trục đi qua tâm quả cầu	$\frac{2}{5}MR^2$
Quả cầu rỗng bán kính $R$	Trục đi qua tâm quả cầu	$\frac{2}{3}MR^2$

Về phương pháp tính toán, khóa luận tiếp cận bài toán vi phân bậc hai có điều kiện ban đầu: $$\frac{d^2\theta}{dt^2} = f\left(t, \theta, \frac{d\theta}{dt}\right), \quad \theta(t_0) = \theta_0, \quad \omega(t_0) = \omega_0$$

Quy trình giải số được thực hiện bằng cách hạ bậc phương trình vi phân cấp hai về hệ hai phương trình vi phân cấp một với $N = 2$, đặt biến trạng thái $x_1(t) = \theta(t)$ và $x_2(t) = \omega(t)$: $$\frac{dx_1}{dt} = x_2(t)$$ $$\frac{dx_2}{dt} = f(t, x_1(t), x_2(t))$$

Đoạn thời gian khảo sát $[a, b]$ được chia 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 Runge-Kutta bậc 4 (RK4) được lập trình để tính toán giá trị nghiệm tại bước $i+1$ dựa trên bốn hệ số góc đánh giá tại các điểm giữa và cuối bước chia: $$y_{i+1} = y_i + \frac{1}{6}(k_1 + 2k_2 + 2k_3 + k_4)h$$ với: $$k_1 = f(t_i, y_i)$$ $$k_2 = f\left(t_i + \frac{h}{2}, y_i + \frac{h}{2}k_1\right)$$ $$k_3 = f\left(t_i + \frac{h}{2}, y_i + \frac{h}{2}k_2\right)$$ $$k_4 = f(t_i + h, y_i + hk_3)$$

Nguồn dữ liệu tính toán được thiết lập từ các thông số thực nghiệm và mô hình chuẩn hóa trong các tài liệu tham khảo chuyên khảo về tính toán cơ học và đại dương học. Toàn bộ chương trình tính toán số được tác giả hiện thực hóa bằng ngôn ngữ Fortran với chương trình con RK4 và hàm xác định đạo hàm DERIV.

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 các định luật và công thức của cơ học vật rắn. Mở đầu bằng chuyển động quay của chất điểm và vật rắn, chương phân tích động học quay với các biến số vị trí góc $\theta(t)$, vận tốc góc $\omega(t)$ và gia tốc góc $\gamma(t)$. Khái niệm moment quán tính được trình bày tường minh cho cả hệ rời rạc và vật thể liên tục thông qua tích phân khối lượng riêng, khẳng định tính chất cộng của moment quán tính theo nguyên lý chồng chất. Định lý dời trục Steiner-Huygens được diễn giải toán học và tổng hợp thành bảng giá trị cho 5 dạng vật thể đối xứng cơ bản.

Ở phần động lực học vật rắn, chương xây dựng phương trình chuyển động quay dựa trên moment lực $\vec{\tau} = \vec{r} \times \vec{F}$ và moment động lượng $\vec{L} = \vec{r} \times \vec{p}$. Tác giả đúc kết quy trình giải bài toán động lực học vật rắn thành 4 bước:

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

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

Chương 2 trình bày chi tiết cơ sở toán học của phương pháp sai phân số để giải bài toán giá trị ban đầu Cauchy. Tác giả chứng minh cách chuyển đổi một phương trình vi phân cấp hai thành hệ phương trình vi phân cấp một. Khái niệm chuỗi Taylor được sử dụng để phân tích sai số của các thuật toán: ngắt chuỗi Taylor ở bậc nhất dẫn đến phương pháp Euler với sai số cắt cụt địa phương cấp $O(h^2)$; ngắt chuỗi ở bậc hai dẫn đến thuật toán Runge-Kutta bậc 2 (RK2).

Từ việc cân bằng hệ số của hàm gia đa bước, chương dẫn xuất chi tiết công thức tích phân số của phương pháp Runge-Kutta bậc 4 (RK4) đạt bậc chính xác địa phương $O(h^5)$ mà không đòi hỏi phải tính đạo hàm bậc cao. Phần cuối chương trình bày cấu trúc lập trình giải số bằng ngôn ngữ Fortran: sử dụng mảng biến temp(1) tương ứng với tọa độ góc $y$ và temp(2) tương ứng với vận tốc góc $y'$. Trọng tâm tính toán đặt vào việc xây dựng chính xác hàm DERIV cho từng bài toán cơ học cụ thể.

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

Chương 3 triển khai giải số cho 4 bài toán dao động cơ học vật rắn với các mức độ phức tạp tăng dần:

1. Con lắc đơn: Phương trình chuyển động $d^2\theta/dt^2 = -(g/l)\sin\theta$. Khi giả thiết góc dao động nhỏ ($\sin\theta \approx \theta$), cơ hệ dao động điều hòa với nghiệm giải tích $\theta_{ana} = \theta_0 \cos(\sqrt{g/l}t) = \theta_0 \cos(\pi t)$ (với tham số chọn $g/l = \pi^2$). Khóa luận tiến hành giải số RK4 với bước thời gian $tstep = 0.01$ s và so sánh trực tiếp với nghiệm giải tích ở các biên độ ban đầu khác nhau:

   • Tại $\theta_0 = \pi/20$ rad và $\pi/10$ rad: Nghiệm giải số và nghiệm giải tích hoàn toàn trùng khớp, khẳng định tính chính xác của thuật toán trong phạm vi góc nhỏ ($\theta < 10^\circ \approx \pi/18$ rad).
   • Tại $\theta_0 = \pi/6$ rad và $\pi/3$ rad: Xuất hiện sự sai lệch rõ nét giữa nghiệm giải số và nghiệm giải tích. Độ lệch pha và sai số biên độ tăng dần theo thời gian, chứng minh tính chất phi điều hòa khi biên độ góc lớn dù chuyển động vẫn bảo toàn tính tuần hoàn.

2. Con lắc vật lý dạng thanh thẳng: Xét một thanh đồng chất khối lượng $m = 1$ kg, chiều dài $d = 0.25$ m, có trục quay đi qua một đầu thanh. Khoảng cách khối tâm $l_{cm} = d/2 = 0.125$ m. Moment quán tính qua khối tâm $I_{cm} = \frac{1}{12}md^2$, áp dụng định lý Steiner-Huygens tính được moment quán tính đối với điểm treo $I_S = \frac{1}{3}md^2 = 1/48$ kg$\cdot$m$^2$. Lấy gia tốc trọng trường $g = \pi^2$, phương trình vi phân chuyển động chính xác là $d^2\theta/dt^2 = -6\pi^2 \sin\theta$. Chương trình Fortran khai báo hàm DERIV = -59.2176264*sin(temp(1)) với bước thời gian $tstep = 0.01$ s. Kết quả mô phỏng trong 10 giây đầu với $\theta_0 = \pi/10$ rad và $\theta_0 = \pi/6$ rad chỉ ra rằng ở góc lớn $\pi/6$ rad, chu kỳ dao động thực tế bị kéo dài so với chu kỳ điều hòa $T_0 = 2\pi/\sqrt{6}\pi$ s, thể hiện rõ tính chất phi tuyến của trường trọng lực.

3. Con lắc vật lý đặt trong môi trường chất lỏng (antenna dưới đáy đại dương): Khảo sát chuyển động của thiết bị antenna đo tín hiệu điện từ thăm dò hydrocarbon có chiều dài $L = 5$ m, treo tại điểm $M$. Gia tốc góc của hệ được tham số hóa từ thực nghiệm theo phương trình vi phân phi tuyến: $$\frac{d^2\theta}{dt^2} = -c_1 \sin\theta - c_2 \frac{d\theta}{dt}\left|\frac{d\theta}{dt}\right|$$ với các hệ số thực nghiệm $c_1 = 1.0$ (đại diện cho tác dụng kết hợp của trọng lực và lực đẩy Archimedes) và $c_2 = 0.2$ (hệ số lực cản môi trường chất lỏng tỉ lệ với bình phương vận tốc). Điều kiện ban đầu: $\theta_0 = 10^\circ = \pi/18$ rad, $\omega_0 = 0$. Mô phỏng giải số trong thời gian $t = 60$ s với bước thời gian $tstep = 0.01$ s cho thấy biên độ góc $\theta(t)$, vận tốc góc $\omega(t)$ và gia tốc góc $\gamma(t)$ đều suy giảm liên tục theo thời gian, tái hiện chính xác hiện tượng dao động tắt dần dưới tác dụng của lực cản phi tuyến $-c_2 |\omega|\omega$.

4. Con lắc vật lý dưới tác dụng của gió: Khảo sát quả cầu đặc bán kính $R = 0.02$ m, mật độ khối lượng $\rho_m = 7000$ kg/m$^3$, khối lượng $M = \frac{4}{3}\pi R^3 \rho_m$, treo bởi dây nhẹ chiều dài $L = 1$ m. Moment quán tính toàn phần $I = \frac{2}{5}MR^2 + ML^2$. Môi trường không khí có mật độ $\rho_k = 1.293$ kg/m$^3$, hệ số lực cản $D = 12\rho_k R^2$, gia tốc trọng trường $g = 9.8$ m/s$^2$. Lực cản không khí phụ thuộc vận tốc tương đối: $\vec{F}_D = -D|\vec{v} - \vec{u}|(\vec{v} - \vec{u})$. Khóa luận phân tích ba trường hợp:

   • Không có gió ($\vec{u} = \vec{0}$): Phương trình chuyển động có moment cản $\tau_{z,D} = -DL^3 |\omega|\omega$. Khảo sát trong 20 giây với các vận tốc góc ban đầu $\omega_{max} = 2, 6, 10$ rad/s. Tại $\omega_{max} = 2$ và $6$ rad/s, con lắc dao động tắt dần quanh vị trí cân bằng $\theta = 0$. Tại $\omega_{max} = 10$ rad/s, do cơ năng ban đầu rất lớn, quả cầu quay liên tục khoảng 4 vòng ($\theta \approx 8\pi \approx 25.12$ rad) trong 5 giây đầu tiên trước khi chuyển sang dao động tắt dần quanh vị trí cân bằng mới $\theta = 8\pi$.
   • Gió thổi đều theo phương ngang với vận tốc không đổi $u(t) = u_0 = 20$ m/s: Phương trình vi phân xuất hiện biểu thức vận tốc tương đối phức hợp: $$\gamma = -\frac{MgL}{I}\sin\theta - \frac{LD}{I}\sqrt{\omega^2 L^2 + u_0^2 - 2\omega L u_0 \cos\theta}(\omega L - u_0 \cos\theta)$$ Mô phỏng cho thấy quá trình tắt dần diễn ra nhanh hơn đáng kể so với khi không có gió. Vị trí cân bằng của con lắc bị lệch khỏi trục thẳng đứng, không còn nằm tại các vị trí $\theta = 2n\pi$, tương tự như bài toán con lắc đặt trong điện trường đều nằm ngang.
   • Gió thổi điều hòa theo thời gian: Vận tốc gió biến thiên theo quy luật $u(t) = \frac{u_0}{2}[1 + \cos(2\pi t / T)]$ với $u_0 = 20$ m/s, khảo sát hai chu kỳ gió $T = 5$ s và $T = 10$ s tại $\omega_{max} = 8$ rad/s. Sự phụ thuộc tường minh vào biến thời gian $t$ làm biến đổi sâu sắc chu kỳ dao động của con lắc. Con lắc không bị tắt dao động hoàn toàn mà duy trì chuyển động do liên tục nhận năng lượng kích thích tuần hoàn từ luồng gió. Tại chu kỳ $T = 10$ s, ở thời điểm $t = 15$ s biên độ dao động tiệm cận về 0, nhưng đến thời điểm $t = 20$ s biên độ lại tăng vọt trở lại khi vận tốc gió đạt giá trị cực đại.

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

Thông số và đặc tính	Con lắc đơn	Con lắc thanh thẳng	Antenna đáy đại dương	Quả cầu chịu lực gió
Kích thước / Khối lượng	Chất điểm, dây dài $l$	$m = 1$ kg, $d = 0.25$ m	$L = 5$ m	$M$ ($\rho_m = 7000$ kg/m$^3$), $R = 0.02$ m, $L = 1$ m
Moment quán tính ($I$)	$ml^2$	$I_S = \frac{1}{3}md^2 = \frac{1}{48}$ kg$\cdot$m$^2$	Xác định từ thực nghiệm	$I = \frac{2}{5}MR^2 + ML^2$
Lực cản môi trường	Bỏ qua	Bỏ qua	Lực cản chất lỏng: $-c_2 |\omega|\omega$	Lực cản không khí: $-D|\vec{v}-\vec{u}|(\vec{v}-\vec{u})$
Ngoại lực đặc biệt	Trọng lực	Trọng lực	Lực đẩy Archimedes ($c_1 = 1.0$)	Luồng gió: $u_0 = 20$ m/s hoặc $u(t)$ biến thiên
Phương pháp giải	So sánh giải tích & RK4	So sánh giải tích & RK4	Giải số RK4 ($tstep = 0.01$ s)	Giải số RK4 ($tstep = 0.01$ s, $nsteps = 2000$)
Hiện tượng vật lý chính	Mất tính điều hòa ở góc lớn	Chu kỳ tăng khi biên độ lớn	Tắt dần đồng thời $\theta(t), \omega(t), \gamma(t)$	Quay vòng khi kích thích lớn; lệch vị trí cân bằng; kích thích duy trì bởi gió điều hòa

Khóa luận đạt được các kết quả khoa học chính:

1. Đã xây dựng thành công bộ chương trình tính toán số bằng Fortran ứng dụng thuật toán Runge-Kutta bậc 4, giải quyết ổn định hệ phương trình vi phân cấp một tương đương cho các bài toán cơ học vật rắn phi tuyến.
2. Làm rõ giới hạn của phép gần đúng góc nhỏ ($\sin\theta \approx \theta$) thông qua định lượng sai số giữa nghiệm giải tích và nghiệm số ở các góc kích thích lớn ($\pi/6$ rad và $\pi/3$ rad).
3. Mô phỏng định lượng chuyển động dao động tắt dần của mô hình antenna đại dương có xét đồng thời lực đẩy Archimedes và lực cản bình phương vận tốc trong thời gian 60 giây.
4. Phát hiện và mô tả các hiệu ứng động lực học phi tuyến của con lắc chịu tác dụng của gió: hiện tượng quay tròn 4 vòng trước khi dao động khi vận tốc góc ban đầu đạt $10$ rad/s; sự dịch chuyển vị trí cân bằng động khi gió thổi đều $20$ m/s; và hiện tượng phục hồi biên độ dao động tại các thời điểm gió đạt vận tốc cực đại khi chịu tác động của gió điều hòa chu kỳ $T = 10$ s.

Đóng góp tự nêu của tác giả: Khóa luận đã hệ thống hóa phương pháp tiếp cận số hóa cho các mô hình cơ học vật rắn vượt ra ngoài khuôn khổ giải tích truyền thống của giáo trình vật lý đại cương, chứng minh tính hiệu quả của công cụ giải số trong việc xử lý các bài toán cơ học có lực cản phi tuyến và ngoại lực phụ thuộc thời gian.

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

Khóa luận tồn tại một số giới hạn nghiên cứu được tác giả nêu trong văn bản:

• Phạm vi mô hình mới dừng lại ở các bài toán con lắc dao động phẳng hai chiều (2D) quanh một trục quay cố định; chưa mở rộng sang không gian ba chiều (3D).
• Luồng gió tác dụng trong mô hình mới khảo sát trường hợp thổi theo phương ngang cố định (thổi đều hoặc điều hòa); chưa xét đến trường hợp gió thổi theo phương bất kỳ trong không gian hoặc gió đổi hướng liên tục theo thời gian.
• Chương trình tính toán mới tập trung giải quyết các phương trình vi phân bậc hai đưa về hệ hai phương trình vi phân bậc một ($N = 2$); chưa triển khai cho các hệ nhiều bậc tự do hoặc phương trình vi phân bậc cao hơn.

Định hướng phát triển tiếp theo bao gồm việc ứng dụng kết quả tính toán và đồ thị mô phỏng số vào đổi mới phương pháp dạy học môn Vật lý đại cương tại các trường đại học và trường sư phạm; đồng thời phát triển thuật toán giải số cho các hệ dao động phức tạp hơn chịu tác động của gió đa phương và môi trường cản biến tính.

Giá trị tham khảo và ứng dụng

Tài liệu này hữu ích cho:

• Giảng viên và sinh viên ngành Sư phạm Vật lý, Vật lý học: Làm tài liệu tham khảo trực tiếp cho học phần Cơ học thuộc chương trình Vật lý đại cương, minh họa trực quan sự sai lệch giữa gần đúng giải tích góc nhỏ và bản chất dao động phi tuyến ở góc lớn.
• Sinh viên và học viên chuyên ngành Cơ học kỹ thuật, Toán ứng dụng, Vật lý tính toán: Cung cấp tài liệu mẫu về quy trình thiết lập phương trình vi phân chuyển động từ định luật hai Newton, kỹ thuật hạ bậc phương trình vi phân và mã nguồn Fortran chuẩn RK4 để giải các bài toán giá trị ban đầu.
• Người nghiên cứu dao động công trình và mô hình thủy khí: Tham khảo phương pháp thiết lập lực cản phi tuyến tỉ lệ bình phương vận tốc tương đối và phương pháp tham số hóa các hệ số lực đẩy môi trường trong mô hình antenna đại dương.

Phần đáng tham khảo nhất là Chương 2 về thuật toán Runge-Kutta bậc 4 cùng cấu trúc lập trình hàm DERIV trong Fortran, và Mục 3.4 phân tích cơ học chuyển động của con lắc dưới tác dụng của luồng gió biến thiên.

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

1. Khóa luận sử dụng thuật toán giải số nào và viết trên ngôn ngữ lập trình gì? Khóa luận sử dụng phương pháp Runge-Kutta bậc 4 (RK4) để giải hệ hai phương trình vi phân bậc một tương đương với bước thời gian $tstep = 0.01$ s, toàn bộ chương trình tính toán được lập trình bằng ngôn ngữ Fortran.

2. Điểm khác biệt giữa nghiệm giải tích và nghiệm giải số của con lắc đơn thể hiện như thế nào ở các góc biên độ khác nhau? Ở biên độ góc nhỏ ($\theta_0 = \pi/20$ rad và $\pi/10$ rad), nghiệm giải tích $\theta(t) = \theta_0 \cos(\pi t)$ hoàn toàn trùng khớp với nghiệm giải số RK4. Khi tăng biên độ lên góc lớn ($\theta_0 = \pi/6$ rad và $\pi/3$ rad), nghiệm giải tích gần đúng góc nhỏ lệch đáng kể so với nghiệm số, độ lệch tăng dần theo thời gian, thể hiện hệ mất tính chất dao động điều hòa dù vẫn giữ tính tuần hoàn.

3. Trong mô hình antenna dưới đáy đại dương, ý nghĩa vật lý của các hệ số $c_1$ và $c_2$ là gì? Hệ số $c_1 = 1.0$ đại diện cho độ lớn của gia tốc biểu kiến sinh ra bởi lực đẩy Archimedes cùng phương nhưng ngược chiều với trọng lực tác dụng lên antenna; hệ số $c_2 = 0.2$ đặc trưng cho độ lớn của lực cản môi trường chất lỏng tỉ lệ với bình phương vận tốc thông qua số hạng $-c_2 |\omega|\omega$.

4. Hiện tượng gì xảy ra với con lắc chịu lực cản không khí khi kích thích ban đầu với vận tốc góc rất lớn $\omega_{max} = 10$ rad/s trong trường hợp không có gió? Do nhận được cơ năng ban đầu rất lớn, con lắc không dao động ngay mà quay liên tục khoảng 4 vòng tròn tương ứng với góc quét $\theta \approx 8\pi \approx 25.12$ rad trong 5 giây đầu tiên; sau đó, năng lượng bị tiêu tán bởi lực cản không khí và con lắc mới bắt đầu dao động tắt dần quanh vị trí cân bằng mới tại $\theta = 8\pi$.

5. Khi luồng gió thổi theo phương ngang có vận tốc không đổi $u_0 = 20$ m/s, trạng thái cân bằng của con lắc thay đổi như thế nào so với trường hợp không có gió? Khi không có gió, vị trí cân bằng của con lắc luôn nằm ở các góc $\theta = 2n\pi$. Khi có gió thổi đều $u_0 = 20$ m/s, lực đẩy của luồng gió làm vị trí cân bằng của con lắc bị lệch khỏi phương thẳng đứng và không còn nằm tại các góc $\theta = 2n\pi$, tương tự như hiện tượng con lắc đơn dao động trong điện trường đều có phương nằm ngang.

Kết luận

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

Trích đoạn nội dung tài liệu

Chương 1 Lý thuyết cơ học vật rắn Trong chương này, khóa luận sẽ trình bày tóm tắt về lý thuyết cơ học vật rắn. Đầu tiên, các đại lượng động học và trạng thái chuyển động trong chuyển động quay của một chất điểm được trình bày ngắn gọn. Chuyển động quay và chuyển động tịnh tiến có những tính chất và cách tiếp cận tương tự nhau. Để tiếp cận một bài toán cơ học vật rắn, ta thường tiếp cận bằng phương pháp phân tích lực hoặc phương pháp năng lượng.

Ngoài ra, ta cũng có thể tiếp cận dựa trên cơ học Lagrange hay Hamilton. Tuy nhiên, trong giới hạn của khóa luận này, chúng tôi không trình bày phần cơ học Lagrange và Hamilton ở đây. Nội dung chương này được tham khảo từ các tài liệu [1–3]. Chuyển động quay Ta có thể làm một vật rắn chuyển động bằng cách tịnh tiến hoặc quay nó.

Ta đã biết chuyển động tịnh tiến của một chất điểm thường được mô tả bởi vector tọa độ ⃗r trong không gian Cartesian. Trong khi đó, chuyển động quay của vật rắn được mô tả bởi góc θ(t). Cụ thể, trong không gian 2 chiều, nó được mô tả bởi góc θ, điểm quay O và hướng trục quay ⃗k. 3 Trong suốt quá trình quay, góc θ(t) biến đổi theo thời gian.

Vận tốc góc của vật khi đó là dθ ω(t) = .1) dt Nếu vận tốc góc biến đổi theo thời gian, khi đó gia tốc góc của vật được xác định như sau dω d2 θ γ(t) = = 2.2) dt dt Một điểm trên vật quay cách trục quay một khoảng R có vận tốc tiếp tuyến ds dθ v= = R = Rω.3) dt dt Bảng 1.1 thể hiện sự so sánh (tương ứng) giữa hai loại chuyển động quay và chuyển động tịnh tiến. Để giải các bài toán chuyển động quay thì tương tự như trong chuyển động tịnh tiến, ta cần giải phương trình vi phân cấp hai d2 θ = γ(t, θ, ω), (1.4) dt2 với điều kiện đầu: θ(t0 ) = θ0 và ω(t0 ) = ω0. Khi gia tốc góc γ = γ(t) là hàm phụ thuộc tường minh vào thời gian t Bảng 1.1: So sánh chuyển động tịnh tiến và chuyển động quay. Chuyển động Tịnh tiến Quay Vị trí x(t) θ(t) Vận tốc v(t) = dx/dt ω(t) = dθ/dt Gia tốc a(t) = dv/dt = d2 x/dt2 γ(t) = dω/dt = d2 θ/dt2 4 thì ta có thể giải phương trình vi phân bằng tích phân trực tiếp ˆ t ω(t) = ω(t0 ) + γ(t) dt, (1.6) t0 Trong trường hợp đặc biệt, với chuyển động có vận tốc góc không đổi, ta có: θ(t) = θ(t0 ) + ω(t − t0 ).7) Nếu chuyển động có gia tốc góc không đổi, ta có: ω(t) = ω(t0 ) + γ(t − t0 ), (1.9) 2 Tổng quát, chuyển động quay xảy ra quanh một trục được cho bởi vector vận tốc góc ω ⃗.

Khi đó, vận tốc của một điểm trên vật quay tại vị trí ⃗r là ⃗ × ⃗r.10) Gia tốc tại một điểm trên vật quay tại vị trí ⃗r là ⃗a = ⃗γ × ⃗r + ω ⃗ × (⃗ω × ⃗r). Chuyển động quay của vật rắn Một vật rắn có thể chuyển động tịnh tiến hoặc chuyển động quay nhưng khoảng cách giữa hai điểm bất kì trên vật rắn không thay đổi. Đại lượng đặc trưng cho quán tính trong chuyển động quay là moment 5 quán tính. Moment quán tính của một hệ quanh một trục là X IO = mi ρ2i , (1.12) i trong đó mi là khối lượng của chất điểm i và ρi là khoảng cách từ chất điểm i đến trục quay O.

Moment quán tính của vật rắn với khối lượng riêng ρM quanh trục O là ˚ IO = ρM ρ2 dV, (1.13) trong đó ρ là khoảng cách từ phần tử dV đến trục O. Các moment quán tính được cộng theo nguyên lý chồng chất: Moment quán tính của hai hệ A và B quanh trục O là tổng của các moment quán tính của mỗi hệ IO,AB = IO,A + IO,B .14) Khi thay đổi trục quay, moment quán tính được xác định thông qua định lý Steiner-Huygens (định lý trục song song hay định lý dời trục). Định lý Steiner-Huygens: Moment quán tính IO của một vật quanh trục O có thể được xác định từ moment quán tính của vật quanh một trục song song đi qua khối tâm Icm của vật IO = Icm + md2 , (1.15) trong đó d là khoảng cách giữa trục O và khối tâm (hình 1.2 đưa ra các công thức moment quán tính đối với trục quay đi qua khối tâm của một số vật rắn thường gặp có hình dạng đặc biệt. Tương tự như trong chuyển động tịnh tiến, chuyển động quay cũng có cơ năng bao gồm động năng quay và thế năng.

Minh họa định lý Steiner-Huygens. Động năng của một vật rắn quay quanh một trục O cố định là 1 K = IO ω 2 , (1.16) 2 trong đó ω là vận tốc góc quanh trục O. Thế năng của vật rắn trong trọng trường đều là U = mgz, (1.2: Moment quán tính của một số vật thể đồng chất có trục quay đi qua khối tâm. Vật thể Icm Thanh thẳng (1/12)M L2 Vòng dây / Trụ rỗng M R2 Đĩa / Trụ đặc (1/2)M R2 Quả cầu đặc (2/5)M R2 Quả cầu rỗng (2/3)M R2 7 trong đó z là vị trí theo phương thẳng đứng của khối tâm hệ so với mốc thế năng.

Động lực học vật rắn ⃗ tại điểm cách trục quay đi qua O một khoảng ⃗r là Moment của lực F ⃗τ = ⃗r × F⃗ .18) Các moment lực có thể được cộng để xác định moment tổng quanh một trục cho trước. Moment của trọng lực đối với trục quay bất kì tác dụng lên vật rắn là ⃗τ = ⃗r × m⃗g , (1.19) trong đó ⃗r là vị trí khối tâm của vật (vector nối từ trục quay đến khối tâm). Với một vật rắn quay quanh một trục cố định (trục Oz ) đi qua khối tâm, gia tốc tịnh tiến của khối tâm là X F⃗j = m⃗a, (1.20) j gia tốc góc quanh khối tâm là X τz,j = Iz γ, (1.21) j ⃗j là moment lực j tác dụng tại điểm ⃗rcm,j so với trong đó ⃗τj = ⃗rcm,j × F khối tâm, và Iz là moment quán tính của vật quanh trục Oz đi qua khối tâm. 8 Moment động lượng của một chất điểm với động lượng (tịnh tiến) p⃗ tại vị trí ⃗r đối với điểm O là ⃗ = ⃗r × p⃗.22) Moment động lượng của hệ quanh điểm O là X ⃗O = L ⃗ri × mi⃗vi .23) i Định luật hai Newton cho chuyển động quay của hệ quanh điểm O cố định là ⃗O X dL = ⃗ri × F⃗i = ⃗τO .24) dt i Moment động lượng của hệ quanh khối tâm của nó là X ⃗ cm = L ⃗rcm,i × mi⃗vi .25) i Định luật hai Newton cho chuyển động quay của hệ quanh khối tâm của nó là ⃗ cm X dL = ⃗rcm,i × F⃗i = ⃗τcm .26) dt i Moment động lượng của vật rắn quay quanh trục cố định là ⃗ O,z = IO,z ω L ⃗.27) Moment động lượng của vật rắn quay quanh một trục cố định đi qua khối tâm là ⃗ cm,z = Icm,z ω L ⃗.28) Ta có thể tóm tắt quy trình giải một bài toán động lực trong cơ học vật rắn gồm 4 bước [1] 9 1.

Nhận dạng: Xác định góc θ(t) (chọn hệ quy chiếu). Tìm các điều kiện đầu θ(t0 ) và ω(t0 ). Mô hình: Tìm các lực tác dụng lên vật rắn. Xác định moment lực cho mỗi lực.

Áp dụng định luật hai Newton cho chuyển động quay để tìm γ = γ(θ, ω, t). Giải: Ta giải phương trình vi phân d2 θ = γ(θ, ω, t), (1.29) dt2 với các điều kiện đầu. Ta có thể giải bằng các kĩ thuật giải tích hoặc giải số. Phân tích: Kiểm tra tính đúng đắn của θ(t) và ω(t).

Đánh giá quá trình chuyển động. 10 Chương 2 Phương pháp tính toán Chương này sẽ trình bày cách giải một phương trình vi phân với các điều kiện ban đầu. Với các phương trình đơn giản, ta có thể dễ dàng tìm được nghiệm giải tích của chúng. Tuy nhiên một số phương trình vi phân không thể hoặc chưa có lời giải chính xác, khi đó ta cần đến sự tính toán của máy tính với các phương pháp khác nhau nhằm đưa ra nghiệm giải số.

Từ đó, các đại lượng vật lý trong hệ vật rắn có thể được biểu diễn bằng đồ thị. Trong khóa luận này, phương pháp Runge-Kutta bậc 4 được sử dụng để đưa ra nghiệm giải số. Nội dung chương này được tham khảo từ tài liệu [4, 5]. Bài toán Cauchy Một phương trình vi phân bậc hai có thể được viết dưới dạng y ′′ (t) = f (t, y(t), y ′ (t)), (2.1) trong đó ta có thể tìm được hàm y(t) từ các đạo hàm của nó.

Một phương trình vi phân có vô số nghiệm riêng cấu thành nên nghiệm tổng quát, chúng khác nhau bởi một hằng số. Với các hằng số được xác 11 định, ta thu được nghiệm riêng từ nghiệm tổng quát. Khi cho trước giá trị ban đầu y0 của hàm y(t) tại t0 ta sẽ nhận được một nghiệm riêng của phương trình. Bài toán Cauchy còn được gọi là bài toán giá trị ban đầu.

Xét phương trình vi phân cấp hai thỏa 2 điều kiện được biểu diễn như sau  ′′ ′  y (t) = f (t, y(t), y (t))   y(a) = y0 , a ≤ t ≤ b, (2. Bất kì một phương trình vi phân bậc cao luôn có thể đưa được về hệ phương trình vi phân bậc nhất. Phương trình vi phân cấp hai (2.2) có thể đưa về phương trình vi phân cấp 1 bằng cách biến đổi sau y(t) = x1 (t), (2.5) với điều kiện đầu của phương trình (2. Với bất kì một phương pháp giải số nào, ta luôn phải rời rạc hóa các thành phần trong phương trình.

Để tìm nghiệm gần đúng của bài toán (2.2) ta chia đoạn [a, b] thành n đoạn nhỏ bằng nhau với khoảng chia ∆t ≡ h = (b − a)/n. Khi đó các điểm chia là t0 = a, tk = t0 + k · h với k = 0, 1, 2,. Giá trị gần đúng cần tìm của hàm tại điểm tk , được kí hiệu là yk. Phương pháp Runge-Kutta Trong kĩ thuật giải số, phương pháp Runge-Kutta được sử dụng để gần đúng nghiệm của các phương trình vi phân tuyến tính.

Phương pháp này được phát triển từ những năm đầu của thế kỉ trước bởi hai nhà toán học người Đức Carl Runge và Wilhelm Kutta. Phương pháp Euler là phương pháp giải số bậc nhất cho phương trình vi phân toàn phần với một điều kiện đầu. Nó thuộc nhóm phương pháp Runge-Kutta đơn giản với bậc thấp nhất [6]. Ở các bậc cao hơn, ta có các phương pháp Runge-Kutta bậc 2 và bậc 4.