Một mét chìm thêm tạo ra bao nhiêu lực phục hồi?
Heave là chuyển động tịnh tiến theo phương thẳng đứng của vật nổi. Muốn mô hình hóa chuyển động này, trước hết phải biết lực nổi thay đổi ra sao khi thân dịch khỏi vị trí cân bằng. Tài liệu QBlade trình bày độ cứng heave tuyến tính bằng khối lượng riêng, gia tốc trọng trường và diện tích mặt đường nước, đồng thời giới hạn phép tuyến tính hóa khi hình học phần ngập thay đổi đáng kể. [1]
Bài này xét một vật nổi được giữ không nghiêng, chỉ có một bậc tự do thẳng đứng. Các kích thước, khối lượng và hệ số đều là giả định minh họa. Không suy rộng kết quả thành đáp ứng sóng của một sà lan hoặc giàn cụ thể; ở đó còn có ghép heave–pitch, kích động sóng, bức xạ và các lực khác.
Chọn gốc tại trạng thái cân bằng
Lấy y dương xuống dưới, tính từ vị trí cân bằng ban đầu trong nước yên. Khối lượng thực m=1.200.000 kg, nước có ρ=1.025 kg/m³ và g=9,81 m/s². Giả định lực neo tăng thêm bằng không tại gốc, nên thể tích chiếm nước cân bằng là m/ρ=1.170,731707 m³. Nếu có lực kéo neo thường trực, phải đưa nó vào cân bằng tĩnh trước khi đặt gốc.
Diện tích mặt đường nước ban đầu A₀=48 m². Khi y nhỏ và thành thân gần thẳng đứng, thể tích ngập tăng A₀y; lực nổi tăng ρgA₀y hướng lên. Vì y dương xuống, lực đó mang dấu phục hồi âm trong phương trình lực.
C_h=ρgA₀=482.652 N/m.
Giả định hệ giữ vị trí có độ cứng thẳng đứng gia tăng k_m=100.000 N/m, tuyến tính quanh gốc. Độ cứng tổng C=582.652 N/m. Hệ giữ này được mô hình hóa như một phần tử đàn hồi đã biết, không coi mọi kiểu neo ngoài khơi đều tạo ra cùng một độ cứng thẳng đứng.
Khối lượng thêm tham gia quán tính, không tự tăng trọng lượng
Giả định khối lượng thêm heave a₃₃=800.000 kg trong miền chuyển động đang xét. Khối lượng quán tính hiệu dụng M=m+a₃₃=2.000.000 kg. Không đưa M vào phương trình cân bằng lực nổi như thể có thêm 800 tấn hàng nằm trên boong; khối lượng thêm biểu diễn phản ứng quán tính của nước đối với chuyển động.
Bài tập heave và roll của MIT phân biệt khối lượng thân với khối lượng thêm trong tính tần số riêng và phương trình động lực. [2] Với mô hình vô hướng hệ số hằng ở đây, phương trình quanh cân bằng là:
Mÿ+bẏ+Cy=P(t).
P(t) là lực ngoài gia tăng hướng xuống, đơn vị N; b có đơn vị N·s/m. Các dấu chấm biểu diễn đạo hàm theo thời gian. Hạng Cy là lực phục hồi đã chuyển sang vế trái, không phải một lực tiếp tục làm vật chìm xuống.
Chu kỳ và quá trình thả tự do
Tần số góc riêng không cản là ω_n=√(C/M)=0,539746 rad/s; chu kỳ T_n=2π/ω_n=11,640999 s. Giả định tỷ số cản ζ=0,08, khi đó b=2ζ√(MC)=172.718,795735 N·s/m. Tần số có cản ω_d=ω_n√(1−ζ²)=0,538016 rad/s, cho chu kỳ 11,678430 s.
Đẩy vật xuống y₀=0,30 m rồi thả từ trạng thái đứng yên, không có P(t) sau khi thả. Nghiệm tuyến tính đúng điều kiện vận tốc ban đầu bằng không là:
y(t)=y₀e^(−ζω_nt)[cos(ω_dt)+ζ sin(ω_dt)/√(1−ζ²)].
Tại t=3 s, y=0,009736 m, tức vật gần quay lại gốc. Chỉ viết y₀e^(−ζω_nt)cos(ω_dt) sẽ cho vận tốc ban đầu khác không; lỗi nhỏ ở biểu thức ban đầu có thể làm sai thời điểm qua mức tham chiếu. Biên bao suy giảm cũng không đồng nghĩa mọi thời điểm nằm trên đường biên ấy.
Thông số b ở đây là tổng cản tuyến tính giả định. Tài liệu Orcina nêu rằng dữ liệu khối lượng thêm và cản bức xạ có thể phụ thuộc tần số; dùng giá trị hằng cần phù hợp với miền tần số chuyển động. [3] Vì vậy một cặp M, b suy từ một thí nghiệm không tự đại diện cho mọi trạng thái biển.
Mặt đường nước mở rộng khi chìm xuống
Bây giờ giữ nguyên trạng thái gốc nhưng cho diện tích mặt đường nước biến đổi theo quy luật giả định A(y)=48+12y m², với y tính bằng mét và chỉ áp dụng trong khoảng hình học đang xét. Hệ số 12 có đơn vị mét. Độ tăng thể tích phải lấy tích phân, không lấy A(y)y:
ΔV=48y+6y².
Lực phục hồi tổng trở thành R(y)=Cy+qy², với q=ρg×6=60.331,5 N/m². Tại y=1 m, R=642.983,5 N, lớn hơn 60.331,5 N so với mô hình giữ diện tích ban đầu. Dùng A(y)y sẽ tính phần bậc hai gấp đôi giá trị đúng.
Giả định một cơ cấu ngoài tác dụng lực xuống không đổi P=300.000 N nhưng không thêm khối lượng vào vật. Cân bằng mới thỏa qy_e²+Cy_e=P. Lấy nghiệm dương phù hợp miền hình học:
y_e=[√(C²+4qP)−C]/(2q).
Kết quả y_e=0,490023 m. Mô hình tuyến tính cho P/C=0,514887 m, sâu hơn khoảng 24,864 mm. Sai khác này xuất phát từ hình học tăng diện tích mặt đường nước, không phải một hệ số an toàn được tự thêm vào.
Chu kỳ quanh cân bằng mới dùng độ cứng tiếp tuyến
Nếu chỉ dao động rất nhỏ quanh y_e, độ cứng phù hợp là đạo hàm của lực phục hồi:
C_t=C+2qy_e=641.779,676450 N/m.
Giữ M không đổi theo giả thiết, chu kỳ không cản cục bộ là 11,091796 s. Không dùng P/y_e làm độ cứng tiếp tuyến: đó là độ cứng dây cung từ gốc cũ đến cân bằng mới. Nếu lực 300 kN thực ra do bổ sung hàng, khối lượng thực cũng tăng; phải sửa cả M và vị trí cân bằng, không chỉ thay độ cứng.
Sự thay đổi mớn còn có thể làm khối lượng thêm và cản thay đổi. Việc giữ chúng cố định trong phép so sánh này nhằm tách riêng ảnh hưởng của mặt đường nước, chứ không khẳng định thủy động lực độc lập với mớn.
Kiểm năng lượng cho một lần thả không cản
Xét một thí nghiệm giả định khác: bỏ cản và lực ngoài, dùng mô hình R(y)=Cy+qy², rồi thả vật từ y=0,50 m với vận tốc bằng không. Thế năng phục hồi so với gốc là:
U(y)=Cy²/2+qy³/3.
Tại vị trí thả, U=75.345,3125 J. Khi vật lần đầu qua y=0, toàn bộ phần này chuyển thành động năng trong mô hình bảo toàn; tốc độ là √(2U/M)=0,274491 m/s. Nếu dùng lực tuyến tính, kết quả là 0,269873 m/s. Hai số không được đem so với nghiệm có cản ở phần trước như thể cùng một phép thử.
Phép kiểm năng lượng giúp bắt lỗi dấu và tích phân. Nó cũng cho thấy giới hạn của việc gọi một chu kỳ duy nhất cho dao động lớn: lực bậc hai làm hai phía của gốc không đối xứng. Khi mặt boong chạm nước, pontoon nổi khỏi nước, dây neo đổi nhánh hoặc có va đập, quy luật đa thức hiện tại có thể hết hiệu lực.
Đối với hồ sơ kỹ thuật, cần lưu trạng thái cân bằng, mặt đường nước theo mớn, nguồn hệ số thủy động lực và điều kiện của phép thử. Chu kỳ heave chỉ là một đặc trưng của mô hình quanh trạng thái đó; khả năng làm việc ngoài khơi còn phụ thuộc phổ sóng, lực kích động, chuyển động tương đối và giới hạn của thiết bị được phục vụ.
Nguồn tham khảo
[1] Hydrostatic Forces — hydrostatic stiffness matrix — QBlade Documentation. Truy cập ngày 09/11/2026.
[2] Floating Structure Heave and Roll — section 41 — MIT OpenCourseWare, 2.017J. Truy cập ngày 09/11/2026.
[3] Vessel types: Stiffness, added mass and damping — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Heave là chuyển động tịnh tiến theo phương thẳng đứng của vật nổi. Muốn mô hình hóa chuyển động này, trước hết phải biết lực nổi thay đổi ra sao khi thân dịch khỏi vị trí cân bằng. Tài liệu QBlade trình bày độ cứng heave tuyến tính bằng khối lượng riêng, gia tốc trọng trường và diện tích mặt đường nước, đồng thời giới hạn phép tuyến tính hóa khi hình học phần ngập thay đổi đáng kể. [1]
Bài này xét một vật nổi được giữ không nghiêng, chỉ có một bậc tự do thẳng đứng. Các kích thước, khối lượng và hệ số đều là giả định minh họa. Không suy rộng kết quả thành đáp ứng sóng của một sà lan hoặc giàn cụ thể; ở đó còn có ghép heave–pitch, kích động sóng, bức xạ và các lực khác.
Chọn gốc tại trạng thái cân bằng
Lấy y dương xuống dưới, tính từ vị trí cân bằng ban đầu trong nước yên. Khối lượng thực m=1.200.000 kg, nước có ρ=1.025 kg/m³ và g=9,81 m/s². Giả định lực neo tăng thêm bằng không tại gốc, nên thể tích chiếm nước cân bằng là m/ρ=1.170,731707 m³. Nếu có lực kéo neo thường trực, phải đưa nó vào cân bằng tĩnh trước khi đặt gốc.
Diện tích mặt đường nước ban đầu A₀=48 m². Khi y nhỏ và thành thân gần thẳng đứng, thể tích ngập tăng A₀y; lực nổi tăng ρgA₀y hướng lên. Vì y dương xuống, lực đó mang dấu phục hồi âm trong phương trình lực.
C_h=ρgA₀=482.652 N/m.
Giả định hệ giữ vị trí có độ cứng thẳng đứng gia tăng k_m=100.000 N/m, tuyến tính quanh gốc. Độ cứng tổng C=582.652 N/m. Hệ giữ này được mô hình hóa như một phần tử đàn hồi đã biết, không coi mọi kiểu neo ngoài khơi đều tạo ra cùng một độ cứng thẳng đứng.
Khối lượng thêm tham gia quán tính, không tự tăng trọng lượng
Giả định khối lượng thêm heave a₃₃=800.000 kg trong miền chuyển động đang xét. Khối lượng quán tính hiệu dụng M=m+a₃₃=2.000.000 kg. Không đưa M vào phương trình cân bằng lực nổi như thể có thêm 800 tấn hàng nằm trên boong; khối lượng thêm biểu diễn phản ứng quán tính của nước đối với chuyển động.
Bài tập heave và roll của MIT phân biệt khối lượng thân với khối lượng thêm trong tính tần số riêng và phương trình động lực. [2] Với mô hình vô hướng hệ số hằng ở đây, phương trình quanh cân bằng là:
Mÿ+bẏ+Cy=P(t).
P(t) là lực ngoài gia tăng hướng xuống, đơn vị N; b có đơn vị N·s/m. Các dấu chấm biểu diễn đạo hàm theo thời gian. Hạng Cy là lực phục hồi đã chuyển sang vế trái, không phải một lực tiếp tục làm vật chìm xuống.
Chu kỳ và quá trình thả tự do
Tần số góc riêng không cản là ω_n=√(C/M)=0,539746 rad/s; chu kỳ T_n=2π/ω_n=11,640999 s. Giả định tỷ số cản ζ=0,08, khi đó b=2ζ√(MC)=172.718,795735 N·s/m. Tần số có cản ω_d=ω_n√(1−ζ²)=0,538016 rad/s, cho chu kỳ 11,678430 s.
Đẩy vật xuống y₀=0,30 m rồi thả từ trạng thái đứng yên, không có P(t) sau khi thả. Nghiệm tuyến tính đúng điều kiện vận tốc ban đầu bằng không là:
y(t)=y₀e^(−ζω_nt)[cos(ω_dt)+ζ sin(ω_dt)/√(1−ζ²)].
Tại t=3 s, y=0,009736 m, tức vật gần quay lại gốc. Chỉ viết y₀e^(−ζω_nt)cos(ω_dt) sẽ cho vận tốc ban đầu khác không; lỗi nhỏ ở biểu thức ban đầu có thể làm sai thời điểm qua mức tham chiếu. Biên bao suy giảm cũng không đồng nghĩa mọi thời điểm nằm trên đường biên ấy.
Thông số b ở đây là tổng cản tuyến tính giả định. Tài liệu Orcina nêu rằng dữ liệu khối lượng thêm và cản bức xạ có thể phụ thuộc tần số; dùng giá trị hằng cần phù hợp với miền tần số chuyển động. [3] Vì vậy một cặp M, b suy từ một thí nghiệm không tự đại diện cho mọi trạng thái biển.
Mặt đường nước mở rộng khi chìm xuống
Bây giờ giữ nguyên trạng thái gốc nhưng cho diện tích mặt đường nước biến đổi theo quy luật giả định A(y)=48+12y m², với y tính bằng mét và chỉ áp dụng trong khoảng hình học đang xét. Hệ số 12 có đơn vị mét. Độ tăng thể tích phải lấy tích phân, không lấy A(y)y:
ΔV=48y+6y².
Lực phục hồi tổng trở thành R(y)=Cy+qy², với q=ρg×6=60.331,5 N/m². Tại y=1 m, R=642.983,5 N, lớn hơn 60.331,5 N so với mô hình giữ diện tích ban đầu. Dùng A(y)y sẽ tính phần bậc hai gấp đôi giá trị đúng.
Giả định một cơ cấu ngoài tác dụng lực xuống không đổi P=300.000 N nhưng không thêm khối lượng vào vật. Cân bằng mới thỏa qy_e²+Cy_e=P. Lấy nghiệm dương phù hợp miền hình học:
y_e=[√(C²+4qP)−C]/(2q).
Kết quả y_e=0,490023 m. Mô hình tuyến tính cho P/C=0,514887 m, sâu hơn khoảng 24,864 mm. Sai khác này xuất phát từ hình học tăng diện tích mặt đường nước, không phải một hệ số an toàn được tự thêm vào.
Chu kỳ quanh cân bằng mới dùng độ cứng tiếp tuyến
Nếu chỉ dao động rất nhỏ quanh y_e, độ cứng phù hợp là đạo hàm của lực phục hồi:
C_t=C+2qy_e=641.779,676450 N/m.
Giữ M không đổi theo giả thiết, chu kỳ không cản cục bộ là 11,091796 s. Không dùng P/y_e làm độ cứng tiếp tuyến: đó là độ cứng dây cung từ gốc cũ đến cân bằng mới. Nếu lực 300 kN thực ra do bổ sung hàng, khối lượng thực cũng tăng; phải sửa cả M và vị trí cân bằng, không chỉ thay độ cứng.
Sự thay đổi mớn còn có thể làm khối lượng thêm và cản thay đổi. Việc giữ chúng cố định trong phép so sánh này nhằm tách riêng ảnh hưởng của mặt đường nước, chứ không khẳng định thủy động lực độc lập với mớn.
Kiểm năng lượng cho một lần thả không cản
Xét một thí nghiệm giả định khác: bỏ cản và lực ngoài, dùng mô hình R(y)=Cy+qy², rồi thả vật từ y=0,50 m với vận tốc bằng không. Thế năng phục hồi so với gốc là:
U(y)=Cy²/2+qy³/3.
Tại vị trí thả, U=75.345,3125 J. Khi vật lần đầu qua y=0, toàn bộ phần này chuyển thành động năng trong mô hình bảo toàn; tốc độ là √(2U/M)=0,274491 m/s. Nếu dùng lực tuyến tính, kết quả là 0,269873 m/s. Hai số không được đem so với nghiệm có cản ở phần trước như thể cùng một phép thử.
Phép kiểm năng lượng giúp bắt lỗi dấu và tích phân. Nó cũng cho thấy giới hạn của việc gọi một chu kỳ duy nhất cho dao động lớn: lực bậc hai làm hai phía của gốc không đối xứng. Khi mặt boong chạm nước, pontoon nổi khỏi nước, dây neo đổi nhánh hoặc có va đập, quy luật đa thức hiện tại có thể hết hiệu lực.
Đối với hồ sơ kỹ thuật, cần lưu trạng thái cân bằng, mặt đường nước theo mớn, nguồn hệ số thủy động lực và điều kiện của phép thử. Chu kỳ heave chỉ là một đặc trưng của mô hình quanh trạng thái đó; khả năng làm việc ngoài khơi còn phụ thuộc phổ sóng, lực kích động, chuyển động tương đối và giới hạn của thiết bị được phục vụ.
Nguồn tham khảo
[1] Hydrostatic Forces — hydrostatic stiffness matrix — QBlade Documentation. Truy cập ngày 09/11/2026.
[2] Floating Structure Heave and Roll — section 41 — MIT OpenCourseWare, 2.017J. Truy cập ngày 09/11/2026.
[3] Vessel types: Stiffness, added mass and damping — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
