1. Nhiều dây cùng giữ một vật nổi
Tải ngang trên một hệ neo không được chia đều chỉ theo số dây. Hướng dây, độ dốc đặc tính lực–khoảng cách và lực căng trước đều ảnh hưởng phân phối. Tài liệu Ansys Aqwa xây dựng độ cứng từ đạo hàm tải–chuyển vị và chuyển các đóng góp về cùng hệ trục trước khi tìm cân bằng. [1]
Bài này xét một điểm nối chung chỉ được chuyển ngang x, y; cao độ và các góc quay bị giữ cố định. Đây là mô hình hai bậc tự do để kiểm phép lắp ma trận. Một vật nổi thật có nhiều vị trí fairlead, chuyển động đứng và quay sẽ cần mô hình đầy đủ hơn. MoorPy phân biệt độ cứng của toàn bộ bậc tự do với độ cứng sau khi cho các bậc tự do khác tự cân bằng. [2]
2. Hai nguồn đóng góp vào độ cứng
Một neo nằm tại aᵢ; điểm nối tại r. Khoảng ngang Rᵢ = |aᵢ − r|. Giả sử độ lớn lực ngang Hᵢ chỉ phụ thuộc Rᵢ ở cao độ đang giữ. Véc tơ lực lên điểm nối hướng về neo:
Fᵢ = Hᵢ(Rᵢ)(aᵢ − r)/Rᵢ
Một dịch chuyển nhỏ vừa đổi Hᵢ, vừa đổi hướng véc tơ. Đặt kr = dH/dR là độ cứng dọc hướng hình chiếu ngang, kg = H/R là phần độ cứng hình học do đổi hướng. Nếu θ là góc véc tơ nối tới neo, đo ngược chiều kim đồng hồ từ +x, đóng góp của dây là:
Kxx = kr cos²θ + kg sin²θ
Kyy = kr sin²θ + kg cos²θ
Kxy = Kyx = (kr − kg)sinθ cosθ
Các hệ số tính bằng kN/m khi H dùng kN và R dùng m. Đây là đạo hàm của véc tơ lực vừa nêu, với quy ước ΔF = −KΔr. Thành phần kg vẫn tồn tại dù độ lớn lực gần không đổi, vì một dây đang căng bị xoay sẽ đổi hướng lực. Nếu chỉ chiếu kr rồi bỏ kg, ta đã bỏ một phần cơ học của chính mô hình hướng tâm này.
3. Hệ ba dây giả định
Chọn ba neo cách gốc 100 m, ở các góc 0°, 120°, 240°. Tọa độ tương ứng là (100; 0), (−50; 86,60254), (−50; −86,60254) m. Ở r = 0, mỗi dây có H₀ = 200 kN. Ba véc tơ lực bằng nhau và cách nhau 120° nên cân bằng khi chưa có tải ngoài.
Các độ cứng kr giả định lần lượt 20, 30, 10 kN/m. Chúng đại diện ba đặc tính khác nhau quanh cùng lực căng trước; không phải số liệu một hệ neo thực. Với R₀ = 100 m, tất cả kg = 2 kN/m. Sau khi cộng ba ma trận, ta nhận:
Kxx = Kyy = 33 kN/m
Kxy = Kyx = −8,660254 kN/m
Định thức bằng 1.014 (kN/m)², dương. Hai trị riêng là 24,339746 và 41,660254 kN/m, nên mô hình có độ cứng phục hồi dương theo mọi hướng ngang gần gốc. Điều này chỉ kiểm mô hình hai chiều ở trạng thái ban đầu; không kết luận về ổn định đứng, lắc, tải đỉnh hoặc trạng thái đứt dây của một công trình.
Hai phần tử chéo chính bằng nhau nhưng hệ vẫn ghép x–y. Dây ở 120° cứng hơn dây ở 240° làm hệ số ngoài đường chéo khác không. Vì vậy “độ cứng x bằng độ cứng y” không có nghĩa dịch chuyển sẽ luôn cùng hướng với tải.
4. Giải chuyển vị do tải ngoài
Đặt tải giả định Px = 120 kN, Py = 40 kN tại điểm nối. Mô hình tuyến tính cho:
33Δx − 8,660254Δy = 120
−8,660254Δx + 33Δy = 40
Nghiệm Δx = 4,246953 m và Δy = 2,326657 m. Nếu bỏ hệ số ghép và tính riêng 120/33, 40/33, sẽ nhận 3,636364 m và 1,212121 m, khác đáng kể. Không thể coi ma trận như hai lò xo độc lập khi các phần tử ngoài đường chéo chưa được chứng minh nhỏ.
Nếu chỉ bỏ độ cứng hình học kg nhưng giữ phần chiếu kr, ma trận có hai phần tử chéo bằng 30 và phần tử ghép vẫn −8,660254. Khi đó nghiệm là Δx = 4,783527 m, Δy = 2,714219 m. Đây là sai khác do bỏ một thành phần độ cứng, không phải do tuyến tính hóa khác thời điểm.
Phép kiểm đơn giản nhất là nhân lại K với chuyển vị để khôi phục tải. Tiếp đó kiểm các lực ban đầu có tổng bằng không. Nếu một cấu hình đã lệch cân bằng ngay từ đầu, chỉ giải KΔr = P mà bỏ lực dư sẽ dẫn tới một vị trí sai dù phép nhân ma trận hoàn toàn đúng.
5. Kiểm bằng mô hình hình học phi tuyến
Để kiểm ảnh hưởng đổi hướng, định nghĩa một mô hình giả định có đặc tính lực theo khoảng ngang tuyến tính:
Hᵢ = 200 + krᵢ(Rᵢ − 100)
Chỉ dùng trong vùng cả ba Hᵢ còn dương. Ở mỗi vị trí r, tính Rᵢ và hướng thực của dây, rồi giải:
P + ΣFᵢ(r) = 0
Mặc dù H theo R tuyến tính, véc tơ tổng lực vẫn phi tuyến theo x, y vì cả khoảng cách và hướng đều thay đổi. Đây chưa phải nghiệm catenary đầy đủ; nó là mô hình hướng tâm được định nghĩa rõ để kiểm phép tuyến tính hóa.
Với tải 120; 40 kN, nghiệm cân bằng là x = 4,225159 m, y = 2,472398 m. Các khoảng ngang là 95,806748; 100,091202; 104,281890 m. Lực ngang ba dây lần lượt 116,134963; 202,736053; 242,818896 kN. Cộng các véc tơ theo hướng mới cho đúng (−120; −40) kN, cân bằng tải ngoài.
Đưa nghiệm tuyến tính (4,246953; 2,326657) vào mô hình hình học này còn lực dư khoảng (−2,1811; +4,7233) kN. Đây là phép kiểm có ý nghĩa hơn chỉ thấy dịch chuyển “có vẻ nhỏ”. Tài liệu Aqwa mô tả lặp vị trí để giảm mất cân bằng, đồng thời lưu ý hệ phi tuyến có thể vượt quá nghiệm mong muốn nếu bước lặp không được kiểm soát. [3]
6. Độ nhạy khi tăng tải và mất dây
Nhân đôi tải thành (240; 80) kN. Tuyến tính hóa tại gốc cho x = 8,493906 m, y = 4,653315 m. Mô hình hình học phi tuyến cho x = 8,438758 m, y = 5,283954 m. Các lực tương ứng khoảng 34,271645; 204,167240; 288,954378 kN. Dây thứ nhất đã giảm lực đáng kể nhưng vẫn dương trong kịch bản này.
Nếu tải tiếp tục tăng làm một Hᵢ tính ra âm, dây không được biến thành thanh chịu nén. Phải chuyển sang nhánh chùng theo mô hình đã định nghĩa, giải lại cân bằng và cập nhật độ cứng. Không nên kéo dài đường tuyến tính ban đầu qua điểm mất căng rồi coi kết quả là lực neo hợp lệ.
Mất một dây cũng làm thay đổi cả lực ban đầu lẫn ma trận. Ví dụ bỏ dây góc 120° tại gốc khiến tổng lực hai dây còn lại bằng (100; −173,2051) kN trước khi cộng tải môi trường. Vì vậy tình huống mất dây không được tính chỉ bằng cách xóa một ma trận rồi dùng nguyên tải ngoài; phải đưa mất cân bằng do lực căng trước vào bài toán và tìm một trạng thái mới. Bài này chưa tính tải quá độ ngay sau sự cố.
7. Áp dụng vào mô hình vật nổi
Với fairlead đặt xa trọng tâm, lực dây còn gây mômen; một góc quay nhỏ dịch chuyển từng fairlead khác nhau. Aqwa lưu ý việc chuyển ma trận về tâm vật phải xét thêm thay đổi mômen khi vật quay. [1] Không được ghép các độ cứng tại nhiều điểm bằng phép cộng trực tiếp nếu chúng chưa cùng gốc và cùng bậc tự do.
Hồ sơ tính nên giữ vị trí neo/fairlead, lực căng trước, đặc tính H(R) hoặc mô hình dây, các tọa độ bị khóa, ma trận và lực dư sau kiểm lại. Chuyển vị ngang chỉ là một đầu ra; lực từng dây và điều kiện còn căng cần được kiểm cùng lúc. Như vậy kết quả giải hệ sẽ gắn với cấu hình cơ học cụ thể thay vì chỉ là một phép chia tổng tải cho “độ cứng neo”.
Nguồn tham khảo
[1] Aqwa Theory Manual 2024 R2 — The Stiffness Matrix — Ansys. Truy cập ngày 09/11/2026.
[2] MoorPy Usage — System Equilibrium and Stiffness Matrices — MoorPy developers, National Renewable Energy Laboratory. Truy cập ngày 09/11/2026.
[3] Aqwa Theory Manual 2024 R2 — Iteration Towards Equilibrium — Ansys. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Tải ngang trên một hệ neo không được chia đều chỉ theo số dây. Hướng dây, độ dốc đặc tính lực–khoảng cách và lực căng trước đều ảnh hưởng phân phối. Tài liệu Ansys Aqwa xây dựng độ cứng từ đạo hàm tải–chuyển vị và chuyển các đóng góp về cùng hệ trục trước khi tìm cân bằng. [1]
Bài này xét một điểm nối chung chỉ được chuyển ngang x, y; cao độ và các góc quay bị giữ cố định. Đây là mô hình hai bậc tự do để kiểm phép lắp ma trận. Một vật nổi thật có nhiều vị trí fairlead, chuyển động đứng và quay sẽ cần mô hình đầy đủ hơn. MoorPy phân biệt độ cứng của toàn bộ bậc tự do với độ cứng sau khi cho các bậc tự do khác tự cân bằng. [2]
2. Hai nguồn đóng góp vào độ cứng
Một neo nằm tại aᵢ; điểm nối tại r. Khoảng ngang Rᵢ = |aᵢ − r|. Giả sử độ lớn lực ngang Hᵢ chỉ phụ thuộc Rᵢ ở cao độ đang giữ. Véc tơ lực lên điểm nối hướng về neo:
Fᵢ = Hᵢ(Rᵢ)(aᵢ − r)/Rᵢ
Một dịch chuyển nhỏ vừa đổi Hᵢ, vừa đổi hướng véc tơ. Đặt kr = dH/dR là độ cứng dọc hướng hình chiếu ngang, kg = H/R là phần độ cứng hình học do đổi hướng. Nếu θ là góc véc tơ nối tới neo, đo ngược chiều kim đồng hồ từ +x, đóng góp của dây là:
Kxx = kr cos²θ + kg sin²θ
Kyy = kr sin²θ + kg cos²θ
Kxy = Kyx = (kr − kg)sinθ cosθ
Các hệ số tính bằng kN/m khi H dùng kN và R dùng m. Đây là đạo hàm của véc tơ lực vừa nêu, với quy ước ΔF = −KΔr. Thành phần kg vẫn tồn tại dù độ lớn lực gần không đổi, vì một dây đang căng bị xoay sẽ đổi hướng lực. Nếu chỉ chiếu kr rồi bỏ kg, ta đã bỏ một phần cơ học của chính mô hình hướng tâm này.
3. Hệ ba dây giả định
Chọn ba neo cách gốc 100 m, ở các góc 0°, 120°, 240°. Tọa độ tương ứng là (100; 0), (−50; 86,60254), (−50; −86,60254) m. Ở r = 0, mỗi dây có H₀ = 200 kN. Ba véc tơ lực bằng nhau và cách nhau 120° nên cân bằng khi chưa có tải ngoài.
Các độ cứng kr giả định lần lượt 20, 30, 10 kN/m. Chúng đại diện ba đặc tính khác nhau quanh cùng lực căng trước; không phải số liệu một hệ neo thực. Với R₀ = 100 m, tất cả kg = 2 kN/m. Sau khi cộng ba ma trận, ta nhận:
Kxx = Kyy = 33 kN/m
Kxy = Kyx = −8,660254 kN/m
Định thức bằng 1.014 (kN/m)², dương. Hai trị riêng là 24,339746 và 41,660254 kN/m, nên mô hình có độ cứng phục hồi dương theo mọi hướng ngang gần gốc. Điều này chỉ kiểm mô hình hai chiều ở trạng thái ban đầu; không kết luận về ổn định đứng, lắc, tải đỉnh hoặc trạng thái đứt dây của một công trình.
Hai phần tử chéo chính bằng nhau nhưng hệ vẫn ghép x–y. Dây ở 120° cứng hơn dây ở 240° làm hệ số ngoài đường chéo khác không. Vì vậy “độ cứng x bằng độ cứng y” không có nghĩa dịch chuyển sẽ luôn cùng hướng với tải.
4. Giải chuyển vị do tải ngoài
Đặt tải giả định Px = 120 kN, Py = 40 kN tại điểm nối. Mô hình tuyến tính cho:
33Δx − 8,660254Δy = 120
−8,660254Δx + 33Δy = 40
Nghiệm Δx = 4,246953 m và Δy = 2,326657 m. Nếu bỏ hệ số ghép và tính riêng 120/33, 40/33, sẽ nhận 3,636364 m và 1,212121 m, khác đáng kể. Không thể coi ma trận như hai lò xo độc lập khi các phần tử ngoài đường chéo chưa được chứng minh nhỏ.
Nếu chỉ bỏ độ cứng hình học kg nhưng giữ phần chiếu kr, ma trận có hai phần tử chéo bằng 30 và phần tử ghép vẫn −8,660254. Khi đó nghiệm là Δx = 4,783527 m, Δy = 2,714219 m. Đây là sai khác do bỏ một thành phần độ cứng, không phải do tuyến tính hóa khác thời điểm.
Phép kiểm đơn giản nhất là nhân lại K với chuyển vị để khôi phục tải. Tiếp đó kiểm các lực ban đầu có tổng bằng không. Nếu một cấu hình đã lệch cân bằng ngay từ đầu, chỉ giải KΔr = P mà bỏ lực dư sẽ dẫn tới một vị trí sai dù phép nhân ma trận hoàn toàn đúng.
5. Kiểm bằng mô hình hình học phi tuyến
Để kiểm ảnh hưởng đổi hướng, định nghĩa một mô hình giả định có đặc tính lực theo khoảng ngang tuyến tính:
Hᵢ = 200 + krᵢ(Rᵢ − 100)
Chỉ dùng trong vùng cả ba Hᵢ còn dương. Ở mỗi vị trí r, tính Rᵢ và hướng thực của dây, rồi giải:
P + ΣFᵢ(r) = 0
Mặc dù H theo R tuyến tính, véc tơ tổng lực vẫn phi tuyến theo x, y vì cả khoảng cách và hướng đều thay đổi. Đây chưa phải nghiệm catenary đầy đủ; nó là mô hình hướng tâm được định nghĩa rõ để kiểm phép tuyến tính hóa.
Với tải 120; 40 kN, nghiệm cân bằng là x = 4,225159 m, y = 2,472398 m. Các khoảng ngang là 95,806748; 100,091202; 104,281890 m. Lực ngang ba dây lần lượt 116,134963; 202,736053; 242,818896 kN. Cộng các véc tơ theo hướng mới cho đúng (−120; −40) kN, cân bằng tải ngoài.
Đưa nghiệm tuyến tính (4,246953; 2,326657) vào mô hình hình học này còn lực dư khoảng (−2,1811; +4,7233) kN. Đây là phép kiểm có ý nghĩa hơn chỉ thấy dịch chuyển “có vẻ nhỏ”. Tài liệu Aqwa mô tả lặp vị trí để giảm mất cân bằng, đồng thời lưu ý hệ phi tuyến có thể vượt quá nghiệm mong muốn nếu bước lặp không được kiểm soát. [3]
6. Độ nhạy khi tăng tải và mất dây
Nhân đôi tải thành (240; 80) kN. Tuyến tính hóa tại gốc cho x = 8,493906 m, y = 4,653315 m. Mô hình hình học phi tuyến cho x = 8,438758 m, y = 5,283954 m. Các lực tương ứng khoảng 34,271645; 204,167240; 288,954378 kN. Dây thứ nhất đã giảm lực đáng kể nhưng vẫn dương trong kịch bản này.
Nếu tải tiếp tục tăng làm một Hᵢ tính ra âm, dây không được biến thành thanh chịu nén. Phải chuyển sang nhánh chùng theo mô hình đã định nghĩa, giải lại cân bằng và cập nhật độ cứng. Không nên kéo dài đường tuyến tính ban đầu qua điểm mất căng rồi coi kết quả là lực neo hợp lệ.
Mất một dây cũng làm thay đổi cả lực ban đầu lẫn ma trận. Ví dụ bỏ dây góc 120° tại gốc khiến tổng lực hai dây còn lại bằng (100; −173,2051) kN trước khi cộng tải môi trường. Vì vậy tình huống mất dây không được tính chỉ bằng cách xóa một ma trận rồi dùng nguyên tải ngoài; phải đưa mất cân bằng do lực căng trước vào bài toán và tìm một trạng thái mới. Bài này chưa tính tải quá độ ngay sau sự cố.
7. Áp dụng vào mô hình vật nổi
Với fairlead đặt xa trọng tâm, lực dây còn gây mômen; một góc quay nhỏ dịch chuyển từng fairlead khác nhau. Aqwa lưu ý việc chuyển ma trận về tâm vật phải xét thêm thay đổi mômen khi vật quay. [1] Không được ghép các độ cứng tại nhiều điểm bằng phép cộng trực tiếp nếu chúng chưa cùng gốc và cùng bậc tự do.
Hồ sơ tính nên giữ vị trí neo/fairlead, lực căng trước, đặc tính H(R) hoặc mô hình dây, các tọa độ bị khóa, ma trận và lực dư sau kiểm lại. Chuyển vị ngang chỉ là một đầu ra; lực từng dây và điều kiện còn căng cần được kiểm cùng lúc. Như vậy kết quả giải hệ sẽ gắn với cấu hình cơ học cụ thể thay vì chỉ là một phép chia tổng tải cho “độ cứng neo”.
Nguồn tham khảo
[1] Aqwa Theory Manual 2024 R2 — The Stiffness Matrix — Ansys. Truy cập ngày 09/11/2026.
[2] MoorPy Usage — System Equilibrium and Stiffness Matrices — MoorPy developers, National Renewable Energy Laboratory. Truy cập ngày 09/11/2026.
[3] Aqwa Theory Manual 2024 R2 — Iteration Towards Equilibrium — Ansys. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
