Vì sao phải biết bước sóng khi tính công trình ngoài khơi?
Thông tin “sóng cao 2 m, chu kỳ 8 giây” chưa đủ để xác định cách tải phân bố lên một kết cấu. Bước sóng quyết định sự khác pha giữa hai vị trí cách nhau theo phương truyền sóng và tham gia vào phân bố vận tốc nước theo chiều sâu. Nếu nhập nhầm bước sóng nước sâu cho một vùng có độ sâu hữu hạn, sai số có thể đi tiếp vào bài toán tải, chuyển động và tương tác giữa các bộ phận công trình.
Bài này xem một thành phần sóng điều hòa trên đáy phẳng, không có dòng chảy nền. Tài liệu sóng tuyến tính của MIT thiết lập quan hệ giữa tần số góc, số sóng và độ sâu từ các điều kiện biên của bài toán mặt thoáng. [1] Phần tính số bên dưới là ví dụ tự xây dựng nhằm kiểm tra phương pháp, không phải bộ thông số khí tượng hải văn cho một mỏ dầu khí.
Quan hệ tán sắc và các ký hiệu
Công thức cần giải là:
ω² = g × k × tanh(kh)
Trong đó ω = 2π/T là tần số góc, đơn vị rad/s; T là chu kỳ, s; k = 2π/L là số sóng, rad/m; L là bước sóng, m; h là độ sâu nước, m; g là gia tốc trọng trường, m/s². Tích kh không có thứ nguyên. Trong phép tính, radian được xem là không thứ nguyên, nhưng vẫn nên ghi rad/s và rad/m để tránh nhầm với tần số Hz hoặc nghịch đảo bước sóng.
Khi đã biết T và h, k xuất hiện cả bên ngoài lẫn bên trong hàm tanh. Vì vậy không thể lấy k = ω²/g rồi gọi đó là nghiệm tổng quát. Giá trị này chỉ là giới hạn nước sâu. Tài liệu MIT về nghiệm quan hệ tán sắc phân biệt giới hạn nước sâu, nước nông và vùng trung gian cần giải phương trình đầy đủ. [2]
Có thể đặt x = kh và μ = ω²h/g để đưa bài toán về x tanh(x) = μ. Cách này thuận lợi cho kiểm tra vì chỉ còn một phương trình không thứ nguyên. Khi x > 0, đạo hàm của vế trái là tanh(x) + x/cosh²(x), luôn dương. Do đó, với μ > 0 có một nghiệm dương duy nhất. Tính đơn điệu cũng giúp phương pháp chia đôi hoạt động ổn định nếu khoảng ban đầu bao lấy nghiệm.
Cách giải có thể kiểm tra lại
Đặt F(k) = gk tanh(kh) − ω². Chọn hai giá trị dương ka và kb sao cho F(ka) < 0 và F(kb) > 0. Lấy trung điểm, tính dấu của F rồi giữ nửa khoảng còn chứa nghiệm. Tiếp tục đến khi khoảng đủ nhỏ theo mục đích tính. Đây là thủ tục tìm nghiệm số cho ví dụ; không phụ thuộc một phần mềm thương mại cụ thể.
Nếu dùng Newton, công thức lặp có thể viết k mới = k − F(k)/F′(k), với F′(k) = g[tanh(kh) + kh/cosh²(kh)]. Cần kiểm tra k luôn dương và kết quả cuối thực sự thỏa phương trình. Việc phần mềm báo “hội tụ” chỉ cho biết tiêu chí lặp đã đạt; một phép đổi đơn vị sai ở đầu vào vẫn có thể tạo nghiệm hội tụ rất đẹp của một bài toán khác.
Nên kiểm tra hai sai số riêng. Sai số phương trình có thể dùng ε = |gk tanh(kh) − ω²|/ω². Sai số dữ liệu lại đến từ chu kỳ, độ sâu và giả thiết vật lý. Đẩy ε xuống cực nhỏ không làm giảm bất định của số liệu khảo sát. Với bài viết hướng dẫn, giữ nhiều chữ số trong tính nội bộ rồi làm tròn kết quả trình bày thường hợp lý hơn in ra một bước sóng dài hàng chục chữ số.
Ví dụ giả định ở độ sâu 20 m
Lấy T = 8 s, h = 20 m và g = 9,81 m/s². Khi đó ω = 0,785398 rad/s và ω² = 0,616850 s⁻². Giải phương trình cho k xấp xỉ 0,07076243 rad/m; tích kh bằng 1,41525. Bước sóng L = 2π/k = 88,79267 m, làm tròn 88,79 m. Vận tốc pha c = L/T = 11,09908 m/s.
Kiểm tra ngược bằng số sóng chưa làm tròn: 2π/√[gk tanh(kh)] phải trả lại chu kỳ 8 s. Đây là phép kiểm tra độc lập với việc chỉ lấy L chia T. Nếu kiểm tra ngược không đạt, cần xem lại hàm lượng giác hyperbolic, độ sâu và hệ đơn vị trước khi sử dụng kết quả cho vận tốc hạt nước.
Theo xấp xỉ nước sâu, L∞ = gT²/(2π) = 99,92384 m. Giá trị này lớn hơn nghiệm hữu hạn khoảng 12,54%. Theo xấp xỉ nước nông, Lsh = T√(gh) = 112,05713 m, lớn hơn khoảng 26,20%. Cả hai công thức đều đơn giản, nhưng ở bộ dữ liệu này việc chọn tùy ý một giới hạn tạo sai số rõ rệt. Dấu hiệu cần chú ý là kh không đủ lớn để tanh(kh) gần một và cũng không đủ nhỏ để thay tanh(kh) bằng kh.
Độ nhạy theo chu kỳ và độ sâu
Giữ h = 20 m, tăng T từ 8 lên 10 s cho nghiệm L = 121,23691 m, tăng khoảng 36,54%, trong khi chu kỳ chỉ tăng 25%. Nếu thay bằng công thức nước sâu, L tăng theo T²; ở độ sâu hữu hạn sự phụ thuộc không còn đúng hoàn toàn dạng bình phương. Vì vậy không nên nội suy bước sóng đơn giản theo tỷ lệ chu kỳ khi khoảng thay đổi đáng kể.
Giữ T = 8 s nhưng giảm độ sâu từ 20 xuống 10 m, nghiệm trở thành L = 70,89835 m. Vận tốc pha giảm còn 8,86229 m/s. Đây là so sánh hai trạng thái trong mô hình tuyến tính, chưa phải mô phỏng toàn bộ quá trình sóng đi qua đáy biển biến đổi. Muốn mô tả tuyến truyền thực, còn cần phương, địa hình và các cơ chế biến đổi năng lượng tương ứng.
Để đánh giá một khoảng bất định nhỏ hơn, có thể thử T = 7,2 và 8,8 s tại h = 20 m. Hai bước sóng lần lượt khoảng 75,37 và 101,96 m. Như vậy sai khác ±10% về chu kỳ tạo khoảng bước sóng không đối xứng quanh giá trị cơ sở. Trong hồ sơ phân tích tải, nên giữ lại các trường hợp này thay vì nén mọi bất định vào một bước sóng trung bình.
Giới hạn và cách chuyển kết quả sang mô hình tải
Quan hệ trên thuộc lý thuyết sóng biên độ nhỏ. Nó không tự chứng minh một trường hợp sóng lớn, gần vỡ hay có dòng mạnh vẫn được mô tả phù hợp. Orcina lưu ý các lý thuyết sóng điều hòa đều cần thận trọng khi sóng vỡ hoặc gần vỡ. [3] Bài tính không đưa ra một ngưỡng chấp nhận kỹ thuật để thay cho tiêu chí của dự án.
Trước khi truyền k sang mô hình tải, cần thống nhất hệ cao độ: h được đo từ mực nước nào, có tính triều và nước dâng trong trường hợp đang xét hay không. Cần phân biệt vận tốc truyền đỉnh sóng c với vận tốc hạt nước tác dụng lên thanh kết cấu. Cùng một giá trị 11,10 m/s trong ví dụ, nếu bị nhập nhầm thành vận tốc dòng để tính lực cản sẽ dẫn đến một bài toán hoàn toàn khác.
Kết quả nên lưu cùng T, h, g, phương pháp tìm nghiệm và sai số kiểm tra. Khi bộ thông số khí tượng hải văn thay đổi, chỉ cần chạy lại chuỗi tính có truy vết ấy. Một phép tính tán sắc ngắn nhưng rõ giả thiết sẽ hữu ích hơn một bảng bước sóng không nói được nó thuộc độ sâu và chu kỳ nào.
Nguồn tham khảo
[1] 13.012 Reading 6: Linear Free Surface Waves — Massachusetts Institute of Technology — A. H. Techet. Truy cập ngày 09/11/2026.
[2] Solution of the Dispersion Relationship — Massachusetts Institute of Technology. Truy cập ngày 09/11/2026.
[3] Waves: Breaking waves — Orcina — OrcaFlex documentation. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Thông tin “sóng cao 2 m, chu kỳ 8 giây” chưa đủ để xác định cách tải phân bố lên một kết cấu. Bước sóng quyết định sự khác pha giữa hai vị trí cách nhau theo phương truyền sóng và tham gia vào phân bố vận tốc nước theo chiều sâu. Nếu nhập nhầm bước sóng nước sâu cho một vùng có độ sâu hữu hạn, sai số có thể đi tiếp vào bài toán tải, chuyển động và tương tác giữa các bộ phận công trình.
Bài này xem một thành phần sóng điều hòa trên đáy phẳng, không có dòng chảy nền. Tài liệu sóng tuyến tính của MIT thiết lập quan hệ giữa tần số góc, số sóng và độ sâu từ các điều kiện biên của bài toán mặt thoáng. [1] Phần tính số bên dưới là ví dụ tự xây dựng nhằm kiểm tra phương pháp, không phải bộ thông số khí tượng hải văn cho một mỏ dầu khí.
Quan hệ tán sắc và các ký hiệu
Công thức cần giải là:
ω² = g × k × tanh(kh)
Trong đó ω = 2π/T là tần số góc, đơn vị rad/s; T là chu kỳ, s; k = 2π/L là số sóng, rad/m; L là bước sóng, m; h là độ sâu nước, m; g là gia tốc trọng trường, m/s². Tích kh không có thứ nguyên. Trong phép tính, radian được xem là không thứ nguyên, nhưng vẫn nên ghi rad/s và rad/m để tránh nhầm với tần số Hz hoặc nghịch đảo bước sóng.
Khi đã biết T và h, k xuất hiện cả bên ngoài lẫn bên trong hàm tanh. Vì vậy không thể lấy k = ω²/g rồi gọi đó là nghiệm tổng quát. Giá trị này chỉ là giới hạn nước sâu. Tài liệu MIT về nghiệm quan hệ tán sắc phân biệt giới hạn nước sâu, nước nông và vùng trung gian cần giải phương trình đầy đủ. [2]
Có thể đặt x = kh và μ = ω²h/g để đưa bài toán về x tanh(x) = μ. Cách này thuận lợi cho kiểm tra vì chỉ còn một phương trình không thứ nguyên. Khi x > 0, đạo hàm của vế trái là tanh(x) + x/cosh²(x), luôn dương. Do đó, với μ > 0 có một nghiệm dương duy nhất. Tính đơn điệu cũng giúp phương pháp chia đôi hoạt động ổn định nếu khoảng ban đầu bao lấy nghiệm.
Cách giải có thể kiểm tra lại
Đặt F(k) = gk tanh(kh) − ω². Chọn hai giá trị dương ka và kb sao cho F(ka) < 0 và F(kb) > 0. Lấy trung điểm, tính dấu của F rồi giữ nửa khoảng còn chứa nghiệm. Tiếp tục đến khi khoảng đủ nhỏ theo mục đích tính. Đây là thủ tục tìm nghiệm số cho ví dụ; không phụ thuộc một phần mềm thương mại cụ thể.
Nếu dùng Newton, công thức lặp có thể viết k mới = k − F(k)/F′(k), với F′(k) = g[tanh(kh) + kh/cosh²(kh)]. Cần kiểm tra k luôn dương và kết quả cuối thực sự thỏa phương trình. Việc phần mềm báo “hội tụ” chỉ cho biết tiêu chí lặp đã đạt; một phép đổi đơn vị sai ở đầu vào vẫn có thể tạo nghiệm hội tụ rất đẹp của một bài toán khác.
Nên kiểm tra hai sai số riêng. Sai số phương trình có thể dùng ε = |gk tanh(kh) − ω²|/ω². Sai số dữ liệu lại đến từ chu kỳ, độ sâu và giả thiết vật lý. Đẩy ε xuống cực nhỏ không làm giảm bất định của số liệu khảo sát. Với bài viết hướng dẫn, giữ nhiều chữ số trong tính nội bộ rồi làm tròn kết quả trình bày thường hợp lý hơn in ra một bước sóng dài hàng chục chữ số.
Ví dụ giả định ở độ sâu 20 m
Lấy T = 8 s, h = 20 m và g = 9,81 m/s². Khi đó ω = 0,785398 rad/s và ω² = 0,616850 s⁻². Giải phương trình cho k xấp xỉ 0,07076243 rad/m; tích kh bằng 1,41525. Bước sóng L = 2π/k = 88,79267 m, làm tròn 88,79 m. Vận tốc pha c = L/T = 11,09908 m/s.
Kiểm tra ngược bằng số sóng chưa làm tròn: 2π/√[gk tanh(kh)] phải trả lại chu kỳ 8 s. Đây là phép kiểm tra độc lập với việc chỉ lấy L chia T. Nếu kiểm tra ngược không đạt, cần xem lại hàm lượng giác hyperbolic, độ sâu và hệ đơn vị trước khi sử dụng kết quả cho vận tốc hạt nước.
Theo xấp xỉ nước sâu, L∞ = gT²/(2π) = 99,92384 m. Giá trị này lớn hơn nghiệm hữu hạn khoảng 12,54%. Theo xấp xỉ nước nông, Lsh = T√(gh) = 112,05713 m, lớn hơn khoảng 26,20%. Cả hai công thức đều đơn giản, nhưng ở bộ dữ liệu này việc chọn tùy ý một giới hạn tạo sai số rõ rệt. Dấu hiệu cần chú ý là kh không đủ lớn để tanh(kh) gần một và cũng không đủ nhỏ để thay tanh(kh) bằng kh.
Độ nhạy theo chu kỳ và độ sâu
Giữ h = 20 m, tăng T từ 8 lên 10 s cho nghiệm L = 121,23691 m, tăng khoảng 36,54%, trong khi chu kỳ chỉ tăng 25%. Nếu thay bằng công thức nước sâu, L tăng theo T²; ở độ sâu hữu hạn sự phụ thuộc không còn đúng hoàn toàn dạng bình phương. Vì vậy không nên nội suy bước sóng đơn giản theo tỷ lệ chu kỳ khi khoảng thay đổi đáng kể.
Giữ T = 8 s nhưng giảm độ sâu từ 20 xuống 10 m, nghiệm trở thành L = 70,89835 m. Vận tốc pha giảm còn 8,86229 m/s. Đây là so sánh hai trạng thái trong mô hình tuyến tính, chưa phải mô phỏng toàn bộ quá trình sóng đi qua đáy biển biến đổi. Muốn mô tả tuyến truyền thực, còn cần phương, địa hình và các cơ chế biến đổi năng lượng tương ứng.
Để đánh giá một khoảng bất định nhỏ hơn, có thể thử T = 7,2 và 8,8 s tại h = 20 m. Hai bước sóng lần lượt khoảng 75,37 và 101,96 m. Như vậy sai khác ±10% về chu kỳ tạo khoảng bước sóng không đối xứng quanh giá trị cơ sở. Trong hồ sơ phân tích tải, nên giữ lại các trường hợp này thay vì nén mọi bất định vào một bước sóng trung bình.
Giới hạn và cách chuyển kết quả sang mô hình tải
Quan hệ trên thuộc lý thuyết sóng biên độ nhỏ. Nó không tự chứng minh một trường hợp sóng lớn, gần vỡ hay có dòng mạnh vẫn được mô tả phù hợp. Orcina lưu ý các lý thuyết sóng điều hòa đều cần thận trọng khi sóng vỡ hoặc gần vỡ. [3] Bài tính không đưa ra một ngưỡng chấp nhận kỹ thuật để thay cho tiêu chí của dự án.
Trước khi truyền k sang mô hình tải, cần thống nhất hệ cao độ: h được đo từ mực nước nào, có tính triều và nước dâng trong trường hợp đang xét hay không. Cần phân biệt vận tốc truyền đỉnh sóng c với vận tốc hạt nước tác dụng lên thanh kết cấu. Cùng một giá trị 11,10 m/s trong ví dụ, nếu bị nhập nhầm thành vận tốc dòng để tính lực cản sẽ dẫn đến một bài toán hoàn toàn khác.
Kết quả nên lưu cùng T, h, g, phương pháp tìm nghiệm và sai số kiểm tra. Khi bộ thông số khí tượng hải văn thay đổi, chỉ cần chạy lại chuỗi tính có truy vết ấy. Một phép tính tán sắc ngắn nhưng rõ giả thiết sẽ hữu ích hơn một bảng bước sóng không nói được nó thuộc độ sâu và chu kỳ nào.
Nguồn tham khảo
[1] 13.012 Reading 6: Linear Free Surface Waves — Massachusetts Institute of Technology — A. H. Techet. Truy cập ngày 09/11/2026.
[2] Solution of the Dispersion Relationship — Massachusetts Institute of Technology. Truy cập ngày 09/11/2026.
[3] Waves: Breaking waves — Orcina — OrcaFlex documentation. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
