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

Việt Nam sở hữu bờ biển dài trên 3.200 km cùng nguồn lợi thủy hải sản phong phú, đưa kinh tế biển trở thành một trong những ngành kinh tế mũi nhọn. Để đáp ứng yêu cầu khai thác xa bờ, đội tàu cá Việt Nam đã có sự chuyển biến mạnh mẽ về công suất: từ mức bình quân 15 CV (năm 1985) nâng lên khoảng 90 CV – 150 CV vào giai đoạn 2006 và định hướng đóng mới hàng loạt tàu có công suất từ 200 CV trở lên. Tàu thủy là công trình nổi hoạt động trong môi trường sóng gió phức tạp; việc đảm bảo tính năng an toàn chống chìm lật và ổn định thân tàu là yêu cầu kỹ thuật sống còn nhằm giảm thiểu tai nạn trên biển.

Trong lý thuyết và thực tiễn thiết kế tàu, bài toán tính toán cánh tay đòn ổn định ($l_\theta$) giữ vai trò quyết định. Đa số các phương pháp truyền thống (phương pháp trực tiếp, phương pháp Viện sĩ Crưlôp...) tiếp cận theo hướng bài toán thuận (chỉ kiểm tra ổn định sau khi đã hoàn thiện tuyến hình tàu), phụ thuộc vào việc đo đạc thủ công trên bản vẽ đường hình lý thuyết tàu (ĐHLTT), dẫn đến khối lượng tính toán lớn, sai số tích tụ nhiều và khó tự động hóa. Đề tài được thực hiện nhằm nghiên cứu, so sánh hai hướng tiếp cận giải bài toán nghịch hiện đại của Giáo sư V.G. Vlaxôv và Phó Giáo sư - Tiến sĩ Nguyễn Quang Minh, kết hợp thuật toán Spline để nâng cao độ chính xác tính toán.

Mục tiêu và nhiệm vụ nghiên cứu

  1. Tìm hiểu, tổng hợp cơ sở lý thuyết về tính nổi và tính ổn định tàu thủy;
  2. Nghiên cứu cơ sở giải tích và thuật toán tính tay đòn ổn định của GS. Vlaxôv và PGS.TS. Nguyễn Quang Minh;
  3. Nghiên cứu và xây dựng thuật toán Spline bậc ba nhằm hàm hóa đường hình lý thuyết thân tàu;
  4. Lựa chọn tàu mẫu, tính toán cánh tay đòn ổn định tĩnh theo hai phương pháp trên phần mềm máy tính;
  5. Xây dựng đồ thị ổn định tĩnh (đồ thị Reed), phân tích và so sánh sai số, khối lượng tính toán cũng như khả năng lập trình;
  6. Đánh giá sai số diện tích và mômen của các dạng đường cong hình học giữa thuật toán Spline và phần mềm đồ họa AutoCAD;
  7. Đưa ra nhận xét, kết luận và đề xuất khả năng ứng dụng thực tế cho đội tàu cá Việt Nam.

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

  • Đối tượng nghiên cứu: Phương pháp tính toán cánh tay đòn ổn định tàu thủy của GS. Vlaxôv và mô hình tàu ổn định tương đương của PGS.TS. Nguyễn Quang Minh.
  • Phạm vi nghiên cứu: Bài toán tính tay đòn ổn định tĩnh ($l_\theta$) cho tàu thủy trong điều kiện nước tĩnh.
  • Khảo sát thực tế và thời gian thực hiện: Chuyến đi thực tế thu thập số liệu tại Quảng Ngãi (01 tuần); thời gian thực hiện đề tài từ tháng 08/2006 đến tháng 11/2006 tại Khoa Cơ khí - Trường Đại học Nha Trang.

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

Khung lý thuyết và tiêu chuẩn ổn định

  • Lý thuyết ổn định tĩnh thân tàu: Tàu chịu tác động của mômen nghiêng tĩnh làm lệch phương giữa trọng lực $P$ và lực nổi $D$, hình thành mômen hồi phục: $$M_{hp} = P \cdot l_\theta = P \cdot (l_{hd} - l_{tl})$$ Trong đó:

    • $l_{hd} = y_c \cos\theta + (z_c - z_{co})\sin\theta$: Tay đòn ổn định hình dáng, phụ thuộc vào hình dạng hình học phần chìm của vỏ tàu.
    • $l_{tl} = z_g \sin\theta$: Tay đòn ổn định trọng lượng, phụ thuộc vào phân bố tải trọng và cao độ trọng tâm tàu $z_g$.
    • Ở góc nghiêng nhỏ ($\theta < 10^\circ - 12^\circ$), tay đòn hồi phục xác định qua chiều cao tâm ổn định ban đầu: $l_{hp} = h_0 \cdot \theta$ với $h_0 = r_0 - (z_g - z_c)$.
  • Tiêu chuẩn ổn định:

    • Tiêu chuẩn vật lý (Tiêu chuẩn Liên bang Nga): Đánh giá qua hệ số an toàn $K = \frac{M_{gh}}{M_{ng}} > 1$, trong đó $M_{gh}$ là mômen giới hạn chống lật và $M_{ng} = 0,001 \cdot z \cdot P_g \cdot S$ là mômen gió.
    • Tiêu chuẩn thống kê (IMO 1974 / Phương pháp Rahola): Áp dụng hệ 6 điều kiện an toàn, tiêu biểu gồm $h_0 \ge 0,35\text{ m}$ và $l_{\theta 30} \ge 0,030\text{ m}$.
    • Tiêu chuẩn tương đương cho tàu cá Việt Nam: PGS.TS. Nguyễn Quang Minh chứng minh đối với tàu cá ven bờ các tỉnh phía Nam, tiêu chuẩn $l_{\theta 30} \ge 0,20\text{ m}$ tương đương với toàn bộ hệ tiêu chuẩn IMO 1974, từ đó dẫn xuất tiêu chuẩn hình học đơn giản: $\frac{B}{H} \ge \frac{29 - 2B}{10}$.

Phương pháp nghiên cứu và giải thuật tính toán

Phương pháp Bản chất giải thuật Ưu điểm Hạn chế
Phương pháp trực tiếp Đo đạc tọa độ từ ĐHLTT, tính tích phân diện tích ($\omega$), thể tích ($V$) và mômen ($M_V$) bằng công thức hình thang. Trực quan, quen thuộc trong tính toán truyền thống. Khối lượng tính toán lớn, sai số đo đạc cao, khó lập trình tự động.
Phương pháp Viện sĩ Crưlôp Dựa trên định lý Euler về đường nước đẳng tích nghiêng vô cùng bé ($\Delta\theta = 5^\circ - 10^\circ$), xác định quỹ tích tâm nổi $C(y_c, z_c)$ qua bán kính $r_0 = \frac{J_x}{V}$. Cơ sở toán học chặt chẽ, chính xác ở góc nghiêng nhỏ. Sai số tích tụ qua nhiều bước lặp, quy trình thủ công phức tạp ở góc nghiêng lớn.
Phương pháp GS. Vlaxôv Tiếp cận bài toán nghịch bằng cách xấp xỉ tay đòn hình dáng $l_{hd}(\theta)$ theo chuỗi Fourier dạng đa thức lượng giác lẻ: $l_{hd} = a_1\sin\theta + a_2\sin 2\theta + a_3\sin 4\theta + a_4\sin 6\theta$. Rút gọn thành: $l_{hd} = y_{c90} f_1(\theta) + (z_{c90} - z_{co}) f_2(\theta) + r_0 f_3(\theta) + r_{90} f_4(\theta)$. Khối lượng tính toán giảm mạnh; chỉ cần xác định 5 thông số quy đổi tại $\theta = 0^\circ$ và $\theta = 90^\circ$ ($z_{co}, z_{c90}, y_{c90}, r_0, r_{90}$). Phụ thuộc vào các hệ số và bảng tra hàm thực nghiệm $f_i(\theta)$.
Phương pháp PGS.TS. Nguyễn Quang Minh Sử dụng mô hình tàu ổn định tương đương qua phép biến hình afin mở rộng: quy đổi tàu thực 3D về mô hình thân ống phẳng 2D có tiết diện là mặt cắt ngang (MCN) giữa, liên hệ qua hệ thức $Z_{ct} = Z_E$ và $Y_{ct} = Y_E \cdot \frac{2\eta}{y}$. Chuyển bài toán 3D phức tạp về khảo sát hình học MCN giữa; loại bỏ sai số đo đạc thủ công; tối ưu hóa cho lập trình máy tính. Đòi hỏi thuật toán giải tích chính xác để hàm hóa đường biên sườn tàu.
  • Thuật toán Spline bậc ba: Rời rạc hóa đường cong sườn tàu thành $(n-2)$ đoạn cong bậc ba $y = a + bx + cx^2 + dx^3$. Áp dụng điều kiện liên tục hàm, liên tục đạo hàm bậc nhất ($y'$) và đạo hàm bậc hai ($y''$) tại $(n-3)$ điểm nối, kết hợp điều kiện biên tiếp tuyến ban đầu. Thiết lập hệ phương trình ma trận cấp $4(n-2)$: $[A][X] = [B]$ và giải tìm hệ số trên phần mềm MATLAB.

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

Chương 1: Đặt vấn đề

Chương 1 khái quát thực trạng và định hướng phát triển của ngành tàu cá Việt Nam, chỉ rõ tầm quan trọng của tính năng ổn định đối với an toàn hàng hải. Tác giả phân tích ý nghĩa khoa học và thực tiễn của việc chuyển đổi từ bài toán thuận sang bài toán nghịch trong thiết kế tàu. Phần tổng quan hệ thống hóa các phương pháp tính toán ổn định kinh điển (phương pháp trực tiếp, phương pháp Crưlôp), đồng thời xác lập mục tiêu, đối tượng, phạm vi nghiên cứu và kế hoạch thực hiện đề tài.

Chương 2: Cơ sở lý thuyết về hai phương pháp và thuật toán hàm hoá

Chương 2 trình bày chi tiết cơ sở toán học của bản vẽ đường hình lý thuyết tàu, các hệ số hình học ($\alpha, \beta, \delta, \chi, \phi$) và các hệ thống tiêu chuẩn ổn định hiện hành (vật lý, thống kê IMO, tiêu chuẩn Việt Nam). Trọng tâm của chương là phân tích sâu:

  1. Giải thuật xấp xỉ chuỗi Fourier của GS. Vlaxôv;
  2. Cơ sở biến hình afin và mô hình tàu tương đương của PGS.TS. Nguyễn Quang Minh;
  3. Quy trình toán học xây dựng thuật toán Spline bậc ba để tính các tích phân diện tích MCN ($\omega = \int y dz$), diện tích mặt đường nước ($S = \int y dx$), thể tích chiếm nước ($V = \int \omega dx$), mômen tĩnh ($M_{\omega oy}, M_{\omega oz}, M_{Soy}, M_{Sox}$) và mômen quán tính ($J_x$);
  4. Thử nghiệm hàm hóa Spline trên 4 dạng đường cong hình học đặc trưng và kiểm chứng sai số với phần mềm AutoCAD.

Chương 3: Lựa chọn tàu mẫu, tính chính xác tay đòn theo hai phương pháp và so sánh kết quả

Chương 3 tiến hành chọn tàu mẫu thực tế (tàu cá vỏ gỗ/thép), ứng dụng các đoạn chương trình MATLAB để giải ma trận hệ số Spline cho từng đường sườn và đường nước. Tác giả tính toán chi tiết tọa độ trọng tâm, tâm nổi và cánh tay đòn ổn định $l_\theta$ tại các góc nghiêng từ $0^\circ$ đến $90^\circ$ theo cả hai phương pháp (Vlaxôv và Nguyễn Quang Minh). Kết quả tính toán được biểu diễn trên đồ thị Reed để so sánh độ dốc, giá trị cực đại ($l_{max}$) và đánh giá khả năng tự động hóa khi lập trình.


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

Kết quả thử nghiệm và đánh giá sai số thuật toán Spline

Để chứng minh độ chính xác của giải thuật hàm hóa Spline trước khi áp dụng tính toán thủy tĩnh, tác giả đã tiến hành khảo nghiệm trên các dạng đường cong thân tàu điển hình (cho bởi tập hợp 7 điểm gián đoạn từ $Z = 0$ đến $Z = 2100\text{ mm}$, thiết lập hệ ma trận 20 phương trình), sau đó tính tích phân diện tích ($\omega_{sp}$) và mômen tĩnh ($M_{\omega sp}$) để so sánh với giá trị đo trực tiếp chính xác trên AutoCAD ($\omega_{cx}, M_{\omega cx}$):

$$\delta_{dt} = \frac{|\omega_{cx} - \omega_{sp}|}{\omega_{cx}} \cdot 100%; \quad \delta_m = \frac{|M_{\omega cx} - M_{\omega sp}|}{M_{\omega cx}} \cdot 100%$$

Dạng đường cong Tọa độ điểm $Y$ (mm) ứng với $Z = (0, 350, 700, 1050, 1400, 1750, 2100)$ Diện tích xấp xỉ $\omega_{sp}$ ($\text{mm}^2$) Mômen tĩnh xấp xỉ $M_{\omega sp}$ ($\text{mm}^3$) Sai số diện tích $\delta_{dt}$ (%) Sai số mômen $\delta_m$ (%)
Dạng 1 (Mũi quả lê) $(0, 600, 650, 520, 570, 800, 1160)$ $1.325.045$ — $< 0,1%$ $0,011%$
Dạng 2 (Đường cong hai chiều) $(0, 55, 178, 400, 765, 1000, 1100)$ $1.031.103$ — $< 0,1%$ $0,037%$
Dạng 3 (Đường cong một chiều) $(0, 55, 178, 400, 765, 1200, 1700)$ $1.192.996$ — $< 0,1%$ $0,015%$
Dạng 4 (Đường cong mạn/đáy) $(0, 685, 1045, 1285, 1460, 1590, 1700)$ — — — —

Ghi chú: Kết quả kiểm chứng cho thấy sai số hình học của phương pháp hàm hóa Spline so với AutoCAD đều ở mức rất nhỏ ($< 0,05%$), đảm bảo độ tin cậy tuyệt đối cho các phép tính tích phân thủy tĩnh tiếp theo.

Đóng góp chính của đề tài

  1. Về mặt học thuật và phương pháp luận: Làm sáng tỏ bản chất toán học của hai phương pháp giải bài toán nghịch ổn định tàu thủy; chỉ rõ ưu thế của mô hình tàu tương đương (PGS.TS. Nguyễn Quang Minh) trong việc chuyển đổi tích phân không gian 3 chiều về tích phân phẳng trên MCN giữa.
  2. Về mặt kỹ thuật tính toán: Ứng dụng thành công công cụ MATLAB để giải hệ phương trình Spline bậc ba và thực hiện tích phân giải tích các yếu tố hình học tàu, khắc phục triệt để nhược điểm sai số tích tụ và sự cồng kềnh của phương pháp đo đạc đồ giải truyền thống.
  3. Về mặt thực tiễn: Cung cấp cơ sở tính toán và công cụ hỗ trợ thiết kế nhanh, tin cậy cho các kỹ sư đóng tàu khi tính toán kiểm tra ổn định các dòng tàu cá hoạt động tại ngư trường Việt Nam.

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

  • Hạn chế:
    • Nghiên cứu chỉ giới hạn trong bài toán ổn định tĩnh trên mặt nước phẳng tĩnh (nước tĩnh), chưa khảo sát ổn định động dưới tác động liên tục của sóng và gió giật.
    • Số lượng tàu mẫu thử nghiệm còn hạn chế trong khuôn khổ thời gian thực hiện đồ án tốt nghiệp đại học (hơn 03 tháng).
  • Hướng nghiên cứu tiếp:
    • Mở rộng thuật toán Spline và mô hình tàu tương đương để giải quyết bài toán ổn định động và chuyển động lắc ngang của tàu trên sóng;
    • Xây dựng phần mềm tích hợp hoàn chỉnh có giao diện người dùng, tự động hóa từ khâu nhập bảng tọa độ đường hình đến xuất hồ sơ ổn định theo tiêu chuẩn IMO và Đăng kiểm Việt Nam.

Giá trị tham khảo

  • Đối tượng sử dụng: Sinh viên, học viên cao học chuyên ngành Kỹ thuật Tàu thủy, Cơ khí Đóng tàu, Vỏ tàu thủy và Kỹ thuật Hàng hải.
  • Nội dung tham khảo hữu ích:
    • Tài liệu mẫu về phương pháp lập ma trận $[A][X] = [B]$ cho thuật toán Spline bậc ba giải trên MATLAB;
    • Cơ sở lý thuyết và các công thức quy đổi tọa độ tâm nổi của mô hình tàu ổn định tương đương;
    • Bảng tổng hợp công thức tính toán tích phân các yếu tố thủy tĩnh thân tàu ($\omega, S, V, M_{\omega}, M_S, J_x$).

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

1. Tác giả, người hướng dẫn và cơ sở đào tạo của đề tài tốt nghiệp này là ai?

Đề tài do sinh viên Nguyễn Xuân Nam (MSSV: 44D1194, Lớp 44TT, ngành Cơ khí đóng tàu thủy / Cơ khí tàu thuyền, Khoa Cơ khí - Trường Đại học Nha Trang) thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Quang Minh, hoàn thành vào tháng 11 năm 2006.

2. Điểm khác biệt cốt lõi giữa hai phương pháp Vlaxôv, Nguyễn Quang Minh so với các phương pháp tính ổn định truyền thống là gì?

Các phương pháp truyền thống tiếp cận theo bài toán thuận (chỉ kiểm tra ổn định sau khi đã định hình tàu bằng cách đo đạc thủ công trên bản vẽ ĐHLTT). Ngược lại, phương pháp của GS. Vlaxôv và PGS.TS. Nguyễn Quang Minh tiếp cận theo bài toán nghịch, đưa ra các mô hình xấp xỉ giải tích và mô hình ổn định tương đương, cho phép thiết kế chủ động các kích thước thân tàu thỏa mãn tiêu chuẩn ổn định ngay từ khâu nhiệm vụ thiết kế.

3. Điều kiện toán học để thiết lập hệ phương trình Spline bậc ba trong đề tài là gì?

Hệ phương trình Spline bậc ba $y = a + bx + cx^2 + dx^3$ được xây dựng dựa trên:

  • Điều kiện đường cong đi qua các điểm nút gián đoạn đã cho;
  • Điều kiện trơn và liên tục tại các điểm nối: đạo hàm bậc nhất bằng nhau ($y'{1} = y'{2}$) và đạo hàm bậc hai bằng nhau ($y''{1} = y''{2}$);
  • Điều kiện biên về tiếp tuyến/hệ số góc ($k$) tại điểm đầu xuất phát.

4. Tiêu chuẩn ổn định quy đổi nào được đề xuất cho đội tàu cá phía Nam Việt Nam trong tài liệu?

Dựa trên nghiên cứu của PGS.TS. Nguyễn Quang Minh được trích dẫn trong đề tài, đối với đội tàu cá ven bờ các tỉnh phía Nam, tiêu chuẩn tay đòn ổn định tại góc nghiêng $30^\circ$ là $l_{\theta 30} \ge 0,20\text{ m}$ được xem là tương đương với toàn bộ hệ 6 điều kiện của tiêu chuẩn IMO 1974, tương ứng với tiêu chuẩn kích thước hình học rút gọn là $\frac{B}{H} \ge \frac{29 - 2B}{10}$.

5. Độ chính xác của thuật toán hàm hóa Spline trong đồ án được kiểm chứng như thế nào?

Tác giả kiểm chứng bằng cách so sánh diện tích ($\omega_{sp}$) và mômen tĩnh ($M_{\omega sp}$) tính từ tích phân các hàm Spline với diện tích chính xác ($\omega_{cx}$) và mômen chính xác ($M_{\omega cx}$) đo trực tiếp trên phần mềm AutoCAD qua 4 dạng đường cong đặc trưng (mũi quả lê, đường cong hai chiều, đường cong một chiều và đường cong mạn/đáy), đạt sai số mômen rất nhỏ (dao động từ $0,011%$ đến $0,037%$).


Kết luận

Đề tài tốt nghiệp của tác giả Nguyễn Xuân Nam đã hệ thống hóa và phân tích sâu sắc cơ sở lý thuyết của hai phương pháp tính toán ổn định tàu thủy hiện đại theo hướng bài toán nghịch của GS. Vlaxôv và PGS.TS. Nguyễn Quang Minh. Việc ứng dụng thành công thuật toán Spline bậc ba kết hợp công cụ tính toán MATLAB đã giải quyết triệt để bài toán hàm hóa đường hình lý thuyết thân tàu, mang lại độ chính xác cao và giảm thiểu đáng kể khối lượng tính toán thủ công. Kết quả nghiên cứu cung cấp tài liệu tham khảo có giá trị học thuật và thực tiễn cao cho công tác thiết kế, kiểm tra tính ổn định của đội tàu cá Việt Nam.