Bài viết Tính toán định tâm hệ trục tàu thủy bằng phương pháp phần tử hữu hạn tập trung nghiên cứu xây dựng chương trình tính toán định tâm hệ trục theo tải trọng gối đỡ bằng phương pháp phần tử hữu hạn nhằm mục đích nội địa hóa phần mềm tính toán.
TẠP CHÍ ISSN: 1859-316X KHOA HỌC CƠNG NGHỆ HÀNG HẢI KHOA HỌC - CÔNG NGHỆ JOURNAL OF MARINE SCIENCE AND TECHNOLOGY TÍNH TỐN ĐỊNH TÂM HỆ TRỤC TÀU THỦY BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN SHIP PROPULSION SHAFTING ALIGNMENT USING FINITE ELEMENT METHOD LÊ ĐÌNH DŨNG1, PHẠM XUÂN DƯƠNG2, CAO ĐỨC THIỆP1* Khoa Máy tàu biển, Trường Đại học Hàng hải Việt Nam Trường Đại học Hàng hải Việt Nam *Email liên hệ: thiepcd@vimaru.edu.vn Tóm tắt Quá trình định tâm hệ trục tàu thủy ảnh hưởng đến hoạt động an toàn tin cậy hệ trục trình khai thác vận hành tàu thủy Hiện nay, có nhiều phương pháp khác sử dụng để tính tốn định tâm hệ trục tàu thủy Phổ biến nay, phương pháp định tâm hệ trục tàu thủy theo tải trọng gối đỡ phương pháp phần tử hữu hạn Trên giới, xuất phần mềm tính tốn định tâm hệ trục phương pháp trên, nhiên để sở hữu phần mềm địi hỏi chi phí cao Chính vậy, nghiên cứu này, tác giả tập trung nghiên cứu xây dựng chương trình tính tốn định tâm hệ trục theo tải trọng gối đỡ phương pháp phần tử hữu hạn nhằm mục đích nội địa hóa phần mềm tính tốn Kết tính tốn thử nghiệm cho tàu cụ thể nghiên cứu cho kết đảm bảo độ xác so sánh với kết tính tốn tính nhà máy đóng tàu Đây sở để phát triển phần mềm nội địa tính tốn định tâm hệ trục tàu thủy tương lai Từ khóa: Hệ trục tàu thủy, định tâm, tải trọng gối đỡ, phần tử hữu hạn Abstract The propulsion shaft alignment affects the safe and reliable operation of the ship shaft system during its exploitation Currently, there are many different methods used to compute the shafting alignment The most popular method is the shafting alignment according to the bearing reaction using the finite element method In the world, these software have appeared to compute shafting alignment, but they require very high costs Therefore, in this study, the authors focus on building a program to compute the shaft alignment according to the bearing reaction using the finite element method The specific ship's computation results are compared with the SỐ 71 (8-2022) shipyard's results to validate the program's accuracy This study is the primary step in developing domestic shaft alignment software Keywords: Propulsion, Shaft Alignment, Bearing Reaction, Finite Element Method Mở đầu Q trình tính tốn định tâm hệ trục tàu thủy trình tính tốn liên quan đến xác định đường tâm trục, lắp ráp thành phần hệ trục lên tàu chỉnh theo yêu cầu đưa Với tàu cụ thể, lên đà sửa chữa, tài liệu liên quan tới đường tâm trục tiêu chuẩn chỉnh hệ trục thường sẵn Do đó, đo độ dịch tâm (GAP) gãy khúc (SAG) hai bích nối hai đoan trục thường áp dụng phổ biến [1, 2] Phương pháp đo độ dịch tâm gãy khúc lặp lặp lại cho hệ trục lần sửa chữa lên đà tàu Với tàu biển đóng mới, tùy thuộc vào lực nhà máy đóng tàu mà phương pháp chỉnh hệ trục cách đo dịch tâm gãy khúc thực Tuy nhiên, phương pháp đo dịch tâm gãy khúc tồn nhược điểm áp dụng cho tàu biển đóng [3-5], bao gồm: chưa quan tâm tới trọng lượng thân hệ trục tải trọng tác dụng lên gối đỡ, đặc biệt gối đỡ trục chân vịt Theo phương pháp này, hệ trục tàu thủy thường mơ hình hóa hệ trục thẳng khơng bị biến dạng Do đó, phương pháp dẫn tới kết định tâm thiếu xác Phương pháp định tâm theo tải trọng gối đỡ nghiên cứu nhiều nhà khoa khoa học giới công ty hay tổ chức đăng kiểm [6-9] áp dụng phương pháp đo độ căng (strain gauge measurement) để đưa tải trọng gối đỡ cho loạt tàu Tại Việt Nam, phương pháp tải trọng gối đỡ số tác giả nghiên cứu áp dụng Tuy nhiên, nghiên cứu dựa phần mềm chun dụng có sẵn khó chủ động tính tốn 11 TẠP CHÍ KHOA HỌC - CÔNG NGHỆ Hiện nay, quan đăng kiểm khác giới Đăng kiểm Nhật Bản NK [7], Đăng kiểm Mỹ ABS [6],… đưa vào Quy phạm việc yêu cầu bảng tính định tâm hệ trục phải phê duyệt trước tiến hành lắp ráp hệ trục tàu thủy Trong đó, phương pháp định tâm theo tải trọng gối đỡ có nhiều ưu điểm [10, 11], khuyến khích áp dụng phổ biến lĩnh vực thiết kế sửa chữa tàu thủy Các cơng ty thiết kế tàu thủy nước ngồi thường tự xây dựng phát triển phần mềm tính tốn định tâm hệ trục tàu thủy cho riêng để sử dụng nội mua phần mềm chuyên dụng hãng phát triển phần mềm Tuy nhiên, chi phí để sở hữu phần mềm chuyển giao kỹ sử dụng phần mềm cao, tốn nhiều thời gian không chủ động Ở Việt Nam nay, công ty thiết kế tàu thủy hầu hết khơng tính tốn bảng tính định tâm hệ trục chưa có cơng cụ tính tốn định tâm theo tải trọng gối đỡ Khi cần phải có bảng tính định tâm hệ trục để trình Đăng kiểm phê duyệt Đối với tàu biển có đường kính trục 400mm [10, 11] họ phải th cơng ty nước ngồi giải phần việc dẫn đến chi phí đắt đỏ đồng thời bị động suốt trình thi cơng lắp ráp hệ trục Chính vậy, việc nghiên cứu phương pháp tính tốn định tâm hệ trục phương pháp phần tử hữu hạn đáp ứng yêu cầu thực tiễn quan đăng kiểm tăng cường lực cho ngành đóng tàu Việt Nam thiết kế tàu cỡ lớn Phương pháp nghiên cứu 2.1 Mơ hình tính tốn theo phương pháp phần tử hữu hạn Ở tốn tính nghiệm bền hệ trục tàu thủy, xây dựng mơ hình tính tốn ta cần xét đoạn trục tương quan lắp ghép với thiết bị hệ trục mà không cần xét đến trục khuỷu động diesel Tuy nhiên, tốn tính định tâm hệ trục, độ co bóp trục khủy ảnh hưởng lớn đến độ dịch tâm gãy khúc toàn hệ trục nên cần phải đưa trục khuỷu vào mơ hình tính Khi trục khuỷu tính tốn quy đổi thành đoạn trục thẳng tựa nhiều gối đặt tâm cổ trục động Do vậy, mơ hình tính định tâm hệ trục tàu thủy mơ tả Hình ISSN: 1859-316X KHOA HỌC CÔNG NGHỆ HÀNG HẢI JOURNAL OF MARINE SCIENCE AND TECHNOLOGY Phản lực độ dịch chỉnh gối đỡ xác định theo công thức sau: (1) R = A T Trong đó: R = R1R2 .R9 T = 1 Để rời rạc hóa mơ hình cho Hình thành mơ hình phần tử hữu hạn sẵn sàng cho tính tốn sử dụng phần tử dầm chịu uốn nút, gọi phần tử dầm Bernoulli-Euler biểu diễn Hình Tại nút phần tử có bậc tự gọi chuyển vị nút, độ võng v(x) góc xoay (x) Hình Phần tử dầm chịu uốn Bernoulli-Euler Với quy ước chiều chuyển vị đại lượng phần tử dầm Bernoulli-Euler [2] trình bày Hình 2, ta thành lập phương trình phần tử dạng ma trận sau: Keqe = Pe (2) T Trong đó, we = v v véc tơ 1 2 chuyển vị nút phần tử, Ke ma trận độ cứng phần tử, Pe véc tơ tải phần tử, tính sau: T K e = B BDdV Ve T T Pe = N g e dV + N Pe dS Ve Se với N = N1 N N N (3) (4) gọi ma trận hàm dạng Hermite, có dạng sau: Hình Mơ hình tính toán định tâm hệ trục tàu thủy 12 SỐ 71 (8-2022) TẠP CHÍ ISSN: 1859-316X KHOA HỌC CƠNG NGHỆ HÀNG HẢI KHOA HỌC - CÔNG NGHỆ JOURNAL OF MARINE SCIENCE AND TECHNOLOGY N =1− x N =3 L 3 x x x + 2 + ; N2 = x − L L L L 3 x x x − ;N = − + L L L x (5) Trong B = N gọi ma trận tính biến dạng Thực tích phân cơng thức (3) (4) với hàm dạng (5) thu ma trận độ cứng phần tử dầm Bernoulli-Euler sau: T T K e = B BDdV = E B BdFdx Ve LF = 12 EI 0 L dx 6L 4L 0 6L −6 L L 12 −6 L 2 4L −12 (6) Bước cần xây dựng phương trình tồn cục hay toàn kết cấu sử dụng nguyên lý bảo tồn lượng phương pháp Galerkin Đó hệ phương trình đại số tuyến tính bậc nhất, viết dạng ma trận sau: q − P = (7) Phương trình (7) chứa véc tơ ẩn chuyển vị nút toàn kết cấu q Tuy nhiên, ta áp dụng nguyên lý để thành lập phương trình chưa áp đặt điều kiện biên cho kết cấu Do ma trận suy biến không tồn ma trận nghịch đảo Ta phải đưa thêm vào điều kiện biên để đảm bảo hệ có nghiệm Phương trình (7) phương trình đại số tuyến tính có nhiều phương pháp giải khác như: phương pháp xác phương pháp lặp Những phương pháp trình bày cách hệ thống, chi tiết tài liệu liên quan đến đại số tuyến tính phương pháp tính Trong hầu hết gói phần mềm phân tích phần tử hữu hạn, việc tìm lời giải hệ phương trình tuyến tính lập trình sẵn hàm, thủ tục tiêu chuẩn có sẵn để sử dụng 2.2 Sơ đồ thuật tốn tính định tâm hệ trục theo tải trọng gối đỡ Để tính tốn định tâm hệ trục theo tải trọng gối đỡ phương pháp phần tử hữu hạn cần xây dựng chương trình tính mơi trường Matlab với tên gọi MATSHAL Sơ đồ thuật tốn tính tốn định tâm SỐ 71 (8-2022) thể Hình Theo đó, khối hình chữ nhật biểu thị bước cơng việc nhiệm vụ quan trọng cần thực hiện, khối hình elip thể liệu biểu diễn ma trận, tệp m-file đồ thị, hình thoi thể bước kiểm tra theo liệu kiểu Boolean (đúng-sai), mũi tên biểu thị luồng vận động chương trình Chức khối Nhập liệu để nhập liệu cần thiết cho chương trình tính dạng tệp văn bao gồm liệu hình học, kết cấu, tính vật liệu hệ trục cần tính tốn Khối Tính độ dịch chuyển gối đỡ ban đầu (offset) thực tính tốn độ dịch chuyển gối cần hiệu chỉnh ban đầu gối đỡ hệ trục để đảm bảo theo tiêu chí tải trọng gối đỡ Khối Tính phản lực gối ban đầu thực tính phản lực gối đỡ với giả thiết gối đỡ có độ dịch chuyển ban đầu tính Khối Xác định độ hiệu chỉnh gối đỡ thực xác định hiệu chỉnh cần thiết gối để đảm bảo tải trọng phân bố gối đỡ đạt theo tiêu chí đặt Khối Phản lực gối đỡ sau hiệu chỉnh thực tính phản lực gối đỡ sau hiệu chỉnh độ lệch gối theo giá trị hiệu chỉnh gối tính Khối Kiểm tra điều kiện tải trọng gối đỡ thực kiểm tra phản lực gối đỡ thỏa mãn điều kiện tải trọng gối đỡ, đạt yêu cầu tiến hành xuất liệu Hình Sơ đồ thuật tốn tính định tâm hệ trục theo tải trọng gối đỡ Toàn q trình tính tốn, bao gồm tính độ lệch gối, phản lực gối đỡ, độ hiệu chỉnh gối, mô men uốn, lực cắt, ứng suất,… gối đỡ toàn chiều dài hệ trục thực với việc sử dụng lý thuyết dầm Becnulli-Euler phương pháp 13 TẠP CHÍ ISSN: 1859-316X KHOA HỌC CÔNG NGHỆ HÀNG HẢI KHOA HỌC - CÔNG NGHỆ JOURNAL OF MARINE SCIENCE AND TECHNOLOGY phần tử hữu hạn cho toán dầm nhiều nhịp nhằm đảm bảo cho kết tin cậy, xác Tính tốn thử nghiệm mơ hình tàu 45.000DWT 3.1 Mơ hình hóa hệ trục tàu 45.000DWT Nhằm mục đích kiểm nghiệm chương trình tính tốn định tâm hệ trục theo tải trọng gối đỡ phương pháp phần tử hữu hạn trình bày trên, trường hợp tính tốn cho tàu 45.000DWT sử dụng để kiểm tra độ xác chương trình tính Hệ trục tàu chở dầu/hóa chất 45.000DWT, tốc độ tối đa đạt 15,3knots, chạy biển cấp không hạn chế kiểm tra giám sát đóng phân cấp đăng kiểm NK Cấu trúc bố trí chung tàu thể Hình Dựa vào tài liệu thiết kết hệ trục tàu, cần tiến hành xây dựng mơ hình hệ trục bao gồm kích thước, khối lượng phục vụ cho q trình tính định tâm hệ trục Mơ hình bố trí hệ trục mơ tả Hình Để tính tốn định tâm hệ trục tàu bước cần xây dựng sơ đồ hệ lực tác dụng lên hệ trục Đối với hệ trục khảo sát bao gồm 01 trục chân vịt, 01 trục trung gian 01 trục khuỷu Số lượng gối đỡ: tổng số gối đỡ 09, bao gồm gối đỡ trục chân vịt, 01 gối đỡ trục trung gian, 06 gối đỡ trục khuỷu để đảm bảo yêu cầu số lượng gối đỡ bên máy khơng nhỏ [6, 7] Sơ đồ ngoại lực tác dụng lên hệ trục thể Hình Trong tính tốn định tâm hệ trục theo tải trọng gối đỡ cần phải tính tốn tải trọng gối đỡ gối đỡ Hình Bố trí chung tàu 45000DWT Hình Mơ hình bố trí hệ trục Hình Sơ đồ ngoại lực tác dụng lên hệ trục 14 SỐ 71 (8-2022) TẠP CHÍ ISSN: 1859-316X KHOA HỌC CÔNG NGHỆ HÀNG HẢI KHOA HỌC - CÔNG NGHỆ JOURNAL OF MARINE SCIENCE AND TECHNOLOGY nằm đường tâm lý thuyết trạng thái dằn tàu, chân vịt nằm hoàn toàn nước với ký hiệu TT0 (trạng thái 0) Từ lập ma trận hệ số ảnh hưởng xác định gối đỡ có phản lực âm tức có chiều hướng từ lên Có trạng thái tàu hệ trục cần tính tốn định tâm theo gợi ý quan đăng kiểm bao gồm: Trạng thái (ký hiệu TT1) tàu xét trạng thái dằn, chân vịt nằm hoàn toàn nước gối đỡ điều chỉnh theo trạng thái nóng; trạng thái (ký hiệu TT2) xét tàu trạng thái dằn, chân vịt nằm hoàn toàn nước gối đỡ điều chỉnh theo trạng thái lạnh; trạng thái (ký hiệu TT3) tương ứng với trạng thái lắp đặt hệ trục tàu trạng thái chân vịt nửa chìm nước gối đỡ điều chỉnh theo trạng thái lạnh Trong q trình tính tốn cần xem xét lưu ý đến thay đổi chiều cao gối đỡ nhiệt độ thay đổi gối đỡ máy Thơng thường giá trị cung cấp hồ sơ máy gối đỡ khơng cịn có tác dụng đỡ trục Bên cạnh đó, hệ trục vị trí gối đỡ không tiếp xúc với gối đỡ nửa bạc Do đó, hệ trục hoạt động gây dao động, tăng ứng suất, gây phá hủy gối đỡ nhanh chóng Rõ ràng, điều không thỏa mãn yêu cầu thiết kế lắp ráp hệ trục cần thay đổi khoảng offset để đảm bảo yêu cầu vận hành hệ trục 3.2 Kết tính tốn kiểm chứng Trước tiên, tính tốn tải trọng tác dụng lên gối đỡ xét gối đỡ nằm thẳng hàng đường tâm lý thuyết tương ứng với độ dịch chỉnh offset = ngoại lực tác dụng lên hệ trục trạng thái dằn tàu tương ứng với chân vịt ngập hoàn toàn nước biển Biểu đồ mô men uốn, lực cắt độ võng hệ trục TT0 mơ tả Hình Giá trị tải trọng tác dụng lên gối đỡ trường hợp gối đỡ nằm đường tâm lý thuyết thể Bảng Bảng Bảng giá trị tải trọng gối đỡ hệ trục trạng thái gối đỡ nằm đường tâm lý thuyết Bearing Position (m) No Force Reaction (T) 2,46 29,562 Hình Biểu đồ nội lực hệ trục TT0 5,42 -1,773 10,685 8,467 15,619 13,259 16,374 -2,995 17,224 10,496 18,074 9,660 18,924 8,197 19,774 11,089 Để thành lập ma trận ảnh hưởng hệ trục điều kiện không kể đến tải trọng thân hệ trục tải trọng tác dụng cần tiến hành nâng gối đỡ từ số đến số khoảng 0,1mm lấy giá trị tải trọng gối đỡ Ma trận ảnh hưởng phản ánh độ cứng hệ trục sở để tính tốn xác định khoảng dịch chuyển gối đỡ cho trạng thái định tâm Theo bảng giá trị tải trọng gối đỡ khác phân bố tải trọng lên hệ trục không vị trí gối đỡ Xuất phản lực gối đỡ gối đỡ số số mang giá trị “âm” điều có nghĩa phản lực gối đỡ có chiều từ xuống, tức SỐ 71 (8-2022) Như trình bày có trạng thái hệ trục cần thiết phải tính tốn định tâm gối đỡ để đảm bảo điều kiện phân bố tải cần dịch chuyển offset Bảng tổng hợp độ dịch chuyển gối đỡ định tâm trạng thái khác 15 TẠP CHÍ ISSN: 1859-316X KHOA HỌC CƠNG NGHỆ HÀNG HẢI KHOA HỌC - CÔNG NGHỆ JOURNAL OF MARINE SCIENCE AND TECHNOLOGY Bảng Bảng tổng hợp độ dịch chuyển offset cho Bảng Bảng tổng hợp phản lực gối đỡ đối gối đỡ đối trạng thái khác trạng thái khác Bearing Offset (m) Bearing Force Reaction (T) No TT1 TT2 TT3 No TT1 TT2 TT3 0,078 0,078 0,078 26,577 1,825 28,243 1,272 28,049 1,141 -1,1 -1,96 -1,196 -2,2 -1,196 -2,2 7,769 5,978 8,155 4,103 8,172 4,089 -1,96 -1,96 -2,2 -2,2 -2,2 -2,2 3,760 10,181 5,303 10,140 5,303 0,140 -1,96 -1,96 -2,2 -2,2 -2,2 -2,2 9,872 8,160 9,971 8,160 9,981 8,172 -1,96 -2,2 -2,2 11,121 11,125 11,122 Bảng Bảng tính tải trọng gối đỡ Bảng Bảng so sánh kết tính tốn phản lực gối đỡ trước thêm vật nặng MATSHAL nhà máy đóng tàu Force Reaction (T) Bearing Position No (m) Force Reaction (T) Bearing No SHIP MATSHAL Error (%) BUILD COM 2,46 29,629 26,577 26,547 0,11 1,8245 1,85 1,37 7,7692 7,852 1,05 5,42 -2,360 10,685 11,955 15,619 -2,091 5,978 5,939 0,65 3,7602 3,62 3,87 16,374 9,625 17,224 10,144 10,18 10,012 1,67 9,9716 10,359 3,73 18,074 9,970 18,924 8,160 8,1603 8,34 2,15 19,774 11,121 11,121 11,34 1,93 Sau xác định giá trị khoảng dịch chuyển gối đỡ offset cần tiến hành tính tốn kiểm tra lại phản lực gối đỡ Bảng tổng hợp phản lực gối đỡ trạng thái định tâm khác Để xác định giá trị SAG GAP (độ gãy khúc độ dịch tâm) cho trình lắp ráp hệ trục vị trí bích nối trục chân vịt trục trung gian trạng thái tự bảng tính phản lực gối đỡ thể Bảng Theo Bảng 4, tải trọng gối đỡ số có giá trị âm (-2.360 T) dựa vào giá trị phải thêm tải trọng đặt bích nối để đảm trục chân vịt tiếp xúc với gối trục: Chọn thử giá trị tải đặt lên 2,1 Sau đặt thêm tải trọng phụ thêm này, phản lực gối đỡ số dương với giá trị 0,572T đạt yêu cầu Khi độ võng góc xoay bích nối trục chân vịt 0,213951mm 0,000176rad Đối với trục trung gian, cần bổ sung thêm gối đỡ giả hai vị trí mơ hình tính trạng thái tự 16 Tại gối đỡ trục trung gian độ võng 0,100962mm điều chỉnh gối đỡ giả cho vị trí trục gối trục đạt giá trị -1,196mm Khi độ võng góc xoay bích nối trục trung gian đầu nối với trục chân vịt nối với trục khuỷu -1,933867mm, 0,000265rad -1,855155mm, 0,000245rad Đối với trục khuỷu độ võng góc xoay bích nối trục khuỷu là: 2,206939mm -0,000014rad Từ giá trị ta xác định giá trị SAG, GAP bích nối trục chân vịt trục trung gian: GAP | 0, 213951 ( 1, 933867) | 1, 72 SAG 880x(0, 000176 0, 000265) 0, 078 Tính giá trị SAG, GAP bích nối trục khuỷu động trục trung gian: GAP | 1, 855155 2, 206939 | 0, 35 SAG 900x(0, 000265 0, 000014) 0, 22 SỐ 71 (8-2022) TẠP CHÍ ISSN: 1859-316X KHOA HỌC CÔNG NGHỆ HÀNG HẢI JOURNAL OF MARINE SCIENCE AND TECHNOLOGY Để kiểm tra đánh giá phương pháp tính định tâm theo tải trọng gối đỡ sử dụng chương trình MATSHAL xây dựng, kết tính tốn so sánh với kết tính nhà máy đóng tàu Nhà máy đóng tàu thuê phần mềm chuyên dụng quan đăng kiểm phê duyệt quan đăng kiểm NK Với kết cấu hệ trục tính trường hợp hệ trục trạng thái nóng tàu trạng thái dằn với giá trị offset So sánh tải trọng gối đỡ, mơ men uốn, góc xoay tính phần mềm MATSHAL nhà máy đóng tàu thể Bảng Hình Hình So sánh kết tính tốn MATSHAL nhà máy đóng tàu Qua biểu đồ so sánh ta thấy kết tính tốn định tâm hệ trục chương trình tính tốn định tâm MATSHAL phù hợp với kết tính tốn đề xuất nhà máy đóng tàu Sai số hai trường hợp tính khơng 4%, đảm bảo độ xác tính tin cậy phần mềm MATSHAL Kết luận Bài báo xây dựng chương trình tính MATSHAL thực tính tốn số liệu cần thiết cho trình định tâm hệ trục tàu thủy theo tiêu chí tải trọng gối đỡ phương pháp phần tử hữu hạn cho tàu thực tế, tàu chở hóa chất 45.000DWT Tác giả tiến hành so sánh, phân tích số liệu tính tốn từ chương trình vừa xây dựng với kết từ nhà SỐ 71 (8-2022) KHOA HỌC - CƠNG NGHỆ máy đóng tàu phê duyệt quan đăng kiểm NK để khẳng định tính xác, khả thi chương trình tính Cụ thể, biến thiên số liệu hai chương trình tính thống sai số giá trị nhỏ 4% Vì vậy, chương trình MATSHAL sử dụng thực tiễn tính toán định tâm hệ trục tàu thủy TÀI LIỆU THAM KHẢO [1] Đặng Hộ (1985), Thiết kế trang trí động lực tàu thủy, Tập 1+2, NXB Giao thông vận tải [2] Võ Như Cầu, Tính kết cấu theo phương pháp phần tử hữu hạn, NXB Xây dựng, Hà Nội [3] Nguyễn Mạnh Thường (2011), Tính tốn tải trọng gối trục chong chóng tàu thủy, Tạp chí Giao thơng vận tải, Số tháng 6/2011, Hà Nội [4] Phạm Quốc Việt (2015), Nghiên cứu ứng dụng phần mềm Lira để tính toán kiểm tra định tâm lắp ráp hệ trục đóng tàu, Đề tài KHCN cấp Cơ sở, Trường Đại học Hàng hải Việt Nam [5] Phan Trung Kiên (2013), Nghiên cứu phương pháp tính tốn định tâm hệ trục tàu thủy theo tải trọng gối đỡ phần mềm SAP, Luận văn Thạc sỹ Trường Đại học Hàng hải Việt Nam [6] Batra, A., K Shankar, and S Swarnamani (2007), Propulsion shaft alignment measurements on warships afloat and alignment solution using multi-objective optimisation, Proceedings of IMarEST-Part A-Journal of Marine Engineering and Technology, pp.39-49 [7] SHIRAKI, D., et al (2005), Advanced Calculation Method for Marine Propulsion Shafting Alignment, ClassNK technical bulletin [8] Cowper, B., A DaCosta, and S Bobyn (1999), Shaft alignment using strain gages: Case studies Marine Technology and SNAME News, Vol.36(02): pp.77-91 [9] Vulić, N., A Šestan, and V Cvitanić (2006), Shafting alignment calculation and validation criteria Trogir/Split [10] Lê Đình Dũng (2021), Các phương pháp tính toán định tâm hệ trục tàu thủy theo tải trọng gối đỡ, Chuyên đề Tiến sĩ số 1, Trường ĐH Hàng hải Việt Nam [11] Chang, M H., Juang, S H (2014), Theoretical Design Study on Shafting Aligment Calculation for High Speed Craft Ngày nhận bài: Ngày nhận sửa: Ngày duyệt đăng: 08/8/2022 17/8/2022 22/8/2022 17 ... q trình tính định tâm hệ trục Mơ hình bố trí hệ trục mơ tả Hình Để tính tốn định tâm hệ trục tàu bước cần xây dựng sơ đồ hệ lực tác dụng lên hệ trục Đối với hệ trục khảo sát bao gồm 01 trục chân... nghiên cứu phương pháp tính toán định tâm hệ trục phương pháp phần tử hữu hạn đáp ứng yêu cầu thực tiễn quan đăng kiểm tăng cường lực cho ngành đóng tàu Việt Nam thiết kế tàu cỡ lớn Phương pháp nghiên... cứu 2.1 Mơ hình tính tốn theo phương pháp phần tử hữu hạn Ở tốn tính nghiệm bền hệ trục tàu thủy, xây dựng mơ hình tính tốn ta cần xét đoạn trục tương quan lắp ghép với thiết bị hệ trục mà không