Phổ sóng cần đi qua đáp ứng của đúng điểm đang quan tâm
Một biển có cùng chiều cao sóng có nghĩa có thể tạo chuyển động boong khác nhau khi năng lượng tập trung ở các tần số khác nhau. RAO mô tả đáp ứng theo một tần số và hướng sóng; nó gồm biên độ và pha. Orcina nhấn mạnh phải biết quy ước trục, hướng sóng, gốc pha và đơn vị của RAO trước khi ghép dữ liệu. [1]
Bài này tính một ví dụ giả định với ba dải tần, rồi chuyển từ heave tại gốc sang chuyển động đứng của điểm lệch khỏi gốc. Các con số không đại diện cho một tàu thật. Mô hình dùng hệ tuyến tính, trạng thái biển dừng và các thành phần tần số khác nhau không tương quan.
Vì sao phải bình phương môđun?
Nếu một thành phần sóng có biên độ a và hàm truyền phức H, biên độ đáp ứng là |H|a. Phương sai của một thành phần sin có biên độ a là a²/2; do đó phương sai đáp ứng nhân với |H|². Tài liệu waveresponse của 4Subsea viết phổ đầu ra dưới dạng tích H nhân liên hợp phức của H và phổ đầu vào, rồi tích phân theo hướng nếu có phổ định hướng. [2]
S_z(f)=|H_z(f)|²S_η(f).
Ở đây f tính bằng Hz, S_η có đơn vị m²/Hz và H_z có đơn vị m/m. Không dùng H² theo nghĩa bình phương một số phức, vì kết quả có thể phức và không phải mật độ phương sai. Cũng không chỉ nhân phổ với |H|, vì phổ là đại lượng bậc hai theo biên độ.
Ba ô phổ có bề rộng khác nhau
Giả định phổ sóng chỉ có năng lượng trong ba dải: 0,04–0,06 Hz; 0,08–0,12 Hz; và 0,18–0,22 Hz. Mật độ không đổi trong từng dải lần lượt bằng 5; 10; 2,5 m²/Hz. Bề rộng tương ứng 0,02; 0,04; 0,04 Hz, nên phương sai sóng của từng dải là 0,10; 0,40; 0,10 m².
Tổng phương sai m₀η=0,60 m², cho tham số phổ H_m0=4√m₀η=3,098387 m. Việc đặt mật độ bằng không trong các khoảng giữa dải là một phần định nghĩa ví dụ; không nội suy nối các điểm giữa dải rồi gọi đó là cùng một phổ.
Giả định H_z không đổi trên từng dải và lần lượt bằng 1; 0,8−0,6i; 0,2−0,1i. Môđun bình phương là 1; 1; 0,05. Phương sai heave tại gốc là:
m₀z=0,10×1+0,40×1+0,10×0,05=0,505 m².
Độ lệch chuẩn là 0,710634 m. Nếu bỏ bề rộng dải và chỉ cộng ba mật độ phổ, kết quả còn đơn vị m²/Hz và không phải phương sai; căn của nó cũng không phải độ lệch chuẩn chuyển vị.
Tịnh tiến điểm tính phải giữ pha heave và pitch
Chọn trục x dọc tàu, y ngang và z lên trên theo hệ thuận tay phải; pitch dương quanh trục y. Một điểm P tại x=20 m, y=0 có chuyển vị đứng góc nhỏ z_P=z−20θ. AQWA mô tả chuyển RAO từ trọng tâm đến một điểm khác bằng quan hệ vật rắn và véc tơ vị trí của điểm ấy. [3]
Giả định RAO pitch θ, tính bằng rad/m, trên ba dải lần lượt là 0,010; −0,015+0,020i; 0,005−0,003i. Tất cả dùng cùng gốc pha với H_z. RAO của điểm P phải được cộng dưới dạng phức trước khi bình phương:
H_P=H_z−20H_θ.
Ba giá trị thu được là 0,8; 1,1−1,0i; 0,1−0,04i. Môđun bình phương tương ứng 0,64; 2,21; 0,0116. Phương sai tại P bằng 0,064+0,884+0,00116=0,94916 m²; độ lệch chuẩn là 0,974248 m.
Chỉ riêng dải giữa đã đóng góp 0,884 m², lớn hơn cả phương sai heave tại gốc. Lý do không phải đã tạo thêm năng lượng sóng: chuyển động quay làm điểm cách gốc 20 m dịch đứng nhiều hơn, và quan hệ pha trong dải này làm hai thành phần tăng cường nhau tại P.
Kiểm bằng hiệp phương sai
Phương sai pitch trong ví dụ là 0,0002634 rad². Hiệp phương sai heave–pitch tính bằng tổng phần thực của H_z nhân liên hợp H_θ, với trọng số phương sai sóng từng dải, bằng −0,00847 m·rad. Vì vậy:
Var(z−xθ)=Var(z)+x²Var(θ)−2xCov(z,θ).
Tại x=20 m, hai phần phương sai riêng cộng thành 0,61036 m², còn số hạng chéo thêm 0,3388 m², đưa tổng về đúng 0,94916. Bỏ số hạng chéo tương đương giả định heave và pitch không tương quan; giả định ấy sai ở đây vì cả hai cùng do một trường sóng kích động.
Nếu dùng điểm đối xứng x=−20 m, dấu của số hạng chéo đổi. Phương sai chỉ còn 0,27156 m², độ lệch chuẩn 0,521114 m. Hai điểm có cùng khoảng cách đến gốc vẫn có thể có mức chuyển động khác nhau do dấu vị trí và pha của các bậc tự do.
Độ nhạy với lỗi pha và đơn vị góc
Nếu chỉ lưu biên độ heave và pitch rồi giả sử chúng luôn cùng pha hoặc luôn ngược pha, có thể tạo một kết quả tăng hoặc giảm giả. Chuyển toàn bộ RAO sang một gốc pha mới bằng cùng một hệ số phức có môđun bằng một không đổi phương sai của một điểm. Nhưng đổi pha chỉ cho một bậc tự do sẽ thay số hạng chéo và làm thay đổi bài toán.
RAO pitch trong ví dụ đã là rad/m. Nếu số 0,020 thực ra là độ/m mà bị dùng như rad/m, phần chuyển vị quay bị phóng đại 180/π lần; phương sai riêng do quay bị phóng đại bình phương tỷ số đó. Phép đổi phải diễn ra trước bước cộng với heave, không chỉ sửa nhãn đơn vị ở kết quả cuối.
Tích phân theo hướng là cộng phương sai của thành phần độc lập
Xét một dải tần riêng có phương sai sóng 0,20 m², chia 60% và 40% cho hai hướng độc lập. Giả định môđun RAO của một đáp ứng vô hướng là 1 và 0,3 theo hai hướng. Phương sai đáp ứng là 0,12×1²+0,08×0,3²=0,1272 m².
Cách lấy biên độ trung bình 0,6×1+0,4×0,3 rồi bình phương và nhân 0,20 cho 0,10368 m², không đúng với các hướng độc lập đã định nghĩa. Nếu hai thành phần có quan hệ pha kết hợp xác định, phải xét phổ chéo phù hợp; không được vừa gọi chúng độc lập vừa cộng biên độ như hai tín hiệu đồng bộ.
Đổi từ Hz sang radian mỗi giây
Đặt ω=2πf. Để bảo toàn phương sai, S_η,ω(ω)=S_η,f(f)/(2π), đồng thời dω=2πdf. Dải đầu có mật độ 5/(2π) m²/(rad/s) và bề rộng 2π×0,02 rad/s, nên tích vẫn bằng 0,10 m². Đổi nhãn trục tần số mà không đổi mật độ sẽ làm sai tích phân một hệ số 2π.
Tương tự, khi nội suy bảng RAO cần giữ đúng trục tần số, hướng sóng và trạng thái tải. Ba dải hằng trong bài cho phép tính chính xác bằng tổng; dữ liệu thực có đỉnh hẹp cần lưới đủ mịn để tích phân hội tụ. Tăng số điểm ở phần phổ gần bằng không không bù được việc bỏ qua một đỉnh đáp ứng hẹp.
Độ lệch chuẩn chưa phải cực đại trong ca làm việc
Với P, đại lượng 4√m₀P=3,896994 m là một độ cao chuyển động theo quy ước phổ. Không được gọi đó là biên độ một phía hoặc cực đại chắc chắn. Orcina phân biệt biên độ có nghĩa, cực đại theo thời lượng và chu kỳ đáp ứng, đồng thời nêu các giả thiết phổ hẹp, Gaussian và tuyến tính khi dùng công thức thống kê của báo cáo. [4]
Để đánh giá thao tác cẩu hoặc đi lại trên boong, cần tiếp tục tính đại lượng liên quan như gia tốc, chuyển động tương đối hay vượt ngưỡng, cùng thời lượng và điều kiện áp dụng. Bảng kết quả nên giữ điểm tính, quy ước pha, trạng thái tải, miền tích phân và phần năng lượng bị cắt. Nhờ đó người đọc có thể kiểm vì sao một con số thay đổi khi đổi vị trí hoặc hướng tàu.
Nguồn tham khảo
[1] Vessel theory: RAOs and phases — Orcina. Truy cập ngày 09/11/2026.
[2] Calculate response spectrum — waveresponse documentation — 4Subsea. Truy cập ngày 09/11/2026.
[3] AQWA Theory Manual 2024 R2 — section 4.4 Response Amplitude Operators — ANSYS. Truy cập ngày 09/11/2026.
[4] Vessel response reports — spectral response — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một biển có cùng chiều cao sóng có nghĩa có thể tạo chuyển động boong khác nhau khi năng lượng tập trung ở các tần số khác nhau. RAO mô tả đáp ứng theo một tần số và hướng sóng; nó gồm biên độ và pha. Orcina nhấn mạnh phải biết quy ước trục, hướng sóng, gốc pha và đơn vị của RAO trước khi ghép dữ liệu. [1]
Bài này tính một ví dụ giả định với ba dải tần, rồi chuyển từ heave tại gốc sang chuyển động đứng của điểm lệch khỏi gốc. Các con số không đại diện cho một tàu thật. Mô hình dùng hệ tuyến tính, trạng thái biển dừng và các thành phần tần số khác nhau không tương quan.
Vì sao phải bình phương môđun?
Nếu một thành phần sóng có biên độ a và hàm truyền phức H, biên độ đáp ứng là |H|a. Phương sai của một thành phần sin có biên độ a là a²/2; do đó phương sai đáp ứng nhân với |H|². Tài liệu waveresponse của 4Subsea viết phổ đầu ra dưới dạng tích H nhân liên hợp phức của H và phổ đầu vào, rồi tích phân theo hướng nếu có phổ định hướng. [2]
S_z(f)=|H_z(f)|²S_η(f).
Ở đây f tính bằng Hz, S_η có đơn vị m²/Hz và H_z có đơn vị m/m. Không dùng H² theo nghĩa bình phương một số phức, vì kết quả có thể phức và không phải mật độ phương sai. Cũng không chỉ nhân phổ với |H|, vì phổ là đại lượng bậc hai theo biên độ.
Ba ô phổ có bề rộng khác nhau
Giả định phổ sóng chỉ có năng lượng trong ba dải: 0,04–0,06 Hz; 0,08–0,12 Hz; và 0,18–0,22 Hz. Mật độ không đổi trong từng dải lần lượt bằng 5; 10; 2,5 m²/Hz. Bề rộng tương ứng 0,02; 0,04; 0,04 Hz, nên phương sai sóng của từng dải là 0,10; 0,40; 0,10 m².
Tổng phương sai m₀η=0,60 m², cho tham số phổ H_m0=4√m₀η=3,098387 m. Việc đặt mật độ bằng không trong các khoảng giữa dải là một phần định nghĩa ví dụ; không nội suy nối các điểm giữa dải rồi gọi đó là cùng một phổ.
Giả định H_z không đổi trên từng dải và lần lượt bằng 1; 0,8−0,6i; 0,2−0,1i. Môđun bình phương là 1; 1; 0,05. Phương sai heave tại gốc là:
m₀z=0,10×1+0,40×1+0,10×0,05=0,505 m².
Độ lệch chuẩn là 0,710634 m. Nếu bỏ bề rộng dải và chỉ cộng ba mật độ phổ, kết quả còn đơn vị m²/Hz và không phải phương sai; căn của nó cũng không phải độ lệch chuẩn chuyển vị.
Tịnh tiến điểm tính phải giữ pha heave và pitch
Chọn trục x dọc tàu, y ngang và z lên trên theo hệ thuận tay phải; pitch dương quanh trục y. Một điểm P tại x=20 m, y=0 có chuyển vị đứng góc nhỏ z_P=z−20θ. AQWA mô tả chuyển RAO từ trọng tâm đến một điểm khác bằng quan hệ vật rắn và véc tơ vị trí của điểm ấy. [3]
Giả định RAO pitch θ, tính bằng rad/m, trên ba dải lần lượt là 0,010; −0,015+0,020i; 0,005−0,003i. Tất cả dùng cùng gốc pha với H_z. RAO của điểm P phải được cộng dưới dạng phức trước khi bình phương:
H_P=H_z−20H_θ.
Ba giá trị thu được là 0,8; 1,1−1,0i; 0,1−0,04i. Môđun bình phương tương ứng 0,64; 2,21; 0,0116. Phương sai tại P bằng 0,064+0,884+0,00116=0,94916 m²; độ lệch chuẩn là 0,974248 m.
Chỉ riêng dải giữa đã đóng góp 0,884 m², lớn hơn cả phương sai heave tại gốc. Lý do không phải đã tạo thêm năng lượng sóng: chuyển động quay làm điểm cách gốc 20 m dịch đứng nhiều hơn, và quan hệ pha trong dải này làm hai thành phần tăng cường nhau tại P.
Kiểm bằng hiệp phương sai
Phương sai pitch trong ví dụ là 0,0002634 rad². Hiệp phương sai heave–pitch tính bằng tổng phần thực của H_z nhân liên hợp H_θ, với trọng số phương sai sóng từng dải, bằng −0,00847 m·rad. Vì vậy:
Var(z−xθ)=Var(z)+x²Var(θ)−2xCov(z,θ).
Tại x=20 m, hai phần phương sai riêng cộng thành 0,61036 m², còn số hạng chéo thêm 0,3388 m², đưa tổng về đúng 0,94916. Bỏ số hạng chéo tương đương giả định heave và pitch không tương quan; giả định ấy sai ở đây vì cả hai cùng do một trường sóng kích động.
Nếu dùng điểm đối xứng x=−20 m, dấu của số hạng chéo đổi. Phương sai chỉ còn 0,27156 m², độ lệch chuẩn 0,521114 m. Hai điểm có cùng khoảng cách đến gốc vẫn có thể có mức chuyển động khác nhau do dấu vị trí và pha của các bậc tự do.
Độ nhạy với lỗi pha và đơn vị góc
Nếu chỉ lưu biên độ heave và pitch rồi giả sử chúng luôn cùng pha hoặc luôn ngược pha, có thể tạo một kết quả tăng hoặc giảm giả. Chuyển toàn bộ RAO sang một gốc pha mới bằng cùng một hệ số phức có môđun bằng một không đổi phương sai của một điểm. Nhưng đổi pha chỉ cho một bậc tự do sẽ thay số hạng chéo và làm thay đổi bài toán.
RAO pitch trong ví dụ đã là rad/m. Nếu số 0,020 thực ra là độ/m mà bị dùng như rad/m, phần chuyển vị quay bị phóng đại 180/π lần; phương sai riêng do quay bị phóng đại bình phương tỷ số đó. Phép đổi phải diễn ra trước bước cộng với heave, không chỉ sửa nhãn đơn vị ở kết quả cuối.
Tích phân theo hướng là cộng phương sai của thành phần độc lập
Xét một dải tần riêng có phương sai sóng 0,20 m², chia 60% và 40% cho hai hướng độc lập. Giả định môđun RAO của một đáp ứng vô hướng là 1 và 0,3 theo hai hướng. Phương sai đáp ứng là 0,12×1²+0,08×0,3²=0,1272 m².
Cách lấy biên độ trung bình 0,6×1+0,4×0,3 rồi bình phương và nhân 0,20 cho 0,10368 m², không đúng với các hướng độc lập đã định nghĩa. Nếu hai thành phần có quan hệ pha kết hợp xác định, phải xét phổ chéo phù hợp; không được vừa gọi chúng độc lập vừa cộng biên độ như hai tín hiệu đồng bộ.
Đổi từ Hz sang radian mỗi giây
Đặt ω=2πf. Để bảo toàn phương sai, S_η,ω(ω)=S_η,f(f)/(2π), đồng thời dω=2πdf. Dải đầu có mật độ 5/(2π) m²/(rad/s) và bề rộng 2π×0,02 rad/s, nên tích vẫn bằng 0,10 m². Đổi nhãn trục tần số mà không đổi mật độ sẽ làm sai tích phân một hệ số 2π.
Tương tự, khi nội suy bảng RAO cần giữ đúng trục tần số, hướng sóng và trạng thái tải. Ba dải hằng trong bài cho phép tính chính xác bằng tổng; dữ liệu thực có đỉnh hẹp cần lưới đủ mịn để tích phân hội tụ. Tăng số điểm ở phần phổ gần bằng không không bù được việc bỏ qua một đỉnh đáp ứng hẹp.
Độ lệch chuẩn chưa phải cực đại trong ca làm việc
Với P, đại lượng 4√m₀P=3,896994 m là một độ cao chuyển động theo quy ước phổ. Không được gọi đó là biên độ một phía hoặc cực đại chắc chắn. Orcina phân biệt biên độ có nghĩa, cực đại theo thời lượng và chu kỳ đáp ứng, đồng thời nêu các giả thiết phổ hẹp, Gaussian và tuyến tính khi dùng công thức thống kê của báo cáo. [4]
Để đánh giá thao tác cẩu hoặc đi lại trên boong, cần tiếp tục tính đại lượng liên quan như gia tốc, chuyển động tương đối hay vượt ngưỡng, cùng thời lượng và điều kiện áp dụng. Bảng kết quả nên giữ điểm tính, quy ước pha, trạng thái tải, miền tích phân và phần năng lượng bị cắt. Nhờ đó người đọc có thể kiểm vì sao một con số thay đổi khi đổi vị trí hoặc hướng tàu.
Nguồn tham khảo
[1] Vessel theory: RAOs and phases — Orcina. Truy cập ngày 09/11/2026.
[2] Calculate response spectrum — waveresponse documentation — 4Subsea. Truy cập ngày 09/11/2026.
[3] AQWA Theory Manual 2024 R2 — section 4.4 Response Amplitude Operators — ANSYS. Truy cập ngày 09/11/2026.
[4] Vessel response reports — spectral response — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
