Khối lượng thêm không phải một số luôn dùng được ở mọi tần số
Trong mô hình vật nổi, nước tạo phản lực khi thân gia tốc và khi thân phát sóng ra môi trường. Tài liệu WEC-Sim biểu diễn phần này bằng khối lượng thêm và cản bức xạ phụ thuộc tần số; khi tính trong miền thời gian, có thể dùng lực phụ thuộc lịch sử vận tốc. [1] Vì vậy một con số khối lượng thêm được lấy tại chu kỳ nào phải đi kèm chính chu kỳ đó.
Bài viết dùng một hệ một bậc tự do và một mô hình nhớ toán học do tác giả đặt ra. Tất cả hệ số là giả định; chúng không được trích từ thử nghiệm hay mô hình nhiễu xạ của một thân tàu. Mục tiêu là kiểm tính nhất quán giữa các cách tính, không cung cấp hệ số để thiết kế một vật nổi cụ thể.
Phương trình điều hòa và ba ý nghĩa của chu kỳ
Với chuyển vị z và quy ước điều hòa e^(iωt), đặt m là khối lượng thân, A(ω) là khối lượng thêm, C là độ cứng và B(ω)+b₀ là tổng cản. Độ cứng động phức có phần thực và phần ảo:
D_R=C−ω²[m+A(ω)].
D_I=ω[B(ω)+b₀].
Biên độ dưới lực điều hòa có biên độ F₀ là F₀/√(D_R²+D_I²). Tần số làm D_R bằng không là một dấu mốc cân bằng phần quán tính và độ cứng. Đỉnh biên độ còn phụ thuộc D_I thay đổi theo tần số. Dao động tự do lại do các nghiệm đặc trưng trong miền thời gian quyết định. Ba kết quả có thể gần nhau nhưng không được định nghĩa là một.
Orcina lưu ý rằng cách dùng khối lượng thêm và cản hằng phải chọn hệ số phù hợp miền tần số chuyển động, còn miền rộng cần dữ liệu phụ thuộc tần số. [2] Việc giải ω²[m+A(ω)]=C là một bài toán ẩn, không phải chỉ thay A tại tần số bất kỳ vào công thức căn.
Dựng một mô hình có nhớ để tránh ghép hệ số tùy ý
Giả định m=1.000.000 kg, A_∞=1.000.000 kg, C=500.000 N/m và b₀=20.000 N·s/m. Chọn hai tham số của hàm nhớ là k_r=150.000 N/m và λ=0,5 s⁻¹. Hàm nhớ theo độ trễ t≥0 được đặt:
K_r(t)=k_r(1−λt)e^(−λt).
Lực nhớ trong vế trái phương trình là tích phân của K_r(t−τ) nhân vận tốc tại τ, từ 0 đến t; phần quán tính tức thời là (m+A_∞)z̈. Mô hình giả định vật không chuyển động trước t=0, còn điều kiện chuyển vị ban đầu được khai báo riêng.
Biến đổi Laplace của hàm nhớ là K̂_r(s)=k_rs/(s+λ)². Vì lực nhớ nhân vận tốc và nằm ở vế trái, đóng góp vào độ cứng động là sK̂_r(s). Với đầu vào điều hòa ổn định:
D(iω)=C−(m+A_∞)ω²+iωb₀+iωK̂_r(iω).
Tách phần thực và ảo, rồi so với D_R và D_I ở trên, thu được cặp hệ số:
A(ω)=A_∞+k_r(λ²−ω²)/(λ²+ω²)².
B(ω)=2k_rλω²/(λ²+ω²)².
Đây là hai hệ số suy ra từ cùng một hàm nhớ, không chọn hai bảng độc lập cho tiện khớp kết quả. B luôn không âm và tiến về không ở hai giới hạn tần số; mô hình là một minh họa nhân quả có tiêu tán. Các tính chất này chưa chứng minh nó là đáp ứng của một hình dạng thân thực tế.
Orcina mô tả quan hệ giữa khối lượng thêm, cản và đáp ứng xung, đồng thời cảnh báo dữ liệu không nhất quán có thể cho kết quả không đáng tin cậy. [3] Phần mô hình giả định ở đây được xây dựng để người đọc tự kiểm quan hệ ấy, không thay cho dữ liệu thủy động lực đã xác minh.
Tính các hệ số tại ba tần số
Tại ω=0,2 rad/s, A=1.374.554,102259 kg và B=71.343,638526 N·s/m. Tại ω=0,5 rad/s, A=1.000.000 kg và B=150.000 N·s/m. Tại ω=1 rad/s, A=928.000 kg và B=96.000 N·s/m. Ở tần số thấp A lớn hơn A_∞; sau khi ω vượt λ, A thấp hơn A_∞ rồi tiến về giới hạn ấy từ phía dưới. Không nên áp đặt một đường nội suy đơn điệu chỉ vì tên gọi có chữ khối lượng.
Ở tần số bằng không, A(0)=1.600.000 kg. Nếu dùng số này như hằng số cho mọi chuyển động, công thức 2π√[(m+A)/C] cho 14,327867 s. Nếu dùng A_∞ thì được 12,566371 s. Hai phép thế số đều thực hiện đúng số học nhưng chưa giải quyết câu hỏi hệ số nào phù hợp chuyển động.
Giải điều kiện cân bằng phần thực
Đặt y=ω². Sau khi thay hàm A vào D_R=0 và nhân với (λ²+y)², phương trình của ví dụ được phân tích thành:
(0,25−y)(2.000.000y²+850.000y+125.000)=0.
Biến y có đơn vị s⁻²; các hệ số số học trong biểu thức đã mang các đơn vị tương ứng của bộ dữ liệu SI. Thừa số bậc hai luôn dương khi y≥0, nên nghiệm vật lý duy nhất trên miền này là y=0,25 s⁻². Do đó ω_R=0,5 rad/s và T_R=12,566371 s.
Tại đúng tần số này A bằng A_∞, nên T_R tình cờ bằng phép tính dùng A_∞. Điều đó do bộ số giả định đã chọn, không chứng minh có thể luôn dùng giới hạn tần số vô hạn để tìm chu kỳ.
Đỉnh đáp ứng dưới lực điều hòa
Giữ biên độ lực F₀=100.000 N không đổi theo tần số. Tại ω_R, biên độ chuyển vị là 1,176471 m. Tối thiểu hóa D_R²+D_I² trên miền tần số dương cho đỉnh tại ω≈0,495696 rad/s, chu kỳ 12,675490 s, biên độ 1,181545 m. Đỉnh hơi lệch vì phần cản động cũng thay đổi theo ω.
Đây là đáp ứng theo lực, chưa phải RAO theo biên độ sóng. Với lực kích động sóng phụ thuộc tần số, cả tử số cũng thay đổi; vị trí đỉnh có thể dịch thêm. Biên độ hơn một mét chỉ là nghiệm của mô hình tuyến tính giả định và không tự bảo đảm thân thực còn nằm trong miền tuyến tính.
Dao động tự do cần giữ trạng thái nhớ
Đặt lực ngoài bằng không và tìm nghiệm dạng e^(st). Sau khi nhân phương trình với (s+λ)², đa thức đặc trưng là:
2.000.000s⁴+2.020.000s³+1.170.000s²+505.000s+125.000=0.
Các nghiệm là khoảng −0,046329561±0,497848945i s⁻¹ và −0,458670439±0,199051321i s⁻¹. Cặp đầu tắt chậm hơn, có chu kỳ thành phần dao động 12,620666 s; cặp thứ hai có chu kỳ 31,565655 s nhưng phần thực âm lớn nên suy giảm nhanh. Điều kiện ban đầu quyết định biên độ của mỗi thành phần.
Do đó bản ghi tự do ban đầu không nhất thiết là một đường sin có một biên bao duy nhất. Chỉ sau khi thành phần tắt nhanh nhỏ đi, phép khớp một cặp dao động có thể trở nên hợp lý. Chu kỳ 12,620666 s của thành phần tắt chậm khác cả chu kỳ tại D_R=0 và chu kỳ ở đỉnh cưỡng bức.
Kiểm dấu và giới hạn của mô hình nhớ
Tại t=4 s, K_r(t)=−20.300,292485 N/m. Một đoạn âm của hàm nhớ không có nghĩa hệ đang tạo năng lượng từ không khí: dấu của hàm theo thời gian phải được xét cùng tích phân và quan hệ tần số. Trong mô hình này B(ω) không âm; bỏ đoạn âm của K_r sẽ thay đổi cả A lẫn B, chứ không chỉ làm phép tính “an toàn hơn”.
Cũng không cộng thêm A(ω) vào m+A_∞ khi đã dùng đầy đủ lực nhớ, vì đó là tính trùng phản ứng của nước. Khi chuyển sang thân thật, cần giữ trục tham chiếu, mớn, độ sâu, miền tần số và phép chuẩn hóa hệ số. Một phép kiểm hữu ích là tái tạo đáp ứng điều hòa từ mô hình thời gian rồi đối chiếu với bảng A, B đã nhập. Sự phù hợp của ba cách biểu diễn giúp phát hiện lỗi triển khai; độ đúng vật lý vẫn cần dữ liệu thủy động lực hoặc thử nghiệm phù hợp hình học và trạng thái đang đánh giá.
Nguồn tham khảo
[1] WEC-Sim Theory — sinusoidal response and convolution integral formulation — WEC-Sim Documentation / NREL and Sandia. Truy cập ngày 09/11/2026.
[2] Vessel types: Stiffness, added mass and damping — Orcina. Truy cập ngày 09/11/2026.
[3] Vessel theory: Impulse response and convolution — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong mô hình vật nổi, nước tạo phản lực khi thân gia tốc và khi thân phát sóng ra môi trường. Tài liệu WEC-Sim biểu diễn phần này bằng khối lượng thêm và cản bức xạ phụ thuộc tần số; khi tính trong miền thời gian, có thể dùng lực phụ thuộc lịch sử vận tốc. [1] Vì vậy một con số khối lượng thêm được lấy tại chu kỳ nào phải đi kèm chính chu kỳ đó.
Bài viết dùng một hệ một bậc tự do và một mô hình nhớ toán học do tác giả đặt ra. Tất cả hệ số là giả định; chúng không được trích từ thử nghiệm hay mô hình nhiễu xạ của một thân tàu. Mục tiêu là kiểm tính nhất quán giữa các cách tính, không cung cấp hệ số để thiết kế một vật nổi cụ thể.
Phương trình điều hòa và ba ý nghĩa của chu kỳ
Với chuyển vị z và quy ước điều hòa e^(iωt), đặt m là khối lượng thân, A(ω) là khối lượng thêm, C là độ cứng và B(ω)+b₀ là tổng cản. Độ cứng động phức có phần thực và phần ảo:
D_R=C−ω²[m+A(ω)].
D_I=ω[B(ω)+b₀].
Biên độ dưới lực điều hòa có biên độ F₀ là F₀/√(D_R²+D_I²). Tần số làm D_R bằng không là một dấu mốc cân bằng phần quán tính và độ cứng. Đỉnh biên độ còn phụ thuộc D_I thay đổi theo tần số. Dao động tự do lại do các nghiệm đặc trưng trong miền thời gian quyết định. Ba kết quả có thể gần nhau nhưng không được định nghĩa là một.
Orcina lưu ý rằng cách dùng khối lượng thêm và cản hằng phải chọn hệ số phù hợp miền tần số chuyển động, còn miền rộng cần dữ liệu phụ thuộc tần số. [2] Việc giải ω²[m+A(ω)]=C là một bài toán ẩn, không phải chỉ thay A tại tần số bất kỳ vào công thức căn.
Dựng một mô hình có nhớ để tránh ghép hệ số tùy ý
Giả định m=1.000.000 kg, A_∞=1.000.000 kg, C=500.000 N/m và b₀=20.000 N·s/m. Chọn hai tham số của hàm nhớ là k_r=150.000 N/m và λ=0,5 s⁻¹. Hàm nhớ theo độ trễ t≥0 được đặt:
K_r(t)=k_r(1−λt)e^(−λt).
Lực nhớ trong vế trái phương trình là tích phân của K_r(t−τ) nhân vận tốc tại τ, từ 0 đến t; phần quán tính tức thời là (m+A_∞)z̈. Mô hình giả định vật không chuyển động trước t=0, còn điều kiện chuyển vị ban đầu được khai báo riêng.
Biến đổi Laplace của hàm nhớ là K̂_r(s)=k_rs/(s+λ)². Vì lực nhớ nhân vận tốc và nằm ở vế trái, đóng góp vào độ cứng động là sK̂_r(s). Với đầu vào điều hòa ổn định:
D(iω)=C−(m+A_∞)ω²+iωb₀+iωK̂_r(iω).
Tách phần thực và ảo, rồi so với D_R và D_I ở trên, thu được cặp hệ số:
A(ω)=A_∞+k_r(λ²−ω²)/(λ²+ω²)².
B(ω)=2k_rλω²/(λ²+ω²)².
Đây là hai hệ số suy ra từ cùng một hàm nhớ, không chọn hai bảng độc lập cho tiện khớp kết quả. B luôn không âm và tiến về không ở hai giới hạn tần số; mô hình là một minh họa nhân quả có tiêu tán. Các tính chất này chưa chứng minh nó là đáp ứng của một hình dạng thân thực tế.
Orcina mô tả quan hệ giữa khối lượng thêm, cản và đáp ứng xung, đồng thời cảnh báo dữ liệu không nhất quán có thể cho kết quả không đáng tin cậy. [3] Phần mô hình giả định ở đây được xây dựng để người đọc tự kiểm quan hệ ấy, không thay cho dữ liệu thủy động lực đã xác minh.
Tính các hệ số tại ba tần số
Tại ω=0,2 rad/s, A=1.374.554,102259 kg và B=71.343,638526 N·s/m. Tại ω=0,5 rad/s, A=1.000.000 kg và B=150.000 N·s/m. Tại ω=1 rad/s, A=928.000 kg và B=96.000 N·s/m. Ở tần số thấp A lớn hơn A_∞; sau khi ω vượt λ, A thấp hơn A_∞ rồi tiến về giới hạn ấy từ phía dưới. Không nên áp đặt một đường nội suy đơn điệu chỉ vì tên gọi có chữ khối lượng.
Ở tần số bằng không, A(0)=1.600.000 kg. Nếu dùng số này như hằng số cho mọi chuyển động, công thức 2π√[(m+A)/C] cho 14,327867 s. Nếu dùng A_∞ thì được 12,566371 s. Hai phép thế số đều thực hiện đúng số học nhưng chưa giải quyết câu hỏi hệ số nào phù hợp chuyển động.
Giải điều kiện cân bằng phần thực
Đặt y=ω². Sau khi thay hàm A vào D_R=0 và nhân với (λ²+y)², phương trình của ví dụ được phân tích thành:
(0,25−y)(2.000.000y²+850.000y+125.000)=0.
Biến y có đơn vị s⁻²; các hệ số số học trong biểu thức đã mang các đơn vị tương ứng của bộ dữ liệu SI. Thừa số bậc hai luôn dương khi y≥0, nên nghiệm vật lý duy nhất trên miền này là y=0,25 s⁻². Do đó ω_R=0,5 rad/s và T_R=12,566371 s.
Tại đúng tần số này A bằng A_∞, nên T_R tình cờ bằng phép tính dùng A_∞. Điều đó do bộ số giả định đã chọn, không chứng minh có thể luôn dùng giới hạn tần số vô hạn để tìm chu kỳ.
Đỉnh đáp ứng dưới lực điều hòa
Giữ biên độ lực F₀=100.000 N không đổi theo tần số. Tại ω_R, biên độ chuyển vị là 1,176471 m. Tối thiểu hóa D_R²+D_I² trên miền tần số dương cho đỉnh tại ω≈0,495696 rad/s, chu kỳ 12,675490 s, biên độ 1,181545 m. Đỉnh hơi lệch vì phần cản động cũng thay đổi theo ω.
Đây là đáp ứng theo lực, chưa phải RAO theo biên độ sóng. Với lực kích động sóng phụ thuộc tần số, cả tử số cũng thay đổi; vị trí đỉnh có thể dịch thêm. Biên độ hơn một mét chỉ là nghiệm của mô hình tuyến tính giả định và không tự bảo đảm thân thực còn nằm trong miền tuyến tính.
Dao động tự do cần giữ trạng thái nhớ
Đặt lực ngoài bằng không và tìm nghiệm dạng e^(st). Sau khi nhân phương trình với (s+λ)², đa thức đặc trưng là:
2.000.000s⁴+2.020.000s³+1.170.000s²+505.000s+125.000=0.
Các nghiệm là khoảng −0,046329561±0,497848945i s⁻¹ và −0,458670439±0,199051321i s⁻¹. Cặp đầu tắt chậm hơn, có chu kỳ thành phần dao động 12,620666 s; cặp thứ hai có chu kỳ 31,565655 s nhưng phần thực âm lớn nên suy giảm nhanh. Điều kiện ban đầu quyết định biên độ của mỗi thành phần.
Do đó bản ghi tự do ban đầu không nhất thiết là một đường sin có một biên bao duy nhất. Chỉ sau khi thành phần tắt nhanh nhỏ đi, phép khớp một cặp dao động có thể trở nên hợp lý. Chu kỳ 12,620666 s của thành phần tắt chậm khác cả chu kỳ tại D_R=0 và chu kỳ ở đỉnh cưỡng bức.
Kiểm dấu và giới hạn của mô hình nhớ
Tại t=4 s, K_r(t)=−20.300,292485 N/m. Một đoạn âm của hàm nhớ không có nghĩa hệ đang tạo năng lượng từ không khí: dấu của hàm theo thời gian phải được xét cùng tích phân và quan hệ tần số. Trong mô hình này B(ω) không âm; bỏ đoạn âm của K_r sẽ thay đổi cả A lẫn B, chứ không chỉ làm phép tính “an toàn hơn”.
Cũng không cộng thêm A(ω) vào m+A_∞ khi đã dùng đầy đủ lực nhớ, vì đó là tính trùng phản ứng của nước. Khi chuyển sang thân thật, cần giữ trục tham chiếu, mớn, độ sâu, miền tần số và phép chuẩn hóa hệ số. Một phép kiểm hữu ích là tái tạo đáp ứng điều hòa từ mô hình thời gian rồi đối chiếu với bảng A, B đã nhập. Sự phù hợp của ba cách biểu diễn giúp phát hiện lỗi triển khai; độ đúng vật lý vẫn cần dữ liệu thủy động lực hoặc thử nghiệm phù hợp hình học và trạng thái đang đánh giá.
Nguồn tham khảo
[1] WEC-Sim Theory — sinusoidal response and convolution integral formulation — WEC-Sim Documentation / NREL and Sandia. Truy cập ngày 09/11/2026.
[2] Vessel types: Stiffness, added mass and damping — Orcina. Truy cập ngày 09/11/2026.
[3] Vessel theory: Impulse response and convolution — Orcina. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
